《数值分析》matlab实验报告-梯形、辛普森求积公式.doc

文档介绍：东莞理工学院
《数值分析》
实验报告
实验名称：牛顿插值法
系 别：计算机学院
专 业：2013级信息与计算科学
班 级：1班
姓 名：
学 号：
实验日期：
实验内容
用不同数值方法计算积分。
算法说明
梯形求积公式算法：
将积分区间划分为n等份，步长分点为 。
积分。
辛普森求积公式算法：
其中h为步长。
Matlab软件程序清单
梯形求积公式TiXing_quad(a,b,h)：
function t = TiXing_quad(a,b,h)
%a为积分下界，b为积分上界，h为步长。
format long
x = a:h:b;
y = sqrt(x).*log(x);
y(1) = 0;
t = 0;
for k=1:(b-a)/h,
t=t+y(k)+y(k+1);
end
t=t*h/2;
辛普森求积公式Sinpson_quad(a,b,h)：
function s=Sinpson_quad(a,b,h)
format long
x=a:h:b;
y=sqrt(x).*log(x);
z=sqrt(x+h/2).*log(x+h/2);
y(1)=0;
s=0;
for k=1:(b-a)/h,
s=s+y(k)+y(k+1)+4*z(k);
end
s=s*h./6;
运行结果
真值I=-4/9=-
次数
I（梯形公式）
I（辛普森公式）
50
-
-
100
-
-
200
-
-
分析与思考