下载此文档

编写Matlab程序实现复合梯形和龙贝格数值积分.doc


文档分类:高等教育 | 页数:约4页 举报非法文档有奖
1/4
下载提示
  • 1.该资料是网友上传的,本站提供全文预览,预览什么样,下载就什么样。
  • 2.下载该文档所得收入归上传者、原创者。
  • 3.下载的文档,不会出现我们的网址水印。
1/4 下载此文档
文档列表 文档介绍
《数值分析》课程实验报告
所属学院
电气与信息工程学院
专业
学号
实验名称
编写Matlab程序实现复合梯形和龙贝格数值积分
姓名
一、实验目的
;
;

二、实验内容
1. 用Matlab软件编程实现复合梯形数值积分

三、实验步骤
1. 用Matlab软件编程实现复合梯形数值积分
(1)开启软件平台MATLAB,编程并建立M文件
(2)执行M文件trapz(f,a,b,n),观察结果:
2. 用Matlab软件编程实现龙贝格数值积分
编程及结果:
function romberg1(a,b,eps,j)
t=zeros(2,4); t(2,1)=(b-a)*(f0(a)+f0(b))/2
for k=1:j
t(1,:)=t(2,:); s=0;h=(b-a)/(2^k);
for i=1:2^(k-1)
s=s+f0(a+(2*i-1)*h)
end;
t(2,1)=t(1,1)/2+s*h;
for m=1:min(3,k)
t(2,m+1)=(4^m*t(2,m)-t(1,m))/(4^m-1);
end;
k, t
if (abs(t(2,4)-t(1,4))<eps) & (k>3)
break
end;
end
t(2,4)
function f =f0(x)
f=sin(x)/x;
>> a=,b=1,eps=,j=3
a =
-011
b =
1
eps =
-007
j =
3
>> romberg1(a,b,eps,j)
t =
0 0 0 0
0 0 0
s =

k =1
t =
0 0 0
0 0
s =

编写Matlab程序实现复合梯形和龙贝格数值积分 来自淘豆网www.taodocs.com转载请标明出处.

非法内容举报中心
文档信息
  • 页数4
  • 收藏数0 收藏
  • 顶次数0
  • 上传人ocxuty74
  • 文件大小49 KB
  • 时间2018-05-27