Решение линейных разностных уравнений с постоянными коэффициентами

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

Содержание:

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

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

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

Приращение функции при переходе от значения xi к значению xi+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 и т. д.

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

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

Откуда:

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Возьмем функцию и убедимся, что она будет решением уравнения (7.55). Поскольку , тогда . Подставим yt и yt-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.
Общим решением однородного уравнения является функция

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

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

Решим эту систему уравнений относительно A1 и A2:

Откуда

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Разностные уравнения для чайников

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

$$ a_0 y(x) + a_1 y(x+1) + a_2 y(x+2)=f(x), \\ a_0 y(x) + a_1 y(x+1) + a_2 y(x+2)+ a_3 y(x+3)=f(x). $$

Здесь $a_i$ — постоянные коэффициенты, $f(x)$ — правая часть (неоднородность уравнения), $y(x)$ — искомая неизвестная функция.

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

Примеры решений разностных уравнений

Задача 1. Решить разностное уравнение: $y(x+2)-4y(x+1)+4y(x)=2^x$

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

Задача 3. Решить разностное уравнение третьего порядка

$$ y(x+3)-16y(x+2)+83y(x+1)-140y(x)=0, \quad y(0)=3, y(1)=18, y(2)=120. $$

Задача 4. Найти частное решение однородного разностного уравнения:

Помощь с разностными уравнениями

Если вам нужна помощь с решением задач и контрольных по дифференциальным и разностным уравнениям (и другим разделам математического анализа), обращайтесь в МатБюро. Стоимость подробной консультации от 100 рублей , оформление производится в Word, срок от 1 дня.

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

Ниже разберем способы, как решить линейные однородные и неоднородные дифференциальные уравнения порядка выше второго, имеющих постоянные коэффициенты. Подобные уравнения представлены записями y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = 0 и y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = f ( x ) , в которых f 0 , f 1 , . . . , f n — 1 — являются действительными числами, а функция f ( x ) является непрерывной на интервале интегрирования X .

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

Общее решение ЛОДУ и ЛДНУ

Мы зададим формулировку двух теорем, показывающих, какого вида общих решений ЛОДУ и ЛНДУ n -ого порядка следует искать.

Общим решением y 0 ЛОДУ y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = 0 на интервале
X (коэффициенты f 0 ( x ) , f 1 ( x ) , . . . , f n — 1 ( x ) непрерывны на X ) будет линейная комбинация
n линейно независимых частных решений ЛОДУ y j , j = 1 , 2 , . . . , n , содержащая произвольные постоянные коэффициенты C j , j = 1 , 2 , . . . , n , то есть y 0 = ∑ j = 1 n C j · y j .

Общим решением y ЛНДУ y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = f ( x ) на интервале X (коэффициенты f 0 ( x ) , f 1 ( x ) , . . . , f n — 1 ( x ) непрерывны на X ) и функцией f ( x ) будет являться сумма y = y 0 + y

, где y 0 — общее решение соответствующего ЛОДУ y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = 0 , а y

— некоторое частное решение исходного ЛНДУ.

Итак, общее решение линейного неоднородного дифференциального уравнения, содержащего постоянные коэффициенты y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = f ( x ) , нужно искать, как y = y 0 + y

— некоторое его частное решение, а y 0 = ∑ j = 1 n C j · y j – общее решение соответствующего однородного дифференциального уравнения y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = 0 .

В первую очередь рассмотрим, как осуществлять нахождение y 0 = ∑ j = 1 n C j · y j — общее решение ЛОДУ n -ого порядка с постоянными коэффициентами, а потом научимся определять частное решение y

линейного неоднородного дифференциального уравнения n -ого порядка при постоянных коэффициентах.

Алгебраическое уравнение n -ого порядка k n + f n — 1 · k n — 1 + . . . + f 1 · k + f 0 = 0 носит название характеристического уравнения линейного однородного дифференциального уравнения n -ого порядка, содержащего постоянные коэффициенты, записи y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = 0 .

Возможно определить n частных линейно независимых решений y 1 , y 2 , . . . , y n исходного ЛОДУ, исходя из значений найденных n корней характеристического уравнения k 1 , k 2 , . . . , k n .

Методы решения ЛОДУ и ЛНДУ

