Simpson公式误差.pptVIP

博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 * 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 Simpson公式的误差 格林公式中曲线积分处理 右矩形公式应用 求积分的蒙特卡罗方法 龙贝格外推计算公式 《数值分析》 21 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 S[ xk] = I[ xk ] ( k = 0,1,2,3) ? 代数精度为3 2/16 I[1]=b – a I[x]=(b2 – a2)/2 I[x2]=(b3 – a3)/3 I[x3]=(b4 – a4)/4 S[1]=b – a S[x]=(b2 – a2)/2 S[x2]=(b3 – a3)/3 S[x3]=(b4 – a4)/4 与 I[x4] 不等 R[ f ]=I[ f ] – S[ f ] 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 x0=a, x1=(a+b)/2, x2= b. 三次多项式H(x)满足 H(x0)=f(x0), H(x1)=f(x1), H(x2)=f(x2), H’(x1)=f’(x1) 插值误差 3/16 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 令 h =(b – a)/2, x = x0+ t h ,则 令 u = t – 1, 则 Simpson公式的求积误差: 4/16 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 令h=(b – a)/(2m) , 取 xk=a+kh ( k =0,1,···,2m) 复合Simpson公式 5/16 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 格林公式的离散化——多边形面积计算 令 Q=x, P= –y 得区域 D 的面积计算公式 将? 分划, ?j : (xj，yj)?(xj+1，yj+1) ( j = 1,2,···,n) t∈( 0，1 ) 参数方程: D 6/16 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 ? 面积公式 顶点按逆时针排列,且 (xn+1，yn+1)=(x1，y1) 7/16 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 OPj = [xj , yj ], OPj+1 = [xj+1 , yj +1] P6 = P1 当OPj 到 OPj+1 为逆时针时,取正号; 当OPj 到 OPj+1 为顺时针时,取负号; 8/16 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 右矩形公式及其应用 a b Cornu螺线 ( -4.5≤ s ≤4.5 ) 记 sk = 0.1k，(k = 1，2，…，45). 9/16 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 x0=0;y0=0;h=0.1;s=h; x(1)=h*cos(.5*s*s); y(1)=h*sin(.5*s*s); for k=2:45 s=s+h; x(k)=x(k-1)+h*cos(.5*s*s); y(k)=y(k-1)+h*sin(.5*s*s); end plot([0,x],[0,y]) 10/16 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 蒙特卡罗法求积分 N=input(input N:=); f=inline(x.^3./(exp(x)-1)); t=0.1:.1:10;y0=f(t); plot(t,y0),hold on P=rand(N,2); x=5*P(:,1); y=1.5*P(:,2); n=sum(y=f(x)); q=7.5*n/N I=find(y=f(x)); plot(x(I),y(I),g.) N =2000: q= 4.8975, 4.9256, 4.7550, 4.9800 ····· 11/16 博学之，审问之，慎寺之，明辩之，笃行之。精心整理，欢迎收藏 蒙特卡罗法求体积 计算位于锥面 z2 = x2 + y2 上方与球面 x2 + y2 + (z – 1)2 = 1 内部区域的体积 锥形体可表示为 for k=1:7 P=rand(10000,3); x=2*P(:,1)-1; y=2*P(:,2)-1; z=2*P(:,3); RR=x.^2+y.^2;R=sqrt(RR); q(k)=8*sum(z=Rz=1+sqrt(1-RR))/10000; end 3.1224 3.1368 3.1680 3.1768 3.1848 3.0632 3.1192 12/16 *