Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
лекции матан.pdf
Скачиваний:
85
Добавлен:
31.05.2015
Размер:
1.16 Mб
Скачать

Получили систему двух линейных уравнений с двумя неизвестными a и b. Очевидно, что система имеет определенное решение, и что при найденных значениях a и b функция S (a, b) имеет минимум.

Пример 3.7. Пусть на основании эксперимента получены четыре значения искомой функции у= (х) при четырех значениях аргумента (n= 4), которые записаны в таблице:

Таблица 2

 

х

1

2

3

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

3

4

2,5

 

0,5

 

 

 

 

 

 

 

 

 

 

 

Решение. Будем искать функцию в виде линейной функции y=ax+b.

Составляем выражение S(a, b) :

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

S(a,b) [ yi (axi b)]2 .

 

 

 

 

 

 

 

 

i 1

 

 

Чтобы составить систему (3.16) для определения коэффициентов a и b

предварительно вычисляем

 

 

 

4

21,

4

 

4

 

4

 

10.

y x

x 2

39, x 11,

y

i

 

i i

 

 

i

 

 

i

 

 

i 1

 

i 1

 

 

i 1

i 1

 

 

Система (3.16) принимает вид 21 39a 11b 0,

10 11a 4b 0.

Решая эту систему, находим a и b: a =- 26/35, b = 159/35. Искомая прямая (рис. 16) есть y 2635 x 15935 .

у

 

 

 

 

 

(х,a,b)

у1

х

Рис. 16

4. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ

4.1. Основные понятия и определения

Определение. Дифференциальным уравнением называется уравнение,

связывающее независимую переменную х , искомую функцию у= f(x) и ее производные y , y ,..., y(n) .

Символически дифференциальное уравнение можно записать в виде:

30

F(x, y, y ) 0.

F(x, y, y , y ,..., y (n) ) 0.

Если искомая функция у= f(x) есть функция одной переменной, то уравнение называется обыкновенным.

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

Так, например, уравнение y 2xy 2 7 0 есть уравнение первого порядка.

Определение. Решением или интегралом дифференциального уравнения называется всякая функция у= f(x), которая, будучи подставлена в уравнение, превращает его в тождество.

Дифференциальное уравнение первого порядка имеет вид

(4.1)

Если это уравнение можно разрешить относительно y , то его можно за-

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

Теорема. Если в уравнении

у = f(x, у)

(4.2)

функция f(х, у) и ее частная производная

f

по у непрерывны в некоторой об-

y

 

 

ласти D на плоскости Оху, содержащей некоторую точку (х0; у0), то существует единственное решение этого уравнения у= (х), удовлетворяющее ус-

ловию у= у0 при х= х0.

Геометрический смысл теоремы заключается в том, что существует и притом единственная функция у= (x) , график которой проходит через точ-

ку (х0; у0).

Условие, что при х= х0 функция у должна равняться заданному числу у0, называется начальным условием. Оно часто записывается в виде y x x0 y0 .

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

у= (х, С),

(4.3)

которая зависит от одной произвольной постоянной

С и удовлетворяет сле-

дующим условиям:

 

а) она удовлетворяет дифференциальному уравнению при любом кон-

кретном значении постоянной С;

 

б) каково бы ни было начальное условие у= у0

при х= х0, можно найти

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

31

В процессе отыскания общего решения часто получается соотношение

вида

Ф (х, у, С) = 0

(4.4)

не разрешенное относительно у. Равенство (4.4) называется общим интегралом дифференциального уравнения.

Определение. Частным решением называется любая функция у= (х, С0), которая получается из общего решения у= (х, С), если произвольной постоянной С придать определенное значение. Соотношение

Ф(х, у, С0)=0 называется в этом случае частным интегралом.

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

4.2.Уравнения с разделяющимися переменными

Рассмотрим дифференциальное уравнение вида

dydx f1 (x) f 2 ( y),

где правая часть есть произведение функции, зависящей только функцию, зависящую только от у. Преобразуем его следующим (предполагая, что f 2 ( y) 0) ,

1

dy f1

(x)dx. .

 

f 2 ( y)

 

 

(4.5)

от х, на образом,

(4.6)

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

f 21( y) dy f1 (x)dx C .

Дифференциальное уравнение (4.6) называется уравнением с разделенными переменными, а приводящееся к нему уравнение (4.5) – уравнением с

разделяющимися переменными.

Пример 4.1. Дано уравнение

 

dy

 

y

. Разделим переменные

