Файл: азастан республикасы білім жне ылым министрлігі абай атындаЫ аза лтты педагогикалы университеті.docx
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 25.10.2023
Просмотров: 222
Скачиваний: 2
ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
-жатық бағана, rn-тік бағана және 1/ n-сан. Енді А мен А-1 матрицаларын бір-біріне көбейітсек, онда
.
Осыдан
An-1pn-1+unqn=E (5.1)
vnpn-1+annqn=0 (5.2)
An-1rn+un/ n=0 (5.3)
vnrn+ann/ n=1. (5.4)
(5.3) теңдігінен rn=An-1-1un/ n ,
ал (5.4) теңдігінен
nn=ann-vnAn-1-1un (5.5)
белгісіздерін табамыз.
Ал (5.1) теңдігінен
pn-1=An-1-1-An-1-1unqn. (5.6)
Енді (5.2) және (5.5) формулаларының негізінде
vnAn-1-1-vnAn-1-1unqn+annqn=vnAn-1-1-(ann- n)qn+annqn =
vnAn-1-1+ nqn=0.
болғандықтан qn=-vnAn-1-1/ n
pn-1=An-1-1+An-1-1unvnAn-1-1/ n.
, (5.7)
мұнда .
Бұл әдіс арқылы А-1 матрицасын табу ретімен
(а11), матрицаларының кері матрицасын көмкеру әдісі арқылы біртіндеп табу арқылы іске асырылады. Есептеу схемасы мынандай:
( n санын екі жолмен табу, есептеу барысының дұрыстығын қадағалап отыру үшін қажет.)
4) bіk=cіk+ іn nk/ n , (і,k -1) bіn= іn/ n ; bnk= nk/ n , (і,k -1)
dnn=1/ n.
Көмкеру әдісін, теңдеулер жүйесін шешкен кезде, Аn-1-1 белгілі болғанда қолдану тиімді. Мұндай жағдайлар Б.Г.Галеркин немесе В.Ритц әдісімен математикалық физика немесе механика есептерін (n-1) координатты функцияны қолданғандағы дәлдік қанағаттандырмаған, ал (n) координатты функцияны қолданғандағы дәлдік қанағаттанарлық болғанда кездеседі.
Теңдеулер жүйесін дәл әдіспен шешудің басқа да жолдары бар. Олардың кейбір түрлерімен Ө.Сұлтанғазин мен С.Атанбаевтың ,,Есептеу әдістерінің қысқаша теориясы” атты оқулығының 1-кітабынан танысуға болады.
Лекция 11. Теңдеулер жүйесін қарапайым итерациялық әдістермен шешу жолдары.
Итерациялық процесстерді құру принципі.
Көп жағдайда теңдеулер жүйесін дәл әдіске қарағанда итерациялық әдіспен шешкен тиімдірек. Себебі:
Енді осы итерациялық процесті құру схемасымен танысайық.
Айталық ,
(6.1)
теңдеулер жүйесі берілсін және делік. Осы теңдеулер жүйесін шешу үшін
(6.2)
формуласы бойынша ,( -итерацияның
-қадамдағы жуықтауы) векторлар тізбегін табамыз.
Х0-бастапқы кез-келген вектор, Н1, Н2,... матрицалар тізбегі .
(6.2)-формуласын итерациялық әдіс дейміз.
{Н (к)} матрицалар тізбегінің түріне қарай (6.2) формуласынан әртүрлі итерациялық процестерді аламыз. Егер ягни -дан тәуелсіз болса, онда итерациялық процесті “стационар” , ал -дан тәуелді болса “стационар емес” деп, ал болса, онда қадамды итерациялық процесс дейміз.
Итерациялық процестер үшін (6.1) теңдеулер жүйесінің шешуі Х*-қозғалмайтын нүкте болып табылады, яғни деп алсақ, онда (6.2) формула бойынша тапқан келесі векторымыз да болады.
Керісінше -қозғалмайтын нүкте болатын кез-келген мынандай итерациялық процесті
, (6.3)
мұнда -матрицалар тізбегі, -векторлар тізбегі, (6.2) түрінде жазуға болады.
Шынында да үшін
болғандықтан
мұнда .
Енді итерациялық процесстің жинақталуын қарастыру үшін (6.2) теңдігін векторынан алып тастайық, ягни
. (6.4)
Осы формуланы былайша жазсақ
онда -векторлар тізбегінің векторына жинақталуы үшін
матрицасының нөлге ұмтылуы қажетті және жеткілікті, ал матрицасының кез-келген нормасының нөлге ұмтылуы жеткілікті шарт екенің көреміз.Ал -векторлар тізбегінің векторына ұмтылу жылдамдығы матрицасына байланысты болғандықтан, оны тандау әртүрлі итерациялық әдістерге әкеледі.
Енді итерациялық әдістердің жинақтылығын қалай сипаттауға болады, немесе қай уақытта итерацмялық процессті тоқтату керектігіне тоқталайық. Теңдеулер жүйесінің дәлдігін қателік векторы арқылы бағалау мүмкін болмағандықтан көбіне -ауытқу векторы арқылы бағалайды. Шынында да, егер (C-const) болса, онда
.
Осыдан болғанда .
екенің көреміз. Ал
(6.6)
болғандықтан
.
Осыдан мынандай тұжырым жасауға болады: егер (C-const) яғни (6.1) орнықты болса, онда шегінен шегі келіп шығады. Ал матрицасының шарттылығы нашар болған жағдайда бұл тұжырым орындалмауы да мүмкін.
Сондықтан итерациялық процессті , алдын-ала берілген дәлдік - бойынша теңсіздігі орындалған кезде тоқтату ертерек болуы мүмкін.
Енді (6.4) теңдігін былайша жазайық: .
Осыдан .
Бұл теңсіздік ,(6.1) теңдеулер жүйесі орнықты болған жағдайда, итерациялық процесті
(6.7)
теңсіздігі орындалғанда тоқтатуға болатындығын көрсетеді.
Егер А оң анықталған симметриялы матрица болған жағдайда қателік вектордың шамасын қателік функциясы деп аталатын
функцияның шамасымен де анықтауға болады. Шынында да А оң анықталған болғандықтан және болғанда .
болғандықтан, оның мәнін табу мүмкін емес. Бірақ
(6.8)
функционалының қателік функциясынан айырмашылығы тұрақты сан болғандықтан және екеуі де өздерінің ең кіші мәндеріне болғанда ие болатындығын ескерсек, Ф(х) функционалына минимум мән беретін векторды жүйенің шешуі деп қарастыруға болады. Сондықтан
(6.9)
болған жағдайда да итерациялық процесті тоқтатуға болады.
Енді итерациялық процестердің қарапайым түрлерімен танысайық.
Біртіндеп жуықтау әдісі.
Берілген
теңдеулер жүйесін былайша жазайық
, (6.10)
мұнда .
Енді (6.10) теңдеулер жүйесін шешу үшін мынандай итерациялық әдісті пайдаланайық
, (6.11)
мұнда -өзіміздің еркіміз бойынша алынатын кез-келген бастапқы вектор. (6.11) итерациялық әдісін (6.2) формуласы түрінде жазсақ
,
онда екенің және бұл стационар итерациялық әдіс екенің көреміз.
Бұл әдістің жинақтылығын анықтау үшін, математикалық индукцияны қолдану арқылы, (6.11) теңдігінен алынған
(6.12)
және теңдігінен (6.11) теңдігін алып тастағанда алынған
(6.13)
формулаларын қарастырайық.
Енді осы формулаларды қолдана отырып әдістің жинақтылығы туралы бірнеше теоремаларды берейік.
6.1-теорема: Кез-келген -де біртіндеп жуықтау әдісінің жинақты болуы үшін матрицасының меншікті мәндерінің модулі бірден кіші болуы қажетті және жеткілікті.
Бұл теореманың шартын тексеру қиын болғандықтан, шарттары оңай анықталатын теоремаларды қарастырайық.
6.2-теорема: Біртіндеп жуықтау әдісі жинақталуы үшін болуы жеткілікті.
.
Осыдан
An-1pn-1+unqn=E (5.1)
vnpn-1+annqn=0 (5.2)
An-1rn+un/ n=0 (5.3)
vnrn+ann/ n=1. (5.4)
(5.3) теңдігінен rn=An-1-1un/ n ,
ал (5.4) теңдігінен
nn=ann-vnAn-1-1un (5.5)
белгісіздерін табамыз.
Ал (5.1) теңдігінен
pn-1=An-1-1-An-1-1unqn. (5.6)
Енді (5.2) және (5.5) формулаларының негізінде
vnAn-1-1-vnAn-1-1unqn+annqn=vnAn-1-1-(ann- n)qn+annqn =
vnAn-1-1+ nqn=0.
болғандықтан qn=-vnAn-1-1/ n
pn-1=An-1-1+An-1-1unvnAn-1-1/ n.
Сонымен
, (5.7)
мұнда .
Бұл әдіс арқылы А-1 матрицасын табу ретімен
(а11), матрицаларының кері матрицасын көмкеру әдісі арқылы біртіндеп табу арқылы іске асырылады. Есептеу схемасы мынандай:
-
( 1,n ,..., n-1,n)=-An-1-1un , ( n,1 ,..., n,n-1)=-vnAn-1-1 , -
.
( n санын екі жолмен табу, есептеу барысының дұрыстығын қадағалап отыру үшін қажет.)
4) bіk=cіk+ іn nk/ n , (і,k
dnn=1/ n.
Көмкеру әдісін, теңдеулер жүйесін шешкен кезде, Аn-1-1 белгілі болғанда қолдану тиімді. Мұндай жағдайлар Б.Г.Галеркин немесе В.Ритц әдісімен математикалық физика немесе механика есептерін (n-1) координатты функцияны қолданғандағы дәлдік қанағаттандырмаған, ал (n) координатты функцияны қолданғандағы дәлдік қанағаттанарлық болғанда кездеседі.
Теңдеулер жүйесін дәл әдіспен шешудің басқа да жолдары бар. Олардың кейбір түрлерімен Ө.Сұлтанғазин мен С.Атанбаевтың ,,Есептеу әдістерінің қысқаша теориясы” атты оқулығының 1-кітабынан танысуға болады.
Лекция 11. Теңдеулер жүйесін қарапайым итерациялық әдістермен шешу жолдары.
Итерациялық процесстерді құру принципі.
Көп жағдайда теңдеулер жүйесін дәл әдіске қарағанда итерациялық әдіспен шешкен тиімдірек. Себебі:
-
Егер итерациялық процесс жылдам жинақталатын болса, онда теңдеулер жүйесін шешуге жіберілген арифметикалық амалдар саны мен оны есептеуге жіберілген уақыт үнемделеді. -
Итерациялық әдіс өзін-өзі түзетін отыратындықта, есептеу кезінде жіберілген қателер жалпы теңдеулер жүйесінің шешуіне әсерін тигізбейді. -
Итерация әдісінің алгоритмі онша күрделі болмағандықтан ЕЭМ-ға программа оңай құрылады. -
Теңдеулер жүйесінің матрицасының элементтерінің көпшілігі нөлге тең болған жағдайда итерация әдісінің тиімділігі арта түседі.
Енді осы итерациялық процесті құру схемасымен танысайық.
Айталық ,
(6.1)
теңдеулер жүйесі берілсін және делік. Осы теңдеулер жүйесін шешу үшін
(6.2)
формуласы бойынша ,( -итерацияның
-қадамдағы жуықтауы) векторлар тізбегін табамыз.
Х0-бастапқы кез-келген вектор, Н1, Н2,... матрицалар тізбегі .
(6.2)-формуласын итерациялық әдіс дейміз.
{Н (к)} матрицалар тізбегінің түріне қарай (6.2) формуласынан әртүрлі итерациялық процестерді аламыз. Егер ягни -дан тәуелсіз болса, онда итерациялық процесті “стационар” , ал -дан тәуелді болса “стационар емес” деп, ал болса, онда қадамды итерациялық процесс дейміз.
Итерациялық процестер үшін (6.1) теңдеулер жүйесінің шешуі Х*-қозғалмайтын нүкте болып табылады, яғни деп алсақ, онда (6.2) формула бойынша тапқан келесі векторымыз да болады.
Керісінше -қозғалмайтын нүкте болатын кез-келген мынандай итерациялық процесті
, (6.3)
мұнда -матрицалар тізбегі, -векторлар тізбегі, (6.2) түрінде жазуға болады.
Шынында да үшін
болғандықтан
мұнда .
Енді итерациялық процесстің жинақталуын қарастыру үшін (6.2) теңдігін векторынан алып тастайық, ягни
. (6.4)
Осы формуланы былайша жазсақ
онда -векторлар тізбегінің векторына жинақталуы үшін
матрицасының нөлге ұмтылуы қажетті және жеткілікті, ал матрицасының кез-келген нормасының нөлге ұмтылуы жеткілікті шарт екенің көреміз.Ал -векторлар тізбегінің векторына ұмтылу жылдамдығы матрицасына байланысты болғандықтан, оны тандау әртүрлі итерациялық әдістерге әкеледі.
Енді итерациялық әдістердің жинақтылығын қалай сипаттауға болады, немесе қай уақытта итерацмялық процессті тоқтату керектігіне тоқталайық. Теңдеулер жүйесінің дәлдігін қателік векторы арқылы бағалау мүмкін болмағандықтан көбіне -ауытқу векторы арқылы бағалайды. Шынында да, егер (C-const) болса, онда
.
Осыдан болғанда .
екенің көреміз. Ал
(6.6)
болғандықтан
.
Осыдан мынандай тұжырым жасауға болады: егер (C-const) яғни (6.1) орнықты болса, онда шегінен шегі келіп шығады. Ал матрицасының шарттылығы нашар болған жағдайда бұл тұжырым орындалмауы да мүмкін.
Сондықтан итерациялық процессті , алдын-ала берілген дәлдік - бойынша теңсіздігі орындалған кезде тоқтату ертерек болуы мүмкін.
Енді (6.4) теңдігін былайша жазайық: .
Осыдан .
Бұл теңсіздік ,(6.1) теңдеулер жүйесі орнықты болған жағдайда, итерациялық процесті
(6.7)
теңсіздігі орындалғанда тоқтатуға болатындығын көрсетеді.
Егер А оң анықталған симметриялы матрица болған жағдайда қателік вектордың шамасын қателік функциясы деп аталатын
функцияның шамасымен де анықтауға болады. Шынында да А оң анықталған болғандықтан және болғанда .
Сонымен қатар
болғандықтан, оның мәнін табу мүмкін емес. Бірақ
(6.8)
функционалының қателік функциясынан айырмашылығы тұрақты сан болғандықтан және екеуі де өздерінің ең кіші мәндеріне болғанда ие болатындығын ескерсек, Ф(х) функционалына минимум мән беретін векторды жүйенің шешуі деп қарастыруға болады. Сондықтан
(6.9)
болған жағдайда да итерациялық процесті тоқтатуға болады.
Енді итерациялық процестердің қарапайым түрлерімен танысайық.
Біртіндеп жуықтау әдісі.
Берілген
теңдеулер жүйесін былайша жазайық
, (6.10)
мұнда .
Енді (6.10) теңдеулер жүйесін шешу үшін мынандай итерациялық әдісті пайдаланайық
, (6.11)
мұнда -өзіміздің еркіміз бойынша алынатын кез-келген бастапқы вектор. (6.11) итерациялық әдісін (6.2) формуласы түрінде жазсақ
,
онда екенің және бұл стационар итерациялық әдіс екенің көреміз.
Бұл әдістің жинақтылығын анықтау үшін, математикалық индукцияны қолдану арқылы, (6.11) теңдігінен алынған
(6.12)
және теңдігінен (6.11) теңдігін алып тастағанда алынған
(6.13)
формулаларын қарастырайық.
Енді осы формулаларды қолдана отырып әдістің жинақтылығы туралы бірнеше теоремаларды берейік.
6.1-теорема: Кез-келген -де біртіндеп жуықтау әдісінің жинақты болуы үшін матрицасының меншікті мәндерінің модулі бірден кіші болуы қажетті және жеткілікті.
Бұл теореманың шартын тексеру қиын болғандықтан, шарттары оңай анықталатын теоремаларды қарастырайық.
6.2-теорема: Біртіндеп жуықтау әдісі жинақталуы үшін болуы жеткілікті.