К.Ю. Богачёв - Практикум на ЭВМ. Методы приближения функций (1133845), страница 19
Текст из файла (страница 19)
+,1+6B1H xi > xk :95 i > k 5 yj > yl :95 j > l , -+ 0,7 9.23767//D79.2/+,-5 f (x1Y : : : Y xi y1Y : : : Y yj ) 0 (6) 0D45,6;F-,; -+6B1+ 47972 2/.47/5; @H/1<55 f 0 -+41.A (xi yj ) 2 .95 m = n 0D9.?7/57 (6) >+?/+ 2.:5,.-B 0 0537Pf (x y) =Xi+j n+1x 2.3.3.f (x1 Y : : : Y xi y1Y : : : Y yj );1iY;1 jYk=1 l=1(x ; xk )(y ; yl )-$>!@ !& 7;.,,>+-95> +:5,.//HF 0DE7 2.3.4H 65/78/+8 5/-79:+6;<55 0 4.,-/+> ,6H4.7 l = 1 0 (2), -.7. :95C65?7/57 65/78/+8 @H/1<578.-$>!@ $ )!"C+2/.45> 079E5/D -97H=+6B/51. 47972 (x1 y1), (x2 y2), (x3 y3).
.>-97CH7-,; :+,-9+5-B :95C65?.FGHF @H/1<5; 053.Pf (x y) = x + y + H3+067-0+9;FGHF H,6+05;>Pf (xk yk ) = f (xk yk ) k = 1 2 3:(7)+I@@5<57/-D , , >+=H- CD-B :+6H47/D 52 H9.0/7/5; :6+,1+,-5, :9+A+3;G78 47972 -95 -+415 (xk yk f (xk yk )), k = 1 2 3. 3/.1+, 36; 3.6B/78E7=+526+?7/5; /.> H3+C/77 007,-5 C.25,/D7 @H/1<55 2.3.45 65/78/+8 5/-79:+6;<55 5 :+,-9+5-B Pf 1.1 5A 65/78/HF 1+>C5/.<5F.C+2/.45> x ; u y ; v L(u1 v1 ) (u2 v2) (x y) = u ; u1 v ; v1 = (x ; u1)(v2 ; v1) ; (y ; v1 )(u2 ; u1) (8)2121{ 65/78/.; @H/1<5;, +C9.G.FG.;,; 0 /H6B /. :9;>+8, :9+A+3;G78 47972 -+415(u1 v1) 5 (u2 v2). 0737> -.1?7 @H/1<551 (x y) = L(x2 y2) (x3 y3) (x y)2 (x y) = L(x1 y1) (x3 y3) (x y)3 (x y) = L(x1 y1) (x2 y2) (x y)-.7. x ; x y ; y 1 (x y) = x ; x2 y ; y2 = (x ; x2 )(y3 ; y2) ; (y ; y2)(x3 ; x2) x3 ; x 2 y3 ; y 2 2 (x y) = x ; x1 y ; y1 = (x ; x1 )(y3 ; y1) ; (y ; y1)(x3 ; x1) x3 ; x 1 y3 ; y 1 3 (x y) = x ; x1 y ; y1 = (x ; x1 )(y2 ; y1) ; (y ; y1)(x2 ; x1)2121..
!"# $B $Bx2. '&+-.- *3 C+'. -.1?7116'i(x y) = i((xx yy)) i = 1 2 3:(9)i i i>77>:'i(xj yj ) = ij i j = 1 2 3:"H/1<55 'i(x y), i = 1 2 3 0+2B>7> 0 1.47,-07 C.25,/DA 0 2.3.47 65/78/+85/-79:+6;<55, -.7. :95C65?.FGHF @H/1<5F Pf CH37> 5,1.-B 0 05373XPf (x y) = i'i(x y)i=1=37 1+I@@5<57/-D i +:97376;F-,; 52 H,6+058 (7). .1 1.1 >.-95<. (I.1.3)2.3.45 65/78/+8 5/-79:+6;<55 { 735/54/.;, -+ i = f (xi yi), i = 1 2 3 5 :95C65?.FG.; @H/1<5; 5>77- 053Pf (x y) =3Xi=1f (xi yi)'i(x y):/.47/5; :95C65?.FG78 @H/1<55 Pf /.
,-+9+/7 -97H=+6B/51. (/.:95>79,S(x1 y1) (x2 y2)]) +:97376;F-,; -+6B1+ 2/.47/5;>5 @H/1<55 f 0 -+41.A, 67?.G5A /. I-+8 ,-+9+/7 (0 /.E7> :95>797 { (x1 y1) 5 (x2 y2)). 673+0.-76B/+,7,65 -97H=+6B/515 1 5 2 5>7F- +CGHF ,-+9+/H, -+ 2/.47/5; :95C65?.FG5A@H/1<58, :+,-9+7//DA 0 1.?3+> 52 /5A, ,+0:.3.F- 0 -+41.A I-+8 +CG78 ,-+9+/D. J-+ :+20+6;7- ,-9+5-B 52 :95C65?.FG5A @H/1<58 Pk f , +:973767//DA0 -97H=+6B/51.A k , 1H,+4/+-65/78/HF .::9+1,5>.<5F Pf , +:973767//HF /.N-95./=H6;<55 S k , :+ :9.056H: Pf (x y) = Pk f (x y) :95 (x y) 2 k .k=1-$>!@ $) !;7+,-9+7//HF 0DE7 :95C65?.FGHF @H/1<5; Pf , +:973767//HF /. -95./=HN6;<55 D0 = S k , >+?/+ :+6H45-B -.1?7 ,673HFG5> ,:+,+C+>.k=15$#!.
C+2/.45> 47972 V1 (D 0 ) 65/78/+7 :9+,-9./,-0+ /7:979D0/DA0 D0 @H/1<58, ;06;FG5A,; 65/78/D>5 0 1.?3+> -97H=+6B/517 i D0 =SN .kk=10737> C.25, :9+,-9./,-0. V1 (D0):(1)0'(1)i 2 V1 (D ) 'i (xj yj ) = ij i j = 1 2 : : : M1 (10)=37 f(xj yj )gj=1 2 ::: M1 { >/+?7,-0+ 079E5/ -97H=+6B/51+0 k D0 . "H/1<55(1)'(1)i 4.,-+ /.2D0.F- >. /.65-547,15 @H/1<5F 'i >+?/+2.3.-B ,673HFG5> +C9.2+>. H,-B k , l = 1 2 : : : ni { -97H=+6B/515, 5>7FG57 -+41H (xi yi) 0 1.47,-07 079E5/D. C+2/.45> 47972 'k (x y) 65/78/HFll..
!"# $B $Bx2. '&+-.- *3 C+'117@H/1<5F, +:973767//HF 0 -97H=+6B/517 k , 9.0/HF 1 0 -+417 (xi yi) 5 0 0+,-.6B/DA 079E5/.A k . .157 @H/1<55 >D :+,-9+565 9./77, ,>. (9). +=3.ll8>< 'k (x y) 7,65 (x y) 2 k(1)'i (x y) = >n7,65 (x y) 2= S k: 0l=1llil,;1.; g 2 V1(D0) 9.,16.3D0.7-,; :+ C.25,H 52 @H/1<58 f'(1)i (x y )gi=1 2 ::: M1,673HFG5> +C9.2+>:g(x y) =M1Xi=1g(xi yi)'(1)i (x y ):78,-05-76B/+, 0 ,56H (10) 5>77> 36; gb(x y) = P g(xi yi)'(1)i (x y ): gb(xj yj ) =i=1g(xj yj ), j = 1 2 : : : M1 .
673+0.-76B/+, gb(x y) ;06;7-,; 0 1.?3+> -97H=+6B/517 k 65/78/+8 5/-79:+659HFG78 @H/1<578. ,56H 1+9971-/+,-5 2.3.4565/78/+8 5/-79:+6;<55 5/-79:+6;<5+//D8 >/+=+467/ 735/,-07//7/. 673+0.N-76B/+, g = gb 0 1.?3+> -97H=+6B/517 k , 5, 2/.45-, 0 D0 = S k .k=1.,,>+-95> 2.3.4H /.A+?37/5; :95C65?.FG78 @H/1<55 Pf 2 V1(D0), ,+0:.3.FG78 , @H/1<578 f 0 -+41.A f(xj yj )gj=1 2 ::: M1 :M1Pf (xj yj ) = f (xj yj ) j = 1 2 : : : M1 :(11)"H/1<5F Pf CH37> 5,1.-B 0 0537Pf (x y) =M1Xi=1i'(1)i (x y )=37 1+I@@5<57/-D i +:97376;F-,; 52 H,6+058 (11). ,56H (10) >.-95<. (I.1.3)2.3.45 65/78/+8 5/-79:+6;<55 { 735/54/.;.
+I-+>H i = f (xi yi), i =1 2 : : : M1 5 :95C65?.FG.; @H/1<5; 5>77- 053Pf (x y) =x 2.3.4.M1Xi=1f (xi yi)'(1)i (x y ):-$>!@ %#$ 7;.,,>+-95> +:5,.//HF 0DE7 2.3.4H 65/78/+8 5/-79:+6;<55 0 4.,-/+> ,6H4.7 l = 2 0 (2), -.7. :95C65?7/57 10.39.-547,1+8 @H/1<578... !"# $B $Bx2. '&+-.- *3 C+'118-$>!@ $ )!"C+2/.45> 079E5/D -97H=+6B/51. 47972 (x1 y1), (x2 y2), (x3 y3), ,79735/D 7=+ ,-+9+/ { 47972(x4 y4) (x1 2 y1 2) (x5 y5) (x2 3 y2 3) (x6 y6) (x1 3 y1 3)=37 xk l = (xk + xl )=2, yk l = (yk + yl )=2.
.> -97CH7-,; :+,-9+5-B :95C65?.FGHF@H/1<5F 053. (2) H3+067-0+9;FGHF H,6+05;>Pf (xk yk ) = f (xk yk ) k = 1 2 : : : 6:(12)+,-9+5> C.25,/D7 @H/1<55 2.3.45 65/78/+8 5/-79:+6;<55.0737> @H/1<551 (x y) = L(x2 y2) (x3 y3) (x y)2 (x y) = L(x1 y1) (x3 y3)(x y)3 (x y) = L(x1 y1) (x2 y2) (x y)4 (x y) = L(x1 2 y1 2) (x1 3 y1 3) (x y)5 (x y) = L(x1 2 y1 2 ) (x2 3 y2 3) (x y) 6 (x y) = L(x1 3 y1 3) (x2 3 y2 3) (x y)(=37 @H/1<55 L +:973767/D 0 (8)), .
-.1?7i(x y) = i((xx yy)) i = 1 2 : : : 6:i i i+6+?5>'1(x y) = 1 (x y)4(x y) '2(x y) = 2 (x y)5(x y)'3(x y) = 3 (x y)6(x y) '4(x y) = 1 (x y)2(x y)'5(x y) = 2 (x y)3(x y) '6(x y) = 1 (x y)3(x y):(13)>77>:'i (xj yj ) = ij i j = 1 2 : : : 6:"H/1<55 'i(x y), i = 1 2 : : : 6 0+2B>7> 0 1.47,-07 C.25,/DA 0 2.3.47 65/78/+85/-79:+6;<55, -.7. :95C65?.FGHF @H/1<5F Pf CH37> 5,1.-B 0 05376XPf (x y) = i'i(x y)i=1=37 1+I@@5<57/-D i +:97376;F-,; 52 H,6+058 (12). .1 1.1 >.-95<. (I.1.3)2.3.45 65/78/+8 5/-79:+6;<55 { 735/54/.;, -+ i = f (xi yi), i = 1 2 : : : 6 5:95C65?.FG.; @H/1<5; 5>77- 053Pf (x y) =..6Xi=1f (xi yi)'i(x y): !"# $B $Bx2.
'&+-.- *3 C+'119/.47/5; :95C65?.FG78 @H/1<55 Pf /. ,-+9+/7 -97H=+6B/51. (/.:95>79,S(x1 y1) (x2 y2)]) +:97376;F-,; -+6B1+ 2/.47/5;>5 @H/1<55 f 0 -+41.A, 67?.G5A /. I-+8 ,-+9+/7 (0 /.E7> :95>797 { (x1 y1), (x2 y2) 5 (x1 2 y1 2)).
673+0.-76B/+, 7,65 -97H=+6B/515 1 5 2 5>7F- +CGHF ,-+9+/H, -+ 2/.47/5; :95C65?.FG5A @H/1<58, :+,-9+7//DA 0 1.?3+> 52 /5A, ,+0:.3.F- 0 -+41.A I-+8+CG78 ,-+9+/D. J-+ :+20+6;7- ,-9+5-B 52 :95C65?.FG5A @H/1<58 Pk f , +:973767//DA 0 -97H=+6B/51.A k , 1H,+4/+-10.39.-547,1HF .::9+1,5>.<5F Pf ,N+:973767//HF /. -95./=H6;<55 S k , :+ :9.056H: Pf (x y) = Pk f (x y) :95k=1(x y) 2 k .-$>!@ $) !;7+,-9+7//HF 0DE7 :95C65?.FGHF @H/1<5; Pf , +:973767//HF /. -95./=HN6;<55 D0 = S k , >+?/+ :+6H45-B -.1?7 ,673HFG5> ,:+,+C+>.k=15$#!.
C+2/.45> 47972 V2 (D 0 ) 65/78/+7 :9+,-9./,-0+ /7:979D0/DA0 D0 @H/1<58, ;06;FG5A,; 10.39.-547,15>5 0 1.?3+> -97H=+6B/517 i D0 =SN .kk=10737> C.25, :9+,-9./,-0. V2 (D0):(2)0'(2)i 2 V2 (D ) 'i (xj yj ) = ij i j = 1 2 : : : M2 (14)=37 f(xj yj )gj=1 2 ::: M2 { >/+?7,-0+ 079E5/ 5 ,973/5A -+471 ,-+9+/ -97H=+6B/51+0 k D0 ./.65-547,15 @H/1<5F '(2)i >+?/+ 2.3.-B ,673HFG5> +C9.2+>. H,-B k ,l = 1 2 : : : ni { -97H=+6B/515, -.157, 4-+ (xi yi) 2 k .
C+2/.45> 47972 'k (x y)10.39.-547,1HF @H/1<5F, +:973767//HF 0 -97H=+6B/517 k , 9.0/HF 1 0 -+417(xi yi) 5 0 0 +,-.6B/DA 079E5/.A 5 ,79735/.A ,-+9+/ k . .157 @H/1<55 >D:+,-9+565 9./77, ,>. (13). +=3.lllll8>< 'k (x y) 7,65 (x y) 2 k(2)'i (x y) = >n7,65 (x y) 2= S k: 0llil=1l,;1.; g 2 V2(D0) 9.,16.3D0.7-,; :+ C.25,H 52 @H/1<58 f'(2)i (x y )gi=1 2 ::: M2,673HFG5> +C9.2+>:M2Xg(x y) = g(xi yi)'(2)i (x y ):i=178,-05-76B/+, 0 ,56H (14) 5>77> 36; gb(x y) = P g(xi yi)'(2)i (x y ): gb(xj yj ) =i=1g(xj yj ), j = 1 2 : : : M2 . 673+0.-76B/+, gb(x y) ;06;7-,; 0 1.?3+> -97H=+6B/517 k 10.39.-547,1+8 5/-79:+659HFG78 @H/1<578.
,56H 1+9971-/+,-5 2.3.45M2.. !"# $B $Bx2. '&+-.- *3 C+'12065/78/+8 5/-79:+6;<55 5/-79:+6;<5+//D8 >/+=+467/ 735/,-07//7/. 673+0.N-76B/+, g = gb 0 1.?3+> -97H=+6B/517 k , 5, 2/.45-, 0 D0 = S k .k=1.,,>+-95> 2.3.4H /.A+?37/5; :95C65?.FG78 @H/1<55 Pf 2 V2(D0), ,+0:.3.FG78 , @H/1<578 f 0 -+41.A f(xj yj )gj=1 2 ::: M2 :Pf (xj yj ) = f (xj yj ) j = 1 2 : : : M2 :(15)"H/1<5F Pf CH37> 5,1.-B 0 0537Pf (x y) =M2Xi=1i'(2)i (x y )=37 1+I@@5<57/-D i +:97376;F-,; 52 H,6+058 (15). ,56H (14) >.-95<. (I.1.3)2.3.45 65/78/+8 5/-79:+6;<55 { 735/54/.;.
+I-+>H i = f (xi yi), i =1 2 : : : M2 5 :95C65?.FG.; @H/1<5; 5>77- 053Pf (x y) =x 2.3.5.M2Xi=1f (xi yi)'(2)i (x y ):-$>!@ > 7;.,,>+-95> +:5,.//HF 0DE7 2.3.4H 65/78/+8 5/-79:+6;<55 0 4.,-/+> ,6H4.7 l = 3 0 (2), -.7. :95C65?7/57 1HC547,1+8 @H/1<578.-$>!@ $ )!"C+2/.45> 079E5/D -97H=+6B/51. 47972 (x1 y1), (x2 y2), (x3 y3), -+415,376;G57 7=+ ,-9+/D /. -95 9.0/D7 4.,-5, { 47972(x4 y4) (x1 2 1 y1 2 1) (x5 y5) (x1 2 2 y1 2 2) (x6 y6) (x2 3 1 y2 3 1)(x7 y7) (x2 3 2 y2 3 2) (x8 y8) (x1 3 1 y1 3 1) (x9 y9) (x1 3 2 y1 3 2)=37 xk l s = ((3 ; s)xk + sxl )=3, yk l s = ((3 ; s)yk + syl )=3, <7/-9 -;?7,-5 -97H=+6B/51.
{ 47972(x10 y10) (x1 2 3 y1 2 3)=37 x1 2 3 = (x1 + x2 + x3 )=3, y1 2 3 = (y1 + y2 + y3)=3. .> -97CH7-,; :+,-9+5-B:95C65?.FGHF @H/1<5F 053. (2) H3+067-0+9;FGHF H,6+05;>Pf (xk yk) = f (xk yk ) k = 1 2 : : : 10:(16)+,-9+5> C.25,/D7 @H/1<55 2.3.45 65/78/+8 5/-79:+6;<55... !"# $B $Bx2. '&+-.- *3 C+'1210737> @H/1<551 (x y) = L(x2 y2) (x3 y3 )(x y)2 (x y) = L(x1 y1) (x3 y3)(x y)3 (x y) = L(x1 y1) (x2 y2 )(x y)4 (x y) = L(x1 2 1 y1 2 1 ) (x1 3 1 y1 3 1) (x y)5 (x y) = L(x1 2 2 y1 2 2 ) (x1 3 2 y1 3 2) (x y) 6 (x y) = L(x1 2 2 y1 2 2 ) (x2 3 1 y2 3 1) (x y)7 (x y) = L(x1 2 1 y1 2 1 ) (x2 3 2 y2 3 2) (x y) 8 (x y) = L(x1 3 2 y1 3 2 ) (x2 3 2 y2 3 2) (x y)9 (x y) = L(x1 3 1 y1 3 1 ) (x2 3 1 y2 3 1) (x y)(=37 @H/1<55 L +:973767/D 0 (8)), .
-.1?7i(x y) = i((xx yy)) i = 1 2 : : : 10:i i i+6+?5>'1(x y) = 1(x y)4(x y)5(x y) '2 (x y) = 2(x y)6(x y)7(x y)'3(x y) = 3(x y)8(x y)9(x y) '4 (x y) = 1(x y)2(x y)5(x y)'5(x y) = 1(x y)2(x y)7(x y) '6 (x y) = 2(x y)3(x y)7(x y)(17)'7(x y) = 2(x y)3(x y)9(x y) '8 (x y) = 1(x y)3(x y)5(x y)'9(x y) = 1(x y)3(x y)9(x y) '10 (x y) = 1(x y)2(x y)3(x y):>77>:'i(xj yj ) = ij i j = 1 2 : : : 10:"H/1<55 'i (x y), i = 1 2 : : : 10 0+2B>7> 0 1.47,-07 C.25,/DA 0 2.3.47 65/78/+85/-79:+6;<55, -.7. :95C65?.FGHF @H/1<5F Pf CH37> 5,1.-B 0 0537Pf (x y) =10Xi=1i'i(x y)=37 1+I@@5<57/-D i +:97376;F-,; 52 H,6+058 (16).