dy

 

 

dx

.

 

 

 

 

 

 

 

dx

 

x

 

 

 

 

 

 

 

 

 

 

 

y

 

 

x

Интегрируя, находим

 

dy

dx

C, т.е.

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

x

ln

y

ln

x

ln

C

ln

y

ln

C / x

;

 

 

 

 

 

 

 

 

 

 

 

 

отсюда получаем общее решение: у = С/x.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4.2. Дано уравнение

 

(1 x) ydx (1 y)xdy 0.

Разделяя

 

пере-

менные, находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

32

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x

 

 

 

 

 

1 y

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

dx

 

dy 0,

 

1 dx

 

 

 

1 dy 0,

 

 

 

 

 

 

x

 

 

 

 

 

y

 

 

 

 

x

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

интегрируя, получаем ln

 

x

 

x ln

 

y

 

y C

или

ln

 

xy

 

x y C ; последнее со-

 

 

 

 

 

 

 

 

 

 

отношение есть общий интеграл уравнения.

 

 

 

 

 

 

 

4.3. Линейные уравнения первого порядка

Определение. Линейным уравнением первого порядка называется урав-

нение линейное относительно неизвестной функции и ее производной. Оно имеет вид

 

dy

P(x) y Q(x) ,

(4.7)

 

 

 

dx

 

где P(x) и Q(x) - заданные непрерывные функции или постоянные.

 

Решение линейного уравнения (4.7) будем искать в виде произведения

двух функций от х:

 

у= u (x) v (x) .

(4.8)

Одну из этих функций можно взять произвольной, другая определится на основании уравнения (4.7). Дифференцируя обе части равенства (4.8), находим

 

 

 

 

 

dy

u

dv

v

du

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

dx

dx

dx

 

 

 

 

Подставляя полученное выражение производной в уравнение (4.7), бу-

дем иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dv

v

du

Puv Q

 

 

dv

 

du

 

 

u

 

 

или u

 

Pv v

 

Q .

(4.9)

dx

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

dx

 

 

Выберем функцию v такой, чтобы выражение в скобках обратилось в

нуль

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dv

Pv 0 .

 

 

 

