1612042534-1f807c74b48d50003337997455390972 (542297), страница 27
Текст из файла (страница 27)
(НГУ)ЛЕКЦИЯ 16Новосибирск, 2018 г.21 / 24Если для уравнений возмущённого движения можно найти функциюЛяпунова, обладающую по сравнению с обычной устойчивостьюнекоторыми дополнительными свойствами, то движение будетасимптотически устойчивым.Теорема Ляпунова. (об асимптотической устойчивостиневозмущённого движения)Если дифференциальные уравнения возмущённого движения таковы,что можно найти знакоопределённую функцию V (t, x), допускающуюбесконечно малый высший предел, производная по времени которойV 0 (t, x), взятая в силу дифференциальных уравнений возмущённогодвижения, является знакоопредёленной функцией противоположного сV знака, то невозмущённое движение — устойчиво асимптотически.Батяев Е. А.
(НГУ)ЛЕКЦИЯ 16Новосибирск, 2018 г.22 / 24Итак, теоремы Ляпунова дают достаточные условия устойчивостидвижения. Применение этих теорем требует знания функции ЛяпуноваV (t, x), обладающей вполне определенными свойствами для обычнойили асимптотической устойчивости.Однако, общих методов построения таких функций нет и вопрос с еенахождением является совершенно открытым. Известно лишь что дляустойчивых движений такая функция заведомо существует.Это обстоятельство ограничивает практическое применение теорем.Тем не менее значение теорем Ляпунова не исчерпывается тем, чтоона дает средство для творческого решения задач об устойчивостидвижений. С помощью них устанавливается ряд других важныхтеоретических результатов.Батяев Е. А.
(НГУ)ЛЕКЦИЯ 16Новосибирск, 2018 г.23 / 24Во многих практически важных случаях функцию Ляпунова можнопостроить если известны первые интегралы уравнений возмущённогодвижения. Пусть U1 , . . . , Up – не зависящие от времени первые интегралыуравнений возмущённого движения. Без ограничения общности можносчитать, что функции Uj (x1 , . . . , xm ) (j = 1, . . . , p) обращаются в ноль вначале координат x1 = . . . = xm = 0. Пусть в общем случае ни одна из Uj неявляется знакоопределенной.
Будем искать функцию Ляпунова в виде:V =pX(λj Uj + µj Uj2 )j=1где λj , µj – неопределенные постоянные. Понятно, что V будет первыминтегралом уравнений возмущённого движения (т.е. V 0 ≡ 0 на любомрешении дифференциальных уравнений возмущённого движения). Если намудастся подобрать λj , µj так, чтобы разложение V в ряд по возмущениямначиналось со знакоопределенной квадратичной формы, то все условиятеоремы Ляпунова об устойчивости движения будут выполнены.Батяев Е. А. (НГУ)ЛЕКЦИЯ 16Новосибирск, 2018 г.24 / 24ТЕОРЕТИЧЕСКАЯ МЕХАНИКА2 СЕМЕСТРЛЕКЦИЯ 17УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМТЕОРЕМА ЛЯПУНОВА ОБАСИМПТОТИЧЕСКОЙ УСТОЙЧИВОСТИНЕВОЗМУЩЕННОГО ДВИЖЕНИЯ ПОПЕРВОМУ (ЛИНЕЙНОМУ) ПРИБЛИЖЕНИЮКРИТЕРИИ АСИМПТОТИЧЕСКОЙУСТОЙЧИВОСТИЛектор: Батяев Евгений АлександровичБатяев Е. А.
(НГУ)ЛЕКЦИЯ 17Новосибирск, 2018 г.1 / 25Линейная системаПусть дифференциальные уравнения возмущенного движения –линейны и имеют постоянные коэффициенты, т.е. образуютлинейную систему:mXdxk (t)= Xk (t, x1 , . . . , xm ) =aki xi (t)dti=1(k = 1, . . . , m)Или вводя вектор x̄(t) и постоянную матрицу A получимпредставление линейной системы в матричном видеx1a11 . . .
a1mdx̄(t) . ... ..= Ax̄(t)x̄ = .. A = .... dtxmam1 . . . ammБатяев Е. А. (НГУ)ЛЕКЦИЯ 17Новосибирск, 2018 г.2 / 25Т.к. это система линейных обыкновенных дифференциальныхуравнений, то ищем общее решение в видеx̄ = ūeλtПодставляя это выражение в уравнении и сокращая на eλt ,получим соотношения, связывающие искомые величинывектора-столбца ū = (u1 , . . . , un ) и λ:(A − λI)ū = 0где I – единичная матрица. Нас интересует нетривиальноерешение данной однородной системы линейных алгебраическихуравнений , т.е.
|ui | 6= 0 хотя бы для одного номера i = 1, . . . , mnXили kūk =|ui |2 > 0, значит определитель системы (матрицыi=1A − λI) должен быть равен нулю:det(A − λI) = 0характеристическое или вековое уравнениеБатяев Е. А. (НГУ)ЛЕКЦИЯ 17Новосибирск, 2018 г.3 / 25Характеристическое уравнение det(A − λI) = 0 являетсяалгебраическим уравнением m-ой степени относительно λ.Корни характеристического уравнения (λ) называются –характеристическими числами матрицы A.Вектор-столбец ū 6= 0, удовлетворяющий вместе с числом λсистеме линейных однородных алгебраических уравнений(A − λI)ū = 0называется – собственным вектором матрицы A,соответствующим характеристическому числу λ.Итак, в каждом частном решении исходной линейной системыобыкновенных дифференциальных уравнений, имеющим видx̄ = ūeλt , λ – характеристическое число матрицы A, а ū –соответствующий собственный вектор.Батяев Е.
А. (НГУ)ЛЕКЦИЯ 17Новосибирск, 2018 г.4 / 25В общем случае, если про матрицу A ничего не известно, токорни характеристического уравнения могут быть ивещественными и комплексными, кратными и не кратными.Если, например, все корни различные λk (k = 1, . . . , m), тогдакаждому λk соответствует собственный вектор ūk и частноерешение вида ūk eλk t .
Линейная комбинация этих решений сновабудет решением:mXx̄ =Ck ūk eλk tk=1где Ck – постоянные.Легко показать, что все ūk – линейно независимы, а это решениеявляется общим, т.е. удовлетворяющим уравнениям и начальнымданным путём соответствующего подбора коэффициентов Ck .Батяев Е. А. (НГУ)ЛЕКЦИЯ 17Новосибирск, 2018 г.5 / 25В самом же общем случае, корни характеристического уравнениямогут быть кратными. В курсах по теории дифференциальныхуравнений доказывается, что в таком случае, формула для решенияможет усложниться. Пусть например, какой-либо корень λj , кратностиkj , у которого имеется только nj (nj < kj ) линейно независимыхсобственных векторов ūj . Если nj = kj , т.е. хоть корень λj и кратный,но из системы (A − λj I)ū = 0 можно найти столько независимыхсобственных векторов ū какова кратность корня, что возможно, еслиrang(A − λj I) = m − kj , и обычный вид общего решения сохранится.Иначе недостающие решения для λj можно искать в виде векторногомногочлена относительно t степени rj = kj − nj (метод Эйлера)´³r(r )(0)(1)x̄j = ūj + ūj · t + .
. . + ūj j · trj · eλj t = Pj j (t) · eλj t(1)(r )Чтобы найти векторы присоединенные векторы ūj , . . ., ūj j –необходимо подставить данное решение в исходную системудифференциальных уравнений. Приравняв коэффициенты подобныхчленов (одной степени по t) в левой и правой частях системы,получим уравнения для нахождения этих векторов.Батяев Е.
А. (НГУ)ЛЕКЦИЯ 17Новосибирск, 2018 г.6 / 25r(r )(1)Члены данного полинома Pj j (t): ūj · t, . . . , ūj j · trj –неограниченно возрастающие со временем, называются – вековые.Таким образом, в общем случае произвольное решение системыдифференцируемых уравнений определяется формулой:mm³´XX(0)(1)(r )x̄ =Ck x̄k =Ck ūk + ūk · t + . . . + ūk k · trk · eλk tk=1k=1Для дальнейшего анализа решения на устойчивость, вспомнимизвестный факт из математического анализа:lim tn e−µt = 0t→∞∀ n, µ > 0Рассмотрим самый общий вид корня характеристического уравнения –характеристического числа – комплексного:λk = µk + iνkт.е. µk = Re(λk ) – вещественная часть λk . Тогда решение примет вид:m³´X(0)(1)(r )x̄ =Ck ūk + ūk · t + . .
. + ūk k · trk · eµk t · (sin νk t + i cos νk t)k=1Батяев Е. А. (НГУ)ЛЕКЦИЯ 17Новосибирск, 2018 г.7 / 25Рассмотрим различные случаи.1. Пусть все µk < 0. Тогда имеет место оценка:m³´X(0)(1)(r )|x̄| <|Ck | |ūk eµk t | + |ūk · teµk t | + . . . + |ūk k · trk eµk t | −−−→ 0t→∞k=1значит lim |x̄(t)| = 0, и нулевое решение линейной системыt→∞дифференциальных уравнений возмущенного движения, аследовательно, невозмущённое движение – асимптотически устойчиво.2. Если хотя бы один µk > 0, то очевидно, что соответствующее λkчастное решение, а значит и общее решение линейной системыдифференциальных уравнений – будет неограниченно возрастать.Следовательно, в этом случае невозмущённое движение – неустойчиво.3. Пусть нет ни одного корня с положительной вещественной частью,но есть корни с нулевыми вещественными частями.
Тогда вопрос обустойчивости невозмущённого движения определяется наличиемвековых членов соответствующего частного решения: если они есть, торешение неустойчиво, если их нет – то устойчиво. При этом речь вэтом случае идет – об обычной (не асимптотической) устойчивости.Батяев Е. А.
(НГУ)ЛЕКЦИЯ 17Новосибирск, 2018 г.8 / 25Линейное приближениеПредставим исходную систему уравнений возмущенного движения(вообще говоря, нелинейных) в матричном виде:dx̄(t)dxk (t)= X̄(t, x̄) ⇔= Xk (t, x1 , . . . , xm ) (k = 1, . . . , m)dtdtОтносительно X̄ = (X1 , . . . , Xm ) – полагаем, что они аналитичны.Тогда раскладывая каждую в ряд по возмущениям в окрестностиначала отсчета (нулевого решения) – в ряд Маклорена – получим:∂ X̄1 ∂ 2 X̄(t, 0) · x̄ +(t, 0) · x̄ · x̄ + .
. .∂ x̄2 ∂ x̄∂ x̄Первое слагаемое X̄(t, 0) ≡ 0, что соответствует невозмущённомудвижению x̄ = 0. Обозначим матрицу во втором слагаемом и суммуостальных членов ряда начиная со второго порядка относительно x̄,следующим образом:X̄(t, x̄) = X̄(t, 0) +A(t) =∂ X̄(t, 0)∂ x̄Батяев Е. А. (НГУ)f̄ (t, x̄) =ЛЕКЦИЯ 171 ∂ 2 X̄(t, 0) · x̄ · x̄ + . .
.2 ∂ x̄∂ x̄Новосибирск, 2018 г.9 / 25Таким образом получим систему следующего вида:dx̄(t)= A(t)x̄ + f̄ (t, x̄)dtf̄ (t, x̄) – тоже аналитичная функция.Отбрасывая все нелинейные члены f̄ (t, x̄), получим системулинейных обыкновенных дифференциальных уравнений:dx̄(t)= A(t)x̄dtназываемую – линейным приближением – длянелинейной системы.Батяев Е. А. (НГУ)ЛЕКЦИЯ 17Новосибирск, 2018 г.10 / 25Различают:периодический случай:X̄(t, x̄) = X̄(t + τ, x̄)⇒A и f̄ − периодические, с периодом τ .X̄ = X̄(x̄)стационарный случай:⇒A − постоянная, f̄ = f̄ (x̄);В обоих этих случаях Ляпуновым было установлено:• из асимптотической устойчивости линейного приближения⇒ асимптотическая устойчивость нелинейной системы;• из неустойчивости линейного приближения⇒ неустойчивость нелинейной системы.В случае обычной устойчивости линейного приближения сделатьзаключение об устойчивости нелинейной системы нельзя.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
