Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

СеверинВП ДУ высших порядков

.pdf
Скачиваний:
45
Добавлен:
02.02.2015
Размер:
2 Мб
Скачать
☆

именно положим

 

y C1(x) y1 C2 (x) y2 ... Cn (x) yn ,

(6.6)

где C1(x) , C2 (x) , …, Cn (x) уже не произвольные постоянные, как в (6.5), а неизвестные пока функции. Подставляя выражение (6.6) в не-

однородное уравнение

L[ y] f (x) ,

получим для нахождения

C1(x) ,

C2 (x) , …,

Cn (x)

лишь одно уравнение. Для однозначного определе-

ния этих функций следует наложить на них

 

n 1 дополнительных

условий. Покажем, как это сделать самым удобным способом.

 

 

 

Дифференцируя равенство (6.6), имеем

 

 

 

 

 

 

 

 

y C (x) y

C (x) y

... C (x) y

C

(x) y C

(x) y

2

... C

(x) y

n

.

1

1

2

2

 

 

n

 

n

1

1

 

2

 

n

 

 

Потребуем теперь, чтобы не только

y имело такой же вид, как u , но

чтобы и y

выглядело так же,

как и u , то есть, чтобы было

 

 

 

 

 

 

y C (x) y C (x) y ... C (x) y .

 

 

(6.7)

 

 

 

 

1

1

2

 

 

2

n

 

 

n

 

 

 

 

 

Для этого нужно потребовать выполнение равенства

 

 

 

 

 

 

 

 

C (x) y C

(x) y

2

... C (x) y

n

0 .

 

 

 

 

 

 

 

 

1

1

2

 

 

 

n

 

 

 

 

 

 

 

Это и есть первое из n 1 дополнительных условий.

 

 

 

 

 

Дифференцируя равенство (6.7), получим

 

 

 

 

 

 

 

 

y C (x) y C (x) y

... C (x) y

C

(x) y C

(x) y

... C

(x) y .

1

1

2

2

 

 

n

 

n

1

1

 

2

 

2

n

 

n

 

Потребуем, чтобы y выглядело так же, как и

u , то есть, чтобы было

y C (x) y

C (x) y ... C

n

(x) y .

 

1

1

2

2

 

n

Для этого нужно положить

 

 

 

 

 

C

(x) y

C

(x) y

... C (x) y

0 .

1

1

2

2

n

 

n

 

Это и есть второе дополнительное условие. Аналогично, требуя, чтобы y имело тот же вид, что и u , получим третье условие

71

 

C (x) y C

(x) y

... C

(x) y 0 .

 

 

1

1 2

 

2

n

n

 

 

 

Продолжая этот процесс дифференцирования, на

n 1

требуем, чтобы было

 

 

 

 

 

 

 

 

 

y(n 1)

C (x) y (n 1)

C

2

(x) y (n 1) ... C

n

(x) y (n 1) .

 

1

1

 

 

2

 

 

n

Тогда придем к последнему из дополнительных условий

шаге по-

(6.8)

C (x) y(n 2) C (x) y

(n 2)

... C

(x) y(n 2)

0 .

 

1

1

 

2

 

2

 

 

n

 

n

 

Дифференцируя (6.8), имеем

 

 

 

 

 

 

 

y(n) C (x) y(n) C

2

(x) y(n) ... C (x) y(n)

 

 

 

1

1

 

 

2

 

n

n

 

 

C (x) y(n 1) C

(x) y(n 1)

... C (x) y(n 1) .

 

 

1

1

2

 

 

2

 

 

n

n

 

Итак, если выполняется система равенств

 

 

 

 

 

 

(x) y2

 

 

 

 

(x) yn

0,

 

 

C1

(x) y1 C2

... Cn

 

 

 

C

(x) y

C

(x) y

... C

(x) y

0,

 

 

 

1

1

2

2

 

 

 

n

n

 

 

 

 

C

(x) y C

(x) y

... C

(x) y

0,

 

 

 

1

1

2

2

 

 

 

