книги / Численные методы. Ч. 4
.pdfАлгоритм решения
Рассмотрим процедуру решения системы линейных алгебраических урав нений методом Гаусса в общем случае.
Система уравнений А х - f в развернутой (компонентной) записи пред ставляется в виде
ап х, +аа -х2+а,} х, +...+ аы ■х„=/„
а21 ■х1+а22-х2+агз-х3 +...+ а2т-х„= f 2,
( и )
°31 •*! + а 12 ' Х2 + °33 ‘ *3 + - + ^ п , - Х т = / з ,
°ml *1 + а ш2 - Х2 + ° т} - Х3 + - + а я т -Хт = / „ .
Пусть аи * 0 . Тогда первое уравнение системы (1.1) можно поделить на этот коэффициент,
1*1 + С \ 2 ' Х2 +Св-*3 +--.+q m' хт =^1 = — ,
а \\
где с12= а12/а м>с1з = а\ъ! а\\>•••» с\т - а\т/ап тема (1.1) преобразуется к виду
х2 + С13 • Xз +■■■+Сы-Х*
■Хг+а» *3
Л»)
•*з
• С помощью этого уравнения сис
II
II |
= / з '-*2Г.У|> |
|
|
м |
-«31-Тр |
.0) •*2+ а»3 *з +•••+«- • * . = л ои |
~ а/п1 ’ У\• |
||||
|
аи |
|
' с1;’i,j = 2, т . |
В полученной |
|
можно выделить подсистему m -1 |
линейных уравнений с т -1 неизвестными |
||||
величинами: |
|
|
|
|
|
°22 |
*2 + “23 |
*3 |
“2/я |
Хт~ J2 |
» |
* к |
-Х2 + а И |
*3 + - + О з 2 |
= / з (1). |
а ^ х 1+а ^ .х 3+...+ а ^ . хт= й '\
Пусть теперь а^ * 0. Поделим первое уравнение новой системы на этот
коэффициент: |
|
1 ■хг +си -х3+... +с2т-хя = у1 = Д ^ = ^ 2 |
3<|> — |
°22 |
а322. |
где с23 =с$Цс$1 , си = С2г} 1 ап > с2ш= аг11а 72 С помощью этого соотноше ния уравнения системы (1.1) преобразуются к виду
'1 • X, + с12 • х2 + с13 • х3 +...+ с1ш-xm= yv
О • X, +1 • х2 + с23 • х3 +...+ с2т • хт = у 2,
0-х, +0-х2 +а£> -Ху+...+а'£ ■хт= f { v = f<l>- а £ ■у 2,
0-х, + 0 x 2+ a% -x} +...+ a(JJ,-х„ = f 2>= f l - а™ -у2. |
|
|
Здесь обозначено: afj* = a fj - |
-c2j , i , j = 3,m. |
|
В результате преобразований получена подсистема т - 2 |
уравнений с |
|
т - 2 неизвестными: |
|
|
a<l>-xi+ -+ aZ > .xm=f?>, |
|
|
°% -x> + -+ °l2- x m= f i 2) |
|
|
Предполагая, что в первом уравнении последней системы |
* 0 , делим |
|
это уравнение на этот коэффициент: |
|
|
/ (2)
1*х3 +.. •+С2т*Хт = У} —"(2)" ‘ азз
Снова выполняются операции по понижению порядка системы алгебраи ческих уравнений, и так далее, до тех пор, пока вся система уравнений (1.1) не будет преобразована к виду
1 • Xj + с12 • х2 + с13 • х3 + |
...+ с Хт - х т - у |
Ь х 2+с2з-х3+... |
+ с2т-хт = У2 > |
1 *х3 +... |
+ с3т -хт = у 2, |
Процедура получения матрицы такого вида носит название «прямого хо да» метода Гаусса. Очевидным условием для успешного выполнения прямого
хода является |
* 0, у = 1,/я . |
«Обратный ход» метода позволяет определить искомые величины:
'Хт=Ут'
Х т- 1 — У я»—1 “ ^т-\тХ т '
, Х т -2 ~ Ут-1 ~ ^т-2тХ т ~~ ^т -2т -\Х т-\ >
= У1~ Z c'1клк- к=2
Вычитание строк в методе Гаусса (образование линейных комбинаций уравнений) не изменяет значения определителя матрицы. В результате выпол нения всех преобразований метода Гаусса определитель исходной матрицы может быть вычислен с использованием формулы
det(^) = n ^ - ' ) У=1
Таким образом, сохраняя значения коэффициентов, расположенных после преобразования уравнений на главной диагонали (до операции деления коэф фициентов строки на первый ненулевой элемент), можно вычислить определи тель исходной матрицы.
Пусть a pqi p,q = l,m - коэффициенты обратной матрицы А~х Согласно
|
т |
h |
/ - j |
определению, |
^ aik 'a kj = , где 8/у =< |
- символ Кронекера. Теперь |
|
q-й столбец |
к=1 |
I0’1 ф J |
|
обратной матрицы можно рассматривать как результат ре |
|||
шения системы линейных алгебраических уравнений вида |
|||
|
а 1(ч) |
|
|
|
«V»; |
*—< |
|
|
|
||
|
а р(я) |
|
\<п |
|
^ т(ч), |
|
Р(т) t |
ТаКим образом, для нахождения обратной матрицы необходимо решить т систем Линейных алгебраических уравнений с правыми частями, определен ными социальным образом. При этом матрицу коэффициентов следует преоб разовать лишь один раз, но одновременно преобразовывать т правых частей всех сиСтем уравнений.
Выполнение расчетов
Пусть система, состоящая из 15 линейных алгебраических уравнений, за дана с помощью табл. 1.1. Решение этой системы уравнений с помощью вы числительной программы, реализующей алгоритм метода Гаусса, приведено в табл. 1.2.
Компоненты обратной матрицы А~1 приведены в табл. 1.3. Детерминант матрицы А равен 5,71183*1017, детерминант матрицы А~] равен 1,75075-КГ18
Для проверки правильности определения обратной матрицы А~1 целесо образно вычислить произведение матриц А и А~1, которое по определению должно давать единичную матрицу:
1 |
0 |
0 |
О |
0 |
1 |
0 |
О |
0 |
0 |
1 |
О . |
0 |
0 |
0 |
1 |
Результат перемножения матриц А и А 1показан в табл. 1.4. Кроме того,
det(^ • А~х)= det(i4)det(>4"1)= det(£) = 1.
Перемножение приведенных выше значений определителей матриц А и
А"1дает
d e t^ d e t^ " 1)= 5,71183 • 1017 • 1,75075-Ю"18 =0,999999.
Это свидетельствует о хорошей точности проведенных вычислений. Поскольку точное решение поставленной задачи - вектор х - неизвестно,
определить погрешность полученного решения 8х = х - х невозможно. В пред положении, что величины А и / введены в память компьютера без погрешно стей, вычислим невязку т| = Ах - / , получаемую в результате подстановки в
исходную задачу приближенного решения х вместо точного решения х. Вы
численные значения невязки |
приведены в табл. 1.2. |
Учитывая далее, что |
|
Т| = Ах —f |
- Ах - Ах - А(х - дс)= АЪх, |
можно оценить погрешность получаемого решения с помощью выражения
8X = / T V ||8x||< |^'1||ri||,
где в качестве норм вектора и матрицы берутся, например,
н - т » ы , |
<|2 ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 1.1 |
|
|
|
Коэффициенты матрицы а,у и правая частьf |
системы линейных алгебраических уравнений A x = f |
|
|||||||||||||
|
|
у = 1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
f |
/= |
1 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
205,8 |
|
2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
222,8 |
|
3 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
153,2 |
|
4 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
209,7 |
|
5 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
119,0 |
|
6 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
126,9 |
|
7 |
0,0 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
45,5 |
|
8 |
0,0 |
0,0 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
0,0 |
0,0 |
100,4 |
|
9 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
0,0 |
-56,7 |
|
10 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
0,0 |
-108,9 |
|
11 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
0,0 |
66,2 |
|
12 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
0,6 |
61,3 |
|
13 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
1,2 |
124,5 |
|
14 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
3,5 |
-18,1 |
|
15 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,6 |
1,2 |
3,5 |
16,0 |
-15,7 |
|
|
|
В соответствии |
с |
выбранным |
||
|
Решение системы уравнений |
определением норм (1.2) для рас |
|||||
|
сматриваемой |
системы |
уравнений |
||||
|
и невязка уравнения |
получены значения: |
|
|
|||
|
на полученном решении |
|л|| = 1,4925-10-5, |
|||||
У |
Решение х, |
Невязка Л/ |
|||||
|
|
|
|
||||
1 |
10,0 |
1,09673-10-6 |
|
И | = 26,6, |
|
||
2 |
10,000001 |
1,34706-10‘5 |
||/Г,|| = 0,10373. |
||||
3 |
4,0000005 |
8,20160-10-6 |
|||||
4 |
9,999999 |
-1,08659-Ю"5 |
Следовательно, |
|
погрешность |
||
5 |
3,0 |
6,84261-10"6 |
найденного решения не превышает |
||||
6 |
6,000001 |
1,49250-10"5 |
||5х|| < 1,5482-1 O'* |
||||
7 |
—1,490116-10-7 |
9,53674-10'7 |
|||||
|
|
|
|
||||
8 |
7,0000005 |
1,01328-Ю-5 |
Числом обусловленности MAi ха |
||||
9 |
-4,000001 |
—1,16348-Ю-5 |
рактеризующим устойчивость реше |
||||
10 |
-8,0 |
3,81470-10"* |
ния системы |
линейных |
алгебраиче |
||
И |
5,0 |
6,79493-10-6 |
ских уравнений по отношению к от |
||||
12 |
2,0000005 |
4,33922-10-* |
носительному |
возмущению правой |
|||
13 |
8,0 |
-3,24249-10"6 |
части 5/, называется величина, опре |
||||
14 |
-3,0 |
-4,17233-Ю"7 |
деляемая выражением |
|
|||
15 |
-1,0000001 |
-1,90735-10-* |
Ч |
, = М |
' 1 4 |
||
|
|
|
Для рассматриваемой задачи, с
учетом полученных значений, МА= 2,7592.
Выводы
1.Разработана программа решения системы линейных алгебраических уравнений методом Гаусса.
2.С помощью этой программы найдено решение заданной системы урав нений (см. табл. 1.2).
3.Построена обратная матрица (см. табл. 1.3).
4. Вычислены определители исходной det(.4) — 5,71183*10 и обратной detOT1) = 1,75075-10-18 матриц.
5.Определено значение числа обусловленности МА = 2,7592. Поскольку это значение близко к 1, решение системы уравнений является устойчивым по отношению к возмущению правой части.
6.Погрешность решения системы уравнений не превышает 1,548*10
7.Для решения заданной системы линейных алгебраических уравнений
методом Гаусса на компьютере с процессором Intel® Pentium® 4 (тактовая час тота 2,2 ГГц, объем оперативной памяти 512 Мбайт) требуется 1,3*10 5 с.
PNRPU
0,0657 |
-0,0139 |
-0,0017 |
-0,0012 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
-0,0139 |
0,0687 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
-0,0017 |
-0,0136 |
0,0687 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
-0,0012 |
-0,0014 |
-0,0136 |
0,0688 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0014 |
-0,0014 |
-0,0136 |
0,0688 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0014 |
-0,0014 |
-0,0136 |
0,0688 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0014 |
-0,0014 |
-0,0136 |
0,0688 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0014 |
-0,0014 |
-0,0136 |
0,0688 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0014 |
-0,0014 |
-0,0136 |
0,0688 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0014 |
-0,0014 |
-0,0136 |
0,0688 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0014 |
-0,0014 |
-0,0136 |
0,0688 |
-0,0136 |
-0,0014 |
-0,0014 |
0,0009 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0014 |
-0,0014 |
-0,0136 |
0,0688 |
-0,0136 |
-0,0014 |
-0,0012 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 -0,0014 -0,0014 -0,0136 |
0,0687 |
-0,0136 |
-0,0017 |
|||
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0014 |
-0,0014 |
-0,0136 |
0,0687 |
-0,0139 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0 |
0,0009 |
-0,0012 |
-0,0017 |
-0,0139 |
0,0657 |
1,0 |
-92-10* |
3,0 КГ9 |
1,410"'° |
4,2-Ю"10 |
1,1-Ю"10-2,2-10"" -1,2-10'" -2,410'12 |
1,710'12 2,710"13 3,010Г13 1,810т14 8,710'15 -7,310т16 |
||||||||
-6,410* |
1,0 |
2,010"* |
1Д10"9 |
2,3-10"9 -6,910"10-6,810"" -6,910"" |
1,010'" |
1,910'12 2,310"12 -22-КГ12 1,2-10"'3 -8,410т’4 62-10Г’5 |
||||||||
-5,4 ЮГ* |
2,610"* |
1,0 |
-1,610"* |
1,910"9 |
-2,010"9 |
1,1-10"9 -1,0 Ю"10 |
2,310'" |
3,410'" 7,010"12-5,1-10"13 -4,510~13 |
4,810Г13 |
3,7- 10Г14 |
||||
- 1,110 * |
4,5 10* |
6,1-10* |
1,0 |
7,610"9 |
-2 ,5 10"9 -2,010"9 |
4,810"" |
9,410"” |
3,510"” -7,610"13 |
8,410~12 |
1,510"'2 |
1,0 ИГ12-9,610Г13 |
|||
-4,6-10-10 |
3,710* |
зд-кг9 |
-2,410"* |
1,0 |
-1,510т9 |
-7,7-10"9 -5,910Г10 |
9,810"’° |
3,010"’° -4,610"” -5,810'" |
1,410"” |
1,610"’2 -4,010Г14 |
||||
-2,410'10 |
3,01ОТ9 |
1,510"9 |
-5,010"9 |
-3,010"* |
1,0 |
2,610"9 -4,9109 |
-2,010"9 |
7,910"'° 2,910"'° 23-10"" -9,310'" |
1210"” |
1,510~12 |
||||
3,3-1(Г" |
2,410"10 |
8,3- 10г'° |
-6,710Г9 |
—4,510"9 |
7310Г* |
1,0 |
-4,7 Ю"9 |
-1,410"9 -2,010"9 |
1,310"9 -92-10"" -1,1-10"” -82-10"" |
9,1-10"’2 |
||||
1,0 КГ12-3,410"'2 -8,810'" |
-9,310"'° |
- 13-ю"9 -3,4109 |
5.510"9 |
1,0 |
-12-10* |
—4,2-10"9 -2,0109 -5,1-10"'° |
1,2-10"'° -4,610"” |
3,310г" |
||||||
2,4 КГ" -6,510"" -4,010"'2 |
4,610"" |
- п о ю 9 |
4 7 - Ю"9 |
33-109 |
1,0 ИГ* |
1,0 |
1,2-10* -1,610"* |
1,610"9 |
2,310"'° |
2,91СГ10-8,810"” |
||||
7,01(Г'2 1,010"" -72-10"” -7,71<Г" |
2,510"'° 2,81(Г10 -1 ,4 109 |
-9 ,9 10"9 |
2,510* |
1,0 |
9,9109 |
-53- Ю"9 -1,310"9 2,3410"’° |
12-ltr'0 |
|||||||
-1,1-1<Г'3 6,1-10"12 |
1,010г" —1,5 Ю"10 |
5,310"12 |
2,0 Ют'0 |
1,6109 |
1,410Г10 -4,910'9 |
3,310"* |
1,0 |
2,410"* |
-1,010* |
2310* |
4,1-10"’° |
|||
4,01(Г|Э-1,1-1ОТ12 |
2,71СГ12-2,010"" -1,510"" -1,0 Ю'30 |
3,010Г10 |
1,810"’° |
-3,2-10"9 |
-2,810"9 |
62-109 |
1,0 |
6,310* |
-1,410* |
7,610"” |
||||
-1,9-Ю-'3 5,710"'5-8,710"'3 |
3,М0"12 |
8,410"12-1,3 юг" |
4,410‘12 |
1,2-10"’° |
8,410"'° |
- З. ИО"9 -1,610* 2,810* 1,0 |
13Ю"* |
-7210* |
||||||
-3,Я0"'4 -3,1-10"'4 -3,2-1ОТ13 -1,Ы<Г12 |
5,510"12-5,910г'2 -1,310т" 1,810"12 |
3,2- Ю"10 |
8,510"10 -4,710 9 |
5,810"” |
1,410"* |
1,0 |
1,1-10* |
|||||||
-2,8-1(Г14 |
2,810"'4 -1,710'13 |
0,0 |
9,010"13 9,0910"'3 |
0,0 |
0,0 |
-72-10"12 |
23-10"’° |
1,810* |
0,0 |
-1,1-10* |
1,410"* |
1,0 |
7.7.2.Метод квадратного корня
Задание. Для системы линейных алгебраических уравнений А х - f с сим метричной матрицей А:
-разработать вычислительную программу, реализующую метод квад ратного корня;
-найти решение этой системы уравнений;
-построить обратную матрицу А~1;
- вычислить определители det(i4) и d et^"1);
-оценить погрешность нахождения решения системы уравнений;
-подсчитать значение числа обусловленности МА матрицы А;
-оценить быстродействие вычислительной программы.
Алгоритм решения
Метод квадратного корня предназначен для решения систем линейных ал гебраических уравнений вида А х - f z симметричной матрицей коэффициен
тов atj = ajh i,j = 1,m . Метод основан на разложении матрицы коэффициентов
А в произведение:
A = S TDS,
где S - верхняя треугольная матрица с положительными значениями на главной диагонали; D - диагональная матрица со значениями +1 или -1. Соотношения
для вычисления диагональных значений матриц SnD: |
|
||||
Г |
i-l |
\ |
„ |
Jz}у-' |
__ |
djj = sign |
ai j - L sljdkk |
> aj j = s 2jdjj + H sh dM' |
J = ^ m > |
||
|
|
) |
|
A=lЫ1 |
|
|
|
2jj |
^LiSkjdkk , |
j = \ , m . |
|
|
|
|
k=\ |
|
|
Наддиагональные элементы матрицы S определяются по формуле |
|||||
|
a iу ” S |
s ki^kks kj |
___ |
|
|
|
SU = ------ Z Z --------- > |
= |
|
||
|
|
*iiu ii |
|
|
|
Если ввести обозначения |
у = Sx9 Z = D y, то решение системы линейных |
алгебраических уравнений методом квадратного корня S TDSx = / можно рас сматривать как последовательность трех процессов: 1) решения системы урав нений S Tz = / , то есть вычисление решения z системы уравнений с нижней
треугольной матрицей S r ; 2) решения системы уравнений Dy = z , то есть вы числения решения у системы уравнений с диагональной матрицей D; 3) реше
ния системы уравнений Sx = у , то есть определения искомого решения из сис темы уравнений с верхней треугольной матрицей S.
Выполнение расчетов
Пусть система 15 линейных алгебраических уравнений задана, как и ранее, с помощью табл. 1.1. Решение этой системы уравнений с помощью вычисли тельной программы, реализующей алгоритм метода квадратного корня, приве дено в табл. 1.5.
Компоненты вектора невязки т|, получаемой в результате подстановки
приближенного |
решения х вместо точного решения х в исходную |
задачу, |
||
^ ^ Лх |
приведены в табл. 1.5. Определенная с помощью выражения (1.2) |
|||
Норма невязки |
|
|
|
|
|
|
||л|| = 1,8573-10-5 |
|
|
Компоненты обратной матрицы |
приведены в табл. 1.6. В соответствии |
|||
с выбранным |
определением норм |
(1.2) получены значения |
норм |
||л|| = 1,8573 10"5, ||у4_1|| = 0,10373. Это позволяет оценить погрешность найден
ного решения, которая не превышает ||5х|| ^ 1,9266 • 10"6
Таблица 1.5
Решение системы уравнений и невязка уравнения
на полученном решении
/ |
Решение х, |
Невязка т|/ |
1 |
10,0 |
-6,9141410^ |
2 |
9,9999990 |
-1,85728-Ю-5 |
3 |
3,9999998 |
-5,2452 М О'7 |
4 |
10,0 |
3,90693-10-6 |
5 |
3,0000002 |
5,60284-10"6 |
6 |
5,9999995 |
-6,43730-10-6 |
7 |
0,0000003 |
4,19617-Ю"4 |
8 |
7,0 |
-2,98023-10-6 |
9 |
^1,0 |
2,86102-10-7 |
10 |
-8,0000010 |
-1,08719-10'5 |
11 |
5,0 |
7,70092-10-6 |
12 |
2,0000005 |
7,00951-10'* |
13 |
8,0000010 |
1,02043-10'5 |
14 |
-3,0000005 |
-4,64916-1 O'6 |
15 |
-1,0000002 |
-3,81470-Ю-6 |
Выводы
1.Разработана программа реше ния системы линейных алгебраиче ских уравнений с симметричной мат рицей методом квадратного корня.
2.С помощью этой программы
найдено решение заданной системы уравнений.
3.Построена обратная матрица (см. табл. 1.6).
4.Невязка решения системы
уравнений не превышает 1,8573-10"5 5. Для решения заданной систе мы линейных алгебраических урав
нений методом квадратного корня на компьютере с процессором Intel® Pentium® 4 (тактовая частота 2,2 ГГц, объем оперативной памяти 512 Мбайт) требуется 2,4-10-5 с.