Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
LibreTexts - Ukrayinska

8: Метод Гаусса-Зайделя

Цілі навчання

Прочитавши цю главу, ви повинні мати можливість:

  1. вирішувати множину рівнянь методом Гауса-Зайделя,
  2. визнати переваги і підводні камені методу Гаусса-Зайделя, і
  3. визначити, за яких умов завжди сходиться метод Гаусса-Зайделя.

Навіщо потрібен інший метод для вирішення безлічі одночасних лінійних рівнянь?

У певних випадках, наприклад, коли система рівнянь велика, ітераційні методи розв'язання рівнянь більш вигідні. Методи ліквідації, такі як гаусова елімінація, схильні до великих похибок округлення для великого набору рівнянь. Ітераційні методи, такі як метод Гаусса-Зайделя, дають користувачеві контроль над помилкою округлення. Крім того, якщо фізика задачі добре відома, початкові здогадки, необхідні в ітераційних методах, можуть бути зроблені більш розумно, що призводять до більш швидкої збіжності.

Який алгоритм для методу Гауса-Зайделя? З огляду на загальний набірn рівнянь іn невідомих, ми маємо

a11x1+a12x2+a13x3+...+a1nxn=c1

a21x1+a22x2+a23x3+...+a2nxn=c2

                   

an1x1+an2x2+an3x3+...+annxn=cn

Якщо діагональні елементи ненульові, кожне рівняння переписується на відповідне невідоме, тобто перше рівняння переписується зx1 лівого боку, друге рівняння переписується зx2 лівої сторони і так далі наступним чином

x1=c1a12x2a13x3a1nxna11

x2=c2a21x1a23x3a2nxna22

xn1=cn1an1,1x1an1,2x2an1,n2xn2an1,nxnan1,n1

xn=cnan1x1an2x2an,n1xn1ann

Ці рівняння можуть бути переписані у вигляді підсумовування як

x1=c1nj=1j1a1jxja11

x2=c2nj=1j2a2jxja22

xn1=cn1nj=1jn1an1,jxjan1,n1

xn=cnnj=1jnanjxjann

Значить для будь-якого рядуi,

xi=cinj=1jiaijxjaii,i=1,2,,n.

Тепер, щоб знайтиxi 's, один припускає початкове припущення дляxi' s, а потім використовує переписані рівняння для обчислення нових оцінок. Пам'ятайте, що завжди використовуються найновіші оцінки для обчислення наступних оцінок,xi. Наприкінці кожної ітерації обчислюється абсолютна відносна приблизна похибка для кожноїxi як

|a|i=|xnewixoldixnewi|×100

деxnewi - нещодавно отримане значенняxi, іxoldi є попереднім значеннямxi.

Коли абсолютна відносна приблизна похибка для кожногоxi менше, ніж попередньо заданий допуск, ітерації зупиняються.

Приклад 1

Швидкість вгору ракети наведена в три різні часи в наступній таблиці.

Таблиця 1: Дані про швидкість проти часу.

Час,t(s) Швидкість,v(m/s)
\ (t (s)\) ">5 \ (v (м/с)\) ">106.8
\ (t (s)\) ">8 \ (v (м/с)\) ">177.2
\ (t (s)\) ">12 \ (v (м/с)\) ">279.2

Дані про швидкість наближаються поліномом як

v(t)=a1t2+a2t+a3,5t12

Знайти значення заa1, a2,and a3 допомогою методу Гаусса-Зайделя. Припустімо початкове припущення рішення як

[a1a2a3]=[125]

і провести дві ітерації.

Рішення

Многочлен проходить через три точки даних,(t1,v1),(t2,v2),and(t3,v3) звідки з наведеної вище таблиці.

t1=5,v1=106.8

t2=8,v2=177.2

t3=12,v3=279.2

Вимога, щоv(t)=a1t2+a2t+a3 проходить через три точки даних дає

v(t1)=v1=a1t21+a2t1+a3

v(t2)=v2=a1t22+a2t2+a3

v(t3)=v3=a1t23+a2t3+a3

Підстановка даних(t1,v1),(t2,v2),and(t3,v3) дає

a1(52)+a2(5)+a3=106.8

a1(82)+a2(8)+a3=177.2

a1(122)+a2(12)+a3=279.2

або

25a1+5a2+a3=106.8

64a1+8a2+a3=177.2

144a1+12a2+a3=279.2

Коефіцієнтиa1,a2,anda3 для вищевказаного виразу задаються

