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

Министерство образования и науки Российской Федерации Московский государственный технологический университет «СТАНКИН»

А.В. Боголюбов, О.К. Иванова

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

для студентов 3 семестра

Москва 2007

УДК 517 (075.8)

Задачи и контрольные вопросы по математике для студентов 3 семестра.

Контрольные задания /Боголюбов А.В., Иванова О.К. под ред. А.В. Боголюбова, А.Г. Елькина, Н.Н. Холщевниковой. – М.: МГТУ «Станкин», «Янус-К», 2007. – 100 с.

Содержит задачи и контрольные вопросы по курсу математики для третьего семестра и включает следующие разделы: дифференциальные уравнения и системы, числовые и функциональные ряды, операционное исчисление. Во всех параграфах даются необходимые краткие теоретические сведения – определения, формулировки теорем, формулы, – а также приводятся примеры решений типовых задач. Задачи сопровождаются ответами. Главы 8 и 10 написаны Боголюбовым А.В., глава 9 – Ивановой О.К.

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

Рис. 8. Библ. 7 назв.

© МГТУ «Станкин», 2007

ОГЛАВЛЕНИЕ

 

Предисловие......................................................................................................

5

Глава 8

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

6

§ 1.

Основные понятия. Уравнения 1-го порядка..................................

6

1.Основные понятия……..…………...……...……………………..6

2.Уравнения 1-го порядка. Графический метод построения интегральных кривых (метод изоклин)…….....……...…...……7

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

4.Однородные уравнения…………..…….......….….……………12

5.Линейные уравнения. Уравнения Бернулли…....…...………..13

6.Особые решения. Уравнения Клеро…………....…...………...15

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

1-го порядка…………...……….……………………………….16

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

19

1.Основные понятия…………………………………...…………19

2.Уравнения, допускающие понижение порядка……...……….21

3.Линейные уравнения……………………...……………………23

4.Однородные линейные уравнения

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

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

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

§ 3. Системы дифференциальных уравнений

......................................33

1.Основные понятия……….……………..………………………33

2.Сведение уравнения n-го порядка к нормальной системе.

Метод исключения……….………………...…………………..35 3. Собственные значения и собственные векторы матрицы.......36

4.Решение однородных линейных систем с постоянными коэффициентами методами линейной алгебры…........………38

§ 4. Элементы теории устойчивости

....................................................42

1.Основные понятия…………………………...…………………42

2.Классификация точек покоя однородной линейной системы 2-го порядка с постоянными коэффициентами......…..………42

Глава 9 . Ряды и их применение……………........……………………….47

§1. Числовые ряды…………………………………………………….47

1.Сходимость ряда…………...…………….……………………..47

2.Признаки сходимости для рядов

сположительными членами…….………….…………………..51

3.Абсолютная сходимость и условная сходимость.

3

Признак Лейбница….……………...…………………………...58

§2. Функциональные ряды……………………………………………63

1.Область сходимости функционального ряда….....…………...63

2.Равномерная сходимость……………...…..…………………...65

§3. Степенные ряды………………………...………………………....68

1.Область сходимости степенного ряда………...……………....68

2.Разложение функций в ряд Тейлора (Маклорена)…...………71

3.Применение степенных рядов…...………...………………….75

§4. Ряды Фурье………………………...………………………………79

Глава 1 0 . Операционное исчисление………..…………………….……84

§1. Преобразование Лапласа………...………………………………..84

1.Определение и свойства преобразования Лапласа…………...84

2.Восстановление оригинала по изображению……....………....87

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

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

Ответы......................................................................................

91

Список литературы……………………....……….………………………….101

4

ПРЕДИСЛОВИЕ

Настоящее пособие предназначается для студентов третьего семестра втузов, изучающих курс общей математики в объеме 270–280 аудиторных часов, что соответствует программам по математике подготовки бакалавров техники и технологий и инженеров-экономистов на дневных факультетах МГТУ «Станкин». Содержит следующие разделы: дифференциальные уравнения и системы, числовые и функциональные ряды, операционное исчисление.

Каждый параграф начинается с краткого теоретического введения, в котором приводятся определения и обозначения математических понятий, математические факты, формулировки теорем и утверждений, а также формулы, необходимые для решения задач. Задачам, предлагаемым для самостоятельного решения, предшествуют подробно разобранные типовые примеры. Начало решений разобранных примеров и задач помечается знаком 3, конец решения – знаком 4.

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

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

5

Глава 8

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

§1. Основные понятия. Уравнения 1-го порядка

1.Основные понятия. (Обыкновенным) дифференциальным уравнением назы-

вается равенство вида

 

) = 0 ,

(1)

F( x, y, y , , y

(n)

 

 

 

связывающее между собой независимую переменную x, искомую функцию y(x) и ее производные y, y′′, , y(n) до производной некоторого порядка n включительно.

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

Решением (частным решением) уравнения (1) на интервале I называется всякая функция y = ϕ(x) , при подстановке которой в это уравнение вместе с ее производны-

ми, уравнение обращается в тождество относительно x I . (В более широком смысле функцию, обладающую указанным свойством на любом открытом множестве, в практике также называют решением). Если решение дифференциального уравнения задается неявно уравнением Φ(x, y) = 0 , то последнее равенство называют интегралом (ча-

стным интегралом) уравнения (1).

График всякого решения дифференциального уравнения (или кривая на плоскости xy, заданная его интегралом) называется интегральной кривой этого уравнения.

Пример 1. Показать, что кривая, заданная уравнением y = 3x sin 2x , является интегральной кривой дифференциального уравнения y′′+ 4 y =12x .

3Имеем y = 3x sin 2x , y′ = 3 2cos 2x , y′′ = 4sin 2x . Подставив y′′ и y в данное уравнение, получим 4sin 2x + 4(3x sin 2x) =12x или 12x =12x , т.е. тождество от-

носительно x R. Следовательно, функция

y = 3x sin 2x

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

ее график – интегральной кривой данного дифференциального уравнения.4

 

 

 

 

 

 

 

8.1.

Каков

порядок

 

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

 

 

уравнения

y

4

y

′′

 

′′′

 

 

 

 

 

 

 

 

 

 

 

 

 

+ sin(xy ) = 0 ?

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.2.

Каков порядок дифференциального уравнения xy

3

 

′′ 4

= y

?

 

 

 

 

 

( y )

 

 

 

 

 

8.3.

Является

ли

функция

y = 2x ln x

решением

уравнения

xy′= 2x + y ?

 

 

 

y = xex

 

 

 

 

 

 

 

 

 

 

 

8.4.

Является

ли

функция

решением

 

уравнения

y′′= 2( y′− y) ?

 

 

x2 + 2 y 3 = 0

 

 

 

 

 

 

 

 

 

 

 

8.5.

Является

ли

парабола

интегральной

кривой

уравнения (2 y + x 3) y′′+ xy′ = 0?

xy = 5 интегральной кривой уравнения

 

 

 

 

8.6.

Является ли гипербола

xy′+ y = 0 ?

 

 

 

 

 

 

 

 

 

 

 

6

 

Задачи повышенной сложности

 

8.7.

Показать, что уравнение ln C(x2 + y2 ) = 2arctg

y

, где C > 0 , яв-

x

 

 

 

 

ляется интегралом дифференциального уравнения ( x y) y′ = x + y .

8.8.

Показать, что линия y 1 x2 x

1 y2 = C

является инте-

гральной кривой уравнения 1 x2 y′= 1y2

для любого C [–1, 1].

2. Уравнения 1-го порядка. Графический метод построения интегральных кривых (метод изоклин). Дифференциальное уравнение 1-го порядка, разрешенное относительно производной, имеет вид:

y′ = f (x, y) .

(2)

Начальным условием для уравнения (2) называется равенство вида

 

y(x0 ) = y0

(3)

(используется также запись y |x=x0 = y0 ), где x0 , y0 – заданные числа (начальные значе-

ния). Задача отыскания решения уравнения (2), удовлетворяющего заданному начальному условию (3), называется задачей Коши для этого уравнения. Геометрически задача Коши состоит в отыскании интегральной кривой уравнения (2), проходящей через заданную точку (x0 , y0 ) на плоскости xy.

Теорема (существования и единственности решения задачи Коши для уравне-

ния (2)). Если функция f (x, y) непрерывна вместе с частной производной fy в некото-

рой области D R2 , то для любой точки (x0 , y0 ) D задача Коши для дифференциаль-

ногоуравнения(2) с начальным условием (3) имеетипритом единственноерешение.

Функция y = ϕ(x,C) называется общим решением уравнения (2), если:

1)при любом допустимом значении параметра C она является решением этого уравнения;

2)любое частное решение уравнения (2) (кроме, быть может, отдельных реше-

ний) представимо в виде y = ϕ(x,C0 ) при некотором значении C0 этого параметра (т.е. для любой точки (x0 , y0 ) D (см. теорему) найдется такое значение C0 , что

y0 = ϕ( x0 ,C0 ) ).

Уравнение Φ(x, y,C) = 0 , определяющее общее решение уравнения (2) неявно, называют общим интегралом этого уравнения.

Пример 2. Является ли функция y = Cx2 общим решением дифференциаль-

ного уравнения y′ = 2xy ? Решить задачу Коши для этого уравнения при начальном ус-

ловии y(1) = −2 .

3При любом допустимом значении C ( C 0 ) данная функция является решением данного уравнения:

 

 

2

2 Cx2

 

y =

Cx

= 2 Cx 2 Cx = x

2 Cx = 2 Cx ( x 0) ,

 

y

т.е. первое условие в определении общего решения выполнено. Второе условие, однако, не выполняется, т. к. при y0 < 0 и x0 0 не существует решения уравнения

7

y

0

=

Cx 2

относительно C (в семействе парабол y = Cx2

нет кривых, проходящих в

 

 

0

 

 

нижней полуплоскости!), в то время как по теореме существования и единственности (условия которой выполняются как для правой полуплоскости D1 ( x > 0 ), так и для ле-

вой D2 ( x < 0 )) через любую точку (x0 , y0 ) ,

где x0 0 , проходит единственная инте-

гральная кривая. Таким образом, функция

y =

Cx2

общим решением не является.

Общим решением в данном случае

является,

очевидно, функция y = Cx2 , C R

(ветви этих парабол уже целиком заполняют и правую и левую полуплоскости). Для решения задачи Коши, положив x =1 , y = −2 , из уравнения 2 = C 12 найдем C = −2 , и, таким

образом, получимискомоечастноерешение y = −2x2 ( x > 0 ). Ответ: нет; y = −2x2 .4

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

8.9.

xy′ =1+ xex .

8.10. y′= 6cos3x 2x .

8.11.

Показать, что функция

y = sin x 1 +Cesin x является общим

решением уравнения y′+ y cos x = sin x cos x .

8.12. Является ли функция y = Cx22 общим решением уравнения xy′+2 y = 0 ?

8.13.Является ли равенство y(1) = 2 начальным условием для уравнения y′= f ( x, y) ?

8.14.Является ли равенство y(1) =1 начальным условием для уравнения y′= x2 y2 ?

8.15. Найти функцию, принадлежащую семейству y = 2 +C cos x и удовлетворяющую условию y(0) = −1.

8.16. Найти кривую семейства y =1 +xCx , проходящую через точку

(–1, 1).

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

8.17. y′= 2x 3

x ; y(1) = 0 .

8.18. (1 + x2 ) y′ = 2x ; y(0) = −2 .

Если G – область определения функции

f (x, y) в дифференциальном уравнении

(2), то в каждой точке

(x, y) координатной

плоскости, где (x, y) G , равенство

tg α = f (x, y) определяет направление касательной к проходящей через эту точку инте-

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

8

Для построения поля направлений целесообразно выбирать точки в множестве G не случайно, а выделять сразу все точки, в которых направление поля одинаково. Множество всех точек плоскости, в которых направление поля, определяемое уравнением (2), одно и то же, называется изоклиной. Всякая изоклина задается уравнением f (x, y) = k , где k – некоторое число, и обычно является некоторой кривой. Построив

несколько изоклин для различных значений k и указав в точках каждой изоклины соответствующее направление поля, определяемое равенством tg α = k , получим поле на-

правлений. Этот способ построения поля направлений, а затем и интегральных кривых, носит название «метод изоклин».

Пример 3. Для уравнения y′ = y + 2 x2 найти изоклину, проходящую через точку (1, 0) , и выяснить направление поля в ее точках.

3Уравнение семейства изоклин имеет вид y + 2 x2 = k .

Положив x = −1, y = 0 , найдем k =1, и, следовательно, уравнение искомой изоклины y + 2 x2 =1, т.е. изоклиной является парабола y = x2 1. Поскольку k =1, то направление поля в точках изоклины образует угол α = arctg1 = 45° с осью x. Ответ:

y= x2 1; α = 45°.4

8.19.Составить уравнение касательной к интегральной кривой уравнения y′= xy2 + x3 в точке (–1, 1).

8.20.

Под каким углом интегральные кривые уравнения y′= xy +3

пересекают ось x?

8.21.

Построить изоклину уравнения y′= x2 + y2 2 y , проходящую

через точку (–2, 1).

8.22. Построить изоклину уравнения (x y) y′ = x , проходящую через точку (1, –1).

Методом изоклин построить приближенно семейство интегральных кривых следующих дифференциальных уравнений:

8.23.

y′= x + y .

8.24. yy′+ x = 0 .

8.25.

Задачи повышенной сложности

Показать, что уравнение

2x + y 1 = Ce2 yx является общим

интегралом дифференциального уравнения (4x + 2 y 3) y′ = 2x + y +1.

