数值积分方法详解：复合梯形、辛普森和龙贝格算法

数值积分方法详解：复合梯形、辛普森和龙贝格算法

本文将详细讲解三种常用的数值积分方法：复合梯形公式、复合辛普森公式和龙贝格算法，并通过代码示例展示如何使用这些方法计算函数的积分。

1. 定义函数

function y=fun41(x)
y=sqrt(x).*log(x);
end

• 该代码定义了一个名为 fun41 的函数，输入参数为 x，输出结果为 x 的平方根与自然对数的乘积。

2. 复合梯形公式和复合辛普森公式

clear;
clc;
h=0.001; % h为步长, 可分别令h=1,0.1,0.01,0.001
n=1/h;
t=0;
for i=1:n-1
    t=t+fun41(i*h);
end
T=h/2*(0+2*t+fun41(1));
T=vpa(T,10)
% 以上为复合梯形公式

% 以下为复合辛普森公式
s1=0;
s2=0;
for i=0:n-1
    s1=s1+fun41(h/2+i*h);
end
for i=1:n-1
    s2=s2+fun41(i*h);
end
S=h/6*(0+4*s1+2*s2+fun41(1));
S=vpa(S,10)

• 该代码段首先清空命令窗口和工作空间。
• 定义步长 h 为 0.001，并将 n 设为 1/h。
• 初始化 t 为 0，循环计算 fun41 的值并累加到 t 中。
• 最终根据复合梯形公式计算数值积分 T。
• 将 T 保留 10 位小数。
• 初始化 s1 和 s2 为 0，分别循环计算 fun41 的值并累加到 s1 和 s2 中。
• 最终根据复合辛普森公式计算数值积分 S。
• 将 S 保留 10 位小数。

3. 龙贝格算法

clear;
clc;
m=16;
h=1;
T(1)=(0+fun41(1))*h/2
for i=2:m
    h=h/2;
    n=1/h;
    t=0;
    for j=1:2:n-1
        t=t+fun41(j*h);
    end
    T(i)=T(i-1)/2+h*t; %梯形公式
end
for i=1:m-1
    for j=m:i+1
        T(j)=4^i/(4^i-1)*T(j)-1/(4^i-1)*T(j-1);
        % 通过不断的迭代求得T(j)，即T表的对角线上的元素。
    end
end
vpa(T(m),10)

• 该代码段首先清空命令窗口和工作空间。
• 初始化 m 为 16，h 为 1，根据龙贝格算法计算第一项数值积分 T(1)。
• 循环计算 T(i) 的值，每次步长减半，累加 fun41 的值并计算梯形公式，更新 T 表中的元素。
• 根据龙贝格算法迭代公式，不断更新 T 表中对角线上的元素，最终输出 T(m) 的值保留 10 位小数。

总结

本文介绍了三种常用的数值积分方法：复合梯形公式、复合辛普森公式和龙贝格算法。通过代码示例，可以更直观地理解这些方法的具体实现过程。对于复杂的函数，使用这些数值积分方法可以有效地求解其积分值。