数值积分复华辛普森、梯形公式与误差定区间块

最新推荐文章于 2023-02-13 10:06:46 发布

鹦鹉学徒

最新推荐文章于 2023-02-13 10:06:46 发布

阅读量1.4k

点赞数

文章标签： matlab 矩阵

本文链接：https://blog.csdn.net/qq_43442189/article/details/110521415

版权

该博客介绍了如何使用MATLAB实现数值积分，包括梯形公式（T）和复华辛普森公式（FX）。通过函数JF，可以根据用户指定的积分函数、区间和精度要求计算积分。同时，提供了辅助函数Qf用于根据误差要求确定合适的区间划分数量。Qf函数利用泰勒展开和误差分析来估计满足误差要求的n值。

摘要由CSDN通过智能技术生成

%JF数值积分，fun为积分函数，例如fun=@(x)3*x，a、b分别为下上限，n为划分区间，c为所选积分方法（T代表梯形，FX代表复华辛普森）。其中n可由第二个函数通过设定误差R_n求出满足的n。

function [op]=JF(fun,a,b,n,c)
if (nargin<5)
c='T';
end
x1=(b-a)/2;
x2=(b+a)/2;
h=(b-a)./n;
Sum=0;
Sum_1=0;
if c=='T'
op=x1*fun(x2);
elseif c=='X'
op=x1/3*(fun(a)+4*fun(x2)+fun(b));
elseif c=='FT'
for i=1:n-1
x=a+i*h;
Sum=Sum+fun(x);
end
op=h/2*(fun(a)+2*Sum+fun(b));
elseif c=='FX'
for i=0:n-1
x=a+i*h;
x_1=x+1/2*h;
Sum_1=Sum_1+fun(x_1);
if i~=0
Sum=Sum+fun(x);
end
e