Matlab学习：复合梯形求积公式和复合辛普生求积公式

最新推荐文章于 2024-05-07 18:36:03 发布

晞微

最新推荐文章于 2024-05-07 18:36:03 发布

阅读量8.8k

点赞数 16

文章标签：学习 matlab

本文链接：https://blog.csdn.net/qq_43013173/article/details/134483213

版权

复合梯形求积公式：

function t = agui_trapz(fname,a,b,n)
% 为被积函数，a,b分别为下界和上界，n为等分数
% 复合梯形求积公式
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));
end

测试结果：

复合辛普生求积公式：

function s = agui_simpson(fname,a,b,n)
% 为被积函数，a,b分别为下界和上界，n为等分数
% 复合辛普生求积公式
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;

测试结果：