Решение линейных уравнений с помощью метода крамера. Правило Крамера. Метод обратной матрицы. Описание метода Крамера

Выявление и анализ тенденции временного ряда часто производится с помощью его выравнивания или сглаживания. Экспоненциальное сглаживание - один из простейших и распространенных приемов выравнивания ряда. Экспоненциальное сглаживание можно представить как фильтр, на вход которого последовательно поступают члены исходного ряда, а на выходе формируются текущие значения экспоненциальной средней.

Пусть - временной ряд.

Экспоненциальное сглаживание ряда осуществляется по рекуррентной формуле: , .

Чем меньше α, тем в большей степени фильтруются, подавляются колебания исходного ряда и шума.

Если последовательно использовать рекуррентное это соотношение, то экспоненциальную среднюю можно выразить через значения временного ряда X.

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

После появления работ Р. Брауна экспоненциальное сглаживание часто используется для решения задачи краткосрочного прогнозирования временных рядов.

Постановка задачи

Пусть задан временной ряд: .

Необходимо решить задачу прогнозирования временного ряда, т.е. найти

Горизонт прогнозирования, необходимо, чтобы

Для того, чтобы учитывать устаревание данных, введем невозрастающую последовательность весов , тогда

Модель Брауна

Предположим, что D - невелико (краткосрочный прогноз), то для решения такой задачи используют модель Брауна .

Если рассматривать прогноз на 1 шаг вперед, то - погрешность этого прогноза, а новый прогноз получается в результате корректировки предыдущего прогноза с учетом его ошибки - суть адаптации.

При краткосрочном прогнозировании желательно как можно быстрее отразить новые изменения и в то же время как можно лучше «очистить» ряд от случайных колебаний. Т.о. следует увеличивать вес более свежих наблюдений: .

С другой стороны, для сглаживания случайных отклонений, α нужно уменьшить: .

Т.о. эти два требования находятся в противоречии. Поиск компромиссного значения α составляет задачу оптимизации модели. Обычно, α берут из интервала (0,1/3).

Примеры

Работа экспоненциального сглаживания при α=0.2 на данных ежемесячных отчетов по продажам иностранной автомобильной марки в России за период с января 2007 по октябрь 2008. Отметим резкие падения в январе и феврале, когда продажи традиционно снижаются и повышения в начале лета.

Проблемы

Модель работает только при небольшом горизонте прогнозирования. Не учитываются тренд и сезонные изменения. Чтобы учесть их влияние, предлагается использовать модели: Хольта (учитывается линейный тренд) Хольта-Уинтерса (мультипликативные экспоненциальный тренд и сезонность), Тейла-Вейджа (аддетивные линейный тренд и сезонность).

Экспоненциальное сглаживание - способ сглаживания временных рядов, вычислительная процедура которого включает обработку всех предыдущих наблюдений, при этом учитывается устаревание информации по мере удаления от прогнозного периода. Иначе говоря, чем "старше" наблюдение, тем меньше оно должно влиять на величину прогнозной оценки. Идея экспоненциального сглаживания состоит в том, что по мере "старения" соответствующим наблюдениям придаются убывающие веса.

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

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

где I - 2, 3, 4 - периоды сглаживания; 5, - сглаженная величина на период £; У, - фактическое значение уровня на период 1 5, 1 - сглаженное значение на период Ь-Ьг- сглаженное значение тренда на период 1 - сглаженное значение на период I- 1; А и В - сглаживающие константы (числа между 0 и 1).

Сглаживающие константы А и В характеризуют фактор взвешивания наблюдений. Обычно Л, В < 0,3. Так как (1 - А) < 1, (1 - В) < 1, то они убывают по экспоненциальному закону по мере удаления наблюдения от текущего периода I. Отсюда данная процедура получила название экспоненциально сглаживания.

Уравнение добавляется в общую процедуру для сглаживания тренда. Каждая новая оценка тренда получается как взвешенная сумма разности между последними двумя сглаженными значениями (текущая оценка тренда) и предыдущей сглаженной оценки. Данное уравнение позволяет существенно сократить влияние случайных возмущений на тренд с течением времени.