Укажем все существующие варианты и приведем примеры на каждый.

  1. Когда все решения k 1 , k 2 , . . . , k n характеристического уравнения k n + f n — 1 · k n — 1 + . . . + f 1 · k + f 0 = 0 действительны и различны, линейно независимые частные решения будут выглядеть так:
    y 1 = e k 1 · x , y 2 = e k 2 · x , . . . , y n = e k n · x . Общее же решение ЛОДУ n -ого порядка при постоянных коэффициентах запишем как: y 0 = C 1 · e k 1 · x + C 2 · e k 2 · x + . . . + C n · e k n · x .

Пример 1

Задано ЛОДУ третьего порядка, содержащее постоянные коэффициенты y ‘ ‘ ‘ — 3 y » — y ‘ + 3 y = 0 . Определите его общее решение.

Решение

Cоставим характеристическое уравнение и найдем его корни, разложив предварительно многочлен из левой части равенства на множители, используя метод группировки:
k 3 — 3 k 2 — k + 3 = 0 k 2 ( k — 3 ) — ( k — 3 ) = 0 ( k 2 — 1 ) ( k — 3 ) = 0 k 1 = — 1 , k 2 = 1 , k 3 = 3

Ответ: найденные корни являются действительными и различными, значит общее решение ЛОДУ третьего порядка с постоянными коэффициентами запишем как: y 0 = C 1 · e — x + C 2 e x + C 3 · e 3 x .

  1. Когда решения характеристического уравнения являются действительными и одинаковыми ( k 1 = k 2 = . . . = k n = k 0 ) , линейно независимые частные решения линейного однородного дифференциального уравнения n -ого порядка с постоянными коэффициентами буду иметь вид: y 1 = e k 0 · x , y 2 = x · e k 0 · x , . . . , y n = x n — 1 · e k 0 · x .

Общее же решение ЛОДУ будет выглядеть так:
y 0 = C 1 · e k 0 · x + C 2 · e k 0 · x + . . . + C n · x n — 1 · e k 0 · x = = e k 0 · x · C 1 + C 2 · x + . . . + C n · x n — 1

Задано дифференциальное уравнение: y ( 4 ) — 8 k ( 3 ) + 24 y » — 32 y ‘ + 16 y = 0 . Необходимо определить его общее решение.

Решение

Составим характеристическое уравнение заданного ЛОДУ: k 4 — 8 k 3 + 24 k 2 — 32 k + 16 = 0 .

Преобразуем данное характеристическое уравнение, используя формулу бинома Ньютона, оно примет вид: k — 2 4 = 0 . Отсюда мы выделим его четырехкратный корень k 0 = 2 .

Ответ: общим решением заданного ЛОДУ станет: y 0 = e 2 x · C 1 + C 2 · x + C 3 · x 2 + C 4 · x 3

  1. Когда решения характеристического уравнения линейного однородного дифференциального уравнения n -ого порядка при постоянных коэффициентах — различные комплексно сопряженные пары α 1 ± i · β 1 , α 2 ± i · β 2 , . . . , α m ± i · β m , n = 2 m , линейно независимые частные решения такого ЛОДУ будут иметь вид:
    y 1 = e α 1 x · cos β 1 x , y 2 = e α 1 x · sin β 1 x , y 3 = e α 2 x · cos β 2 x , y 4 = e α 2 x · sin β 2 x , … y n — 1 = e α m x · cos β m x , y n = e α m x · sin β m x

Общее же решение запишем так:

y 0 = e α 1 x · C 1 · cos β 1 x + C 2 · sin β 1 x + + e α 2 x · C 3 · cos β 2 x + C 4 · sin β 2 x + . . . + + e α m x · C n — 1 · cos β m x + C n · sin β m x

Задано ЛОДУ четвертого порядка при постоянных коэффициентах y ( 4 ) — 6 y ( 3 ) + 14 y » — 6 y ‘ + 13 y = 0 . Необходимо его проинтегрировать.

Решение

Составим характеристическое уравнение заданного ЛОДУ: k 4 — 6 k 3 + 14 k 2 — 6 k + 13 = 0 . Осуществим преобразования и группировки:

k 4 — 6 k 3 + 14 k 2 — 6 k + 13 = 0 k 4 + k 2 — 6 k 3 + k + 13 k 2 + 1 = 0 k 2 + 1 k 2 — 6 k + 13 = 0