(4.10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

Разделяя переменные в полученном уравнении относительно v, находим

dvv Pdx .

Интегрируя, получим

ln C1

ln v

 

Pdx , или

v C1e Pdx .

 

Так как достаточно какого-нибудь отличного от нуля решения уравнения

(4.10), то за функцию v (x)можно взять v e Pdx , где

Pdx - какая-нибудь

первообразная. Очевидно, что v (x) 0.

 

 

 

Подставляя

найденное значение v (x)

в уравнение (4.9), получим

 

du

du

 

Q( x)

, откуда u

Q( x)

C .

 

v(x)

 

Q( x) или

 

 

 

 

dx

 

dx

dx

v( x)

v( x)

 

Подставляя u и v в формулу (4.8), окончательно получим

33

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q( x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y v( x)

 

 

 

 

 

 

 

 

 

dx C .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v( x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 4.3. Решить уравнение

 

dy

 

 

 

2

 

y (x 1)3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

x

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

 

 

 

dv

 

 

 

du

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Пусть y=uv, тогда

 

 

 

 

u

 

 

v

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

dx

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя

выражение

 

 

 

dy

в

 

 

 

исходное

 

уравнение,

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

dv

v

du

 

2

 

uv (x 1)3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dv

 

 

2

 

 

 

v

du

(x 1)

3

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

(4.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для определения v

решим уравнение

dv

 

 

2

 

v 0 , т.е.

 

dv

 

 

 

2dx

, откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

x 1

 

 

 

v

x 1

 

 

 

2 ln

 

x 1

 

или

v (x 1)2 .

Подставляя выражение функции v в уравнение

ln

v

 

 

(4.11),

 

 

 

получаем

для

определения

u

 

уравнение

x 1 2

du

(x 1)3 , или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

du

 

(x

1) , откуда

u

(x 1) 2

C . Следовательно, общий интеграл заданного

 

dx

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнения будет иметь вид

y

(x 1)

4

C(x 1) 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученное семейство является общим решением. Каково бы ни было начальное условие (х0, у0), где х0 -1, всегда можно так подобрать С, чтобы соответствующее частное решение удовлетворяло заданному начальному условию. Например, частное решение, удовлетворяющее условию у0=3 при

х0=0, найдем следующим образом:

3

(0 1)4

C(0 1)2 , C 5 / 2. Следовательно,

2

 

 

 

 

 

 

 

 

 

искомое частное решение таково: y

( x 1)

4

 

5

(x 1)

2 .

Однако, если начальное

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

условие (х0, у0) выбрать так, что х0 = -1, то мы не найдем частного решения, удовлетворяющего этому условию. Это объясняется тем, что при х0 = -1

2

 

 

 

 

 

 

 

 

функция P( x)

 

разрывна и, следовательно, не удовлетворяет условиям

x 1

теоремы существования и единственности.

 

 

 

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

дифференциального уравнения

y ytgx cos2 x , удовлетворяющее начальному условию у(0) =1.

Решение. Положим y=uv, тогда

dy

u

dv

v

du

.

 

 

 

 

 

 

dx

 

dx

 

dx

34

dv

 

 

du

 

2

 

u

 

vtgx

v

 

cos

 

x . Определим v так, чтобы выражение в скобках обра-

 

 

 

dx

 

 

dx

 

 

 

тилось в нуль.

Тогда

dv

vtgx , разделяя переменные, получим

dv

 

sin x

dx ,

 

v

 

 

 

 

 

 

 

 

dx

 

 

cos x

интегрируя уравнение, найдем ln v ln cos x или v cos x

 

 

 

 

 

Для определения u имеем уравнение

 

 

 

 

 

du

 

du

x ; u cos xdx sin x C .

 

 

 

 

cos x

 

cos2 x ,

 

 

cos

 

 

 

 

dx

 

dx

 

 

 

 

Умножив u на v, получим общее решение y cos x(sin x c) . Используя начальное условие у(0) =1, найдем 1= сos0

(sin0 +C), откуда С=1. Искомое частное решение будет иметь вид y cos x(sin x 1) .

4.4. Дифференциальные уравнения высших порядков

Как уже было сказано выше, дифференциальное уравнение n-го порядка можно записать в виде

F(x, y, y , y ,..., y (n) ) 0

или, если его можно разрешить относительно n-й производной,

 

 

 

 

 

 

 

 

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

(4.12)

Для уравнений, разрешенных относительно производной, имеет место

теорема о существовании и единственности решения.

 

Теорема. Если в уравнении

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

функция

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

и

ее частные производные по аргументам

y, y , y ,..., y (n 1)

непрерывны

в

некоторой

области, содержащей

значе-

ния x x0, y y0, y

 

 

 

 

(n 1)

y0

(n 1)

, то существует и притом единственное

 

y0,....,y

 

 

решение y y(x)

уравнения, удовлетворяющее условиям

 

 

y0 ,

 

 

y

 

y0 ,

... ,

y (n 1)

 

y0

(n 1) .

 

y

x x0

 

 

x x0

x x0

 

Эти условия называются начальными условиями.

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

y (x, C1 , C2 ,...,Cn ) ,

зависящая от nпроизвольных постоянных C1 , C2 ,...,Cn и такая, что :

а) она удовлетворяет уравнению

при любых значениях постоян-

ных C1 , C2 ,...,Cn ;

 

 

 

 

 

 

 

 

б) при заданных начальных условиях

 

 

 

 

y

 

x x0 y0 , ... ,

y (n 1)

 

 

y0

(n 1)

y

x x0 y0 ,

 

 

x x0

 

 

 

 

 

35

 

 

постоянные C1 , C2 ,...,Cn можно подобрать так, что функция y (x, C1 , C2 ,...,Cn ) будет удовлетворять этим условиям.

Соотношение вида Ф(x, C1 , C2 ,...,Cn ) 0 называется общим интегралом дифференциального уравнения.

Всякая функция, получающаяся из общего решения при конкретных значениях постоянных C1 , C2 ,...,Cn , называется частным решением.

4.5. Линейные однородные уравнения. Определения и общие свойства

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

производных

 

 

 

 

(n)

, т.е. имеет вид

 

 

 

y , y ,..., y

 

 

 

 

 

 

 

 

a

0

y (n) a y (n 1) ... a

n

y f (x) ,

(4.13)

 

 

 

 

 

1

 

 

где a0 , a1 ,...,an

и

f ( x)

- заданные функции от х или постоянные, причем

a0 0 для всех значений х из той области, в которой рассматривается уравнение (4.13). Будем предполагать, что функции a0 , a1 ,...,an и f (x) непрерывны при всех значениях х, причем коэффициент a0 1 (если он не равен 1, все члены уравнения надо поделить на него). Функция f ( x) , стоящая в правой части уравнения, называется правой частью уравнения.

Если f (x) 0, то уравнение называется линейным неоднородным или

уравнением с правой частью. Если f ( x) 0, то уравнение имеет вид

 

y (n) a

y (n 1) ... a

n

y 0

(4.14)

1

 

 

 

и называется линейным однородным или уравнением без правой части. Установим некоторые основные свойства линейных однородных урав-

нений.

Теорема 1. Если у1 и у2 два частных решения линейного однородного уравнения второго порядка

y a1 y a2 y 0 , (4.15)

то у1 + у2 есть также решение этого уравнения.

Теорема 2. Если у1 есть решение линейного однородного уравнения второго порядка (4.15) и С – постоянная, то Су1 есть также решение этого уравнения.

Определение. Два решения уравнения (4.15) у1 и у2 называются линейно независимыми на отрезке а,b , если их отношение на этом отрезке не яв-

ляется постоянным, т.е. если

y1

const .

y2

 

 

В противном случае решения называются линейно зависимыми. Иными словами, два решения у1 и у2 называются линейно зависимыми на отрезке

36

а,b , если существует такое постоянное число , что

y1

. В этом случае у1

y2

 

 

= у2.

Теорема 3. Если у1 и у2 два линейно независимых решения уравнения

(4.15), то

у=С1 у1 + С2 у2,

где С1 и С2 – произвольные постоянные, есть общее решение этого уравнения.

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

Имеем линейное однородное уравнение второго порядка

y py

q y

0 ,

(4.16)

где p и q – постоянные действительные числа. Чтобы найти общий интеграл этого уравнения достаточно, как следует из теоремы 3, найти два линейно независимых частных решения.

Будем искать частные решения в виде

y ekx , где k = const;

Тогда y kekx ; y k 2 ekx .

Подставляя полученные выражения производных в уравнение (4.16), находим

ekx (k 2 pk q) 0.

 

Так как ekx 0, то значит

 

k 2 pk q 0.

(4.17)

Следовательно, если k будет удовлетворять уравнению (4.17), то e kx будет решением уравнения (4.16). Уравнение (4.17) называется характеристическим уравнением по отношению к уравнению (4.16). Характеристическое уравнение есть квадратное уравнение, имеющее два корня; обозначим их через k1 и k2. При этом

k1

 

p

 

p 2

q,

k2

 

p

 

p 2

q .

 

4

 

4

 

2

 

 

 

2

 

 

Возможны следующие случаи:

1.k1 и k2 – действительные и притом не равные между собой числа;

2.k1 и k2 – действительные равные числа;

3. k1 и k2 – комплексные числа ( k1 i ,

k2 i , где

и

действительные числа, а i – называют мнимой единицей: i2= 1.).

 

Рассмотрим каждый случай отдельно.

 

 

1. Корни характеристического уравнения действительны и различны:

k1 k2. В этом случае частными решениями будут функции

 

y1 ek1x ,

y2 ek2 x .

 

 

 

37

 

 

Эти решения линейно независимы, так как

y1

 

ek1x

e

(k1

k2 ) x

const .

y2

ek2 x

 

 

Следовательно, общее решение имеет вид

 

 

 

 

 

 

 

y C1ek1x C2 ek2 x .

 

 

 

 

 

 

(4.18)

Пример 4.5. Найти общее решение уравнения y y 2 y 0 .

Решение. Составим характеристическое уравнение k 2 k 2 0 . Находим корни характеристического уравнения:

k1,2

 

1

 

1

2; k1

1;

k2 2.

2

4

 

 

 

 

 

 

Общее решение имеет вид y C1e x C2 e 2x .

2. Корни характеристического уравнения действительные и равные. В

этом случае k1= k2.

Одно частное решение y1 ek1x получается на основании предыдущих рассуждений. Нужно найти второе частное решение. Линейно независимое с

первым (функция ek2 x тождественно равна e k1x и поэтому не может рассматриваться в качестве второго частного решения).

Будем искать второе частное решение в виде y2 u(x)ek1x , где u(x) - неизвестная функция, подлежащая определению.

Дифференцируя, находим

 

 

k1x

k1ue

k1x

e

k1x

(u

 

uk1 ) ,

 

 

y2

u e

 

 

 

 

 

y

u ek1x

2k u ek1x k 2

uek1x ek1x (u

2k u uk 2 ).

2

 

 

1

 

 

1

 

 

 

1

1

Подставляя

выражения

производных

в

уравнение (4.16), получаем

ek1x [u (2k1 p)u (k12 pk1 q)u] 0 .

Так как k1 – кратный корень характеристического уравнения, то

k1= k2 = - p/2 или 2k1 p,

2k1 p 0 .

Следовательно, для

того чтобы найти u(x) , надо решить уравнение

ek1x u 0, или u 0, интегрируя уравнение, получаем u=Ax+B. В частности, можно положить А=1, В=0. Таким образом, в качестве второго частного решения можно взять

y2 xek1x .

 

 

 

 

 

 

Это решение линейно независимо с первым,

так как

y2

x const . Поэтому

y1

 

 

 

 

 

 

 

 

общим решением будет функция

 

 

 

 

 

 

 

 

y C ek1x C

2

xek1x ek1x (C

xC

2

) .

(4.19)

1

 

1

 

 

 

 

38

Пример 4.6. Найти общее решение уравнения y 4 y 4 y 0 .

Решение. Составим характеристическое уравнение k 2 4k 4 0 . Находим корни характеристического уравнения: k1 k2 2.

Общее решение имеет вид y C1e2x C2 xe2x .

3. Корни характеристического уравнения комплексные.

Так как комплексные корни входят попарно сопряженными, то обозна-

чим k1 i ,

k2 i , где p / 2,

 

q

p 2

.

4

 

 

 

 

 

 

 

 

 

Частные решения можно записать в виде

 

 

 

 

y

1

e( i ) x ,

y

2

e( i ) x .

(4.20)

 

 

 

 

 

 

 

 

Это комплексные функции действительного аргумента, удовлетворяющие дифференциальному уравнению (4.16).

Очевидно, что если какая-либо комплексная функция действительного аргумента y u(x) iv(x) удовлетворяет уравнению (4.16), то этому уравнению удовлетворяют функции u(x) и v( x) . Эти функции являются решениями уравнения (4.16).

Перепишем комплексные решения (4.20) в виде суммы действительной

и мнимой части:

y1 e x cos x ie x sin x,

y2 e x cos x ie x sin x.

Тогда частными решениями уравнения (4.16) будут действительные

 

 

 

e x

cos x,

 

 

 

 

 

 

 

e x sin x.

 

 

функции y

 

y

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функции

y1

,

 

 

y2

линейно независимы, так как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e x cos x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

1

 

 

ctg x

const.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e x sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

корней характеристического уравнения имеет вид

 

y C

 

 

C

 

 

 

 

C e x cos x C

 

e x

sin x

 

 

y

1

2

y

2

2

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y e x (C1 cos x C2 sin x) ,

(4.21)

где С1 и С2 - произвольные постоянные.

Важным частным случаем решения (4.21) является случай, когда корни характеристического уравнения мнимые. Это имеет место, когда в уравнении

(4.16) р=0, и оно имеет вид y qy 0 .

 

 

 

 

 

 

 

Корни характеристического уравнения k1,2 i

 

q i . Решение (4.21)

принимает вид y C1 cos x C2 sin x .

 

 

 

Пример 4.7. Найти общее решение уравнения

y 2 y 5y 0 и частное

решение, удовлетворяющее начальным условиям y x 0 0, y x 0 1.

39

Решение. Составим характеристическое уравнение

k 2 2k 5 0 и

найдем его корни:

 

 

 

 

 

 

 

 

 

k1 1 2i, k2 1 2i .

 

 

Общее решение имеет вид

y e x (C1 cos2x C2 sin 2x) .

Найдем частное решение, удовлетворяющее заданным начальным усло-

виям, определим

С1 и С2.

На основании первого условия находим:

0= e 0 (C cos(2 0) C

2

sin(2 0)),

откуда

С

.

Найдем

производную

1

 

 

1

 

 

 

y e x (2C2 cos2x C2 sin 2x) . Из второго условия получим 1=2 С2, т.е. С2=1/2. Таким образом, искомое частное решение имеет вид

y 12 e x sin 2x .

4.7. Неоднородные линейные уравнения второго порядка

Пусть имеем неоднородное линейное уравнение второго порядка

y a1 y a2 y f (x).

(4.22)

Структура общего решения такого уравнения определяется следующей теоремой:

Теорема. Общее решение неоднородного уравнения представляется как сумма какого-нибудь частного решения этого уравнения у* и общего реше-

ния y соответствующего однородного уравнения y a1 y a2 y 0.

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

Пусть имеем уравнение

y py qy f (x) ,

(4.23)

где p и q - действительные числа.

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

1. Пусть правая часть уравнения (4.23) представляет собой произведе-

ние показательной функции на многочлен, т.е. имеет вид

 

f x Pn (x)e x ,

(4.24)

где Pn ( x) - многочлен n й степени. Тогда возможны случаи.

 

а) Число не является корнем характеристического

уравнения

k 2 pk q 0.

 

В этом случае частное решение нужно искать в виде

 

y Qn (x)e x ( A0 x n A1 x n 1 .. An ) .

(4.25)

40

 

Действительно, подставляя решение (4.25) в уравнение (4.23) и сокращая все члены на множитель e x , будем иметь:

Qn (x) (2 p)Qn (x) ( 2 p q)Qn (x) P(x).

(4.26)

Qn (x) - многочлен степени n, Qn ( x) - многочлен степени n – 1,

Qn (x) - многочлен степени n – 2. Таким образом, слева и справа от знака равенства стоят многочлены n-й степени. Приравнивая коэффициенты при одинаковых степенях х, получим систему n+1 уравнений для определения неизвестных коэффициентов A0 , A1 , A2 ,..., An .

б) Число есть простой корень характеристического уравнения.

