数学实验3.docxVIP

实验3 插值与数值积分实验目的：1．掌握用MATLAB 计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，对三种插值结果进行初步分析；2．掌握用MATLAB 及梯形公式、辛普森公式、高斯公式计算数值积分；3．通过实例学习用插值和数值积分解决实际问题。实验内容：实验预备1、编制计算拉格朗日插值的M文件，如下：function y=lagr(x0,y0,x)n=length(x0);m=length(x);for i=1:m z=x(i); s=0;for k=1:n p=1;for j=1:nif j~=k p=p*(z-x0(j))/(x0(k)-x0(j));endend s=p*y0(k)+s;end y(i)=s;end2、编制辛普森公式的.m文件，并命名为simp.m文件。如下：function s=simp(x,y)n=length(x);s=0;m=(n-1)/2;for k=1:(m-1)s=s+4*y(2*k+2)+2*y(2*k+1);ends=s+y(1)+y(n);s=s*(x(2)-x(1))/3;end实验题目【问题1】表3.7 中给出的x,y 数据位于机翼剖面的轮廓线上，y1 和y2 分别对应轮廓的上下线。假设需要得到x 坐标每改变0.1 时的y 坐标。试完成加工所