Из полученного результата несложно записать две пары комплексно сопряженных корней k 1 , 2 = ± i и k 3 , 4 = 3 ± 2 · i .

Ответ: общее решение заданного линейного однородного дифференциального уравнения n -ого порядка с постоянными коэффициентами запишется как:
y 0 = e 0 · C 1 · cos x + C 2 · sin x + e 3 x · C 3 · cos 2 x + C 4 · sin 2 x = = C 1 · cos x + C 2 · sin x + e 3 x · C 3 · cos 2 x + C 4 · sin 2 x

  1. Когда решения характеристического уравнения — это совпадающие комплексно сопряженные пары α ± i · β , линейно независимыми частными решениями линейного однородного дифференциального уравнения n-ого порядка с постоянными коэффициентами будут записи:
    y 1 = e α · x · cos β x , y 2 = e α · x · sin β x , y 3 = e α · x · x · cos β x , y 4 = e α · x · x · sin β x , … y n — 1 = e α · x · x m — 1 · cos β x , y n = e α · x · x m — 1 · sin β x

Общим решением ЛОДУ будет:

y 0 = e α · x · C 1 · cos β x + C 2 · sin β x + + e α · x · x · C 4 · cos β x + C 3 · sin β x + . . . + + e α · x · x m — 1 · C n — 1 · cos β x + C n · sin β x = = e α · x · cos β x · C 1 + C 3 · x + . . . + C n — 1 · x m — 1 + + e α · x · sin β x · C 2 + C 4 · x + . . . + C n · x m — 1

Задано линейное однородное дифференциальное уравнение с постоянными коэффициентами y ( 4 ) — 4 y ( 3 ) + 14 y » — 20 y ‘ + 25 y = 0 . Необходимо определить его общее решение.

Решение

Составим запись характеристического уравнения, заданного ЛОДУ, и определим его корни:

k 4 — 4 k 3 + 14 k 2 — 20 k + 25 = 0 k 4 — 4 k 3 + 4 k 2 + 10 k 2 — 20 k + 25 = 0 ( k 2 — 2 k ) 2 + 10 ( k 2 — 2 k ) + 25 = 0 ( k 2 — 2 k + 5 ) 2 = 0 D = — 2 2 — 4 · 1 · 5 = — 16 k 1 , 2 = k 3 , 4 = 2 ± — 16 2 = 1 ± 2 · i

Таким образом, решением характеристического уравнения будет двукратная комплексно сопряженная пара α ± β · i = 1 ± 2 · i .

Ответ: общее решение заданного ЛОДУ: y 0 = e x · cos 2 x · ( C 1 + C 3 · x ) + e x · sin 2 x · ( C 2 + C 4 · x )

  1. Встречаются различные комбинации указанных случаев: некоторые корни характеристического уравнения ЛОДУ n -ого порядка с постоянными коэффициентами являются действительными и различными, некоторые — действительными и совпадающими, а какие-то — комплексно сопряженными парами или совпадающими комплексно сопряженными парами.

Пример 5

Задано дифференциальное уравнение y ( 5 ) — 9 y ( 4 ) + 41 ( 3 ) + 35 y » — 424 y ‘ + 492 y = 0 . Необходимо определить его общее решение.

Решение

Составим характеристическое уравнение заданного ЛОДУ: k 5 — 9 k 4 + 41 k 3 + 35 k 2 — 424 k + 492 = 0 .

Левая часть содержит многочлен, который возможно разложить на множители. В числе делителей свободного члена определяем двукратный корень k 1 = k 2 = 2 и корень k 3 = — 3 .

На основе схемы Горнера получим разложение: k 5 — 9 k 4 + 41 k 3 + 35 k 2 — 424 k + 492 = k + 3 k — 2 2 k 2 — 8 k + 41 .

Квадратное уравнение k 2 — 8 k + 41 = 0 дает нам оставшиеся корни k 4 , 5 = 4 ± 5 · i .

Ответ: общим решением заданного ЛОДУ с постоянными коэффициентами будет: y 0 = e 2 x · C 1 + C 2 x + C 3 · e — 3 x + e 4 x · C 4 · cos 5 x + C 5 · sin 5 x