[25516481144121][a1a2a3]=[106.8177.2279.2]

Переписування рівнянь дає

a1=106.85a2a325

a2=177.264a1a38

a3=279.2144a112a21

Ітерація #1

Задано початкове припущення вектора розв'язку як

[a1a2a3]=[125]

отримуємо

a1=106.85(2)(5)25=3.6720

a2=177.264(3.6720)(5)8=7.8150

a3=279.2144(3.6720)12(7.8510)1=155.36

Абсолютна відносна приблизна похибка для кожногоxi тоді дорівнює

|a|1=|3.672013.6720|×100=72.76%

|a|2=|7.851027.8510|×100=125.47%

|a|3=|155.365155.36|×100=103.22%

В кінці першої ітерації оцінка вектора розв'язку дорівнює

[a1a2a3]=[3.67207.8510155.36]

і максимальна абсолютна відносна приблизна похибка дорівнює125.47.

Ітерація #2

Оцінка вектора розв'язку в кінці Ітерації #1 дорівнює

[a1a2a3]=[3.67207.8510155.36]

Тепер отримуємо

a1=106.85(7.8510)(155.36)25=12.056

a2=177.264(12.056)(155.36)8=54.882

a3=279.2144(12.056)12(54.882)1=798.34

Абсолютна відносна приблизна похибка для кожногоxi тоді дорівнює

|a|1=|12.0563.672012.056|×100=69.543%

|a|2=|54.882(7.8510)54.882|×100=85.695%

|a|3=|798.34(155.36)798.34|×100=80.540%

В кінці другої ітерації оцінка вектора розв'язку дорівнює

[a1a2a3]=[12.05654.882798.54]

і максимальна абсолютна відносна приблизна похибка дорівнює85.695%.

Проведення більшої кількості ітерацій дає наступні значення вектора розв'язку та відповідні абсолютні відносні наближені похибки.

Ітерація a1 |a|1% a2 |a|2% a3 |a|3%
1 \ (a_ {1}\) ">3.672 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">72.767 \ (a_ {2}\) ">—7.8510 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">125.47 \ (a_ {3}\) ">—155.36 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">103.22
2 \ (a_ {1}\) ">12.056 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">69.543 \ (a_ {2}\) ">—54.882 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">85.695 \ (a_ {3}\) ">—798.34 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">80.54
3 \ (a_ {1}\) ">47.182 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">74.447 \ (a_ {2}\) ">—255.51 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">78.521 \ (a_ {3}\) ">—3448.9 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">76.852
4 \ (a_ {1}\) ">193.33 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">75.595 \ (a_ {2}\) ">—1093.4 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">76.632 \ (a_ {3}\) ">—14440 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">76.116
5 \ (a_ {1}\) ">800.53 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">75.85 \ (a_ {2}\) ">—4577.2 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">76.112 \ (a_ {3}\) ">—60072 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">75.963
6 \ (a_ {1}\) ">3322.6 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">75.906 \ (a_ {2}\) ">—19049 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">75.972 \ (a_ {3}\) ">—249580 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">75.931

Як видно з наведеної вище таблиці, оцінки рішення не сходяться до істинного рішення

a1=0.29048

a2=19.690

a3=1.0857

Наведена вище система рівнянь, схоже, не сходиться. Чому?

Ну, підводним каменем більшості ітераційних методів є те, що вони можуть або не можуть сходитися. Однак рішення певного класу системи одночасних рівнянь завжди сходиться за допомогою методу Гаусса-Зайделя. Цей клас системи рівнянь - це де матриця коефіцієнтів[A] в[A][X]=[C] діагонально домінуюча, тобто

|aii|nj=1ji|aij| for all i

|aii|>nj=1ji|aij|for at least one i

Якщо система рівнянь має матрицю коефіцієнтів, яка не є домінуючою по діагоналі, вона може сходитися або не сходитися. На щастя, багато фізичних систем, що призводять до одночасних лінійних рівнянь, мають діагонально домінуючу матрицю коефіцієнтів, яка потім забезпечує збіжність для ітераційних методів, таких як метод Гаусса-Зайделя розв'язування одночасних лінійних рівнянь.

Приклад 2

Знайти розв'язок наступної системи рівнянь за допомогою методу Гаусса-Зайделя.

12x1+3x25x3=1

x1+5x2+3x3=28

3x1+7x2+13x3=76

Використовувати

[x1x2x3]=[101]

як початкове вгадати і провести дві ітерації.

Рішення 1

матриця коефіцієнтів

