Система линейных уравнений

Система линейных уравнений Реферат Студентки группы 1113а Лариной Натальи Владиславовны научный руководитель Винокурова Татьяна Васильевна, к. м. н доцент кафедры математического анализа Владивосток 2007 Содержание Глава 1. Критерий совместимости….3 Глава 2. Метод Гаусса….5 Глава 3. Формулы Крамера… …11

Глава 4. Матричный метод… 14 Список литературы…15 Глава 1 Критерий совместимости Система линейных уравнений имеет вид: a11x1 + a12x2 + + a1nxn = b1 a21x1 + a22x2 + + a2nxn = b2 (5.1) … am1x2 + am2x2 + + amnxn = bm Здесь аij и bi (i = ; j = ) - заданные, а xj - неизвестные действительные числа. Используя понятие произведения матриц, можно переписать систему (5.1) в виде:

AX = B, (5.1) где A = (аij) - матрица, состоящая из коэффициентов при неизвестных системы (5.1), которая называется матрицей системы, X = (x1, x2 xn)T, B = (b1, b2 bm)T - векторы-столбцы, составленные соответственно из неизвестных xj и из свободных членов bi Упорядоченная совокупность n вещественных чисел (c1, c2 cn) называется решением системы (5.1), если в результате подстановки этих чисел вместо соответствующих переменных x1, x2 xn каждое уравнение системы обратится в арифметическое тождество; другими словами, если существует

вектор C= (c1, c2 cn)T такой, что AC ≡ B. Система (5.1) называется совместной, или разрешимой, если она имеет по крайней мере одно решение. Система называется несовместной, или неразрешимой, если она не имеет решений. Матрица Ã = образованная путем приписывания справа к матрице A столбца свободных членов, называется расширенной матрицей системы.

Вопрос о совместности системы (5.1) решается следующей теоремой. Теорема Кронекера- Капелли. Система линейных уравнений совместна тогда и только тогда, когда ранги матриц A и Ã совпадают, т.е. r(A) = r(Ã) = r Для множества М решений системы (5.1) имеются три возможности: 1. M = Ø (в этом случае система несовместна);

2. M состоит из одного элемента, т.е. система имеет единственное решение (в этом случае система называется определенной); 3. M состоит более чем из одного элемента (тогда система называется неопределенной). В третьем случае система (5.1) имеет бесчисленное множество решений. Система имеет единственное решение только в том случае, когда r(A) = n. При этом число уравнений - не меньше числа неизвестных (m ≥ n); если m > n, то m-n уравнений

являются следствиями остальных. Если 0 < r < n, то система является неопределенной. Для решения произвольной системы линейных уравнений нужно уметь решать системы, в которых число уравнений равно числу неизвестных так называемые системы крамеровского типа: a11x1 + a12x2 + + a1nxn = b1 a21x1 + a22x2 + + a2nxn = b2 (5.3) … an1x2 + an2x2 + + annxn = bn Системы (5.3) решаются одним из следующих способов:

1. методом Гаусса, или методом исключения неизвестных; 2. по формулам Крамера; 3. матричным методом. Глава 2 Метод Гаусса Исторически первым, наиболее распространенным методом решения систем линейных уравнений является метод Гаусса, или метод последовательного исключения неизвестных. Сущность этого метода состоит в том, что посредством последовательных исключений неизвестных данная

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

Возвратимся матрицам. Рассмотрим теперь прямоугольную матрицу, имеющую m строк и n столбцов. Ее называют матрица размером m на n. А(mхn). Выделим в этой матрице произвольные к строк и к столбцов. Они образуют квадратную матрицу B(kхk) Например: Минором К-ого порядка матрицы А называется определитель квадратной матрицы, получающейся из данной матрицы выделением произвольных к строк и к столбцов.  B =detB- является минором третьего порядка.