8.26. Показать, что функция y = xexxdx является общим решением уравнения xy′− y = xex .

ку y′=

y =

8.27. Найти и построить на плоскости xy область, через каждую точкоторой проходит единственная интегральная кривая уравнения

ln x +ln y .

8.28. Доказать, что функция y = x ln(Cx) +1 является, а функция

x(C +ln x) +1 не является общим решением уравнения y′= x + y +1. x

9

3. Уравнения с разделяющимися переменными. Дифференциальное уравне-

ние 1-го порядка, приводящееся к виду

 

y′ = f (x)g( y) ,

(4)

называется уравнением с разделяющимися переменными. Так как

 

y′ = dydx ,

(5)

то путем деления на g( y) и умножения на dx (при g( y) 0 ) уравнение (4) приводится к виду gdy( y) = f (x)dx (этот переход называется разделением переменных). Интегрируя

обе части этого равенства (левую часть по y, а правую по x), получаем общий интеграл данного уравнения в виде G( y) = F(x) +C , откуда, если это возможно, выражаем об-

щее решение. Кроме этого, если какое-либо число y0 – корень функции g( y) , т.е. g( y0 ) = 0 , то, очевидно, функция y = y0 также является решением уравнения (4).

Замечание. Дифференциальное уравнение 1-го порядка, разрешенное относительно производной (уравнение (2)), может быть задано в дифференциалах:

M (x, y)dx + N (x, y)dy = 0

(6)

(переход от (2) к (6) и обратно осуществляется с помощью (5)). Поскольку в уравнении

(6) переменные x и y равноправны, то решениями этого уравнения наряду с функциями вида y = ϕ(x) называются функции x = ψ( y) , обращающие его в тождество относи-

тельно y на некотором интервале. Уравнение (6) является уравнением с разделяющимися переменными, если M (x, y) = M1( x)M2 ( y) и N (x, y) = N 1( x)N2 ( y) (переменные разделяются путем деления на N 1(x)M2 ( y) ); в этом случае среди прочих решений оно

может иметь решения вида x = x0 , если N1(x0 ) = 0 , и y = y0 , если M2 ( y0 ) = 0 . Пример 4. Решить уравнение

(2 y2 +1) y′ = 2xy .

3Данное уравнение является уравнением с разделяющимися переменными (ср. с (4) при f (x) = 2x , g( y) = 2 y2y+1). Полагая y′ = dydx , разделяем переменные:

(2 y2 +1)

dy

= 2xy

2 y2 +1

dy = 2xdx ,

dx

y

 

 

 

иинтегрируем:

(2 y + 1y)dy = 2xdx y2 + ln | y |= x2 +C (общий интеграл уравнения);

кроме того, при разделении переменных потеряно решение y = 0 , которое не может быть получено из общего интеграла ни при каком значении C. Таким образом, все решения данного уравнения задаются неявно общим интегралом y2 +ln | y |= x2 +C , C R, (вы-

разить

отсюда общее решение в явном виде невозможно) и уравнением y = 0 .

Ответ:

y2 +ln | y |= x2 +C ; y = 0 .4

Пример 5. Найти частное решение уравнения y′+ y = 0 , удовлетворяющее начальному условию y(0) = 3 .

3Найдем сначала общее решение. Разделяем переменные:

10

dydx = −y dyy = −dx .

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

dy

= −

dx ln | y |= −x + ln C ( C > 0 ).

y

 

1

1

 

 

 

Тогда

ln | y |= ln ex + ln C1 ln | y |= ln C1ex | y |= C1ex y = ±C1ex .

Решение y = 0 , которое теряется при разделении переменных, не входит в полученное общее решение, т. к. параметр ±C1 может иметь любые значения, кроме 0. Введем новый параметр C, принимающий также и нулевое значение. Тогда функция y = Cex будет общим решением, включающим и решение y = 0 .

Для нахождения частного решения, положив в общем решении x = 0 , y = 3, из уравнения 3 = Ce0 найдем C = 3 , и, следовательно, y = 3ex .4

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

8.29. 2 yy′+sin x = 0 .

8.31. y′= y tg x .

8.33. (3y2 +1) y′−2x = 0.

8.35.

y′= 2x

1y2 .

8.37.

xydx +

1 x2dy = 0 .

8.30.

yy′ = x .

8.32.

xy′+ 2 y = 0 .

8.34.

( y +1) y′= y .

8.36.

1 x2 y′ = 2 y .

8.38.

2(1 + y2 )xdx (1 + x2 )dy = 0 .

Решить задачи Коши:

 

 

 

 

 

8.39.

(x +1) y′+ xy = 0;

y(0) = 2 .

8.40.

y′= ex+y ; y(0) = 0 .

8.41.

(2xy2 + x)dy +(2x2 y y)dx = 0 ; y(1) =1.

8.42.

y′=

sin x

y(π) = 0.

 

 

 

 

 

;

 

 

 

 

2 +cos y

 

 

 

 

8.43.

xy′ = y y2 ;

y(3) =1.

 

 

8.44.

y′= y2 tg x ; y(1) = 0.

 

 

Задачи повышенной сложности

8.45.

Найти общее решение уравнения

1 x2 y′ =1 + y2 .

8.46.

Решить уравнение y2 cos xdx sin2 xdy = 0 .

Решить задачи Коши:

 

 

 

 

 

8.47.

3 xy′ = 3 + y ;

y(2) = −3.

 

 

8.48.

2sin xdy + ydx = 0 ;

y(π

2

) =1.

 

 

 

 

 

 

 

 

 

 

11

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

 

y

 

y

= g( x)

(7)

 

(или к виду y′ = f (x, y) , где функция f (x, y) является однородной 0-го порядка, т.е. удовлетворяет условию f (λx,λy) = f (x, y) тождественно относительно x, y и λ > 0 ).

С помощью подстановки y = u x (при этом y′ = ux + u ), где u = u(x) – новая

неизвестная функция, однородное дифференциальное уравнение 1-го порядка приводится к уравнению с разделяющимися переменными относительно функции u(x) . Най-

дя все решения этого уравнения и подставив в них u = xy , получим все решения данного

уравнения.

Замечание. Уравнение (6) является однородным, если функции M (x, y) и N (x, y)

однородные

функции одного порядка,

т.е. существует

такое k Z,

что

M (λx,λy) = λk M ( x, y) и N (λx,λy) = λk N (x, y)

тождественно относительно x, y и λ > 0 .

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

 

 

 

 

( x2 + 2 y2 )dx xydy = 0 .

 

 

3Данное

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

т. к. функции

M ( x, y) = x2 + 2 y2

и

N (x, y) = −xy являются однородными 2-го порядка:

 

 

M (λx,λy) = (λx)2 + 2(λy)2 = λ2 (x2 + 2 y2 ) = λ2 M ( x, y) , также и N (λx,λy) = λ2 N ( x, y) (или по другому – приведя уравнение к виду (2), разрешенному относительно производной:

dy

=

x2 + 2 y2

 

(при xy 0 ), получаем, что функция

f (x, y) =

x2 + 2 y2

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

dx

xy

 

 

 

xy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (λx,λy) = f (x, y) ; либо, наконец, f (x, y) =

1 + 2(

y

)2

= g(

y

) , где g(t) =

1 + 2t2

).

вию

x

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

x

 

t

 

 

 

Положим y = ux . Тогда

y

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

= u x + u .

y

 

x2 + 2 y2

 

x2

+ 2(ux)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

xy

,

получим u x + u =

 

 

ux2

,

или,

после преобразований,

уравнение с

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

x dudx =1 +uu2 .

Разделяем переменные:

 

udu

= dx

 

1 + u2

x

и интегрируем:

ln(1 + u2 ) = 2 ln | x | +ln C

(произвольная постоянная для удобства потенцирования представлена в логарифмической форме). Отсюда

1 + u2 = Cx2 .

Возвращаясь к функции y ( u = xy ), получим общий интеграл: x2 + y2 = Cx4 , C > 0 ,

12

откуда можно выразить общее решение y = ±x Cx2 1. Кроме этого, исходное уравнение имеет решение x = 0 , которое теряется при разрешении его относительно y. От-

вет: y = ±x Cx2 1 ( C > 0 ); x = 0 .4

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

 

 

x

 

y

 

 

y′=

y

8.49.

y′ = y

+

 

.

 

 

8.50.

 

.

 

 

x

 

x y

8.51.

y2dx + x2dy = xydy .

 

8.52.

(x y)dx + xdy = 0 .

8.53.

xy′− y = x sin y x .

 

8.54.

y = x( y′−e

y

 

 

) .

 

x

Решить задачи Коши:

 

 

 

 

 

 

 

8.55.

xy′ = y + x tg y x ; y(12) = π4 .

8.56.

xy′ = y ln y x ; y(1) =1.

8.57.

( y +

x2 + y2 )dx xdy = 0; y(1) = 0 .

 

 

 

 

 

8.58.

(x2 + y2 )dy xydx = 0 ;

y(0) = 0,5.

 

 

 

 

 

8.59.

 

 

 

Задачи повышенной сложности

Решить уравнение y2dx = x2 (dy + 2dx) .

Решить задачи Коши:

 

 

 

 

 

 

 

8.60.

(xy′− y) arctg y x = x ;

y(1) = 0 .

 

 

 

 

 

 

 

y′=

y2 2xy x2

 

 

 

 

 

 

 

8.61.

 

; y(1) = −1.

 

 

 

 

 

 

y2 + 2xy x2

 

 

 

 

 

 

8.62. Найти функцию ϕ(t) , такую, чтобы общее решение уравнения xϕ( x y) = xy′− y имело вид y = lnxCx (C 0 ).

5. Линейные уравнения. Уравнения Бернулли. Дифференциальное уравнение

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

 

 

 

 

 

y′+ a( x) y = b(x)

(8)

(т.е. содержит y и yв первой степени).

 

 

 

 

Для решения линейного уравнения положим y = u v , где u и v – функции от x.

Тогда

y

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

 

 

= u v + uv

 

 

 

 

 

 

 

+ a(x)v) = b( x) .

(9)

 

 

 

 

 