n

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C (x) y(n 2) C (x) y(n 2) ... C (x) y(n 2) 0,

1 1 2 2 n n

то

 

y C1 (x) y1 C2 (x) y2 ... Cn (x) yn ,

 

 

 

 

 

y C (x) y

C

2

(x) y

... C

n

(x) y

,

 

 

 

 

 

1

1

 

 

 

2

 

 

 

 

n

 

 

 

 

 

 

y C (x) y C

2

(x) y

... C

n

(x) y ,

 

 

 

 

 

1

1

 

 

 

 

2

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(n 1) C (x) y (n 1) C

2

(x) y (n 1) ... C

n

(x) y (n 1)

,

 

 

1

1

 

 

 

 

 

 

2

 

 

 

 

 

 

n

 

 

y(n) C (x) y(n)

 

C

2

(x) y(n)

... C

n

(x) y(n)

 

 

1

 

1

 

 

 

 

 

2

 

 

 

 

 

 

n

 

 

C

(x) y(n 1) C (x) y(n 1)

... C

(x) y(n 1) .

 

 

1

 

1

 

 

 

 

2

 

2

 

 

 

 

 

n

 

 

n

 

(6.9)

(6.10)

72

Умножая равенства (6.10) соответственно на pn (x) , pn 1(x) ,

pn 2 (x) , …, p1(x) , 1

и подставляя их в уравнение

L[ y] f (x) , полу-

чим

 

 

 

 

 

 

 

C1(x)L[ y1] C2 (x)L[ y2 ] ... Cn (x)L[ yn ]

C

(x) y(n 1) C

(x) y(n 1)

... C

(x) y(n 1)

f (x) .

1

1

2

2

n

n

 

Поскольку y1 , y2 , …,

yn есть решения уравнения L[ y] 0 , то

 

L[ y1] 0 ,

L[ y2 ] 0 , …,

L[ yn ] 0 ,

а, значит,

 

 

 

 

 

 

C (x) y(n 1)

C (x) y(n 1) ... C (x) y(n 1)

f (x) .

1

1

2

2

n

n

 

Это и есть то основное условие, по отношению к которому условия (6.9) являются дополнительными.

Таким образом, если для C1(x) , C2 (x) ,

…, Cn (x) выполнены

условия

 

 

 

 

 

 

 

 

 

 

 

 

 

(x) y1

 

 

(x) y2

 

 

(x) yn

0,

 

 

C1

C2

... Cn

 

 

 

C

(x) y

C

(x) y

... C

(x) y

0,

 

 

 

1

1

 

2

2

 

n

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.11)

 

C

(x) y(n 2)

C (x) y(n 2)

... C (x) y(n 2)

0,

 

1

1

 

 

2

2

 

 

n

n

 

 

C

(x) y(n 1)

C (x) y(n 1)

... C (x) y(n 1) f (x),

 

1

1

 

 

2

2

 

 

n

n

 

то функция (6.6) является решением уравнения L[ y] f (x) . Определителем неоднородной системы уравнений (6.11) является

вронскиан W (x) , а так как функции y1 , y2 , …, yn линейно независи-

мы, то W (x) 0 . Поэтому система (6.11) имеет единственное решение:

C

(x) (x) ,

C

(x)

2

(x) , …,

C

(x)

n

(x) ,

1

 

1

2

 

 

 

n

 

 

где 1(x) ,

2 (x) , …,

n (x)

– некоторые функции, являющиеся ку-

сочно-непрерывными, а, значит, и интегрируемыми. Поэтому получим

73

C1(x) 1(x) C1 , C2 (x) 2 (x) C2 , …, Cn (x) n (x) Cn ,

где C1 , C2 , …, Cn – независимые произвольные постоянные. Подставляя эти функции в равенство (6.6), находим

y [1(x) C1]y1 [2 (x) C2 ]y2 ... [n (x) Cn ]yn . (6.12)

Покажем, что это решение уравнения L[ y] f (x) является его общим решением. Для этого перепишем равенство (6.12) в виде

