matlab 矩形公式求积分

matlab 矩形公式求积分

k=[1 2 3 4 5 6 7];
xk=[0.3 0.5 0.7 0.9 1.1 1.3 1.5];
fxk=[0.3895 0.6598 0.9147 1.1611 1.3971 1.6212 1.8325];


s1=sum(fxk(1:end-1).*(diff(xk))); %矩形
s2=trapz(xk,fxk);           %梯形  


矩形公式就是把将微分的面积元看成是矩形 ΔS(i)=fxk(i)*Δx(i)
梯形公式就是把将微分的面积元看成是梯形 ΔS(i)=(fxk(i)+fxk(i+1))/2*Δx(i)




矩形





梯形





一般梯形的面积计算更为精确

trapz，矩形公式是微积分的定义式，假设每一个积分区间的值是一样的，积分区间和所成的值构成了矩形追问我知道用梯形公式s2=trapz(x,y)
但不知道怎么用“矩阵公式”追答假设积分：sin(x)，那是设置积分区间为0--pi（上下限），积分步长为pi/100
X = 0:pi/100:pi;Y = sin(x);Z = trapz(X,Y)
题解：
X=0.3:0.2:1.5;Y=[0.3895 0.6598 0.9147 1.1611 1.3971 1.6211 1.8325];trapz(X,Y)
结果：1.3730