К.Ю. Богачёв - Практикум на ЭВМ. Методы приближения функций (1133845), страница 18
Текст из файла (страница 18)
'&+-.- *3 C+'111(n + 1) (m + 1) >.-95<H F (xy) = (fij(xy) ) G(x) F (G(y) )t ./.47/5; @H/1<5+/.6+0 Fij = (ij kl(f ))k l=1 2 3 0 (14):0 (xy) (x)(xy)BB fi j fi j fi j+1Fij = BBB fi(yj) fi j fi(yj)+1@ (xy) (x) (xy)fi+1 j fi+1 j fi+1 j+11CCCCCA(15)$)67; $; #& % 5; G*9./7/57 5/@+9>.<55 + @H/1<55 5 77 :95C65?7/55 H3+C/+ +,HG7,-05-B,673HFG5> +C9.2+>. A+3/D7 3.//D7 { 2/.47/5; @H/1<5+/.6+0 (ij kl(f )), i =1 2 : : : n , j = 1 2 : : : m , k l = 1 2 3 4, A9./;-,; 0 47-D97A 00737//DA 0DE7>.-95<.A F 2 Mn m , F (x) 2 Mn+1 m , F (y) 2 Mn m+1 , F (xy) 2 Mn+ m+1 .
.-95<.Fij = (ij kl(f ))k l=1 2 3 :+6H4.7-,; ,+=6.,/+ (15). +I@@5<57/-D :95C65?.FG78@H/1<55 { (ij kl), i = 1 2 : : : n , j = 1 2 : : : m , k l = 1 2 3, 0D45,6;F-,; +35/9.2 5 A9./;-,; 0 4-A >79/+> >.,,507 ; = (ijkl), i = 1 2 : : : n , j = 1 2 : : : m ,k l = 1 2 3 9.2>79. (n m 3 3), =37 ijkl = ij kl .x2. " ,, ,.J-+- >7-+3 5/-79:+6;<55 :+20+6;7- :95C6525-B @H/1<55 >/+=5A :797>7//DA, 2.3.//+8 0 +C6.,-;A ,6+?/+8 @+9>D.
6; :9+,-+-D >D CH37> 9.,,>.-950.-B -+6B1+ @H/1<55 30HA :797>7//DA.x 2.1.-% 6# ! $5!;7H,-B 2.3./. @H/1<5; f : D ! R1 , =37 D R2 { +=9./547//.; +C6.,-B, 3+,-.-+4/+ =6.31+8 =9./5<78, f 2 F (D R1) { /71+-+9+>H @H/1<5+/.6B/+>H:9+,-9./,-0H.973,-.05> +C6.,-B D 0 0537 +CO735/7/5; 1+/74/+=+ 45,6. -97H=+6B/51+0Nk D , /7 5>7FG5A +CG5A 0/H-97//5A -+471: D0 S k D . +,-9+7/57k=1-.1+=+ /.C+9. -97H=+6B/51+0 (/.2D0.7>+=+ +C6.,-5 D ) 0 +CG7>,6H4.7 ;06;7-,; +47/B ,6+?/+8 2.3.478 5 36; 77 97E7/5; ,HG7,-0H7- C+6BE+745,6+ >7-+3+0. 5?7 >D 9.,,>+-95> /7,1+6B1+ ,.>DA :9+,-78E5A 52 /5A..677, /.
1.?3+> -97H=+6B/517 k @H/1<5; f :95C65?.7-,; -7> 565 5/D>>7-+3+>, -.7. ,-9+5-,; @H/1<5; Pk f , H3+067-0+9;FG.; -7> 565 5/D> 65/78/D> H,6+05;> 5/-79:+6;<55. 95C65?.FG57 @H/1<550 ,+,73/5A -97H=+6B/51.A 3+6?/D CD-B ,+=6.,+0./D: 7,65 (x y) 2 k1 T k2 { -+41. +CG78 =9./5<D-97H=+6B/51+0, -+ Pfk1 (x y) = Pfk2 (x y). /.47/57 :95C65?.FG78 @H/1<55.. !"# $B $Bx2. '&+-.- *3 C+'112Pf 0 -+417 (x y) 2 D0 0D45,6;7-,; ,673HFG5> +C9.2+>: +:97376;7-,; k , -.1+7, 4-+ (x y) 2 k , 5 :+ +:973767/5F :+6.=.7-,; Pf (x y) = Pk f (x y).
J-+:+3A+3 ;06;7-,; +C+CG7/57> 1H,+4/+->/+=+467//+8 .::9+1,5>.<55 /. ,6H4.8@H/1<58 >/+=5A :797>7//DA.+2>+?7/ 39H=+8 ,:+,+C :+,-9+7/5; :95C65?.FG78 @H/1<55. . +C6.,-50D +:97376;F- /.C+9 C.25,/DA @H/1<58 'l (x y) 5 :95C65?.FGHF @H/1<5F5GH- 0 0537 5A 65/78/+8 1+>C5/.<55. .25,/D7 @H/1<55 'l (x y) +CD4/+ 0DC59.F- -.1, 4-+CD 5A /+,5-76B ,+379?.6 1.1 >+?/+ >7/BE7 -97H=+6B/51+0 k(+- I-+=+ ,0+8,-0. C.25,/DA @H/1<58 :+E6+ /.20./57 \ %")..:95>79, 0DC59.7-,; -+41. (xl yl) 2 D0 5 C.25,/.; @H/1<5; 'l (x y) :+6.=.7-,; 9.0/+8 >/+=+467/H 2.3.//+8 ,-7:7/5 0 1.?3+> 52 -97H=+6B/51+0 k , -.15A,4-+ (xl yl) 2 k . +,-.6B/DA -97H=+6B/51.A 'l (x y) :+6.=.7-,; 9.0/+8 /H6F.,65 @H/1<5+/.6D l 0 2.3.47 65/78/+8 5/-79:+6;<55 7,-B 2/.47/5; @H/1<550 -+417: l (f ) = f (xl yl), -+ :95 -.1+> :+,-9+7/55 C.25,/DA @H/1<58 ,H?7/57:95C65?.FG78 @H/1<55 Pf /.
-97H=+6B/51 k 5 +:5,.//.; 0DE7 @H/1<5;Pk f ,+0:.3.F-.x 2.2.5>& 5$; $) !;7 >!5?7 >D +:5E7> /7,1+6B1+ :9+,-78E5A ,:+,+C+0 :+,-9+7/5; -95./=H6;<55+=9./547//+8 +C6.,-5 D , 3+,-.-+4/+ =6.31+8 =9./5<78.1) C6.,-B D +=9./547/., :+I-+>H ,HG7,-0H7- :9;>+H=+6B/51 Sa b] Sc d],-.1+8, 4-+ D Sa b]Sc d].
.?3.; 52 ,-+9+/ I-+=+ :9;>+H=+6B/51. 9.2C50.7-,;/. 9.0/D7 4.,-5 -+41.>5 xi = a + ihx , i = 0 1 : : : n , hx = (b ; a)=n , yj = a + jhy ,j = 0 1 : : : m , hy = (d ; c)=m . k ;06;7-,; -97H=+6B/51+> , 079E5/.>5 0-+41.A Qk 1 , Qk 2 , Qk 3 , 7,65 0D:+6/7/+ +3/+ 52 47-D97A H,6+058:Qk 1 = (xi yj ) 2 DQk 1 = (xi yj ) 2 DQk 1 = (xi yj ) 2 DQk 1 = (xi+1 yj ) 2 DQk 2 = (xi+1 yj ) 2 DQk 2 = (xi yj+1) 2 DQk 2 = (xi+1 yj ) 2 DQk 2 = (xi+1 yj+1) 2 DQk 3 = (xi+1 yj+1) 2 DYQk 3 = (xi+1 yj+1) 2 DYQk 3 = (xi yj+1) 2 DYQk 3 = (xi yj+1) 2 D:(1),6+057 k D >+?7- CD-B /7 0D:+6/7/+ 0 :95=9./54/DA -97H=+6B/51.A.+S6H47//.; -95./=H6;<5; CH37- -+4/+ .::9+1,5>59+0.-B +C6.,-B (-.7. D = k ),k7,65 D ;06;7-,; >/+=+H=+6B/51+> , 079E5/.>5 0 -+41.A (xi yj ), ,-+9+/D 1+-+9+=+ 65C+ :.9.6676B/D +,;> 1++935/.-, 65C+ ,+,-.06;F- , /5>5 H=+6 =4.2) D:+6/;7> +:5,.//HF 0 :H/1-7 1) -95./=H6;<5F.
6; :+0DE7/5; -+4/+,-5 .::9+1,5>.<55 +C6.,-5 9.,,>+-95> 0,7 :9;>+H=+6B/515 k 6 D , 079E5/.>5 Qk 1 = (xi yj ), Qk 2 = (xi+1 yj ), Qk 3 = (xi+1 yj+1), Qk 4 = (xi yj+1),-.15>5, 4-+ A+-; CD +3/. 52 -+471 Qk l 2 D . 305=./57> -+471 Qk l ! Q0k l ,l = 1 2 3 4 /. /.5>7/BEHF 0+2>+?/HF 076545/H :+ /.:9.067/5F 1 =9./5<7@D +C6.,-5 D 3+C50.7>,; -+=+, 4-+CD 65C+ 0,7 /+0D7 -+415 Q0k l 2= D , 65C+ 36;.. !"# $B $Bx2. '&+-.- *3 C+'1131.1+8-65C+ -9+815 Q0k l (565 30HA -9+71) CD6+ 0D:+6/7/+ +3/+ 52 H,6+058 (1).97H=+6B/515 , -.15>5 079E5/.>5 Q0k l 2 D 3+C.06;7> 1 -95./=H6;<55. S+6H47//.; -95./=H6;<5; CH37- -+4/+ .::9+1,5>59+0.-B +C6.,-B (-.7. D = k ),k7,65 D ;06;7-,; >/+=+H=+6B/51+> , 079E5/.>5 0 -+41.A (xi yj ).3) 973:+6+?5>, 4-+ >+?/+ :+,-9+5-B :97+C9.2+0./57 (=6.31+7 5 C72 +,+CDA-+471), :7970+3;G77 +C6.,-B G , 36; 1+-+9+8 >D H>77> ,-9+5-B -95./=H6;<5F(/.:95>79, :9;>+H=+6B/51), 0 9.,,>.-950.7>HF +C6.,-B D .
+=3. 079E5/D -97H=+6B/51+0 -95./=H6;<55 +C6.,-5 D >+=H- CD-B :+6H47/D 1.1 +C9.2D 079E5/-95./=H6;<55 +C6.,-5 G :95 I-+> :97+C9.2+0./55.4) 973:+6+?5>, 4-+ +C6.,-B D >+?/+ 9.2972.-B :9;>D>5 65/5;>5 /. :+3+C6.,-5, 0 1.?3+8 52 1+-+9DA >D >+?7> :+,-9+5-B -95./=H6;<5F. ,65 079E5/D -97H=+6B/51+0, 67?.G57 /. 65/5;A 9.2972+0, 0 1.?3+8 52 -95./=H6;<58,++-07-,-07//+ ,+0:.3.F-, -+ +CO735/7/57 I-5A -95./=H6;<58 3.7- -95./=H6;<5F 0,78 +C6.,-5 D .x 2.3.-$>!@ 7 % $ )!".,,>+-95> 2.3.4H :95C65?7/5; @H/1<55 f = f (x y), 2.3.//+8 0 -97H=+6B/517 .x 2.3.1.-% 6# ! $5!;7H,-B 1.?3.; ,-+9+/. -97H=+6B/51. 9.2C5-. /. l 9.0/DA 4.,-78, 5 47972-+415 9.2C57/5; (016F4.; 079E5/D -97H=+6B/51.) :9+0737/D :9;>D7, :.9.6676B/D7 7=+ ,-+9+/.>. C+2/.45> 47972 T >/+?7,-0+, ,+,-+;G77 52 -+471 :797,747/5; I-5AP :9;>DA, 67?.G5A 0 2.>1/H-+> -97H=+6B/517 .
5,6+ I-5A -+4719.0/+ n = l + 1k = (l + 1)(l + 2)=2. H37> +C+2/.4.-B 5A 47972 Q1 (x1 y1), : : : ,k=1Qn(xn yn)..,,>+-95> 2.3.4H :+,-9+7/5; >/+=+467/. Pf ,-7:7/5 lX ijPf (x y) =aij x y (2)i+j l:95/5>.FG7=+ 0 -+41.A Qk (xk yk ) 2.3.//D7 2/.47/5;Pf (xk yk ) = f (xk yk ) k = 1 2 : : : n:(3),6+05; (3) :973,-.06;F- ,+C+8 n 65/78/DA H,6+058 5/-79:+6;<55. +,1+6B1H 45,6+ /75207,-/DA 1+I@@5<57/-+0 aij -.1?7 9.0/+ n , -+ ,++-/+E7/5; (3) +C9.2HF- ,5,-7>H 52 n H9.0/7/58 , n /75207,-/D>5. ,65 2.3.4.
65/78/+8 5/-79:+6;<55 1+9971-/., -.7. I-. ,5,-7>. 9.297E5>., -+ 52 /77 >+?/+/.8-5 1+I@@5<57/-D aij .,1+>D8 >/+=+467/ Pf >+?/+ 0D:5,.-B ;0/+ 5 0 9.,,>+-97//+> 0DE7 +CG7> ,6H4.7. 3/.1+, I-+ 0D9.?7/57 3+,-.-+4/+ ,6+?/+ 5 /. :9.1-517 5,:+6B2H7-,; 9731+. 5?7 >D 0D:5E7> 053 >/+=+467/. Pf 0 ,6H4.7, 1+=3. -97H=+6B/51..
!"# $B $Bx2. '&+-.- *3 C+'114 { :9;>+H=+6B/D8 ,+ ,-+9+/.>5, :.9.6676B/D>5 +,;> 1++935/.-, . -.1?7 36;,6H4.70 l = 1 2 3.x 2.3.2.-$>!@ 7 % 5$; )!" $ )!"7E5> :+,-.067//HF 0DE7 2.3.4H 65/78/+8 5/-79:+6;<55 0 ,6H4.7, 1+=3.-97H=+6B/51 { :9;>+H=+6B/D8 ,+ ,-+9+/.>5, :.9.6676B/D>5 +,;> 1++935/.-.
H,-B 7=+ 079E5/.>5 ;06;F-,; -+415 (a c), (b c), (a d) (-.7. -97H=+6B/51+C9.2+0./ 3767/57> :9;>+H=+6B/51. Sa b] Sc d] 7=+ 35.=+/.6BF). H,-Ba = x1 < x2 < : : : < xn;1 < xn = b c = y1 < y2 < : : : < ym;1 < ym = d5 5207,-/D 2/.47/5; @H/1<55 f 0 -+41.A (xi yj ) 2 . 97CH7-,; :+,-9+5-B:95C65?.FGHF @H/1<5F Pf 053. (2), H3+067-0+9;FGHF H,6+05;>Pf (xi yj ) = f (xi yj ) (xi yj ) 2 :(4)J-. :+,-./+01. ;06;7-,; C+677 +CG78, 47> +:5,.//.; 0 :973D3HG7> 9.23767,1+-+9.; :+6H4.7-,; :95 9.0/+>79/+ 9.,:973767//DA -+41.A xi , yj 5 m = n .6; <7678 0D0+3.
.6=+95->. :95C65?7/5; >D CH37> ,45-.-B, 4-+ 5207,-/D2/.47/5; @H/1<55 f 0 -+41.A (xi yj ) 2= (/.:95>79, CH37> ,45-.-B f 9.0/+8/H6F 0 I-5A -+41.A). +=3. 36; @H/1<55 f :95>7/5>D .6=+95->D .::9+1,5>.<55 -7/2+9/D>5 :9+520737/5;>5 0 :9;>+H=+6B/517 Sa b] Sc d].
.> H3+C/+5,:+6B2+0.-B :95C65?7/57 5/-79:+6;<5+//D> >/+=+467/+> .=9./?. 0 @+9>7BF-+/., -.7. :95C65?.FGHF @H/1<5F (1.8), 1+-+9HF >D 2.:5E7> 0 ,673HFG7>0537i;1 j ;1f (x1Y : : : Y xi y1Y : : : Y yj ) Q Q (x ; xk )(y ; yl )k=1 l=1i j : (x y )2i;1 j ;1P+f (x Y : : : Y x y Y : : : Y y ) Q Q (x ; x )(y ; y )Pf (x y) =Piji j : (x y )2= ii11jjkk=1 l=1l(5)(237,B ,45-.7-,; 2.>1/H-D> >/+?7,-0+>). + :+,-9+7/5F, I-. @H/1<5; H3+067-0+9;7- H,6+05F Pf (xi yj ) = f (xi yj ) 36; 0,7A (xi yj ) 2 Sa b]Sc d].
.>7-5>,4-+;1iY;1 jY(x ; xk )(y ; yl ) = 0k=1 l=1(x y) = (xp yq )36; 0,7A -+471p = 1 2 : : : i ; 1 q = 1 2 : : : j ; 1:673+0.-76B/+, 0-+9.; ,H>>. 0 (5) 9.0/. 0 36; 0,7A (x y) = (xi yj ) 2 (:+,1+6B1H 36; 0,;15A (xi yj ) 2= , (xp yq ) 2 0D:+6/7/+ 65C+ xi > xp , 65C+yj > yq ). +I-+>H @H/1<5;Pf (x y) =Xi j : (x yi..j)2f (x1Y : : : Y xi y1Y : : : Y yj );1iY;1 jYk=1 l=1(x ; xk )(y ; yl )(6) !"# $B $Bx2. '&+-.- *3 C+'115H3+067-0+9;7- H,6+05;> (4), -.7. ;06;7-,; 97E7/57> 2.3.45 65/78/+8 5/-79:+6;<55.