u v +u(v

 

Выберем функцию v(x) так, чтобы выражение в скобке обратилось в нуль, т.е. найдем любое ненулевое решение уравнения с разделяющимися переменными v′ = −a(x)v .

Подставляя эту функцию вместо v в уравнение (9), получаем уравнение с разделяющимися переменными относительно функции u(x) :

uv(x) = b(x) .

13

Находим общее решение этого уравнения u = u(x) +C , после чего, перемножая найденные функции v(x) и u(x) +C , получаем общее решение исходного уравнения (8):

y = v(x)(u(x) +C) .

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

xy′ = cos x y .

3Данное уравнение является линейным:

 

 

 

 

 

 

 

 

y′+

y

=

cos x

.

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x

 

 

 

 

 

Полагая здесь y = uv ,

y

 

 

, получим

 

 

 

 

 

 

 

 

 

= u v + uv

 

 

 

 

 

 

 

 

 

 

 

 

uv

cos x

 

 

 

v

cos x

 

 

 

u v + uv

 

+ x

=

x

или

u v + u(v

 

+ x) =

x

.

(*)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

Найдем какую-нибудь функцию v(x)

( 0 ), решая уравнение v′+ x = 0 :

dv = −v

 

dv

= −dx

ln | v |= −ln | x | (C = 0) v =

1 .

dx

x

 

v

 

 

 

x

 

 

 

 

 

 

 

 

 

 

x

Подставляя v(x) в уравнение (*), получим:

ux= cosx x или dudx = cos x .

Общее решение этого уравнения:

u = sin x +C .

Перемножая теперь u и v, получаем общее решение исходного уравнения:

y = sin x +C .4 x

Аналогично линейному уравнению (с помощью подстановки y = u v ) решается

уравнение Бернулли:

y′+ a(x) y = b(x) yα , где α ≠ 0 , α ≠1

(при α = 0 это уравнение является линейным, а при α =1 – уравнением с разделяющимися переменными). Уравнение Бернулли можно свести к линейному с помощью подстановки z = y1−α . Заметим, что при α > 0 у уравнения Бернулли есть решение y = 0 .

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

 

 

 

 

 

2 y

 

3y

8.63.

y

 

= ex (x +1)2 .

8.64.

y′− x = x .

x +1

8.65.

y′= 3x2ex2 2xy .

8.66.

y′= e3x 2 y .

8.67.

xy′ = y + x2 cos x .

8.68.

xy′ = ex + xy .

Решить задачи Коши:

 

 

8.69.

(1 + x2 ) y′ = 2xy +(1 + x2 )2 ;

y(1) = 4 .

 

 

 

 

1

 

 

 

8.70.

y

+ y tg x = cos x ; y(0) = 0 .

 

 

 

 

 

14

Задачи повышенной сложности

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

y

8.72. (2ey x) y′ =1.

8.71. y′ = x +2 y3 .

8.73. y′+ 4xy = 2xex2

8.75. xy′ = x(2ln x +1)

Решить задачу Коши:

8.76. y′= 2 y + ex x ;

y .

8.74. xy2 y′ = x2 + y3 .

y .

 

y(0) = 0,25 .

 

6. Особые решения. Уравнения Клеро. Решение или интеграл дифференциального уравнения F(x, y, y) = 0 называется особым, если через любую точку соответ-

ствующей (особой) интегральной кривой проходит, кроме нее, еще и другая касающаяся ее интегральная кривая. Для уравнения y′ = f (x, y) , в частности, в каждой точке

особой интегральной кривой нарушается единственность решения задачи Коши. Если однопараметрическое семейство интегральных кривых дифференциального уравнения, определяемое уравнением Φ(x, y,C) = 0 , имеет огибающую, то она является особой ин-

тегральной кривой, уравнения которой (явное или неявное) получаются путем исключения (если это возможно) параметра C из системы

Φ( x, y,C) = 0,

(10)

{Φ′C (x, y,C) = 0

 

(где Φ′C ( x, y,C) = ∂ΦC – частная производная по C функции Φ(x, y,C) ), либо путем ре-

шения этой системы относительно x и y (параметрические уравнения с параметром C). Пример 8. Решить уравнение y′ = 33 y2 .

3Данное уравнение – с разделяющимися переменными. Поэтому dydx = 33 y2 33dyy2 = dx 3 y = x +C

и y = (x +C)3 – его общее решение (т.е. семейство интегральных кривых состоит из кубических парабол). Имеется, однако, решение y = 0 , которое не может быть получе-

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

Коши в точках оси x – нарушение непрерывности частной производной fy = 32y в этих

точках (см. п. 2). Заметим, что решение y = 0 может быть получено в результате нахождения огибающей полученного семейства парабол, т.е. путем исключения параметра C из системы y = (x +C)3 , 0 = 3(x +C)2 . Ответ: y = (x +C)3 ; y = 0 .4

Примером дифференциального уравнения 1-го порядка, имеющим особые ре-

шения, является уравнение Клеро:

y = xy′+ϕ( y) .

15

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

y = Cx (C) ,

(11)

а особые решения согласно (10) определяются системой

y = Cx (C) , 0 = x + ϕ (C) .

Особая интегральная кривая, определяемая этой системой, является огибающей семейства прямых (11), т.е. это семейство состоит из ее касательных.

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

4 y + y2 = 4( x +1) y.

3Выражая отсюда y:

y = (x +1) y′− y42 или y = xy′+ y′− y42 ,

получаем, что данное уравнение есть уравнение Клеро. Следовательно, его общее решение имеет вид

y = Cx +C C42 .

Для нахождения особого решения продифференцируем обе части этого равенства по C: 0 = x +1 C2

и составим систему

 

 

C2

 

y = Cx +C

4

,

 

= x +1 C .

 

 

0

 

 

 

2

 

 

Исключая отсюда C (из второго уравнения выражаем C = 2(x +1) и подставляем в пер-

вое уравнение), получим особое решение y = (x +1)2 . Ответ: y = Cx +C

C2

;

4

y = (x +1)2 .4

 

 

 

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

8.77.

xy2 + 2xy′− y = 0 ; ( y C)2 = 4Cx .

 

 

 

 

 

8.78.

y = y2 xy′+0,5x2 ;

y = 0,5x2 +Cx +C2 .

 

 

 

 

Решить уравнения Клеро:

 

 

 

 

 

 

8.79.

y = xy′+

 

1

 

8.80.

y = xy′−ey.

 

 

 

y.

 

 

 

8.81.

y

= e

yxy

.

 

8.82.

y +

y

= 0 .

 

 

 

 

y (x +1)

 

7. Смешанные задачи. Геометрические и физические задачи, приводящие к решению дифференциальных уравнений 1-го порядка.

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

 

 

 

 

8.83. xy

y = x ctg

y

x .

8.84.

+ y

2

.

 

 

y = y (x +1)

 

16

8.85.

y = xy′+ 1 + y2 .

Решить задачи Коши:

8.87.

y

x

y(0) = 0,5.

= −y2e2 ;

8.89.

y′+ y tg x = −cos x ; y(1) = 0 .

8.91.

y

xy y2

; y(1) = 0,5.

= x2

 

8.86.

y′= 2e2 x+y .

 

8.88.

 

 

 

 

ex

π

 

 

 

 

 

 

 

 

; y( 2) = 0 .

y +y =sin2 x

8.90.

y

=

2x + y

;

 

y(1) = 0 .

 

 

 

x

 

8.92.

y

=

2xy

;

 

y(1) =1.

x2 +4

 

 

Задача из геометрии о нахождении уравнения линии y = y(x) по заданному

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

терминахкоординатточки касания( x и y ) иуглового коэффициента( yили y1).

Пример 10. Найти уравнение линии, проходящей через точку (0, 2), если поверхность, образованная при вращении этой линии вокруг оси x (зеркало прожектора), отражает все лучи, выходящие из начала координат O, параллельно оси x (в положительном направлении).

3Пусть y = y(x) – уравнение искомой линии, M (x, y) – произвольная точка на ней, MA – нормаль в этой точке, A – точка пересечения нормали с осью x . Тогда, ввиду закона отражения из оптики, линия определяется свойством ее нормали: угол OMA между прямой OM и нормалью (угол падения) равен углу OAM между прямой x и нормалью (угол отражения), что эквивалентно равнобедренности треугольника AMO . Выразим равенство OA = OM через x , y и y. Уравнениенормалив точке M имеетвид:

Y y = − y1( X x)

(здесь X и Y – текущие координаты нормали). Полагая в этом уравнении YA = 0 , находим X A = x + yy. Таким образом, OA =| x + yy| . Поскольку OM = x2 + y2 , то получаем равенство | x + yy|= x2 + y2 , представляющее собой дифференциальное уравнение 1-го порядка, которому удовлетворяет искомая функция y = y(x) . Уравнение при-

водится к виду

 

 

 

y′ =

x ±

x2 + y2

,

 

y

 

 

 

т.е. является однородным. Решая его (с помощью подстановки y = ux ), найдем общий интеграл:

y2 = C2 + 2Cx .

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

Для выделения линии, проходящей через данную точку, используем начальное условие y(0) = 2 , откуда C = 2 , и, следовательно, y = 2 1 + x – уравнениеискомойлинии.4

17

8.93. Найти и построить линию, проходящую через точку (12 , 1) и

обладающую следующим свойством: абсцисса точки пересечения касательной в произвольной точке линии с осью x равна квадрату абсциссы точки касания.

8.94. Найти и построить линию, проходящую через точку ( 5, 1) и

обладающую следующим свойством: отрезок нормали к ней в произвольной точке M, заключенный между осями координат, делится точкой M пополам.

8.95.Найти уравнение линии, проходящей через точку (1, 2) и обладающей следующим свойством: произведение абсциссы произвольной точки M линии на абсциссу точки A пересечения нормали MA с осью x равно удвоенному квадрату расстояния от начала координат до точки M.

8.96.Найти уравнение линии, проходящей через точку (12 ,7 4) , если

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

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

dtd (mv) = F(t, v) для прямолинейного и dtd (Iω) = M (t, ω) для вращательного движения

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

Пример 11. В помещении цеха вместимостью 10800 м3 воздух содержит 0,12% углекислоты. Вентиляторы доставляют свежий воздух, содержащий 0,04% углекислоты, со скоростью 1500 м3/мин. Предполагая, что углекислота распределяется по помещению равномерно в каждый момент времени, найти процентное содержание углекислоты через 10 мин после начала работы вентиляторов.

3Пусть x = x(t) – процентное содержание углекислоты через t мин после начала

работы вентиляторов, x + x – в момент времени t + t . За малый промежуток времени t вместе со свежим воздухом в помещение поступило 1500 t 0,04% углекислоты и приблизительно (с точностью до o( t) ) 1500 t x% углекислоты удалено вместе с воздухомпомещения(воздухиз помещения цеха удаляетсястойжескоростью). Такимобразом

x 1500 t 0,04 1500 t x . 10800

Заменяя приращения x и t дифференциалами, получаем дифференциальное уравнение 1-го порядка (с разделяющимися переменными)

36dx = 5(0,04 x)dt .

18

Его общее решение x = 0,04 +Ce365 t . Используя начальное условие x(0) = 0,12 , находим C = 0,08 и, следовательно, зависимость процентного содержания углекислоты в воздухе цеха от времени имеет вид:

x = 0,04(1 + 2e365 t ) .

Отсюда при t =10 получаем x 0,06% .4

8.97.Замедляющее действие трения на диск, вращающийся в жидкости (момент сил трения), пропорционально угловой скорости вращения. Начав вращение со скоростью 5 об/c, по истечении двух минут диск вращается со скоростью 3 об/c. Через какое время от начала движения его угловая скорость уменьшится до 1 об/c?

8.98.Скорость охлаждения тела пропорциональна разности температур тела и окружающей его среды (закон Ньютона). Через какое время температура тела, нагретого до 100°C, понизится до 25°C, если температура помещения равна 20°C и за первые 10 мин тело охладилось до 60°C?

8.99.В дне котла, имеющего форму полушара радиусом R = 43 см, образовалась пробоина площадью S = 0,2 см2. Через какое время вода, наполняющая котел, вытечет из него? (По закону Торричелли вода вытекает

из сосуда через малое отверстие со скоростью v = 0,6 2gh , где h – высота

уровня воды над отверстием, g = 9,81 м/c2 – ускорение свободного падения). 8.100. Пуля, двигаясь со скоростью v0 = 400 м/с, пробивает стену толщиной h = 20 см и вылетает, имея скорость v1 =100 м/с. Полагая силу сопротивления стены пропорциональной квадрату скорости движения пу-

ли, найти время прохождения пули через стену.

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

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

 

 

y

 

 

 

 

) .

(1)

 

 

(n)

= f ( x, y, y , , y

(n1)

 

 

 

 

 

 

 

 

Начальными условиями для уравнения (1) называются n равенств вида

 

 

y(x0 ) = y0 , y(x0 ) = y0, , y(n1) ( x0 ) = y0(n1)

(2)

(или y |x=x = y0 ,

y|x=x

= y0, ,

y(n1) |x=x

0

= y0(n1) ),

где x0 , y0 , y0, , y0(n1)

– заданные

0

 

0

 

 

 

 

 

 

числа (начальные значения).

Задача отыскания решения уравнения (1), удовлетворяющего заданным начальным условиям (2), называется задачей Коши для этого уравнения.

19

Теорема (существования и единственности решения задачи Коши для уравне-

ния (1)). Если функция f ( x, y, y, , y(n1) ) непрерывна вместе с частными производными fy , yf,, y(nf1) в некоторой области D Rn+1 , то для любой точки

(x0 , y0 , y0, , y0(n1) ) D задача Коши для дифференциального уравнения (1) с начальными условиями (2) имеет и притом единственное решение.

Функция y = ϕ(x,C1,C2 , ,Cn ) называется общим решением уравнения (1), если: 1) при любых допустимых значениях параметров C1,C2 , ,Cn она является ре-

шением этого уравнения; 2) любое его частное решение (кроме, быть может, отдельных решений) пред-

ставимо в виде y = ϕ(x,C1,C2 , ,Cn ) при некоторых значениях параметров (т.е. для любой точки (x0 , y0 , y0, , y0(n1) ) D (см. теорему) найдутся такие значения параметровC1,C2 , ,Cn , что функция y = ϕ(x,C1,C2 , ,Cn ) будет удовлетворять начальным условиям (2)).

Уравнение Φ(x, y,C1, ,Cn ) = 0 , определяющее общее решение уравнения (1) неявно, называют общим интегралом этого уравнения.

Пример 1. Показать, что функция y = C1ex +C2ex является общим решением уравнения y′′− y = 0 . Найти частное решение этого уравнения, удовлетворяющее начальным условиям y(0) = 0 , y(0) =1.

3Проверим выполнение условий из определения общего решения.

1)Число произвольных постоянных в данной функции равно 2, что совпадает с порядком уравнения.

2)При любых значениях C1 и C2 эта функция является решением уравнения:

y = C1ex +C2ex y′ = C1ex C2ex y′′ = C1ex +C2ex ;

подставляя y′′ и y в уравнение, получим тождество:

C1ex + C2ex (C1ex + C2ex ) = 0 0 = 0 x R. 3) Каковы бы ни были начальные значения x0 , y0 , y0, система уравнений

{y0 = C1ex0 +C2ex0 , y0′ = C1ex0 C2ex0

имеет, притом единственное решение относительно C1 и C2 , т. к. ее определитель

=

 

ex0

ex0

 

= −2 0 .

 

 

 

 

ex0

ex0

 

 

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

Для нахождения частного решения после подстановки начальных значений x = 0 , y = 0 и y′ =1 в общее решение и его производную получим систему

{0 = C1 +C2 ,

1 = C1 C2 ,

откуда C1 = 12 , C2 = − 12 и, следовательно, искомое частное решение y = ex 2ex .4

20

8.101.

Являются ли равенства

y(2) = −1,

 

= 2 начальными ус-

y

(1)

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

′′

 

 

 

 

 

 

 

 

= f (x, y, y ) ?

 

 

 

 

8.102.

Является ли гипербола

xy = 3 интегральной кривой уравне-

ния yy′′= 2 y2 ?

 

 

 

 

 

 

 

 

 

8.103.

Является ли функция y = Cx2 +C3

общим решением уравне-

ния xy′′ = y?

 

 

 

y = (C1 +C2 )sin x

 

8.104.

Является ли

функция

общим решением

уравнения y′′+ ytg x = 0?

 

 

 

 

 

 

 

 

8.105.

Показать, что функция y = C x3 +C

2

является общим решени-

ем уравнения xy′′−2 y′= 0 .

1

 

 

 

 

 

1

 

 

 

 

 

8.106.

Показать, что функция y =

sin(C x +C

2

) при любых C 0

 

 

 

 

 

1

 

1

 

 

 

 

C1

 

 

 

 

и C2 является решением уравнения yy′′+1 = y2 .

 

 

 

 

2. Уравнения, допускающие понижение порядка.

 

 

 

а) Уравнения вида y(n)

= f (x) . Общее решение этого уравнения получается пу-

тем n-кратного интегрирования (при каждом интегрировании порядок уравнения по-

нижается на 1).

 

 

 

 

 

 

 

 

′′

(не содержащие явно y). С по-

 

 

 

 

 

 

 

 

 

б) Уравнения 2-го порядка вида F(x, y , y ) = 0

