Математические методы исследования устойчивости семейства систем дифференциальных уравнений с последействием тема диссертации и автореферата по ВАК РФ 01.01.09, кандидат физико-математических наук Купцов, Сергей Юрьевич

• Купцов, Сергей Юрьевич
• кандидат физико-математических наук
• 1998

Оглавление диссертации кандидат физико-математических наук Купцов, Сергей Юрьевич

Введение

Обозначения

Глава 1. Некоторые свойства матриц

Глава 2. Оценка радиуса локальной устойчивости системы по параметрам.

Глава 3. Анализ устойчивости гироскопической системы с одним запаздыванием.

Глава 4. Абсолютная устойчивость линейной системы у * и и дифференциальных уравнении с последействием по запаздываниям.

Введение диссертации (часть автореферата) на тему «Математические методы исследования устойчивости семейства систем дифференциальных уравнений с последействием»

За последние 45-50 лет теория систем дифференциальных уравнений с отклоняющимся аргументом нашла многочисленные приложения в самых разнообразных областях механики, физики, биологии, технических и экономических наук. Особенно велики приложения к теории автоматического регулирования и к теории колебаний. В связи с быстрым развитием и расширением областей приложения уравнений с отклоняющимся аргументом естественно возрастает и интерес к ним. За этот период опубликовано, наверное, тысячи работ, в которых рассматривались данные объекты, создавался аппарат их исследования, вводились все новые и новые обобщения. Было бы бессмысленно перечислять все эти работы, так как в каждой из них имеется что-то такое, что не было освещено во всех предыдущих. Однако очевидно, что все вместе они как раз и образуют ту громадную теорию, на которой основано исследование данных математических объектов и их свойств, а также их применение в конкретных областях науки и техники. В данной работе нас, прежде всего, будет интересовать вопрос сохранения линейной системой свойства асимптотической устойчивости по Ляпунову при вариациях параметров данной системы, а именно, вариациях элементов матриц системы и ее запаздываний. Сразу заметим, что в главе 2 нами будут получены явные оценки на вариации этих величин, при выполнении которых данное свойство будет иметь место. Причем, процедуры, которые мы будем использовать, не будут требовать от нас проверки условий Эрмита-Билера (см., например, [12]) или других известных условий, для проверки которых необходимо построение не только самой характеристической функции, соответствующей данной системе уравнений, но и ее производной.

Маленький экскурс в историю начнем с проблемы устойчивости.

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

