Опубликован: 15.03.2007 | Уровень: специалист | Доступ: платный | ВУЗ: Донецкий национальный технический университет
Лекция 9:

Компьютерное моделирование и решение линейных и нелинейных многомерных систем

< Лекция 8 || Лекция 9: 123 || Лекция 10 >
Аннотация: Лекция рассматривает метод и алгоритм решения систем линейных уравнений методом Гаусса

При моделировании экономических задач, таких как задачи управления и планирования производства, определения оптимального размещения оборудования, оптимального плана производства, оптимального плана перевозок грузов (транспортная задача), распределения кадров и др., может быть положена гипотеза линейного представления реального мира.

Математические модели таких задач представляются линейными уравнениями. Если задача многомерна, то ее математическая модель представляется системой линейных уравнений.

Линейные математические модели также используются в нелинейных системах при условии, если эта нелинейная система условно линеаризирована.

В общем виде система линейных уравнений имеет вид:

\left\{ \begin{array}{l}
a_{11}x_1 + a_{12}x_2 + a_{13}x_3 + a_{14}x_4 + \ldots + a_{1n}x_n = b_1,\\ 
a_{21}x_1 + a_{22}x_2 + a_{23}x_3 + a_{24}x_4 + \ldots + a_{2n}x_n = b_2,\\ 
a_{31}x_1 + a_{32}x_2 + a_{33}x_3 + a_{34}x_4 + \ldots + a_{3n}x_n = b_3,\\ 
\ldots\\ 
a_{n1}x_1 + a_{n2}x_2 + a_{n3}x_3 + a_{n4}x_4 + \ldots + a_{nn}x_n = b_n, 
\end{array} \right.

где

aij - коэффициенты при неизвестных системы,

bi - свободные члены,

xj - неизвестные системы,

i=\overline{1,n} - номер строки,

j=\overline{1,n} - номер столбца,

n - порядок системы.

В матричной форме система линейных уравнений имеет вид:

A \cdot \overline X = \overline B,

где

{A} = \left [\begin{array}{ccc} a_{11}  a_{12}  \ldots  a_{1n}\\ 
a_{21}  a_{22}  \ldots  a_{2n}\\ 
\ldots  \ldots  \ldots  \ldots\\ 
a_{n1}  a_{n2}  \ldots  a_{nn}\\ 
\end{array} \right] \text{ –матрица коэффициентов системы порядка } {(n \times n)},
\mathbf{\overline X} = \left[ \begin{array}{ccc} x_1\\ x_2\\ \ldots\\ x_n\\ \end{array} \right] \text{ – вектор неизвестных системы размерностью n},
\mathbf{\overline B} = \left[ \begin{array}{ccc} b_1\\ b_2\\ \ldots\\ b_n\\ \end{array} \right] \text{ – вектор свободных членов размерностью n}.

Численные методы решения систем линейных уравнений (СЛУ) можно разделить на две группы:

  1. точные или прямые методы,
  2. приближенные методы.

Приближенные методы реализуют на ЭВМ нахождение корней с заданной точностью и являются итерационными методами.

Точные методы позволяют получить решение системы за конечное число итераций. К точным методам относятся:

  • правило Крамера,
  • метод Гаусса,
  • метод прогонки.

Решение систем линейных уравнений методом Гаусса

Метод Гаусса является точным методом. Он позволяет получить решение системы за конечное число арифметических действий. В основе метода лежит идея последовательного исключения неизвестных. Метод состоит из двух этапов. На первом этапе (прямой ход) система при помощи последовательного исключения неизвестных приводится к треугольному виду. На втором этапе (обратный ход) из системы треугольного вида последовательно, в обратном порядке, начиная c n-го уравнения, находятся неизвестные системы.

В качестве примера возьмем систему 4 порядка.

\left\{ \begin{array}{l} 
a_{11}^{(0)}x_1 + a_{12}^{(0)}x_2 + a_{13}^{(0)}x_3 + a_{14}^{(0)}x_4 = b_1^{(0)},\\ 
a_{21}^{(0)}x_1 + a_{22}^{(0)}x_2 + a_{23}^{(0)}x_3 + a_{24}^{(0)}x_4 = b_2^{(0)},\\ 
a_{31}^{(0)}x_1 + a_{32}^{(0)}x_2 + a_{33}^{(0)}x_3 + a_{34}^{(0)}x_4 = b_3^{(0)},\\ 
a_{41}^{(0)}x_1 + a_{42}^{(0)}x_2 + a_{43}^{(0)}x_3 + a_{44}^{(0)}x_4 = b_4^{(0)}, 
\end{array} \right. ( 9.1)

Прямой ход. На первом шаге прямого хода (к=1) находим x1 из первого уравнения системы (9.1).

a_{11}^{(0)} - ведущий элемент первой строки.

Если a_{11}^{(0)} \neq 0, то

x_1= - \frac {a_{12}^{(0)}}{a_{11}^{(0)}} x_2 - \frac {a_{13}^{(0)}}{a_{11}^{(0)}} x_3 - \frac {a_{14}^{(0)}}{a_{11}^{(0)}} x_4 + \frac {b_1^{(0)}}{a_{11}^{(0)}}, ( 9.2)

Обозначим:

\frac {a_{12}^{(0)}}{a_{11}^{(0)}} = a_{12}^{(1)}; \frac {a_{13}^{(0)}}{a_{11}^{(0)}} = a_{13}^{(1)}; \frac {a_{14}^{(0)}}{a_{11}^{(0)}} = a_{14}^{(1)}; \frac {b_1^{(0)}}{a_{11}^{(0)}} = b_1^{(1)}. ( 9.3)

Подставляя (9.3) в (9.2), получим

x_1 = -a_{12}^{(1)}x_2 - a_{13}^{(1)}x_3 - a_{14}^{(1)}x_4 + b_1^{(1)}, ( 9.4)

где

a_{1j}^{(1)} = \frac{a_{1j}^{(0)}}{a_{11}^{(0)}}, j=2,3,4 = \overline{(k+1),n}\\ 
b_1^{(1)} = \frac{b_1^{(0)}}{a_{11}^{(0)}}.

Подставляем (9.4) во 2, 3 и 4 уравнение системы (9.1), получим:

\left\{ \begin{array}{l}
(a_{22}^{(0)} - a_{21}^{(0)} \cdot a_{12}^{(1)}) x_2 + (a_{23}^{(0)} - a_{21}^{(0)} \cdot a_{13}^{(1)}) x_3 + (a_{24}^{(0)} - a_{21}^{(0)} \cdot a_{14}^{(1)}) x_4 = b_2^{(0)} - a_{21}^{(0)}\cdot b_1^{(1)},\\ 
(a_{32}^{(0)} - a_{31}^{(0)} \cdot a_{12}^{(1)}) x_2 + (a_{33}^{(0)} - a_{31}^{(0)} \cdot a_{13}^{(1)}) x_3 + (a_{34}^{(0)} - a_{31}^{(0)} \cdot a_{14}^{(1)}) x_4 = b_3^{(0)} - a_{31}^{(0)}\cdot b_1^{(1)},\\ 
(a_{42}^{(0)} - a_{41}^{(0)} \cdot a_{12}^{(1)}) x_2 + (a_{43}^{(0)} - a_{41}^{(0)} \cdot a_{13}^{(1)}) x_3 + (a_{44}^{(0)} - a_{41}^{(0)} \cdot a_{14}^{(1)}) x_4 = b_4^{(0)} - a_{41}^{(0)}\cdot b_1^{(1)}, 
\end{array} \right.

Обозначив коэффициенты при неизвестных полученной системы через a_{ij}^{(1)}, а свободные члены через b_i^{(1)} перепишем полученную систему:

\left\{ \begin{array}{l}
a_{22}^{(1)} x_2 + a_{23}^{(1)} x_3 + a_{24}^{(1)} x_4 = b_2^{(1)},\\ 
a_{32}^{(1)} x_2 + a_{33}^{(1)} x_3 + a_{34}^{(1)} x_4 = b_3^{(1)},\\ 
a_{42}^{(1)} x_2 + a_{43}^{(1)} x_3 + a_{44}^{(1)} x_4 = b_4^{(1)}. 
\end{array} \right. ( 9.5)

где

a_{ij}^{(1)} = a_{ij}^{(0)} - a_{il}^{(0)} \cdot a_{lj}^{(1)},\\ 
b_i^{(1)} =b_i^{(0)} -a_{il}^{(0)} b_l^{(1)},\\ i=\overline{(k+1),n}; j=\overline{(k+1),n}.

Таким образом, в результате выполнения первого шага прямого хода исходная система (9.1) n-го порядка преобразована к совокупности уравнения (9.4) и системы линейных уравнений (9.5), порядок которой равен n-1.

На втором шаге прямого хода (к=2) из первого уравнения системы (9.5) находим x2.

a_{22}^{(1)} -ведущий элемент первой строки системы (9.5).

Если a_{22}^{(1)} \neq 0, то из первого уравнения системы (9.5) имеем:

x_2= - a_{23}^{(2)}x_3 - a_{24}^{(2)}x_4 + b_2^{(2)}, ( 9.6)

где

a_{2j}^{(2)} = a_{2j}^{(1)}/ a_{22}^{(1)},\\ 
b_2^{(2)} = b_2^{(1)}/ a_{22}^{(1)}.\\ j=3,4=\overline{(k+1),n}

Подставив выражение (9.6) во второе и третье уравнения системы (9.5), получим новую систему линейных уравнений, порядок которой равен n-2.

\left\{ \begin{array}{l}
a_{33}^{(2)}x_3 + a_{34}^{(2)}x_4 = b_3^{(2)},\\ 
a_{43}^{(2)}x_3 + a_{44}^{(2)}x_4 = b_4^{(2)}. 
\end{array} \right. ( 9.7)

где

a_{ij}^{(2)} = a_{ij}^{(1)} - a_{i2}^{(1)} \cdot a_{2j}^{(2)}, \\b_i^{(2)} = b_i^{(1)} - a_{i2}^{(1)}b_2^{(2)},\\ 
i=\overline{(k+1),n}; j=\overline{(k+1),n}.

Таким образом, в результате выполнения второго шага прямого хода исходная система (9.1) преобразована к совокупности уравнений (9.4), (9.6) и системы линейных уравнений (9.7),порядок которой равен n-2.

< Лекция 8 || Лекция 9: 123 || Лекция 10 >
Равиль Султанов
Равиль Султанов

В уравнениях движения кривошипно-шатунного механизма вместо обозначения радиуса кривошипа "r" ошибочно записан символ "γ" (гамма).

P.S. Может быть это слишком очевидно, но не упомянуто, что угол поворота кривошипа φ считается малым.

Александр Никитин
Александр Никитин

Добрый день.

В расчете параметра Т4 xi суммируется с величиной h/2 ?

Yusupov Ozod
Yusupov Ozod
Узбекистан, Samar
Владимир Ленчицкий
Владимир Ленчицкий
Россия, Губкинский