Таким образом, мы рассмотрели основные случаи, когда возможно определить y 0 — общее решение ЛОДУ n -ого порядка с постоянными коэффициентами.

Следующее, что мы разберем – это ответ на вопрос, как решить линейное неоднородное дифференциальное уравнение n -ого порядка с постоянными коэффициентами записи y ( n ) + f n — 1 · y ( n — 1 ) + . . . + f 1 · y ‘ + f 0 · y = f ( x ) .

Общее решение в таком случае составляется как сумма общего решения соответствующего ЛОДУ и частного решения исходного ЛНДУ: y = y 0 + y

. Поскольку мы уже умеем определять y 0 , остается разобраться с нахождением y

, т.е. частного решения ЛНДУ порядка n с постоянными коэффициентами.

Приведем все способы нахождения y

согласно тому, какой вид имеет функция f ( x ) , находящаяся в правой части рассматриваемого ЛНДУ.

    Когда f ( x ) представлена в виде многочлена n -ой степени f ( x ) = P n ( x ) , частным решением ЛНДУ станет: y

= Q n ( x ) · x γ . Здесь Q n ( x ) является многочленом степени n , а r – указывает, сколько корней характеристического уравнения равно нулю.
Когда функция f ( x ) представлена в виде произведения многочлена степени n и экспоненты f ( x ) = P n ( x ) · e α · x , частным решением ЛНДУ второго порядка станет: y

= e α · x · Q n ( x ) · x γ . Здесь Q n ( x ) является многочленом n —ой степени, r указывает, сколько корней характеристического уравнения равно α .
Когда функция f ( x ) записана как f ( x ) = A 1 cos ( β x ) + B 1 sin ( β x ) , где А 1 и В 1 – числа, частным решением ЛНДУ станет запись y

= A cos β x + B sin β x · x γ . Здесь где А и В являются неопределенными коэффициентами, r – указывает, сколько комплексно сопряженных пар корней характеристического уравнения равно ± i β .
Когда f ( x ) = e α x · P n ( x ) sin β x + Q k x cos β x , то y

= e α 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 = m a x ( n , k ) .

Коэффициенты, которые неизвестны, определяются из равенства y

( n — 1 ) + . . . + f 1 y

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

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

Пусть нам заданы y j , j = 1 , 2 , . . . , n — n линейно независимые частные решения соответствующего ЛОДУ, тогда, используя различные вариации произвольных постоянных, общим решением ЛНДУ
n -ого порядка с постоянными коэффициентами будет запись: н = ∑ j = 1 n C j ( x ) · y j . В нахождении производных функций C j ( x ) , j = 1 , 2 , . . . , n поможет система уравнений:

∑ j = 1 n C j ‘ ( x ) · y j = 0 ∑ j = 1 n C j ‘ ( x ) · y ‘ j = 0 ∑ j = 1 n C j ‘ ( x ) · y » j = 0 … ∑ j = 1 n C j ‘ ( x ) · y j ( n — 2 ) = 0 ∑ j = 1 n C j ‘ ( x ) · y j ( n — 1 ) = 0

а собственно функции C j ( x ) , j = 1 , 2 , . . . , n найдем при последующем интегрировании.

Задано ЛНДУ с постоянными коэффициентами: y ‘ ‘ ‘ — 5 y » + 6 y ‘ = 2 x . Необходимо найти его общее решение.

Решение

Составим характеристическое уравнение: k 3 — 5 k 2 + 6 k = 0 . Корни данного уравнения: k 1 = 0 , k 2 = 2 и k 3 = 3 . Таким образом, общим решением ЛОДУ будет запись: y 0 = C 1 + C 2 · e 2 x + C 3 · e 3 x , а частные линейно независимые решения это: y 1 = 1 , y 2 = e 2 x , y 3 = e 3 x .

Варьируем произвольные постоянные: y = C 1 ( x ) + C 2 ( x ) · e 2 x + C 3 ( x ) · e 3 x .

Чтобы определить C 1 ( x ) , C 2 ( x ) и C 3 ( x ) , составим систему уравнений:

