复合梯形公式课件.pptVIP

复合梯形公式;椭圆周长计算:;x = a sinφ cosθ y = b sinφ sinθ z = c cos φ ;积分和式的计算: f (x)∈C[a, b]单增, 令 h = (b – a)/n, xk = a + kh, (k = 0,1,2,···,n); 4.8612e+004 4.8660e+004 4.8683e+004 4.8695e+004 4.8701e+004 4.8704e+004 4.8706e+004 4.8707e+004 4.8707e+004;数值求积公式的一般形式(机械求积公式);梯形公式: 取A0 =A1 =(b – a )/2;;插值型求积公式: 在 [a，b]上取 a≤ x0 x1 x2 …… xn≤b 作Lagrange插值;插值型求积公式的余项 ;例3取 x0 =a, x1 =0.5(a+b), x2 = b ,则 h=0.5(b – a );定义: 若一个求积公式对f(x)= xi(i=0,1,...,m)能精确成立，但对f(x)= xm+1不精确成立，则称该公式具有m次代数精度。 令机械求积公式对f(x)= xi(i=0,1,...,n)精确成立，那么得线性方程组 当节点xk (k=0,1,...,n)给定且互异时，系数Ak可由上式确定。;定理：(n+1)个节点的求积公式为插值型的充要条件是该公式至少有n次代数精度.;解: 取f(x)= 1, x, x2 令求积公式准确成立 ;取等距结点xj = a + jh时,插值型求积公式称为Newton-Cotes公式 ;复合梯形求积公式 ;;取;f=inline(sqrt((7782.5*sin(x)).^2cos(x).^2)); t=0.25*pi*(f(0)+f(pi/2)); n=1;h=pi/2;e=1;k=0; while e0.01 s=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+1 end 4*t;f=inline(sqrt((7782.5*sin(x)).^2cos(x).^2)); T=4*quad(f, 0, pi/2) Ans=4.8707e+004