8: Метод Гаусса-Зайделя
Прочитавши цю главу, ви повинні мати можливість:
- вирішувати множину рівнянь методом Гауса-Зайделя,
- визнати переваги і підводні камені методу Гаусса-Зайделя, і
- визначити, за яких умов завжди сходиться метод Гаусса-Зайделя.
Навіщо потрібен інший метод для вирішення безлічі одночасних лінійних рівнянь?
У певних випадках, наприклад, коли система рівнянь велика, ітераційні методи розв'язання рівнянь більш вигідні. Методи ліквідації, такі як гаусова елімінація, схильні до великих похибок округлення для великого набору рівнянь. Ітераційні методи, такі як метод Гаусса-Зайделя, дають користувачеві контроль над помилкою округлення. Крім того, якщо фізика задачі добре відома, початкові здогадки, необхідні в ітераційних методах, можуть бути зроблені більш розумно, що призводять до більш швидкої збіжності.
Який алгоритм для методу Гауса-Зайделя? З огляду на загальний набірn рівнянь іn невідомих, ми маємо
a11x1+a12x2+a13x3+...+a1nxn=c1
a21x1+a22x2+a23x3+...+a2nxn=c2
⋮ ⋮
an1x1+an2x2+an3x3+...+annxn=cn
Якщо діагональні елементи ненульові, кожне рівняння переписується на відповідне невідоме, тобто перше рівняння переписується зx1 лівого боку, друге рівняння переписується зx2 лівої сторони і так далі наступним чином
x1=c1−a12x2−a13x3……−a1nxna11
x2=c2−a21x1−a23x3……−a2nxna22
⋮
⋮
xn−1=cn−1−an−1,1x1−an−1,2x2……−an−1,n−2xn−2−an−1,nxnan−1,n−1
xn=cn−an1x1−an2x2−……−an,n−1xn−1ann
Ці рівняння можуть бути переписані у вигляді підсумовування як
x1=c1−∑nj=1j≠1a1jxja11
x2=c2−∑nj=1j≠2a2jxja22
⋮
xn−1=cn−1−∑nj=1j≠n−1an−1,jxjan−1,n−1
xn=cn−∑nj=1j≠nanjxjann
Значить для будь-якого рядуi,
xi=ci−∑nj=1j≠iaijxjaii,i=1,2,…,n.
Тепер, щоб знайтиxi 's, один припускає початкове припущення дляxi' s, а потім використовує переписані рівняння для обчислення нових оцінок. Пам'ятайте, що завжди використовуються найновіші оцінки для обчислення наступних оцінок,xi. Наприкінці кожної ітерації обчислюється абсолютна відносна приблизна похибка для кожноїxi як
|∈a|i=|xnewi−xoldixnewi|×100
деxnewi - нещодавно отримане значенняxi, іxoldi є попереднім значеннямxi.
Коли абсолютна відносна приблизна похибка для кожногоxi менше, ніж попередньо заданий допуск, ітерації зупиняються.
Швидкість вгору ракети наведена в три різні часи в наступній таблиці.
Час,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,5≤t≤12
Знайти значення за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.8−5a2−a325
a2=177.2−64a1−a38
a3=279.2−144a1−12a21
Ітерація #1
Задано початкове припущення вектора розв'язку як
[a1a2a3]=[125]
отримуємо
a1=106.8−5(2)−(5)25=3.6720
a2=177.2−64(3.6720)−(5)8=−7.8150
a3=279.2−144(3.6720)−12(−7.8510)1=−155.36
Абсолютна відносна приблизна похибка для кожногоxi тоді дорівнює
|∈a|1=|3.6720−13.6720|×100=72.76%
|∈a|2=|−7.8510−2−7.8510|×100=125.47%
|∈a|3=|−155.36−5−155.36|×100=103.22%
В кінці першої ітерації оцінка вектора розв'язку дорівнює
[a1a2a3]=[3.6720−7.8510−155.36]
і максимальна абсолютна відносна приблизна похибка дорівнює125.47.
Ітерація #2
Оцінка вектора розв'язку в кінці Ітерації #1 дорівнює
[a1a2a3]=[3.6720−7.8510−155.36]
Тепер отримуємо
a1=106.8−5(−7.8510)−(−155.36)25=12.056
a2=177.2−64(12.056)−(−155.36)8=−54.882
a3=279.2−144(12.056)−12(−54.882)1=−798.34
Абсолютна відносна приблизна похибка для кожногоxi тоді дорівнює
|∈a|1=|12.056−3.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.056−54.882−798.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|≥n∑j=1j≠i|aij| for all i
|aii|>n∑j=1j≠i|aij|for at least one i
Якщо система рівнянь має матрицю коефіцієнтів, яка не є домінуючою по діагоналі, вона може сходитися або не сходитися. На щастя, багато фізичних систем, що призводять до одночасних лінійних рівнянь, мають діагонально домінуючу матрицю коефіцієнтів, яка потім забезпечує збіжність для ітераційних методів, таких як метод Гаусса-Зайделя розв'язування одночасних лінійних рівнянь.
Знайти розв'язок наступної системи рівнянь за допомогою методу Гаусса-Зайделя.
12x1+3x2−5x3=1
x1+5x2+3x3=28
3x1+7x2+13x3=76
Використовувати
[x1x2x3]=[101]
як початкове вгадати і провести дві ітерації.
Рішення 1
матриця коефіцієнтів
[A]=[123−51533713]
є діагонально домінуючим як
|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=1−3x2+5x312
x2=28−x1−3x35
x3=76−3x1−7x213
Припускаючи початкове припущення
[x1x2x3]=[101]
Ітерація #1
x1=1−3(0)+5(1)12=0.50000
x2=28−(0.50000)−3(1)5=4.9000
x3=76−3(0.50000)−7(4.9000)13=3.0923
Абсолютна відносна приблизна похибка в кінці першої ітерації дорівнює
|∈a|1=|0.50000−10.50000|×100=100.00%
|∈a|2=|4.9000−04.9000|×100=100.00%
|∈a|3=|3.0923−13.0923|×100=67.662%
Максимальна абсолютна відносна приблизна похибка дорівнює100.00%
Ітерація #2
x1=1−3(4.9000)+5(3.0923)12=0.14679
x2=28−(0.14679)−3(3.0923)5=3.7153
x3=76−3(0.14679)−7(3.7153)13=3.8118
В кінці другої ітерації абсолютна відносна приблизна похибка дорівнює
|∈a|1=|0.14679−0.500000.14679|×100=240.61%
|∈a|2=|3.7153−4.90003.7153|×100=31.889%
|∈a|3=|3.8118−3.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]
Задано систему рівнянь
3x1+7x2+13x3=76
x1+5x2+3x3=28
12x1+3x2−5x3=1
знайти рішення за допомогою методу Гауса-Зайделя. Використовувати
[x1x2x3]=[101]
як початкова здогадка.
Рішення 1
Переписуючи рівняння, отримуємо
x1=76−7x2−13x33
x2=28−x1−3x35
x3=1−12x1−3x2−5
Припускаючи початкове припущення
[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]=[3713153123−5]
не є домінуючим по діагоналі, як
|a11|=|3|=3≤|a12|+|a13|=|7|+|13|=20
Значить, метод Гауса-Зайделя може сходитися, а може і не сходитися.
Однак це той самий набір рівнянь, що і попередній приклад, і який сходився. Єдина відмінність полягає в тому, що ми обмінялися першим і третім рівнянням один з одним, і це зробило матрицю коефіцієнтів не домінуючою по діагоналі.
Тому не виключено, що система рівнянь може бути зроблена діагонально домінуючою, якщо обмінюватися рівняннями один з одним. Однак це можливо не для всіх випадків. Наприклад, наступний набір рівнянь
x1+x2+x3=3
2x1+3x2+4x3=9
x1+7x2+x3=9
неможливо переписати, щоб зробити матрицю коефіцієнтів по діагоналі домінуючою.
Вікторина методу Гауса-Зайделя
Квадратна матриця[A]n×n є домінуючою по діагоналі, якщо
(А)|aii|≥n∑j=1i≠j|aij|,i=1,2,...,n
(B)|aii|≥n∑j=1i≠j|aij|,i=1,2,...,n і|aii|>∑nj=1i≠j|aij|, для будь-якогоi=1,2,...,n
(C)|aii|≥n∑j=1|aij|,i=1,2,...,n і|aii|>∑nj=1|aij|, для будь-якогоi=1,2,...,n
(D)|aii|≥n∑j=1|aij|,i=1,2,...,n
Використовуючи[x1,x2,x3]=[1,3,5] як початкову здогадку, значення[x1,x2,x3] після трьох ітерацій у методі Гаусса-Зайделя для
[127315127−11][x1x2x3]=[2−56]
є
(А)[−2.8333−1.4333−1.9727]
(Б)[1.4959−0.90464−0.84914]
(С)[0.90666−1.0115−1.0243]
(D)[1.2148−0.72060−0.82451]
Для того, щоб наступна система рівнянь,
2x1+7x2−11x3=6x1+2x2+x3=−57x1+5x2+2x3=17
сходиться за допомогою методу Гаусса-Зайделя, можна переписати вищевказані рівняння наступним чином:
(А)[27−11121752] [x1x2x3]=[6−517]
(Б)[75212127−11] [x1x2x3]=[17−56]
(С)[75212127−11] [x1x2x3]=[6−517]
(D) Рівняння не можна переписати у формі для забезпечення збіжності.
Для[127315127−11] [x1x2x3]=[227−2] та використання[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
Алгоритм розв'язання методу Гаусса-Зайделя[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 |
Термістори вимірюють температуру, мають нелінійний вихід і цінуються для обмеженого діапазону. Отже, коли виготовляється термістор, виробник постачає криву опору проти температури. Точне уявлення кривої, як правило, задається
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
Вправа методу Гауса-Зайделя
У системі рівнянь[A][X]=[C], якщо[A] діагонально домінантний, то метод Гауссейдала-Зайделя
- завжди сходиться
- може сходитися або не сходитися
- завжди розходиться
- Відповідь
-
A
У системі[A][X]=[C] рівнянь якщо не[A] є домінуючою по діагоналі, то метод Гаусса-Зайделя
- Завжди сходиться
- Може сходитися або не сходитися
- Завжди розходиться.
- Відповідь
-
Б
У системі рівнянь[A][X]=[C], якщо[A] вона не домінуюча по діагоналі, система рівнянь завжди може бути переписана, щоб зробити її домінуючою по діагоналі.
- Правда
- Помилковий
- Відповідь
-
Б
Розв'яжіть наступну систему рівнянь методом Гауса-Зайделя
12x1+7x2+3x3=2x1+5x2+x3=−52x1+7x2−11x3=6
Проведіть 3 ітерації, обчислите максимальну абсолютну відносну приблизну похибку в кінці кожної ітерації і виберіть[x1x2x3]=[135] як початкову здогадку.
- Відповідь
-
[x1x2x3]=[0.90666−1.0115−1.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%]
Розв'яжіть наступну систему рівнянь методом Гауса-Зайделя
12x1+7x2+3x3=2
x1+5x2+x3=−5
2x1+7x2−11x3=6
Проведіть 3 ітерації, обчислите максимальну абсолютну відносну приблизну похибку в кінці кожної ітерації та виберіть[x1x2x3]=[135] як початкову здогадку.
- Відповідь
-
[x1x2x3]=[0.90666−1.0115−1.0243]
Вирішіть наступну систему рівнянь методом Гауса-Зайделя
x1+5x2+x3=5
12x1+7x2+3x3=2
2x1+7x2−11x3=6
Проведіть 3 ітерації, обчислите максимальну абсолютну відносну наближену похибку в кінці кожної ітерації, і вибрати[x1x2x3]=[135] в якості початкової здогадки.
- Відповідь
-
[x1x2x3]=[−1163.71947.61027.2]
[|∈a|1|∈a|2|∈a|3]=[89.156%89.139%89.183%]