МИНОБРНАУКИ РОССИИ Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Мордовский государственный университет им. Н. П. Огарёва» В. Д. Бочкарева Математика для студентов географических направлений. Метод Гаусса решения систем линейных уравнений Учебно-методическое пособие Саранск 2012 СИСТЕМЫ ЛИНЕЙНЫХ УРАВНЕНИЙ. МЕТОД ГАУССА Сведения из теории Многие задачи естествознания своими моделями имеют системы линейных уравнений с несколькими неизвестными. Линейным уравнением с n неизвестными x1 , x 2 , , xn называется уравнение вида: a1 x1 a2 x2 an xn b , где a1 , a2 , …, an , b R (мы рассматриваем только вещественные уравнения, т.е. уравнения в которых коэффициенты, свободный член и значения неизвестных являются вещественными числами). Решением данного уравнения называется упорядоченный набор n вещественных чисел c1 , c 2 , …, c n , удовлетворяющих этому уравнению. Другими словами, упорядоченный набор =( c1 , c 2 , …, c n ) вещественных чисел называется решением уравнения a1 x1 a2 x2 an xn b , если числовое равенство c1 x1 c2 x2 cn xn b истинно. Системой m линейных уравнений с n неизвестными x1 , x 2 , , xn называется система вида: a11 x1 a12 x2 a1n xn b1 , a x a x a x b , 22 2 2n n 2 (1) 21 1 где aik R , bi R , , am1 x1 am 2 x2 amn xn bm , значения неизвестных – вещественные числа. Заметим, нумерация коэффициентов двойная. Например, a 25 (читается “а два пять”). Первый индекс указывает номер уравнения, а второй индекс – номер неизвестного, при котором стоит коэффициент. Например, a 25 означает, что это коэффициент из второго уравнения при пятом неизвестном. Свободные члены имеют только один индекс – номер уравнения, в котором этот свободный член находится. Например, b7 означает, что это свободный член седьмого уравнения. Число уравнений m может равняться числу неизвестных n . В этом случае система называется квадратной n – го порядка. Число m может быть меньше n , число m может быть больше n В этом случае система называется прямоугольной. По поведению свободных членов системы подразделяют на два типа: однородные и неоднородные. Система линейных уравнений называется однородной, если все ее свободные члены равны нулю. Система линейных уравнений называется неоднородной, хотя бы один из ее свободных членов отличен от нуля. Решением системы (1) называется такой упорядоченный набор n вещественных чисел c1 , c 2 , …, c n , который удовлетворяет каждому уравнению системы. Другими словами, упорядоченный набор =( c1 , c 2 , …, c n ) называется решением системы (1), если система числовых равенств a11c1 a12 c2 a1n cn b1 , a c a c a c b , 22 2 2n n 2 (2) 21 1 – истина. , am1c1 am 2 c2 amn cn bm , Решить систему – это значит найти множество ее решений. Если множество решений системы не пустое, то система называется совместной. Если множество решений системы пустое, то система называется несовместной. Если длина множества решений системы равна 1 (т.е. система имеет только одно решение), то система называется совместной и определенной. Если длина множества решений системы больше единицы (т.е. система имеет хотя бы два решения), то система называется совместной, но неопределенной. Решить линейную систему (найти множество ее решений) можно различными способами. В данной теме мы рассмотрим метод, который называется методом Гаусса или методом последовательного исключения неизвестных. Метод Гаусса базируется на применении следующих свойств решений линейной системы уравнений. Напомним, две системы линейных уравнений называются равносильными (эквивалентными), если множества их решений совпадают. Теорема 1. Если в системе поменять местами два уравнения, то система переходит в эквивалентную ей систему. Теорема 2. Если в системе какое-нибудь уравнение умножить (разделить) на число, отличное от нуля, то получится система, эквивалентная заданной. Теорема 3. Если в системе отбросить нулевое уравнение, то система перейдет в эквивалентную ей систему. Теорема 4. Если в системе есть два одинаковых уравнения, и одно из них отбросить, то система перейдет в эквивалентную систему. Теорема 5. Если в системе есть два пропорциональных уравнения, и одно из них отбросить, то система перейдет в эквивалентную систему. Теорема 6. Если в системе к какому-нибудь уравнению прибавить другое уравнение этой же системы, умноженное на любое число, то получится система, эквивалентная заданной. Все доказательства указанных теорем проводятся на основании определений решения линейного уравнения с n неизвестными и решения системы линейных уравнений. Заметим, что запись системы зависит от числа неизвестных, значений коэффициентов и свободных членов, но не зависит от названия неизвестных, а только от их номеров (неизвестные менять местами нельзя). Поэтому систему удобнее записывать в виде таблицы (матрицы), состоящей из коэффициентов и свободных членов системы, где вместо знаков «=» будет стоять вертикальная черта, отделяющая свободные члены от коэффициентов. Именно, вместо стандартной записи системы (1) мы будем иметь ее запись в виде следующей таблицы: a 11 a 12 a 1n b1 a 21 a 22 a 2n b 2 . a m1 a m2 a mn b n 2x1 x 2 3x 3 1, Например, система x1 - x 2 5, будет иметь следующую x 2x x 0. 2 3 1 1 3 1 2 таблицу: 1 1 0 5 . 1 2 1 0 Заметим, если в каком-то уравнении системы пропущена запись неизвестного, то это значит, что у этого неизвестного коэффициент равен нулю. Если известно количество неизвестных системы и их обозначение, то по таблице системы можно восстановить стандартную запись системы. Например, известно, что система линейных уравнений имеет неизвестные а, b, с и имеет таблицу: 1 1 1 2 0 1 0 2 . 0 0 1 1 1a - 1b 2c 1, В стандартном виде система имеет вид: 0a 1b 0c 2, или 0a 0b 1c -1. a - b 2c 1, 2, b c -1. Итак, теперь систему мы будем называть матрицей (таблицей), уравнение – строкой этой матрицы. Например, вместо «второе уравнение системы» мы будем говорить «вторая строка матрицы». Тогда перечисленные выше теоремы 1–6 будут звучать на языке «матрица, строки» и могут быть сформулированы следующим образом: Матрицу системы мы имеем право преобразовать: 1) менять в ней местами строки; 2) умножать делить все элементы строки на число, отличное от нуля; 3) отбрасывать нулевую строку; 4) из двух одинаковых строк одну отбрасывать; 5) из двух пропорциональных строк одну отбрасывать; 6) к строке матрицы прибавлять другую строку этой же матрицы, умноженное на любое число. При этом мы будем получать матрицу системы, эквивалентной исходной. Заметим, чтобы умножить строку на число, надо элемент (число) этой строки умножить на это число; чтобы к строке прибавить другую строку, надо к каждому элементу данной строки прибавить соответствующий элемент другой строки. Метод Гаусса предполагает выполнение специальных стандартных шагов преобразования исходной системы для получения эквивалентной системы наиболее простой конструкции. Рассмотрим применение метода Гаусса на конкретном примере. Пусть дана система линейных уравнений: x1 x 2 x 3 1, 2x 1 x 2 x 3 2, x x 2 x 1. 2 3 1 1 1 1 1 Запишем систему в виде матрицы: 2 1 1 2 . 1 1 2 1 1 шаг. Среди строк таблицы выбираем одну (любую), в которой содержится хотя бы один отличный от нуля коэффициент. В нашем примере в качестве такой строки подойдет и первая, и вторая, и третья. Отметим первую, обведем ее в рамку. Эту сроку назовем рабочей строкой. 1 1 1 1 2 1 1 2 . 1 1 2 1 Среди коэффициентов (чисел слева от вертикальной черты) отметим любой отличный от нуля. В нашем примере это любой из чисел 1, 1, 1. Пусть это будет второй коэффициент. Этот элемент назовем ведущим. Заключим в рамочку весь столбик, в котором находится ведущий элемент. 1 1 1 1 2 1 1 2 . 1 1 2 1 Теперь с помощью рабочей строки изменим все остальные с помощью преобразований 1–6) так, чтобы против ведущего элемента в отмеченном столбце получились нули. 1 1 1 1 1 1 1 1 -1 1 0 2 0 . -1 2 1 1 2 1 1 2 1 0 0 1 0 Здесь мы ко второй строке прибавили рабочую, умноженную на (-1), а к третьей строке прибавили рабочую, умноженную на (-1). Смотрим, не получилась ли в результате нулевая строка. Если получилась, то ее отбросим. Смотрим, не получились ли две одинаковые строки. Если «да», то одну из них отбросим. Смотрим, не получились ли две пропорциональные строки. Если «да», то одну из них отбросим. Смотрим, не имеют ли все элементы какой-либо строки общий множитель, отличный от нуля. Если «да», то эту строку на этот множитель разделим. На этом заканчивается первый шаг. 2 шаг. Среди строк последней таблицы, которые не были рабочими, отмечаем рабочую и в ней отмечаем ведущий элемент. С помощью преобразований 1–6) в столбце, в котором находится ведущий элемент, напротив ведущего элемента получаем нули. В нашем примере: 1 1 1 1 1 1 0 1 -1 1 0 0 0 . 1 0 2 0 -2 0 0 0 0 1 0 1 0 Из строк полученной матрицы отбрасываем нулевую (если она возникла), одну из двух одинаковых строк (если они возникли), одну из пропорциональных строк (если они возникли), произведем деление на общий множитель всех элементов строки (если такая возникла). На этом заканчивается второй шаг. И так далее до тех пор, пока все строки таблицы не побывают рабочими. В нашем примере: 1 1 1 1 1 1 1 1 1 1 0 1 -1 -1 -1 1 0 0 0 -1 1 0 2 0 2 1 1 2 -2 1 1 2 1 0 0 0 0 1 0 1 0 0 1 0 1 1 0 0 0 . 0 0 1 0 В итоге могут возникнуть следующие три случая. Случай 1. Среди строк итоговой матрицы найдется строка, в которой все элементы равны нулю, а свободный член отличен от нуля. Например, 0 0 0 5 . Восстановив по этой строке уравнение, мы получим: 0x1 0x 2 0 x 3 5 . Это уравнение явно не имеет решений. Т. е., в таком случае система несовместна. Случай 2. Итоговая таблица имеет столько же строк, сколько в системе имеется неизвестных. И ни одна из строк не показывает на несовместность системы. Например, 0 1 5 0 0 2 0 3 . 0 0 3 0 В этом случае по таблице восстанавливаем стандартную запись системы. В этой системе каждое уравнение будет содержать только по одному неизвестному с коэффициентом, отличным от нуля. Все остальные коэффициенты будут равны нулю. Например: 5 x1 0 x 2 0 x 3 1, 0 x1 2 x 2 0 x 3 3, 0 x 0 x 3 x 0. 1 2 3 5 x1 1, Обычно, члены с нулевыми коэффициентами не пишут: 2 x 2 3, 3x 0. 3 Разделим каждое уравнение получившейся системы на отличный от 1 x1 5 3 нуля коэффициент этого уравнения: x 2 . Получилась система 2 x 3 0 простейшего вида, эквивалентная исходной. Последняя система явно имеет 1 3 только одно решение ; ; 0 . 5 2 Обратимся к первоначальному примеру. 1 1 1 2 1 1 1 1 2 0 1 0 1 0 0 0 0 1 1 0 ; 0 1 2 1 -1 1 1 1 1 1 1 0 1 -1 -1 1 0 0 0 -1 1 0 2 0 -2 0 0 0 0 1 0 1 0 x 2 1 x1 0 . Система имеет только одно решение =(0; 1; 0), x 0 3 x1 x 2 x 3 1, значит и исходная система 2x 1 x 2 x 3 2, имеет только одно решение x x 2 x 1. 2 3 1 =(0;1;0). Случай 3. Итоговая таблица имеет строк меньше, чем число неизвестных. x1 x 2 x 3 x 4 1, Например: 2x1 x 2 3 x 3 x 4 2, 3x 2x 4 x 2x 3. 2 3 4 1 1 1 1 1 1 1 1 1 1 1 -1 1 1 1 1 1 2 1 3 1 2 1 0 2 0 1 1 0 2 0 1 -2 3 2 4 2 3 1 0 2 0 1 отбрасываем 0 1 1 1 0 . Последняя таблица имеет две строки, а неизвестных 1 0 2 0 1 системы четыре. Среди столбиков коэффициентов последней таблицы отмечаем столько, сколько строк в этой таблице. Именно, отмечаем такие столбики, которые после перестановки могут задать таблицу, в которой все элементы главной диагонали (элементы, стоящие по диагонали с верхнего угла к правому нижнему) были отличны от нуля, а все остальные элементы – нули: c1 0 , c 0 . i ck 0 В нашем примере этими столбиками могут быть первый и второй. Если 1 0 в уме их поменять местами, то получится 0 1 Неизвестные, коэффициенты которых попали в отмеченные столбики, называют главными. В нашем примере это х1 и х2. Все остальные неизвестные называются свободными. Теперь по последней таблице восстанавливаем стандартную запись системы так, чтобы слева от знаков “=” были члены только с главными неизвестными, а справа – только свободные члены и члены со свободными неизвестными. При этом: 1) если коэффициент главного неизвестного в уравнении равен нулю, то этот член писать не будем; 2) если коэффициент свободного неизвестного равен нулю, то этот член будем писать обязательно; 3) нумерация неизвестных в каждом уравнении слева от знака “=” возрастающая, справа – тоже возрастающая. В нашем примере: х1, х2 – главные неизвестные, х3, х4 – свободные неизвестные. x 2 0 x 3 x 4 , x1 1 2 x 3 0 x 4 . Если свободным неизвестным придать конкретные числовые значения, то согласно последней системе мы найдем соответствующие числовые значения главных неизвестных и тем самым найдем решение последней системы, а, следовательно, заданной. Например, 1 = (1-6+0; 0+3-5; 3; 5) = (-5; -2; 3; 5.) Это так называемое «частное решение» системы. Итак, частное решение системы получается при конкретных (частных) числовых заданиях свободных неизвестных систем. Т. к. свободных неизвестных у системы по крайней мере одно, и оно может принять любое числовое значение из R , то частных решений будет бесконечно много. Чтобы описать (задать) бесконечное множество, нужно указать характеристическое свойство его элементов. Это можно сделать следующим образом. Обозначим множество решений нашей системы М. Тогда М = { = (1-2с1+0с2; 0+с1-с2; с1 ;с2), с1, с2 R }. Решение системы записанное в виде = (1-2с1+0с2; 0+с1-с2; с1 ;с2) = (1-2с1; с1-с2;с1 ; с2). Называют общим решением системы. Другими словами, общим решение системы с бесконечным множеством решений называют решение, в котором свободным неизвестным приписаны значения, обозначенные буквами, в зависимости от которых произведена запись главных неизвестных. Заметим, т. к. выбор главных неизвестных, а следовательно, и свободных, может происходить неоднозначно, то вид записи общего решения тоже неоднозначен. Например, в нашем случае в качестве главных неизвестных можно взять х1 и х4. тогда свободными неизвестными будут х2 и х3. Итоговая система принимает вид: x 4 0 x 2 x 3 , x1 1 0 x 3 2 x 3 . И общее решение системы будет иметь вид: = (1+0с1-2с2; с1 ;с2; 0-с1+с2) = (1-2с2; с1; с2; -с1+с2), с1, с2 R . Лабораторная работа № 14. Метод Гаусса решения систем линейных уравнений Вопросы к работе 1. Какая система называется системой линейных уравнений? 2. Какие типы систем линейных уравнений? 3. Что такое решение системы n уравнений? 4. Что значит «решить» систему? 5. На какие свойства эквивалентных систем опирается метод Гаусса? 6. Как записать систему линейных уравнений в виде таблицы? 7. Какие преобразования таблицы системы линейных уравнений мы имеем право делать? 8. В чем заключается выполнение шага метода Гаусса? Сколько шагов в методе Гаусса? 9. В каком случае система линейных уравнений будет несовместна? 10.В каком случае система линейных уравнений будет иметь только одно решение? 11.Как найти общее решение неопределенной системы линейных уравнений? 12.Что такое частное решение неопределенной системы? Образцы решения заданий 1. Решить систему методом Гаусса. 3 х1 х2 х3 4 а) 2 х1 5 х2 3х3 17 х1 х2 х3 0 Решение. Записываем систему в виде таблицы и преобразуем эту таблицу: 3 1 1 4 2 5 3 17 1 1 1 0 0 0 1 0 1 1 1 0 0 1 0 1 0 2 0 1 0 1 0 0 0 1 0 0 4 4 4 0 1 1 1 0 1 1 1 7 1 17 0 7 1 17 0 8 0 16 1 1 0 1 1 1 0 1 0 0 1 3 2 1 Записываем систему в стандартном виде: х3 3 х 2 2 ; система имеет только одно решение 1 х1 α=(1,2,3). Проверим, удовлетворяет ли это решение исходной системе. 3 1 2 3 4 2 1 5 2 3 3 17 1 2 3 0 4 4(и ) 17 17(и ) Ответ: α=(1,2,3). 0 0(и ) зх1 х2 5 б) 2 х1 х2 х3 0 х1 х3 1 Решение. Записываем систему в виде таблицы: 3 1 0 5 3 1 0 5 0 0 0 4 2 1 1 0 3 1 0 1 3 0 0 1 0 1 1 1 0 1 1 1 0 1 1 1 восстанавливаем по таблице первое уравнение: 0х1 0х2 0х3 4 . Это уравнение решений не имеет. Значит и система решений не имеет. Следовательно, заданная система несовместима. 2 х1 х2 3х3 2 х4 4 в) 3х1 х2 х3 2 х4 6 х1 х2 7 х3 2 х4 2 Записываем таблицу заданной системы и преобразуем ее: 2 1 3 2 4 1 0 4 0 2 1 0 4 0 2 1 0 4 0 2 3 1 1 2 6 2 0 8 0 4 1 1 7 2 2 0 1 11 2 0 1 1 7 2 2 1 1 7 2 2 Главные неизвестные: х1 , х2 . Свободные неизвестные: х3 , х4 . Восстанавливаем систему: х1 2 4 х3 0 х4 ; х2 0 11х3 2 х4 х1 2 4 х3 0 х4 х2 0 11х3 2 х4 Общее решение системы: 2 4c1 0c2 ; 0 11c1 2c2 ; c1 ; c2 2 4c1 ; 11c1 2c2 ; c1 ; c2 , где c1 , c2 R 2(2 4c1) (11c1 2c2) 3c1 2c2 4 Проверка: 3(2 4c1) (11c1 2c2) c1 2c2 6 (2 4c1) (11c1 2c2) 7c1 2c2 2 4 8c1 11c1 2c2 3c1 2c2 4 6 12c1 11c1 2c2 c1 2c2 6 ; 2 4c1 11c1 2c2 7c1 2c2 2 4 4 6 6 2 2 Ответ: (2 4с1 ;11с1 2с2 ; с1 ; с2) . где c1 , c2 R Упражнения Решить систему линейных уравнений: 3х1 2 х2 х3 5 3) 2 х1 3х2 х3 1 2 х1 х2 3х3 11 х1 х2 х3 х4 1 6 2) х1 2 х2 3х3 4 х4 3х1 2 х2 х3 5 х4 2 х 5 х х 8 х 5 2 3 4 1 х1 2 х2 х3 4 4) 3х1 5 х2 3х3 1 4 х1 3х2 4 х3 0 5 х1 8 х2 х3 4 5) х1 2 х2 3х3 1 2 х1 3х2 2 х3 9 7 х1 4 х2 х3 13 6) 2 х1 3х2 х3 10 9 х1 х2 1 х1 х2 2 х3 х4 1 1) 2 х1 х2 х3 х4 5 х1 2 х2 х3 4 Индивидуальное задание. Решить систему. 2 х1 х2 3х3 7 1) а) 2 х1 3х2 х3 1 3х1 2 х2 х3 6 2 х1 х2 2 х3 3 2) а) х1 х2 2 х3 4 4 х1 х2 4 х3 3 3х1 х2 х3 12 3) а) х1 2 х2 4 х3 6 5 х1 х2 2 х3 3 2 х1 х2 3х3 4 4) а) х1 3х2 х3 11 х1 2 х2 2 х3 7 3х1 2 х2 4 х3 12 5) а) 3х1 4 х2 2 х3 6 2 х1 х2 х3 9 8 х1 3х2 6 х3 4 6) а) х1 х2 х3 2 4 х1 х2 3х3 5 4 х1 х2 3х3 9 7) а) х1 х2 х3 2 8 х1 3х2 6 х3 12 2 х1 3х2 4 х3 33 8) а) 7 х1 5х2 24 4 х1 11х3 39 2 х1 7 х2 3х3 х4 6 б) 3х1 5х2 2 х3 2 х4 4 9 х1 4 х2 х3 7 х4 2 2 х1 3х2 5 х3 7 х4 1 б) 4 х1 6 х2 2 х3 3х4 2 2 х1 3х2 11х3 15 х4 1 3х1 4 х2 х3 2 х4 3 б) 6 х1 8 х2 2 х3 5х4 7 9 х1 12 х2 3х3 10 х4 13 3х1 2 х2 5 х3 4 х4 2 б) 6 х1 4 х2 4 х3 3х4 3 9 х1 6 х2 3х3 2 х4 4 2 х1 х2 3х3 7 х4 5 б) 6 х1 3х2 х3 4 х4 7 4 х1 2 х2 14 х3 31х4 18 9 х1 3х2 5 х3 6 х4 4 б) 6 х1 2 х2 3х3 х4 5 3х1 3х2 3х3 14 х4 8 х1 5 х2 4 х3 3х4 1 б) 2 х1 х2 2 х3 х4 0 5 х1 3х2 8 х3 х4 1 х1 х2 х3 х4 х5 6 б) х1 3х2 х3 3х4 3х5 8 х1 х2 3х4 х5 4 2 х1 3х2 4 х3 12 9) а) 7 х1 5х2 х3 33 4 х1 х3 7 3х1 х2 х3 2 х4 4 б) х1 х2 х3 2 х4 1 х1 х2 х3 6 х4 6 х1 4 х2 х3 6 10) а) 5 х2 4 х3 20 3х1 2 х2 5 х3 22 х1 х2 2 х3 х4 0 б) 2 х1 х2 2 х3 х4 2 х1 2 х2 4 х3 2 х4 2 1. 2. 3. 4. 1. 2. 3. 4. 5. Задания для самоконтроля Можно ли в таблице системы линейных уравнений преобразовать столбики? Почему? Может ли быть однородная система линейных уравнений несовместимой? Может ли линейная система, в которой число уравнений меньше числа неизвестных иметь только одно решение? Может ли линейная система, в которой число уравнений больше числа неизвестных иметь только одно решение? Литература (Боровиков А. Н. Математическая геология – ее методика или методология? // Пути познания Земли, -М., 1971). Воробьев Н. Н. Роль теории игр в математизации знаний // Методологические проблемы кибернетики: Материалы к Всесоюзной конференции. Т. 1. – М., 1970) Гохман В. М., Гуревич Б. Л., Саушкин Ю. Г. Проблемы метагеографии //Математика в экономической географии. Вопросы географии. Ст.77. – М.: Мысль, 1968; Гохман В.М., Минц А.А., Преображенский В.С. Системный подход в географии // Теоретическая география. Вопросы географии. Сб.88. -М.: Мысль, 1971; Гуревич Б. М., Саушкин Ю. Г. Математический метод в географии // Вести московского университета. Серия 5, География. 1966 №1). Гохман В. М.. Гуревич Б.Л., Саушкин Ю. Г. Проблемы метагеографии // Математика в экономической географии”. Вопросы географии. Сб. 77. – М.: Мысль, 1968; Гохман В. М., Минц А.А., Преображенский В. С. Системный подход в географии. // Теоретическая география. Вопросы географии. Сб.88. – М.: Мысль, 1971; Согава В. Б. Определение некоторых понятий и терминов физической географии: Доклад института географии Сибири и Дальнего Востока. Вып. 3. – Иркутск, 1963; Согава В. Б. Структурно-динамическое ландшафтоведение и географические проблемы будущего: Доклад института географии Сибири и Дальнего Востока. Вып. 16. – Иркутск, 1967). 6. Саушкин Ю.Г. Смирнов А. М. Роль ленинских идей в развитии теоретической географии // Вести Московского университета. Серия 5 География 1970. №1 7. Гохман В.М., Гуревич Б.Л., Саушкин Ю.Т. Проблемы метагеографии // Математика в экономической географии. Вопросы географии сб.77-М.: Мысль, 1968 8. Гуревич Б.Л. Саушкин Ю.Г. Математический метод в географии // Вестн. Моск. ун-та. Сер. 5 Геогр 1966. №1