y C1 y1 C2 y2 ... Cn yn 1(x) y1 2 (x) y2 ... n (x) yn . (6.13)

Первая группа слагаемых справа есть общее решение однородного уравнения L[ y] 0 . Вторая же группа получается из функции (6.12),

если положить в ней С1 С2 ... Сn 0 . Но функция (6.12) удовле-

творяет уравнению L[ y] f (x) при любых C1 , C2 , …, Cn , поэтому вторая группа слагаемых соответствует частному решению v этого уравнения. Следовательно, равенство (6.13) можно переписать в виде

y u v .

Отсюда и следует, что функция (6.13), или, что то же самое, функция (6.12) есть общее решение неоднородного линейного дифференциального уравнения (6.1).

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

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

П р и м е р 6 . 1 . Найти общее решение уравнения y y tg x .

Р е ш е н и е . Характеристическое уравнение 2 1 0 имеет корни 1,2 i . Поэтому фундаментальная система решений соответ-

ствующего однородного уравнения состоит из функций:

74

y1 cos x , y2 sin x .

Запишем общее решение однородного уравнения u C1 cos x C2 sin x .

Тогда общее решение исходного неоднородного уравнения ищем как

y C1(x) cos x C2 (x)sin x .

При n 2 система (6.11) примет вид:

 

 

 

 

(x) y2 0,

C1(x) y1

C2

 

C (x) y

C

(x) y

f (x).

 

1

1

2

2

 

Для данного уравнения получим:

 

C

(x) cos x C

(x) sin x 0,

 

1

2

 

 

C

(x) sin x C (x) cos x tg x.

 

 

1

2

 

Эта система имеет решение:

 

 

C (x)

sin2 x

,

 

C

(x) sin x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

cos x

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируя эти равенства, получим:

 

 

 

 

 

 

 

 

 

 

 

 

cos 2 x 1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1(x)

 

 

dx

cos x

 

 

 

dx sin x ln

tg

 

 

 

 

C1

,

cos x

 

 

 

4

 

 

 

 

 

 

 

 

 

cos x

 

 

 

2

 

 

 

 

 

 

 

 

C2 (x) cos x C2 .

 

 

 

 

 

 

Запишем общее решение

 

 

 

 

 

 

 

 

 

 

 

 

y sin x cos x ln

 

 

 

 

 

x

 

 

 

C1 cos x cos x sin x C2 sin x .

 

 

 

 

 

 

 

 

tg

4

 

 

 

cos x

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно получим

75

y C1 cos x C2 sin x ln

 

 

x

 

 

tg

 

 

 

cos x .

4

2

 

 

 

 

 

Практическое занятие

Т е м а з а н я т и я .

Метод вариации произвольных постоянных. Ц е л ь з а н я т и я .

Изучить метод вариации произвольных постоянных для решения неоднородных линейных дифференциальных уравнений.

А у д и т о р н о е з а н я т и е .

1. Решить дифференциальные уравнения:

 

y 3y

2 y

 

1

 

 

2) y 2y y

3e x

 

 

1)

 

 

;

x 1 ;

 

 

 

 

ex

1

 

 

 

 

 

 

 

 

 

 

 

 

 

3)

y y

2

 

;

 

 

 

 

4) x3 ( y y) x2 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos3 x

 

 

 

 

 

 

 

 

 

 

С а м о с т о я т е л ь н а я р а б о т а .

 

 

 

 

 

 

1. Решить дифференциальные уравнения:

 

 

 

 

 

1)

y 2 y y

ex

;

 

 

 

2) y y

1

 

;

 

 

 

 

 

 

 

x

 

 

 

 

 

sin x

 

 

3)

y 4y 2tg x .

 

 

 

 

 

 

 

 

 

 

Контрольные вопросы

1.Какое уравнение называется неоднородным линейным дифференциальным уравнением n -го порядка?

2.Привести общий вид неоднородного линейного дифференциального уравнения n -го порядка.

3.Сформулировать общие теоремы о неоднородных линейных дифференциальных уравнениях.

4.Доказать общую теорему 1 о неоднородных линейных диффе-

76

ренциальных уравнениях.

5.Доказать общую теорему 2 о неоднородных линейных дифференциальных уравнениях.

6.Сформулировать и доказать теорему об общем решении неоднородного линейного дифференциального уравнения n -го порядка.

7.Какой вы знаете общий метод интегрирования неоднородного линейного дифференциального уравнения?

8.Обосновать метод вариации произвольных постоянных для интегрирования неоднородного линейного дифференциального уравнения.

9.Вывести метод вариации произвольных постоянных для интегрирования неоднородного линейного дифференциального уравнения.

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

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

12.Записать вид общего решения неоднородного линейного дифференциального уравнения n -го порядка, которое находится методом вариации произвольных постоянных.

77

7.НЕОДНОРОДНЫЕ ЛИНЕЙНЫЕ УРАВНЕНИЯ

СПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ

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

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

7.1. Уравнения с правой частью специального вида

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

y(n) p y(n 1)

p

2

y(n 2)

... p

n 1

y p

n

y f (x) ,

(7.1)

1

 

 

 

 

 

 

где p1 , p2 , …, pn 1 , pn

– вещественные числа.

 

 

 

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

78

виду уравнения (7.1).

Поскольку общее решение соответствующего однородного уравнения (5.1) мы уже умеем находить на основании теории и методов раздела 5, то получение общего решения уравнения (7.1) в силу теоремы подраздела 6.2 сводится к нахождению какого-нибудь частного его решения. Оказывается, что эта задача решается весьма просто при некоторых специальных видах функции f (x) в правой части уравнения (7.1). Далее рассматриваются следующие случаи.

1. Уравнения с правой частью в виде полинома.

В этом случае правая часть уравнения (7.1) представима в виде

f (x) Pm (x) ,

где Pm (x) – многочлен степени m с известными коэффициентами. 2. Правая часть в виде произведения полинома и экспоненты. В данном случае правая часть уравнения (7.1) имеет вид

f (x) Pm (x)eax ,

где Pm (x) – многочлен степени m с известными коэффициентами, a

–некоторое число.

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

В этом случае правая часть уравнения (7.1) имеет один из видов:

f (x) P (x)eax sinbx ,

f (x) P (x)eax cos bx ,

 

m

m

 

f (x) [P(x)sinbx R(x) cos bx]eax ,

где Pm (x) –

многочлен степени

m с известными коэффициентами,

P(x) и Q(x)

– многочлены известных степеней с известными коэф-

фициентами,

a и b – некоторые числа.

Для всех этих случаев на основании значений корней характеристического многочлена уравнения (7.1) и вида его правой части можно

79

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

методом неопределенных коэффициентов.

7.2. Уравнения с правой частью в виде полинома

Пусть правая часть неоднородного линейного дифференциального уравнения с постоянными коэффициентами (7.1) имеет вид

f (x) Pm (x) ,

где Pm (x) – многочлен степени m , то есть

P (x) b xm b xm 1 ... b

x b .

 

m

0

 

1

 

 

m 1

 

 

m

 

Тогда уравнение (7.1) примет вид

 

 

 

 

 

 

 

y(n) p y(n 1)

p

2

y(n 2)

... p

n 1

y p

n

y P (x) .

(7.2)

1

 

 

 

 

 

m

 

Рассмотрим два случая.

В п е р в о м с л у ч а е предположим, что в уравнении (7.1) pn 0 . Тогда частное решение уравнения (7.2) можно искать в виде

v Qm (x) ,

(7.3)

где Qm (x) – некоторый многочлен степени m . Действительно, в этом случае, положим

v A0 xm A1xm 1 A2 xm 2 ... Am 1x Am ,

где A0 , A1 , …, Am – неизвестные коэффициенты. Дифференцируя это равенство, получим:

v mA0 xm 1 (m 1)A1xm 2 (m 2)A2 xm 3 ... Am 1 ,

80