Question

matlab计算曲线积分

如图，一组数据在matlab中绘制成曲线如图，怎么求图中阴影各部分的积分？求各部分而不是总和。

Answer 1

图中阴影各部分的积分可以用trapz梯形数值积分函数来求解。求解时，应注意每个低峰的数值范围，可以通过for循环语句找出y<=0的x、y值，用trapz（x，y）求出各低峰的面积。例如找到某一低峰的x、y值，其积分面积

>> X =[3.1416 3.4558 3.7699 4.0841 4.3982 4.7124 5.0265 5.3407 5.6549 5.969 6.2832];

>> Y =[0 -0.30902 -0.58779 -0.80902 -0.95106 -1 -0.95106 -0.80902 -0.58779 -0.30902 0];

>> abs(trapz(X,Y))

ans =       1.9835  %计算结果