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

Lesn3_Integraly

.pdf
Скачиваний:
6
Добавлен:
11.05.2015
Размер:
2.02 Mб
Скачать

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

Представим уравнение в виде y(n) f (x, y, y ,..., y(n 1) ) . Получим y(n) q(x) p1(x) y(n 1) ... pn (x) y . Следовательно,

f q(x) p1(x) y(n 1) ... pn (x) y .

Тогда

f

p (x) ,

f

p

(x) , … ,

f

 

(n 1)

p (x) ;

 

y

n

y

 

n 1

 

 

y

 

1

функции f ,

 

 

 

 

определены и непрерывны на

f y , f y , … , f y(n 1)

множестве D = (a; b) R n.

 

 

 

 

 

 

 

 

 

По теореме Коши для любой точки M (x

, y

0

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

 

 

 

 

 

 

0

 

 

0

0

где x0 (a;b) , а остальные координаты произвольны, дифферен-

циальное уравнение имеет единственное решение y(x), удовлетворяющее начальному условию в точке M. Каждое из этих решений является n раз дифференцируемой на (a; b) функцией.

Рассмотрим оператор

A: Ln (a;b) L(a;b) , определяемый

равенством

 

 

A( y) y(n) p (x) y(n 1) ... p (x) y .

(4)

1

n

 

Этот оператор является линейным. Действительно, для

любых функций u(x), v(x) из Ln (a;b) и любых чисел

, R

в силу линейности операции дифференцирования имеем:

A( u v) ( u v)(n) p1(x)( u v)(n 1) ... pn (x)( u v)

( u(n) v(n) ) p1(x)( u(n 1) v(n 1) ) ... pn (x)( u v)

[u(n) p1(x)u(n 1) ... pn (x)u] [v(n) p1(x)v(n 1) ... pn (x)v]

A(u) A(v) .

Используя оператор A, линейное уравнение (3) можно записать кратко A( y) q(x).

сследуем теперь линейное однородное уравнение:

y(n) p (x) y(n 1)

... p

(x) y 0 .

(5)

1

n

 

 

При условии, что функции

p1(x), p2 (x), ... , pn (x), q(x)

непрерывны на интервале (a; b), оно имеет на этом интервале бесконечное множество решений. Рассмотрим сначала вопрос о линейной зависимости этих решений на интервале (a; b).

110

§3. Линейные однородные дифференциальные уравнения

Теорема 2. Пусть y1, y2 ,..., yn - система решений уравнения (5), линейно независимая на интервале (a; b). Тогда

W[ y1,..., yn ](x) 0 x (a;b) .

Доказательство опустим.

Теоремы 1 и 2 позволяют решать вопрос о линейной независимости решений линейного однородного уравнения (5).

Следствие. Пусть x0 (a; b) - произвольная точка. Решения y1, y2 ,..., yn линейного однородного уравнения линейно независимы на интервале (a; b) тогда и только тогда, когда

W[ y1,..., yn ](x0 ) 0 .

Рассмотрим в линейном пространстве Ln (a;b) подмноже-

ство L всех решений линейного однородного уравнения (5). Пусть y1, y2 L - произвольные его решения. Для них на

интервале (a; b) выполняются тождества A( y1) 0 , A( y2 ) 0 . Из этих тождеств для любого числа следует:

 

A( y1 y2 ) A( y1) A( y2 ) 0 ,

A( y1) A( y1) 0 .

Тождества означают, что функции

y1 + y2,

y1

являются

решениями уравнения (5), то есть

y1 + y2 L, y1

L.

 

Таким образом, доказано следующее утверждение.

Теорема 3. Множество L всех решений линейного однородного

 

уравнения на интервале

(a;

b)

относительно сложения

 

 

решений и умножения решения на число образует под-

 

пространство линейного пространства Ln (a;b)

всех n

 

раз дифференцируемых функций.

 

 

Следствие. Если y1(x), y2 (x),...,yn (x) - решения линейного од-

нородного уравнения (5) и с1, с2,…, сn R, то функция n

y(x) ci yi (x) является решением этого уравнения. i 1

Выясним теперь структуру общего решения линейного однородного уравнения.

111

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

3. Структура общего решения линейного однородного дифференциального уравнения

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

Определение 4. Система решений y1, y2 ,..., yk линейного одно-

родного уравнения порядка n на интервале (a; b) назы-

вается фундаментальной системой решений, если эта система линейно независима на данном интервале и k = n.

Теорема 4. (О структуре общего решения).

Пусть y1, y2 ,..., yn - фундаментальная система решений

линейного однородного уравнения (5). Тогда общее решение этого уравнения имеет вид:

n

 

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

(6)

i 1

 

Доказательство сводится к проверке свойств общего решения.

1. Пусть

c* ,c* ,...,c* - произвольные действительные

 

1

2

n

 

 

 

 

 

 

числа. Тогда

согласно

следствию

из

теоремы 3 функция

 

n

 

 

 

 

 

 

 

 

y(x, c1* , c2* ,...,cn* ) ci* yi (x)

является решением уравнения (5).

 

i 1

 

 

 

 

 

 

 

2. Пусть

u(x) - некоторое решение уравнения (1) на интер-

вале (a; b). Найдем такие числа

c* ,c* ,...,c* R, чтобы на (a; b)

 

 

 

 

 

 

 

1

2

n

выполнялось равенство u(x) y(x,c* ,c* ,...,c* ) .

 

 

 

 

 

 

 

1

2

n

Возьмем произвольную точку

x0 (a; b) и вычислим зна-

чения функции u(x)

и ее производных в этой точке.

 

 

u( x0 ) u0 ;

 

 

 

 

u (x

0

) u ;

 

 

 

 

 

 

 

 

0

 

 

 

 

................

 

(*)

 

 

 

 

 

 

 

 

 

 

 

 

 

(n 1)

 

 

(n 1)

.

 

 

u

 

 

(x0 ) u0

112

§3. Линейные однородные дифференциальные уравнения

Составим систему уравнений относительно C1,...,Cn :

y(x0 ,C1,...,Cn ) u0 ;

 

 

 

 

 

 

 

 

 

 

y (x

0

,C ,...,C

n

) u ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

................................

 

 

 

 

 

 

 

 

(**)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(n 1) (x

0

,C ,...,C

n

 

) u(n 1) ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

Перепишем систему следующим образом:

 

 

 

C1 y1 (x0 ) C2 y2 (x0 ) ... Cn yn (x0 ) u0 ;

 

 

C y (x

0

) C

2

y

(x

0

) ... C

n

y

(x

0

) u

;

 

 

1 1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

n

 

 

0

 

 

 

................................................................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C y(n 1)

(x

0

) C

2

y(n 1) (x

0

) ... C

n

y(n 1) (x

0

) u(n 1) ;

1 1

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

n

 

0

Определитель системы равен определителю Вронского

W[ y1,..., yn ](x0 ) . Так как

 

 

y1, y2 ,..., yn

-

 

фундаментальная систе-

ма решений уравнения (5), то по теореме 2

W[ y1,..., yn ](x0 ) 0 .

Тогда система (**)

 

имеет единственное решение

c* ,...,c* . Его

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

n

можно найти, например, по правилу Крамера.

 

 

 

 

 

 

 

Рассмотрим функцию

 

y(x,c*

,...,c* ) . Она тоже

является

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

n

 

 

 

 

 

 

 

 

решением уравнения (5). Из систем равенств (*) и (**) следует, что решения u(x) и y(x,c1* ,...,cn* ) удовлетворяют одному и тому же начальному условию задачи Коши. Тогда по теореме Коши на интервале (a; b) выполняется равенство u(x) y(x,c1* ,...,cn* ) .

Доказанные свойства функции (6) означают, что она является общим решением линейного однородного уравнения (5).

Согласно доказанной теореме нахождение общего решения линейного однородного уравнения сводится к нахождению фундаментальной системы решений этого уравнения. Всегда ли она существует?

Теорема 5. Для любого линейного однородного уравнения высшего порядка существует хотя бы одна фундаментальная система решений.

113

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

Доказательство. Рассмотрим линейное однородное уравнение порядка n, имеющее решения на интервале (a; b). Возьмем произвольную точку x0 (a; b). Запишем для этого уравнения n начальных условий.

 

 

 

0

)

 

 

 

 

 

 

 

 