мощью замены y

= z(x) , y

′′

= z

 

 

 

 

 

 

порядок уравнения понижается до 1-го: F(x, z, z ) = 0 .

Если z = z(x,C1)

 

– общее решение последнего уравнения, то, интегрируя равенство

y′ = z(x,C1) , получаем общее решение данного уравнения: y = z(x,C1)dx +C2 .

 

 

 

 

 

 

 

 

 

′′

(не содержащие явно x). С по-

в) Уравнения 2-го порядка вида F( y, y , y ) = 0

 

 

 

 

 

′′

 

 

dp

 

 

мощью замены y

= p( y) ,

 

y

= p dy данное уравнение приводится к уравнению 1-го

 

 

 

порядка (относительно функции

p = p( y) ): F( y, p, p dp) = 0 . Если p = p( y,C ) – общее

 

 

 

 

 

 

 

 

 

 

dy

1

 

 

 

 

 

 

 

 

 

 

 

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

y′ = p( y,C1) ,

получаем общий

интеграл

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

dy = p( y,C1)

 

 

dy

 

 

 

 

 

dx

 

 

= dx P( y,C1) = x +C2

(а также решения вида y = y0 , если p( y0 ,C1) = 0 ).

p( y,C )

1

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

xy′′+ 2 y′ = 6x .

 

 

 

 

 

3Данное уравнение 2-го порядка и не содержит y.

Положим y′ = z ; тогда

y′′ = zи уравнение принимает вид

 

2

 

 

 

 

xz′+ 2z =

6x , или z′+ x z = 6 . Это линейное уравне-

ние 1-го порядка. Решая его с

помощью подстановки

z = uv ,

найдем v =

1

,

 

 

 

 

 

 

 

 

 

x2

u = 2x3 +C1 ,

и общее решение z = 2x +C1 . Следовательно, y′ = 2x + C1 , откуда, интег-

 

 

 

 

x2

 

x2

 

 

рируя, получим общее решение данного уравнения y = x2 C1

+C2 .4

 

 

 

 

 

 

 

 

x

 

 

 

 

21

Пример 3. Найти частное решение уравнения yy′′ = y2 , удовлетворяющее начальным условиям y(0) = y(0) = −1 .

3Уравнение 2-го порядка и не содержит x. Положим y′ = p( y) ; тогда y′′ = p dpdy

и уравнение принимает вид p2 dpdy = y2 . Это уравнение с разделяющимися переменны-

ми. Разделяя переменные и интегрируя, найдем его общий интеграл p3 = y3 +C1 . Используя начальные значения y = −1, y′ = p = −1 (при x = 0 ), получим C1 = 0 . Следовательно, y′ = y . Решая это уравнение с разделяющимися переменными, найдем его общее решение y = C2ex , а из условия y(0) = −1 следует C2 = −1 . Итак, искомое частное решение есть y = −ex .4

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

8.107. (1 + x2 )2 y′′+2x = 0 .

8.109. (2x +1) y′′+2 y′ = 0 .

8.111. (1 +ex ) y′′−ex y′ = 0 . 8.113. 2 yy′′=1 + y2 . 8.115. y′′ = 4 y3 y3 .

8.117. ( y 2) y′′= 3y2 .

Решить задачи Коши:

8.108. y′′= 6x +4sin 2x .

8.110. xy′′= 3y.

8.112. y′′+ y2 = 0 . 8.114. ( y 1) y′′ = 2 y2 .

8.116. y′′+12 y2 y3 = 0 . 8.118. yy′′+ y2 = 0 .

8.119.

x

3

y

′′′

= 2 ;

 

y(1) =

1,

 

 

 

′′

=

3 .

 

 

 

 

y (1) =

1, y

(1)

8.120.

2

 

 

xy

′′

 

= 3;

 

y(1) = −1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y (1) = 0 .

 

 

 

8.121. (x 6) y

′′

=

2( y

3) ;

y(7)

= 22 ,

 

 

= 0.

 

 

 

y (7)

8.122. (x 1) y

′′

= y

 

 

 

 

 

 

 

 

 

 

= 2 .

 

 

 

 

 

 

+2 ; y(3) =1, y (3)

8.123. (1 +4x

2

) y

′′

2 y

2

 

2 = 0 ;

y(0) =1,

y

 

 

 

 

 

 

 

(0) = 0 .

8.124.

y

′′

= e

x

y

2

;

 

y(0)

= 2 ,

 

 

 

 

 

 

′′

 

 

 

 

 

y

(0) = −1.

 

 

 

 

8.125.

y

+e

2 y

=

0 ; y(0)

= 0 ,

=1.

 

 

 

 

 

 

 

 

 

 

 

y

(0)

 

 

 

8.126. (1 + y

2

) y

′′

= 2 yy

2

;

y(0) =

0 ,

 

 

 

 

 

 

 

 

 

 

y (0) =1.

Задачи повышенной сложности

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

8.127. x3 y′′+ x2 y′−1 = 0 .

8.129. y′′+ ytg x = sin 2x .

8.131. yy′′+ y3 = y2 . 8.133. 2 yy′′−3y2 = 4 y2 .

Решить задачи Коши:

8.128. xy′′− y′ = x2ex . 8.130. x( y′′−2) = 2 y. 8.132. yy′′ = y2 + y.

8.134. yy′′− y2 = y2 y.

8.135.

y′′tg y = 2 y2 ;

 

y(0) = π2 ,

y(0) = −1.

8.136.

xy

′′

ln

y

x

;

y(1) = 0 ,

= e .

 

= y

 

y (1)

22

8.137. Тело, находящееся от центра Земли на расстоянии H, падает на Землю из состояния покоя под действием силы притяжения с ускорением, обратно пропорциональным квадрату его расстояния от центра Земли (закон Ньютона). Радиус Земли равен R, ускорение на поверхности равно g. Пренебрегая сопротивлением атмосферы, определить, через сколько времени оно упадет на Землю.

8.138. Тело массой m движется прямолинейно под действием постоянной силы F. Найти скорость движения тела и пройденный им путь как функции времени, если в начальный момент они оба равны нулю, а сопротивление среды пропорционально квадрату скорости с коэффициентом пропорциональности k.

3. Линейные уравнения. Линейным дифференциальным уравнением n-го поряд-

ка называется уравнение вида

y(n) + a1(x) y(n1) + … + an1(x) y′+ an (x) y = b(x) ,

(3)

где коэффициенты a1( x), a2 ( x), , an ( x) , а также правая часть b(x) предполагаются

непрерывными функциями на некотором интервале I. Все решения этого уравнения определены на интервале I.

Если b(x) = 0 x I , то уравнение (3) имеет вид

y(n) + a1(x) y(n1) + … + an1( x) y′+ an (x) y = 0

(4)

и называется однородным; в противном случае уравнение (3) называется неоднородным. Множество решений однородного линейного дифференциального уравнения

(о.л.д.у.) обладает следующими свойствами:

1) если y1(x) и y2 (x) – какие-нибудь два решения уравнения (4), то их сумма y1 + y2 также есть решение этого уравнения;

2) если y(x) – какое-нибудь решение уравнения (4) и C – любое число, то их произведение Cy также есть решение этого уравнения.

Совокупность (система) из m функций y1( x), y2 ( x), , ym ( x) называется линейно зависимой на интервале I, если существуют числа α1, α2 ,, αm , не все равные нулю

и такие, что

α1 y1 + α2 y2 +…+ αm ym = 0 x I ;

если же это равенство возможно только в случае α1 = α2 =…= αm = 0 , то система на-

зывается линейно независимой. Две функции y1 и

y2 ( y2 0 ), в частности, тогда и

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

y1

= const .

 

 

 

 

 

 

 

y2

 

 

 

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

системы функций y1(x), y2 (x), ,

ym ( x) называется определитель

 

 

 

 

 

 

 

y1 y2

ym

 

 

 

 

 

 

′ ′

 

 

W (x) =

y1 y2

ym

 

.

 

 

 

y1(m1) y2(m1) ym(m1)

 

23

Если хотя бы в одной точке x0 I имеем W (x0 ) 0 , то система функций линейно неза-

висима на I.

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

Теорема (о структуре общего решения о.л.д.у.). Если y1, y2 , , yn – какая-

нибудь фундаментальная система решений уравнения (4), то его общее решение имеет вид:

y = C1 y1 +C2 y2 + … +Cn yn ,

где C1,C2 , ,Cn – произвольные постоянные.

Для данного неоднородного уравнения (3) однородное уравнение (4) (с той же левой частью) называется соответствующим однородным уравнением.

Множество решений неоднородного линейного дифференциального уравнения (н.л.д.у.) обладает следующими свойствами:

1)если y1(x) – какое-нибудь решение неоднородного уравнения (3), а y0 ( x) – какое-нибудь решение соответствующего однородного уравнения (4), то их сумма y1 + y0 есть решение уравнения (3);

2)если y1(x) и y2 (x) – какие-нибудь два решения уравнения (3), то их разность y1 y2 есть решение соответствующего однородного уравнения (4).

Теорема (о структуре общего решения н.л.д.у.). Общее решение н.л.д.у. равно сумме любого частного решения этого уравнения и общего решения соответствующего о.л.д.у.

Если правая часть н.л.д.у (3) есть сумма нескольких функций b( x) = b1(x) + b2 (x) + … + bk (x)

и yi (x) – какие-нибудь частные решения уравнений

y(n) + a1(x) y(n1) + … + an1(x) y′+ an (x) y = bi (x)

( i =1, 2,, k ) соответственно, то их сумма y1 + y2 + … + yk есть решение уравнения (3)

(принцип суперпозиции решений).

Пример 4. Проверить, что функции x и cos x являются решениями уравнения

(1 + x tg x) y′′− xy′+ y = 0 ,

и найти его решение, удовлетворяющее начальным условиям y(0) = −1, y(0) = 2 . 3Имеем y = x , y′ =1 , y′′ = 0 . Подставив в уравнение, получаем тождество:

(1 + x tg x) 0 x 1 + x = 0 0 = 0 .

Аналогично проверяется вторая функция. Данное уравнение является однородным линейным дифференциальным уравнением 2-го порядка (к виду (4) оно приводится путем деления на коэффициент при y′′; один из интервалов непрерывности коэффициентов, содер-

жащий точку x = 0 , – I = (π2, π2) ). Так как решения x и cos x линейно независимы

( cosx x const ), то они образуют фундаментальную систему решений данного уравнения, и по теоремео структуреобщегорешения о.л.д.у. его общеерешение имеет вид

y = C1x +C2 cos x .

24

После подстановки начальных значений x = 0 , y = −1 и y′ = 2 в общее решение и его производную y′ = C1 C2 sin x найдем C1 = 2 , C2 = −1 , и, следовательно, искомое частное решение y = 2x cos x .4

Исследовать на линейную зависимость следующие системы функций:

8.139.

8.142.

8.145.

ex , xex .

8.140.

2 x , 3ln 2x .

8.141.

ex ,

ex+1 .

x , ln x .

8.143.

x , 2x , x2 .

8.144.

ex ,

xex , x2ex .

sin x , cos x , sin 2x .

 

8.146.

1, sin2 x , cos 2x .

8.147. Функции sin x и cos x являются решениями о.л.д.у. Будет ли разность cos x sin x решением этого уравнения?

8.148. Образуют ли функции ex и ex фундаментальную систему решений уравнения y′′− y = 0 ?

8.149. Образуют ли функции sin 2x и sin x cos x фундаментальную систему решений уравнения y′′cos 2x + 2 ysin 2x = 0 ?

8.150. Функции x2 и x3 образуют фундаментальную систему решений о.л.д.у. Найти решение этого уравнения, удовлетворяющее начальным условиям: y(1) =1, y(1) = 0 .

8.151. Функции ex , e2 x , e3x – решения о.л.д.у. 3-го порядка. Найти его частное решение, удовлетворяющее начальным условиям: y(0) = 6 , y(0) =14 , y′′(0) = 36 .

8.152. Составить о.л.д.у. 2-го порядка, решениями которого являются функции ex и xex .

8.153. Составить о.л.д.у. 2-го порядка, решениями которого являются функции x3 и x4 .

8.154. Дано н.л.д.у. 2-го порядка, причем y = ln x – его частное ре-

шение, а y = tg x и y = cos x

– частные решения соответствующего одно-

родного уравнения. Найти общее решение данного уравнения.

 

4. Однородные линейные уравнения с постоянными коэффициентами. Об-

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

 

 

 

y(n) + a1 y(n1) + … + an1 y′+ an y = 0 ,

(5)

где a1, a2 , , an – некоторые числа.

 

Многочлен n-й степени от переменной λ

 

 

 

L(λ) = λn + a1λn1 + … + an1λ + an ,

 

полученный из

левой

части

уравнения (5) путем замены производных

y(k )

( k = 0, 1, 2, , n ;

y(0) = y

по определению) искомой функции степенями λk , называет-

ся характеристическим многочленом, а алгебраическое уравнение n-й степени

 

 

 

λn + a1λn1 + … + an1λ + an = 0

(6)

характеристическим уравнением для уравнения (5).

 

25

Теорема (о характеристическом уравнении). Функция eλx тогда и только то-

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

Решения уравнения (5), соответствующие различным корням характеристического уравнения, линейно независимы; кроме того, если λ – действительный корень кратности r >1, то наряду с функцией eλx решениями дифференциального уравнения являются функции

