x1+C12x2+ . . . +C1nxn=y1
x2+C23x3+…+C2nxn=y2 (4.11)
----------------------
xn=yn
теңдеулер жүйесін аламыз.
Қорыта келгенде (4.11) теңдеулер жүйесінің коэффициенттері мен оң жағы мына формулалар арқылы табылатының көреміз:
akj(0)=akj (k,j=1,2,…,n),
( j=k+1,k+2,…,n; k=1,…,n).
аіj(k)=aіj(k-1)- aіk(k-1)Ckj ( і,j=k+1,k+2,…,n; k=1,…,n),
fk(0)=fk , , (k=1,…,n),
fі(k)=fі(k-1)-aіk(k-1) yk (і=k+1,k+2,…,n.),
Ал (4.11) теңдеулер жүйесінің шешуі -х1,х2,…,хn мына формула арқылы табылады:
хn=yn , (і=n-1,n-2,…,1). (4.13)
Гаусс әдісімен теңдеулер жүйесін есептеу кезінде қолданылған арифметикалық амалдар саны -
N=m(m2+3m-1)/3 екенің есептеу онша қиын емес.
Гаусс әдісімен теңдеулер жүйесін шешкен кезде (4.12) формуласындағы бас элемент -акк(к-1) өте аз шама болған жағдайда дөңгелектеу қателері өсіп, есептеу дәлдігі төмендеп кетуі мүмкін. Сондықтан, кейбір жағдайда, бас элемент ретінде max|aіk(k-1)| (і=k,n) коэфициентін алады, яғни |хіk| коэфициенті ең үлкен теңдеуді к-шы теңдеумен орын ауыстыру арқылы есептеуді жалғастырады. Бұл әдіс Гаусстың тік жол боынша бас элементті таңдау әдісі деп аталады.
Егер akk -ны max|aіk(k-1)| (і=k+1,n) коэфициентімен алмастырсақ, яғни хк белгісізін хі белгісізімен орын алмастыру арқылы есептесе , онда бұл әдісті - Гаусстың жатық жол бойынша бас элементті таңдау әдісі дейді.
Лекция 8. Қуалау әдісі
Бұл әдіс көбінесе үш диогналды теңдеулер системасын шешуге арналған. Айталық,
(4.19)
(і=1,...,n), a1=cn=0
теңдеулер жүйесі берілсін. Бұл теңдеулер жүйесінің шешулерінің ара-қатынасы
і=xі+1 і+1+yі+1 (і=n,...,1), (4.20)
болсын делік. xі , yі белгісіздерін табу үшін (4.20) теңдігін (4.19) теңдігіне қойсақ, онда мынандай теңдеулер системасын аламыз:
.
Осы теңдеуден
. (4.21)
Жоғарыда көрсетілген (4.20) мен (4.21) формулаларын салыстыра отырып
, (4.22)
(мұнда і=1,...,n) формуласын аламыз. Егер x1=y1= n+1=0 десек , онда xі , yі ( і=1,...,n) белгісіздерін табу арқылы , (4.20) формуласын қолданып, теңдеулер жүйесінің шешуін табамыз.
Енді осы әдістің орнықты болуы үшін аі 0, сі 0,
|bі|>|aі|+|cі| ( і=1,...,n) (4.23)
шарттарының жеткілікті екенің көрсетейік.
Шынында да х1=0 болғандықтан |x2|=|c1|/|b1|<1.
Енді бұл шарт кез-келген і үшін |xі|<1 деп алып |x1+і|<1 теңсіздігі де орындалатынын көрсетейік. (4.23) формуладан теңсіздігін алуға болатындықтан .
Осыдан (4.22) формулаларының бөлімі нөлге тең еместігін көреміз және , яғни
(і=1,...,n). (4.24)
(4.23) шарты (4.19) жүйесінің шешуінің бар болуын және оның біреу ғана болатындығын қамтамасыз етіп қоймай, жүйенің орнықты екенің де көрсетеді.
Шынында да (4.20) формуласы бойынша і=і0+1 болғанда і+1-дің орнына есептелінді делік. Онда і=і0 болғанда формуласының орнына формуласын аламыз және .
Осыдан болатындықтан жіберілген қатенің өспейтіндігін көреміз.
Лекция 9. Матрицаның анықтауышын есептеу.
Гаусс әдісін матрицаның анықтауышын есептеу үшін де қолдануға болады. Біз оның бір ғана түрін қарастырамыз.
Айталық,
анықтауышы берілген және а11 0 болсын. Осыдан
, мұнда b11=a1і/a11 .
Енді бірінші жатық жолды әр төменгі жатық жолдың бірінші элементіне көбейтіп және сол жатық жолдан алып тастап отырсақ, онда
.
Егер а22.1 0 болса жоғарыдағы қолданған әдісті қалған
(n-1) өлшемді анықтауышқа қолданамыз. Осы әдісті аяғына дейін жалғастырсақ, онда
.
Егер аіі.і-1=0 немесе нөлге жақын болған жағдайда, алдын ала жатық жолдарды немесе тік жолдарды алмастыру арқылы аіі.і-1 0 матрицаға түрлендіруге болады. Ең дұрысы матрицаның бірінші жатық жолының бірінші элементі, оның барлық элементтерінен, абсолют шамасы бойынша ең үлкен болғаны тиімді. Анықтауышты есептеген кезде қолданылған көбейту мен бөлудің саны- .
Гаусс әдісі бойынша кері матрицаны табу.
Бізге ерекше емес
А=[aіj] (і,j=1,2,...,n) (4.31)
матрицасы берілсін.
Оның кері матрицасы
А-1=[хіj (і,j=1,2,...,n) (4.32)
болсын.
Онда АА-1 =Е теңдігі орындалатындықтан А мен А-1 матрицаларын көбейту арқылы мына теңдеулер жүйесін аламыз:
(і,j=1,2,...,n), (4.33)
мұнда
Теңдеулер жүйесінің оң жағы n бағанадан тұратындықтан Гаусс әдісін барлығына паралель қолдану тиімді, яғни
xkj (k,j=1,2,...,n) мына формулалар арқылы табылады:
akj(0)=akj, (k,j=1,...,n), (j=k+1,k+2,...,n; k=1,...,n),
(і,j=k+1,k+2,...,n; k=1,2,...,n-1). (4.34)
(4.35)
(j=k+1,k+2,...,n; і,k=1,2,...,n)
(і=n-1,n-2,...,1). (4.36)
Лекция 10. Кері матрицаны табу әдістері.
Гаусс әдісін қолдану арқылы кері матрицаны табу жолын өткен параграфта қарастырдық. Енді кері матрицаны табудың басқа жолдарын қарастырайық.
1>1>
Достарыңызбен бөлісу: |