MATLAB函数trapz(x, y, n), 其中y是x的积分, 使用梯形法则逼近函数y = f(x)的积分, 并且n(可选)与维度n进行积分。
句法
Z = trapz(Y)
Z = trapz (X, Y)
Z = trapz (..., dim)
例1
使用MATLAB函数trapz(x, y)估算积分成本
并与精确值进行比较, 评估n = 5和n = 10时的百分比误差。
解
确切的值从
对于使用梯形法则的逼近, 我们让x5代表n = 5的行向量, x10代表n = 10的向量, 即∆x = 1/5和∆x = 1/10。相应的值表示为y5和y10, 曲线下的面积分别表示为area5和area10。
创建以下脚本
x5=linspace (1, 2, 5);
x10=linspace (1, 2, 10);
y5=1. /x5; y10=1. /x10;
area5=trapz (x5, y5), area10=trapz (x10, y10)
MATLAB显示以下结果:
area5 =
0.6970