Минором второго порядка является, например определитель Сами элементы матрицы можно рассматривать как миноры первого порядка. Какие-то из миноров равны нулю, какие-то нет. Рангом матрицы называется наибольший из порядков ее миноров, отличных от нуля. Если ранг А обозначаемый r (A) равен r, то это означает, что в А имеется хотя бы один отличный от нуля минор порядка r, но всякий минор, порядка больше чем r, равен

нулю. Например, найдем ранг матрицы 1. Проверяем минор 4 порядка  A =0, т.к. матрица содержит нулевой столбец. 2. Проверяем миноры 3 порядка Итак, процесс вычислений миноров прекращаем, поскольку миноров 4 порядка, не равных нулю нет, а минор 3-го порядка найден. Значит r(A)=3. Мы уже рассматривали методы решения систем линейных алгебраических уравнений, но только в тех случаях,

когда матрица системы – квадратная, то есть число уравнений равно числу неизвестных. Рассмотрим теперь самый простой и употребительный способ решения систем линейных уравнений – метод Гаусса. Рассмотрим его на простейшем примере, решая систему трех уравнений с тремя неизвестными. Мы хотим исключить х1 из всех уравнений, кроме первого. Для этого мы должны вычесть из второго уравнения первое, умноженное на 4, а к третьему прибавить первое,

умноженное на 5. На втором шаге исключения мы не трогаем первое уравнение. Другие два уравнения содержат два неизвестных х2 и х3 и к ним можно применить ту же процедуру исключения. Для этого к третьему уравнению прибавляем второе, умноженное на 3. Далее наши действия очевидны. Из третьего уравнения х3=-1, подставляя это значение во второе уравнение, получаем х2=-3 и наконец, из первого уравнения получаем х1=2.

Этот простой процесс называется простой подстановкой. Таким образом, процесс решения системы линейных алгебраических уравнений по методу Гаусса состоит из двух этапов. Первый этап (прямой ход метода) – система приводится к треугольному виду. Второй этап (обратный ход) – неизвестные определяются последовательно, начиная с последнего неизвестного и кончая первым. Аналогично, эту идею последовательного исключения можно применить и в случае матрицы

А(mxm) размера больше 3х3. Без ограничения общности можно считать, что в нашей системе ведущий элемент a11 0 первого шага (иначе просто переставим уравнение). На первом шаге мы просто исключим х1 из всех уравнений, начиная со второго, для чего из второго уравнения почленно вычтем первое, умноженное на а21/а11, из третьего почленно вычтем первое, помноженное на а31/а11 и т.д Тогда система заменится эквивалентной системой:

Коэффициенты при неизвестных и свободные члены в последних m-1 уравнениях системы, определяются формулами: Таким образом, на первом шаге уничтожаются все коэффициенты, лежащие под первым ведущим элементом a11 0 На втором шаге уничтожаются элементы, лежащие под вторым ведущим элементом а22(1) (если a22(1) 0) Продолжая этот процесс и дальше, мы, наконец, на (m-1) шаге приведем исходную систему к треугольной системе. Матрица этой системы имеет вид: На этом прямой ход метода

Гаусса заканчивается. Второй этап – обратный ход, заключается в решении треугольной системы. Из последнего уравнения находим xm. По найденному xm из (m-1) уравнения находим xm-1. Затем по xm-1 и xm из (m-2) уравнения находим xm-2. Процесс продолжаем, пока не найдем x1 из первого уравнения. Если у нас число уравнений меньше числа неизвестных, то мы придем не к треугольной системе, а к ступенчатой.

Так как прямой ход метода Гаусса прервется, когда уравнения закончатся, а неизвестные еще останутся. В таком случае в каждом уравнении системы перенесем все члены с неизвестными xk+1,… xm в правую часть. Придавая неизвестным xk+1,… xm (называемым свободными) произвольные значения, получим треугольную систему из которой последовательно найдем все остальные неизвестные (называемые базисными). Так как произвольные значения можно придавать любыми способами, система будет иметь бесчисленное множество

