数学建模培训matlab基础讲义office9版摘要.ppt

依据积分中值定理， 就是说，底为 而高为 的矩形面积恰恰等于所求曲边梯形的面积。 取 内若干个节点 处的高度 ，通过加权平均的方法生成平均高度 ，这类求积公式称机械求积公式： 式中 称为求积节点， 称为求积系数，亦称伴随节点的权。 d=pi/1000; t=0:d:3*pi; y=fun(t); re1s=cumsum(y)*d; re1=re1s(end); %显示积分值的最后两个结果 re2=trapz(y)*d; disp([re1 re2]); return; %---------- function y=fun(t) y=exp(-0.5*t).*sin(t+pi/6); return; 输出： cumsum--0.9016 trapz--0.9008 精确值：0.9008407878 梯形求积比矩形求积精度高。 note：节点数与精度有什么关系？ 3）变步长辛普生法 基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数的调用格式为： [I,n]=quad(fname