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

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

1. 定义函数

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

定义函数 fun41(x)，输入参数为 x，输出参数为 y，y 的值为 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，并根据 h 计算出积分间隔的个数 n。然后，使用循环累加 fun41(i*h) 的值到 t 中，并使用复合梯形公式计算数值积分 T，最后使用 vpa 函数将 T 保留 10 位有效数字。

同理，代码还使用循环累加 fun41(h/2+i*h) 和 fun41(i*h) 的值到 s1 和 s2 中，并使用复合辛普森公式计算数值积分 S，最后使用 vpa 函数将 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 和初始步长 h，并使用梯形公式计算了 T(1) 的值。

然后，使用循环对 h 进行二分，并使用梯形公式计算 T(i) 的值。最后，使用双重循环对 T 数组进行迭代，通过递推公式不断更新 T(j) 的值，最终得到 T(m) 的值，即龙贝格算法计算的数值积分结果。

总结

本文详细讲解了三种常用的数值积分方法：复合梯形公式、复合辛普森公式和龙贝格算法，并用MATLAB代码示例演示了如何使用这些方法计算数值积分。希望本文能够帮助大家理解这些数值积分方法，并将其应用到实际问题中。