значений. В решении следующего примера не будем выписывать каждую систему, а ограничимся лишь преобразованиями над матрицами: и Такая модификация метода называется методом Жордана-Гауcса. Решение систем линейных алгебраических уравнений методом Жордана – Гаусса 1 шаг: а) первую строку не меняем б) из второй вычитаем первую, умноженную на 2 в) третью не меняем, т.к. там неизвестное х1 и так отсутствует.

2 шаг: а) вторую строку делим на - 4 б) из третьей строки вычитаем новую вторую (поделенную на -4). 3 шаг: делим третью строку на (-7/4). Последней матрице соответствует система: или х3 = -2 + 10/7х4 + 3/7х5 Глава 3 Формулы Крамера Метод Крамера состоит в том, что мы последовательно находим главный определитель системы (5.3), т.е. определитель матрицы А  = det (ai j) и n вспомогательных определителей  i (i= ), которые получаются

из определителя  заменой i-го столбца столбцом свободных членов. Формулы Крамера имеют вид:   x i =  i ( i = ), (5.4) Из (5.4) следует правило Крамера, которое дает исчерпывающий ответ на вопрос о совместности системы (5.3): если главный определитель системы отличен от нуля, то система имеет единственное решение, определяемое по формулам: x i =  i / . Если главный определитель системы  и все

вспомогательные определители  i = 0 (i= ), то система имеет бесчисленное множество решений. Если главный определитель системы  = 0, а хотя бы один вспомогательный определитель отличен от нуля, то система несовместна. Пример. Решить методом Крамера систему уравнений: x1 + x2 + x3 + x4 = 5, x1 + 2x2 - x3 + 4x4 = -2, 2x1 - 3x2 - x3 - 5x4 = -2, 3x1 + x2 +2x3 + 11 x4 = 0. Решение.

Главный определитель этой системы: значит, система имеет единственное решение. Вычислим вспомогательные определители  i ( i = ), получающиеся из определителя  путем замены в нем столбца, состоящего из коэффициентов при xi, столбцом из свободных членов: Отсюда x1 =  1/ = 1, x2 =  2/ = 2, x3 =  3/ = 3, x4 =  4/ = -1, решение системы - вектор

С=(1, 2, 3, -1)T. Глава 4 Матричный метод Если матрица А системы линейных уравнений невырожденная, т.е. det A  0, то матрица А имеет обратную, и решение системы (5.3) совпадает с вектором C = A1B. Иначе говоря, данная система имеет единственное решение. Отыскание решения системы по формуле X=C, C=A1B называют матричным способом решения системы,

или решением по методу обратной матрицы. Пример. Решить матричным способом систему уравнений: x1 - x2 + x3 = 6, 2x1 + x2 + x3 = 3 x1 + x2 +2x3 = 5 Решение. Обозначим: Тогда данная система уравнений запишется матричным уравнением AX=B. Поскольку, то матрица A невырождена и поэтому имеет обратную: Для получения решения X мы должны умножить вектор-столбец

B слева на матрицу A: X = A1B. В данном случае и, следовательно, Выполняя действия над матрицами, получим: x1 = 1/5(16+33-2& #61655;5) = 1/5 (6+9-10) = 1 x2 = 1/5 (-36 +13 - 15) = 1/5 (- 18 + 3 + 5) = -2 x3 = 1/5 (16 - 23 + 35) = 1/5 (6 -6 + 15) = 3

Итак, С = (1, -2, 3)T. Список литературы 1. Г.И. Кручкович. “Сборник задач по курсу высшей математике”, М. “Высшая школа”, 1973 год. 2. В.СШипачев. “Высшая математика”, М. “Высшая школа”, 1985 год. 3. Б.М. Владимирский, А.Б. Горстко, Я.М. Ерусалимский. «Математика. Общий курс», Ст Петербург, «Лань», 2002 год.