原创力文档- 1、本文档共7页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
WORD格式--可编辑--专业资料
问题背景描述：（P66）表3.6给出的x，y数据位于机翼断面的轮廓线上，Y1和
Y2分别对应轮廓的上下线。假设需要得到x坐标每改变0.1时的y坐标。试玩成加工
所需数据，画出曲线，求加工断面的面积。
表3.6机翼断面轮廓线上的数据
x035791112131415
Y01.82.22.73.03.12.92.52.01.6
1
Y01.21.72.02.12.01.81.21.01.6
2
1、实验目的：
a.掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法，改变节
点的数目，对三种插值结果进行初步分析。
b.掌握用MATLAB及梯形公式、辛普森公式计算数值积分。
c.通过实例学习用插值和数值积分解决实际问题。
2、实验原理与数学模型：
插值模型
3、实验所用软件及版本：
MATLAB6.1
4、实验内容：
A、用MATLAB计算拉格朗日插值的方法：
x=[035791112131415];
y1=[01.82.22.73.03.12.92.52.01.6];
y2=[01.21.72.02.12.01.81.21.01.6];
plot(x,y1,x,y2)
xx=0:0.5:15;
yy1=interp1(x,y1,xx,spline);
yy2=interp1(x,y2,xx,spline);
plot(xx,yy1,xx,yy2)
trapz(xx,yy1)-trapz(xx,yy2)
利用数值积分可算出机翼面积=11.3053（平方米）
机翼断面曲线如下：
--完整版学习资料分享
WORD格式--可编辑--专业资料
B、用MATLAB计算分段线性插值的方法：
x0=[0,3:2:11,12:15];
y1=[01.82.22.73.03.12.92.52.01.6];
y2=[01.21.72.02.12.01.81.21.01.6];
x=0:0.1:15;
y1_in=interp1(x0,Y1,x);
y2_in=interp1(x0,Y2,x);
[x,y1_in,y1_sp,y2_in,y2_sp]
subplot(2,1,1),plot(x,y1_in,x,y2_in,b),title(interp)
trapz(x,y1_in)-trapz(x,y2_in)
利用数值积分可算出机翼面积S=10.7500（平方米）
机翼断面曲线如下机翼断面曲线如下：
--完整版学习资料分享
WORD格式--可编辑--专业资料
C、用MATLAB计算三次样条插值的方法：
x0=[0,3:2:11,12:15];
y1=[01.82.22.73.03.12.92.52.01.6];
y2=[01.21.72.02.12.01.81.21.01.6];
x=0:0.1:15;
y1_sp=spline(x0,Y1,x);
y2_sp=spline(x0,Y2,x);
[x,y1_in,y1_sp,y2_in,y2_sp]
subplot(2,1,2),plot(x,y1_sp,x,y2_sp,b),title(spline)
trapz(x,y1_sp)-trapz(x,y2
文档评论（0）