复化Simpson求积公式 | matlab

最新推荐文章于 2021-11-10 17:02:22 发布

Sany 何灿

最新推荐文章于 2021-11-10 17:02:22 发布

阅读量1.1w

点赞数 8

CC 4.0 BY-SA版权

分类专栏：数值计算

本文链接：https://blog.youkuaiyun.com/SanyHo/article/details/107251949

数值计算专栏收录该内容

44 篇文章

订阅专栏

本文介绍了一种使用辛普森法则进行数值积分的方法，详细解释了如何通过将积分区间等分为n份并利用特定公式计算积分值。通过具体实例演示了如何在MATLAB中实现该算法，并给出了sin函数在[0,π]区间的积分结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

function T=simpson(f_name, a, b, n)
% f_name为要求的定函数y=f(x)所在的程序文件名
% a为积分下限
% b为积分上限
% n为积分区间[a,b]划分成小区间的等份数

h = (b-a) / n;
x = a + (0:n) * h;
f = feval(f_name, x);
N = length(f) - 1;
if N == 1
    fprintf('Data has only one interval')
    return;
end
if N == 2
    T = h / 3 * (f(1) + 4 * f(2) + f(3));
    return;
end
if N == 3
    T = 3 / 8 * h * (f(1) + 3 * f(2) + 3 * f(3) + f(4));
    return;
end
T = 0;
if 2 * floor(N / 2) == N
    T = h / 3 * (2 * f(N-2) + 2 * f(N-1) + 4 * f(N) + f(N+1));
    m = N - 3;
else
    m = N;
end
    T = T + (h/3) * (f(1) + 4 * sum(f(2:2:m)) + 2 * f(m+1));
if m > 2
    T = T + (h/3) * 2 * sum(f(3:2:m));
end

function y=simpson_f(x)
y = sin(x)
end

T = simpson('simpson_f', 0, pi, 4)

%=> y = 0    0.7071    1.0000    0.7071    0.0000
%=> T = 2.0046