【计算方法数值分析】复化梯形公式、复化辛普森公式和龙贝格数值积分

1、复化梯形公式

%复化梯形公式
function t=agui_trapz(fname,d2fname,a,b,e)
%fname为被积函数，d2fname为函数fname的二阶导数，a,b分别为下界和上届，e为精度
y=abs(feval(d2fname,a:1e-5:b));
m=max(y);
h=abs(sqrt(12*e(b-a)./m));
n=ceil((b-a)/h)
h=(b-a)/n;
fa=feval(fname,a);
fb=feval(fname,b);
f=feval(fname,a+h:h:b-h+0.001*h);
t=h*(0.5*(fa+fb)+sum(f));

2、复化辛普森公式

%复化辛普森公式
function s=agui_simpson(fname,d4fname,a,b,e)
%fname为被积函数，d4fname为函数fname的四阶导数，a,b分别为下界和上届，e为精度
y=abs(feval(d4fname,a:1e-5:b));
m=max(y);
h=abs((2800*e(b-a)./m).^(0.25));
n=ceil((b-a)/h)
h=(b-a)/n;
fa=feval(fname,a);
fb=feval(fname,b);
s=fa-fb;
x=a;
for i=1:n;x=x+h/2;s=s+4*feval(fname,x);x=x+h/2;s=s+2*feval(fname,x);
end
s=s*h/6;

3、龙贝格数值积分法

%龙贝格数值积分法
function r=agui_rbg(fname,a,b)
%fname为被积函数，a,b分别为下界和上界
e=1e-6;
i=1;j=1;h=b-a;
T(i,1)=h/2*(feval(fname,a)+feval(fname,b));
T(i+1,1)=T(i,1)/2+sum(feval(fname,a+h/2:h:b-h/2+0.001*h))*h/2;
T(i+1,j+1)=4^j*T(i+1,j)/(4^j-1)-T(i,j)/(4^j-1);
while abs(T(i+1,i+1)-T(i,i))>ei=i+1;h=h/2;T(i+1,1)=T(i,1)/2+sum(feval(fname,a+h/2:h:b-h/2+0.001*h))*h/2;for j=1:iT(i+1,j+1)=4^j*T(i+1,j)/(4^j-1)-T(i,j)/(4^j-1);end
end
T
r=T(i+1,j+1);
end

结果分析

1.由表格比较可得出，无论是复化梯形公式、复化辛普森公式还是龙贝格积分公式，都有着较高的精度，相较之下，龙贝格积分公式精度最高。
2.数值求积的误差除了与被积函数有关之外，还与积分区间的长度(b-a)有关，积分区间越小，则求积公式的截断误差也越小。而在对积分区间作同样的分割的条件下，复合辛普生求积公式比复合梯形公式的计算精度高。
3.在计算速度方面，从表中可看出，复化梯形公式的等分数要比其它两个大得多，且从计算结果上可看出复化辛普生公式也比复化梯形公式的收敛速度快得多。
4.从计算公式复杂性而言，龙贝格数值积分公式的计算量最少。

【计算方法数值分析】复化梯形公式、复化辛普森公式和龙贝格数值积分相关推荐

复化辛普森公式的误差matlab,求两个matlab的程序，分别利用复化梯形公式和辛普森公式求解误差函数erf(x)。...
main.m g=@(t)exp(-t.^2); x=linspace(0,5,100); y1=zeros(1,100); y2=zeros(1,100); for i = 1:100 y1(i)= ...
Python实现梯形公式、辛普森公式、复化梯形公式、复化辛普森公式
数值分析:梯形公式 .辛普森公式.复化梯形公式.复化辛普森公式 Python实现梯形公式 .辛普森公式.复化梯形公式.复化辛普森公式数值求积公式概念梯形公式与辛普森公式梯形公式与辛普森公式的余项 ...
数值分析复化求积matlab,MATLAB数值分析实验二(复合梯形、辛普森和龙贝格求积,以及二重积分计算等)...
1.理解如何在计算机上使用数值方法计算定积近似值; 2.学会复合梯形.复合Simpson和龙贝格求积分公式的编程与应用. 3.探索二重积分在矩形区域的数值积分方法. 佛山科学技术学院实验报告 ...
数值分析复化梯形公式matlab,数值分析复化梯形公式,复化Simpson公式MATLAB程序
<数值分析复化梯形公式,复化Simpson公式MATLAB程序>由会员分享,可在线阅读,更多相关<数值分析复化梯形公式,复化Simpson公式MATLAB程序(1页珍藏版)>请 ...
辛普森复合求积公式matlab,MATLAB数值分析实验二(复合梯形、辛普森和龙贝格求积,以及二重积分计算等).doc...
[摘要]佛山科学技术学院实验报告课程名称数值分析实验项目数值积分专业班级机械工程姓名余红杰学号 2111505010 指导教师陈剑成绩日期月日一.实验目的 ...
数值计算大作业：数值积分（梯形、辛普森与龙贝格方法在Matlab实现）
作为研究生的入门课,数值计算的大作业算是所有研究生开学的重要编程作业. 针对数值积分的编程,我把梯形.辛普森与龙贝格方法在MATLAB中编写的程序放在文章最后了,需要的同学自取. PS:附录中的程序并 ...
【数值分析实验MATLAB】数值积分：梯形公式、辛普森公式、复化梯形公式、复化辛普森公式、龙贝格算法、自适应求积方法
(1)Matlab代码: f1: if x==0y=1; else%y=exp(-x^2);y=(sin(x))/x; end 梯形公式: a=0; b=1; tx=(b-a)/2*(f1(a)+f1 ...
复化梯形公式，Newton-Cotes公式，变量代换后的复化梯形公式，Gauss-Legendre公式，Gauss-Jacobi公式插值积分的精确度比较
1.问题分别计算积分 Ic=∫01cos⁡xxdx=1.809048475800...I_c=\int_0^1\frac{\cos{x}}{\sqrt{x}}dx=1.809048475800... ...
求解函数定积分，梯形公式、复化梯形公式、复合辛普森公式求解定积分近似值程序
以函数f(x)=sin(x)/x为例,求解其在[0,1]区间的定积分. <span style="font-family:SimSun;font-size:12px;"> ...

【计算方法数值分析】复化梯形公式、复化辛普森公式和龙贝格数值积分

【计算方法数值分析】复化梯形公式、复化辛普森公式和龙贝格数值积分

【计算方法数值分析】复化梯形公式、复化辛普森公式和龙贝格数值积分相关推荐

最新文章

热门文章