梯形求积

用trapz函数可以求梯形数值积分。

【例21】

例如， $\int_{0}^{\pi} {\operatorname{sin}{(}x)dx}$ 的精确值为2。下面用trapz函数在均匀间隔的网格上求该积分的数值近似。

code.matlab

>> X=0:pi/100:pi;
>> Y=sin(X);
>> Z=trapz(X,Y)
Z=
1.9998

对于非均匀间隔的网格，则命令行为

code.matlab

>> X=sort(rand(1,101)*pi);
>> Y=sin(X);
>> Z=trapz(X,Y)
Z=
    1.9981