C ‘ 1 ( x ) · y 1 + C ‘ 2 ( x ) · y 2 + C ‘ 3 ( x ) · y 3 = 0 C ‘ 1 ( x ) · y ‘ 1 + C ‘ 2 ( x ) · y ‘ 2 + C ‘ 3 ( x ) · y ‘ 3 = 0 C ‘ 1 ( x ) · y » 1 + C ‘ 2 ( x ) · y » 2 + C ‘ 3 ( x ) · y » 3 = 2 x ⇔ C ‘ 1 ( x ) · 1 + C ‘ 2 x · e 2 x ‘ + C ‘ 3 ( x ) · y 3 = 0 C ‘ 1 ( x ) · 1 ‘ + C ‘ 2 x · e 2 x ‘ + C ‘ 3 ( x ) · e 3 x ‘ = 0 C ‘ 1 ( x ) · 1 ‘ ‘ + C ‘ 2 x · e 2 x ‘ ‘ + C ‘ 3 ( x ) · e 3 x ‘ ‘ = 2 x ⇔ C ‘ 1 ( x ) · 1 + C ‘ 2 x · e 2 x + C ‘ 3 ( x ) · e 3 x = 0 C ‘ 1 ( x ) · 0 + C ‘ 2 ( x ) · 2 e 2 x + C ‘ 3 ( x ) · 3 e 3 x = 0 C ‘ 1 ( x ) · 0 + C ‘ 2 ( x ) · 4 e 2 x + C ‘ 3 ( x ) · 9 e 3 x = 2 x

Решаем, используя метод Крамера:

∆ = 1 e 2 x e 3 x 0 2 e 2 x 3 e 3 x 0 4 e 2 x 9 e 3 x = 18 e 2 x · e 3 x — 12 e 2 x · e 3 x = 6 e 5 x ∆ C 1 ‘ ( x ) = 0 e 2 x e 3 x 0 2 e 2 x 3 e 3 x 2 x 4 e 2 x 9 e 3 x = e 5 x · 2 x ⇒ C ‘ 1 ( x ) = ∆ C 1 ‘ ( x ) ∆ = e 5 x · 2 x 6 e 5 x = 1 6 · 2 x ∆ C 2 ‘ ( x ) = 1 0 e 3 x 0 0 3 e 3 x 0 2 x 9 e 3 x = — 3 e x · 2 x ⇒ C ‘ 2 ( x ) = ∆ C 2 ‘ ( x ) ∆ = — 3 e 3 x · 2 x 6 e 5 x = — 1 2 · e — 2 x · 2 x ∆ C 3 ‘ ( x ) = 1 e 2 x 0 0 2 e 2 x 0 0 4 e 2 x 2 x = 2 e 2 x · 2 x ⇒ C ‘ 3 ( x ) = ∆ C 3 ‘ ( x ) ∆ = 2 e 2 x · 2 x 6 e 5 x = 1 3 · e — 3 x · 2 x

Интегрируем C ‘ 1 ( x ) = 1 6 · 2 x с помощью таблицы первообразных, а
C ‘ 2 ( x ) = — 1 2 · e — 2 x · 2 x и C ‘ 3 ( x ) = 1 3 · e — 3 x · 2 x при помощи метода интегрирования по частям, получим:
C 1 ( x ) = 1 6 · ∫ 2 x d x = 1 6 · 2 x ln 2 + C 4 C 2 ( x ) = — 1 2 · ∫ e — 2 x · 2 x d x = — 1 2 · e — 2 x · 2 x ln 2 — 2 + C 5 C 3 ( x ) = 1 3 · ∫ e — 3 x · 2 x d x = 1 3 · e — 3 x · 2 x ln 2 — 3 + C 6

Ответ: искомым общим решением заданного ЛОДУ с постоянными коэффициентами будет:

y = C 1 ( x ) + C 2 ( x ) · e 2 x + C 3 ( x ) · e 3 x = = 1 6 · 2 x ln 2 + C 4 + — 1 2 · e — 2 x · 2 x ln 2 — 2 + C 5 · e 2 x + + 1 3 · e — 3 x · 2 x ln 2 — 3 + C 6 · e 3 x

где C 4 , C 5 и C 6 – произвольные постоянные.


источники:

http://www.matburo.ru/ex_ma.php?p1=marazn

http://zaochnik.com/spravochnik/matematika/differentsialnye-uravnenija/linejnye-differentsialnye-uravnenija-vysshih-por-1/