Если бы в этом случае частное решение мы стали искать в форме (4.25), то в равенстве (4.26) слева получился бы многочлен (n-1)степени, так как ко-

эффициент при Qn (x) , т.е.

( 2 p q) равен нулю, а многочлены

Qn ( x)

и

Q

 

(x) имеют степень, меньшую n. Следовательно, ни при каких A , A , A ,..., A

 

n

 

0

1 2

n

равенство (4.26) не было бы тождеством. Поэтому в рассматриваемом случае частное решение нужно искать в виде многочлена (n +1)-й степени, но без свободного члена (так как свободный член этого многочлена исчезнет при дифференцировании):

y xQn (x)e x .

(4.27)

в) Число есть двукратный корень характеристического уравнения. То-

гда в результате подстановки в дифференциальное уравнение

функ-

ции Qn (x)e x степень многочлена понизится на две единицы.

 

Следовательно, в левой части равенства (4.26) остается многочлен (n-2)- й степени. Для того чтобы в результате подстановки получить многочлен

степени n, следует частное решение искать в виде произведения e x на многочлен (n+2)-й степени. При этом свободный член этого многочлена и член в первой степени исчезнут при дифференцировании; поэтому их можно не включать в решение.

Итак, в случае, когда есть двукратный корень характеристического

уравнения, частное решение можно брать в форме

 

