复合梯形公式matlab
1、可以利用matlab的trapz函数命令x=0:0.00001:1;%x用来储存积分点y=(x+1).*sin(x);%y用来求解积分点x处的函数值i=trapz(x,y)i?=0.7608663730793?验证该问题的解析解syms?xy=(x+1)*sin(x);%被积函数表达式ii=int(y,0,1)ii?=sin(1)?-?2*cos(1)?+?1?%ii即为该被积函数的解析解ii_e=eval(ii)?ii_e?=0.760866373071617?%ii的数值解%可以看出梯形求积公式在步长等于0.00001的情况下。
2、数值积分的解与解析解的数值能达到小数点后11位保持一致。
![](http://yyk.iask.sina.com.cn/pic/fimg/160992423928368984881.jpg)