MATLAB实验报告模板范例插值与数值积分含程序结果分析可直接套用.docVIP

实验插值与数值积分

【实验目的】

1．掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，对三种插值结果进行初步分析；

2．掌握用MATLAB及梯形公式、辛普森公式计算数值积分；

3．通过实例学习用插值和数值积分解决实际问题。

【实验内容】

3．用梯形公式和辛普森公式计算由下表数据给出的积分。

k

1

2

3

4

5

6

7

xk

0.3

0.5

0.7

0.9

1.1

1.3

1.5

yk

0.3895

0.6598

0.9147

1.1611

1.3971

1.6212

1.8325

已知该表数据为函数y=x+sinx/3所产生，将计算值与精确值作比较。

【解】:手工分析怎样求解这题。

【计算机求解】:怎样设计程序？流程图？变量说明？能否将某算法设计成具有形式参数的函数形式？

【程序如下】:

%原始数据

x=[0.3:0.2:1.5];

y=[0.38950.65980.91471.16111.39711.62121.8325];

z=0.5*(1.5^2-0.3^2)+1/3*(cos(0.3)-cos(1.5))%精确值

z1=trapz(x,y)%这是一个MATLAB内部的梯形公式，你能否写一个这样的函数？然后在主%程序中调用？

k=length(y);

y1=[y(2:2:k-1)];s1=sum(y1);

y2=