资源描述:
《用MATLAB解插值和曲线拟合问题.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、注:上机作业文件夹以自己的班级姓名学号命名,文件夹包括如下上机报告和Matlab程序。上机报告模板如下:佛山科学技术学院上机报告课程名称数学应用软件上机项目用MATLAB解插值和曲线拟合问题专业班级一.上机目的熟练掌握多种插值方法:线性插值,三次样条插值,三次插值和最近邻点插值(linear、spline、cubic、nearest)Interp1(x,y,cx,’method’)在MATLAB的线性最小二乘拟合中,用得较多的是多项式拟合,其命令为:A=polyfit(x,y,m)二.上机内容1、在飞机的机翼加工时,由于机翼尺寸很大,通常在图纸上只能标出部分关键点的数
2、据.某型号飞机的机翼上缘轮廓线的部分数据如下。用插值法求x每改变0.1时的y值,画出图形表示。x04.749.051938577695114133152171190y05.238.111.9716.1517.116.3414.6312.166.697.033.9902、已知观测数据点如表所示:xy0-0.4470.11.9780.23.280.36.160.47.080.57.340.67.660.79.560.89.480.99.3111.2用3次多项式函数拟合这些数据点,画出图形。3、教材习题1,用griddata插值函数,三次插值。4、教材习题4,初始值为V0=
3、0,=0.05三.上机方法与步骤第1题:要用插值法求x每改变0.1时的y值,可以选用程序y=interp1(xheng,yzhong,x,'spline');xheng=[04.749.051938577695114133152171190];yzhong=[05.238.111.9716.1517.116.3414.6312.166.697.033.990];x=0:0.1:190;y=interp1(xheng,yzhong,x,'spline');plot(xheng,yzhong,'+',x,y,xheng,yzhong,'r:')xlabel('xhen')
4、,ylabel('yzhon')截图:第2题:要用3次多项式函数拟合这些数据点,并画出图形,选用程序A=polyfit(x,y,3),z=polyval(A,x);plot(x,y,'k+',x,z,'r')编程求解即可。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,y,'k+',x,z,'r')截图:第3题:用griddata插值函数,三次插值,选用cK=griddata(T,P,K,cT,cP,'cubic')编程
5、。T=[68688787106106140140];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(T,P,K,cT,cP,'cubic')截图:第4题:要确定t,v,选用x=lsqcurvefit('curvefun4',x0,t,v),f=curvefun4(x,t)functionf=curvefun4(x,t)f=10-(10-x(1))*exp(
6、-t/x(2));t=[0.51234579];v=[6.366.487.268.228.668.999.439.63];x0=[0,0.05];x=lsqcurvefit('curvefun4',x0,t,v)f=curvefun4(x,t)截图:四.上机结果第1题:截图:第2题:A=Columns1through316.0758-33.924529.3246Column4-0.6104截图:第3题:cK=0.0724截图:第4题:x=5.55773.5002f=Columns1through56.14906.66167.49138.11478.5832Column
7、s6through88.93539.39879.6604截图: