О частных решениях неоднородных обыкновенных линейных дифференциальных уравнений и их систем Текст научной статьи по специальности «Математика»

Вестник СамГУ — Естественнонаучная серия. 2005. №5(39).

УДК 517.9

91

О ЧАСТНЫХ РЕШЕНИЯХ НЕОДНОРОДНЫХ ОБЫКНОВЕННЫХ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ И ИХ СИСТЕМ1

© 2005 В.Г. Цирулик, Д.В.Цирулик2

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

1. Предварительные сведения

Многие задачи теории и практики приводят к необходимости интегрирования неоднородных обыкновенных линейных дифференциальных уравнений (ОЛДУ) и систем таких уравнений. Известно достаточно много методов, восходящих к работам Л. Эйлера, Ж. Лагранжа и О. Коши, позволяющих находить частные решения некоторых классов уравнений и систем в замкнутом виде. Укажем, к примеру, на методы неопределенных коэффициентов, метод исключения и метод вариации произвольных постоянных [1-10], метод функции отклика [1, 9], операционный метод [5]. В прикладных дисциплинах применяется еще метод комплексных амплитуд [3] и операторный метод [10].

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

1 Представлена доктором физико-математических наук профессором Л.М. Берковичем.

2Цирулик Владимир Григорьевич ([email protected]), Цирулик Дмитрий Владимирович ([email protected]), кафедра высшей математики Таганрогского государственного радиотехнического университета. 347922, Россия, г. Таганрог, пер. Некрасовский, 44.

п

Обозначим через Ьп = 2 х)ЬІ обыкновенный линейный дифференци-

І=0

альный оператор с переменными коэффициентами, где п — порядок оператора, 6 = ^ — оператор дифференцирования по переменной х.

Рассмотрим уравнение

Ьп (у) = /(х). (1.1)

В основе рассматриваемого метода лежит следующее утверждение [11]. Теорема 1.1. Пусть правая часть уравнения (1.1) представлена в виде

Г ___

/(х) = 2 сіиі(х), гДе не все числа Сі = 0, а функции щ(х), (і = 1 , г) образуют І=1

фундаментальную систему (ФСР) решений ОЛДУ Хг(у) = 0. Тогда, если оператор ХгЬп допускает факторизацию

ХгЬп = МП + Г-дУд, (1.2)

в которой операторы Ьп и Уд (0 < q) взаимно просты справа, то при д = = г частное решение уравнения (1.1) является линейной комбинацией ФСР уравнения Уд(у) = 0. Если же д < г, то возможно, что частное решение имеет указанный вид.

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

2. Частные решения коммутативно факторизуемых уравнений

Пусть ш(х) — достаточное число раз дифференцируемая вещественная функция. Обозначим через П = f любую ее первообразную, а через д

5д(П) = ^ яП — многочлен от П с, вообще говоря, комплексными коэф-І=0

фициентами. Через ш(х)6 обозначим оператор дифференцирования с умножением на функцию со(х) слева — со(х)^, тогда (соб)'—его і-я степень. Расмотрим дифференциальное уравнение

Ьт(у) = £ аі(шбУу = Бд(П), (2.1)

І=0

где аі — комплексные числа, а оператор Ьт несущественно отличается от общего вида коммутативно факторизуемых операторов [12]. В частности, при ш(х) = 1 уравнение (2.1)—с постоянными коэффициентами, при ш(х) = = х — уравнение Эйлера и т.д.

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

Лемма 2.1. Если 8д(П) — многочлен степени д, то (шб)д8д(П) Ф 0, а (шб)(д+;")5д(П) = 0 при ] > 0.

Доказательство. Действительно, шбП = 1, (шб)2П2 = 2!. Допустим, что (шб)д-1Пд-1 = (д - 1)!. Тогда (шб)дПд = (шб)д-1 (шбПд) = д(шб)д-1 Пд-1 = = д! .

т

Многочлен Н(Х) = 2 аіXі назовем характеристическим многочленом од-І=0

нородного уравнения Ьт(у) = 0.

Лемма 2.2. Справедливо тождество

т

Ьт(ие}'°а) = ех°а ^ -Я®(Х0)м®, (2.2)

І=0 І!

где Н(і)(Х0) — значения производных характеристического многочлена при X = Х0. Тождество выводится аналогично [1, с. 83].

Утверждения, аналогичные следующей теореме, имеются в [1—9], но там их доказательства проводятся методом неопределенных коэффициентов. Дадим доказательство, приводящее к иному методу отыскания частного решения, применив теорему 1.1 в том случае, когда правая часть уравне-

т

ния (1.1) является многочленом степени д от П, а оператор Ьт = 2 аі (шб)'.

І=0

В этом случае можно положить Хд+1 = (шб)д+1, и тождество (1.2) превратится в условие коммутирования операторов.

Теорема 2.1. Дифференциальное уравнение (2.1), где аі є С, І = 0,т, а0 Ф 0 и ат Ф 0, имеет частным решением многочлен ид(П) степени д, удовлетворяющий треугольной системе уравнений

((оЬуьт(у) = (шб)г5д(П), / = 0^, (2.3)

рассматриваемой как алгебраическая относительно переменных у, шбу, ..., (шб)ду.

Доказательство. Проведем доказательство индукцией по степени многочлена 8 д(П). При д = 0 утверждение очевидно. При д = 1 искомое частное решение можно определить согласно лемме 2.1 из алгебраической относительно переменных у, шбу системы: а1 шбу + а0у = 81(П), а0шбу = шб81(П). Допустим, что утверждение справедливо при д = к, то есть система

(со6)'Ьш(у) = (со6)'5д(П), / = ОД (2.4)

как алгебраическая относительно переменных у, шбу, ... , (шб)ку определяет многочлен степени к и его производные до к порядка включительно. Положим д = к + 1 и рассмотрим уравнение

Ьт(у) = 8 к+1(П). (2.5)

Применив к этому уравнению оператор шб, получим по лемме 2.1

!т(шбу) = S k (П). (2.6)

Согласно предположению индукции, уравнение (2.6) имеет частным решением многочлен

шбу = yk

степени &, удовлетворяющий вместе со своими производными системе (2.4), причем (шб)2у = у'1 , ..., (шб)&+1 у = у®.

Выразим из уравнения (2.5) у = ^(5^+1(П) - «хсобу - ... - ат((Х>Ъ)ту. Если заменить в этом выражении производную (соб)1у на (соб)'_1_у£, г = 1, иа,

ТО полученный многочлен у = тт-(5^+1(П) - а,1Ук - ... - йш(Сй6)Ш_1_У£ степени

и0

& +1 должен удовлетворять уравнению (2.5) по построению, при этом многочлен со6_у = ^(ш65£+1(П)-а1ш6_у£-...-яш(со6)ш_у£ совпадает с многочленом Уk, удовлетворяющим вместе со своими производными системе (2.4). Таким образом, многочлен у^+\ = ^(5^+1(П) - «хсобу - ... - ат((Х>Ъ)ту удовлетворяет

алгебраической системе ((х>Ь)1Ьп(у) = (со6)'5^+1(П), г = 0Д+1.

Замечание. Теорема 2.1 дает следующий способ отыскания полиномиального частного решения в рассматриваемом случае: применяя к исходному уравнению оператор шб с учетом тождества (шб)]у = 0, ] > q, составим и решим систему линейных уравнений (2.3) как алгебраическую относительно переменных у, шбу, ..., (шб)9у. Очевидно, что эта неоднородная система всегда является треугольной с отличными от нуля элементами на диагонали.

Следствие. Дифференциальное уравнение

т

Lm(y) = ^ а-(шб)'у = Р5(П),

1=г

где 0 ^ г < т, аг ф 0, а Р5(П) — многочлен степени 5, вообще говоря, с комплексными коэффициентами, имеет частным решением многочлен у = = Пгй5(П) степени 5 + г, где <25(П) — многочлен степени 5.

Доказательство. Данное уравнение заменой (шб)у = w сводится к уравнению (2.1) с отличным от нуля младшим коэффициентом.

Рассмотрим теперь уравнение с правой частью специального вида

Lm(y) = еаП(Рп(П) с»(ЬП) + Qk(П) 8ш(ЬП)), (2.7)

где а,-, а и Ь — вещественные числа, Рп(П), 2&(П) — многочлены с вещественными коэффициентами степени п и & соответственно.

Для отыскания частного решения уравнения (2.7) заменим его ком-плексифицированным уравнением, введя для этого новую неизвестную комплекснозначную функцию вещественного аргумента г = у + V где у, V — вещественные функции вещественного аргумента, ,2 = -1.

Уравнение

1т(г) = е%0П(Рп(П) - 12&(П)), (2.8)

где Х.0 = а + іЬ, называется комплексификацией уравнения (2.7).

Отделяя в уравнении (2.8) вещественную и мнимую части, легко убедиться, что его вещественная часть равносильна исходному вещественному уравнению.

Теорема 2.2. Уравнение (2.7) имеет частное решение вида

у = Пг ехр(аП)^е(и5(П)) cos(ЬП) - 1т(и5(П)) sm(ЬП)), (2.9)

где г — кратность числа Хо = а+іЬ как корня характеристического многочлена оператора Ьт, и5(П) — многочлен степени 5 = тах(п, к) с комплексными коэффициентами, который при г = 0 находится из системы (2.3), а при г > 0 совпадает с многочленом 25(П) из следствия; Re(Us(П)), 1т(и5(П)) — вещественная и мнимая части многочлена и5(х).

Доказательство. Сначала комплексифицируем уравнение (2.7), перейдя к уравнению Ьт(г) = Рп(П) - і2к(П). Воспользовавшись тождеством (2.2), выполним в этом уравнении замену переменной

г = ыеХоП, (2.10)

приведя его к виду

Ьт(и) = > -Я^ХпХсобУм =

т і

(и) = ^ — Я®(Х0)(со6)'м = Т,(Ц),

I!

1=Г

где коэффициенты Ы(1)(к0), вообще говоря, комплексные числа, Т5(П) = = Рп(П) - 12к(П).

Если в последнем уравнении г ф 0, то, сделав еще одну замену

(шб)ги = w,

приведем его к виду,

т-г

Lт-г(™) = ^ Б;(Х0)(шб)^ = Т3(П),

1=0

где ^0(^0) ф 0.

По теореме 2.1 последнее уравнение имеет частным решением многочлен w = Ж5(П). Интегрируя при г > 0 дифференциальное уравнение (шб)ги = Ж5(П) при нулевых произвольных постоянных, получим и = Пг25(П). Если г = 0, то и = Ж5(П). В силу подстановки (2.10) частное комплексное решение исходного уравнения имеет вид г = ПгУ5(П)ехр(.0П). Вещественное частное решение исходного уравнения находится по формуле у = \(г + г).

Приведем пример, иллюстрирующий обсуждаемый метод.

Пример 1. Найти частное решение уравнения х2у" + 3ху' + 2у = = х-1 [(1п2 х + 1) СОв (1п х) + (1п х - 2) 8Ш (1п х)].

Решение. Переписав уравнение в виде (хб)2у + 2хбу + 2у = ехр(- 1п х)Х Х[(1п2 х + 1)сов(1п х) + (1п х - 2)8т(1п х)], видим, что ш(х) = х, П = 1п х.

Комплексифицировав последнее уравнение, получим (хб)2г + 2хбг + 2г = = ехр(Хо1п х)(1п2 х + 1 - 11п х + 2г), где Хо = -1 + г. Выполнив замену переменной г = иеХо1пх, найдем (хб)2и + 21хЬы = 1п2 х + 1 - г 1п х + 2г. Положив

в последнем уравнении хби = w, получим xбw + 2г^ = 1п2 х + 1 - г 1п х + 2г. Частное решение этого уравнения является многочленом второй степени относительно 1п х. Поэтому, действуя на уравнение для w слева оператором хб с учетом условия (xб)Jw = 0 при ] > 2, составим систему

xбw + 2г^ = 1п2 х + 1 - г 1п х + 2г,

(хб)2w + 2ixбw = 21п х - г,

2г(хб)^ = 2,

из которой найдем т = ^ Х - ^ + 1- Проинтегрировав уравнение хЬи = т

при нулевых произвольных постоянных, получим и = 1пх^~г^ х - ^ + 1^.

В данном случае многочлен С/гСпх) = ^ х - ^ + 1. Поэтому по формуле

(2.9) вещественное частное решение исходного уравнения имеет вид у =

= х-1 1п X ^СОЗ(1п х) + + з)

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

3. Отыскание частных решений неоднородных систем

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

т

Цу) = ^МыЬУУ = 5 ,(£2), (3.1)

г=о

где А{ — пХп матрицы с комплексными элементами, у = (у\(х), ...,у„(х)) — столбец неизвестных, 8ц(П) = (8Ц1(П),..., 8 Цп(П)) — столбец правых частей, в котором 8Чк(П) — заданные многочлены степени Цк, число ц = шах(ц1,...,цп) назовем весом вектора (П).

Заметим, что в [3] дано решение аналогичных систем с постоянными коэффициентами, но с правыми частями более общего вида, методом исключения, а также изложен метод их нормализации. С некоторыми общими методами интегрирования линейных систем можно познакомиться, например, по монографии [13].

т

Матричный многочлен И(к) =2 А{к1 назовем характеристическим мно-

г=о

гочленом системы (3.1).

Лемма 3.1. Если (П) — столбец веса ц, то (шб)98ц(П) = 5 Ф 0, где

5 — постоянный вектор, (со6)|'9+^59(П) = 5 = 0 при ] > 0.

Доказательство аналогично доказательству леммы 2.1.

Теорема 3.1. Если матрица Ао обратима, то при условии совместности система дифференциальных уравнений (3.1) имеет частным решением столбец

Ф9(П) = (ФР1(П),...,ФРп(П)) (3.2)

веса ц, удовлетворяющий треугольной системе матричных уравнений

(шб)гЬ(у) = (шб)г Бд(П), г = о, ц, (3.3)

рассматриваемой как алгебраическая относительно неизвестных у шбу, ..., (шб)ц У.

Доказательство аналогично доказательству теоремы 2.1.

Пример 2. Найти частное решение системы

Аз/" + Агу” + Аі/ + Ас) = 5 2(х),

где -43 = ( 1 4 )’ Аг = ( 2 -2 )’ А‘ = ( 0 "/ )’ М = ( 2 5 )’ Ш) =

= (2х + 1, х2 - 2).

Решение. Согласно теореме 3.1, данная система имеет частным решением столбец веса два. Поэтому дифференцируя дважды исходную систему, получим следующую алгебраическую систему для отыскания частного решения:

Му" + Му' + Му = 52(П),

А і/' + АоУ = 65 2(П),

А0/' = б252(П).

Найденное из последней системы частное решение таково:

- / 1 2 14 13 1 2 „ 10

У = \-^х + т* _2х+ —

\ 3 3 9 3 9

Следствие. Если система

т

Цу) = ^М((0&Уу = БЧ(Ц), (3.4)

где о ^ г < т, 8ц(П) — столбец веса ц, совместна, а матрица Аг обратима, то она имеет частное решение вида

у = П Ф, (П), (3.5)

где Фц(П) — столбец веса ц.

Доказательство утверждения состоит в приведении системы (3.4) заменой (шб)г_у = V к виду (3.1) с последующим применением теоремы 3.1 и интегрированием системы (шб)г_у = V при нулевых произвольных постоянных.

Рассмотрим теперь систему более общего вида

L(y) = exp(aQ.)[S q(Q.) cos(bQ) + Rk(Q) sin(^n)], (3-6)

в котороой коэффициенты оператора L являются вещественными n X n матрицами, а и b — вещественные числа, Sq(H) и Rк(П) — столбцы многочленов с вещественными коэффициентами веса q и к соответственно. Назовем комплексификацией системы (3.6) следующую систему:

L{z) = ехр(Х0П)[59(П) - Шк(Щ (3.7)

где z = у + iv — новая неизвестная вектор-функция, _у и v — вещественные вектор-функции вещественного аргумента, Ло = а + ib. Очевидно, что вещественная часть последней системы совпадает с системой (3.6).

Теорема 3.2. Если система (3.6) совместна и при Ло = а + ib матрицы

Н^(ко), i = О, г — 1 нулевые, а матрица Н^Г\Хо) обратима, то система имеет

частное решение вида

у = £1Г ехр(аП)[Ие(С/^(П)) cos (6П) - 1т(С/^(П)) sin(&n))], (3-8)

где US(Q.) — столбец веса s = max(g, к), который при г = О находится из системы вида (3.3), а при r > 0 совпадает со столбцом Ф^П), отыскиваемым, как и в следствии; Re(Us(Q.)) и 1т(С/^(П)) — вещественная и мнимая части столбца Us(0).

Доказательство аналогично доказательству теоремы 2.2: комплексифи-цировав систему (3.6), выполним в полученной системе замену переменной z = иехр(ЛоП), приведя ее к виду

т -

Ш) = -Я®(Л0)(соб)гм = ТЖ), (3.9)

i=r

где коэффициенты Н^(ко), вообще говоря, комплексные матрицы, TS(D.) =

= р)-гЩ

Если в последней системе r Ф 0, то, сделав еще одну замену

(шб)ги = 1У,

приведем ее к виду

m-r

Lm-r(u) = BidoXubYw = Т,(Ц), i=0

где матрица Во(Ло) обратима по предположению. Далее по теореме 3.1 или по следствию получаем требуемое.

Заключение

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

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

Литература

[1] Степанов В.В. Курс дифференциальных уравнений. М.: Гостехиздат, 1950. 468 с.

[2] Матвеев Н.С. Методы интегрирования обыкновенных дифференциальных уравнений. М.: Высшая школа, 1963. 531 с.

[3] ПонтрягинЛ.С. Обыкновенные дифференциальные уравнения. М.: Наука, 1965. 331 с.

[4] КарташевА.П., Рождественский Б.Л. Обыкновенные дифференциальные уравнения и вариационное исчисление. М: Наука, 1976. 255 с.

[5] ФедорюкМ.В. Обыкновенные дифференциальные уравнения. М.: Наука, 1980. 350 с.

[6] Лизоркин П.И. Курс дифференциальных и интегральных уравнений с дополнительными главами анализа. М.: Наука, 1981. 381 с.

[7] ТихоновА.Н., ВасильеваА.Б., СвешниковА.Г. Дифференциальные уравнения. М.: Наука, 1985. 231 с.

[8] БибиковЮ.Н. Курс обыкновенных дифференциальных уравнений. М.: Высшая школа, 1991. 300 с.

[9] Боярчук А.К., Головач Г.П. Справочное пособие по высшей математике. Дифференциальные уравнения в примерах и задачах. Ч. 5. М.: Эди-ториал, 2001. 383 с.

[10] ЭльсгольцЛ.Э. Обыкновенные дифференциальные уравнения. СПб.: Лань, 2002. 218 с.

[11] ЦируликВ.Г. О методе неопределенных коэффициентов для неоднородных обыкновенных линейных дифференциальных уравнений // Диф-ференц. уравнения. 2003. №39(11). С. 1571-1573.

[12] Беркович Л.М. Факторизация и преобразования дифференциальных уравнений. Методы и приложения. М.: НИЦ ’’Регулярная и хаотическая динамика”, 2002. 464 с.

[13] Гантмахер Ф.Р. Теория матриц. М.: Наука, 1988. 548 с.

Поступила в редакцию 22/ VI /2005; в окончательном варианте—14/^1/2005.

SOLUTIONS OF INHOMOGENEOUS ORDINARY DIFFERENTIAL EQUATIONS AND SYSTEMS3

© 2005 V.G. Tsiroulik, D.V. Tsiroulik4

In the paper the method of finding solutions of some class of inhomo-geneous ordinary differential equations and systems with right-hand side of a special form is proposed. It is proved that for equations and systems of equations of this class the proposed method considerably reduces the amount of calculations in comparison with the known methods especially in the resonance case when right-hand side contains the polynomial of a high degree.

Paper received 22/ VI /2005. Paper accepted 14/VII/2005.

3Communicated by Dr. Sci. (Phys. & Math.) Prof. L.M. Berkovich.

4Tsiroulik Vladimir Grigorievich ([email protected]), Tsiroulik Dmitri Vladimirovich ([email protected]), Dept. of Mathematics, Taganrog State Radio Engineering University, Taganrog, 347922, Russia.