К.Ю. Богачёв - Практикум на ЭВМ. Методы приближения функций (1133845), страница 15
Текст из файла (страница 15)
'&+-.- ')3 &)(-x 1.3.94!)$ $5!;7 6$& 5$6%#;2 -7+97>D 1 0D-71.7- ,673HFG58 .6=+95-> 0D45,67/5; 5,1+>DA 1+I@@5<57/-+0 ; = (i j )(i=1 ::: n j=1 ::: m) @H/1<55 Pf .C+2/.45> fkl = kl (f ), fl = (f1l : : : fnl)t { 071-+9 52 I67>7/-+0 fkl , F =Sf1 : : : fm ] { n m >.-95<., ,+,-.067//.; 52 071-+9+0 fk . +=3. @+9>H6. (5)H-079?3.7-, 4-+; = A;1 F (B t);1:.,,>+-95> >.-95<HT = A;1 F = A;1Sf1 : : : fm] = SA;1 f1 : : : A;1fm ]: ,56H 972H6B-.-.
(I.1.4) -7+97>D (I.1.1) 7,65 b = (b1 : : : bn)t { :9+520+6B/D8071-+9, -+ 071-+9c = A;1 b 7,-B 071-+9 1+I@@5<57/-+0 5/-79:+659HFG78 @H/1n<55 u = P ciui -.1+8, 4-+ k (u) = bk , k = 1 2 : : : n . 673+0.-76B/+, k -8i=1,-+6C7< >.-95<D T 7,-B 071-+9 1+I@@5<57/-+0 5/-79:+659HFG78 @H/1<55,:+,-9+7//+8 :+ k ->H ,-+6C<H fk >.-95<D F ..,,>+-95> >.-95<Ht; = T (B t);1 = B ;1 T t :,65 +C+2/.45-B T t = St1 : : : tn] { ,-+6C<D ti 2 Rm >.-95<D T t , -+ B ;1T =SB ;1 t1 : : : B ;1tn].
,56H 972H6B-.-. (I.1.4) -7+97>D (I.1.1) 7,65 b = (b1 : : : bm )t{ :9+520+6B/D8 071-+9, -+ 071-+9c = B ;1b 7,-B 071-+9 1+I@@5<57/-+0 5/m-79:+659HFG78 @H/1<55 u = P cj uj -.1+8, 4-+ l(u) = bl , l = 1 2 : : : m .j =1673+0.-76B/+, l -.; ,-9+1. >.-95<D ; 7,-B 071-+9 1+I@@5<57/-+0 5/-79:+659HFG78 @H/1<55, :+,-9+7//+8 :+ l -+8 ,-9+17 tl >.-95<D T ..-95<D T 5 ; +CD4/+ A9./;-,; /.
>7,-7 >.-95<D F . I-+> ,6H4.7 +:5,.//D8 0DE7 .6=+95-> >+?7- CD-B ,@+9>H659+0./ 0 ,673HFG7> 0537.1) 6; 0,7A k = 1 2 : : : n 2.>7/5-B k -8 ,-+6C7< fk >.-95<D F /. 071-+9nk 1+I@@5<57/-+0 5/-79:+659HFG78 @H/1<55 u = P kiui -.1+8, 4-+i=1i(fk ) = fik = i(u), i = 1 2 : : : n .2) 6; 0,7A l = 1 2 : : : m 2.>7/5-B l -HF ,-9+1H tl >.-95<D F /. 071-+9ml 1+I@@5<57/-+0 5/-79:+659HFG78 @H/1<55 v = P lj vj -.1+8, 4-+j =1j (tl ) = tjl = j (u), j = 1 2 : : : m ...
!"# $B $Bx1. '&+-.- ')3 &)(-x 1.4.95-$)$; $!67; !)$H,-B :+3:9+=9.>>. INTERP X(LU, N, ALPFA) :+ 071-+9HLU = (i(u)) (i=1 ::: N)2/.47/58 @H/1<5+/.6+0 i /. /71+-+9+8 @H/1<55 u 2 F (X R1) ,-9+5- 071-+9 ALPFA = (i ) (i=1 ::: N) 1+I@@5<57/-+0 5/-79:+659HFG78 @H/1<55 PX u =PN u :+ /71+-+9+>H >7-+3H 1 65/78/+8 5/-79:+6;<55 @H/1<58 52 F (X R1).i=1 i iH,-B :+3:9+=9.>>. INTERP Y(MV, M, BETA) :+ 071-+9HMV = (j (v)) (j=1 ::: M)2/.47/58 @H/1<5+/.6+0 i /.
/71+-+9+8 @H/1<55 v 2 F (Y R1) ,-9+5- 071-+9 BETA = (i) (j=1 ::: M) 1+I@@5<57/-+0 5/-79:+659HFG78 @H/1<55 PY v =PM v :+ /71+-+9+>H >7-+3H 2 65/78/+8 5/-79:+6;<55 @H/1<58 52 F (Y R1).j =1 j jH,-B f 2 F (X Y R1). C+2/.45> N M >.-95<HNF = ( (i j )(f ) )(i=1 ::: N j=1 ::: M),1+>.; N M >.-95<. 1+I@@5<57/-+0 GAMMA = ; @H/1<55 PfNXMX(Pf )(x y) =i j ui(x)vj (y)i=1 j =10D45,6;7-,; ,673HFG5> .6=+95->+>1. 6; i = 1 N 0D:+6/;-BINTERP X(i{8 ,-+6C7< NF, N, i{8 ,-+6C7< GAMMA)2.
6; j = 1 M 0D:+6/;-BINTERP Y(j{; ,-9+1. GAMMA, M, j{; ,-9+1. GAMMA),65 >7-+3 1 565 2 -97CHF- 36; ,0+7=+ :9+0737/5; +C9.G7/5; /71+-+9+8 >.-95<D, -+ I-. >.-95<. CH37- +3/+8 5 -+8 ?7 :95 1.?3+> 0D2+07 :+3:9+=9.>>DINTERP X 5 INTERP Y. 673+0.-76B/+, :95 :790+> +C9.G7/55 1 :+3:9+=9.>>7 INTERP X 565 INTERP Y >+?/+ :+,-9+5-B, /.:95>79, LU{9.26+?7/57 I-+8>.-95<D 5 ,+A9./5-B I-+ 9.26+?7/57 3+ ,673HFG7=+ 0D2+0. :+3:9+=9.>>D 0+0,:+>+=.-76B/+> >.,,507. 2.05,5>+,-5 +- 5,:+6B2H7>+=+ ;2D1. :9+=9.>>59+0./5; 65C+ :H/1- 1 65C+:H/1- 2 I-+=+ .6=+95->. >+?7- CD-B +,HG7,-067/ C72 3+:+6/5-76B/DA 0D45,65-76B/DA 2.-9.- /. 0D2+0 :+3:9+=9.>> 5/-79:+6;<55 :+ :797>7//+8 X 565 Y...
!"# $B $Bx1. '&+-.- ')3 &)(-96.:95>79, /. ;2D17 FORTRAN, 0 1+-+9+> >.,,50D A9./;-,; :+ ,-+6C<.>, :790D8 :H/1- .6=+95->. >+?/+ 2.:5,.-B 0 0537DO I=1,NCALL INTERP X( NF(1,I), N, GAMMA(1,I) )ENDDO. ;2D17 C, 0 1+-+9+> >.,,50D A9./;-,; :+ ,-9+1.>, 0-+9+8 :H/1- .6=+95->. >+?/+ 2.:5,.-B 0 0537for( j=0Y j < MY j++ ) INTERP Y( GAMMASj], M, GAMMASj] )Yx 1.5.$5!;7 6$& 5$6%#; % ! )! 55$7H,-B f = f (x y) { @H/1<5; 30HA :797>7//DA, 2.3.//.; /.
:9;>+H=+6B/517 Sa b] Sc d]. H,-B -.1?7 2.3./D -+415 a = x1 < x2 < : : : < xn = b 5c = y1 < y2 < : : : < ym = d , 5 5207,-/D 2/.47/5; f (xi yj ), i = 1 2 : : : n ,j = 1 2 : : : m @H/1<55 f 0 -+41.A (xi yj ). 9+37>+/,-959H7> :95>7/7/57 -7+97>D 1 36; 9.,,>+-97//DA 9./77 >7-+3+0 :95C65?7/5;.x 1.5.1.$5!;7& )! )$@H,-B 1.?3D8 52 >7-+3+0 :95C65?7/5; :+ :797>7//D> x 5 y ;06;7-,; .::9+1,5>.<578 >/+=+467/+> .=9./?. (,>.
x I.3). 95C65?.FG.; @H/1<5; (:+1.?3+8 52 :797>7//DA) 0D45,6;7-,; :+ @+9>H67 (I.3.1).0737> 0,7 H4.,-0HFG57 0 :+,-./+017 2.3.45 65/78/+8 5/-79:+6;<55 -7/2+9/D>5 :9+520737/5;>5 +CO71-D. .25,/D7 @H/1<55:n x;xYkui(x) = i = 1 2 : : : nk=1 k6=i xi ; xkm y;yYl j = 1 2 : : : m:vj (y) =l=1 l6=j yj ; yl"H/1<5+/.6Di (u) = u(xi) i = 1 2 : : : n j (v) = v(yj ) j = 1 2 : : : m:(7).-95<D A 0 (2) 5 B 0 (3) { 735/54/D7. .25, 0 :9+,-9./,-07 @H/1<58 { -7/2+9/DA :9+520737/58:nm x;x y;yYYklwij (x y) = ui(x)vj (y) =x;xy;yk jlk=1 k6=i l=1 l6=j ii = 1 2 : : : n j = 1 2 : : : m:.25, 0 :9+,-9./,-07 @H/1<5+/.6+0 { -7/2+9/DA :9+520737/58:ij (f ) = i(j (f )) = f (xi yj ) i = 1 2 : : : n j = 1 2 : : : m:.. !"# $B $Bx1.
'&+-.- ')3 &)(-97+ -7+97>7 1 :95C65?.FG.; @H/1<5; 5>77- 053 (4), =37 1+I@@5<57/-D ij0D45,6;F-,; 0 (5). ,56H A;1 = I , B ;1 = I 5>77>: ij = ij (f ). 5-+=7:95C65?.FG.; @H/1<5; 5>77- 053Pf =x 1.5.2.n XmXi=1 j =1f (xi yj )m x;x y;yYkl:x;xy;yk jlk=1 k6=i l=1 l6=j inY6#!& $6 #!; 7 ) 5$&.23767//D7 9.2/+,-5, 00737//D7 0 x I.8, ;06;F-,; (:+ 3+1.2.//+>H 0 x I.8) 65/78/D>5 /7:979D0/D>5 @H/1<5+/.6.>5. D >+?7> 9.,,>+-97-B 5A -7/2+9/+7:9+520737/57:f (xiY : : : Y xi+k yj Y : : : Y yj+l) == (f (xiY : : : Y xi+k ))(yj Y : : : Y yj+l) = (f (yj Y : : : Y yj+l))(xiY : : : Y xi+k ):.07/,-0+ 237,B ,:9.073650+ 0 ,56H :973,-.067/5; (I.4.2). J-+ +:973767/57 I1050.67/-/+ ,673HFG7>H ./.6+=H +:97376/5; (I.8.1).5$#!.
H,-B 2.3./ /.C+9 -+471 xi : : : xi+k , yj : : : xj +l (/7 +C;2.-76B/+ 9.2654/DA). ) @H/1<55 f /H670+=+ :+9;31. /.2D0.7-,; f (xi yj ) = f (xi yj ), 9.23767//+8 9.2/+,-BF (k + l)-=+ :+9;31. /.2D0.7-,;f (xiY : : : Y xi+k yj Y : : : Y yj+l)8 k>1 @ f (x y Y : : : Y y ) 7,65 x = : : : = x 5 f 2 C (k+l)>j +lii+k>k! @xk i j>>(f (xiY : : : Y xr;1Y xr+1Y : : : Y xi+k yj Y : : : Y yj+l);>>>;f (xiY : : : Y xs;1Y xs+1Y : : : Y xi+k yj Y : : : Y yj+l))=(xs ; xr)>>><7,65 ,HG7,-0HF- xs 6= xr => l1 @ f (x Y : : : Y x y ) 7,65 y = : : : = y 5 f 2 C (k+l)>>i+k jjj +l>l! @yl i>>(f (xiY : : : Y xi+k yj Y : : : Y yr;1Y yr+1Y : : : Y yj+l);>>>;f (xiY : : : Y xi+k yj Y : : : Y ys;1Y ys+1Y : : : Y yj+l))=(ys ; yr )>>:7,65 ,HG7,-0HF- ys 6= yr x 1.5.3.$5!;7; $ ! "<H,-B 1.?3D8 52 >7-+3+0 :95C65?7/5; :+ :797>7//D> x 5 y ;06;7-,; .::9+1,5>.<578 >/+=+467/+> .=9./?., 1+-+9D8 0D45,6;7-,; :+ @+9>H67 BF-+/.
(,>. x I.5). 95C65?.FG.; @H/1<5; (:+ 1.?3+8 52 :797>7//DA) 0D45,6;7-,; :+ @+9>H67 (I.5.2)... !"# $B $Bx1. '&+-.- ')3 &)(-980737> 0,7 H4.,-0HFG57 0 :+,-./+017 2.3.45 65/78/+8 5/-79:+6;<55 -7/2+9/D>5 :9+520737/5;>5 +CO71-D. .25,/D7 @H/1<55:ui(x) =vj (y) =iY;1(x ; xk ) i = 1 2 : : : nk=1jY;1l=1(y ; yl)j = 1 2 : : : m(237,B ,45-.7-,;, 4-+ :9+520737/57 9.0/+ 1, 7,65 079A/58 :97376 >7/BE7 /5?/7=+). "H/1<5+/.6D 2.3.F-,; 0 (7).
3/.1+, H3+C/77 9.,,>+-97-B 39H=HF ,5,-7>H@H/1<5+/.6+0, 1+-+9.; 3.7- -H ?7 :95C65?.FGHF @H/1<5F. ,56H (I.10.1) 01.47,-07 @H/1<5+/.6+0 >+?/+ 02;-Bi(u) = u(x1Y : : : Y xi) i = 1 2 : : : nj (v) = v(y1Y : : : Y yj ) j = 1 2 : : : m:6; I-+8 ,5,-7>D @H/1<5+/.6+0 >.-95<D A 0 (2) 5 B 0 (3) { 735/54/D7. .25,0 :9+,-9./,-07 @H/1<58 { -7/2+9/DA :9+520737/58:;1iY;1 jYwij (x y) = ui(x)vj (y) =(x ; xk )(y ; yl)k=1 l=1i = 1 2 : : : n j = 1 2 : : : m:.25, 0 :9+,-9./,-07 @H/1<5+/.6+0 { -7/2+9/DA :9+520737/58:ij (f ) = i(j (f )) = f (x1 Y : : : Y xi y1Y : : : Y yj )i = 1 2 : : : n j = 1 2 : : : m:+ -7+97>7 1 :95C65?.FG.; @H/1<5; 5>77- 053 (4), =37 1+I@@5<57/-D ij0D45,6;F-,; 0 (5).
,56H A;1 = I , B ;1 = I 5>77>: ij = ij (f ). 5-+=7:95C65?.FG.; @H/1<5; 5>77- 053;1n XmiY;1 jYXPf (x y) =f (x1Y : : : Y xi y1Y : : : Y yj )(x ; xk )(y ; yl )(8)i=1 j =1k=1 l=1D45,67/57 9.23767//DA 9.2/+,-78 +,HG7,-06;7-,; ,+=6.,/+ x 1.4 :95>7/7/57>.6=+95->. 52 x I.5 :+ 1.?3+8 52 :797>7//DA x 5 y .x 1.5.4.$5!;7; "#%@ =;" )!H,-B 1.?3D8 52 >7-+3+0 :95C65?7/5; :+ :797>7//D> x 5 y ;06;7-,; 5/-79:+6;<578 "305?HG5>5,;" >/+=+467/.>5 (,>.
x I.6) ,-7:7/5 k 5 l ,++-07-,-07//+. J-5 >7-+3D ;06;F-,; /. 1.?3+> 52 +-9721+0 Sxi xi+1 ], i = 1 2 : : : n ;k ; 1, Sxn;k xn] 5 Syj yj+1], j = 1 2 : : : m ; l ; 1, Sym;l ym] .::9+1,5>.<578 5/-79:+6;<5+//D> >/+=+467/+> .=9./?., 2.:5,.//D> 0 @+9>7 BF-+/..+I-+>H :95C65?.FG.; @H/1<5; ,-9+5-,;, -.1 ?7, 1.1 0 :973D3HG7> 9.23767... !"# $B $Bx1.
'&+-.- ')3 &)(-x 1.5.5.$5!;7; $&996!H,-B 1.?3D8 52 >7-+3+0 :95C65?7/5; :+ :797>7//D> x 5 y ;06;7-,;.::9+1,5>.<578 >/+=+467/+> .=9./?. , 19.-/D>5 H26.>5 (,>. x I.9). 7+97>. I.9.1 :+20+6;7- 9.,,>.-950.-B -.1HF .::9+1,5>.<5F ./.6+=54/+ +:5,.//+80DE7 C72 19.-/DA H26+0.x 1.5.6.6!@ 5 )! B>&C%H,-B 1.?3D8 52 >7-+3+0 :95C65?7/5; :+ :797>7//D> x 5 y ;06;7-,; .::9+1,5>.<578 >/+=+467/.>5 7CDE70. (,>. x I.13). 95C65?.FG.; @H/1<5;(:+ 1.?3+8 52 :797>7//DA) 0D45,6;7-,; :+ @+9>H67 (I.13.10).0737> 0,7 H4.,-0HFG57 0 :+,-./+017 2.3.45 65/78/+8 5/-79:+6;<55 -7/2+9/D>5 :9+520737/5;>5 +CO71-D. .25,/D7 @H/1<55:!bTi(x) (x) = Ti 2x ; (b + a) i = 0 1 : : : n ; 1 b;a !(d + c) j = 0 1 : : : m ; 1Tbj(y) (y) = Tj 2y ;d ;c(9)"H/1<5+/.6D 2.3.F-,; 0 (I.13.4):DE(x)i(u) = Tbi(x) u n i = 0 1 : : : n ; 1DE(y)j (v) = Tbi(y) v m j = 0 1 : : : m ; 1=37hf ginx( )=nXi=1f (xi)g(xi)hf gimy( )=mXj =1f (yj )g(yj ):6; I-+8 ,5,-7>D @H/1<5+/.6+0 >.-95<D A 0 (2) 5 B 0 (3) { 35.=+/.6B/D7(,>.
(I.13.9)):A = diag n n2 : : : n2 2 Mn B = diag m m2 : : : m2 2 Mm :.25, 0 :9+,-9./,-07 @H/1<58 { -7/2+9/DA :9+520737/58:wij (x y) = ui(x)vj (y) = Tbi(x)(x)Tbj(y) (y)i = 0 1 : : : n ; 1 j = 0 1 : : : m ; 1:.25, 0 :9+,-9./,-07 @H/1<5+/.6+0 { -7/2+9/DA :9+520737/58:n XmXxy)ij (f ) = i (j (f )) =wij (xi yj )f (xi yj ) hwij f i(nm(10)i=1 j =1i = 0 1 : : : n ; 1 j = 0 1 : : : m ; 1..