在MATLAB中，如何利用cumsum函数求积分

在MATLAB中，如何利用cumsum函数求积分

在MATLAB语言中，采用矩形求积法求解积分由来实现。其使用格式如下。

·对向量x,cumsum(x)命令返回一个向量，该向量的第N个元素是x的前N个元素的和。

·对矩阵x,cumsum(x)命令返回一个和x同型的矩阵，该矩阵的列即为对x的每一列的积累和。

·对N维数组x，cumsum(x)命令从第一个非独立数组开始操作。

·cumsum(x,DIM)命令中，参数DIM指明是从第一个非独立维开始。

例题 利用cumsum函数求积分

>> x1=[1 2 3 4 5 6 7 8 9]

x1 =

1　2　3　4　5　6　7　8　9

>> cumsum(x1)

ans =

1　3　6　10　15　21　28　36　45

>> x2=[1 2 3;4 5 6;7 8 9]

x2 =

1　2　3

4　5　6

7　8　9

>> cumsum(x2)

ans =

1　2　3

5　7　9

12　15　18

>> cumsum(x2,1)

ans =

1　2　3

5　7　9

12　15　18

>> cumsum(x2,2)

ans =

1　3　6

4　9　15

7　15　24

>> cumsum(x2,3)

ans =

1　2　3

4　5　6

7　8　9

>> t=0:0.1:10;

>> x=sin(t);

>> y=cumsum(x)*0.1;

>> plot(t,x,'r-',t,y,'k*')

运行结果如图5-3所示。

图5-3 矩形求积图像

从图5-3可以看出，求出的积分曲线和余弦曲线形状相同，这与理论计算的结果是相符合的。