资源描述:
《应用软件9求插值问题和曲线拟合问题》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、注意:上机作业文件夹以自己的姓名学号命名,文件夹中包括如下上机报告和Matlab程序。佛山科学技术学院上机报告课程名称数学应用软件上机项口求插值问题和曲线拟合问题专业班级姓名学号一.上机目的1•了解各种插值的用法,一维插值的基本格式是interpl(xzyzex,'method*);二维插值的基本格式是interp2(x,y,z,ex,cy,'method*)和插值基点为散乱节点griddata(x,y,z,ex,cy,'method*);2.了解各种插值函数的method选择方法;3.了解曲线拟合的程序用法,线性最小乘拟合的程序polyval(A,x),非线性最小乘拟合的程序ls
2、qcurvefit和lsqnonlin。二.上机内容1、在飞机的机翼加工时,由于机翼尺寸很大,通常在图纸上只能标岀部分关键点的数据.某型号飞机的机翼上缘轮廓线的部分数据如下。用插值法求兀每改变0」时的y值,画出图形表示。x04.749.051938y05.23&111.9716.155717」7616.349514.6311413315217112.166.697.031903.9902、已知观测数据点如表所示:x00.10.20.30.40.50.60.70.80.91y・0.4471.9783.286.167.087.347.669.569.489.311.2用3次多项式函数拟
3、合这些数据点,画出图形。3、教材习题1,用griddata插值函数,三次插值。4、教材习题4,初始值为V°=0,T=0.05三.上机方法与步骤给出相应的问题分析及求解方法,并写出Matlab程序,并有上机程序显示。1、在飞机的机翼加工时,由于机翼尺寸很大,通常在图纸上只能标出部分关键点的数据.某型号飞机的机翼上缘轮廓线的部分数据如下。用插值法求兀每改变0.1吋的y值,画出图形表示。X04.749.051938577695114133152171190y05.238.111.9716.1517.116.3414.6312.166.697.033.990问题分析:它是关于一维插值,利用
4、程序interpl定义插值,method运用三次样条插值spline,再运用plot画出插值和题中所给的数据。matlab程序:x=[04.749.051938577695114133152171190];y=[05.238.111.9716.1517.116.3414.6312.166.697.033.990];xl=0:0.1:190;yl=interpl(x,y,xl,'spline');plot(x,y,'*',xl,yl,'b',x,y,'r:')xlabel('x'),ylabel('y')2、己知观测数据点如表所示:x00.10.20.30.40.50.60.70.8
5、0.91丿-0.447L9783.286.167.087.347.669.569.489.311.2用3次多项式函数拟合这些数据点,画111图形。问题分析:它是关于曲线拟合的线性拟合,利用程序polyfit定义多项式,再运用plot画11!插值和题中所给的数据。matlab程序: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,yz'k+1,x,zz'r')Editor-F:lchazhi2.m*FileEditTextGoCell
6、ToolsDebugDesktopWindowHelp:dad回b目1.0+1.1箕癖逆Q,1一x=0:0.1:1:1一y=[~0.4471.9783.286.167.087.347.669.569.489.311.2]:2-A=polyfit(x?y?3)3一z=polyval(A,x);4—plot(x,yjk+'',x,z,?r?)3、教材习题1,用griddata插值函数,三次插值。问题分析:它是关于插值基点为散乱节点,利用程序griddata定义多项式,method运用三次插值cubic,再运用plot画出插值和题屮所给的数据。matlab程序:T=[6868878710
7、6106140140];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(TzP,K,cT,cP,1cubic1)叼Editor・F:lchazhi3.mFileEditTextGoCellToolsDebugDesktopWindowHelp:♦回「回