y( x

 

 

1

 

0

 

 

0

 

 

( x0 )

 

 

 

 

 

 

 

y

 

0

1

 

0

...............

 

...

...

...

....

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n 1)

( x0 )

 

 

 

 

 

 

 

y

 

 

 

0

0

 

1

Согласно теореме Коши для каждого i-го условия существует единственное решение yi(x) уравнения, которое удовлетворяет этому начальному условию. Получаем n решений y1(x), y2(x), …, yn(x) рассматриваемого уравнения.

Рассмотрим определитель Вронского для этих решений y1(x), y2(x), … , yn(x), вычисленный в точке x0. Очевидно, W[ y1,..., yn ](x0 ) | E | 1 0 . Это означает, что система решений

y1(x), y2(x), … , yn(x) линейно независима на интервале (a; b), то есть является фундаментальной системой решений уравнения.

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

Следствие. Линейное пространство L решений линейного однородного уравнения порядка n является n-мерным.

Действительно, согласно теореме 5 в пространстве решений L существует фундаментальная система, состоящая из n решений. По теореме 4 она является базисом пространства L.

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

114

§3. Линейные однородные дифференциальные уравнения

Лекция 13

4.Множество комплексных чисел

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

ножество комплексных чисел является расширением множества действительных чисел. Действительные числа интерпретируются как точки координатной прямой. При этом между числами и точками прямой устанавливается взаимно однозначное соответствие. "Свободных мест" на координатной прямой нет. Поэтому комплексные числа можно попытаться интерпретировать как точки на координатной плоскости. Каждая такая точка характеризуется парой действительных чисел - декартовых координат. Поэтому каждое комплексное число должно полностью определяться упорядоченной парой действительных чисел.