xeλx , x2eλx , , xr1eλx ,

причем все эти r решений также линейно независимы, а каждой паре комплексно сопряженных корней λ = α±iβ кратности s соответствует s пар действительных линейно

независимых решений вида

eαx cosβx, xeαx cosβx, , xs1eαx cosβx , eαx sinβx, xeαx sinβx, , xs1eαx sinβx .

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

В частности, вопрос о решении уравнения 2-го порядка

 

 

 

 

 

 

 

y′′+ a1 y′+ a2 y = 0

 

 

 

 

(7)

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

 

 

 

 

 

 

а) характеристическое уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ2 + a λ + a

2

= 0

 

 

 

 

(8)

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

имеет 2 различных действительных корня λ1

и λ2 (дискриминант D > 0 ) – фундамен-

тальную систему образуют решения eλ1x , eλ2x ,

и, следовательно, согласно теореме о

структуре общего решения о.л.д.у. (п. 3) оно имеет вид y = C eλ1x +C eλ2x ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

б)

характеристическое уравнение имеет двукратный действительный корень λ

( D = 0 , λ = λ

2

) eλx

и xeλx

– фундаментальная система решений и y = eλx (C +C

2

x)

 

1

 

 

 

 

 

 

 

 

 

 

 

1

 

– общее решение;

 

 

 

 

 

 

 

 

 

 

 

 

 

в)

корни характеристического уравнения λ1,2 = α±iβ

комплексно сопряжены

( D < 0 )

 

 

eαx cosβx

и

eαx sinβx

фундаментальная

система

решений

 

и

y = eαx (C cosβx +C

2

sinβx) – общее решение.

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

y′′− 4 y = 0 .

 

 

 

 

 

 

 

3Характеристическое уравнение λ2 4 = 0 имеет корни λ = −2 ,

λ

2

= 2 . Фун-

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

даментальную систему решений образуют функции e2 x и e2 x . Следовательно, общее решение имеет вид y = C1e2 x +C2e2 x .4

Пример 6. Найти общее решение уравнения y′′+ 2 y′+5y = 0 .

3Характеристическое уравнение λ2 + 2λ +5 = 0 имеет корни λ1,2 = −1 ± 2i . Следовательно, функции ex cos 2x и ex sin 2x ( α = −1, β = 2 ) образуют фундаментальную систему решений, а общее решение имеет вид y = ex (C1 cos 2x +C2 sin 2x) .4

26

Пример 7. Найти решение уравнения

y(4) + y′′ = 0 ,

удовлетворяющее начальным условиям y(0) = y′′′(0) = −1 , y(0) =1, y′′(0) = 0 . 3Характеристическое уравнение λ4 + λ2 = 0 имеет двукратный корень λ1,2 = 0

и комплексно сопряженные корни λ3,4 = 0 ±i . Поэтому фундаментальная система ре-

шений имеет

вид

1,

x,

cos x , sin x ,

откуда

общее решение –

y = C1 +C2 x +C3 cos x +C4 sin x . Для определения

произвольных постоянных

найдем

производные

 

 

 

 

 

 

 

y′ = C2 C3 sin x +C4 cos x ,

y′′ = −C3 cos x C4 sin x , y′′′ = C3 sin x C4 cos x

и используем

начальные

условия. Получаем:

1 = C1 +C3 ,

1 = C2 +C4 ,

0 = −C3 ,

1 = −C4 , откуда C1 = −1,

C2 = 0 , C3 = 0 и C4 =1 .

Следовательно, искомое частное ре-

шение имеет вид y = sin x 1.4

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

8.155.

y′′−6 y′+9 y = 0 .

 

 

 

8.156.

y′′−4 y′+5y = 0 .

8.157.

3y′′−2 y′−8 y = 0 .

 

 

 

8.158.

4 y′′+ 4 y′+ y = 0 .

8.159.

y′′+6 y′+13y = 0 .

 

 

 

8.160.

y′′−5y′+ 4 y = 0 .

8.161.

y′′−3y = 0.

 

 

 

 

8.162.

y′′+ 2 y = 0 .

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

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

 

 

 

 

 

 

8.163. λ1,2 = 2 ±i .

 

 

 

 

8.164. λ1 = −1, λ2 = 2 .

Решить задачи Коши:

 

 

 

 

 

 

 

 

 

8.165.

y

′′

+ 2 y

+ y = 0 ; y(0) =

 

 

 

 

 

 

 

 

′′

2 , y (0) = −5 .

 

 

 

 

8.166.

y

+ 2 y

= 0 ;

y(0) = −1,

 

 

 

 

 

 

 

 

′′

 

 

y (0) = −2 .

 

 

 

 

 

8.167.

y

+ 4 y = 0 ;

 

 

 

= −2 .

 

 

 

 

 

 

y(0) = 0 , y (0)

 

 

 

 

 

8.168.

y

′′

4 y

+ 4 y = 0 ; y(0) = 0 ,

 

 

 

 

 

 

 

y (0) =1.

 

 

 

 

8.169.

Найти интегральную кривую дифференциального уравнения

y′′− y = 0 , касающуюся в точке (0, 0) прямой y = x .

 

π

 

8.170.

Решить задачу Коши: y

′′

+ y = 0 ; y(

π

2) =

2) = −2 .

 

 

3, y (

 

 

 

 

 

Задачи повышенной сложности

 

 

 

8.171. Решить уравнение 2 y′′′− y′′−2 y′+6 y = 0 ,

если известно, что

функция 2ex sin x является его решением.

 

 

 

 

 

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

 

 

 

 

 

8.172.

y′′′−2 y′′−3y′+10 y = 0 .

 

 

 

8.173.

y(4) + 2 y′′+ y = 0.

27

8.174. Показать, что общее решение уравнения x + ω2 x = 0 может быть представлено в виде x = Asin(ωt ) или x = Acos(ωt + ϕ) , где A 0

и ϕ – произвольные постоянные.

Решить задачи Коши:

 

 

 

 

8.175.

x +3x +3x + x = 0 ;

x(0) =1, x(0) = x(0) = −1.

8.176.

y

′′′

y

 

= −1,

′′

=1.

 

 

= 0 ; y(0) = 3 , y (0)

y (0)

8.177. Узкая длинная трубка вращается с постоянной угловой скоростью ω вокруг перпендикулярной к ней вертикальной оси. Шарик, находящийся внутри трубки, скользит по ней без трения. Найти закон движения шарика относительно трубки (расстояние шарика от оси вращения как функцию времени), если в начальный момент шарик находился на оси вращения и имел начальную скорость v0 .

8.178. Материальная точка массой m движется прямолинейно под действием силы притяжения к неподвижному центру, пропорциональной расстоянию от точки до центра с коэффициентом пропорциональности k. Найти закон движения точки, если в начальный момент расстояние от точки до центра равно a, а ее скорость равна нулю.

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

Согласно теореме о структуре общего решения н.л.д.у. (п. 3) общее решение н.л.д.у. с постоянными коэффициентами

y(n)

+ a1 y(n1) + … + an1 y′+ an y = b(x)

(9)

записывается в виде y = y0

+ y , где y0 – общее решение соответствующего однородно-

го уравнения (5), а y – любое частное решение уравнения (9). Для отыскания y в об-

щем случае можно воспользоваться методом вариации произвольных постоянных (за-

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

Ограничимся рассмотрением случая n = 2 . Пусть y1 и

y2 – фундаментальная

система решений уравнения (7). Будем искать частное решение уравнения

y′′+ a1 y′+ a2 y = b( x)

(10)

в виде y = C1(x) y1 +C2 (x) y2 (этот вид получается, если в выражении общего решения

соответствующего однородного уравнения произвольные постоянные считать функциями от x – отсюда название метода). Для того чтобы функция y была решением

уравнения (10), достаточно, чтобы функции C1(x) и C2 (x) удовлетворяли системе уравнений

{C1y1 +C2y2 = 0, C1y1′ +C2y2′ = b(x).

Решив эту линейную систему относительно C1и C2, далее путем интегрирования найдем C1(x) и C2 (x) .

28

Пример 8. Найти общее решение уравнения y′′− 2 y′+ y = x2e+x 1.

3Найдем сначала общее решение y0 соответствующего однородного уравнения y′′− 2 y′+ y = 0 .

Характеристическое уравнение

λ2 2λ +1 = 0 имеет двукратный корень λ

=1, по-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2

этому y = ex

и y

2

= xex

– фундаментальная система решений, откуда

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

0

 

= C ex

+C

2

xex .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

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

y

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

ции произвольных постоянных. Полагаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = C (x)ex

+C (x)xex .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Найдя y1′ = ex

и y2′ = (1 + x)ex , составим систему уравнений относительно C1

и C2:

 

 

 

 

 

 

 

x

+C

xe

x

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

C e

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ex

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1 + x)e = 2

 

 

,

 

 

 

 

 

 

 

C1e

 

+C2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решив которую,

получим

C1′ = −

 

 

x

 

 

 

 

,

C2

=

 

 

1

 

 

.

Путем интегрирования этих ра-

x2

 

+1

x2 +1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

венств находим C (x) = −ln

x2 +1,

C

2

(x) = arctg x (для нахождения частного решения

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

достаточно взять по одной из первообразных функций C1и C2). Следовательно,

 

 

 

 

 

y = −ex ln

 

 

 

x2 +1 + xex arctg x .

 

Складывая в заключение y0

и y , получаем общее решение данного уравнения

 

 

 

 

y = ex (C +C

2

x ln

 

x2 +1 + x arctg x) .4

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

8.179. λ2 +λ −1 = 0 ; y = ex (x 3) .

8.180. λ2 −λ −2 = 0;

y = 2x 1.

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

 

 

 

 

 

 

 

 

′′

 

2

 

 

′′

 

 

 

ex

 

8.181. y

+ y +ctg

x = 0 .

8.182. y

2 y

+ y =

4 x2 .

 

 

 

 

В случае функции b(x) следующих двух специальных видов для нахождения част-

ного решения уравнения (10) (уравнения (9) в общем случае) целесообразно использовать другой метод – метод неопределенных коэффициентов, который не требует операции интегрирования. Обозначим через λ1, λ2 – корнихарактеристическогоуравнения(8).

29

I. Пусть b( x) = eax P(x) , где a R, а P(x) – многочлен (в частности, b(x) = Aeax или b(x) = P(x) ). Тогда:

а) если a ≠ λ1 и a ≠ λ2 (т.е. число a не является корнем характеристического уравнения), то частное решение ищется в виде y = eaxQ(x) (т.е. в виде, аналогичном функции b( x) ),

б) если λ ≠ λ

2

и

a = λ

или a = λ

2

(a – простой корень), то y = xeaxQ(x) ,

1

 

 

1

 

 

в) если a = λ = λ

2

(a – двукратный корень), то y = x2eaxQ( x) ,

1

 

 

 

 

 

 

где Q(x) – многочлен той же степени, что и P(x) , записанный в общем виде, т.е. с не-

определенными (буквенными) коэффициентами (напомним, в частности, что общий вид многочлена нулевой степени – Q(x) = A , A 0 ). Коэффициенты многочлена Q(x)

находят после подстановки функции y в данное уравнение (10) из условия обращения этого уравнения в тождество.

II.

Пусть

b( x) = eax (P1(x) cos bx + P2 ( x) sin bx) , где

a,b R, а

P1(x) и P2 (x)

многочлены (в частности, b(x) = Acos bx + Bsin bx , A, B R). Тогда:

 

 

а) если комплексное число a + bi не является корнем характеристического урав-

нения, то частное решение ищется в виде y = eax (Q (x) cos bx +Q (x) sin bx) ,

 

 

 

1

2

 

 

б)

если

a + bi является корнем характеристического

уравнения,

то

y = xeax (Q (x) cos bx +Q ( x) sin bx) .

 

 

 

 

1

2

 

 

 

Здесь Q1( x) и Q2 ( x) – многочлены одинаковой степени, равной наибольшей из степеней многочленов P1(x) и P2 (x) , записанные в общем виде, т.е. с неопределенными (буквенными) коэффициентами. Коэффициенты многочленов Q1(x) и Q2 ( x) находят после подстановки функции y в данное уравнение (10) из условия обращения этого урав-

нения в тождество.

Пример 9. Найти общее решение уравнения y′′+ y′− 2 y = ex (6x + 2) .

3Найдем общее решение y0 соответствующего однородного уравнения y′′+ y′− 2 y = 0 .

Характеристическое уравнение λ2 + λ −2 = 0

имеет корни

λ = −2

, λ

2

=1

, поэтому

 

 

1

 

 

 

e2 x и ex – фундаментальная система решений, откуда y0 = C1e2 x +C2ex .

Частное решение y исходного уравнения найдем методом неопределенных коэффициентов. Так как правая часть уравнения – вида I ( a =1 , P(x) = 6x + 2 – многочлен 1-й степени) и число a совпадает с корнем λ2 (случай б)), то частное решение будем искать в виде

y = xex ( Ax + B) = ex ( Ax2 + Bx) .

Для отыскания коэффициентов A и B найдем производные y, y′′ и подставим y , yи y′′ в исходное уравнение. Получим (после сокращения на ex ) тождество

( Ax2 + Bx + 4 Ax + 2B + 2 A) + ( Ax2 + Bx + 2 Ax + B) 2( Ax2 + Bx) 6x + 2 ,

или

6Ax +(2 A + 3B) 6x + 2 .

30

Приравнивая коэффициенты при одинаковых степенях x в обеих частях этого тождества, получим систему уравнений для определения неизвестных A и B:

