ВУЗ:
Составители:
⎪
⎪
⎪
⎪
⎩
⎪
⎪
⎪
⎪
⎨
⎧
=++++−−−−
=++++−−−−
=+++++++
=+++++++
=+++++++
n
nnnnnnnnnn
nnnn
nnnnnnnnnnn
nnnn
nnnn
JUgUgUgUbUbUb
JUgUgUgUbUbUb
JUbUbUbUgUgUg
JUbUbUbUgUgUg
JUbUbUbUgUgUg
"'...""'...''
....
"'...""'...''
'"...""'...''
.....
'"...""'...
''
'"...""'...''
"22112211
1"12121111212111
212112211
222221212222121
112121111212111
Формализовать задачу, т.е. коэффициенты и неизвестные
записать в общем виде
⎪
⎪
⎩
⎪
⎪
⎨
⎧
=+++
=+++
=+++
nnnnnn
nn
nn
cUaUaUa
cUaUaUa
cUaUaUa
...
....
...
...
2211
22222121
11212111
(20)
Глобально алгоритм разделяется на 2 этапа: прямой ход и
обратный. Задача прямого хода – свести систему
уравнений к треугольной типа
⎪
⎪
⎩
⎪
⎪
⎨
⎧
=
=++
=+++
nnnn
nn
nn
bUb
bUbUb
bUbUbU
....
...
...
22222
112121
(21)
Обратный ход – нахождение решения подстановкой
решений снизу
Общее количество действий
3
n≈ – много!
Последовател
ьность
действий при
выполнении
алгоритма
единственног
о деления
Прямой ход
Действие 1.
преобразовать матрицу n×n таким образом,
чтобы диагональные элементы были доминирующими,
т.е. по абсолютной величине наибольшими в строке. При
этом столбик переменных перемещается полностью
Пример: -2u
1
+3u
2
-6u
3
+u
4
=-6 -6u
3
-2u
1
+3u
2
+u
4
=-6
5u
1
- u
2
+3u
3
=5 3u
3
+5u
1
–u
2
=5
u
1
+7u
2
-2u
4
=1 u
1
+7u
2
-2u
4
=1
-u
1
-2u
2
+3u
3
-4u
4
=-8 3u
3
-u
1
-2u
2
-4u
4
=-8
Действие 2. Прямой ход.
Шаг 1. Исключаем из всех уравнений, начиная со
второго, первый член, для чего
а) разделить все коэффициенты и свободный член на
коэффициент при первом уравнении (на –6);
б) получившееся в результате шага 1а уравнение
умножить на коэффициент перед первым членом
i-того
уравнения с обратным знаком и сложить эти два
уравнения. В результате действий первый член в
i-том
уравнении обратится в 0. (2: 1а
×-3, сложить 1а и рез-т)
Шаг 2. Исключить из всех уравнений, начиная с 3-го 2-й
член, для чего
а) уравнение 1а переписать без изменений
б) для укороченной на 1-й член системы проделать
операции шага 1
Шаг 3 и далее. Повторять операции шага 2, пока не
придем к треугольной системе уравнений.
Действие 3
. Обратный ход.
Схема
Жордана
=алгоритм без обратного хода:
действие1, шаг 1 те же. Начиная с шага 2 проводятся
дополнительные операции с уравнениями от1 до
i-1, в
которых проводятся преобразования таким образом,
чтобы коэффициент столбца, соответствующий
преобразуемому уравнению обратился в 0. Операций еще
больше, поэтому редко применяется.
Достоинства
и недостатки
метода
Гаусса
(Жордана)
Достоинства
: простота, наглядность
Недостатки: возможна большая погрешность, вызванная
1) округлением результатов расчетов (диагональный
элемент близок к нулю – избегаем, но на следующих
шагах возможно)
2) неточностью задания исходных данных. Всегда
присутствует погрешность в задании параметров
системы. Как правило, погрешность полученных
решений бывает соизмерима с погрешностью задания
параметров системы, т.е. коэффициентов матрицы.
Однако
могут быть случаи, когда погрешность исх.
данных дает очень большую погрешность решения.
Причина –
плохая обусловленность матрицы
коэффициентов, признаком кторой является малость
определителя матрицы относительно оценки Адамара
⎧ g11U '1 + g12U ' 2 +... + g1nU ' n +b11U "1 +b12U "2 +... + b1nU "n = J 1 ' коэффициент при первом уравнении (на –6); ⎪ б) получившееся в результате шага 1а уравнение ⎪ g 21U '1 + g 22U ' 2 +... + g 2 nU ' n +b21U "1 +b22U "2 +... + b2nU "n = J 2 ' умножить на коэффициент перед первым членом i-того ⎪..... ⎪ уравнения с обратным знаком и сложить эти два ⎨ g n1U '1 + g n 2U ' 2 +... + g nnU ' n +bn1U "1 +bn12U "2 +... + bnnU "n = J n ' уравнения. В результате действий первый член в i-том ⎪− b U ' −b U ' −... − b U ' + g U " + g U " +... + g U ' = J " ⎪ 11 1 12 2 1n n 11 1 12 2 1n "n 1 уравнении обратится в 0. (2: 1а×-3, сложить 1а и рез-т) ⎪.... Шаг 2. Исключить из всех уравнений, начиная с 3-го 2-й ⎪ член, для чего ⎩− bn1U '1 −bn 2U ' 2 −... − bnnU ' n + g n1U "1 + g n 2U "2 +... + g nnU '"n = J "n а) уравнение 1а переписать без изменений Формализовать задачу, т.е. коэффициенты и неизвестные б) для укороченной на 1-й член системы проделать записать в общем виде ⎧a11U 1 + a12U 2 + ... + a1nU n = c1 операции шага 1 ⎪ ⎪a 21U 1 + a 22U 2 + ... + a 2nU n = c 2 Шаг 3 и далее. Повторять операции шага 2, пока не ⎨ (20) ⎪.... придем к треугольной системе уравнений. ⎪⎩a n1U 1 + a n 2U 2 + ... + a nnU n = c n Действие 3. Обратный ход. Схема =алгоритм без обратного хода: Глобально алгоритм разделяется на 2 этапа: прямой ход и Жордана действие1, шаг 1 те же. Начиная с шага 2 проводятся обратный. Задача прямого хода – свести систему дополнительные операции с уравнениями от1 до i-1, в уравнений к треугольной типа которых проводятся преобразования таким образом, ⎧U1 + b12U 2 + ... + b1nU n = b1 ⎪ чтобы коэффициент столбца, соответствующий ⎪ b22U 2 + ... + b2 nU n = b2 преобразуемому уравнению обратился в 0. Операций еще ⎨ (21) ⎪.... больше, поэтому редко применяется. ⎪⎩ bnnU n = bn Достоинства Достоинства: простота, наглядность Обратный ход – нахождение решения подстановкой и недостатки Недостатки: возможна большая погрешность, вызванная решений снизу метода 1) округлением результатов расчетов (диагональный Гаусса элемент близок к нулю – избегаем, но на следующих Общее количество действий ≈ n 3 – много! (Жордана) шагах возможно) Последовател Прямой ход 2) неточностью задания исходных данных. Всегда ьность Действие 1. преобразовать матрицу n×n таким образом, присутствует погрешность в задании параметров действий при чтобы диагональные элементы были доминирующими, системы. Как правило, погрешность полученных выполнении т.е. по абсолютной величине наибольшими в строке. При решений бывает соизмерима с погрешностью задания алгоритма этом столбик переменных перемещается полностью параметров системы, т.е. коэффициентов матрицы. единственног Пример: -2u1+3u2-6u3+u4=-6 -6u3-2u1+3u2+u4=-6 Однако могут быть случаи, когда погрешность исх. о деления 5u1- u2+3u3 =5 3u3+5u1 –u2 =5 данных дает очень большую погрешность решения. u1+7u2 -2u4=1 u1+7u2-2u4=1 Причина – плохая обусловленность матрицы -u1-2u2+3u3-4u4=-8 3u3 -u1 -2u2 -4u4=-8 коэффициентов, признаком кторой является малость Действие 2. Прямой ход. определителя матрицы относительно оценки Адамара Шаг 1. Исключаем из всех уравнений, начиная со второго, первый член, для чего а) разделить все коэффициенты и свободный член на
Страницы
- « первая
- ‹ предыдущая
- …
- 38
- 39
- 40
- 41
- 42
- …
- следующая ›
- последняя »