Будем рассматривать комплексные числа как многочлены первой степени a + bi переменной i. Переменная i называется мнимой единицей и обладает особым свойством

i 2 = –1.

Действительные числа a и b называются соответственно,

действительной и мнимой частями комплексного числа.

Обозначения: z = a + bi, a = Re z, b = Im z;

C - множество всех комплексных чисел z.

Комплексные числа складываются и умножаются как многочлены первой степени переменной i. При записи результата в стандартном виде учитывается равенство i 2 = –1.

Отметим особо пары комплексных чисел вида z = a + bi и z a bi . Они называются комплексно сопряженными числами.

Их сумма и произведение являются действительными числами. ассмотрим свойства корней многочленов с действитель-

115

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

ными коэффициентами.

Любое квадратное уравнение az2 bz c 0 с действительными коэффициентами имеет два корня, которые находятся по обычной формуле корней квадратного уравнения

z

b

b2 4ac

 

.

 

 

 

1,2

 

2a

 

 

 

 

 

 

Пример 4.

 

 

 

 

Решим квадратное уравнение

z2 2z 5 0 .

По формуле корней квадратного уравнения получаем:

z1,2

 

2

 

 

4 4 1 5

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

2

 

16

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 4

 

 

 

 

 

2 4 i2

 

 

 

 

 

 

 

 

 

 

 

z

 

1

 

 

2 4i

;

 

 

 

 

 

1,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z1,2

1 2i .

 

 

 

 

 

 

 

 

 

 

 

 

 

Как и в случае действительных корней, имеет место раз-

ложение квадратного трехчлена:

 

 

 

 

 

 

 

 

 

 

 

z 2 2z 5 (z ( 1 2i))( z ( 1 2i)) ( z z

)( z z

2

) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

В общем случае любой многочлен

 

 

 

 

 

 

 

 

 

 

P (z) zn a zn 1 ... a

 

z a

,

 

 

 

 

 

n

 

1

 

 

 

 

n 1

 

 

n

 

 

 

 

n > 1, с действительными коэффициентами

ai

единственным

образом раскладывается на множители вида

(z z0 )s , где все

числа z0 C

различны и s

.

 

Сумма всех показателей

s

равна степени

n

многочлена. Разложение такого вида называют

стандартным разложением многочлена на множители.

 

Если стандартное разложение многочлена

Pn(z) содержит

множитель (z z0 )s , то число

z0 является корнем этого множи-

теля и, вместе с тем, корнем многочлена

 

 

Pn(z). Показатель

s

называется кратностью корня z0. Если

s = 1, то корень назы-

вается простым.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если число

z0 = a + bi -

корень многочлена Pn(z), то и со-

пряженное число z0 a bi является корнем этого многочлена.

116

§3. Линейные однородные дифференциальные уравнения

5.Линейные однородные уравнения

спостоянными коэффициентами

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

Определение 5. Дифференциальное уравнение вида

 

y(n) a y(n 1) ... a y 0 ,

(7)

 

1

n

 

 

где a1, …, an ,

называется линейным однородным

 