y x 2Qn (x)e x .

(4.28)

Пример 4.8. Найти общее решение уравнения y 4 y 3y x. Решение. Общее решение соответствующего однородного уравнения

имеет вид

 

C e x C

 

xe 3x .

 

y

2

 

 

 

 

1

 

 

Так как

правая

часть данного неоднородного уравнения имеет вид

xe0x (т.е. вид

P (x)e0x ), причем 0 не является корнем характеристического

 

 

 

1

 

 

 

уравнения, то частное решение будем искать в форме

y Q (x)e0x , т.е. поло-

 

 

 

 

 

 

1

жим y A0 x A1 . Подставляя это выражение в заданное уравнение, будем иметь 4A0 3( A0 x A1 ) x .

41

 

Приравнивая

коэффициенты

при одинаковых степенях х, получим

3A0 1,

4A0 3A1

0 ,

откуда,

A0 1/ 3,

A1 4 / 9.

Следовательно,

y

1

 

4

. Общее решение y

 

y будет

 

 

 

 

x

y

 

 

 

 

 

 

 

 

 

 

 

3

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y C1e x C2 e 3x

1

x

4

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

9

 

 

Пример 4.9. Найти общее решение уравнения. y 9 y (x 2 1)e3x .

Решение. Общее решение будет иметь вид y y y .