[A]=[12351533713]

є діагонально домінуючим як

|a11|=|12|=12|a12|+|a13|=|3|+|5|=8

|a22|=|5|=5|a21|+|a23|=|1|+|3|=4

|a33|=|13|=13|a31|+|a32|=|3|+|7|=10

і нерівність строго більше, ніж для хоча б одного ряду. Значить, розчин повинен сходитися за методом Гауса-Зайделя.

Переписуючи рівняння, отримуємо

x1=13x2+5x312

x2=28x13x35

x3=763x17x213

Припускаючи початкове припущення

[x1x2x3]=[101]

Ітерація #1

x1=13(0)+5(1)12=0.50000

x2=28(0.50000)3(1)5=4.9000

x3=763(0.50000)7(4.9000)13=3.0923

Абсолютна відносна приблизна похибка в кінці першої ітерації дорівнює

|a|1=|0.5000010.50000|×100=100.00%

|a|2=|4.900004.9000|×100=100.00%

|a|3=|3.092313.0923|×100=67.662%

Максимальна абсолютна відносна приблизна похибка дорівнює100.00%

Ітерація #2

x1=13(4.9000)+5(3.0923)12=0.14679

x2=28(0.14679)3(3.0923)5=3.7153

x3=763(0.14679)7(3.7153)13=3.8118

В кінці другої ітерації абсолютна відносна приблизна похибка дорівнює

|a|1=|0.146790.500000.14679|×100=240.61%

|a|2=|3.71534.90003.7153|×100=31.889%

|a|3=|3.81183.09233.8118|×100=18.874%

Максимальна абсолютна відносна приблизна похибка дорівнює240.61%. Це більше, ніж значення, отримане100.00% нами в першій ітерації. Рішення розходиться? Ні, оскільки ви проводите більше ітерацій, рішення сходиться наступним чином.

Ітерація a1 |a|1% a2 |a|2% a3 |a|3%
1 \ (a_ {1}\) ">0.5 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">100 \ (a_ {2}\) ">4.9 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">100 \ (a_ {3}\) ">3.0923 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">67.662
2 \ (a_ {1}\) ">0.14679 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">240.61 \ (a_ {2}\) ">3.7153 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">31.889 \ (a_ {3}\) ">3.8118 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">18.874
3 \ (a_ {1}\) ">0.74275 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">80.236 \ (a_ {2}\) ">3.1644 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">17.408 \ (a_ {3}\) ">3.9708 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">4.0064
4 \ (a_ {1}\) ">0.94675 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">21.546 \ (a_ {2}\) ">3.0281 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">4.4996 \ (a_ {3}\) ">3.9971 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">0.65772
5 \ (a_ {1}\) ">0.99177 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">4.5391 \ (a_ {2}\) ">3.0034 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">0.82499 \ (a_ {3}\) ">4.0001 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">0.074383
6 \ (a_ {1}\) ">0.99919 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">0.74307 \ (a_ {2}\) ">3.0001 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">0.10856 \ (a_ {3}\) ">4.0001 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">0.00101

Це близько до точного вектора розв'язку

[x1x2x3]=[134]

Приклад 3

Задано систему рівнянь

3x1+7x2+13x3=76

x1+5x2+3x3=28

12x1+3x25x3=1

знайти рішення за допомогою методу Гауса-Зайделя. Використовувати

[x1x2x3]=[101]

як початкова здогадка.

Рішення 1

Переписуючи рівняння, отримуємо

x1=767x213x33

x2=28x13x35

x3=112x13x25

Припускаючи початкове припущення

[x1x2x3]=[101]

наступні шість ітераційних значень наведені в таблиці нижче.

Ітерація a1 |a|1% a2 |a|2% a3 |a|3%
1 \ (a_ {1}\) ">21 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">95.238 \ (a_ {2}\) ">0.8 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">100 \ (a_ {3}\) ">50.68 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">98.027
2 \ (a_ {1}\) ">—196.15 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">110.71 \ (a_ {2}\) ">14.421 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">94.453 \ (a_ {3}\) ">—462.30 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">110.96
3 \ (a_ {1}\) ">1995 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">109.83 \ (a_ {2}\) ">—116.02 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">112.43 \ (a_ {3}\) ">4718.1 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">109.8
4 \ (a_ {1}\) ">—20149 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">109.9 \ (a_ {2}\) ">1204.6 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">109.63 \ (a_ {3}\) ">—47636 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">109.9
5 \ (a_ {1}\) ">2.0364×105 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">109.89 \ (a_ {2}\) ">—12140 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">109.92 \ (a_ {3}\) ">4.8144×105 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">109.89
6 \ (a_ {1}\) ">—2.0579×106 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {1}\%\) ">109.89 \ (a_ {2}\) ">1.2272×105 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {2}\%\) ">109.89 \ (a_ {3}\) ">—4.8653×106 \ (\ ліворуч |\ in_ {a}\ праворуч |_ {3}\%\) ">109.89

