К.Ю. Богачёв - Практикум на ЭВМ. Методы приближения функций (1133845), страница 12
Текст из файла (страница 12)
!"# # $#x19. '&+-.- 0x 19.8.3.*-0* 70.-\%&" )$&76!%;J-+- >7-+3 5,:+6B2H7-,;, 1+=3. /51.1+8 3+:+6/5-76B/+8 5/@+9>.<55 + :95C65?.7>+8 @H/1<55 f /7-.+=3. :+-97CH7> +- :95C65?.FG78 @H/1<55 Pf(Pf )00(x1 ) = 0(Pf )00(xn) = 0-.7.P100(x1 ) = 0 Pn00;1(xn) = 0:2 H9.0/7/58 (15) :95 f 00(x1 ) = 0 5 f 00 (xn) = 0 :+6H4.7>2d1 + d2 = 3f (x1Y x2)dn;1 + 2dn = 3f (xn;1Y xn)J-5 H9.0/7/5; 5,:+6B2HF-,; -.1 ?7, 1.1 H9.0/7/5; (15).,65 H :95C65?.7>+8 @H/1<55 f 00(x1 ) 6= 0 565 f 00 (xn) 6= 0, -+ :9+520+6 00DC+97 =9./54/DA H,6+058 :950+35- 1 :+0DE7//+8 :+=97E/+,-5 0 :95=9./54/DA H26.A (:+9;31. O(h2)), 4-+ +,+C7//+ :6+A+ +-9.?.7-,; /. :95C65?7/55,:6.8/.>5, =37 0 ,56H 7=+ /76+1.6B/+,-5 I-. :+=97E/+,-B CH37- 0+ 0,7A -+41.A+-9721..x 19.8.4.!% \ %;6!" % 5$)$&6!J-+- >7-+3 5,:+6B2H7-,;, 1+=3. /51.1+8 3+:+6/5-76B/+8 5/@+9>.<55 + :95C65?.7>+8 @H/1<55 f /7-.+=3.
:+-97CH7> +- :95C65?.FG78 @H/1<55 PfP1 P2 Pn;2 Pn;1(16)-.7. >/+=+467/ P1 P2 5/-79:+659H7- f /. Sx1 x3 ] :+ -+41.> x1 x2 x3 :P1(x1 ) = f (x1 ) P1(x2) = f (x2 ) P1(x3 ) = f (x3 ) 5 P10 (x3 ) = d3 , . >/+=+467/ Pn;2 Pn;1 5/-79:+659H7- f /. Sxn;2 xn] :+ -+41.> xn;2 xn;1 xn : Pn;1(xn;2) =f (xn;2) Pn;1(xn;1) = f (xn;1) Pn;1(xn) = f (xn) 5 Pn0 ;1(xn;2 ) = dn;2 . 673+0.-76B/+, >/+=+467/ P1 P2 ,-9+5-,; 1.1 5/-79:+6;<5+//D8 >/+=+467/ -97-B78 ,-7:7/5, :95/5>.FG58 0 -+41.A x1 x2 x3 2/.47/5; f (x1) f (x2) f (x3) 55>7FG58 0 -+417 x3 :9+520+3/HF, 9.0/HF d3 Y >/+=+467/ Pn;1 Pn;2 ,-9+5-,; 1.1 5/-79:+6;<5+//D8 >/+=+467/ -97-B78 ,-7:7/5, :95/5>.FG58 0 -+41.Axn;2 xn;1 xn 2/.47/5; f (xn;2) f (xn;1) f (xn) 5 5>7FG58 0 -+417 xn;2 :9+520+3/HF, 9.0/HF dn;2 ..15> +C9.2+>, 0>7,-+ n ; 1 >/+=+467/.
Pi i = 1 2 : : : n ; 1 5 n :.9.>7-9+0di i = 1 2 : : : n 5>77> n ; 3 >/+=+467/. P1 P2 P3 P4 : : : Pn;3 Pn;2 Pn;1 5n ; 4 :.9.>7-9. d3 d4 : : : dn;2 .7> ,.>D> :9+C67>. =9./54/DA H,6+058 0 >7-+3.A, 526+?7//DA 0 x 19.4,x 19.5 5 x 19.6, ,/;-., :+,1+6B1H :.9.>7-9D di i = 3 4 : : : n ; 2 +:97376;F-,;:+ :950737//D> 0 >7-+3.A @+9>H6.> (,>. (7), (8), (9))... !"# # $#x19.
'&+-.- 0*-0* 70.-77 >7-+37 :95C65?7/5; 1HC547,15>5 ,:6.8/.>5 (x 19.7) 0+2>+?/D 30. :+3A+3..1. 2 H,6+05; Pf 2 C (2) (Sa b]) :+6H4.7>Pi00;1(xi ) = Pi00(xi ) i = 3 4 : : : n ; 2{ n;4 H9.0/7/5; +-/+,5-76B/+ n;4 /75207,-/DA d3 d4 : : : dn;2 . 9.0/7/5; 36;i = 4 5 : : : n ; 3 ,+0:.3.F- , 0D:5,.//D>5 9./77 H9.0/7/5;>5 (12), H9.0/7/5;36;P100(x3 ) P200(x3 ) = P300(x3) Pn00;3(xn;2) = Pn00;2(xn;2 ) Pn00;1(xn;2)52>7/;-,; 0 ,56H 52>7/50E7=+,; 053.
>/+=+467/+0 P1 P2 , Pn;2 Pn;1 .2. H37> 2.3.0.-B P1 P2 : : : Pn;1 1.1 9./BE7 /.C+9+> d1 d2 : : : dn . D:+6/7/5; H,6+05; (16) 3+CO7>,;, :+-97C+0.0 /7:979D0/+,-5 (Pf )000 0 -+41.A x2 5xn;1 , -.7.P1000(x2 ) = P2000(x2 ) Pn000;2(xn;1 = Pn000;1(xn;1)(17)(H,6+05; (16) 5 (17) I1050.67/-/D 36; >/+=+467/+0 -97-B78 ,-7:7/5, :+,1+6B1H0 ,56H (10) H P1 5 P2 , Pn;2 5 Pn;1 ,+0:.3.F- 0,7 :9+520+3/D7 +- /H670+=+ 3+-97-B7=+ :+9;31. 016F45-76B/+).++-/+E7/5; (17) 3.F- 30. /73+,-.FG5A 36; 2.>D1./5; ,5,-7>D (12) H9.0/7/5;.2 :973,-.067/5; >/+=+467/. Pi 0 0537 (2) /.A+35>, 4-+ ,++-/+E7/5; (17)I1050.67/-/D(18)c4 1 = c4 2 c4 n;2 = c4 n;1.) .,,>+-95> :790+7 H9.0/7/57 (18). :+>+GBF (3) /.A+35>d1 + d2 ; 2f (x1 Y x2) = d2 + d3 ; 2f (x2 Y x3) (x2 ; x1 )2(x3 ; x2 )2565(x3 ; x2 )2d1 + (x3 ; x2 )2d2 ; (x2 ; x1 )2d2 ; (x2 ; x1)2 d3 == 2f (x1Y x2 )(x3 ; x2 )2 ; 2f (x2 Y x3)(x2 ; x1 )2-.7.(x3 ; x2 )2d1 + (x3 ; x1 )(x3 ; 2x2 + x1 )2 d2 ; (x2 ; x1 )2d3 == 2f (x1Y x2 )(x3 ; x2 )2 ; 2f (x2 Y x3)(x2 ; x1 )2,65 3+C.05-B I-+ H9.0/7/57 1 ,5,-7>7 (12), -+ /.9HE5-,; 77 -97A35.=+/.6B/+,-B.
+I-+>H :95C.05> 1 I-+>H H9.0/7/5F :790+7 H9.0/7/57 ,5,-7>D (12)(-.7. :95 i = 2):(x3 ; x2 )d1 + 2(x3 ; x1)d2 + (x2 ; x1 )d3 == 3f (x1 Y x2)(x3 ; x2 ) + 3f (x2Y x3 )(x2 ; x1 ).. !"# # $#x19. '&+-.- 0*-0* 70.-78H>/+?7//+7 /. (x2 ; x1 ), :+6H45>(x3 ; x2 )(x3 ; x2 + x2 ; x1 )d1 + (x3 ; x1 )(x3 ; 2x2 + x1 + 2(x2 ; x1 ))d2 == f (x1 Y x2)(x3 ; x2 )(2(x3 ; x2 ) + 3(x2 ; x1 )) + f (x2 Y x3)(x2 ; x1 )2565(x3 ; x2 )d1 + (x3 ; x1 )d2 =2= f (x1Y x2 )(x3 ; x2 )(2x3 +xx2;;x3x1 ) + f (x2 Y x3)(x2 ; x1 ) :31J-+ H9.0/7/57 3+C.06;7-,; 1 ,5,-7>7 (12) 0 1.47,-07 :790+=+.C) .,,>+-95> 0-+9+7 H9.0/7/57 (18).
:+>+GBF (3) /.A+35>dn;2 + dn;1 ; 2f (xn;2Y xn;1) = dn;1 + dn ; 2f (xn;1Y xn) (xn;1 ; xn;2)2(xn ; xn;1 )2565(xn ; xn;1)2 dn;2 + (xn ; xn;1)2 dn;1 ; (xn;1 ; xn;2)2dn;1 ; (xn;1 ; xn;2 )2dn == 2f (xn;2Y xn;1)(xn ; xn;1)2 ; 2f (xn;1Y xn)(xn;1 ; xn;2 )2-.7.(xn ; xn;1)2 dn;2 + (xn ; xn;2)(xn ; 2xn;1 + xn;2)2 dn;1 ; (xn;1 ; xn;2)2 dn == 2f (xn;2Y xn;1)(xn ; xn;1)2 ; 2f (xn;1Y xn)(xn;1 ; xn;2 )2,65 3+C.05-B I-+ H9.0/7/57 1 ,5,-7>7 (12), -+ /.9HE5-,; 77 -97A35.=+/.6B/+,-B. +I-+>H 0D4-7> I-+ H9.0/7/57 52 :+,673/7=+ H9.0/7/5; ,5,-7>D (12)(-.7.
:95 i = n ; 1):(xn ; xn;1)dn;2 + 2(xn ; xn;2 )dn;1 + (xn;1 ; xn;2)dn == 3f (xn;2Y xn;1)(xn ; xn;1) + 3f (xn;1Y xn)(xn;1 ; xn;2 )H>/+?7//+=+ /. (xn ; xn;1 ), :+6H45>(xn ; xn;2)(2(xn ; xn;1 ) ; (xn ; 2xn;1 + xn;2))dn;1++(xn;1 ; xn;2 )(xn ; xn;1 ; (xn;1 ; xn;2 ))dn =2= f (xn;2Y xn;1)(xn ; xn;1 ) + f (xn;1 Y xn)(xn;1 ; xn;2 )(3(xn ; xn;1 )+2(xn ; xn;2 ))565(xn ; xn;2)dn;1 + (xn;1 ; xn;2 )dn =f(xn;2 Y xn;1 )(xn ; xn;1 )2 + f (xn;1 Y xn )(xn;1 ; xn;2 )(3xn ; xn;1 ; 2xn;2 )=:x ;xnn;2J-+ H9.0/7/57 3+C.06;7-,; 1 ,5,-7>7 (12) 0 1.47,-07 :+,673/7=+... !"# # $#x19. '&+-.- 0x 19.8.5.*-0* 70.-5!!"&6! % 5$)$&796!J-+- >7-+3 5,:+6B2H7-,;, 1+=3. /51.1+8 3+:+6/5-76B/+8 5/@+9>.<55 + :95C65?.7>+8 @H/1<55 f /7-.0737> 3+:+6/5-76B/D7 -+415 x0 < x1 5 xn+1 > xn 5 2/.47/5; :95C65?.7>+8 @H/1<55 f (x0 ) 5 f (xn+1) 0 I-5A -+41.A.7> ,.>D> :9+C67>.
=9./54/DA H,6+058 0 >7-+3.A, 526+?7//DA 0 x 19.4 5x 19.6, ,/;-., :+,1+6B1H :.9.>7-9D di i = 1 2 : : : n -7:79B +:97376;F-,; :+:950737//D> 0 >7-+3.A @+9>H6.> (,>. (7), (9)). >7-+37 :95C65?7/5; 1HC547,15>5 ,:6.8/.>5 (x 19.7) 0+2>+?/D /7,1+6B1+:+3A+3+0.1. 0. 3+:+6/5-76B/DA H9.0/7/5; 7,-BP1(x0 ) = f (x0 )Pn;1(xn+1) = f (xn+1):2 :973,-.067/5; (2) >/+=+467/.
Pi 5 053. (3) 7=+ 1+I@@5<57/-+0 053/+, 4-+ 0I-5 H9.0/7/5; 0A+3;- ,++-07-,-07//+ d1 d2 5 dn;1 dn . 673+0.-76B/+, I-5 H9.0/7/5; >+?/+ 3+C.05-B 1 ,5,-7>7 (12) 0 1.47,-07 :790+=+ 5 :+,673/7=+ H9.0/7/58, ,+A9./7/57> 77 -97A35.=+/.6B/+=+ 053..2. /+=+467/ P1 5/-79:+659H7- f /. Sx0 x2] :+ -+41.> x0 x1 x2 : P1(x0 ) =f (x0 ) P1(x1 ) = f (x1 ) P1(x2 ) = f (x2) 5 P10 (x2 ) = d2 , . >/+=+467/ Pn;1 5/-79:+659H7- f /. Sxn;1 xn+1] :+ -+41.> xn;1 xn xn+1 : Pn;1(xn;1 ) = f (xn;1) Pn;1(xn) =f (xn) Pn;1(xn+1) = f (xn+1) 5 Pn0 ;1(xn;1) = dn;1 . 673+0.-76B/+, >/+=+467/ P1,-9+5-,; 1.1 5/-79:+6;<5+//D8 >/+=+467/ -97-B78 ,-7:7/5, :95/5>.FG58 0-+41.A x0 x1 x2 2/.47/5; f (x0 ) f (x1) f (x2) 5 5>7FG58 0 -+417 x2 :9+520+3/HF,9.0/HF d2 Y >/+=+467/ Pn;1 ,-9+5-,; 1.1 5/-79:+6;<5+//D8 >/+=+467/ -97-B78,-7:7/5, :95/5>.FG58 0 -+41.A xn;1 xn xn+1 2/.47/5; f (xn;1) f (xn) f (xn+1)5 5>7FG58 0 -+417 xn;1 :9+520+3/HF, 9.0/HF dn;1 ..15> +C9.2+>, 0>7,-+ n :.9.>7-9+0 di i = 1 2 : : : n 5>77> 5 n ; 2 :.9.>7-9.
d2 d3 : : : dn;1 .2 H,6+05; Pf 2 C (2) (Sa b]) :+6H4.7>Pi00;1(xi ) = Pi00(xi )i = 2 3 : : : n ; 1{ n;2 H9.0/7/5; +-/+,5-76B/+ n;2 /75207,-/DA d2 d3 : : : dn;1 . 9.0/7/5; 36;i = 3 4 : : : n ; 2 ,+0:.3.F- , 0D:5,.//D>5 9./77 H9.0/7/5;>5 (12), H9.0/7/5;36;P100(x2 ) = P200(x2 ) Pn00;2(xn;1 ) = Pn00;1(xn;1)52>7/;-,; 0 ,56H 52>7/50E7=+,; 053. >/+=+467/+0.3. D45,6;7> :95C65?7/57 1 f 0(x1 ) 5 f 0(xn ) ./.6+=54/+ :95C65?7/5F 1HC547,15>5 >/+=+467/.>5 7,,76; (,>.x 19.4).
9H=5>5 ,6+0.>5, >D 0D45,6;7> d1 5dn :+ @+9>H6.> (7) :95 i = 1 5 i = n . J-5 H9.0/7/5; 3+C.06;F-,; 1 ,5,-7>7 (12)0 1.47,-07 :790+=+ 5 :+,673/7=+ H9.0/7/58 , ,+A9./7/57> 77 -97A35.=+/.6B/+=+053.... !"# # $#x20. '&+-.- &+* &+,80/.47/5; f (x0) 5 f (xn+1 ) >+=H- /.A+35-B,; 52 ,++C9.?7/58 :795+354/+,-5,47-/+,-5 5 -.:.. ,65 /51.1+8 3+:+6/5-76B/+8 5/@+9>.<55 + :95C65?.7>+8@H/1<55 f /7-, -+ >+?/+ 5,:+6B2+0.-B I1,-9.:+6;<5F, /.:95>79, 65/78/HF:x0 = x1 ; (x2 ; x1 )f (x0 ) = f (x1 ) ; (f (x2) ; f (x1))xn+1 = xn + (xn ; xn;1 ) f (xn+1) = f (xn ) + (f (xn) ; f (xn;1 )):x 19.8.6.G$5!;7; % 5$)$&6!J-+- >7-+3 5,:+6B2H7-,;, 1+=3. /51.1+8 3+:+6/5-76B/+8 5/@+9>.<55 + :95C65?.7>+8 @H/1<55 f /7-.