уравнением с постоянными коэффициентами.

 

 

спомним: линейное однородное уравнение

y a1 y 0 с

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

имеет общее

решение y Ce a1x . По аналогии с этим будем искать решения

уравнения (7) в виде

y ekx , где

k R.

 

 

 

 

 

Подставив в уравнение (7) функцию

y

и ее производные

y kekx ,

y k 2ekx ,..., y(n)

k nekx , получим равенство

 

 

 

ekx (k n a k n 1 ... a

k a ) 0 .

 

 

 

 

 

1

 

n 1

n

 

 

 

 

 

Так как ekx 0 , то функция

 

y ekx

является решением

уравнения (7) тогда и только тогда, когда

k

является решением

алгебраического уравнения

k n a k n 1 ...

a

k a 0 .

 

 

 

 

 

1

 

 

 

n 1

n

 

Определение 6. Алгебраическое уравнение

 

 

 

 

 

 

 

k n a k n 1

... a

 

k a

0

 

(8)

 

 

 

 

 

 

 

1

 

n 1

n

 

 

 

 

 

называется характеристическим уравнением диффе-

 

ренциального уравнения (7).

 

 

 

 

 

 

 

Алгебраическое уравнение (8) может иметь как действи-

тельные

корни k

кратности

s,

 

так и комплексные

i

кратности t. При этом сумма кратностей всех корней равна n.

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

117

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

уравнения (7).

1. Для каждого действительного корня k кратности s в фундаментальную систему включаем набор из s функций:

ekx , xekx , x2ekx ,..., xs 1ekx .

2. Для каждой пары i комплексных корней кратности t в фундаментальную систему включаем 2 набора по t функций:

 

e x sin x,

xe x sin x, ..., xt 1e x sin x ;

 

e x cos x,

xe x cos x, ..., xt 1e x coss x .

Оформим это соответствие в виде таблицы:

 

 

 

 

 

 

Корни

Кратности

 

 

Функции ФСР

 

 

 

 

 

 

k

s

 

ekx ,

xekx , x2ekx ,...,

xs 1ekx

+ i

 

 

e x sin x,

xe x sin x, ...,

xt 1e x sin x

i

t

 

e x cos x,

xe x cos x, ...,

xt 1e x coss x

 

 

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

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

аким образом, получаем

алгоритм решения линейного однородного уравнения

спостоянными коэффициентами

1.По дифференциальному уравнению (7) составляем характеристическое уравнение (8).

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

3.По таблице, в соответствии с корнями и их кратностями, составляем фундаментальную систему решений.

4.Записываем общее решение (6) этого уравнения.

118

§3. Линейные однородные дифференциальные уравнения

Пример 5.

Решим уравнение y(5) 9y 0 .

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

1.Записываем характеристическое уравнение k5 9k3 0 .

2.Находим его корни и их кратности:

k3 (k 2 9) 0 ;

k 3 (k 3)(k 3) 0 ;

k1 0,

k2 3,

 

 

k3 3 ;

 

 

 

 

 

 

 

 

 

s = 3,

s = 1,

 

 

 

s = 1.

 

 

 

 

 

 

 

 

 

3. По корням составляем фундаментальную систему решений:

k 0,

s 3,

1, x, x2 ;

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k2 3, s 1,

e 3x ;

 

 

 

 

 

 

 

 

 

k3 3,

s 1,

e3x .

 

 

 

 

 

 

 

 

 

4) Записываем общее решение дифференциального уравнения:

y C 1 C x C x2

C e 3x C e3x .

 

 

 

 

 

1

2

 

 

3

4

5

 

 

 

 

 

 

 

Пример 6.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решим уравнение y 8y 0 .

 

 

 

 

 

 

 

 

 

1. Записываем характеристическое уравнение

 

 

 

k3 8 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Находим его корни:

 

 

 

 

 

 

 

 

 

 

 

 

 

k3 23

0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(k 2)(k 2 2k 4) 0 ;

 

 

 

 

 

 

 

 

 

k1 2 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 4 4 1 4

 

 

 

 

 

 

 

 

 

 

k

2,3

 

1 3 1 3i .

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. По корням составляем фундаментальную систему решений:

k 2,

s 1,

e2x ;

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

2,3

1 3i,

 

t 1,

 

e x cos 3x,

e x sin 3x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4. Записываем общее решение дифференциального уравнения: y C1e2x C2e x cos 3x C3e x sin 3x .

119

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]