Вычислительная математика и структура алгоритмов. Воеводин (2006) (1186027), страница 11
Текст из файла (страница 11)
7 # . #& #–,.86% #1/ $ .#& (# #" +"#&( #"%/ #"#$ %(# %#$# # #/ * # #" , # #$#/$&* $1/ [1].(# & $# (# # $ * ! *, #"& $&$1$( #"% &. #&$1$ ( #"* #"! % & $$ # (# #, 0% $#/.* +#*#"/*&* $1/ (*$ # % &%( #"#$& $ 0 . 4#- # !1,### $&)#1 '## (# $* # $ &!" * " +"#.!1 ## , $ #& ' ' % #- # !(# #1&!"& 1%!"#*, $ %*"!& ( $.# &% '#( #",%!0, *(#'"#,* #* # #.
/ 1&!# ( ),%!0- #, $(##( ' % # #. /. 5 #' %& $ 0 ## $$!0, #"!0 # $!0 &)0 (% '##* "##$.%! !"##$ $ 1 # !$ $(# %#$$#%0, #".4 $ & (#"#,#1 # !$(#'#". 4!$ ( #$ X % 'D. 4 &(# #1 , # "#$#&(# #0 $ #'D# # #/ !(# #/ (# $##$$$- * # *%/*- # $ #* +! ) /. #(# $##"! $ '!"#& #.# ,% *& 1 $ ##"!#"## (#$ #. 6! )#"! '#!"#& # #1,$ #, %$ % * # *. 1 # # # #& #–( $#& #"# %+! ) .#(!, # $ ( #$ X +#$ %! / #/ # ! K,'# #$ (#$ * ## # #"# #' %#$#/ " ( ( # # /. %#$- # # !'!&(,"# (## ! ' $( #$ X. &$# #"(# #"## # $&!0, :%+! ) $$ #x % #'D, #"!+! )#"! '.%* (# #* #(# $# , # #(#( & 0 $ #(#/ # ! , ( ! #$ , #$ .# !*#&$ #x.# #(#/ # ! '.- #!$"! ' !# # ,(& 1 , *# # #/##/ ## . 4# ##' %#"#* " +% $ 0., # $#"#" #%' $ 0(#%!+ ". ##"#,%$& !"#$# (# %!! ! , ##*, # #" #.
' !0(#&'#$#"#( #) #*( / (0."#,( & $/ $ #(# $ & ,#"# #$.'!( # $#/ # #'#/ % '# # "!#(# #1&!".870'#/ (#' /"#(#($,(#/ '# #'#',*/ * %$ # .$ * #1 # $% ,(,/+! ) #( $ 0,$ #,#$( & 0,( $ 0,$ #" ( ( # # /, #' %!0, *"#(# #"# # ! . 5*( & $/ ("#,0, "# ( $/ " +.(# % $(, #" # ##*% *"##$! & $$ #(!0 * !.*#&%& )#* #(/"##$,$ "& $ #* ! *%! $ &(#&*#& , #' %( #$ X, #'D, #(# #"# # ! K (# $#( '## . ##' % (#$ "&*#&.(# #, 0 $$ "#- #(#* ! #$ * ( &!( $0'!0 ( #"! #1 # $# $ #1# ! " %&!) #$. 5 # #1 (#& % , 1$ " &( #$# X#'D.#- #"#, $#% #1 #$ &( #"" %&! ) #$(#& % $ ,# (#-$ & # ! $##" *! * " + % ( #$) /, # # # $##', - # ! # # # # #1 #, #1 # ! #1$# - #( #$##/ %# .4"#$## # , #!0! & !# # #1 # # !, $# ##, #! #1" + ( - #0 $% , !0! ! !.$#,!0 1% , !0! ! !0"#.&(#" +"#% ( #$) / $ (#&*#& ,( #$##/%#% $% & /! .0#((&( ($ .
4 &(# #1 ,#1& /(# %!/"#'!& #(" +# , (# # # # !%! #1 #%$#! !#"#"#.# # ( #,*%!(# %#$!#"& & (#(#&*#& , * #' %#$ ( *#' ,- #()"##$.*#&% " ##/() , # $ ( $ $ . #(%! ! & $ 0$"#,(# %!0,$% #*, $ #,"!& ( $. #% # #$! $%!0*() 0,$$ #& #"# % $ 1 /. * #&! / $ '! #1% 1#/)$ # , #$#/()!& ##- # %" & .$ % #%$#%&!0,& . %$#, # #&%!& # / $# (*#$, # #'# $ # #.+#) # #/(#/! ! , $ %# ! $& #/# #"#/ # #$* % #1. 4 & #1*$(#(#&*#&/ # #$ .#( #&(#'!0( )%#' "#*.#' !1&* ! #1 #**- ( * # ( 0#"# #' %#$. (#1##& 1$(#$#% #1 # (##'*#& # .88%!" ##& + ) #$ , # % $.! ##(# #"## ! .4#'( #"##"&$# , "# *$&% ) &(#/ ( & #1#/ *#&, # #%&* (#/ $*49'("1#& 1$:,,.'#$/ # ##& + )(#%$#! # (% .#/ & , – - # # #"#"#. # #1 #. $#( # # $.(,,!!,""( )##"",,,-,-## #(" +#$ #*"##$.'!" +$ ( #$ *)/'!($*(#$ 0, +! ) .
# $#1 # (#% #$ [1].'!&% ," + $ (#&*#& , * ( #$ *#/ %# .'#% ,'!& #( #, #'#1 # ' # " # (# %+#) # !0! ! !*"##$( #& ##$$#% #1'#/ * %$ # .## $$#( &" +"#!1 # ! % $&!" , $ %#/ #& **1$*#& * &* (# $ ." + .## #(# % $'!&! , " $#' %# , %-%#1 #* %#' 1. ##"&'!& # !(# - #"# ( $ .3"5 "$. 4!.% &$( # %$ &A = BC &$!* $ &*) B, C (# & n. !) A, B, C#'#%* aij, bik, ckj. #"#( &0 #( ) ! #1)a ij =nb ik c kj , i , j = 1, 2 ,, n.k =15+# !$# ##(# %!0&(# & $#"#$#$) A.(# ' ##( & 0"##& #%#,#( &(# !#$( # %$ &/ bikckj(#& % # !. & #%, # $ #$ & (% $/.1# # ! $! %# #' 0"#- '# (# &( '#& #$ i, j.89#( )#1(# &!#$%!!. 4!%% ) +# ! :! #1$ (#- $ $( $#&*- # ##' 1/ $ '0# #, # $#& # !# ! 1&!0, /"#aij(0) = 0,aij(k) = aij(k –1) + bikckj,aij =i, j, k = 1, 2, …, n,aij(n).#$ $ # ! % (% ( '#& #$ i, j.
& # (# & ! k(%,- #& 1 ( '(# $## 1 n.4# #(" + - #"# "#. !&, #$ ." +## $$!0 #( )$ & a + bc. 8 #'$ #$$" +"##( $&(#!& # , $ ." +%(# "( # %$# #. 4/ (# #' *(# #1(#& % $+#% (.#(#!"# !0. ! $*#( #$## &i, j, k. # $ ) #!%.&1 ≤ i, j, k ≤ n (#$ ." + .% ! % ( ,!& #!' &$ # ,# $ $ . !## &i, j, k &k >1 '!&( & $%!$ (##( ) , ## $$!0, / $ .## &i, j, k – 1.90ijk. 9.1.+().#1! ## # ( # #.( &n2$ %*1&! #'#/ (#&" +#$.1& / (#&" + #& 1 n $ .( & $#& (! ,(# #1/(# # k.$ ##1 & , $ " + # 1# 1 (% , # # / '$ # ! % $% (.n = 2 " + "#( & $. 9.1.+(#/'#/ * %$ # , $$( $ 0, * $ # #$ # # * #"!'$%, (, ## &$ # .#$ # / i, j & 0#'#',%$, $ # $ # k –#"!0.)( )0 %"$. 4!(." '*/ * ! $/ Ax = b$ #1/!"# #/) / A (# & n# #'#/ (#$ .
'#%% aij, bi,xi ) , ( $#/$ # - .. 4 &(# #1( ,#) A$!"#& "#, $1. #"&+"#i −1aij x j 2 ≤ i ≤ n.x1 = b1, xi = bi —j =15 % ((# !#$% (#( $.(# $#%0&!0, #' %#!1&"###& #%(j. ## $,#,#( &,(# $#$!0, / "##1 #91xi(0) = bi,xi(j) = xi(j — 1) — aijxj(j — 1), i = 1, 2, …, n, j = 1, 2, …, i — 1,xi = xi(i — 1).(## & $ *#$#( )"#$ & a — bc.(!* %/& #$ i, j.##( ) # !, $ 0i, j (# #( $ $.$#/## (#!"# !0 ## &!0! .
"# 1 (#$ !%( 2≤ i ≤ n, 1 ≤ j ≤ i – 1 $ ." + , ## $$!0,#( )a – bc.% ! $ %1&! #( ), (# #" + "#, $ 0 $ $"#1 $ .,$# % !0, $$#& &* aij bi. 5 # " + &!n=5( & $. 9.2$ .*!" #$ $ .$#. 9.2 *#&$ ### &i = 1, j = 0.+(#/'#/ * %$ # , $$( $ 0, * $ # #$ # # * #"!i &'$%,(, ## &$ # .#$##'#',!0 %$!, $ # $ # j –#"!0.ji1b123j41ib123b2b2a21a21b3b3a32a31a32a31b4b4a42a41a43a42a41a52a51a53. 9.2.92a43b5b5(##' 1!#$! / #.#1 # (#!' $04$#&a54+ &a52a51!"#*a53a54.!#$ $(#- #(# $#(# #'!#$, # $ '## (# $#%0 &j'&, $##', "#$# ,$ - # $ '# (#$ & #*- '# (!, $, #"##'#/ (#$!#$(#j. ## $$!0, / "##$:xi(i) = bi,xi(j) = xi(j+1) — aijxj(j), i = 1, 2, …, n, j = i — 1, i — 2, …, 1,xi = xi(1).n=5 ( & $.
9.2 ( $ .($ *"# " + &!!" #$ $ .*#&$ ### &i = 1, j = 1.4% ,$ ., ## $$!0, #( )( a – bc, (#!*# '#/- '!& (#/ +#,#' !1 $, #($ 1!#1*# # #& $ .. '3- #, #$ $ ." +( $#. 9.2 1#& # (! . 5 #(! (# % (!.+ "## # #& !/ !0%$!, (#'#', !0.( $ 0, / $ #( $$ # i.4# !/%!!& # ' # #1 & ./ $#, #'#* "#( & %&.#& #/#/ 1 % & .(# ##& **1*#& * +# ! * $ # #.#*$/ - $ $. '"##$ . # #$#%*% )#& #( #) #* # ( 0*,'!0$ (##$#"##( ) / ! #1$#$#"# #'3(.!"#*#'"#& 1 - $ $#%$#. '# # !".," + #'# * "##$ ()(##0& !" # & !" ."#%#$ $(#/ $#/,0, / n ! $* ( #) # #$, # ( $ /"##1 #%#$% $( #(# ) ## n,$ # #/ – # # % $( #(# ) ## n 2.( $#!&% " !1#( #) # #$' %0,5, $# $ # #!# ' %0.#' %# ,"#, #$ . # #$#%% )"#'#$*" # ( 0*, #"! '()(#%#%% )(* # ( 0*.)$ ) + " -!' *! $0 "$ %"$." * % &/ #/" ' , $#%0, * (.! $/#/+ %##&,# $% &./ *" '*! $/' # #-&$!*& "#) .
4- # $ & "#' #( & $ 0 & "#) , & "#' # — &$!*& "#) . !&& #( , #)$$#/!"# #/. 4!#' #/ (# m (# ' # #$, $ / n.,#!/ * " '*! $/ Au = f&!0, "# $ & :93B1U1ΟD1 B 2D2B3ΟDm−1BmU2U3UmF1F2=F3,Fm"&d 1kb1kBk =e 2k b3kΟDk =,Οu1ku 2ku nk.' ##-&$!*& "##/,d 3ke n − 1, k bnkUk =Οd 2kΟe1k b2 kd nkf1k,Fk =f 2k.f nk#( &!#:–1Uk = Bk (Fk — Dk –1Uk –1), 1 ≤ k ≤ m,(# #1U0 = 0, D0 = 0.##( ) X = B –1(F — DY) $$ # X (# $ #F, Y) B, D. 4# # " + "#,,# 1&% "# $ .## $$! /##( ) .$ & #, #'!&,(# % #. 9.3. # .%$ .&!"!(#&$ 0 ,# #( )$ 0#1( +#) .94B1D0=0B2D1F1U1F2U2BDm-m1U0=0. 9.3.(#(#%#$$#" + &" +#($( & $.' # '&"#.) /.
& #&!(#&$!*& "##/) $ # #$' #%! $ & #,#-$ #FUmm#-& "##/.#* #("#$) /, # #$#"#($.41#( ) $ # .1&$!*& "##/. % - * + #$ #1 #$ $#& #$#% #1 #*# #. "#($$#/' # #-&$!*& "##/#/ $ $#& ' ' ( 1& $.!0% ("#.' # #. ! # $$ &*#'#%/#$:uik = (fik – ei – 1, kui – 1, k — di, k — 1ui, k — 1)bikk = 1, 2, …, m, i = 1, 2, …, n.4 - # ( &(# ", # ui0 = u0k = e0k = di0 = 0 #$ #/ , (# !, $!, &$#/ $$ * i, k. - #/ % (#( )u = b –1(f – ex – dy),f, e, x, y, b, d.(# ## #$$! u (# $" +"##(#!"# !0 .
!, !%# # #/0) ### &i, k. # $ !%.& 1 ≤ i ≤ n, 1 ≤> k ≤ m(#$ ." +'!&, #### $$!0 #( ) u.'!&! % $$ ., (# $ 0,$*#&&! $%# # * "!#$.% ! (#!0 % ( ,!& #!' &$ # ,# $ $ . !## &i, k '!&! ( & $%!$ (##( ) /,## $$!0, * $ .## &i –1, k i, k – 1.#+#) , #'*#&&% ) #( )## &i, k, $$*#& #/!1&% )# # - #/ #( ) .95ki. 9.4.+&' ##-&$!*& "#.#/ki. 9.5. 7 !#'#',#/ (#/ +#.!n = 5, m = 9 " + "#( & $.
9.4. % "#&! ,# $#($#% #1#1 &" +"#(#($.+(#/'# #'#',*/ *%$ # , $$( $ 0, * $ # #$ # # * #"! '$%## &$ # ./%$( $ 0,$ # # (1, 1)$#"#/.. 9.5 $ - # 1 " + (!#'$ &" !(( $ . .1&%* " !(( ## $$! #& #/#$ ." + , ( & $#"#. 9.3. % - *! #$ $ & #,#' %#*# #.#($/"##1( $$($/(!& # ! !(#( ) /.96'"$- )* $ ! $'" " .$#/ % && #& #'!/.u (y, z), "&' !" )#"# ! $.#( #( #$#& #&. 4!∂u∂ 2u=, 0 ≤ z, 0 ≤ y ≤ 1∂y∂z 2u(0, z) = ϕ(z), u(y, 0) = u0(y), u(y, 1) = u1(y).4# ## (#$ #!0(!. "# h (# z$ '# $ 'u j (i) − u j (i−1)τ4!"#%!#(# #i, j. 4" +"#y, z $ %h2τh,u j−1(i−1) − 2u j (i−1) + u j+1(i−1)=$ ## $u j ( i ) = u (ji −1) +.
"# τ (# y. 4 &(# #1$*2$+#! #/( u (ji−−11) − 2 u (ji −1) + u (ji+−11) ).$$ &((#!"#!0! ## &i, j ## #.y = τ2i, z = hj.4#$ 1& / !% ) ### $$!0, /ω = a (1 −$ (##/ &%τ = T/6 " +"##'(# #1$"* %/.%$ # .$#"! '$%, (/%$. /#/ #(.)$ .! " +'!&2ττ) +( b + c ),2hh 2* %/ "!#$ a, b, c.!h = 1/8,( & $#'# *! * 9.6.").,$# % !0, $$#&*&*+(#/'# #'#',*/ ** ( $ 0, * $ # #$ $## &i ,j, $ #(1, 1) (1, – 1).0#"% *% &( $ 0, $ # # (1, 0).97iijj.