资源描述：
《用MATLAB解插值和曲线拟合问题》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、注：上机作业文件夹以自己的班级姓名学号命名，文件夹包括如下上机报告和Matlab程序。上机报告模板如下：佛山科学技术学院上机报告课程名称数学应用软件上机项目用MATLAB解插值和曲线拟合问题专业班级一.上机目的熟练掌握多种插值方法：线性插值，三次样条插值，三次插值和最近邻点插值(linear、spline.cubicsnearest)Interpl(x,y,ex,'method')在MATLAB的线性最小二乘拟合中，用得较多的是多项式拟合，其命令为：A=polyfit(x,y,m)二.上机内容1、在飞机的机翼加工时，由于机翼尺寸很大，通常在
2、图纸上只能标出部分关键点的数据.某型号飞机的机翼上缘轮廓线的部分数据如下。用插值法求x每改变0」时的y值，画出图形表示。x04.749.051938577695114133152171190y05.2311.9716.1517」16.3414.63」12」66.697.033.9902、已知观测数据点如表所示：x00.10.20.30.40.50.60.70.80.91丿-0.4471.9783.286.167.087.347.669.569.489.311.2用3次多项式函数拟合这些数据点，画出图形。3^教材习题1,用griddata插值
3、函数，三次插值。4、教材习题4,初始值为V()=0,T=0.05三.上机方法与步骤第1题：要用插值法教每改变0」时的y值，可以选用程序y=interpl(xheng,yzhong,x,'spline1);xheng=[O4.749.051938577695114133152171190];yzhong=[05.238.111.9716.1517.116.3414.6312.166.697.033.990];x=0:0.1:190;y=interpl(xheng,yzhong,x,fsplinef);plot(xheng,yzhong,1+f
4、zx,y,xheng,yzhong,1r:1)xlabel(1xhen1),ylabel(1yzhon1)截图：第2题：要用3次多项式函数拟合这些数据点，并画出图形，选用程序A=polyfit(xzyz3),z=polyval(A,x);plot(x,y,*k+*,xzz,'r*)编程求解即可。x=0:0.1:1;y=[-0.4471.9783.286.167.087.347.669.569.489.311.2];A=polyfit(x,yz3)z=polyval(A,x);plot(x,y,'k+',x,z,*r1)截图：FileEdit
5、TextGoCellToolsDebugDesktopWindowHelp%c：t，增cl1.02-3_4_5一x=0:0.1:1;y=[-0.4471.9783.286.167.087.347.669.569.489.311.2]：A=polyfit(x,y,3)z=polyval(A,x):plot(x,yjk+，,x,z/r')f•••••feyijiyil.mx[sanciduoxiangshi2・mx]第3题：用griddata插值函数，三次插值，选用cK=griddata(T,P,K,cT,cP,*cubic*)编程。T=[68
6、688787106106140140];P=[9.798113.3249.007813.3559.791814.2779.656312.463];K=[0.08480.08970.07620.08070.06960.07530.06110.0651];cT=99;cP=10.3;cK=griddata(TzPzK,cTzcP,*cubic*)截图：'Editor・C:User5AdministratorDesktopxishuk3.m-口FileEditTextGoCellToolsDebugDesktopWindowHelp珞■越
7、町色輪旨*41♦♦他
8、IB■鳞FW□唔:♦回厂回•目b目1.0+-r1.1x恋疫Q1-T=[68688787106106140140];2-P=[9.798113.3249.007813.3559.791814.2779.656312.463];3-K=[0.08480.08970.07620.08070.06960.07530.06110.0651];4一cT=99;5-cP=10.3;6-cK=griddata(T,P,K,cT,cP,cubic7)8第4题：要确定t,v,选用x=lsqcurvefit(^functionf=curve
9、fun4(x^t)一f=10-(10-x(l))*exp(-t/x(2)):curvefun41,x0zt,v),f=curvefun4(xzt)functionf=curvef
