资源描述：
《matlab实现一元函数积分.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、一元函数的数值积分l函数1trapz功能复合梯形求积公式格式T=trapz(Y)%用等距复合梯形求积公式近似计算Y的积分。若Y是一向量，则trapz(Y)为Y的积分；若Y是一矩阵，则trapz(Y)为Y的每一列的积分。T=trapz(X,Y)%用复合梯形求积公式计算Y在X点上的积分，X和Y可以是负数。若X为一列向量，Y为矩阵，且size(Y,1)=length(X)，则trapz(X,Y)通过Y的第一个非单元集方向进行计算。T=trapz(…,dim)%沿着dim指定的方向对Y进行积分。若参量中包含
2、X，则应有length(X)=size(Y,dim)。例1计算的数值积分。计算可得。准备工作：>>X=0:pi/100:pi;>>Y=sin(X);函数调用：>>T=trapz(X,Y)可得>>T=1.9998或者>>T=pi/100*trapz(Y)可得>>T=1.9998例2计算复积分>>z=exp(i*pi*(0:100)/100);>>trapz(z,1./z)可得>>ans=0.0000+3.1411iproduceZ=1.9998l函数2quad、quadl、quad8功能数值定积分，自
3、适应复合Simpson求积公式。格式q=quad(fun,a,b)%近似地从a到b计算函数fun的数值积分，误差为10-6。若给fun输入向量x，应返回向量y，即fun是一单值函数。q=quad(fun,a,b,tol)%用指定的绝对误差tol代替缺省误差。tol越大，函数计算的次数越少，速度越快，但结果精度变小。q=quad(fun,a,b,tol,trace,p1,p2,⋯)%将可选参数p1,p2,…等传递给函数fun(x,p1,p2,…)，再作数值积分。若tol=[]或trace=[]，则用缺
4、省值进行计算。[q,n]=quad(fun,a,b,⋯)%同时返回函数计算的次数n[q,n]=quadl(fun,a,b,⋯)%用高精度进行计算，效率可能比quad更好。[q,n]=quad8(fun,a,b,⋯)%该命令是将废弃的命令，用quadl代替。例3计算积分方法1：首先写一个计算被积函数函数值的函数文件myfun.mfunctiony=myfun(x)y=1./(x.^3-2*x-5);然后调用函数计算数值积分>>Q=quad(@myfun,0,2)可得>>Q=-0.4605方法2：F=@
5、(x)1./(x.^3-2*x-5);%将被积函数用匿名函数来表达Q=quad(F,0,2);注：注意以上两种方法quad中被积函数的调用方式！
