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

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

1. 定义函数

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

定义函数fun41(x)，用于计算sqrt(x) * log(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，并计算n=1/h；
使用for循环计算复合梯形公式，将每个点的函数值相加；
根据复合梯形公式的公式，计算数值积分T；
使用vpa函数将T保留10位小数；
使用for循环计算复合辛普森公式，将每个点的函数值相加；
根据复合辛普森公式的公式，计算数值积分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=16，h=1，以及T(1)=h/2*(0+fun41(1))，即第一次迭代的梯形公式；
使用for循环，不断对T进行迭代，计算出T表的对角线上的元素；
使用vpa函数将最终的数值积分保留10位小数。

通过以上三种方法，可以有效地求解无法直接求解的定积分。选择哪种方法取决于积分函数的性质以及精度要求。龙贝格算法通常具有更高的精度，但计算量也更大。