复合梯形公式、复合辛普森公式&.doc

文档介绍：复合梯形公式、复合辛普森公式&(精选.)
复合梯形公式、复合辛普森公式&(精选.)
1 / 4
复合梯形公式、复合辛普森公式&(精选.)
宁夏师范学院数学与计算机科学学院
《数值分析》实验报告
复合梯形公式、复合辛普森公式&(精选.)
复合梯形公式、复合辛普森公式&(精选.)
1 / 4
复合梯形公式、复合辛普森公式&(精选.)
宁夏师范学院数学与计算机科学学院
《数值分析》实验报告
实验序号： 4 实验项目名称： 复合梯形公式、 复合辛普森公式
学 号 姓 名 专业、班级
实验地点 指导教师 时 间 2013年11月6日
一、实验目的及要求
掌握复合梯形公式与复合辛普森公式的基本思想。
编程实现用复合梯形公式与复合辛普森公式求积分。
熟悉 matlab 软件的使用。二、实验设备（环境）及要求
1、环境要求：
硬件：一般要求 486 以上的处理器、 16MB以上内存、足够的的硬盘可用空间 ( 随安装组
件的多少而定 ) ；
软件： MATLAB编程软件。
三、实验内容及步骤
计算积分值： I= 1 sin( x) dx .
x
复合梯形公式：
源程序：
function y=fx(x)
y=sin(x)./x;
function T_n=fht(a,b,n)
h=(b-a)/n;
for k=0:n
x(k+1)=a+k*h;
if x(k+1)==0
x(k+1)=10^(-10);
end
end
T_1=h/2*(fx(x(1))+fx(x(n+1)));
for i=2:n
F(i)=h*fx(x(i));
end
T_2=sum(F);
T_n=T_1+T_2;
复合辛普森公式：
1 / 4word.
复合梯形公式、复合辛普森公式&(精选.)
复合梯形公式、复合辛普森公式&(精选.)
4 / 4
复合梯形公式、复合辛普森公式&(精选.)
源程序：
function y=f(x)
y=sin(x)./x;
function S_n=S_P_S(a,b,n)
h=(b-a)/n;
for k=0:n
x(k+1)=a+k*h;
x_k(k+1)=x(k+1)+1/2*h;
if (x(k+1)==0)|(x_k(k+1)==0)
x(k+1)=10^(-10);
x_k(k+1)=10^(-10);
end
end
S_1=h/6*(f(x(1))+f(x(n+1)));
for i=2:n
F_1(i)=h/3*f(x(i));
end
for j=1:n
F_2(j)=2*h/3*f(x_k(j));
end
S_2=sum(F_1)+sum(F_2);
S_n=S_1+S_2;
四．实验结果与数据处理
：
T_1=fht(0,1,2) T