Найдем общее решение соответствующего однородного уравнения. Со-

ставим

характеристическое

уравнение

и

найдем

его

корни

k 2 9 0,

k

1,2

3i.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

общее

решение

соответствующего

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

 

 

C1 cos3x C2 sin 3x .

 

 

 

 

 

 

y

 

 

 

 

 

 

 

Правая часть

данного

неоднородного

уравнения (x 2 1)e3x

имеет

вид P2 (x)e3x . Так как коэффициент 3 в показателе степени не является корнем характеристического уравнения, то частное решение будем искать в форме y Q2 (x)e3x , т.е. положим y ( Ax2 Bx C)e3x . Подставляя это выражение в заданное уравнение, будем иметь

[9( Ax2 Bx C) 6(2Ax B) 2A 9( Ax2 Bx C)]e3x (x 2 1)e3x

Сокращая на e 3x и приравнивая коэффициенты при одинаковых степенях х, получим

18A 1,

12A 18B 0 ,

2A 6B 18C 1,

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

A 1/ 18,

B 1/ 27.,

C 5 / 81.

Следовательно,

 

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

 

 

1

 

2

 

1

 

 

5

 

3x

 

 

 

 

 

 

 

 

 

 

 

 

 

будет y

 

 

 

x

 

 

 

 

