Введение в теорию устойчивости движения. Д. Р. Меркин (774242), страница 42
Текст из файла (страница 42)
Тогда из неравенств (7.42) найдем ,гВ ~Ь вЂ” А, ЗЬ+4ар Ь + а (7 44 Если границы а, А, Ь, В функций а ((, л, в) и р ((, х, и) удовлетворяют неравенствам (7.42) для всех л, й, (,;ь то будут выполнены условия теоремы Барбапгкпа — Красовского '). В атом случае невозмущенное движение л = О и й = О будет асимптотически устойчиво в целом, т. е. пРи любых начальных возмУщениЯх лз и йа в). й 7.5.
Устойчивость линейных систем с периодическими коэффициентами Исследование устойчивости движения многих систем, встречающихся в рааличных технических задачах, часто сводится к анализу линейных дифференциальных уравнений с периодическими коэффициентами. В матричной форме эти уравнения могут быть записаны так (см. 3 5.2, формула (5А9а)): (7.45) х = Р (()ль В атом уравнении л — матрица-столбец или вектор хг (7.46) ') Теорема Барба~лина — Красовского сформулирована в $2.8 лля автономных свстем.
Пряфуякцин у(г), не зывнсящей явно от зреызян г в удовлетворяющей условию (2.16), зга теорема остаегся справедливой н в тоы случае, когда производная )г, завися явно от времеви, является определеяно-огрвцагельвой функцией з смысле Ляпунова. ") Для случая, когда функции гг к [) заввсят только ог времени Г, условия (7.44) были получены другим методом В. М. Старящаским [47). Пряввдзниый здесь вывод опубликован в работе [39[.
232 гл ч11. ъстОнчнвость неавтОномных систем а Р (2) — квадратная матрица раа(1) ... р „(1) ~ (1.47) РН)=-! Будем считать, что все элементы ргз (Т) матрицы Р (2), а следовательно, и сама матрица Р являются периодическими функциями времени 2 одного и того же периода Т, так что для любого момента времени 1 справедливо равенство Р (~ + Т) = Р (2). (7.48) Для получения критериев устойчивости таких систем кратко остановимся на некоторых общих вопросах теории линейных дифференциальных уравнений с периодическими коэффициентами, принадлежащей Флокс(г")оа)пей).
Совокупность и линейно невависимых решений уравнения (7.45) ха1 Х1 аа (7.49) а Хаа называется фундаментальной системой решений этого уравнения, а матрица Х (2) = Ц Х„..., Х„Ц (7.50) или, что то же самое, хв хм ха а х'и Х(2) = (7.51) хт хт "хьа называется фундаментальной матрнцей.
Здесь, как и в дальнейшем, первый индекс элемента хь, обозначает номер функции, а второй — номер решения. Общее решение х (2) уравнения (7.45) определяется обычной формулой общего решения линейного однородного дифференциального уравнения: (7.52) Х(2)=-С1Х1+ С1Х +... + С„Х„, 1 7.ь системы с певиодическими коэФФициентлми 233 .где С„..., ф— произвольные' постоянные, определяемые из начальных условий. В матричной форме общее решение (7.52) имеет вид ш (1) = Х (1).С, (7.53) где С вЂ” матрица-столбец с (7.54) с„ Не нарушая общности, можно считать, что фундаментальная система решений удовлетворяет следующим начальным условиям: (1, й=у, х„;(О)= ~ ~0, -дУ или в матричной форме Х(0) =Е, (7.55) где Š— единичная матрица 1 О...
О О 1... О О О... 1 Обозначим через Д (1) определитель фундаментальной матрицы Д (1) = се1 Х (1). ' .56) Пользуясь равенством (7.55), получим Д (О) = сеС Х (О) = ое1 Е = 1. (7.57) В теории дифференциальных уравнений доказывается следующая формула Лнувилля: с 5'Оь+-+ ..М' Д (1) =. Д (О) ео Учитывая равенство (7.57), найдем, что при 1 =- Т т ) О„+, -ер,„ы Д (7') — е" (7.58) 234 гл, чп, устойчивость нклзтономньгл снетки Если в каком-нибудь решении хт заменить г на г + Т, то в силу периодичности матрицы Р (~) мы снова получим решение, так как вектор хт (г + Т) будет по-прежнему удовлетворять уравнению (7.45), если ему удовлетворял вектор хз (т). Полученное решение не будет совпадать с первоначальным решением ха (Г), но, как всякое решение уравнения (7.45), оно может быть получено из общего решения (7.53) выбором соответствующей матрицы-столбца С.
Обозначив зту матрицу через Аю найдем х„(г + Т) = Х (г).А, (я = 1, 2,..., и). Отсюда видно, что фундаментальная матрица Х (с + Т) решений хт (г + Т), хт (г + Т),..., х (С + Т) имеет вид Х (~+ Т) = Х (~)А, (7.59) где А — постоянная матрица: ап ам ° а„ ам а„... а,„ А =((Ап ..., А„()= (7.60) а,о а„~ По предположению фундаментальная система решений х, (С),..., х„(Ю) удовлетворяет начальным условиям (7,55).
Позтому, положив в равенстве (7.59) й = О, найдем Х (Т) = Х (0)А = ЕА = А. Следовательно, если известна фундаментальная матрица Х (г), то матрица А определится из последнего равенства: тн (Т) ... и „(Т) (7.6() А=-Х(Т)= Ую (Т) ° хам (Т) Покажем далее, что существует решение х (г), удовлетворяющее следующей зависимости: х(~+ Т) = рх(~), (7.62) где р — некоторое постоянное число (такое решение называется нормальным). Действительно, любое решение уравнения (7.45) можно получить из общего решения (7.53). Поэтому если нормальное решение существует, то должна существовать такая постоянная матрица- С е5.
системы с пегиодическими кОэФФициентАми 235 столбец (7.63) для которой будет справедливо равенство (С) = Х (С)(). Так как по предположению х (С) удовлетворяет равенству (7.62), то, имея) в виду, что ю (С+ Т) = Х (С+ Т)(3, получим Х (С + Т)р = рХ (г)(С или, пользуясь р««енством (7.59)„ Х (С)А() = рХ (С)(). Группируя члены, найдем Х (С) (А — рЕ)р = О. Учитывая, что это равенство должно выполняться при всех С, получим (А — рЕ)() Это матричное уравнение, являются матрица-столбец () и скалярным уравнениям (а„— р)(), + а,ар, + аа7р7 + (а„— р)ра + а„7р +а~,К+...+ в котором неизвестными число р, эквивалентно в ...
+а,„(7„= О, ... +«,„Р„=О, (а — р) р„= О. ам — р а7а 1« ам ам — р с)еС (А — рЕ) = = О. (7.64) а«7 ааа ° ° ° а„„вЂ” р Каждому корню ра этого характеристического уравнения (для простоты будем предполагать, что среди корней нет кратных) отвечает свое решение ю„(С), удовлетво- Для того чтобы система и алгебраических однородных уравнений относительно р„..., р„имела решение, отличное от нуля, необходимо и достаточно, чтобы определитель системы равнялся нулю: 236 Гл.
Тп. устОЙчиВОсть нелитономссьсх систем ряющее условию (7.62). В результате получим и нормальных решений х, (1),..., х„(1), удовлетворяющих условию (7,62); х, (1 + Т) =- р,г, (1),..., х„(1+ Т) = р х„(1), (7.65) Эта система решений линейно независима и может быть принята за фундаментальную систему решений.
Покажем теперь, что нормальные решения имеют вид х» (1) = е"к'срк (1) ()с = $,..., и), (7.66) где срс (1) — периодическая матрица-столбец периода Т ) св, (с) ~ фк (1) = . ~, фк (1 + Т) = фк (1), (7.67) ч р)! а ак — постоянные числа, называемые хараитериетичеекими показателями, вычисляются по формуле с) а»= —,, 1пр». (7»=1,..., и). (7.68) 1 х(1)=. ~ Скхк(1). с.=с (7.69) с) Напомним определение логарифма коьшлексного числа (корни уравнения (7.64) могут быть комплексными числами): 1н р = 1п ( р ) + с агх р, где ) р ) есть модуль р и агп р обовначает аргумент р.
Действительно, внесем в равенства (7.66) значение 1, равное 1+ Т. Получим х» (1 + Т) = е »П'тзфк (1+ Т). Пользуясь равенствами (7.67) и (7.68), последовательно найдем (' )т х»(1+ Т) =е ге'с*фа(1-1- Т)=е"к'е ' к срк(1)= = Ркхк (1)» т. е. равенства (7.65). Теперь можно перейти к исследованию устойчивости движения. Будем считать, что за фундаментальчую систему решений уравнения (7.45) принята система нормальных решений хс (1),..., х„(1), удовлетворяющих условиям (7.65). Общее решение уравнения (7.45) запишем в форме (7.52) З еа системы с пеукодическими кОэФФициентАми 237 Вектор х (») определяет изображающую точку М, а слагаемые Ськг (») являются составляющими его.
Через период Т положение изображающей точки М будет определяться равенством х (» + Т) = лг Сгеек (» + Т). з=г Составляющие этого вектора, согласно (7.б5), равны Сьх„(» + Т) = ркСкх„(»). (7.70) Ото»ода ) Сею„(» + Т)(=-) р„() Скос„(») ~. Это равенство показывает, что если все ~ рь ~ 1, то через период Т модули всех составляющих вектора х (»+ Т) уменьшатся и, следовательно, изображающая точка М приблизится к началу координат; если модуль хотя бы одного корня рь больше единицы, то через период Т соответствующая составляющая Сглаз (» + Т) вектора ю (» + Т) увеличится по модулю и изображающая точка М начнет отдаляться от начала координат; наконец, если среди корней характеристического уравнения имеются такие, модули которых равны единице, то модули соответствующих составляющих вектора ж (» — , 'Т) останутся беа изменения. Из этого вытекают следующие условия устойчивости системы, возмущенное двигкение которой описывается линейными дифференциальными уравнениями с периодическими коэффициентами.
Если модули всех корней характериспгического уравнения (7.б»ь) меньше единицы, то невозмущенное движение хг =... = х„= О асимппютически устойчиво. Если среди корней характеристического уравнения имеется хотя бы один, модуль которого больше единицы, то невозмущенное движение неустойчиво, Если среди корней характеристического уравнения имеются такие, модули которых равны единиг»е, а модули остальных корней меньше единицы, то нсвозмущенное движение устойчиво, хотя и не асимптотически.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
