Filippov (531376), страница 16
Текст из файла (страница 16)
При отклонении груза на расстонние х, пружина действует на него с силой кх, направленной к положению равновесия. Сила трения равна 7 = сопз1 и направлена в сторону, противоположную скорости 107 ГС 17. сразовап плоскость Указание. При малых колебаниях считать вгаз з. Изменение длины маятника происходит мгновенно (скачком), при этом угол отклонения маятника и его момент количества движения относительно оси не испытывают скачков. Начертить на фазовой плоскости траектории систем 1040 — 1046, записанных в полярных координатах, и исследовать, имеются ли предельные циклы. 1040.
с'т = т(1 — тз), с11 1041. —" = т(т — 1Нт — 2), ~ = 1. с11 сМ 1042. — т = с(1 — с)з. сЫ вЂ” = 1. сМ 1043. — = юптс с1т сЫ 1044. —" = т()т — 1) — (т — 2! — 2т+ 3), т = 1. с1ьс с11 1045. — = та1п-. бт .1 с1с 1046. — = т(1 — с) з1п с1т 1 Ф 1 — с' 1047*. При каких условиях система груза. При 1 = 0 груз находится на расстоннии сс от положения равновесия н имеет нулевую скорость. Вывести уравнение движения груза. Приняв т = 2, 1с = 2, 1 = 1, 6 = 5, изобразить движение груза на фазовой плоскости. 1039. Изобразить на фазовой плоскости малые колебания мантника переменной длины, считая, что при движении маятника вверх его длина равна 1, а при движении вниз равна Л > 1. Во сколько раз увеличитсн амплитуда за одно полное колебание? (Пример: раскачка качелей.) 108 З 18. Зависимость решения от начальных условий где функция 1(г) непрерывна, имеет предельный цикл? При каких условиях этот цикл устойчив? Неустойчив? Полуустойчив? 1048*.
При каких значениях постоннной а система с[~о с[г 2 — = 1, — = (г — 1)(о + з[п ~р) сМ ' сМ имеет устойчивый предельный цикл? Неустойчивый? Для уравнений 1049 — 1052 с помощью изоклин построить траектории на фазовой плоскости и исследовать особые точки. По чертежу сделать заключение о поведении решений при 1 -+ +ос и о возможности существовании замкнутых траекторий.
1049 -.+ .3 х+ . О 1050. х+ (яз — 1)об+ го = О. 1051. т, +х — 2вгс18х+ з: = О. 1052. У+ 2* — х+х = О. 1053*. Длн уравнения х + 2ах — Ьзкпх + х = О (О < а < 1, Ь > 0) построить траектории на фазовой плоскости и найти точки, в которых предельный цикл пересекает ось О:с. Указание. Найти зависимость между абсциссами двух последовательных пересечений траектории с осью Ох. 1054.
Показать, что уравнение х+ Г(х?) +х = О, где функцин г' непрерывна и г'(у) > 0 при у > О. Г(у) < О при у < О. не может иметь предельных циклов на фазовой плоскости. Указание. Исследовать знак полной производной л,(хе+у ). 1055*. Пусть 1(х, у) и Д, )'„' непрерывны, с"(О, 0) < О, а при хз+ уз > Ьз имеем 1(х, у) > О.
Доказать, что уравнение х, + Ят,. х)х -~- т, = О имеет периодическое решение х(с) ф О. Указание. Перейти на фазовую плоскость и исследовать знак полной производной й(х~ + у ). Построить кольцо, из которого не может выйти ни одна траектория. Применить теорему 21 из [3). Ц18.
Зависимость решения от начальных услоаий 109 З 18. ЗАВИСИМОСТЬ РЕШЕНИЯ ОТ НАЧАЛЬНЫХ УСЛОВИЙ И ПАРАМЕТРОВ. ПРИБЛИЖЕННОЕ РЕШЕНИЕ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ !. Рассмотрим систему в векторной записи — =йг ) е)х с)4 ЦУ(4, у) — У(й х)Ц < ЛЦу — хЦ. (2) Через Ц Ц обозначается любая из обычно применнемых норм вектора: ы = 'т Рт .. +) .г ЦхЦ = (хг)+ ... Ч- )х„) или ЦхЦ = шах )х,Ц Пусть х(Г) — решение системы (1), а у(Г) — вектор-функция. удовлетворнющан неравенствам — У(Г, у)! < у. Цу(П) — х(б)Ц < б. !. г1у Тогда имеет место оценка Цх(Г) — у(1)Ц < бежи -~- — (еь'Ц вЂ” 1) .
(8) Это неравенство можно применять для грубой оценки ошибки приближенного решении у(Г) системы (1), а также дли оценки сверху разности решении х(1) системы (1) и решения у(Г) системы ллг = 8(й у), если Цу(Г, у) — )(Г, у)Ц < у. 2. Если в системе уравнений е)г — =)е(1эхш ...,х„,р), 4=1, ...,и (4) гЕсли в выпуклой по х области имеем — ~ < а (П б = 1, ...,а), то оь в этой области выполнено условие Липшица с й = па. где х = (хг, ..., х ), 1 = (А, ..., 1„). Пусть в рассматриваемой области вектор-функции г непрерывне па й х и удовлетворнет усло- вию Липшица по х 110 518. Зависильосшь решения от нач льнах условий с начальными условиями х,(0) = а,(р), ь = 1..... п (5) р является параметром, функции 7, и а, (1 = 1.....
и) непрерывны и имеют непрерывные производные по хы ..., х„, р, то решение имеет непрерывную производную по параметру р. Производные — '-"' = ип ь = 1, ..., и удовлетворнют линейной системе уравнений й, " ВУь аУь (б) и начальным условиям и;(О) = а',(р), 1 = 1, ..., и.. Значении производных — 7ь и Ось в формуле (б) берутсн при хь = хд(1), ..., х„= = х„(1), где хд(1), ..., х„(1) — решение системы (4) с начальными условинми (5). В частности, если положить аь(р) = р, а;(р) = сопзФ при ь ф й и считать, что все функции 1ы ..., 1 не зависят от р, то из предыдущего утверждения будет следовать, что для системы (4) с начальными условиями х;(О) = аа ь = 1, ..., п производные -~- = о „ = и; (ь = 1, ..., и) от компонент решении хм ..., х,„по начальному условию аь существуют и удовлетворяют системе уравнений пи, д1ь а ~ьо — '=7 *и„ь=1, ...,и, у=з и начальным условиям и;(О) = 0 при ь ~ 5, ал(0) = 1.
3. Если в (4) и (5) функции ~; и а, имеют непрерывные производные по хм ..., х„, 1ь (вблизи значения р = О) до порндка пь включительно, то решение тоже имеет непрерывные производные по р до порядка т и, следовательно. разлагается по степеннм параметра р по формуле Тейлора: х(1) = оо(1) -Ь ро,(С) -Ь р'оз(1) + ... + ршэ, (1) -Ь о(рш). (7) Здесь х и о; и-мерные вектор-функции.
Чтобы найти функции о,(1), можно разложить правые части в (4) и (5] по степеннм р, подставить туда разложение (7) и приравннть коэффициенты при одинаковых степенях р. Получим систему дифференциальных уравнений, из которой последовательно определнются оо(1), о1(т), ... В случае, когда ть и а; — - аналитические функции от хы ... ..., х„, р, решение х(1) разлагается в сходящийся при малых р степенной ряд по р (в силу теоремы об аналитической зависимости 318.
Зависимость решения от начальных условий 111 решения от параметра, см. (4), гл. 1, 2 6). Коэффициенты этого ряда совпалают с коэффициентами разложении (7). Изложенный метод можно использовать для отыскания решении дифференциального уравнения при малых р в тех случаях, когда при р = О уравнение решаетсн известными методами. П ример. Разложить по степеням параметра р решение задачи х=х +2рс, х(1) = — 1. (8) ищем решение в виде х(1) = оо(г) + рог(1) + р ог(1) + . Подставлян это в (8) и приравнивая коэффициенты при одинаковых степенях р, получаем систему ее = оо ее(1) = — 1, ог = 2ооог + 21, ог(1) = О, ег = 2соог г ог, ег(1) = О, Из первого уравнения и начального условия находим ео(1) = — 1 Подставляя это во второе уравнение, получаем чи = -21 ги -~-21, ог(1) = О.
Отсюда е,(1) =1 — 1 '. Подставлян найденные ео и чи в третье уравнение, получаем ег = — 21 ог + (1 — 1 ), ог(1) = О. Решив это линейное уравнение н воспользовавшись начальным условием, найдем ог(1) = — ' — г + з, — г. Следовательно, решение задачи (8) имеет вид / 1'1,/1 2 8 х(1)= — — +р(1 — — )+р ( — — — + — — — )+о(р). сг/ (3 С 31г сг) Это разложение можно продолжить дальше тем же способам. Аналогичным методом можно получать разложения по степеням параметра периодических решений нелинейных уравнений„ в частности, уравнений вида (О) х+ а х = рт(С х. х. И), где функция 7 периодическая по й Переходить от уравнения 2-го порядка н системе при этом ненужно. Произвольные постоянные, 112 218.
Зависилгоста решен я от начальных условий возникающие при отыскании оо(1), ог(1), ..., определяются уже не из начальных условий. а иг условий периодичности (см. [4), гл. 2. 5 8). В случае, когда правая часть (9) не зависит от 1, период решения х(1) заранее не известен. Тогда в уравнении (9) надо перейти от 1 к новому независимому переменному т = 1(1+ Ьгр+ Ьгр + ... ) и искать решения х(т) периода 2о/а. коэффициент Ьг обычно определяется из условия существования периодического решения длн иг(т), и т.д.