{62AA +=36,B = 2,

откуда A =1, B = 0 . Итак, y = x2ex , и, следовательно, общее решение данного уравне-

ния имеет вид

y = y0 + y = C1e2 x + ex (C2 + x2 ) .4

Пример 10. Найти решение уравнения

y′′− 4 y′+ 4 y = 8sin 2x ,

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

3Характеристическое уравнение соответствующего однородного уравнения λ2 4λ + 4 = 0 имеет двукратный корень λ1,2 = 2 , поэтому общее решение однородного

уравнения есть y0 = e2 x (C1 +C2 x) .

Частное решение данного уравнения будем искать в виде y = Acos 2x + B sin 2x ,

т. к. правая часть вида II ( a = 0 , b = 2 , P1( x) = 0 , P2 (x) = 8 – многочлен 0-й степени) и число a + bi = 2i не является корнем характеристического уравнения (случай а)). Найдя y, y′′ и подставив y , y, y′′ в исходное уравнение, получим

(4 Acos 2x 4B sin 2x) 4(2 Asin 2x + 2B cos 2x) + 4( Acos 2x + Bsin 2x) 8sin 2x ,

или

8Asin 2x 8B cos 2x 8sin 2x ,

 

 

откуда 8A =1, 8B = 0 , и, следовательно, A =1,

B = 0 и y = cos 2x .

 

 

Общее решение данного уравнения будет

y = y

0

+ y = e2 x (C +C

2

x) + cos 2x .

 

 

1

 

Для нахождения C1 и C2 продифференцируем общее решение:

 

 

y′ = e2 x (2C1 +C2 + 2C2 x) 2sin 2x ,

 

 

и воспользуемся начальными условиями: 1 = C1 +1 ,

1 = 2C1 +C2 .

Отсюда C1 = 0 ,

C2 = −1 . Итак, искомым частным решением является функция y = cos 2x xe2 x .4 Пример 11. Найти общее решение уравнения y′′+ y = 2(cos x 5e3x ) . 3Характеристическое уравнение λ2 +1 = 0 имеет корни λ1,2 = ±i . Следователь-

но, общее решение соответствующего однородного уравнения – y0 = C1 cos x +C2 sin x . Правая часть данного уравнения является суммой двух функций – b1(x) = 2 cos x и b2 ( x) = −10e3x . Для нахождения частного решения воспользуемся принципом супер-

позиции решений (п. 3): если y1 и y2 – частные решения уравнений y′′+ y = b1( x) и y′′+ y = b2 ( x) соответственно, то y = y1 + y2 – частное решение данного уравнения. Так

как

y1 следует искать в виде y1 = x( Acos x + B sin x) ( a = 0 ,

b =1 , число a + bi = i явля-

ется

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

P1(x) = 2

 

– многочлен 0-й степени,

P (x) = 0 , т.е. имеет место случай II, б)), а y

2

– в виде y

2

= Ce3x ( a = −3 не является

2

 

 

 

 

корнем характеристического уравнения, P(x) = −10 – многочлен 0-й степени, т.е. случай I, а)), то частное решение будем искать в виде

y = x( Acos x + B sin x) +Ce3x .

Найдя y′′ и подставив y и y′′ в исходное уравнение, получим тождество

31

2 Asin x+2Bcosx+x(AcosxBsin x)+9Ce3x +x( Acosx+Bsin x)+Ce3x 2cosx10e3x

или

2B cos x 2 Asin x +10Ce3x 2 cos x 10e3x .

Приравнивая коэффициенты при одинаковых функциях в обеих частях этого тождества, получим систему уравнений для определения неизвестных A, B и C:

2B = 2 , 2 A = 0 , 10C = −10 ,

откуда A = 0 , B =1, C = −1. Итак, y = x sin x e3x , и, следовательно, общее решение y = C1 cos x +C2 sin x + x sin x e3x .4

Написать, в каком виде следует искать частное решение данного н.л.д.у. методом неопределенных коэффициентов:

8.183.

y′′−8 y′+16 y = (1x)e4 x

8.184.

y′′+5y′+4 y = 3ex .

8.185.

y′′−7 y′= (x 1)2 .

 

 

8.186.

y′′+ 2 y′+5y = ex (x cos 2x 3sin 2x) .

y′′+4 y = −5cos 2x .

8.187.

y′′+9 y = 3x sin x .

8.188.

8.189.

y′′−4 y′ = 2cos2 2x .

8.190.

y′′−2 y′+ y = x +sin x .

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

8.191.

y′′−4 y′+3y = e3x .

8.192.

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

8.193.

2 y′′− y′= 3x2 12x .

8.194.

8.195.

y′′+4 y′+5y = 4e3x .

8.196.

8.197.

y′′−2 y′+ y = 2ex .

8.198.

8.199.

y′′−5y′+6 y = 78sin 3x .

8.200.

8.201.

y′′−6 y′+10 y =18sin x 12cos x .

8.202.

y′′+ y = −9sin 2x .

8.203.

8.204.

y′′+4 y =8 6cos x .

 

y′′−2 y′+ 2 y = 5sin x .

y′′+ 4 y′+ 4 y = 9ex . y′′− y = ex .

y′′+9 y +3x 9 = 0 . y′′−2 y′+ y = 2sin x .

y′′+3y′− 4 y = 5e4x36xex .

Решить задачи Коши:

 

 

 

 

 

 

8.205.

y

′′

+ y

2 y

= 3

6x ;

 

 

 

′′

 

y(0) = y (0) =1.

 

8.206.

y

+4 y =13e

 

 

;

y(0)

= 2 ,

 

 

 

3x

y (0) = 3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.207.

y

′′

4 y = 4e

 

 

;

y(0) = −2 ,

 

 

 

2 x

y (0) = −3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.208.

y

′′

y

= 2 2x ;

y(0)

=1,

 

 

′′

 

y (0) = 0 .

 

8.209.

y

+ y = 3cos 2x ;

 

= 0 .

 

 

y(0) = y (0)

 

8.210.

y

′′

+2 y

+ y

= 2(cos x sin x) ;

= 2 .

 

 

y(0) =1, y (0)

Задачи повышенной сложности

Найти линейное д. у., общее решение которого имеет вид:

8.211. y = C1 +e2 x (C2 cos 2x +C3 sin 2x) +sin 2x .

32

8.212.

y = C +C

2

x +3x2

+(C +C

4

x)ex +2cos x .

 

 

 

 

 

 

 

1

 

 

 

 

 

3

 

 

 

 

 

 

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

 

 

 

 

8.213.

y′′′+ y′′− y′− y = 4(1 ex ) .

 

 

 

 

 

8.214.

y(4) + y′′= x2 + x .

 

 

 

 

 

 

 

 

Решить задачи Коши:

 

 

 

 

 

 

 

 

 

8.215.

y

 

 

y =8e

 

;

 

y(0) = −1,

 

′′

= 5,

′′′

= 6 .

(4)

x

 

y (0) =

2 , y (0)

y (0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.216.

y

′′

+4 y =8cos

 

x ;

 

 

 

= −6 .

 

 

 

 

2

y(0) =1, y (0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.217. Тяжелый однородный трос переброшен через гладкий гвоздь так, что с одной стороны свисает часть его длиной 8 м, а с другой стороны

– часть длиной 10 м. За какое время трос соскользнет с гвоздя? Принять g = 9,81м/c2.

8.218. Электрическая цепь состоит из последовательно соединенных источника тока с ЭДС E(t) = E0 sin ωt , катушки индуктивности индуктив-

ностью L и конденсатора емкостью C, причем ω=

1

(случай резонанса).

LC

 

 

Найти ток I в цепи как функцию времени t, если I t=0 = I t=0 = 0 .

§ 3. Системы дифференциальных уравнений

1. Основные понятия. Система дифференциальных уравнений вида

 

y1′ = f1(x, y1, y2 , , yn ),

 

y2′ = f2 (x, y1, y2

, , yn ),

(1)

 

 

 

, , yn ),

 

yn′ = fn (x, y1, y2

 

связывающая независимую переменную x, n искомых функций y1( x), y2 ( x), , yn ( x) ,

их производные 1-го порядка, и разрешенная относительно этих производных, называ-

ется нормальной системой.

Число n, равное сумме порядков уравнений системы (совпадающее для системы

(1) с числом искомых функций и числом уравнений), называется порядком системы. Решением (частным решением) системы (1) на интервале I называется всякая

совокупность (система) n функций y1 = ϕ1( x), y2 = ϕ2 ( x), , yn = ϕn ( x) , при подстанов-

ке которых в систему вместе с их производными каждое уравнение системы обращается в тождество относительно x I .

Равенства вида

y1( x0 ) = y10 , y2 ( x0 ) = y20 , , yn (x0 ) = yn0

(2)

(используется также запись y1 |x=x0 = y10 , y2 |x=x0 = y20 , , yn |x=x0 = yn0 ), где x0 , y10 , y20 , , yn0

– заданные числа (начальные значения), называются начальными условиями для систе-

мы (1). Задача отыскания решения системы (1), удовлетворяющего заданным начальным условиям (2), называется задачей Коши для этой системы.

33

Теорема (существования и единственности решения задачи Коши для системы

(1)). Если функции f1, f2 , , fn из нормальной системы (1) непрерывны вместе с част-

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

fi

по переменным y1, y2 , , yn , i, j =1, 2, , n , в некоторой об-

 

 

y j

ласти D Rn+1 , то для любой точки ( x0 , y10 , y20 , , yn0 ) D задача Коши для системы

(1) с начальными условиями (2) имеет и притом единственное решение.

Совокупность n функций

y1 = ϕ1( x,C1,C2 , ,Cn ), y2 = ϕ2 ( x,C1,C2 , ,Cn ), , yn = ϕn (x,C1,C2 , ,Cn ) (3)

называется общим решением системы (1), если:

1) при любых допустимых значениях параметров C1,C2 , , Cn она является ре-

шением этой системы; 2) любое решение системы (кроме, быть может, отдельных решений) может

быть получено из (3) при некоторых значениях параметров (т.е. для любой точки ( x0 , y10 , y20 , , yn0 ) D (см. теорему) найдутся такие значения параметров C1,C2 , , Cn , при которых функции (3) будут удовлетворять начальным условиям (2)).

Пример 1. Показать, что совокупность функций

y = C1ex +C2e3x , z = 2C1ex 2C2e3x

является общим решением системы

{zy==zy4zy, .

Найти частное решение этой системы, удовлетворяющее начальным условиям y(0) = −1, z(0) = 2 .

3Проверим выполнение условий из определения общего решения.

1)Число произвольных постоянных в функциях данной совокупности равно 2, что совпадает с порядком системы.

2)При любых значениях C1 и C2 функции y и z образуют решение системы:

y = C1ex +C2e3x ,

 

 

y′ = −C1ex + 3C2e3x ,

