А.К. Боярчук, Г.П. Головач - Дифференциальные уравнения в примерах и задачах (1109000), страница 68
Текст из файла (страница 68)
х х+2хз=о. ° Переходя к системе Оху, положив х = у, получаем 3 а=у, у= -2х. Отсюда почленным делением уравнений, а затем интегрированием находим Рис. 55 семейство траекторий на фазовой плоскости: у'+ х = С. Каждая траектория представляет собой леформированную окружность (рис. 55). Очевидно, точка (О, 0) — центр, м ЗО8 Гл. б.
Устойчивость и 4сазовые траепории 654. й+2хз — 2х = О. И Полагая х = у, приходим к уравнению г(у 2(х хз) с(х общее решение которого имеет вид: г = сс'сгг -*'г. Рвс. 57 Ряс. 5Ь Поскольку траектории симметричны относительно обеих координатных осей, то далее считаем х > О, у > О. Если в (1) положим х = у = О, то получим С = О. Следовательно, кривая =г гг- г проходит через начало координат (рис. 56). При С > 0 все траектории проходят выше этой кривой (рис. 57); причем в точке х = 0 производная у' = О, а при у = 0 она не ограничена (точнее, у' — -оо при у +0). Далее, при С < 0 из неравенства х(2-х»-С г г следует, что Ряс.
59 с С ~ +~~ гс г -с. Значит, при уменьшении С область существования семейства траекторий сужается и при С = — 1 траектории вырождаются в точку (1, 0) (рис, 58). Наконец, зеркально отобразив кривые, изображенные на рис. 58, относительно оси Ох, а затем полученную картину — относительно оси Оу, будем иметь полную картину траекторий (рис. 59). Очевидно, точки (1, 0) и (-1, 0) — центры, а точка (О, 0) — селло. и 655. х — 2*+ х + 1 = О. и Перейдя на фазовую плоскость, имеем г(у 2с — х — 1 х=у, у=2* — х — 1 (1) с(х у Интегрируя уравнение (1), получаем семейство траекторий: 2*ы у = — — х' — 2х+С.
1п2 (2) Далее, используя обычные методы математического анализа, строим картину интегральных кривых (2) (рис. бО). Отметим, что кривой, проходящей через точку (1, 0), соответствует значение ф С = 3 — —. 1п2 309 Замкнутым кривым, охватывающим начало координат, соответствуют значения С, определяемые неравенством 2 4 — — < С < 3 — —. 1п2 1п2 » "зо )»(о Зл 7л з з' 0 з Т 7л За Рлс. Ог Значения С для остальных кривых указаны на рис. 60. и 3 Рас. Оз 656. х+ 2созх — 1 = О.
М Из системы х = у, у = 1 — 2 соз х получаем Цу 1 — 2созх з уз = 2х — 4йпх+ С. с(х у р(х) = 2х — 4з(пх (рис. 62), а затем картину семейства ОЗ(х) = 2х — 4мпа+ С путем параллельного переноса графика функции р (рис.63). Теперь построим семейство кривых (5»(х) = 2х — 4з(па+ С > 0 (рис. 64). Далее строим семейство кривых у = ~)/ф»(х) (рис.
65) (точнее, на рис. 65 изобралсена только часть зтопз семейства в окрестности двух особых точек ( — $, 0) и (О, ~т) — центра и седла соответственно). Для получения всей картины семейства траекторий на фазовой плоскости Оху следует картину, изображенную на рис. 65, периодически (с периодом 2»г) продолжить яак влево, так и вправо относительно ее первоначального положения.
Тогда получим следующую картину (рис. 66). 1ь Решив конечную систему у = О, 1 — 2созх = О, находим особые точки: М» Я+2йа.,О), л» (-~у+ 2йзг, О), где (с б х. далее, нетрудно установить, что 1»х» — седла, а 117» — центры (рис. 61). Для построения картины семейства (1) сначала строим кривую Гл. 6. Устейчнвесп и фазовые траектории 657.
х + 2х+ 5х = О. м Сделав замену х = у, приходим к линейной системе х=у, у=-2у — 5х с особой точкой (О, О). Поскольку корни характеристического уравнения этой системы равны -1З:26 то особая точка — устойчивый фокус. Положив в системе х = 1, у = О, находим вектор фазовой скорости е = (О; — 5), с помощью которого мы, учитывая устойчивость фокуса, устанавливаем направление закручивания траекторий на фазовой плоскости (рис.67). и 658. й+*+2х- '=О.
м Полагая в данном уравнении х = у, получаем систему дифференциальных уравнений: 2 в=у, у=в — 2х — у, из которой следует, по точки (О, 0), (2, 0) — особые. Обычное исследование их показывает, по точка (О, 0) — устойчивый фокус, а точка (2, 0) — седло, причем прямые у2 = х — 2, уг = -2х+ 4 Ряс. 42 являются касательными к интегральным кривым, входящим в него (рис,бе). Далее, судя по знаку производной х — 2х — у 2 у = у устанавливаем грубую картину интегральных кривых (рис.
69). Заметим, что на параболе у= а — 2х интегральные кривые достигают экстремальных значений; ось Ох они пересекают под прямым углом, а ось Оу — под углом 45'. Используя эти данные, строим картину интеграяьных кривых (рис.70). и 659. х + хз — х + 1 = О. м Из уравнения * — у — 1 = О 2 2 2 2 Ну х — у — 1 Рве. бв у (е) находим области возрастания и убывания интегральных кривых на фазовой плоскости Оху (рис.71).
Кривые 311 $3.Фазоваа плоскость служат касательными к интегральным кривым, входящим в него. Р .го Особые точки показаны на рис. 72. Таким образом, исходя из сказанного н рисунков 71, 72, строим фазовый портрет (рис. 73). )ь Рас. 71 Ряс. 72 Ра . тЗ 660. х+5х — 41п = О. х +1 2 < Исключив параметр ! из системы дифференциальных уравнений х +1 *=У~ у=4)п 5у 2 придем к уравнению г(у 4!и — * — 2+ — — 5у (1) г!х у с особыми точками: (-1, 0) н (1, 0). Решая неравенства 4 1п — *+- — 5р сО, у устанавливаем области монотонности интегральных кривых уравнения (1) на фазовой плоскости (рис.
74). Кривую Рас. 74 4 х+1 5 2 интегральные кривые пересекают под нулевым углом, а кривую у = 0 — под прямым. Теперь исследуем особые точки. пересекаются интегральными кривыми под нулевым, а ось Ох — под прямым углом. Далее, известным способом находим две особые точки: ( — 1, 0), (О, 1). Поскольку корни характеристического уравнения, соответствующего укороченной линейной системе Р'=О, й=-24, у мнимы, то особая точка (-1, 0) мотает быть для системы г х=у, у=х — у — 1 либо фокусом, либо центром.
Принимая во внимание то, что прн замене у на -у уравнение (4) вида не меняет, Убеждаемся, что х рассматриваемая особая точка является центром для указанной системы. Точка (1, 0) является седлом, а прямые дг — — ьг2(х — 1)г уг = -ггг2(х — 1) Гл. б. Устойчивость и базовые траеаторви 2!г Полагая в системе (е) х = 1+1, у = 0 и отбрасывая нелинейные члены, получаем укороченную систему: ь=б, 6=4ч — 50. Поскольку корни характеристического уравнения -5 х т/41 Л, 7 — (Л, ге 0,70; Лг м -5,7), то особая точка (1, 0) — седло.
Прямые у = Л,(х — 1), у = Л,(х — 1) тис 75 являются касательными к интегральным кривым, входящим в эту точку. Аналогично, полагая в (з) и = = — 1 + (, у = О, получаем укороченную систему: 4=0 6= 44 — 57) характеристический определитель которой имеет нули Л,=-1, Л7=-4.
Следовательно, особая точка (-1, 0) — узеа. Из уравне- ( -4 — 5 — Л|) (е,) находим, что все интегральные кривые, проходящие через узел, касаются прямой у = — х — 1, а из уравнения ( -4 -5 — Л~) (ег) следует, по все указанные кривые пересекаются интегральной кривой, проходящей через узел и имеющей в нем касатезшную у = -4х — 4. Окрестности особых точек изображены на рис. 75. Наконец, учитывая все полученные данные, строим полную фазовую картину (рис.
76). ~ 661. х = 4 — 4х — 2у, у = ху. и Из системы Гас. 76 хр >О 4 — 4х — 2у устанавливаем области монотонности интегральных кривых (рис. 77). Для более четкою представления о поведении интегральных кривых рассмотрим знак второй производной 4У(у — У1)(У вЂ” Уз) (4 — 4х — 2У)з з з где у, = 2 — х+ от, уз = 2 — х — х 7 (х > 0). если х ( О, то 4У (у + у(2х — 4)+(1 — х)(4+ х )) (4 4х гу)з Решая неравенства зО прях>0 (4 — 4х — 2и)з 4 — 4х — 2У=О, ау=О находим особые точки: (1, 0), (О, 2). Известным способом (см.
выше) устанавливаем, что (1, 0) — седло, а точка (О, 2) — вырожденный узел; Гас. 77 т причем в седло входят интегральные кривые под углом о = — 4-. Далее, ось Оу интегральные кривые пересекают горизонтально, а прямую 4-4х-2У = 0 — вертикально. Решая неравенства и ++++ ++++ ++++ ++ ч+ ++++ 2 +++++ +++++ +++в+ +++ +++ +++ +++ +++ +++ +++ +++ — — — +++ ./. + Ъ~ — — — // + + ++ — 'к т+ ++++ 1 ++++ ><О при х< 0, (4 4 2У)з находим области выпуклости интегральных кривых (Рис.78). Заметим, что на прямой 4 — 4х — 2у = 0 и на кривых У = Уг, у = уг часть траекторий меняет направление выпуклости.
Теперь проследим за интегральными кривыми, проходящими через особые точки. Через седло проходит прямая у = О и кривая, пересекающая ось Ох лод углом а = — а/с!8 з. Эта кривзя на параболе у = у, имеет точку перегиба и пересекает ось Оу горизонтально, а прямую 4 — 4х — 2У = 0 — вертикально.
Наконец, под углом 135' к оси Ох она входит в узел. Нижняя часп, ее (при у < 0) проходит под параболой у = у„поскольку; +++ +++ + + .~- +е+ +++ +++ +++ + ++ сг Ъ е+ х+ г +++ +++ м — а ш ху у = 4 — 4х — 2У х на -х (ради удобства), перейдем к интегральному уравне- нию кривой (/у): у(х) = е + / (у(1) <й (х ) )0).