Ви можете бачити, що це рішення не сходиться і матриця коефіцієнтів не є домінуючою по діагоналі. матриця коефіцієнтів

[A]=[37131531235]

не є домінуючим по діагоналі, як

|a11|=|3|=3|a12|+|a13|=|7|+|13|=20

Значить, метод Гауса-Зайделя може сходитися, а може і не сходитися.

Однак це той самий набір рівнянь, що і попередній приклад, і який сходився. Єдина відмінність полягає в тому, що ми обмінялися першим і третім рівнянням один з одним, і це зробило матрицю коефіцієнтів не домінуючою по діагоналі.

Тому не виключено, що система рівнянь може бути зроблена діагонально домінуючою, якщо обмінюватися рівняннями один з одним. Однак це можливо не для всіх випадків. Наприклад, наступний набір рівнянь

x1+x2+x3=3

2x1+3x2+4x3=9

x1+7x2+x3=9

неможливо переписати, щоб зробити матрицю коефіцієнтів по діагоналі домінуючою.

Вікторина методу Гауса-Зайделя

Вікторина 1

Квадратна матриця[A]n×n є домінуючою по діагоналі, якщо

(А)|aii|nj=1ij|aij|,i=1,2,...,n

(B)|aii|nj=1ij|aij|,i=1,2,...,n і|aii|>nj=1ij|aij|, для будь-якогоi=1,2,...,n

(C)|aii|nj=1|aij|,i=1,2,...,n і|aii|>nj=1|aij|, для будь-якогоi=1,2,...,n

(D)|aii|nj=1|aij|,i=1,2,...,n

Вікторина 2

Використовуючи[x1,x2,x3]=[1,3,5] як початкову здогадку, значення[x1,x2,x3] після трьох ітерацій у методі Гаусса-Зайделя для

[12731512711][x1x2x3]=[256]

є

(А)[2.83331.43331.9727]

(Б)[1.49590.904640.84914]

(С)[0.906661.01151.0243]

(D)[1.21480.720600.82451]

Вікторина 3

Для того, щоб наступна система рівнянь,

2x1+7x211x3=6x1+2x2+x3=57x1+5x2+2x3=17

сходиться за допомогою методу Гаусса-Зайделя, можна переписати вищевказані рівняння наступним чином:

(А)[2711121752] [x1x2x3]=[6517]

(Б)[7521212711] [x1x2x3]=[1756]

(С)[7521212711] [x1x2x3]=[6517]

(D) Рівняння не можна переписати у формі для забезпечення збіжності.

Вікторина 4

Для[12731512711] [x1x2x3]=[2272] та використання[x1x2x3]=[121] як початкового припущення значення[x1x2x3] знаходять в кінці кожної ітерації як

Ітерація # x1 x2 x3
1 \ (x_ {1}\) ">0.41667 \ (x_ {2}\) ">1.1167 \ (x_ {3}\) ">0.96818
2 \ (x_ {1}\) ">0.93990 \ (x_ {2}\) ">1.0184 \ (x_ {3}\) ">1.0008
3 \ (x_ {1}\) ">0.98908 \ (x_ {2}\) ">1.0020 \ (x_ {3}\) ">0.99931
4 \ (x_ {1}\) ">0.99899 \ (x_ {2}\) ">1.0003 \ (x_ {3}\) ">1.0000

На яке перше число ітерації ви б довіряли принаймні 1 значній цифрі у своєму рішенні?

(А)1

(Б)2

(С)3

(D)4

Вікторина 5

Алгоритм розв'язання методу Гаусса-Зайделя[A][X]=[C] наведено наступним чином при використанніnmax ітерацій. Початкове значення[X] зберігається в[X].

(A) Суб

Зайдель(n,a,x,rhs,nmax)

Дляk=1 Комуnmax

Дляi=1 Комуn

Дляj=1 Комуn

Якщо (i<>j) Тоді

Сума = Сума +a(i,j)x(j)

кінець

Наступнийj

x(i)=(rhs(i)Sum)/a(i,i)

