复合梯形公式.ppt

《复合梯形公式.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、从积分和式到求积公式插值型求积公式求积公式的代数精度复合梯形公式MATLAB求积分命令第五章数值积分与数值微分椭圆周长计算:x=acosty=bsint0≤t≤2椭圆积分x=asinφcosθy=bsinφsinθz=ccosφD={(θ，φ)

2、0≤θ≤2π,0≤φ≤π}思考题:椭球面的面积计算椭球面积的积分表达式？对二重积分的计算问题？三维体积的离散数据计算？积分和式的计算:f(x)∈C[a,b]单增,令h=(b–a)/n,xk=a+kh,(k=0,1,2,···,n)近似计算abf(x)4.8612e+0044.8660e+0044.8683e+0044

3、.8695e+0044.8701e+0044.8704e+0044.8706e+0044.8707e+0044.8707e+0044.8803e+0044.8755e+0044.8731e+0044.8719e+0044.8713e+0044.8710e+0044.8709e+0044.8708e+0044.8708e+0044.8708e+004a=7782.5c=972.5P.170人造卫星的轨道长度计算数值求积公式的一般形式(机械求积公式)R[f]为数值求积公式余项,x0,x1,···,xn为求积结点;A0,A1,···,An为求积系数.矩形公式:取A0

4、=(b–a)特别地，时分别称为左矩公式，中矩公式，右矩公式。梯形公式:取A0=A1=(b–a)/2abSimpson公式取A0=A3=(b–a)/6，A1=2(b–a)/3，插值型求积公式:在[a，b]上取a≤x0

5、xm+1不精确成立，则称该公式具有m次代数精度。令机械求积公式对f(x)=xi(i=0,1,...,n)精确成立，那么得线性方程组当节点xk(k=0,1,...,n)给定且互异时，系数Ak可由上式确定。定理：(n+1)个节点的求积公式为插值型的充要条件是该公式至少有n次代数精度.例4确定公式使代数精度尽可能高.类似有:Simpson公式具有3阶代数精度梯形公式代数精度为1例.矩形公式代数精度为0解:取f(x)=1,x,x2令求积公式准确成立A1=0求积公式具有至少2阶代数精度容易验证,对f(x)=x3求积公式式不能准确成立.因此这一公式只具有2次代数精度取等距

6、结点xj=a+jh时,插值型求积公式称为Newton-Cotes公式定理:当n为偶数时,n阶Newton-Cotes公式至少有(n+1)阶代数精确度。Newton-Cotes公式代数精度至少为n复合梯形求积公式将积分区间[a,b]n等分.令h=(b-a)/n.xj=a+jh取递推，得给定允许误差界ε>0,当时，结束计算并以T2n作为定积分的近似值.f=inline('sqrt((7782.5*sin(x)).^2+59621550*cos(x).^2)');t=0.25*pi*(f(0)+f(pi/2));n=1;h=pi/2;e=1;k=0;whilee>0

7、.01s=0.5*(t+h*sum(f(.5*h:h:pi/2)));e=abs(s-t);t=s;n=2*n;h=h/2;k=k+1end4*tans=4.8707e+004,(循环次数k=2)复合梯形公式计算将程序第二行改为t=0.5*pi*f(pi/2)用积分部分和式计算,循环次数k=13f=inline('sqrt((7782.5*sin(x)).^2+59621550*cos(x).^2)');T=4*quad(f,0,pi/2)Ans=4.8707e+004MATLAB求定积分命令quad(‘fun’,a,b)高阶求积分命令q=quad8('fun

8、',a,b)重积分计算命令dblquad('fun',inmin,inmax,outmin,outmax)此课件下载可自行编辑修改，此课件供参考！部分内容来源于网络，如有侵权请与我联系删除！