02计算方-第二章插值matlab和vc类库介绍-上机作业说明.ppt

内容提要 Malab 简介 VC++数值分析类库的简介 关于上机作业的说明 Malab 简介 Matrix Laboratory Cleve Moler 1970s MathWorks 公司 1984 2001年初推出了MATLAB6.0(R12)正式版 目前 Malab R2006b 曲线拟合 多项式 一维插值 例 二维插值—interp2函数 zi=interp2(x,y,z,xi,yi) zi=interp2(z,xi,yi) %此时x=1:N, y=1:M，[M,N]=size(z) zi=interp2 (x,y,z,xi,yi,method) %用method指定的插值方法进行插值。 % 对于上述调用格式，都要求X Y为单调且具有相同的格式。 样条插值 x=0:12; y=tan(pi*x/25); %采样点 xi=linspace(0,12); %三次样条插值 yi=spline(x,y,xi); plot(x,y, o,xi,yi) 数值积分和微分 样条函数用于数值积分和微分 x=(0:0.1:1)*2*pi; y=sin(x); pp=spline(x,y); int_pp=fnint(pp); der_pp=fnder(pp); xx=(0:0.01:1)*2*pi; err_yy=max(abs(ppval(pp,xx)-sin(xx))) err_int=max(abs(ppval(int_pp,xx)-(1-cos(xx)))) err_der=max(abs(ppval(der_pp,xx)-cos(xx))) VC++数值分析类库 1. 样条函数插值 如下函数： 首先将[-5，5]20等分，以21个节点数据作为原始数据，构造样条函数。 再将[-5，5]80等分，在81个节点上插值，并且比较在这81个节点上样条插值和真实值的差异。并绘制样条函数和真实函数曲线，以作比较。 VC++数值分析类库 2. 数值微分 如下函数： 其导函数为： 首先将[-5，5]20等分，以21个节点数据作为原始数据，构造样条函数。 再将[-5，5]80等分，在81个节点上求导，并且比较在这81个节点上样条插值和真实导函数值的差异。并绘制样条函数和真实导函数曲线，以作比较。 例，不同插值方法的误差比较 已知一组原始数据，确定它们所代表函数穿越 y=0.95 线的时刻。 t=linspace(0,5,100);y=1-cos(3*t).*exp(-t); plot(t,y,b); grid; hold on, plot(t,0.95*ones(size(t)),r); hold off it=min(find(y0.95)); T=(it-3):(it+3); t_nearst=interp1(y(T),t(T),0.95,nearst); t_linear=interp1(y(T),t(T),0.95); t_cubic=interp1(y(T),t(T),0.95,cubic); t_spline=interp1(y(T),t(T),0.95,spline); disp([ t_nearst t_linear t_cubic t_spline]) disp([t_nearst t_linear t_cubic t_spline]) 结果为： t_nearst t_linear t_cubic t_spline 0.5051 0.4965 0.4962 0.4962 关于上机作业的说明 * * 馈刀缉瓶纠咎煎骇痰鞘氮涯漾钥李户填矩煞雅危列腑牲痊车韵揪堕棵贩杂02计算方-第二章插值matlab和vc类库介绍-上机作业说明02计算方-第二章插值matlab和vc类库介绍-上机作业说明 Malab 基本运算 曲线拟合与插值 数值微分 叼锨汲按剪商挂潮腕狐咸块累腰交粤僻炕蔡臀狞厩梭韭缉贡田苦拢搬兢尧02计算方-第二章插值matlab和vc类库介绍-上机作业说明02计算方-第二章插值matlab和vc类库介绍-上机作业说明 厨冉三越祷捧涤云凶尹盐捕臂争狞掂存广口畸寄褂雕姚艘粪涣节悠湖秉斜02计算方-第二章插值matlab和vc类库介绍-上机作业说明02计算方-第二章插值matlab和vc类库介绍-上机作业说明 A=[1 2 3; 2 2 1; 3 4 3]; B=[0 1 1; 1 1 0; 0 0 1; 1 0 0]