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

本文将详细介绍三种常用的数值积分方法，并提供相应的MATLAB代码示例。

1. 定义函数

首先，定义一个函数 fun41，输入参数为 x，输出为 y，其中 y = sqrt(x) * log(x)。

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

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 和 n = 1/h，将区间 [0, 1] 分成 n 个小区间。
对于复合梯形公式，使用循环计算 t 的值，最后代入公式计算出 T 的值。
对于复合辛普森公式，同样使用循环计算 s1 和 s2 的值，最后代入公式计算出 S 的值。
vpa 函数用于保留 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)

代码解释：

首先，设置初始步长 h = 1，计算第一次迭代的 T(1) 值。
接下来使用循环，每次将步长缩小一半，利用梯形公式计算出新的 T 值。
然后通过迭代求得 T 表的对角线上的元素。
最终输出 T 表的最后一个元素的值，即 vpa(T(m), 10)。

总结

本文详细讲解了三种常用的数值积分方法：复合梯形公式、复合辛普森公式和龙贝格算法，并提供了相应的MATLAB代码示例。希望本文能够帮助读者更好地理解和运用这些方法。