К.Ю. Богачёв - Практикум на ЭВМ. Методы приближения функций (1133845), страница 4
Текст из файла (страница 4)
1.?3+> +-97217 Sxi xi+1 ] i = 1 : : : n ; k +/. ,+0:.3.7.. !"# # $#x6. '&+-.- "(0:-" 4*+19, >/+=+467/+> Pi ,-7:7/5 k , -.15>, 4-+ Pi(xj ) = f (xj ) j = i : : : i + k (-.7. Pi{ 5/-79:+6;<5+//D8 >/+=+467/ .=9./?. ,-7:7/5 k , :+,-9+7//D8 :+ -+41.>xi : : : xi+k ). . +-97217 Sxn;k xn] 5/-79:+659HFG.; @H/1<5; Pf ,+0:.3.7- ,Pn;k { 5/-79:+6;<5+//D> >/+=+467/+> .=9./?. ,-7:7/5 k , :+,-9+7//D> :+-+41.> xn;k : : : xn .+ 5/-79:+6;<5+//+8 @+9>H67 BF-+/.Pi(x) =i+kXj =i(x ; xi ) : : : (x ; xj;1)f (xiY : : : Y xj ):(1)+I@@5<57/-D Pi , i = n ; k n ; k ; 1 : : : 1 0D45,6;F-,; , 5,:+6B2+0./57>-.C65<D 9.23767//DA 9.2/+,-78:f (x1 ) nf (x1 Y x2)n=f (x2 ) nf (x1Y x2 Y x3): : :n=f (x2 Y x3): : : f (x1Y : : : Y xk+1)n=f (x3 ) nf (x2Y x3 Y x4): : :==f (x3 Y x4):::n=f (x4 ) n:::............:::f (xn;k;1)=nf (xn;k;1Y xn;k )=nf (xn;k ) =nf (xn;k;1Y xn;k Y xn;k+1) : : :n=f (xn;k Y xn;k+1)n: : : f (xn;k;1Y : : : Y xn;1 )......f (xn;k Y xn;k+1Y xn;k+2) : : :=n......: : : f (xn;k Y : : : Y xn)=f (xn;2Y xn;1Y xn): : :=f (xn;1) nf (xn;1Y xn) =f (xn) =/.4.67 0D45,6;F-,; 1+I@@5<57/-D >/+=+467/.
Pn;k :f (xn;k ) f (xn;k Y xn;k+1) : : : f (xn;k Y : : : Y xn):2. k(k +1) 0D45-./58 5 k(k +1)=2 3767/58 (:+,1+6B1H Pn;k { 5/-79:+6;<5+//D8>/+=+467/, :+,-9+7//D8 :+ k +1 -+417 xn;k : : : xn ). .-7> -.C65<. 2.:+6/;7-,;:+ 35.=+/.65: 0D45,6;F-,;f (xn;k;1) f (xn;k;1Y xn;k ) : : : f (xn;k;1Y : : : Y xn;1 ):+ @+9>H6.>f (xn;k;1Y : : : Y xn;k+i) = f (xn;k Y : : : Y xn;xk+i) ; ;f (xxn;k;1Y : : : Y xn;k+i;1) n;k+in;k;1i = 0 : : : k ; 1.. !"# # $#x7. . &4<20' 70+3 1'2. 2k 0D45-./58 5 k 3767/58 (I-H :9+<73H9H /.2D0.F- .6=+95->+> ,305=. :973,-.067/5; >/+=+467/.
0 @+9>7 BF-+/.), 5 3.677f (xn;k;2) f (xn;k;2Y xn;k;1) : : : f (xn;k;2Y : : : Y xn;2)::::::::::::f (x1 )f (x1Y x2 ): : : f (x1Y : : : Y xk+1):,7=+ /.3+ 2.:+6/5-B n ; k ; 1 35.=+/.678. . I-+ :+-97CH7-,; 2k(n ; k ; 1)0D45-./58 5 k(n ; k ; 1) 3767/58. CG57 2.-9.-D /. 0D45,67/57 0,7A 1+I@@5<57/-+0: k(k +1)+2k(n;k ;1) = 2kn;k(k +1) 0D45-./58 5 k(k +1)=2+k(n;k ;1) =kn ; k(k + 1)=2 3767/58.&! 6; $5!$ <= 7 Pf % x +,HG7,-06;7-,; :+ ,673HFG7>H .6=+95->H.
:97376;7> i 52 H,6+05; x 2 Sxi xi+1 ].,65 i > n ; k , -+ :+6.=.7> i = n ; k . +=3. (Pf )(x) = Pi(x). D45,67/57 2/.47/5; >/+=+467/. (1) ,-7:7/55 k 0 -+417 x +,HG7,-06;7-,; :+ ,A7>7 +9/79.Pi(x) = f (xi) + (x ; xi)(f (xiY xi+1 ) + (x ; xi+1 )(f (xiY xi+1Y xi+2 ) + : : : ++(x ; xi+k;2 )(f (xiY : : : Y xi+k;1) + (x ; xi+k;1)f (xiY : : : Y xi+k ) |) :{z: :})k;1, 2.-9.-+8 k 0D45-./58, k ,6+?7/58 5 k H>/+?7/58.7. ( " #x I-+> 9.23767 >D :+6H45> +<7/15 :+=97E/+,-5 5/-79:+6;<55 , :+>+GBF5/-79:+6;<5+//+=+ >/+=+467/. .=9./?., 2.:5,.//+=+ 0 @+9>7 BF-+/..f 2 C (Sa b]) fxi gni=1 { Sa b] , x 2 R1f (x) = Ln(x) + (x ; x1) : : : (x ; xn )f (x1Y : : : xnY x)Pn (x ; x ) : : : (x ; x )f (x Y : : : Y x ) . Ln (x) =1i;11i$ 1.
(1)i=1. 72 +=9./547/5; +CG/+,-5 >+?/+ ,45-.-B, 4-+ x 6= xj j =1 : : : n , :+,1+6B1H 0 -+41.A xj :+=97E/+,-B 5/-79:+6;<55 9.0/. /H6F 5 I-+,++-07-,-0H7- (1).95 n = 1 9.07/,-0+ (1) 5>77- 053f (x) = f (x1) + (x ; x1)f (x1 Y x)5 ,:9.073650+ 0 ,56H +:973767/5; (4.1) 9.23767//+8 9.2/+,-5.H,-B 9.07/,-0+ (1) 3+1.2./+ :95 n = k :f (x) = Lk (x) + (x ; x1 ) : : : (x ; xk )f (x1 Y : : : xk Y x):(2)6!"%..
!"# # $#x7. . &4<' 70+3 1'21 ,56H 67>>D (4.3)f (x1 Y : : : Y xk Y xk+1) :f (x1 Y : : : Y xk Y xk+1 x) = f (x1 Y : : : Y xk Y xx) ;;xk+1-,F3. 5>77>f (x1 Y : : : Y xk Y x) = f (x1Y : : : Y xk Y xk+1) + (x ; xk+1)f (x1 Y : : : Y xk Y xk+1 x):(3)+3,-.05> (3) 0 (2), :+6H45>f (x) = Lk (x)+(x ; x1) : : : (x ; xk )f (x1 Y : : : Y xk+1)+(x ; x1) : : : (x ; xk+1)f (x1 Y : : : Y xk Y xk+1 x):(4) ,56H 67>>D 4.5 ,:9.073650+ :973,-.067/57 (4.3)Lk (x) + f (x1 Y : : : Y xk+1)(x ; x1) : : : (x ; xk ) = Lk+1 (x)5 :+-+>H (4) :95/5>.7- 053f (x) = Lk+1(x) + (x ; x1 ) : : : (x ; xk+1 )f (x1Y : : : Y xk Y xk+1 x)-.7. 9.07/,-0+ (1) ,:9.073650+ 36; n = k + 1. 7> ,.>D> -7+97>.
3+1.2./. :+5/3H1<55.f 2 C (n) (Sa b]) fxi gni=1 { Sa b] , x 2 Sa b] = (x) 2 Sy1 y2] ( y1 = y1(x) =minfx1 : : : xn xg , y2 = y2(x) = maxfx1 : : : xn xg) , $ 2. (n)f (x) = Ln (x) + (x ; x1 ) : : : (x ; xn) f n!( ) (5)Ln(x) = P (x ; x1 ) : : : (x ; xi;1 )f (x1Y : : : Y xi) .ni=1.
+ -7+97>7 16!"%f (x) = Ln(x) + (x ; x1) : : : (x ; xn )f (x1Y : : : xnY x):(6)+ 67>>7 4.6 ,HG7,-0H7- = (x) 2 Sy1 y2] (=37 y1 = y1(x) = minfx1 : : : xn xg ,y2 = y2(x) = maxfx1 : : : xn xg) -.1+7, 4-+(n)f (x1 Y : : : Y xnY x) = f n!( ) :+3,-.06;; I-+ 0 (6), :+6H4.7> -97CH7>+7 9.07/,-0+ (5)... !"# # $#x8. )+3 )x' '3 0)+228. . H,-B 2.3./ /.C+9 -+471 xi : : : xi+k (/7 +C;2.-76B/+ 9.2654/DA). ) @H/1<55 f /H670+=+ :+9;31.
/.2D0.7-,; f (xi) =f (xi), 9.23767//+8 9.2/+,-BF k -=+ :+9;31. /.2D0.7-,;81>(k)(k )>>< k! f (xi) 7,65 xi = : : : = xi+k 5 f 2 C f (xiY : : : Y xi+k )= > f (xiY : : : Y xr;1Y xr+1Y : : : Y xi+k ) ; f (xiY : : : Y xs;1Y xs+1Y : : : Y xi+k )>xs ; xr>:7,65 ,HG7,-0HF- xs 6= xr (1)(I-+ +:973767/57 1+9971-/+, -.7. /7 2.05,5- +- 0DC+9. -+471 xr 5 xs , 0 ,56H,5>>7-955 9.23767//+8 9.2/+,-5 :+ ,0+5> .9=H>7/-.>, 3+1.2.-76B,-0+ ,>. 67>>H 4.3).5$#!f 2 C (k) , f (xiY : : : Y xi+k )k + 1 xi : : : xi+k .
1. ! -6!"%. 95 k = 0 9.23767//.; 9.2/+,-B f (xi ) ,+0:.3.7- , f (xi ),1+-+9.; /7:979D0/. :+ H,6+05F.H,-B 67>>. 3+1.2./. 36; 9.23767//DA 9.2/+,-78 :+9;31. k ;1. .,,>+-95>:9+520+6B/HF -+41H (xiY : : : Y xi+k ) 2 Rk+1 .,65 ,HG7,-0HF- i s r i + k -.157, 4-+ xs 6= xr , -+ 0 ,56H +:973767/5; (1)9.23767//.; 9.2/+,-B f (xiY : : : Y xi+k ) /7:979D0/. 1.1 65/78/.; 1+>C5/.<5; /7:979D0/DA (:+ :973:+6+?7/5F 5/3H1<55) @H/1<58:; f (xiY : : : Y xs;1Y xs+1Y : : : Y xi+k ) :f (xiY : : : Y xi+k ) = f (xi Y : : : Y xr;1Y xr+1Y : : : Y xi+xk ) ;s xr,65 xi = : : : = xi+k , -+ :+ +:973767/5F (1)f (xiY : : : Y xi+k ) = k1! f (k)(xi ):+,1+6B1H :+ H,6+05F f (k) /7:979D0/., -+ 36; 0,;1+=+ " > 0 ,HG7,-0H7- =(") > 0 -.1+7, 4-+ 36; 0,7A jxi ; j < jf k (xi) ; f k ()j < ":( )( )H,-B (x0i Y : : : Y x0i+k ) 2 Rk+1 { :9+520+6B/.; -+41., -.1.;, 4-+ jx0j ; xij < ,j = i i + 1 : : : i + k . ,56H 67>>D 4.6 ,HG7,-0H7- 2 Sy1 y2] (=37 y1 =minfx0i : : : x0i+k g , y2 = maxfx0i : : : x0i+k g ) -.1+7, 4-+(k)f (x0iY : : : Y x0i+k ) = f k!( ) :..
!"# # $#x9. '&+-.-23'3 0)+.1 1.1 :+ 0DC+9H (x0iY : : : Y x0i+k ) 2 Rk+1 Sy1 y2] (xi ; xi + ), -+jf (xiY : : : Y xi k) ; f (xiY : : : Y xi k)j = k1! jf k (xi) ; f k ()j < k1! "0+0+( )( )4-+ 5 +2/.4.7- /7:979D0/+,-B 9.23767//+8 9.2/+,-5 0 -+417 (xi Y : : : Y xi+k ).x9. H,-B 2.3./D -+415 a = x1 < x2 : : : < xn = b 5 2/.47/5; f (i) (xj ) i =0 : : : mj ; 1 j = 1 : : : n . 97CH7-,; :+,-9+5-B 5/-79:+6;<5+//D8 >/+=+467/nLm = Pf ,-7:7/5 m ; 1 = P mi ; 1 -.1+8, 4-+i=1L(mi) (xj ) = f (i) (xj ) i = 0 : : : mj ; 1 j = 1 : : : n:(1)J-+- >/+=+467/ /.2D0.F- , .
45,6. m1 : : : mn { x1 : : : xn ,++-07-,-07//+..1+8 >/+=+467/ 735/,-07//7/ ,9735 0,7A >/+=+467/+0 ,-7:7/5 m ; 1. 78,-05-76B/+, 7,65 Qm { >/+=+467/ ,-7:7/5 m ; 1, H3+067-0+9;FG58 H,6+05F (1),-+ 5A 9.2/+,-B Lm ; Qm 5>77- -+41H xj /H67> 19.-/+,-5 mj j = 1 : : : n . 67n3+0.-76B/+, >/+=+467/ Lm ; Qm ,-7:7/5 /7 0DE7 m ; 1 5>77- 0,7=+ P mi = mi=1/H678 (, H47-+> 19.-/+,-5).
+I-+>H Lm ; Qm = 0.$ 1. . (1), & Lm (x) =Lm , -mXi=1(x ; y1) : : : (x ; yi;1) f (y1Y : : : Y yi)(2)y1 = : : : = ym1 = x1 ym1 +1 = : : : = ym1 +m2 = x2 : : : yP =1;1 m +1 = : : : = ym = xn:(3)nii. +6+?5>6!"%x(ij") = xj + (i ; 1)" i = 1 : : : mj j = 1 : : : n(4)=37 0 < " < "0 "0 { -.1+7, 4-+ 0,7 -+415 fx(ij") g 9.2654/D 36; 0,7A " 2 (0 "0).+,-9+5> :+ m -+41.> fx(ij") g 5/-79:+6;<5+//D8 >/+=+467/ 0 @+9>7 BF-+/.Lm(") ,-7:7/5 m ; 1 -.1, 4-+L(m") (x"ij ) = f (x"ij ) 36; 0,7A i = 1 : : : mj j = 1 : : : n 0 < " < "0:..(5) !"# # $#x9. '&+-.-24'3 0)+ ,56H 735/,-07//+,-5 -.1+=+ >/+=+467/.(")Lm (x) =mXi=1") ) f (y (")Y : : : Y y (") )(x ; y1(")) : : : (x ; yi(;11i(6)36; 0,;1+=+ ,:+,+C. :797/H>79+0.-B /.C+9+> y1(") : : : ym(") /.C+9 fx(ij") g .
+I-+>H/.> 3+,-.-+4/+ 3+1.2.-B, 4-+ ,++-/+E7/5; (1) 0D:+6/7/D 36; -+415 x1 .+6+?5>y1(") = x1("1) : : : ym("1) = xm1 1 ym("1) +1 = x(1"2) : : : ym("1) +m2 = x(m")2 2 : : : (7)(")(")( ")(")yP;1 m +1 = x1 n : : : ym = xm n :n=1ini2 (7), (4), (3) 0D-71.7-, 4-+ "limyj(") = yj 36; 0,7A j = 1 : : : m . ,56H /7:97!09D0/+,-5 9.23767//DA 9.2/+,-78 (67>>. 8.1) 0 9.07/,-07 (6) >+?/+ :7978-5 1:97376H :+ " ! 0 (:+6H450E58,; 0 972H6B-.-7 >/+=+467/ +C+2/.4.7> Lm ):mXLm (x) = (x ; y1) : : : (x ; yi;1) f (y1Y : : : Y yi) =i=1m1X= (x ; x1 )i;1f (|x1 Y :{z: : Y x1})+i=1i 9.2mX+(x ; x1 )m1(x ; ym1 +1) : : : (x ; yi;1)f (y1Y : : : Y yi) =i=m1 +1(8)m1 (i;1) x )X1(x ; x1 )i;1+= f(i ; (1)!i=1mX+(x ; x1 )m1(x ; ym1 +1) : : : (x ; yi;1)f (y1Y : : : Y yi)i=m1 +1(0 :+,673/7> 9.07/,-07 5,:+6B2+0./+ +:973767/57 (8.1) 9.23767//+8 9.2/+,-5 ,19.-/D>5 H26.>5).
2 I-+=+ :973,-.067/5; ,673H7L(mi) (x1 ) = f (i) (x1 ) i = 0 : : : m1 ; 1(:+,1+6B1H :9+520+3/D7 3+ :+9;31. m1 ; 1 016F45-76B/+ 0-+9+=+ ,6.=.7>+=+0 (8) 9.0/D /H6F 0 -+417 x1 ). 7> ,.>D> 3+1.2./+, 4-+ >/+=+467/ 053. (2)H3+067-0+9;7- H,6+05;> (1) 0 -+417 x1 .5$#!. .3.0.7>D8 9.07/,-0.>5 (2), (3) >/+=+467/ /.2D0.7-,; 0 @+9>7 BF-+/..$ 2.
f 2 C (m) (Sa b]) (m = max mj ) fxi gni=1 { j =1 ::: n Sa b] , x 2 R1f (x) = Lm (x) + (x ; x1 )m1 : : : (x ; xn)m f (x| 1 Y :{z: : Y x1}Y : : : x| nY :{z: : Y xn}Y x) (9)nLm (x) ..(2), (3).m1mn !"# # $#x9. '&+-.-25'3 0)+. .:5E7> 9.07/,-0+ (9) 0 I1050.67/-/+8 @+9>7f (x) = Lm (x) + (x ; y1) : : : (x ; ym)f (y1Y : : : ymY x)(10)=37 -+415 y1 : : : ym 2.3.F-,; (3), 5 CH37> 3+1.2D0.-B 0>7,-+ (9) 9.07/,-0+ (10).95 I-+> >+?/+ ,45-.-B, 4-+ x 6= yj j = 1 : : : m , :+,1+6B1H 0 -+41.A yj:+=97E/+,-B 5/-79:+6;<55 9.0/. /H6F 5 I-+ ,++-07-,-0H7- (10). .677 3+1.2.-76B,-0+ :+0-+9;7- 3+1.2.-76B,-0+ -7+97>D 7.1.95 m = 1 9.07/,-0+ (9) 5>77- 053(11)f (x) = f (y1) + (x ; y1)f (y1Y x)5 ,:9.073650+ 0 ,56H +:973767/5; (8.1) 9.23767//+8 9.2/+,-5.H,-B 9.07/,-0+ (9) 3+1.2./+ :95 m = k :f (x) = Lm (x) + (x ; y1) : : : (x ; yk )f (y1Y : : : yk Y x):(12) ,56H 67>>D (4.3); f (y1Y : : : Y ykY yk+1) :f (y1Y : : : Y yk Y yk+1 x) = f (y1Y : : : Y yk Y xx) ;yk+1-,F3.