Общее решение неоднородного уравнения теорема

Линейные неоднородные дифференциальные уравнения второго порядка с постоянными коэффициентами

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

Рассмотрим линейное дифференциальное уравнение (ЛНДУ) второго порядка с постоянными коэффициентами вида 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 )

Разностные уравнения

Содержание:

Разностные уравнения

Понятие разницы и разностного уравнения

Если для значений переменной x₁, x₂, x₃, . функция f (x) принимает значения f (x₁), f (x₂), f (x₃) . , то приращения функции составляют f (x₂) – f (x₁), f (x₃) – f (x₂), .

Приращение функции при переходе от значения x_i к значению x_i+1 будем обозначать: В частности можно взять в качестве значения независимых переменных x и x + 1 . Разность Δf (x) = f (x + 1) — f (x) называется первой разностью или разностью первого порядка. Она может рассматриваться в свою очередь как функция от x, а потому и для нее можно определить разницу:

Введем обозначения ΔΔf (x) = Δ 2 f (x), тогда Δ 2 f (x) = f (x + 2) — 2 f (x + 1) + f (x) и называется разностью второго порядка.

Аналогично можно найти разности третьего, четвертого и т. д. порядков.

Определим разности некоторых важнейших функций.

1) Если f (x) = С, где С — постоянная величина, то
Δf (x) = f (x + 1) – f (x) = С – С = 0.

Понятно, что и все разности следующих порядков будут также равняться нулю.

2) Если f (x) = ax + b, то
Δf = Δf (x + 1) — f (x) = a (x + 1) + b — ax — b = a.

Разница первого порядка линейной функции равна постоянной, а все остальные будут равны нулю.

3) Если f (x) = ax 2 + bx + c, то

Поскольку разница первого порядка является линейной функцией, то разница второго порядка — постоянная, а все последующие разности равны нулю.

4) Если f (x) = a x , то

В экономических исследованиях часто встречаются задачи, в которых время t является независимой переменной, а зависимая переменная определяется для времени t, t + 1, t + 2 и т. д.

Обозначим y_t — значение функции y в момент времени t; y_t+1 — значение функции в момент, сдвинутый на одну единицу, например, на следующий час, на следующую неделю и т. д., y_t+2 — значение функции y в момент, сдвинутый на две единицы и т. д.

Очевидно, что

Откуда:

За разность второго порядка, имеем или
поэтому

Аналогично можно доказать, что

Итак, любую функцию

можно представить в виде: (7.50)
и наоборот.

Определение. Уравнение
(7.51)
называется разностным уравнением n-го порядка.

Решить разностное уравнение n-го порядка — это значит найти такую ​​функцию y_t, которая превращает уравнение (7.50) или (7.51) в тождество.

Решение, в котором есть произвольная постоянная, называется общим; решение, в котором постоянная отсутствует, называется частным.

Определение. Уравнение
(7.52)
где a₀, a₁, . a_n — постоянные числа, называется неоднородным разностным
уравнением n-го порядка с постоянными коэффициентами.

Если в уравнении (7.52) f (t) = 0, то уравнение называется однородным разностным уравнением n-го порядка с постоянными коэффициентами:
(7.53)

Уравнение есть однородное разностное уравнение первого порядка с постоянными коэффициентами a и b, а уравнение — неоднородное разностное уравнение второго порядка с постоянными коэффициентами a, b, c.

ТЕОРЕМА 1. Если решениями однородного разностного уравнения (7.53) является y₁ (t) и y₂ (t), то его решением будет также функция y₁ (t) + y₂ (t).

ТЕОРЕМА 2. Если y (t) является решением однородного разностного уравнения (7.53), то его решением будет также функция Ay (t), где А — произвольная постоянная.

ТЕОРЕМА 3. Если y (t) — частное решение неоднородного уравнения (7.52) и y (t, A₁, A₂, . A_n) — общее решение однородного уравнения (7.53), то общим решением неоднородного разностного уравнения будет функция: y (t) + y (t, A₁, A₂, . A_n).

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

Разностные уравнения первого порядка с постоянными коэффициентами

Рассмотрим неоднородное разностное уравнение
(7.54)

Соответствующее ему однородное уравнение будет:
(7.55)

Возьмем функцию и убедимся, что она будет решением уравнения (7.55). Поскольку , тогда . Подставим y_t и y_t-1 в уравнение (7.55):
Итак, является решением уравнения (7.55).