{z = 2C ex

2C

 

e3x

{z′ = −2C ex 6C

e3x ;

 

1

 

 

2

 

 

 

 

 

1

2

 

подставляя в данную систему y, z,

yи z, получим тождества:

C ex + 3C

e3x

= C ex +C

e3x (2C ex 2C e3x ),

1

2

 

 

 

1

 

2

2C

 

1

2

e3x )

{2C ex

6C

e3x = 2C ex

e3x 4(C ex +C

1

 

2

 

 

 

 

1

 

2

 

1

2

 

C ex

+ 3C e3x

= −C ex +3C

e3x ,

x R.

 

1

 

2

e3x

1

 

2

 

{2C ex

6C

= −2C ex 6C e3x

 

1

 

 

2

 

 

 

 

1

 

2

 

 

3) Каковы бы ни были начальные значения x0 , y0 , z0 , система уравнений

{y0 = C1ex0 +C2e3x0 , z0 = 2C1ex0 2C2e3x0

имеет, притом единственное решение относительно C1 и C2 , т. к. ее определитель

=

e

 

x

e

 

3x = −4e2 x0

0 x0 R.

 

 

x0

 

3x0

 

 

 

2e

0

2e

0

 

 

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

34

Для нахождения частного решения после подстановки начальных значений x = 0 , y = −1 и z = 2 в общее решение получим систему

{1 = C1 +C2 , 2 = 2C1 2C2 ,

откуда C = 0 ,

C

2

= −1 и, следовательно, искомое

частное

решение

есть

 

y = −e3x ,

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = 2e3x .4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.219. Показать, что функции y = x +

1

, z

= ex

x

1

 

образуют

x2

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

3

 

 

 

 

 

 

 

 

 

2 y

 

 

 

z′= y + z +

2 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

решение системы y′=1 x

,

 

x

1.

 

 

 

 

 

1

 

 

 

 

 

 

 

8.220.

Показать, что функции

y = C eC2x ,

z = −

 

 

eC2 x образуют

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

C1C2

 

 

 

 

решение системы y′ = −1z ,

 

 

z′ = 1 y

 

 

 

 

 

 

 

 

 

 

 

 

 

при любых C1 ,

C2 ,

C1C2 0 . Найти

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

y(0) =1,

z(0) = 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Сведение уравнения n-го порядка к нормальной системе. Метод исклю-

чения. Дифференциальное уравнение n-го порядка

y

 

 

 

 

 

 

 

 

 

) путем

(n)

= f ( x, y, y , , y

(n1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

введения функций

 

 

 

 

 

 

 

 

 

 

можно свести к нормальной системе n-

y1 = y, y2 = y , , yn = y

(n1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

го порядка:

 

 

 

y

= y

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

= y3

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yn′ = f (x, y1, y2 , , yn ).

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2. Привести к нормальной системе дифференциальное уравнение

 

 

 

 

 

 

 

 

 

 

 

y′′+ xy′− y2 = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

. Тогда z

= y

′′

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

Положим

z = y

 

 

y′ = z,

 

{z′ = y2 xz.

4

Решение нормальной системы n-го порядка методом исключения основано на сведении системы к дифференциальному уравнению n-го порядка относительно одной из неизвестных функций (или, в отдельных случаях, к нескольким дифференциальным уравнениям с одной неизвестной функцией в каждом). Например, при n = 2 , для решения системы

{y1′ = f1(x, y1, y2 ), y2′ = f2 (x, y1, y2 ),

выражаем из первого уравнения y2 = g(x, y1, y1) (или y1 = h(x, y2 , y2) из второго), после чего, находя полную производную dydx2 и подставляя эти выражения во второе уравне-

35

ние (соответственно

dy1

и подставляя в первое), получаем уравнение 2-го порядка от-

 

 

dx

y1 ( y2

 

носительно функции

соответственно). Найдя его общее решение

y1 = ϕ1( x,C1,C2 )

и подставив эту

функцию в выражение y2 = g(x, y1, y1) , получим

y2 = ϕ2 (x,C1,C2 )

и,

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

общее решение системы как совокупность этих

двух функций.

 

 

 

 

Пример 3. Решить методом исключения систему

{zy==zy4zy, .

3Выразим из первого уравнения z: z = y y. Отсюда z′ = y′− y′′. Подставив

эти выражения во второе уравнение системы, получим уравнение 2-го порядка относи-

тельно y:

y′− y′′ = y y′−4 y или y′′− 2 y′−3y = 0 .

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

y = C1ex +C2e3x .

Отсюда y′ = −C1ex + 3C2e3x и

z = y y′ = C1ex +C2e3x (C1ex + 3C2e3x ) = 2C1ex 2C2e3x .

Таким образом,

{y = C1ex +C2e3x , z = 2C1ex 2C2e3x

общее решение данной системы.4

Привести к нормальным системам дифференциальные уравнения:

8.221.

y′′′− xyy′′+ y3 = 0.

8.222. y(4) y2 = 0 .

Решить системы методом исключения:

 

8.223.

y′ = 1z ,

8.224.

y′ = 2 y + z,

 

{z′ = 3y +4z.

 

z′ = 1 y .

 

3. Собственные значения и собственные векторы матрицы. Пусть A – квад-

ратная матрица n-го порядка, X – матрица-столбец размера n ×1, называемая далее (n- мерным) вектором, 0 – нулевой вектор. Число λ (действительное или комплексное) на-

зывается собственным значением, а вектор X 0 (возможно и с комплексными компо-

нентами) – соответствующим собственным вектором матрицы A, если AX = λX .

Алгебраическое уравнение n-й степени det( A −λE) = 0 относительно неизвестной λ, где E – единичная матрица n-го порядка, а det( ) означает определитель, называ-

ется характеристическим уравнением матрицы A. Все собственные значения матрицы A совпадают с корнями ее характеристического уравнения; отсюда следует, что для любой матрицы A существует от одного до n собственных значений. Собственные векторы, соответствующие различным собственным значениям, линейно независимы. Для нахождения всех собственных векторов, соответствующих собственному значению λ, необходимо найти все ненулевые решения однородной системы линейных уравнений

36

( A −λE) X = 0 , где X – вектор из неизвестных; компоненты каждого такого решения,

записанные в столбец, образуют собственный вектор.

Пример 4. Найти собственные значения и какие-нибудь соответствующие собственные векторы матрицы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 0

 

 

 

 

λ

0

 

 

3Составим характеристическоеуравнение матрицы A

( E =

 

 

 

λE

=

 

 

 

0 1

,

 

0

 

):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 −λ

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

det( A −λE) =

 

 

= λ2 2λ −3

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

1 −λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Его корни, т.е. собственные значения матрицы A,

λ = −1 ,

λ

2

 

= 3. Пусть

X

1

=

α1

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α2

 

 

 

– соответствующие собственные векторы. Для нахождения X1 составим од-

X2 =

1

 

 

β2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

нородную систему линейных уравнений

 

 

 

 

 

 

 

 

2 1 α1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( A −λ E) X

 

= ( A + E) X

 

=

 

=

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 α2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{

 

2α − α

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4α1 + 2α2

 

 

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Одно из ненулевых решений этой системы α1 =1 ,

α2 = 2 . Таким образом, один из соб-

ственных векторов, соответствующих собственному значению λ1 , имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X1

= 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Заметим, что общее решение системы α1 = C1 ,

α2 = 2C1 , где C1

– произвольная посто-

янная,

так что равенство X

1

= C

 

1

 

, C 0 ,

определяет все собственные векторы, со-

 

 

 

 

 

 

 

1

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ответствующие собственному значению λ1 .)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично находим X2 :

 

 

 

 

 

 

 

2 1 β

 

 

 

 

 

2β −β

 

= 0,

 

 

 

 

 

 

 

 

 

( A −λ2 E) X2 = ( A 3E) X2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

4

2 β1

= 0 {4β

12β2

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

 

 

 

1

 

 

2

 

1

 

 

 

 

 

β =1,

β

2

= −2 ,

и, следовательно,

 

X

2

=

 

 

.

 

Ответ:

λ = −1 ,

X

1

=

 

;

λ

2

= 3,

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X2 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найти собственные значения и какие-нибудь соответствующие соб-

ственные векторы матриц:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.225.

 

8 2

 

8.226.

 

4

 

5

.

 

 

 

8.227.

2 2

 

 

 

 

 

 

 

 

.

 

 

 

1

 

2

 

 

 

 

 

5 0

.

 

 

 

 

 

 

 

 

 

 

 

 

9 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 3

 

 

 

 

 

 

 

2 1 2

 

 

 

 

 

 

 

4 5 2

 

 

 

 

 

 

8.228.

 

8.229.

 

 

5

3

 

 

3

.

8.230.

 

5 7 3

.

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 1

 

 

 

 

 

 

 

 

1 0

 

2

 

 

 

 

 

 

 

 

 

6

9 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

37

4. Решение однородных линейных систем с постоянными коэффициентами методами линейной алгебры. Однородная линейная система дифференциальных уравнений (о.л.с.д.у.) n-го порядка имеет вид

 

 

x1

= a11(t)x1

+ a12 (t)x2

+ … + a1n (t)xn ,

 

 

 

 

 

 

 

 

x2

= a21

(t)x1

+ a22

(t)x2

+ … + a2n (t)xn ,

 

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= an1(t)x1

+ an2 (t)x2 + … + ann (t)xn ,

 

 

 

 

 

 

 

 

xn

 

 

 

 

 

 

где t – независимая переменная,

x1, x2 , , xn – искомые функции от t,

xi = dxi . Если

ввести матрицы

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a11(t) a12 (t) a1n (t)

 

 

x1

 

 

 

x1

 

 

a21

(t) a22

(t) a2n (t)

 

 

x2

 

и

X

x2

 

,

A(t) =

 

 

 

(матрица системы), X =

 

=

 

 

(t) an2

 

 

 

 

 

 

 

 

 

 

 

 

 

an1

(t) ann (t)

 

 

xn

 

 

 

xn

 

 

то равенство X = A(t) X есть матричная запись системы (4).

Произвольное решение системы (4) можно записать в виде матрицы-столбца (вектор-функции) X (t) . Всякая совокупность из n линейно независимых решений

Xk (t) , k =1, 2, , n , называется фундаментальной системой решений системы (4).

Справедлива

Теорема (о структуре общего решения о.л.с.д.у.). Если X1, X2 , , Xn – какая-

нибудьфундаментальная системарешенийсистемы(4), то ее общее решение имеетвид:

X = C1 X1 +C2 X2 + … +Cn Xn ,

где C1,C2 , ,Cn – произвольные постоянные.

В частном случае систем с постоянными коэффициентами, когда aij (t) = const , т.е. матрица A(t) = A не зависит от t, для отыскания фундаментальной системы реше-

ний может быть использован аппарат собственных значений и собственных векторов. Характеристическое уравнение det( A −λE) = 0 матрицы A о.л.с.д.у. с постоян-

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

 

 

 

X = AX

(5)

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

 

 

 

Теорема (о характеристическом уравнении). Вектор-функция

X (t) =Yeλt

 

,

0

где Y – числовой n-мерный вектор, тогда и только тогда является решением системы

(5), когда λ есть собственное значение матрицы A (корень характеристического уравнения системы), а Y – соответствующий собственный вектор.

Рассмотрим подробно случай n = 2 , т.е. систему

{x = a11x + a12 y, y = a21x + a22 y.

Тогда характеристическое уравнение

a11 − λ

a12

 

= 0

a21

a22 − λ

 

является квадратным и, следовательно, возможны следующие 3 случая:

38

а) характеристическое уравнение имеет 2 различных действительных корня λ1 и λ2 . Тогда, если Y1 и Y2 – какие-нибудь соответствующие собственные векторы, век- тор-функции X1 =Y1eλ1t и X2 =Y2eλ2t образуют фундаментальную систему решений и, следовательно, согласно теореме о структуре общего решения оно имеет вид

X = C1 X1 +C2 X2 ;

б) корни характеристического уравнения λ = α+iβ и λ = α−iβ комплексно со-

пряжены. В этом случае, найдя для собственного значения λ какой-нибудь соответствующий собственный вектор Y (с комплексными компонентами), в качестве фундаментальной системы решений можно взять действительную и мнимую части комплексного решения Yeλt , т.е. X1 = ReYeλt и X2 = ImYeλt ;

в) характеристическое уравнение имеет двукратный действительный корень λ. Тогда одно из решений фундаментальной системы X1 =Yeλt , где Y – один из соответ-

ствующих собственных векторов. Если существует собственный вектор Z , линейно не зависимый от Y – этот случай имеет место, когда A = aE , т.е. система распадается на два независимых уравнения (тогда все ненулевые векторы – собственные), – то другое

решение имеет вид X2 = Zeλt . В противном случае X2 = (Y +tY )eλt

, где Y

– любое ре-

шение системы ( A −λE)Y =Y (другой способ нахождения общего решения системы (5)

в этом случае – положить одну из неизвестных функций равной

eλt (C +C

t) (либо

 

1

2

 

C eλt , если система содержит уравнение с одной неизвестной функцией), а оставшуюся

1

 

 

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

Пример 5. Найти общее решение системы

{xy == xyy4,x.

3Матрица системы

1

1

 

A =

4

.

 

1

Ее собственные значения и соответствующие собственные векторы

λ1 = −1 ,

 

λ2 = 3,

Y =

1

,

Y

=

1

(см. пример

4).

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

вектор-функции

X

1

=

1

et и

1

 

2

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X2 =

 

1

 

– фундаментальная система решений, и общее решение системы имеет вид

 

 

e3t

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X = x

= C X

1

+C

2

X

2

= C

1 et

+C

2

 

1 e3t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

т.е.

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

x = C et +C

e3t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

 

e3t . 4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{y = 2C et

2C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 6. Найти общее решение системы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{xy == x2+xy+, 3y.

 

 

 

 

 

 

 

 

 

 

 

 

3Характеристическое уравнение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 −λ 1

 

 

= λ2

4λ +5 = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

3 −λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

39

имеет комплексно сопряженные корни λ1,2

= 2 ±i . Для нахождения какого-нибудь соб-

ственного вектора Y =

 

α

 

, соответствующего собственному значению λ = 2 +i , име-

 

1

 

ем систему

 

 

 

α2

 

 

 

 

 

 

 

 

 

 

 

 

1 i

1 α

 

 

 

(1 i)α + α = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( A −λE)Y =

2

1 i

α1

 

= 0 или {2α + (11i)α2 = 0.

 

 

 

 

 

 

 

 

2

1

2

Пусть α =1

. Тогда α

 

=1 +i ,

т.е. Y =

1

 

. Итак, данная система имеет комплексное

1

 

2

 

 

 

 

 

1 +i

 

 

 

 

решение вида

 

 

 

 

1

 

 

 

 

 

 

cos t +i sin t

 

 

X =Yeλt =

 

e(2+i)t =

 

 

 

e2t

 

 

 

 

 

1 +i

(cos t sin t) +i(cos t +sin t)

(использована формула Эйлера eit = cos t +i sin t ). В качестве фундаментальной системы решений X1 и X2 возьмем Re X и Im X соответственно:

 

X1

 

cos t

 

 

 

 

=

 

sin t

 

 

 

=

 

e2t , X2

 

 

e2t .

 

 

 

cos t sin t

 

 

 

 

cos t +sin t

 

Итак, общее решение системы

 

 

 

 

 

 

 

 

 

 

 

 

 

x

= C

cos t

 

e2t +C

2

 

sin t

e2t

,

 

 

1

 

 

 

 

 

 

 

 

 

 

или

y

 

cos t

sin t

 

 

 

cos t +sin t

 

x = e2t (C cos t +C

 

sin t),

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

1

+C

 

 

 

(C

C ) sin t). 4

 

 

{y = e2t ((C

2

) cos t +

 

 

 

 

1

 

 

 

 

 

2

1

 

 

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

{x = 2x y, y = 4x + 6 y,

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

2 − λ −−1λ = λ2 8λ +16 = 0 4 6

имеет корень λ = 4 кратности 2. Находим соответствующий собственный вектор:

 

 

 

 

 

 

2 1 α

 

 

 

 

2α − α

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 2

 

α1

 

= 0 {4α

1+ 2α2

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

1

 

 

2

 

 

 

 

 

 

откуда

α =1

,

α

2

= −2

, Y =

1

 

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

X

1

=

1

e4t . Для определения X

2

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

найдем любое решение Y

системы 2 1

β1

=

 

1

,

например, β = 0 ,

β

2

= −1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 2

β2

 

 

 

 

 

 

т.е. Y = 01 . Отсюда

X2 = (Y +tY )eλt = ( 01 + 12 t)e4t = 1 t 2t e4t .

Итак, общее решение данной системы

{x = e4t (C1 +C2t),

y = e4t (2C1 C2 2C2t).

40

Для нахождения частного решения константы C1 и C2 определяем из системы

{10==C21C, 1 C2 ,

полученной в результате подстановки начальных значений t = 0 , x = 0 и y =1 в общее решение, откуда C1 = 0 , C2 = −1 , и, следовательно, искомое частное решение есть

 

 

 

x = −te4t ,

 

 

 

 

 

 

{y = e4t (1 + 2t). 4

 

 

 

Решить характеристические уравнения систем:

 

 

x = −2x y,

 

 

 

x = x 2 y,

8.231. {y = −2 y + x.

 

8.232. {y = 2x 4 y.

Решить системы дифференциальных уравнений:

 

 

x = x +3y,

 

 

x = y,

 

 

x = 4x 5y,

8.233. {y = −x +5y.

8.234.

{y = −2x +3y.

8.235.

{y = 2x 3y.

 

x = 2x 3y,

 

x = 2x 5y,

 

x = 3x 2 y,

8.236. {y = x 2 y.

 

8.237.

{y = 5x 6 y.

8.238.

{y = 4x +7 y.

 

x = −x +2 y,

 

x = x 4 y,

 

 

8.239. {y = −2x 5y.

8.240.

{y = x 3y.

 

 

Решить задачи Коши:

 

 

 

 

 

8.241.

x = −2x y,

y =8x + 2 y; x(0) = 0 ,

y(0) = 2 .

 

8.242.

x = x + 2 y,

y = −x + 4 y;

x(0) = 2 ,

y(0) =1.

 

Задачи повышенной сложности

8.243. Материальная точка массой m движется под действием силы притяжения к неподвижному центру O, пропорциональной расстоянию от точки до центра с коэффициентом пропорциональности k. Движение начинается из точки A на расстоянии a от центра с начальной скоростью v0 , перпендикулярной к отрезку OA. Найти траекторию движения.

8.244. Решить систему: x = x 2 y z, y = −x + y + z, z = x z.

41

§ 4. Элементы теории устойчивости

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

x1

= f1(t, x1, x2 , , xn ),

 

x2

= f2 (t, x1, x2

, , xn ),

(1)

 

 

 

 

= fn (t, x1, x2 , , xn )

 

xn

 

и некоторые начальные условия в точке t0 . Решение ϕ1(t), ϕ2 (t), , ϕn (t) системы (1), удовлетворяющее этим начальным условиям, называется устойчивым по Ляпунову, если для любого ε > 0 существует такое δ > 0 , что для всякого решения x1(t), x2 (t), , xn (t) той же системы, начальные значения которого удовлетворяют неравенствам

| xi (t0 ) −ϕi (t0 ) |< δ, i =1, 2, , n,

(2)

при всех t > t0 выполняются неравенства

| xi (t) −ϕi (t) |< ε, i =1, 2, , n

(т.е. решение изменится сколь угодно мало, если достаточно мало изменятся начальные значения).

Если дополнительно при условии (2)

lim | xi (t) −ϕi (t) |= 0, i =1, 2, , n,

t→+∞

то это решение называется асимптотически устойчивым.

Пример 1. Исследовать на устойчивость решение дифференциального уравнения x =1 x , удовлетворяющего начальному условию x(0) =1 .

3Общее решение уравнения x(t) = Cet +1; частное решение, удовлетворяющее начальному условию, ϕ(t) =1 . Так как | x(t) −ϕ(t) |=| Cet | , то | x(t) −ϕ(t) |<| x(0) −ϕ(0) | при t > 0 . Отсюда, если положить δ = ε, то из | x(0) −ϕ(0) |< δ следует | x(t) −ϕ(t) |< ε при t > 0 . Следовательно, решение ϕ(t) =1 устойчиво. Кроме того, так как

lim | x(t) −ϕ(t) |= lim | Cet |= 0 , то это решение является также асимптотически устой-

t→+∞ t→+∞

чивым.4

Исследовать на устойчивость решения уравнений:

8.245. x = t 1; x(0) = −1 .

8.246. x = t(x 1) ; x(0) =1.

2.Классификация точек покоя однородной линейной системы 2-го порядка

спостоянными коэффициентами. Исследуем на устойчивость точку покоя (нулевое

решение x = 0 , y = 0 ) о.л.с.д.у. 2-го порядка с постоянными коэффициентами

x = a11x + a12 y,

(3)

{y = a21x + a22 y

 

(точкой покоя системы (1) называется решение, состоящее из констант). Ограничимся рассмотрением случая

det A =

a11

a12

0

(4)

 

a21

a22

 

 

42

(при этом условии точка покоя единственна).

Если x = x(t) , y = y(t) – произвольное решение системы (3) (или, вообще, нор-

мальной системы {xy== fg((tt,,xx,, yy),) ), то кривую, параметрически заданную на плоскости

xy этими уравнениями, называют фазовой траекторией. Через каждую точку плоскости xy (фазовой плоскости) проходит единственная фазовая траекториясистемы (3).

Картина фазовых траекторий системы (3) в окрестности точки покоя (0, 0) и направление движений по ним при t → +∞, определяющие устойчивость или неустойчивость точки покоя, существенно зависят от корней λ1 и λ2 характеристического урав-

нения системы. В связи с этим различают следующие случаи ( λ1λ2 0 ввиду (4)):

1) Корни λ1 , λ2 – действительны, λ1 ≠ λ2 (дискриминант характеристического уравнения D > 0 ). Тогда:

а) если λ1 < 0 и λ2 < 0 , то точка покоя асимптотически устойчива и носит на-

звание устойчивый узел (рис. 1; направление движения по всем траекториям, отличным от точки (0,0) , – к точке покоя);

б) λ1 > 0 и λ2 > 0 – точка покоя неустойчива (неустойчивый узел, рис. 1; направление движения – от точки покоя);

в) λ1λ2 < 0 – точка покоя неустойчива (седло, рис. 2).

 

Рис. 1. Узел

Рис. 2.

Седло

2)

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

λ1,2 = α±iβ комплексно сопряжены

( D < 0 ). Тогда:

 

 

а)

если α < 0 , то точка покоя асимптотически устойчива (устойчивый фокус,

рис. 3; направление движения – к точке покоя);

 

 

б)

α > 0 –точка покоя неустойчива (неустойчивый фокус,

рис. 3; направление

движения – от точки покоя);

 

 

в)

α = 0 – точка покоя устойчива (центр, рис. 4).

 

Рис. 3. Фокус

Рис. 4. Центр

43

3) характеристическое уравнение имеет двукратный действительный корень λ ( D = 0 ). Тогда:

а) если λ < 0 , то точка покоя асимптотически устойчива (устойчивый вырож-

денный узел (рис. 5) при A aE и устойчивый дикритический узел (рис. 6) в противном случае; направление движения – к точке покоя);

б) λ > 0 – точка покоя неустойчива (неустойчивый вырожденный узел (рис. 5)

при A aE и неустойчивый дикритический узел (рис. 6) в противном случае; направление движения – от точки покоя).

Рис. 5. Вырожденный узел

Рис. 6. Дикритический узел

Итак, в случае асимптотически устойчивой точки покоя (устойчивый узел, в частности, вырожденный или дикритический, устойчивый фокус) движение по всем фазовым траекториям, отличным от нее самой, при возрастании t направлено к точке покоя. Если же точка покоя неустойчива (неустойчивые узлы, седло, неустойчивый фокус), то это движение происходит в направлении от точки покоя (за исключением двух траекторий в виде открытых лучей с началом в точке (0, 0) и лежащих на одной прямой в случае седла).

Таким образом, если λ1λ2 0 , то при Re λ1 0 и Re λ2 0 точка покоя (0, 0) ус-

тойчива (асимптотически устойчива в случае строгих неравенств); в противном случае

– неустойчива.

Пример 2. Определить характер и исследовать на устойчивость точку покоя системы

{xy == x2+xy+, αy,

взависимости от параметра α ( α ≠ −2 ).

3 Характеристическое уравнение

2 − λ α

 

= λ2 + λ − (α + 2) = 0

 

1 1 − λ

 

 

имеет корни

λ

= −1

2

± 1

2

9 + 4α. Исследуя поведение корней в зависимости от па-

 

1,2

 

 

 

раметра α, получаем:

если α < −9 4 , то корни комплексные, Reλ1,2 = −12 < 0 , и, следовательно, точка покоя – устойчивый фокус;

если α = −9 4 , то λ = −12 – отрицательный корень кратности 2; точка покоя – устойчивый вырожденный узел;

если 9 4 < α < −2 , то корни действительные и отрицательные; точка покоя – ус-

тойчивый узел; если α > −2 , то корни действительные и разных знаков; точка покоя – седло.

44

Таким образом, точка покоя асимптотически устойчива при α < −2 и неустойчива при α > −2 .4

Определить характер точек покоя следующих систем:

x = x +3y,

x = 3x +4 y,

x

8.247. {y = −x 2 y.

8.248. {y = 2x + y.

8.249. {y

x = x +2 y,

x = 3x + y,

x

8.250. {y = −3x + y.

8.251. {y = 2x + y.

8.252. {y

8.253. {xy ==

3yx+xy. ,

8.254. {xy ==x y,

2 y.

 

=2x 5y,

=2x +2 y.

=6x 5y,

=2x 5y.

Задачи повышенной сложности

Чтобы более точно изобразить картину фазовых траекторий в окрестности точки покоя системы (3) в случае узла, седла или вырожденного узла, необходимо прежде всего найти фазовые траектории в виде открытых лучей с началом в точке (0, 0). Направляющими векторами каждой из прямых (двух в случае узла или седла и одной в случае вырожденного узла), на которых лежат эти траектории, являются собственные векторы матрицы системы; направления движения по фазовым траекториям определяется знаком соответствующего собственного значения – к точке покоя при λ < 0 , и от нее – в противном случае. В случае узла остальные фазовые траектории (при добавлении к ним точки (0, 0)) касаются той прямой, для которой направляющий вектор соответствует меньшему по абсолютной величине собственному значению.

Для определения направления, в котором происходит движение по фазовым траекториям вокруг фокуса или центра, достаточно построить в какой-нибудь точке (x, y)

(отличной от начала координат) вектор фазовой скорости v(x, y) = ( x, y) (направляю-

щий вектор касательной к фазовой траектории в этой точке), определяемый по формулам (3). Аналогично определяется направление движения в случае вырожденного узла. Направление закручивания фазовых траекторий в случае фокуса (по часовой стрелке или против) определяется после этого его устойчивостью или неустойчивостью, т.е. по знаку Re λ .

Пример 3. Определить характер точки покоя системы

{x = 5x + y, y = 3x + 3y

инарисовать (схематично) картину ее фазовых траекторий в окрестности начала координат с указанием на них направления движения.

3 Характеристическое уравнение

 

5 −λ 1

 

= λ2 8λ +12 = 0

 

 

 

3 3 −λ

 

 

имеет корни λ1 = 2 , λ2 = 6 . Так как λ1 ,

 

λ2 действительны, λ1 ≠ λ2 и λ1 > 0 , λ2 > 0 , то

точка покоя – неустойчивый узел. Соответствующие собственные векторы Y1 = 13 ,

45

y=–3x

y

y=x

Y2 = 1 . Следовательно, имеются четыре фазо-

 

 

 

1

 

 

 

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

 

 

 

в точке (0,0) – две из них лежат на прямой

 

 

 

y = −3x , проходящей через точку (1, 3) , две дру-

 

 

x

гих – на прямой y = x . Остальные фазовые траек-

 

 

 

тории (при добавлении к ним точки (0, 0) ) каса-

 

 

 

ются в начале координат прямой y = −3x (ее на-

 

 

 

правляющий вектор (1, 3) соответствует мень-

 

 

 

шему собственному значению λ1 ). Направление

 

 

 

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

 

Рис. 7

 

ным от точки (0,0) , ввиду неустойчивости, – от

 

 

 

начала координат (рис. 7).4

Пример 4. Решить аналогичную задачу для системы

y

v

 

x

 

1

Рис. 8

 

{x = x 2 y, y = 4x 3y.

3 Характеристическое уравнение

 

1 −λ −2

 

= λ2 + 2λ +5 =

0

 

 

 

4 3 −λ

 

 

 

имеет комплексно сопряженные корни λ1,2

= −1 ± 2i . Так как

Re λ1,2 < 0 , то точка покоя – устойчивый фокус. Следовательно,

направление движения по всем фазовым траекториям, отличным от точки (0,0) , – к началу координат. Для определения направле-

ния закручивания фазовых траекторий найдем и построим вектор фазовой скорости, например, в точке (1, 0) : при x =1, y = 0 из

данной системы имеем x =1, y = 4 , т.е. v = v(1, 0) = (1, 4) . Таким

образом, закручивание происходитпротивчасовойстрелки(рис. 8).4

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

8.255.

{xy

== xx +

23yy.,

8.256.

{xy

==x3x4+y2. y,

8.257.

{xy

==

3xx24yy. ,

8.258.

{xy

== x2xy,y.

Определить характер и исследовать на устойчивость точку покоя систем в зависимости от параметра α:

x = αx y,

 

x = −3x y,

 

8.259. {y = x + 2 y;

( α ≠ −0,5 ).

8.260. {y = −αx + y;

( α2 3 ).

46