x

 

 

 

e

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

18

 

 

27

 

81

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Общее решение y

 

y будет

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

1

 

 

5

 

3x

 

 

 

 

 

 

 

y C1 cos3x C2 sin 3x

 

 

x

 

 

 

 

x

 

e

 

.

 

 

 

 

 

 

18

 

27

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

81

 

 

Пример 4.10. Найти общее решение уравнения. y 7 y 6 y (x 2)e x .

Решение. Общее решение будет иметь вид y y y .

Найдем общее решение соответствующего однородного уравнения. Со-

ставим

характеристическое

 

уравнение

и

найдем

его

корни

k 2 7k 6 0,

k

1

6,

k

2

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда общее решение соответствующего однородного уравнения

 

 

 

 

 

 

 

 

 

 

 

C e6x C

 

e x .

 

 

 

 

 

 

 

 

 

 

 

 

 

y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

42

 

 

 

 

 

 

Правая часть данного неоднородного уравнения (x 2)e x имеет вид P1 (x)e1x . Так как коэффициент 1 в показателе степени является простым корнем характеристического уравнения, то частное решение будем искать в форме y x( Ax B)e x . Подставляя это выражение в заданное уравнение, будем иметь

[(Ax 2 Bx) (4Ax 2B) 2A 7(Ax 2 Bx) 7(2Ax B) 6(Ax 2 Bx)]ex (x 2)ex .

( 10Ax 5B 2A)ex (x 2)ex.

Сокращая на e x и приравнивая коэффициенты при одинаковых степенях х, получим 10A 1, 5B 2A 2 , откуда A 1/ 10, B 9 / 25. Следова-

 

 

 

1

тельно, частное решение будет

y x

 

 

 

 

 

 

 

 

10

 

9

 

x

 

 

 

 

x

 

Общее решение y y y

 

e

 

.

 

 

 

25

 

 

 

 

 

 

 

 

 

6x

 

x

 

 

1

 

 

9

 

 

x

 

 

y C1e

C2 e

 

x

 

будет

 

 

+ x

 

 

 

 

 

e

 

.

 

 

10

25

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Пусть правая часть уравнения имеет вид f (x) P(x)e x cos x Q(x)e x sin x ,

где P(x), Q(x) - многочлены от х, то форма частного решения определяется

следующим образом:

а) если + i не является корнем характеристического уравнения, то частное решение уравнения (4.23) следует искать в виде

y U(x)e x cos x V (x)e x sin x ,

где U(x), V (x) - многочлены, степень которых равна наивысшей степени многочленов P(x), Q(x) ;

б) если + i является корнем характеристического уравнения, то частное решение уравнения (4.23) следует искать в виде

y x[U(x)e x cos x V (x)e x sin x] .

Следует отметить, что указанные формы частных решений сохраняются и в том случае, когда в правой части уравнения (4.23) один из многочленов P(x), Q(x) тождественно равен нулю, т.е., когда правая часть имеет вид

P(x)e x cos x или Q(x)e x sin x .

Рассмотрим важный частный случай. Пусть правая часть линейного уравнения второго порядка имеет вид:

f (x) M cos x N sin x ,

где M и N – постоянные числа. Тогда:

а) если i не является корнем характеристического уравнения, то частное решение уравнения (4.23) следует искать в виде y Acos x B sin x ;

б) если i является корнем характеристического уравнения, то частное решение уравнения (4.23) следует искать в виде y x[ Acos x B sin x] .

43

Пример 4.11. Найти общее решение линейного неоднородного уравне-

ния y 2 y 5y 2 cos x.

Решение. Общее решение будет иметь вид y y y .

Найдем общее решение соответствующего однородного уравнения. Со-

ставим

характеристическое

 

уравнение

и

найдем

его

корни

k 2 2k 5 0,

