由于在实际应用中有些被积函数很难用公式表示,或者虽然能用公式表示,但很难求得原函数,往往需要对定积分进行近似计算。常用的近似计算方法有矩形法、梯形法和抛物线法等。使用MATLAB中的trapz函数可以用梯形法近似求取定积分的值。
例如
\[\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