Наступнийi

Наступнийj

Кінець Sub

(Б) Суб

Зайдель(n,a,x,rhs,nmax)

Дляk=1 Комуnmax

Дляi=1 Комуn

Сума = 0

Дляj=1 Комуn

Якщо (i<>j) Тоді

Сума = Сума +a(i,j)x(j)

кінець

Наступнийj

x(i)=(rhs(i)Sum)/a(i,i)

Наступнийi

Наступнийk

Кінець Sub

(C) Суб

Зайдель(n,a,x,rhs,nmax)

Дляk=1 Комуnmax

Дляi=1 Комуn

Сума = 0

Дляj=1 Комуn

Сума = Сума +a(i,j)x(j)

Наступнийj

x(i)=(rhs(i)Sum)/a(i,i)

Наступнийi

Наступнийk

Кінець Sub

(D) Суб

Зайдель(n,a,x,rhs,nmax)

Дляk=1 Кому {nmax} $

Дляi=1 Комуn

Сума = 0

Дляj=1 Комуn

Якщо (i<>j) Тоді

Сума = Сума +a(i,j)x(j)

кінець

Наступнийj

x(i)=(rhs(i)Sum)/a(i,i)

Наступнийi

Наступнийk

Кінець Sub

Вікторина 6

Термістори вимірюють температуру, мають нелінійний вихід і цінуються для обмеженого діапазону. Отже, коли виготовляється термістор, виробник постачає криву опору проти температури. Точне уявлення кривої, як правило, задається

1T=a0+a1ln(R)+a2{ln(R)}2+a3{ln(R)}3

деT температура в Кельвіні,R - це опір в Омах, іa0,a1,a2,a3 є константами калібрувальної кривої. З огляду на наступне для термістора

R T
\ (R\) ">Ом \ (Т\) ">C
\ (R\) ">

1101.0

911.3

636.0

451.1

\ (Т\) ">

25.113

30.131

40.120

50.128

значення температури вC для вимірюваного опору900 Ом найбільш майже дорівнює

(А)30.002

(Б)30.473

(С)31.272

(D)31.445

Вправа методу Гауса-Зайделя

Вправа 1

У системі рівнянь[A][X]=[C], якщо[A] діагонально домінантний, то метод Гауссейдала-Зайделя

  1. завжди сходиться
  2. може сходитися або не сходитися
  3. завжди розходиться
Відповідь

A

Вправа 2

У системі[A][X]=[C] рівнянь якщо не[A] є домінуючою по діагоналі, то метод Гаусса-Зайделя

  1. Завжди сходиться
  2. Може сходитися або не сходитися
  3. Завжди розходиться.
Відповідь

Б

Вправа 3

У системі рівнянь[A][X]=[C], якщо[A] вона не домінуюча по діагоналі, система рівнянь завжди може бути переписана, щоб зробити її домінуючою по діагоналі.

  1. Правда
  2. Помилковий
Відповідь

Б

Вправа 4

Розв'яжіть наступну систему рівнянь методом Гауса-Зайделя

12x1+7x2+3x3=2x1+5x2+x3=52x1+7x211x3=6

Проведіть 3 ітерації, обчислите максимальну абсолютну відносну приблизну похибку в кінці кожної ітерації і виберіть[x1x2x3]=[135] як початкову здогадку.

Відповідь

[x1x2x3]=[0.906661.01151.0243][|a|1|a|2|a|3]=[65.001%10.564%17.099%][|a|1|a|2|a|3]=[65.001%10.564%17.099%]

Вправа 5

Розв'яжіть наступну систему рівнянь методом Гауса-Зайделя

12x1+7x2+3x3=2
x1+5x2+x3=5
2x1+7x211x3=6
Проведіть 3 ітерації, обчислите максимальну абсолютну відносну приблизну похибку в кінці кожної ітерації та виберіть[x1x2x3]=[135] як початкову здогадку.

Відповідь

[x1x2x3]=[0.906661.01151.0243]

Вправа 6

Вирішіть наступну систему рівнянь методом Гауса-Зайделя
x1+5x2+x3=5
12x1+7x2+3x3=2
2x1+7x211x3=6
Проведіть 3 ітерації, обчислите максимальну абсолютну відносну наближену похибку в кінці кожної ітерації, і вибрати[x1x2x3]=[135] в якості початкової здогадки.

Відповідь

[x1x2x3]=[1163.71947.61027.2]

[|a|1|a|2|a|3]=[89.156%89.139%89.183%]