Прогнозирование с использованием экспоненциального сглаживания подобно процедуре "наивного" прогнозирования, когда прогнозная оценка на завтра полагается равной сегодняшнему значению. В данном случае в качестве прогноза на один период вперед рассматривается сглаженная величина на текущий период плюс текущее сглаженное значение тренда:

Данную процедуру можно использовать для прогнозирования на любое число периодов, на пример на т периодов:

Процедура прогнозирования начинается с того, что сглаженная величина 51 полагается равной первому наблюдению У, т.е. 5, = У,.

Возникает проблема определения начального значения тренда 6]. Существуют два способа оценки Ьх.

Способ 1. Положим Ьх = 0. Такой подход хорошо работает в случае длинного исходного временного ряда. Тогда сглаженный тренд за небольшое число периодов приблизится к фактическому значению тренда.

Способ 2. Можно получить более точную оценку 6, используя первые пять (или более) наблюдений временного ряда. На их основе гю методу наименьших квадратов решается уравнение У(= а + Ь х г. Величина Ь берется в качестве начального значения тренда.

Методы Крамера и Гаусса – одни из самых популярных методов решения СЛАУ . К тому же, в ряде случаев целесообразно использовать именно конкретные методы. Сессия близка, и сейчас самое время повторить или освоить их с нуля. Сегодня разбираемся с решением методом Крамера. Ведь решение системы линейных уравнений методом Крамера - весьма полезный навык.

Системы линейных алгебраических уравнений

Система линейных алгебраических уравнений – система уравнений вида:

Набор значений x , при котором уравнения системы обращаются в тождества, называется решением системы, a и b – вещественные коэффициенты. Простенькую систему, состоящую из двух уравнений с двумя неизвестными, можно решить в уме либо выразив одну переменную через другую. Но переменных (иксов) в СЛАУ может быть гораздо больше двух, и здесь простыми школьными манипуляциями не обойтись. Что же делать? Например, решать СЛАУ методом Крамера!

Итак, пусть система состоит из n уравнений с n неизвестными.

Такую систему можно переписать в матричном виде

Здесь A – основная матрица системы, X и B , соответственно, матрицы-столбцы неизвестных переменных и свободных членов.

Решение СЛАУ методом Крамера

Если определитель главной матрицы не равен нулю (матрица невырожденная), систему можно решать по методу Крамера.

Согласно методу Крамера, решение находится по формулам:

Здесь дельта – определитель главной матрицы, а дельта x n-ное – определитель, полученный из определителя главной матрицы путем заменой n-ного столбца на столбец свободных членов.

В этом и заключается вся суть метода Крамера. Подставляя найденные по вышеприведенным формулам значения x в искомую систему, убеждаемся в правильности (или наоборот) нашего решения. Чтобы Вы быстрее уловили суть, приведем ниже пример подробного решения СЛАУ методом Крамера:

Даже если у Вас не получится с первого раза, не расстраивайтесь! Немного практики, и Вы начнете щелкать СЛАУ как орешки. Более того, сейчас совершенно необязательно корпеть над тетрадью, решая громоздкие выкладки и исписывая стержень. Можно легко решить СЛАУ методом Крамера в режиме онлайн, лишь подставив в готовую форму коэффициенты. Испробовать онлайн калькулятор решения методом Крамера можно, к примеру, на этом сайте .

А если система оказалась упорной и не сдается, Вы всегда можете обратиться за помощью к нашим авторам, например, чтобы . Будь в системе хоть 100 неизвестных, мы обязательно решим ее верно и точно в срок!

Рассмотрим систему 3-х уравнений с тремя неизвестными

Используя определители 3-го порядка, решение такой системы можно записать в таком же виде, как и для системы двух уравнений, т.е.

(2.4)

если 0. Здесь

Это есть правило Крамера решения системы трех линейных уравнений с тремя неизвестными .

Пример 2.3. Решить систему линейных уравнений при помощи правила Крамера:

Решение . Находим определитель основной матрицы системы

Поскольку 0, то для нахождения решения системы можно применить правило Крамера, но предварительно вычислим еще три определителя:

Проверка:

Следовательно, решение найдено правильно. 

Правила Крамера, полученные для линейных систем 2-го и 3-го порядка, наводят на мысль, что такие же правила можно сформулировать и для линейных систем любого порядка. Действительно имеет место

Теорема Крамера. Квадратная система линейных уравнений с отличным от нуля определителем основной матрицы системы (0) имеет одно и только одно решение и это решение вычисляется по формулам

(2.5)

где  – определитель основной матрицы ,  i определитель матрицы , полученной из основной, заменой i -го столбца столбцом свободных членов .

Отметим, что если =0, то правило Крамера не применимо. Это означает, что система либо не имеет вообще решений, либо имеет бесконечно много решений.

Сформулировав теорему Крамера, естественно возникает вопрос о вычислении определителей высших порядков.

2.4. Определители n-го порядка

Дополнительным минором M ij элемента a ij называется определитель, получаемый из данного путем вычеркивания i -й строки и j -го столбца. Алгебраическим дополнением A ij элемента a ij называется минор этого элемента, взятого со знаком (–1) i + j , т.е. A ij = (–1) i + j M ij .

Например, найдем миноры и алгебраические дополнения элементов a 23 и a 31 определителя

Получаем

Используя понятие алгебраического дополнения можно сформулировать теорему о разложении определителя n -го порядка по строке или столбцу .

Теорема 2.1. Определитель матрицы A равен сумме произведений всех элементов некоторой строки (или столбца) на их алгебраические дополнения:

(2.6)

Данная теорема лежит в основе одного из основных методов вычисления определителей, т.н. метода понижения порядка . В результате разложения определителя n -го порядка по какой-либо строке или столбцу, получается n определителей (n –1)-го порядка. Чтобы таких определителей было меньше, целесообразно выбирать ту строку или столбец, в которой больше всего нулей. На практике формулу разложения определителя обычно записывают в виде:

т.е. алгебраические дополнения записывают в явном виде через миноры.

Примеры 2.4. Вычислить определители, предварительно разложив их по какой-либо строке или столбцу. Обычно в таких случаях выбирают такой столбец или строку, в которой больше всего нулей. Выбранную строку или столбец будем обозначать стрелкой.

2.5. Основные свойства определителей

Разлагая определитель по какой-либо строке или столбцу, мы получим n определителей (n –1)-го порядка. Затем каждый из этих определителей (n –1)-го порядка также можно разложить в сумму определителей (n –2)-го порядка. Продолжая этот процесс, можно дойти до определителей 1-го порядка, т.е. до элементов матрицы, определитель которой вычисляется. Так, для вычисления определителей 2-го порядка придется вычислить сумму двух слагаемых, для определителей 3-го порядка – сумму 6 слагаемых, для определителей 4-го порядка – 24 слагаемых. Число слагаемых будет резко возрастать по мере увеличения порядка определителя. Это означает, что вычисление определителей очень высоких порядков становится довольно трудоемкой задачей, непосильной даже для ЭВМ. Однако вычислять определители можно и по-другому, используя свойства определителей.

Свойство 1 . Определитель не изменится, если в нем поменять местами строки и столбцы, т.е. при транспонировании матрицы :

.

Данное свойство свидетельствует о равноправии строк и столбцов определителя. Иначе говоря, любое утверждение о столбцах определителя справедливо и для его строк и наоборот.

Свойство 2 . Определитель меняет знак при перестановке двух строк (столбцов).

Следствие . Если определитель имеет две одинаковые строки (столбца), то он равен нулю.

Свойство 3 . Общий множитель всех элементов в какой-либо строке (столбце) можно вынести за знак определителя .

Например,

Следствие . Если все элементы некоторой строки (столбца) определителя равны нулю, то и сам определитель равен нулю .

Свойство 4 . Определитель не изменится, если к элементам одной строки (столбца), прибавить элементы другой строки (столбца), умноженной на какое-либо число .

Например,

Свойство 5 . Определитель произведения матриц равен произведению определителей матриц:

Метод Крамера или так называемое правило Крамера – это способ поиска неизвестных величин из систем уравнений. Его можно использовать только если число искомых значений эквивалентно количеству алгебраических уравнений в системе, то есть образуемая из системы основная матрица должна быть квадратной и не содержать нулевых строчек, а также если её детерминант не должен являться нулевым.

Теорема 1

Теорема Крамера Если главный определитель $D$ основной матрицы, составленной на основе коэффициентов уравнений, не равен нулю, то система уравнений совместна, причём решение у неё существует единственное. Решение такой системы вычисляется через так называемые формулы Крамера для решения систем линейных уравнений: $x_i = \frac{D_i}{D}$

В чем заключается метод Крамера

Суть метода Крамера в следующем:

  1. Чтобы найти решение системы методом Крамера, первым делом вычисляем главный определитель матрицы $D$. Когда вычисленный детерминант основной матрицы при подсчёте методом Крамера оказался равен нулю, то система не имеет ни одного решения или имеет нескончаемое количество решений. В этом случае для нахождения общего или какого-либо базисного ответа для системы рекомендуется применить метод Гаусса.
  2. Затем нужно заменить крайний столбец главной матрицы на столбец свободных членов и высчитать определитель $D_1$.
  3. Повторить то же самое для всех столбцов, получив определители от $D_1$ до $D_n$, где $n$ - номер крайнего справа столбца.
  4. После того как найдены все детерминанты $D_1$...$D_n$, можно высчитать неизвестные переменные по формуле $x_i = \frac{D_i}{D}$.

Приёмы для вычисления определителя матрицы

Для вычисления определителя матрицы с размерностью больше чем 2 на 2, можно использовать несколько способов:

  • Правило треугольников, или правило Саррюса, напоминающее это же правило. Суть метода треугольников в том, что при вычислении определителя произведения всех чисел, соединённых на рисунке красной линией справа, записываются со знаком плюс, а все числа, соединённые аналогичным образом на рисунке слева – со знаком минус. B то, и другое правило подходит для матриц размером 3 х 3. В случае же правила Саррюса сначала переписывается сама матрица, а рядом с ней рядом переписываются ещё раз её первый и второй столбец. Через матрицу и эти дополнительные столбцы проводятся диагонали, члены матрицы, лежащие на главной диагонали или на параллельной ей записываются со знаком плюс, а элементы, лежащие на побочной диагонали или параллельно ей - со знаком минус.

Рисунок 1. Правило треугольников для вычисления определителя для метода Крамера

  • С помощью метода, известного как метод Гаусса, также иногда этот метод называют понижением порядка определителя. В этом случае матрица преобразуется и приводится к треугольному виду, а затем перемножаются все числа, стоящие на главной диагонали. Следует помнить, что при таком поиске определителя нельзя домножать или делить строчки или столбцы на числа без вынесения их как множителя или делителя. В случае поиска определителя возможно только вычитать и складывать строки и столбы между собой, предварительно помножив вычитаемую строку на ненулевой множитель. Также при каждой перестановке строчек или столбцов матрицы местами следует помнить о необходимости смены конечного знака у матрицы.
  • При решении методом Крамера СЛАУ с 4 неизвестными, лучше всего будет применять именно метод Гаусса для поиска и нахождения определителей или опредлять детерминант через поиск миноров.

Решение систем уравнений методом Крамера

Применим метод Крамера для системы из 2 уравнений и двумя искомыми величинами:

$\begin{cases} a_1x_1 + a_2x_2 = b_1 \\ a_3x_1 + a_4x_2 = b_2 \\ \end{cases}$

Отобразим её в расширенной форме для удобства:

$A = \begin{array}{cc|c} a_1 & a_2 & b_1 \\ a_3 & a_4 & b_1 \\ \end{array}$

Найдём определитель основной матрицы, также называемый главным определителем системы:

$D = \begin{array}{|cc|} a_1 & a_2 \\ a_3 & a_4 \\ \end{array} = a_1 \cdot a_4 – a_3 \cdot a_2$

Если главный определитель не равен нулю, то для решения слау методом Крамера необходимо высчитать ещё парочку определителей от двух матриц с заменёнными столбцами основной матрицы на строчку свободных членов:

$D_1 = \begin{array}{|cc|} b_1 & a_2 \\ b_2 & a_4 \\ \end{array} = b_1 \cdot a_4 – b_2 \cdot a_4$

$D_2 = \begin{array}{|cc|} a_1 & b_1 \\ a_3 & b_2 \\ \end{array} = a_1 \cdot b_2 – a_3 \cdot b_1$

Теперь найдём неизвестные $x_1$ и $x_2$:

$x_1 = \frac {D_1}{D}$

$x_2 = \frac {D_2}{D}$

Пример 1

Метод Крамера для решения СЛАУ с основной матрицей 3 порядка (3 x 3) и тремя искомыми.

Решите систему уравнений:

$\begin{cases} 3x_1 – 2x_2 + 4x_3 = 21 \\ 3x_1 +4x_2 + 2x_3 = 9\\ 2x_1 – x_2 - x_3 = 10 \\ \end{cases}$

Сосчитаем главный детерминант матрицы пользуясь вышеизложенным под пунктом номер 1 правилом:

$D = \begin{array}{|ccc|} 3 & -2 & 4 \\3 & 4 & -2 \\ 2 & -1 & 1 \\ \end{array} = 3 \cdot 4 \cdot (-1) + 2 \cdot (-2) \cdot 2 + 4 \cdot 3 \cdot (-1) – 4 \cdot 4 \cdot 2 – 3 \cdot (-2) \cdot (-1) - (-1) \cdot 2 \cdot 3 = - 12 – 8 -12 -32 – 6 + 6 = - 64$

А теперь три других детерминанта:

$D_1 = \begin{array}{|ccc|} 21 & 2 & 4 \\ 9 & 4 & 2 \\ 10 & 1 & 1 \\ \end{array} = 21 \cdot 4 \cdot 1 + (-2) \cdot 2 \cdot 10 + 9 \cdot (-1) \cdot 4 – 4 \cdot 4 \cdot 10 – 9 \cdot (-2) \cdot (-1) - (-1) \cdot 2 \cdot 21 = - 84 – 40 – 36 – 160 – 18 + 42 = - 296$

$D_2 = \begin{array}{|ccc|} 3 & 21 & 4 \\3 & 9 & 2 \\ 2 & 10 & 1 \\ \end{array} = 3 \cdot 9 \cdot (- 1) + 3 \cdot 10 \cdot 4 + 21 \cdot 2 \cdot 2 – 4 \cdot 9 \cdot 2 – 21 \cdot 3 \cdot (-1) – 2 \cdot 10 \cdot 3 = - 27 + 120 + 84 – 72 + 63 – 60 = 108$

$D_3 = \begin{array}{|ccc|} 3 & -2 & 21 \\ 3 & 4 & 9 \\ 2 & 1 & 10 \\ \end{array} = 3 \cdot 4 \cdot 10 + 3 \cdot (-1) \cdot 21 + (-2) \cdot 9 \cdot 2 – 21 \cdot 4 \cdot 2 - (-2) \cdot 3 \cdot 10 - (-1) \cdot 9 \cdot 3 = 120 – 63 – 36 – 168 + 60 + 27 = - 60$

Найдём искомые величины:

$x_1 = \frac{D_1} {D} = \frac{- 296}{-64} = 4 \frac{5}{8}$

$x_2 = \frac{D_1} {D} = \frac{108} {-64} = - 1 \frac {11} {16}$

$x_3 = \frac{D_1} {D} = \frac{-60} {-64} = \frac {15} {16}$