«Сандық әдістер» пәнінің оқу-әдістемелік кешені
Осы сияқты екінші және төртінші теңдіктерден
жүктеу/скачать 1.3 Mb.
|
- Бұл бет үшін навигация:
- Сонымен
Осы сияқты екінші және төртінші теңдіктерденN=(D-CA-1B)-1, L=A-1BN. Бұл әдіс қарастырылған кері матрицалар бар болғанда ғана іске асатындығын естен шығармаған жөн. Жоғарыдағы формулаларды былайша өзгертуге болады: N=(D-CA-1B)-1, M=-NCA-1, L=-A-1BN, K=A-1-A-1BM. Немесе K=(A-BD-1C)-1, L=KBD-1, M=-D-1CK, N=D-1-D-1CL. Бұл формулаларда өлшемдері p және q екі матрицаның кері матрицасын табу жеткілікті. Соңғы формулалардан, матрицаның бас диагнолындағы клеткалардың кері матрицасы оңай табылатын жағдайда, клеткаға бөлу әдісі тиімді екенің көреміз. Көмкеру әдісі. Берілген А матрицасын мына түрде жазайық , мұнда . Егер белгілі десек, онда матрицасын мына түрде іздестіреміз: , мұнда біз табуға тиісті pn-1-матрица, qn-жатық бағана, 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), матрицаларының кері матрицасын көмкеру әдісі арқылы біртіндеп табу арқылы іске асырылады. Есептеу схемасы мынандай: ( 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) Осы формуланы былайша жазсақ онда -векторлар тізбегінің векторына жинақталуы үшін матрицасының нөлге ұмтылуы қажетті және жеткілікті, ал матрицасының кез-келген нормасының нөлге ұмтылуы жеткілікті шарт екенің көреміз.Ал -векторлар тізбегінің векторына ұмтылу жылдамдығы матрицасына байланысты болғандықтан, оны тандау әртүрлі итерациялық әдістерге әкеледі. жүктеу/скачать 1.3 Mb. Достарыңызбен бөлісу:
|