k

1

1 2i, k

2

1 2i.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда общее решение соответствующего однородного уравнения

 

 

 

e

x (C

cos2x C

 

 

 

 

 

 

 

 

 

y

2

sin 2x).

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

Правая часть данного неоднородного уравнения 2 cos x , очевидно, что

i =1i не является

корнем характеристического уравнения, частное решение

будем искать в форме y Acos x B sin x , где А и В – постоянные коэффициенты, подлежащие определению.

Найдем производные y :

y

Asin x B cos x, y Acos x B sin x.

 

Подставляя выражения

y и производных в заданное уравнение,

будем

иметь

 

 

 

 

Acos x B sin x 2( Asin x B cos x) 5(cos x B sin x) 2 cos x .

 

Приравнивая коэффициенты при cos x и sin x , получим два уравнения

для определения А и

В:

 

A 2B 5A 2, B 2A 5B 0 ,

откуда

A 2 / 5, B 1/ 5.

решение y y y

Следовательно, частное решение y

2

cos x

1

sin x .Общее

5

5

 

 

 

будет иметь вид

 

 

 

 

y e x (C1 cos2x C2 sin 2x) + 25 cos x 15 sin x .

Пример 4.12. Найти общее решение линейного неоднородного уравне-

ния y 4 y cos2x.

Решение. Общее решение будет иметь вид y y y .

Найдем общее решение соответствующего однородного уравнения. Со-

ставим

характеристическое

уравнение

и

найдем

его

корни

k 2 4 0,

k

1

2i,

k

2

2i.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда общее решение соответствующего однородного уравнения

y C1 cos2x C2 sin 2x.

Правая часть данного неоднородного уравнения cos2x , очевидно, что

i =2 i является корнем характеристического уравнения, частное решение будем искать в форме

y x( Acos x B sin x),

где А и В – постоянные коэффициенты, подлежащие определению. Найдем производные y :

y 2x( Asin 2x B cos2x) ( Acos2x B sin 2x),

44

y 4x( Acos2x B sin 2x) 4( Asin 2x B cos2x).

Подставляя выражения y и производных в заданное уравнение и при-

равнивая коэффициенты при cos2x

и sin 2x , получим два уравнения для оп-

ределения А и В: 4B 1,

4A 0 , откуда A

0, B 1/ 4. Следовательно, ча-

стное решение y

1

x sin x . Общее решение y

 

 

y будет иметь вид

y

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

y C cos2x C

 

sin 2x +

1

x sin x .

 

 

 

2

 

 

 

 

 

1

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

Пример 4.13. Найти частное решение линейного неоднородного урав-

нения

y y 2y cos x 3sin x , удовлетворяющее начальным условиям:

y(0) 1,

y (0) 2.

 

 

 

 

 

 

 

 

Решение. Найдем вначале общее решение, которое будет иметь вид y y y .

Чтобы найти общее решение соответствующего однородного уравнения

 

 

 

 

 

 

 

 

y , составим характеристическое

уравнение

и найдем

его

корни

k 2 k 2 0,

 

k

1

1,

k

2

2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда общее решение соответствующего однородного уравнения

 

 

 

C e x C

 

e 2x .

 

 

 

 

 

 

 

 

 

y

2

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

Правая

 

часть

данного неоднородного

уравнения

имеет

вид

 

f (x) e x (M cos x N sin x) ,

в

данном

случае

0, 1,

i i, M 1, N 3.Так

как число

+ i = i не является корнем

характеристического уравнения, то частное решение будем искать в форме

y ( Acos x B sin x), y Asin x B cos x,

y Acos x B sin x.

Подставим выражения y и производных в заданное уравнение, получим после приведения подобных членов

(B 3A) cos x ( A 3B) sin x cos x 3sin x.

 

Приравняем коэффициенты при cos x

и sin x ,

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

для определения А и В:

B 3A 1,

3B A 3 , откуда A

0,

B 1. Следо-

вательно,

частное решение y sin x . Общее решение y

 

y

будет иметь

y

вид

y C e x C

2

e 2x sin x .

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем С1, С2, используя начальные условия

 

 

 

 

 

 

 

 

 

C e0

C

2

e0

sin 0 1;

 

 

C C

2

1;

 

 

 

 

 

 

 

1

 

 

 

или

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1e0 2C2 e0

cos0 2,

C1 2C2

1 2.

 

 

Отсюда С1=1, С2=0. Искомое частное решение будет y e x sin x .

45