Линейные неоднородные дифференциальные уравнения второго порядка с постоянными коэффициентами
Данная статья раскрывает вопрос о решении линейных неоднородных дифференциальных уравнений второго порядка с постоянными коэффициентами. Будет рассмотрена теория вместе с примерами приведенных задач. Для расшифровки непонятных терминов необходимо обращаться к теме об основных определениях и понятиях теории дифференциальных уравнений.
Рассмотрим линейное дифференциальное уравнение (ЛНДУ) второго порядка с постоянными коэффициентами вида y » + p · y ‘ + q · y = f ( x ) , где произвольными числами являются p и q , а имеющаяся функция f ( х ) непрерывная на интервале интегрирования x .
Перейдем к формулировке теоремы общего решения ЛНДУ.
Теорема общего решения ЛДНУ
Общим решением, находящимся на интервале х , неоднородного дифференциального уравнения вида y ( n ) + f n — 1 ( x ) · y ( n — 1 ) + . . . + f 0 ( x ) · y = f ( x ) с непрерывными коэффициентами интегрирования на x интервале f 0 ( x ) , f 1 ( x ) , . . . , f n — 1 ( x ) и непрерывной функцией f ( x ) равняется сумме общего решения y 0 , которое соответствует ЛОДУ и каким-нибудь частным решением y
, где исходным неоднородным уравнением является y = y 0 + y
Отсюда видно, что решение такого уравнения второго порядка имеет вид y = y 0 + y
. Алгоритм нахождения y 0 рассмотрен в статье о линейных однородных дифференциальных уравнениях второго порядка с постоянными коэффициентами. После чего следует переходить к определению y
Выбор частного решения ЛНДУ зависит от вида имеющейся функции f ( x ) , располагающейся в правой части уравнения. Для этого необходимо рассмотреть отдельно решения линейных неоднородных дифференциальных уравнений второго порядка при постоянных коэффициентах.
Когда f ( x ) считается за многочлен n -ой степени f ( x ) = P n ( x ) , отсюда следует, что частное решение ЛНДУ находим по формуле вида y
= Q n ( x ) · x γ , где Q n ( x ) является многочленом степени n , r – это количество нулевых корней характеристического уравнения. Значение y
является частным решением y
= f ( x ) , тогда имеющиеся коэффициенты, которые определены многочленом
Q n ( x ) , отыскиваем при помощи метода неопределенных коэффициентов из равенства y
Вычислить по теореме Коши y » — 2 y ‘ = x 2 + 1 , y ( 0 ) = 2 , y ‘ ( 0 ) = 1 4 .
Решение
Иначе говоря, необходимо перейти к частному решению линейного неоднородного дифференциального уравнения второго порядка с постоянными коэффициентами y » — 2 y ‘ = x 2 + 1 , которое будет удовлетворять заданным условиям y ( 0 ) = 2 , y ‘ ( 0 ) = 1 4 .
Общим решением линейного неоднородного уравнения является сумма общего решения, которое соответствует уравнению y 0 или частному решению неоднородного уравнения y
, то есть y = y 0 + y
Для начала найдем общее решение для ЛНДУ, а после чего – частное.
Перейдем к нахождению y 0 . Запись характеристического уравнения поможет найти корни. Получаем, что
k 2 — 2 k = 0 k ( k — 2 ) = 0 k 1 = 0 , k 2 = 2
Получили, что корни различные и действительные. Поэтому запишем
y 0 = C 1 e 0 x + C 2 e 2 x = C 1 + C 2 e 2 x .
. Видно, что правая часть заданного уравнения является многочленом второй степени, тогда один из корней равняется нулю. Отсюда получим, что частным решением для y
= Q 2 ( x ) · x γ = ( A x 2 + B x + C ) · x = A x 3 + B x 2 + C x , где значения А , В , С принимают неопределенные коэффициенты.
Найдем их из равенства вида y
Тогда получим, что:
‘ = x 2 + 1 ( A x 3 + B x 2 + C x ) » — 2 ( A x 3 + B x 2 + C x ) ‘ = x 2 + 1 3 A x 2 + 2 B x + C ‘ — 6 A x 2 — 4 B x — 2 C = x 2 + 1 6 A x + 2 B — 6 A x 2 — 4 B x — 2 C = x 2 + 1 — 6 A x 2 + x ( 6 A — 4 B ) + 2 B — 2 C = x 2 + 1
Приравняв коэффициенты с одинаковыми показателями степеней x , получим систему линейных выражений — 6 A = 1 6 A — 4 B = 0 2 B — 2 C = 1 . При решении любым из способов найдем коэффициенты и запишем: A = — 1 6 , B = — 1 4 , C = — 3 4 и y
= A x 3 + B x 2 + C x = — 1 6 x 3 — 1 4 x 2 — 3 4 x .
Эта запись называется общим решением исходного линейного неоднородного дифференциального уравнения второго порядка с постоянными коэффициентами.
Для нахождения частного решения, которое удовлетворяет условиям y ( 0 ) = 2 , y ‘ ( 0 ) = 1 4 , требуется определить значения C 1 и C 2 , исходя из равенства вида y = C 1 + C 2 e 2 x — 1 6 x 3 + 1 4 x 2 + 3 4 x .
y ( 0 ) = C 1 + C 2 e 2 x — 1 6 x 3 + 1 4 x 2 + 3 4 x x = 0 = C 1 + C 2 y ‘ ( 0 ) = C 1 + C 2 e 2 x — 1 6 x 3 + 1 4 x 2 + 3 4 x ‘ x = 0 = = 2 C 2 e 2 x — 1 2 x 2 + 1 2 x + 3 4 x = 0 = 2 C 2 — 3 4
Работаем с полученной системой уравнений вида C 1 + C 2 = 2 2 C 2 — 3 4 = 1 4 , где C 1 = 3 2 , C 2 = 1 2 .
Применив теорему Коши, имеем, что
y = C 1 + C 2 e 2 x — 1 6 x 3 + 1 4 x 2 + 3 4 x = = 3 2 + 1 2 e 2 x — 1 6 x 3 + 1 4 x 2 + 3 4 x
Ответ: 3 2 + 1 2 e 2 x — 1 6 x 3 + 1 4 x 2 + 3 4 x .
Когда функция f ( x ) представляется в виде произведения многочлена со степенью n и экспоненты f ( x ) = P n ( x ) · e a x , тогда отсюда получаем, что частным решением ЛНДУ второго порядка будет уравнение вида y
= e a x · Q n ( x ) · x γ , где Q n ( x ) является многочленом n -ой степени, а r – количеством корней характеристического уравнения, равняющиеся α .
Коэффициенты, принадлежащие Q n ( x ) находятся по равенству y
Найти общее решение дифференциального уравнения вида y » — 2 y ‘ = ( x 2 + 1 ) · e x .
Решение
Уравнение общего вида y = y 0 + y
. Указанное уравнение соответствует ЛОДУ y » — 2 y ‘ = 0 . По предыдущему примеру видно, что его корни равняются k 1 = 0 и k 2 = 2 и y 0 = C 1 + C 2 e 2 x по характеристическому уравнению.
Видно, что правой частью уравнения является x 2 + 1 · e x . Отсюда ЛНДУ находится через y
= e a x · Q n ( x ) · x γ , где Q n ( x ) , являющимся многочленом второй степени, где α = 1 и r = 0 , потому как у характеристического уравнения отсутствует корень, равный 1 . Отсюда получаем, что
= e a x · Q n ( x ) · x γ = e x · A x 2 + B x + C · x 0 = e x · A x 2 + B x + C .
А , В , С являются неизвестными коэффициентами, которые можно найти по равенству y
‘ = e x · A x 2 + B x + C ‘ = e x · A x 2 + B x + C + e x · 2 A x + B = = e x · A x 2 + x 2 A + B + B + C y
‘ ‘ = e x · A x 2 + x 2 A + B + B + C ‘ = = e x · A x 2 + x 2 A + B + B + C + e x · 2 A x + 2 A + B = = e x · A x 2 + x 4 A + B + 2 A + 2 B + C
‘ = ( x 2 + 1 ) · e x ⇔ e x · A x 2 + x 4 A + B + 2 A + 2 B + C — — 2 e x · A x 2 + x 2 A + B + B + C = x 2 + 1 · e x ⇔ e x · — A x 2 — B x + 2 A — C = ( x 2 + 1 ) · e x ⇔ — A x 2 — B x + 2 A — C = x 2 + 1 ⇔ — A x 2 — B x + 2 A — C = 1 · x 2 + 0 · x + 1
Показатели при одинаковых коэффициентах приравниваем и получаем систему линейных уравнений. Отсюда и находим А , В , С :
— A = 1 — B = 0 2 A — C = 1 ⇔ A = — 1 B = 0 C = — 3
Ответ: видно, что y
= e x · ( A x 2 + B x + C ) = e x · — x 2 + 0 · x — 3 = — e x · x 2 + 3 является частным решением ЛНДУ, а y = y 0 + y = C 1 e 2 x — e x · x 2 + 3 — общим решением для неоднородного дифуравнения второго порядка.
Когда функция записывается как f ( x ) = A 1 cos ( β x ) + B 1 sin β x , а А 1 и В 1 являются числами, тогда частным решением ЛНДУ считается уравнение вида y
= A cos β x + B sin β x · x γ , где А и В считаются неопределенными коэффициентами, а r числом комплексно сопряженных корней, относящихся к характеристическому уравнению, равняющимся ± i β . В этом случае поиск коэффициентов проводится по равенству y
Найти общее решение дифференциального уравнения вида y » + 4 y = cos ( 2 x ) + 3 sin ( 2 x ) .
Решение
Перед написанием характеристического уравнения находим y 0 . Тогда
k 2 + 4 = 0 k 2 = — 4 k 1 = 2 i , k 2 = — 2 i
Имеем пару комплексно сопряженных корней. Преобразуем и получим:
y 0 = e 0 · ( C 1 cos ( 2 x ) + C 2 sin ( 2 x ) ) = C 1 cos 2 x + C 2 sin ( 2 x )
Корни из характеристического уравнения считаются сопряженной парой ± 2 i , тогда f ( x ) = cos ( 2 x ) + 3 sin ( 2 x ) . Отсюда видно, что поиск y
будет производиться из y
= ( A cos ( β x ) + B sin ( β x ) · x γ = ( A cos ( 2 x ) + B sin ( 2 x ) ) · x . Неизвестные коэффициенты А и В будем искать из равенства вида y
= cos ( 2 x ) + 3 sin ( 2 x ) .
‘ = ( ( A cos ( 2 x ) + B sin ( 2 x ) · x ) ‘ = = ( — 2 A sin ( 2 x ) + 2 B cos ( 2 x ) ) · x + A cos ( 2 x ) + B sin ( 2 x ) y
» = ( ( — 2 A sin ( 2 x ) + 2 B cos ( 2 x ) ) · x + A cos ( 2 x ) + B sin ( 2 x ) ) ‘ = = ( — 4 A cos ( 2 x ) — 4 B sin ( 2 x ) ) · x — 2 A sin ( 2 x ) + 2 B cos ( 2 x ) — — 2 A sin ( 2 x ) + 2 B cos ( 2 x ) = = ( — 4 A cos ( 2 x ) — 4 B sin ( 2 x ) ) · x — 4 A sin ( 2 x ) + 4 B cos ( 2 x )
Тогда видно, что
= cos ( 2 x ) + 3 sin ( 2 x ) ⇔ ( — 4 A cos ( 2 x ) — 4 B sin ( 2 x ) ) · x — 4 A sin ( 2 x ) + 4 B cos ( 2 x ) + + 4 ( A cos ( 2 x ) + B sin ( 2 x ) ) · x = cos ( 2 x ) + 3 sin ( 2 x ) ⇔ — 4 A sin ( 2 x ) + 4 B cos ( 2 x ) = cos ( 2 x ) + 3 sin ( 2 x )
Необходимо приравнять коэффициенты синусов и косинусов. Получаем систему вида:
— 4 A = 3 4 B = 1 ⇔ A = — 3 4 B = 1 4
= ( A cos ( 2 x ) + B sin ( 2 x ) · x = — 3 4 cos ( 2 x ) + 1 4 sin ( 2 x ) · x .
Ответ: общим решением исходного ЛНДУ второго порядка с постоянными коэффициентами считается
= = C 1 cos ( 2 x ) + C 2 sin ( 2 x ) + — 3 4 cos ( 2 x ) + 1 4 sin ( 2 x ) · x
Когда f ( x ) = e a x · P n ( x ) sin ( β x ) + Q k ( x ) cos ( β x ) , тогда y
= e a x · ( L m ( x ) sin ( β x ) + N m ( x ) cos ( β x ) · x γ . Имеем, что r – это число комплексно сопряженных пар корней, относящихся к характеристическому уравнению, равняются α ± i β , где P n ( x ) , Q k ( x ) , L m ( x ) и N m ( x ) являются многочленами степени n , k , т , m , где m = m a x ( n , k ) . Нахождение коэффициентов L m ( x ) и N m ( x ) производится, исходя из равенства y
Найти общее решение y » + 3 y ‘ + 2 y = — e 3 x · ( ( 38 x + 45 ) sin ( 5 x ) + ( 8 x — 5 ) cos ( 5 x ) ) .
Решение
По условию видно, что
α = 3 , β = 5 , P n ( x ) = — 38 x — 45 , Q k ( x ) = — 8 x + 5 , n = 1 , k = 1
Тогда m = m a x ( n , k ) = 1 . Производим нахождение y 0 , предварительно записав характеристическое уравнение вида:
k 2 — 3 k + 2 = 0 D = 3 2 — 4 · 1 · 2 = 1 k 1 = 3 — 1 2 = 1 , k 2 = 3 + 1 2 = 2
Получили, что корни являются действительными и различными. Отсюда y 0 = C 1 e x + C 2 e 2 x . Далее необходимо искать общее решение, исходя из неоднородного уравнения y
= e α x · ( L m ( x ) sin ( β x ) + N m ( x ) cos ( β x ) · x γ = = e 3 x · ( ( A x + B ) cos ( 5 x ) + ( C x + D ) sin ( 5 x ) ) · x 0 = = e 3 x · ( ( A x + B ) cos ( 5 x ) + ( C x + D ) sin ( 5 x ) )
Известно, что А , В , С являются коэффициентами, r = 0 , потому как отсутствует пара сопряженных корней, относящихся к характеристическому уравнению с α ± i β = 3 ± 5 · i . Данные коэффициенты находим из полученного равенства:
= — e 3 x ( ( 38 x + 45 ) sin ( 5 x ) + ( 8 x — 5 ) cos ( 5 x ) ) ⇔ ( e 3 x ( ( A x + B ) cos ( 5 x ) + ( C x + D ) sin ( 5 x ) ) ) » — — 3 ( e 3 x ( ( A x + B ) cos ( 5 x ) + ( C x + D ) sin ( 5 x ) ) ) = — e 3 x ( ( 38 x + 45 ) sin ( 5 x ) + ( 8 x — 5 ) cos ( 5 x ) )
Нахождение производной и подобных слагаемых дает
— e 3 x · ( ( 15 A + 23 C ) · x · sin ( 5 x ) + + ( 10 A + 15 B — 3 C + 23 D ) · sin ( 5 x ) + + ( 23 A — 15 C ) · x · cos ( 5 x ) + ( — 3 A + 23 B — 10 C — 15 D ) · cos ( 5 x ) ) = = — e 3 x · ( 38 · x · sin ( 5 x ) + 45 · sin ( 5 x ) + + 8 · x · cos ( 5 x ) — 5 · cos ( 5 x ) )
После приравнивания коэффициентов получаем систему вида
15 A + 23 C = 38 10 A + 15 B — 3 C + 23 D = 45 23 A — 15 C = 8 — 3 A + 23 B — 10 C — 15 D = — 5 ⇔ A = 1 B = 1 C = 1 D = 1
Из всего следует, что
= e 3 x · ( ( A x + B ) cos ( 5 x ) + ( C x + D ) sin ( 5 x ) ) = = e 3 x · ( ( x + 1 ) cos ( 5 x ) + ( x + 1 ) sin ( 5 x ) )
Ответ: теперь получено общее решение заданного линейного уравнения:
= = C 1 e x + C 2 e 2 x + e 3 x · ( ( x + 1 ) cos ( 5 x ) + ( x + 1 ) sin ( 5 x ) )
Алгоритм решения ЛДНУ
Любой другой вид функции f ( x ) для решения предусматривает соблюдение алгоритма решения:
- нахождение общего решения соответствующего линейного однородного уравнения, где y 0 = C 1 ⋅ y 1 + C 2 ⋅ y 2 , где y 1 и y 2 являются линейно независимыми частными решениями ЛОДУ, С 1 и С 2 считаются произвольными постоянными;
- принятие в качестве общего решения ЛНДУ y = C 1 ( x ) ⋅ y 1 + C 2 ( x ) ⋅ y 2 ;
- определение производных функции через систему вида C 1 ‘ ( x ) + y 1 ( x ) + C 2 ‘ ( x ) · y 2 ( x ) = 0 C 1 ‘ ( x ) + y 1 ‘ ( x ) + C 2 ‘ ( x ) · y 2 ‘ ( x ) = f ( x ) , а нахождение функций C 1 ( x ) и C 2 ( x ) посредствам интегрирования.
Найти общее решение для y » + 36 y = 24 sin ( 6 x ) — 12 cos ( 6 x ) + 36 e 6 x .
Решение
Переходим к написанию характеристического уравнения, предварительно записав y 0 , y » + 36 y = 0 . Запишем и решим:
k 2 + 36 = 0 k 1 = 6 i , k 2 = — 6 i ⇒ y 0 = C 1 cos ( 6 x ) + C 2 sin ( 6 x ) ⇒ y 1 ( x ) = cos ( 6 x ) , y 2 ( x ) = sin ( 6 x )
Имеем, что запись общего решения заданного уравнения получит вид y = C 1 ( x ) · cos ( 6 x ) + C 2 ( x ) · sin ( 6 x ) . Необходимо перейти к определению производных функций C 1 ( x ) и C 2 ( x ) по системе с уравнениями:
C 1 ‘ ( x ) · cos ( 6 x ) + C 2 ‘ ( x ) · sin ( 6 x ) = 0 C 1 ‘ ( x ) · ( cos ( 6 x ) ) ‘ + C 2 ‘ ( x ) · ( sin ( 6 x ) ) ‘ = 0 ⇔ C 1 ‘ ( x ) · cos ( 6 x ) + C 2 ‘ ( x ) · sin ( 6 x ) = 0 C 1 ‘ ( x ) ( — 6 sin ( 6 x ) + C 2 ‘ ( x ) ( 6 cos ( 6 x ) ) = = 24 sin ( 6 x ) — 12 cos ( 6 x ) + 36 e 6 x
Необходимо произвести решение относительно C 1 ‘ ( x ) и C 2 ‘ ( x ) при помощи любого способа. Тогда запишем:
C 1 ‘ ( x ) = — 4 sin 2 ( 6 x ) + 2 sin ( 6 x ) cos ( 6 x ) — 6 e 6 x sin ( 6 x ) C 2 ‘ ( x ) = 4 sin ( 6 x ) cos ( 6 x ) — 2 cos 2 ( 6 x ) + 6 e 6 x cos ( 6 x )
Каждое из уравнений следует проинтегрировать . Тогда запишем получившиеся уравнения:
C 1 ( x ) = 1 3 sin ( 6 x ) cos ( 6 x ) — 2 x — 1 6 cos 2 ( 6 x ) + + 1 2 e 6 x cos ( 6 x ) — 1 2 e 6 x sin ( 6 x ) + C 3 C 2 ( x ) = — 1 6 sin ( 6 x ) cos ( 6 x ) — x — 1 3 cos 2 ( 6 x ) + + 1 2 e 6 x cos ( 6 x ) + 1 2 e 6 x sin ( 6 x ) + C 4
Отсюда следует, что общее решение будет иметь вид:
y = 1 3 sin ( 6 x ) cos ( 6 x ) — 2 x — 1 6 cos 2 ( 6 x ) + + 1 2 e 6 x cos ( 6 x ) — 1 2 e 6 x sin ( 6 x ) + C 3 · cos ( 6 x ) + + — 1 6 sin ( 6 x ) cos ( 6 x ) — x — 1 3 cos 2 ( 6 x ) + + 1 2 e 6 x cos ( 6 x ) + 1 2 e 6 x sin ( 6 x ) + C 4 · sin ( 6 x ) = = — 2 x · cos ( 6 x ) — x · sin ( 6 x ) — 1 6 cos ( 6 x ) + + 1 2 e 6 x + C 3 · cos ( 6 x ) + C 4 · sin ( 6 x )
Ответ: y = y 0 + y
= — 2 x · cos ( 6 x ) — x · sin ( 6 x ) — 1 6 cos ( 6 x ) + + 1 2 e 6 x + C 3 · cos ( 6 x ) + C 4 · sin ( 6 x )
Структура общего решения системы уравнений
Однородная система линейных уравнений
всегда совместна, так как имеет тривиальное решение . Если ранг матрицы системы равен количеству неизвестных , то тривиальное решение единственное. Предположим, что . Тогда однородная система имеет бесконечно много решений. Заметим, что расширенная матрица однородной системы при элементарных преобразованиях строк приводится к упрощенному виду , т.е. . Поэтому из (5.11) получаем общее решение однородной системы уравнений :
Получим другую форму записи решений однородной системы, которая раскрывает структуру множества решений. Для этого подчеркнем следующие свойства.
Свойства решений однородной системы уравнений
1. Если столбцы — решения однородной системы уравнений, то любая их линейная комбинация также является решением однородной системы.
В самом деле, из равенств следует, что
т.е. линейная комбинация решений является решением однородной системы.
2. Если ранг матрицы однородной системы равен , то система имеет линейно независимых решений.
Действительно, по формулам (5.13) общего решения однородной системы найдем частных решений , придавая свободным переменным следующие стандартные наборы значений (всякий раз полагая, что одна из свободных переменных равна единице, а остальные — равны нулю):
которые линейно независимы. В самом деле, если из этих столбцов составить матрицу, то последние ее строк образуют единичную матрицу. Следовательно, минор, расположенный в последних строках не равен нулю (он равен единице), т.е. является базисным. Поэтому ранг матрицы будет равен . Значит, все столбцы этой матрицы линейно независимы (см. теорему 3.4).
Любая совокупность линейно независимых решений однородной системы называется фундаментальной системой (совокупностью) решений .
Заметим, что фундаментальная система решений определяется неоднозначно. Однородная система может иметь разные фундаментальные системы решений, состоящие из одного и того же количества линейно независимых решений.
Теорема 5.3 об общем решении однородной системы. Если — фундаментальная система решений однородной системы уравнений (5.4), то столбец
при любых значениях произвольных постоянных также является решением системы (5.4), и, наоборот, для каждого решения х этой системы найдутся такие значения произвольных постоянных , при которых это решение удовлетворяет равенству (5.14).
Прямое утверждение теоремы следует из свойства 1 решений однородной системы. Докажем обратное утверждение о том, что любое решение можно представить в виде (5.14). Для этого составим матрицу , приписав к столбцам фундаментальной системы решений столбец
Найдем ранг этой матрицы. Так как первые столбцов линейно независимы, то . Так как каждый из столбцов матрицы является решением системы , то по первой формуле из (5.13) получаем
Следовательно, первая строка матрицы является линейной комбинацией последних строк этой матрицы.
По второй формуле из (5.13) получим, что вторая строка матрицы является линейной комбинацией последних строк этой матрицы, и т.д. По r-й формуле из (5.13) получим, что r-я строка матрицы является линейной комбинацией последних строк этой матрицы. Значит, первые строк матрицы можно вычеркнуть и при этом ранг матрицы не изменится. Следовательно, , так как после вычеркивания в матрице будет всего строк. Таким образом, . Значит, есть базисный минор матрицы , который расположен в первых ее столбцах, а столбец не входит в этот базисный минор. Тогда по теореме о базисном миноре найдутся такие числа , что
Итак, обратное утверждение доказано.
Алгоритм решения однородной системы уравнений
1-5. Выполнить первые 5 пунктов алгоритма Гаусса. При этом не требуется выяснять совместность системы, так как любая однородная система имеет решение (пункт 3 метода Гаусса следует пропустить). Получить формулы (5.11) общего решения, которые для однородной системы будут иметь вид (5.13).
Если ранг матрицы системы равен числу неизвестных , то система имеет единственное тривиальное решение и процесс решения заканчивается.
Если ранг матрицы системы меньше числа неизвестных , то система имеет бесконечно много решений. Структуру множества решений находим в следующих пунктах алгоритма.
6. Найти фундаментальную систему решений однородной системы. Для этого подставить в (5.13) последовательно стандартных наборов значений свободных переменных, в которых все свободные переменные равны нулю, кроме одной, равной единице (см. свойство 2 решений однородной системы).
7. Записать общее решение однородной системы по формуле (5.14).
1. В пункте 6 алгоритма вместо стандартного набора значений свободных переменных можно использовать и другие наборы значений, лишь бы они обеспечивали линейную независимость получаемых частных решений однородной системы.
2. Матрица столбцы которой образуют фундаментальную систему решений однородной системы, называется фундаментальной. Используя фундаментальную матрицу, общее решение (5.14) однородной системы можно записать в виде
3. Если базисный минор матрицы расположен в левом верхнем углу (в первых строках и первых столбцах), то упрощенный вид расширенной матрицы (5.9) однородной системы можно представить в виде блочной матрицы
Тогда блочная матрица размеров является фундаментальной. В этом можно убедиться, используя стандартные наборы значений свободных переменных. Применение блочных матриц может служить вторым способом нахождения фундаментальной системы решений.
Пример 5.4. Найти фундаментальную систему решений и общее решение однородной системы
Решение. 1. Составляем расширенную матрицу системы
2-4. Используя элементарные преобразования над строками матрицы , приводим ее к ступенчатому, а затем и к упрощенному виду (см. решение примера 5.3):
Пункт 3 метода Гаусса пропускаем.
5. Переменные — базисные, а — свободные. Записываем формулу (5.13) общего решения однородной системы
6. Находим фундаментальную систему решений. Так как и , надо подобрать линейно независимых решения. Подставляем в систему стандартные наборы значений свободных переменных:
В результате получили фундаментальную систему решений
7. Записываем общее решение однородной системы по формуле (5.14):
Заметим, что фундаментальную систему решений можно получить, взяв иные наборы значений свободных переменных. Например, и . Тогда получим другую фундаментальную систему решений
Несмотря на различия, обе формулы задают одно и то же множество решений.
Структура общего решения неоднородной системы уравнений
Ранее была выведена формула (5.11) общего решения системы линейных уравнений. Получим другую форму записи, отражающую структуру множества решений.
Рассмотрим неоднородную систему и соответствующую ей однородную систему . Между решениями этих систем имеются связи, выражающиеся следующими свойствами.
Свойства решений неоднородной системы уравнений
1. Разность двух решений и неоднородной системы есть решение однородной системы.
Действительно, из равенств и следует, что .
2. Пусть — решение неоднородной системы. Тогда любое решение неоднородной системы можно представить в виде
В самом деле, для любого решения неоднородной системы разность по свойству 1 является решением однородной системы, т.е. — решение однородной системы.
Теорема 5.4 о структуре общего решения неоднородной системы.
Пусть — решение неоднородной системы, а — фундаментальная система решений соответствующей однородной системы уравнений. Тогда столбец
при любых значениях [i]произвольных постоянных является решением неоднородной системы, и, наоборот, для каждого решения этой системы найдутся такие значения произвольных постоянных , при которых это решение удовлетворяет равенству (5.15).[/i]
Говорят, что общее решение неоднородной системы есть сумма частного решения неоднородной системы и общего решения соответствующей однородной системы.
Доказательство теоремы вытекает из свойств 1, 2 и теоремы 5.3.
Алгоритм решения неоднородной системы уравнений
1-5. Выполнить первые 5 пунктов метода Гаусса решения системы уравнений и получить формулу общего решения неоднородной системы вида (5.11).
6. Найти частное решение неоднородной системы, положив в (5.11) все свободные переменные равными нулю.
7. Записав формулы (5.13) общего решения соответствующей однородной системы, составить фундаментальную систему ее решений. Для этого подставить в (5.13) последовательно стандартных наборов значений свободных переменных, в которых все переменные равны нулю, за исключением одной, равной единице.
8. Записать общее решение неоднородной системы по формуле (5.15).
1. Используя фундаментальную матрицу однородной системы , решение неоднородной системы можно представить в виде
2. Если базисный минор матрицы расположен в левом верхнем углу (в первых строках и первых столбцах), то упрощенный вид расширенной матрицы (5.9) неоднородной системы можно представить в виде блочной матрицы
Тогда блочная матрица оказывается фундаментальной (см. п.3 замечаний 5.3), а столбец является частным решением неоднородной системы (в этом можно убедиться, подставляя в (5.11) нулевой набор свободных переменных). Используя блочные матрицы, общее решение (5 15) неоднородной системы можно представить в виде
где — столбец произвольных постоянных. Полученную формулу можно считать вторым способом решения неоднородной системы.
Пример 5.5. Найти структуру (5.15) общего решения неоднородной системы
Решение. 1-5. Первые 5 пунктов метода Гаусса выполнены при решении примера 5.3, где получены формулы общего решения неоднородной системы:
Переменные — базисные, а — свободные.
6. Полагая , получаем частное решение неоднородной системы .
7. Находим фундаментальную систему решений однородной системы (см. пример 5.4):
8. Записываем по формуле (5.15) общее решение неоднородной системы
Искомая структура множества решений найдена.
Получим формулу общего решения вторым способом , используя п.2 замечаний 5.4. При решении примера 5.3 расширенная матрица системы была приведена к упрощенному виду. Разбиваем ее на блоки:
Записываем частное решение неоднородной системы
и составляем фундаментальную матрицу:
По формуле (5.16) получаем общее решение неоднородной системы, которое преобразуем к виду (5.15):
Структура общего решения линейной неоднородной системы дифференциальных уравнений
достаточно найти общее решение однородной системы уравнений
и частное решение неоднородной системы уравнений (*).
Сумма этих решений представляет собой общее решение неоднородной системы уравнений (*).
В плане практического применения такой метод решения неоднородных систем линейных уравнений не слишком распространен. Однако подход подобного рода является весьма эффективным применительно к решению линейных неоднородных дифференциальных уравнений.
| AX = B | (1) |
может быть представлено в виде суммы общего решения соответствующей однородной системы
| AX = 0 | (2) |
и любого частного решения неоднородной системы (1).
Доказательство. Пусть матрица X 0 является общим решением однородной системы (2):
| AX 0 = 0. | (3) |
Обозначим символическим выражением X ч частное решение неоднородной системы (1):
| AXч = B. | (4) |
Складывая тождества (3) и (4), получам тождество
| A ( X0 + Xч ) = B, | (5) |
справедливое при любых значениях свободных параметров, входящих в общее решение X0 .
Следовательно, матрица X0 + Xч является общим решением матричного уравнения (1).
http://mathhelpplanet.com/static.php?p=struktura-obshchego-resheniya-sistemy-uravnenii
http://portal.tpu.ru/SHARED/k/KONVAL/Sites/Russian_sites/2/17a.htm