e"=Ge""=∑cGv=∑c2v =1|cv+ C2V2+…+ nn 对于相当大的n,有 m)≈~C11 (n+1) ≈n+1 (n+1) 则 入1|=p(G) P椭圆型差分方程的送代送
椭圆型差分方程的迭代法 对于相当大的n ，有 则 + + = + = = = = = − − n n n n n n N i i n i i N i i n i i n n c v c v c v e G e c G v c v 1 2 2 1 2 1 1 1 1 1 ( ) ( 1) 1 1 1 1 1 ( 1) 1 1 1 ( ) e c v e c v n n n+ n+ ( ) ( ) 1 ( 1) G e e n i n i = +
(n+1 entp (n+p-1) 2,(n+P-2)≈.≈/1e (n+p) =(p(G) 要求从第n步迭代到第(n+p)步的误差减少为e(n) 的10-q,即要求 =(p(G)2≤109 P椭圆型差分方程的送代送
椭圆型差分方程的迭代法 要求从第 n 步迭代到第（n+p）步的误差减少为 e（n） 的10–q ，即要求 ( ) 1 2 ( 2) 1 ( 1) 1 (n p) n p n p p n e e e e + + − + − ( ) 1 (n 1) n e e + p p n i n p i G e e ( ( )) ( ) 1 ( ) = + p p q G − 1 = (( )) 10
两边去对数得 lg P(G)(g(G)<0) 或 ≥h10 10 -hP(G (p( p(G)越小,-hnp(G)越大,称 R(G=-hnp(G) 为迭代法的渐近收敛速度。 P椭圆型差分方程的送代送
椭圆型差分方程的迭代法 两边去对数得 或 越小， 越大，称 为迭代法的渐近收敛速度。 plg (G) −q (lg( ) 0) lg ( ) − G G q p ln ( ) ln 10 G q p − 10–q 0 x p ((G)) y (G) − ln (G) R(G) = −ln (G)
二、 Jacobi迭代和 Gauss-Seidel迭代 方程组Ax=b A=D--R n D A 12 R n-In 0 P椭圆型差分方程的送代送
椭圆型差分方程的迭代法 二、Jacobi 迭代和 Gauss-Seidel 迭代 方程组 A = D – L – R Ax = b = n nn n a a a a A 1 11 1 = n n a a a D 22 11 = − − 0 0 1 1 21 n n n a a a L = − − 0 0 1 1 2 1 n n n a a a R
方程组写成(D-L-R)x=b Dx=(l+r)x+b 若an≠0(i=1,2,…,N), Jacobi迭代格式为 D(m+1) =(L+R)x)+b 或x)=D(L+R)x)+Db 即、n+1) ∑a1xm+(=12,…,N J≠l P椭圆型差分方程的送代送
椭圆型差分方程的迭代法 方程组写成 若 ，Jacobi 迭代格式为 或 即 (D − L − R)x = b Dx = (L + R)x +b a 0 (i 1,2, ,N) ii = x D L R x D b (n 1) 1 (n) 1 ( ) + − − = + + Dx L R x b n n = + + ( +1) ( ) ( ) ( 1,2, , ) 1 1 ( 1) ( ) i N a b a x a x i i i N j i j n i j j i i n i = − + = = +