数学应用软件作业9用Matlab解插值问题和曲线拟合问题.doc

注：上机作业文件夹以自己的班级姓名学号命名，文件夹包括如下上机报告和Matlab程序。 上机报告模板如下： 佛山科学技术学院 上 机 报 告 课程名称 数学应用软件 上机项目 用Matlab解插值问题和曲线拟合问题 专业班级 姓 名 学 号 一. 上机目的 本节课我们学习了用Matlab解插值问题和曲线拟合问题，主要有以下内容： 1、了解和掌握用matlab解决各种问题的插值方法。 2、学习用数据拟合函数，并会画出拟合后的图形。 二. 上机内容 1、在飞机的机翼加工时, 由于机翼尺寸很大, 通常在图纸上只能标出部分关键点的数据. 某型号飞机的机翼上缘轮廓线的部分数据如下。用插值法求x每改变0.1时的y值，画出图形表示。 x 0 4.74 9.05 19 38 57 76 95 114 133 152 171 190 y 0 5.23 8.1 11.97 16.15 17.1 16.34 14.63 12.16 6.69 7.03 3.99 0 2、已知观测数据点如表所示： 用3次多项式函数拟合这些数据点，画出图形。 3、教材习题1，用griddata插值函数，三次插值。 4、教材习题4 , 初始值为V0=0，=0.05 三. 上机方法与步骤 给出相应的问题分析及求解方法，并写出Matlab程序，并有上机程序显示截图。 题1： Matlab程序： 建立程序timu1.m如下： xs=[0 4.74 9.05 19 38 57 76 95 114 133 152 171 190]; ys=[0 5.23 8.1 11.97 16.15 17.1 16.34 14.63 12.16 6.69 7.03 3.99 0]; x=0:0.1:190; y=interp1(xs,ys,x,spline); plot(xs,ys,+,x,y,xs,ys,r) xlabel(x),ylabel(y) 题2：Matlab程序： 建立程序timu2.m，如下： x=0:0.1:1; y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.3 11.2]; A=polyfit(x,y,3) z=polyval(A,x); plot(x,y,k+,x,z,r) 题3： Matlab程序： 建立M程序timu3.m,如下： T=[68 68 87 87 106 106 140 140]; P=[9.7981 13.324 9.0078 13.355 9.7918 14.277 9.6563 12.463]; K=[0.0848 0.0897 0.0762 0.0807 0.0696 0.0753 0.0611 0.0651]; CT=99; CP=10.3; CK=griddata(T,P,K,CT,CP,cubic) 题4： 建立主程序timu41.m如下： function f=timu41(x,t) f=10-(10-x(1))*exp(-t/x(2)) 建立程序timu42.m如下： t=[0.5 1 2 3 4 5 7 9]; v=[6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63]; x0=[0,0.05]; x=lsqcurvefit(timu41,x0,t,v) f=timu41(x,t) 上机结果 题1结果： 题2结果： A = 16.0758 -33.9245 29.3246 -0.6104 题3结果： CK = 0.0724 题4结果： x = 5.5577 3.5002 f = 6.1490 6.6616 7.4913 8.1147 8.5832 8.9353 9.3987 9.6604 y x 0 -0.447 0.1 1.978 0.2 3.28 0.3 6.16 0.4 7.08 0.5 7.34 0.6 7.66 0.7 9.56 0.