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

本文将介绍三种常用的数值积分方法：复合梯形公式、复合辛普森公式和龙贝格算法，并用MATLAB代码示例进行说明。

1. 定义函数

首先定义一个函数 fun41(x)，其输入参数为 x，输出为 y=sqrt(x).*log(x)，即 y 等于 x 的平方根乘上 x 的自然对数。

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

2. 使用复合梯形公式和复合辛普森公式计算数值积分

(1) 复合梯形公式

首先清除之前的变量和命令行窗口，将步长 h 设为 0.001，n 为步数，t 为累加器。通过 for 循环计算 t 的值，最后根据复合梯形公式计算数值积分 T，并将其精度保留到小数点后 10 位。

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)

(2) 复合辛普森公式

首先将 s1 和 s2 的初始值设为 0，分别通过 for 循环计算它们的值。最后根据复合辛普森公式计算数值积分 S，并将其精度保留到小数点后 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)

3. 使用龙贝格算法计算数值积分

首先清除之前的变量和命令行窗口，将 m 设为 16，h 设为 1，计算 T(1) 的值。接下来通过 for 循环不断迭代计算 T(i) 的值，直到计算出 T(m)。最后通过双重 for 循环，不断迭代求出 T 表的对角线上的元素，即 T(m) 的值，并将其精度保留到小数点后 10 位。

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)

本文介绍了三种常用的数值积分方法：复合梯形公式、复合辛普森公式和龙贝格算法，并用MATLAB代码示例进行说明，帮助读者理解其原理和应用。希望本文能够帮助读者学习和掌握数值积分方法，并将其应用于实际问题中。