资源描述：
《数值计算在matlab的插值应用》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、例3.16取余弦曲线上11个点的自变量和函数值点作为已知数据，再取41个自变量点。计算插值函数值①打开matlab②在CommandWindow里输入(u>>?,matlab自带的)»x=0:10;y=cos(x);回车，对应在Workspace里创建数据，即11组已知数据那Stack:SelectdatatoplotValue<1x11double><1x11double>MinMax010-0.99...1③计算精确值做比较(“％”为叙述语言，忽略它)»xi=0:.25:10;»yO=cos(xi);%精确值Nam
2、e八Workspace*□0X團F如耳■StacieSelectdatatoplot▼ValueMinMax<1x11double>010<1x41double>010<1x11double>■0.99...1<1x41double>-0.99...1④输入三种插值方式，计算»y1=interp1(x,y,xi);»y2=interpl(x，y，xi,cubic‘)；»y3=interpl(x,y,xi/spline,);%线性插值结果%三次方程式插值结果%样条插值结果Workspace七□，XStack:口)Sele
3、ctdatatoplot▼Name▲ValueMinMax王X<1x11double>010xi<1x41double>010王y<1x11double>•0.99.x1=3y0<1x41double>•0.99一1S<1x41double>•0・99_1Sy2)<1x41double*-0.99^1田y37<1x41double>-1.01^.1⑤绘图(plot)»subplot(2,3,l)；plot(xi,yO,'o',xi,y1);title(*l.与线性插值')Shortcuts劃HowtoAdd劃What'
4、sNewCommandWindowMFigures-Figure1»»»»»»»A»x=0:10;y=cos(x):xi=0:.25:10：y0=cos(xi):%精确值yl=interpl(x,y,xi):%线性插值结果y2=interpl(x,y,xi,cubic"):%三次方程式插值结果y3=interpl(m,y,Mi,'spline'):%样条恬值结果__(sxSplot:plot仗i,y0,、o',xi,yl)5、ndow»»»»»»»»>>“□ox05100-0.Figures・figure1□曰HI◎
6、h
7、铁d紗⑥•必⑤对误丼进行计算，表示线性插值与精确值的误丼»subplot(2,3,4)；plot(xi,yl-yO);title(f差1‘)x=0:10;y=cos(x):xi=0:.25:10:yO=cos(xi):%精确值yl=interpl(x,y,xi):%线性插值结果y2=interpl(x,y,xijcubic'):%三次方程式插值结果y3=interpl(x,y,xi/spline?):%样条插值结果subp
8、lot⑵3,1):plot(xi,yO^J^o",xi,yl):title('1.与线性插值')subplot(2,3,4)^plot7xi?yl-yOr?5it1e('差1')⑥输入所有插值命令，并绘图CommandWindow»x=0:10;y=cos(x):»xi=0:.25:10：»y0=cos(xi)■.%精确值»yl=interpl(x,y,xi):%线性插值结果»y2=interpl(x?y,xi?'cubic-'):%三次方程式插值结果»y3=interpl(x?y,xi,'spline'):%样条插值
9、结果»subplot(2,3,1):plot(xi,yO,o',xi,yl);titleC1.与线性插值')>>subplot(2,3,4):plot(xi?yl-yO);title('差1')»subplot(2,3,2):plot(xi,yO,o',xi,y2,'):title('2.与三次方程式插值')>>subplot(2,3,5):plot(xi,y2-y0);title('差2’)»subplot(2,3,3):plot(xi}y0,o',xi,y3);titleC3.与样条插值')»subplot(2j3
10、j6):plot(xi,y3-y0);titleC差3’)◎IAJ
11、%'響®駁尿•多UILI■口1与线性插值2与三次方程式插值差1差2差3»subplot(2,3,2);plot(xi,y0；o,,xi,y2；-.,);title(,2.与三次方程式插值')»subplot(2,3,5);plot(xi,y2-yO);titleC