Ответ на вопрос о непрерывной и непрерывно-дифференцируемой зависимости решений системы дифференциальных уравнений от значений начальных данных (¿(ьжо) и параметров у впервые был дан независимо друг от друга Пикаром (см. Дарбу [25]) и Бендик-соном [22]. Что вместе с теоремой, известной как теорема Гейне-Кантора, о равномерной непрерывности заданной и непрерывной на компактном множестве функции (см., например, [9]), по-видимому, подтолкнуло А. М. Ляпунова, введя естественное обобщение на тот случай, когда некоторое решение системы определено на всей полуоси £ ^ О, сформулировать для него понятия устойчивости и асимптотической устойчивости [11]. Самим Ляпуновым в этой работе было предложено несколько различных способов решения данной проблемы. В частности, один из методов опирается на понятие положительно определенных функций, которые в некотором смысле определяют "меру" или "степень" отклонения возмущенного движения системы от заданного невозмущенного движения той же системы. Теоремы Ляпунова об устойчивости и неустойчивости были обращены В. И. Зубовым [8], и, таким образом, предложенный А. М. Ляпуновым подход является критериальным.

В [11] была доказана и основная качественная теорема теории устойчивости движений — теорема об исследовании устойчивости по линейному приближению.

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

Еще в 1868 году, то есть за 24 года до написания Ляпуновым [11], вышла работа Т. Максвелла [28], в которой исследовались вопросы, связанные с колебаниями в паровых машинах и, наверное, впервые была показана связь нелинейных колебаний данной системы в окрестности положения равновесия с расположением корней характеристического полинома системы линейного приближения.

Здесь сразу необходимо отметить, что все работы, опубликованные до 1892 года, то есть до защиты А. М. Ляпуновым своей докторской диссертации [11], не содержали решения проблемы законности линеаризации, а в силу этого были быстрее похожи на общие рассуждения, чем на строгие математические выводы (конечно, только в той их части, которая относилась к экстраполяции свойств линейной системы на свойства нелинейной системы).

Однако, работа Максвелла поставила важную для теории линейных дифференциальных уравнений проблему о расположении корней характеристического полинома относительно чисто мнимой оси комплексной плоскости. Надо сказать, что вопрос локализации корней полинома в той или иной области комплексной плоскости интересовал в тот момент многих ученых (Ш. Штурм, О. Коши и др.), и поэтому данная проблема не осталась без внимания. Ее решение было получено Э. Раусом [29], который предложил использовать конечную схему, получившую теперь его имя. Данный подход был с пониманием встречен А. Гурвицем, записавшим в [27] решение в виде известной системы детерминантных неравенств, после чего проблема расположения корней произвольной целой функции комплексной переменной левее чисто мнимой оси получила название проблемы Рауса-Гурвица. Придерживаясь терминологии, которая появилась с того знаменательного момента, далее мы будем называть функции, удовлетворяющие этому свойству, устойчивыми.

Известные методы решения проблемы Рауса-Гурвица впоследствии появились благодаря работам в этой области Льенара, Ши-пара, Михайлова, Найквиста и других прославленных математиков.

Вместе с получением всех этих критериев возникла и новая проблема: сохранит ли функция свойство устойчивости при некоторых вариациях параметров, входящих в ее представление. Пожалуй, первой работой по этой тематике можно считать работу И. А. Выш-неградского по исследованию динамики регулятора Уатта [2].

Позднее исследования по построению областей устойчивости в пространстве параметров получили весьма серьезное развитие в работах многих ученых, отметим работы H.H. Меймана и Н. Г. Чеботарева [12], Ю. И. Неймарка [13] и [14].

Однако можно утверждать, что прорыв в этой области наметился лишь с доказательством В. Л. Харитоновым работы [17], в которой была рассмотрена задача устойчивости интервального полинома, а предложенный в этой работе критерий заключался в проверке устойчивости всего четырех полиномов специального вида, принадлежащих этому семейству. Далее появились более общие и не менее интересные результаты, полученные этим автором [18], [19].

В работе [21] коллективом авторов была доказана так называемая Реберная теорема, которая, дополненная результатами работ [5], [24], и [23], дает условия, позволяющие за конечное число шагов проверить устойчивость произвольного семейства полиномов с линейной зависимостью от неопределенных параметров.

Схожие с рассматриваемыми во второй главе нашего повествования проблемы, но для линейной системы обыкновенных дифференциальных уравнений с постоянными коэффициентами, некоторые из компонент матрицы которой являются независимыми интервальными параметрами, были затронуты в работе [6]. Помимо этого, весьма значительное число работ, среди них можно упомянуть [26], [20], [30] и [15], были посвящены созданию и всевозможным модификациям полных алгоритмов проверки устойчивости таких семейств. Отметим, что класс систем, рассматриваемых в данных работах, достаточно узок, а в предыдущих не совсем ясно прослеживается связь между вариациями коэффициентов характеристической функции и вариациями матриц системы, которым эта функция соответствует. Однако ясно, что при практическом применении результатов, когда речь идет о конкретной системе дифференциальных уравнений, именно, эта зависимость и представляет особый интерес. Ответы на все эти и другие схожие вопросы отражены во второй главе данной диссертации, и автор надеется, что полученные им результаты будут подкреплены дальнейшими исследованиями в этой области.

Заметим, что один из основных результатов, изложенных в диссертации, содержится в первой главе и является обобщением известной теоремы Руше (см., например, [16]) на случай квадратных матриц, коэффициентами которых являются аналитические в некоторой области функции комплексного переменного. На основании данного утверждения, в принципе, построено все дальнейшее изложение. В частности, в главе 3 на примере гироскопической системы рассмотрен критический для главы 2 случай и найдены необходимые и, отдельно, достаточные условия асимптотической устойчивости данной системы по Ляпунову при всех достаточно больших значениях параметра, входящего в ее представление, причем, класс систем, которые не подходят под эти условия обозримо узок.

Здесь необходимо отметить, что понятия устойчивости и асимптотической устойчивости по Ляпунову для решений линейной системы уравнений несколько более общего вида, чем рассматриваемые в главах 2, 3 и 4, где h > О — параметр, а элементы п х n-матрицы G(s) имеют на отрезке s G [—/г, 0] ограниченные вариации, можно сформулировать следующим образом: решение (p(t) данной системы, определенное при всех i ) 0 и удовлетворяющее начальным условиям <¿>(0) = хо, ip(s) = (fo(s) при s £ [—h, 0), где Хо — некоторый вектор, a <^o(s) — некоторая заданная и непрерывная на отрезке [—/г, 0] функция, называется устойчивым по Ляпунову, если по любому достаточно малому значению величины £ > 0 найдется Д > 0, для которого неравенство x(t) — tp(t)| < г будет выполняться для всех t ^ 0 и всех решений x(t) данной системы, для которых функция x(s) непрерывна на отрезке s G [—/г, 0] и выполнено неравенство

Если кроме того, для всех таких х(^): — ¿^+оо>0, то решение <£>(£) называется асимптотически устойчивым по Ляпунову.

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

Известно (см., например, [8]), что устойчивость, асимптотическая устойчивость или неустойчивость по Ляпунову любого решения данной системы эквивалентна, соответственно, устойчивости, асимптотической устойчивости или неустойчивости по Ляпунову решения <£>(/) = 0. Таким образом, можно вести речь об устойчивоо

-h х(0) — xq + max x(s) —

В этой же работе В. И. Зубовым доказано, что данная система будет асимптотически устойчива по Ляпунову тогда и только тогда, когда все корни характеристической функции g(A) = сЫ;Х(А), где находятся на множестве Ле А < 0. При этом элементы матрицы Х(А) являются целыми функциями переменной А, а сама g(A), во-первых, в любой полосе вида а ^ ИеА ^ Ь будет иметь конечное число нулей и, во-вторых, все ее корни будут находиться левее некоторой прямой ЯеА = с. Кроме того, существование некоторого корня Ао у функции g(A), для которого Ие Ао > 0, влечет за собой неустойчивость данной системы по Ляпунову.

Всеми этими утверждениями мы будем пользоваться при доказательстве теорем в главах 2, 3 и 4. При этом не составляет особого труда обобщить результаты данных теорем и на случай систем общего вида, указанного выше. Автор не стал этого делать, поскольку, доказательства этих теорем, были бы полным повторением доказательств тех теорем, которые можно найти в этой диссертации.

В главе 2, вместе с прочими результатами, автором был получен алгоритм для поиска величины запаса устойчивости квазиполинома g(A), соответствующего некоторой асимптотически устойчивой по Ляпунову линейной системе дифференциальных уравнений с конечным числом запаздываний, то есть такого числа со > 0, для которого корни квазиполинома g(A) будут находиться левее и на самой прямой ЯеА = —со- При этом доказано, что построенный метод последовательных приближений будет сходиться, именно, к данному числу Со- Автор не стал углубляться в возможности применения о

-к данного алгоритма для решения других задач, но отметим, что, рассматривая практически тот же алгоритм для семейства функций /(А, ср) = сЫ, (—ХЕ + ег

Еще один из результатов второй главы дает новый критерий проверки устойчивости квазиполинома g(A) и может быть применен на практике, во всяком слз^чае, для систем, имеющих не очень большую размерность. К сожалению, автор не рассмотрел в рамках данной работы проблему, связанную с точностью вычислений, необходимой для получения правильных выводов при прохождении алгоритма, который использует данный критерий. Однако, здесь можно заметить, что все функции, рассматриваемые при прохождении алгоритма, будут вещественны и непрерывны по вещественному параметру и, следовательно, на любом отрезке будут являться равномерно непрерывными функциями. Это дает возможность оценить шаг дробления отрезка, на котором минимизируются данные функции, и получить оценку снизу для точной величины их минимума.

Наконец, в последней главе исследуется проблема сохранения свойства асимптотической устойчивости по Ляпунову некоторой линейной системой дифференциально-разностных уравнений с конечным числом запаздываний, когда данные запаздывания могут варьироваться произвольным образом в пределах множества (0, +оо). Автором доказана теорема, согласно которой решение данного вопроса сводится к исследованию устойчивости одного семейства полиномов специфического вида. Данная задача вырастает из результатов главы 3, в которой от системы потребовалось выполнение подобного свойства в одном из трех частных случаев исследования.

Список используемых обозначений.

Если А — некоторая матрица размерности п х т, то |(А|| — максимальный из модулей элементов А; А* — комплексно сопряженная матрица к матрице А; БрА для квадратной А — сумма элементов главной диагонали. Е — диагональная матрица, с единицами на диагонали.

Тп = \/пР.

Сп — евклидовы пространства размерности п над полем, соответственно, вещественных и комплексных чисел.

Если О — некоторое множество из В? или Сп, то В — его замыкание.

В Кп или Сп: Ке(го), Се(го) обозначают, соответственно, множества {г | \\г - 20|| < г}, {г | \\г - г0\\ = е}, где ||г|| = х = г*, если г £ С1 и неоговорено противное.

С+ — множество точек А ЕС, для которых ИеЛ ^ 0.

Знак ■ по тексту — конец доказательства.

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

Заключение диссертации по теме «Дискретная математика и математическая кибернетика», Купцов, Сергей Юрьевич

Заключение.

После проведения данных исследований хотелось бы подвести некоторую черту под тем, что сделано, и попытаться указать те утверждения, в которых результат может быть улучшен, причем, возможно в значительной степени. Прежде всего это относится к формулировке условий теоремы 1.1. А именно, в ходе доказательства данного утверждения мы использовали оценки, в которых равенство, вообще говоря, может и не достигаться. Таким образом, остается открытым вопрос, существуют ли в пространстве матриц размерности п х п матрицы А и В, для которых: ф 0, с1е1 (А + В) = О, а в условии (1.1) стоит знак равенства. В силу этого замечания все последующие оценки, при получении которых автор ссылается на теорему 1.1 или другие утверждения первой главы, возможно могут быть улучшены даже в общем случае. Однако, принимая во внимание алгоритмы, изложенные во второй главе, данное обстоятельство всегда может быть преодолено.

Заметим, что все результаты второй главы явно или неявно используют результат теоремы 2.1, которая является следствием непрерывной зависимости корней Л уравнения g(A, К) = 0 от координат вектора параметров к при выполнении условий теоремы. Автор привел доказательство данного утверждения, хотя и сомневается, что данная проблема не рассматривалась ранее другими авторами. Здесь можно вспомнить, например, так называемый принцип исключения нуля из [20] для полиномов. Вместе с тем, результаты следствия 2.1 и утверждений первой главы позволяют найти явные оценки на вариации параметров линейной системы дифференциально-разностных уравнений, при выполнении которых система будет сохранять свойство асимптотической устойчивости по Ляпунову, а также дать качественно новый критерий устойчивости ее характеристического квазиполинома и построить метод последовательного приближения для поиска запаса устойчивости данной системы.

В третьей главе автором проведено исследование вопроса асимптотической устойчивости гироскопической системы, при условии существования у системы некоторого ненулевого запаздывания. Заметим, что результаты теорем 3.1 и 3.2 без особого труда могут быть обобщены и на тот случай, когда у данной системы определено не одно, а сразу несколько положительных запаздываний h(H) — (hi(H),., hm(H)), для которых if-1||/i(.ff)|| ^ а0 при Н ^ 1, где «о — достаточно малое число, а сама система записана в виде

Ax(t) + HBx(t) + Cx(t) = D\x(t - hi) + • • • + Dmx(t - hm).

Наконец, в последней главе можно ввести естественные обобщения на те случаи, когда запаздывания системы (4.1) можно разбить на группы соизмеримых запаздываний: h = (hi,., hm) = {aikn,. .,aikPlU.,atqhq,. .,aqkPqq), где ol\, . ,aq — параметры, значения которых могут изменяться в пределах (0,+оо), a kps — заданные натуральные константы. Данная задача может возникнуть при приведении некоторой системы уравнений дифференциально-разностных уравнений /-ого порядка к аналогичной системе уравнений первого порядка. Однако заметим, что в данной главе особый интерес представляет, именно тот факт, что вопрос об абсолютной устойчивости некоторого квазиполинома по т положительным запаздываниям может быть сведен к аналогичной задаче для известного m-параметрического семейства полиномов с неотрицательными параметрами.

Список литературы диссертационного исследования кандидат физико-математических наук Купцов, Сергей Юрьевич, 1998 год

1. Беллман Р. Введение в теорию матриц. — М.: Наука, 1969. — 368 с.

2. Вышнеградский И.А. О регуляторах прямого действия // Известия Санкт-Петербургского технологического института. — 1877. — Т.1. — С. 21-62.

3. Гантмахер Ф.Р. Теория матриц. — М.: Наука, 1967. — 576 с.

4. Жабко А.П., Харитонов В.Л. Методы линейной алгебры в задачах управления. — СПб.: Изд-во С.-Петерб. гос. ун-та, 1993.320 с.

5. Жабко А.П., Харитонов В. Л. Необходимые и достаточные условия устойчивости линейного семейства полиномов // Автоматика и телемеханика. — 1994. — № 10. — С. 125-134.

6. Заде Л. А., Дезоер Л.А. Теория линейных систем: Метод пространства состояний. — М.: Наука, 1970. — 703 с.

7. Зубов В.И. Аналитическая динамика гироскопических систем.

8. Л.: Судостроение, 1970. — 318 с.

9. Зубов В.И. Лекции по теории управления. — М.: Наука, 1975.496 с.

10. Кудрявцев Л.Д. Курс математического анализа. Т.1. — М.: Высш. шк. — 1988. — 712 с.

11. Купцов С.Ю. Об абсолютной устойчивости линейной системы дифференциально-разностных уравнений // Труды XXIX науч. конф. — С.-Петерб. гос. ун-т. Фак-т прикл. мат. проц. управл.1998. — С. 68-69.

12. Ляпунов A.M. Общая задача об устойчивости движения. — М.: ГИТТЛ, 1950. — 472 с.

13. Мейман H.H., Чеботарев Н.Г. Проблема Рауса-Гурвица для полиномов и целых функций. — М.;Л.: Изд-во АН СССР, 1949. — 332 с.

14. Неймарк Ю.И. Устойчивость линеаризованных систем (дискретных и распределенных). — Л.: Изд-во Ленингр. Военно-воздушной Инж. Акад., 1949. — 140 с.

15. Неймарк Ю.И. Структура D-разбиения пространства полиномов и диаграммы Вышнеградского // Докл. АН СССР. — 1948.1. Т.52. — С. 853-856.

16. Поляк Б.Т. Цыпкин Я.З. Частотные критерии робастной устойчивости и апериодичности линейных систем // Автоматика и телемеханика. — 1990. — № 9. — С.45-54.

17. Свешников А.Г., Тихонов А.Г. Теория функций комплексной переменной. — М.: Наука, 1967. — 304 с.

18. Харитонов В.Л. Асимптотическая устойчивость семейства систем линейных дифференциальных уравнений // Дифференц. уравнения. — 1978. — Т.14. — С. 2086-2088.

19. Харитонов В. Л. Проблема Рауса-Гурвица для семейств полиномов и квазиполиномов // Мат. физика. — Киев, 1979. — Вып.26.1. С. 79-89.

20. Харитонов В. Jl. Задача распределения корней характеристического полинома автономной системы // Автоматика и телемеханика. — 1990. — № 3. — С. 42-47.

21. Ackermann J. Does it suffices to check a subset of multilinear parameters in robustness analysis? // IEEE Trans, on Automatic Control. — 1992. — Vol. 37. — P. 487-488.

22. Bartlett A.C., Hollot C.V., Huang Lin. Roots location of an entire polytope of polynomials: It suffices to check the edges // Mathematics of Control, Signals and Systems. — 1988. — Vol. 1. — P. 61-71.

23. Bendixson I. Démonstration de l'existence de l'integrale d'une équation aux derivèes partielles linéaire //Bull. Soc. Math. France. — 1896. — № 24, Vol.3. — P. 220-225.

24. Bialas S. A necessary and sufficient condition for the stability of convex combinations of stable polynomials and matrices // Bull, of Polish Academy of Sciences. — 1985. — Vol. 33. — P. 473-480.

25. Bose N.K. Test of Hurwitz and Schur properties of convex combination of complex polynomials // IEEE Trans, on Automatic Control — 1989. — Vol. 36. — P. 1245-1247.

26. Darboux G. Leçons sur la théorie générale des surfaces. Vol. 5. — Paris: HfGauthier-Villars, 1896. — 363 p.

27. DeGaston R.R.E., Safonov M.G. Exact calculation of the multiloop stability margin // IEEE Trans, on Automatic Control. — 1988. — Vol. 33. — P. 156-171.

28. Hurwitz A. Uber die bedingungen, unter welchen eine gleichung nur wurzeln mit negativen reellen teilen besitzt // Math. Ann. — 1895. — Bd. 46. — S. 273-284.

29. Maxwell J.C. On governers // The Proc. of the Royal Society. — 1868. — № 100. — P. 105-120.

30. Routh E.J. A treatise on the stability of a given state of motion. — London: Macmillan Publ. Co., 1877. — 150 p.

31. Tesi A., Vicino A. Computation of nonconcervative stability perturbation bounds for systems with nonlinearly correlated uncertainties // IEEE Trans, on Automatic Control. — 1990. — Vol. 35. — P. 835-841.