По теореме (2) общее решение однородного разностного уравнения (7.55) является функция , где А — произвольная постоянная.

Пусть — частное решение неоднородного разностного уравнения (7.54). По теореме (3) общим решением неоднородного разностного уравнения (7.54) будет функция

Частное решение найти нетрудно, если f (t) = α, где α — некоторая постоянная. На самом деле, если где u — постоянная. Подставим в уравнение (7.54), имеем: u — au = α, откуда
Итак, общее решение уравнения (7.54) запишем в виде: .

Разностные уравнения второго порядка с постоянными коэффициентами

Пусть задано неоднородное разностное уравнение второго порядка с постоянными коэффициентами:
(7.56)
и соответствующее ему однородное уравнение
(7.57)

Убедимся, что функция будет решением уравнения (7.58). Подставим в уравнение (7.57) (λ ≠ 0), получим Поскольку λ ≠ 0, то поделим на λ t-2 , имеем λ 2 + aλ + b = 0 (7.58)

Это уравнение называется характеристическим уравнением для уравнения (7.57).

Здесь могут иметь место следующие три случая:

1. D = a 2 – 4b > 0, тогда уравнение (7.58) будет иметь два действительных различных корня.
Общее решение уравнения (7.57) запишется в виде:

а общее решение неоднородного уравнения (7.56) запишется так:

2. D = a 2 – 4b = 0, тогда и и

В этом случае однородное уравнение (7.57) примет вид:
(7.59)
Тогда

Легко убедиться, что решением уравнения (7.59) является также функция
Поэтому общим решением уравнения (7.59) является функция а общим решением неоднородного уравнения (7.56) функция

3. D = a 2 – 4b 2 – 5λ + 6 = 0 будет иметь действительные разные корни (D = 25 – 24 = 1 > 0), λ1 =2, λ2 = 3.
Общим решением однородного уравнения является функция

Далее положим, что y_t = y — частное решение неоднородного уравнения, тогда

Таким образом, общим решением неоднородного уравнения является функция Постоянные A₁ и A₂ определим из начальных условий: y₀ = 5, y₁ = 9. Тогда для t = 0 и t = 1 соответственно будем иметь:

Решим эту систему уравнений относительно A₁ и A₂:

Откуда

Итак, — общее решение заданного в условии разностного уравнения.

Примеры применения разностных уравнений в экономических задачах

Пример 1. Пусть некоторая сумма средств выдается под сложный процент p, то к концу t-го года ее размер будет составлять:
Это однородное разностное уравнение первого порядка. Его решением будет функция , где A — некоторая постоянная, которую можно найти из начальных условий.

Если положить y₀ = F , то A = F, откуда

Это известная формула величины фонда F, который выдается под сложный процент.

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

Цена равновесия для данной продукции определяется равенством:
а это разностное уравнение первого порядка.

Положим, что функция спроса определяется формулой а функция предложения — формулой

Цена равновесия запишется: то есть Решением этого уравнения является функция Постоянная A определяется из начальных условий, для t = 0 цена составляет p₀.

Тогда p₀ = A и решением уравнения является функция
Если начальная цена p₀ = 0, то p_t = 0 для всех значений t.

Следовательно, цена не подлежит изменению.

Вообще говоря, функция предложения — возрастающая, а потому b > 0; а функция спроса — убывающая, и поэтому a

Присылайте задания в любое время дня и ночи в ➔

Официальный сайт Брильёновой Натальи Валерьевны преподавателя кафедры информатики и электроники Екатеринбургского государственного института.

Все авторские права на размещённые материалы сохранены за правообладателями этих материалов. Любое коммерческое и/или иное использование кроме предварительного ознакомления материалов сайта natalibrilenova.ru запрещено. Публикация и распространение размещённых материалов не преследует за собой коммерческой и/или любой другой выгоды.

Сайт предназначен для облегчения образовательного путешествия студентам очникам и заочникам по вопросам обучения . Наталья Брильёнова не предлагает и не оказывает товары и услуги.

Структура общего решения системы уравнений

Однородная система линейных уравнений

всегда совместна, так как имеет тривиальное решение . Если ранг матрицы системы равен количеству неизвестных , то тривиальное решение единственное. Предположим, что . Тогда однородная система имеет бесконечно много решений. Заметим, что расширенная матрица однородной системы при элементарных преобразованиях строк приводится к упрощенному виду , т.е. . Поэтому из (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):