WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Санкт-Петербургский государственный университет

На правах рукописи

Шумафов Магомет Мишаустович

СТАБИЛИЗАЦИЯ УПРАВЛЯЕМЫХ ДИНАМИЧЕСКИХ СИСТЕМ

Специальность 05.13.01 – системный анализ, управление и обработка информации (по прикладной математике и процессам управления)

АВТОРЕФЕРАТ

диссертации на соискание ученой степени доктора физико-математических наук

Санкт-Петербург 2012

Работа выполнена на математико-механическом факультете СанктПетербургского государственного университета (СПбГУ).

Научный консультант:

Доктор физико-математических наук, член-корреспондент РАН Леонов Геннадий Алексеевич (Санкт-Петербургский государственный университет)

Официальные оппоненты:

Доктор технических наук, академик РАН Микрин Евгений Анатольевич (Первый заместитель генерального конструктора ОАО РКК «Энергия» им. С.П. Королева) Доктор физико-математических наук, профессор Харитонов Владимир Леонидович (Профессор кафедры теории управления факультета Прикладной математики – процессов управления Санкт-Петербургского государственного университета) Доктор технических наук, член-корреспондент РАН Юсупов Рафаэль Мидхатович (Директор Научно-исследовательского института информатики и автоматизации РАН)

Ведущая организация:

Институт проблем управления им. В.А. Трапезникова РАН

Защита состоится 26 декабря в ____ часов на заседании диссертационного совета Д.212.232.50 по защите диссертаций на соискание ученой степени кандидата наук, на соискание ученой степени доктора наук при СанктПетербургском государственном университете по адресу: 199034, СанктПетербург, В.О., Университетская набережная 7/9. Менделеевский центр.

С диссертацией можно ознакомиться в библиотеке им. М. Горького СанктПетербургского государственного университета.

Автореферат разослан «_____»________________2012г.

Ученый секретарь диссертационного совета доктор физ.-мат. наук, профессор Курбатова И.Г.

АННОТАЦИЯ Диссертационная работа посвящена исследованию вопросов устойчивости и стабилизации динамических систем, определяемых дифференциальными уравнениями. В ней разработан новый, простой, алгоритм стабилизации по состоянию линейных управляемых систем. Доказаны теоремы, дающие эффективно проверяемые, коэффициентные, критерии стабилизируемости линейных управляемых систем обратной связью с запаздыванием.

Получены эффективные частотные критерии устойчивости и стабилизируемости дифференциальных систем с гистерезисными нелинейностями. Исследована задача о переходном процессе и стабилизации динамической системы «машина — регулятор Уатта».

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. Одной из важнейших задач теории управления является задача о стабилизации. Интерес к проблемам стабилизации мотивируется как запросами практики управления, так и формулировками открытых проблем многими известными учеными: Н.Н. Красовским, В.И. Зубовым, У.М. Уонэмом, Д.С. Бернстейном, Р.У.

Брокеттом, Дж. Розенталем, Дж. Виллемсом, В. Блонделем, Е. Сонтагом, М. Видясагаром.

Одной из классических задач стабилизации систем была задача о стационарной стабилизации по состоянию линейной стационарной управляемой системы, т.е. поиск соответствующего линейного стационарного регулятора, обеспечивающего асимптотическую устойчивость замкнутой системы. Эта задача, и более общая задача о размещении собственных чисел матрицы (или полюсов передаточной функции), впервые была сформулирована и решена в работах В.И. Зубова и У.М. Уонэма. После выхода этих работ было написано большое количество статей, где предлагались другие, альтернативные к данным Зубовым и Уонэмом алгоритмы стабилизации. Однако все они достаточно громоздки для вещественного и векторного случая. Предложенные в диссертации алгоритмы стабилизации являются новыми и, на наш взгляд, наиболее простыми.

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

Красовского, В.И. Зубова, Е.А. Микрина, Р.М. Юсупова, Г.А. Леонова, Б.Т. Поляка, П.С. Щербакова, В.Л. Харитонова, В.С. Антончика, Е.Я. Смирнова, А.И. Кирьянена, А.Г.

Александрова, А.В. Назина, М.В.Хлебникова, В.Н.Честнова, И.Я. Каца, Ю.С. Осипова, Е.Л.

Тонкова, В.А. Зайцева, Э.Г. Альбрехта, Г.С. Шелементьева, С.И. Солодушкина, Ю.Ф.

Долгий, И.В. Гайшуна, Д.С. Бернстейна, В. Блонделя, В.Л. Сирмоса, С.Т. Абдаллаха, П. Дорато, К. Григориадиса. Наиболее эффективные методы и алгоритмы стабилизации разработаны для линейных систем управления. Одной из проблем, стимулировавшей немало публикаций, была сформулированная Р. Брокеттом проблема стабилизации линейной стационарной системы с помощью линейной нестационарной обратной связи. Решение проблемы Брокетта в ряде важных для практики случаев дано в работах Г.А. Леонова и Л. Моро, Д. Аэлса. В этих работах построены алгоритмы низкочастотной и высокочастотной стабилизации линейной стационарной системы. Для двумерных и трехмерных стационарных систем показано, как введение нестационарной обратной связи расширяет возможности стационарной стабилизации. Возникает естественный вопрос: существуют ли иные (кроме нестационарных) способы стабилизации линейных стационарных систем, позволяющие расширить возможности стационарной стабилизации? Приходим к основной задаче в постановке Г.А. Леонова: можно ли ограничиваясь лишь линейной стационарной обратной связью стабилизировать линейную стационарную систему введением запаздывания в обратную связь? Каковы возможности линейной стационарной обратной связи с запаздыванием для стабилизации неустойчивых линейных стационарных систем? Хорошо известно, что для некоторых классов уравнений при достаточно малых и достаточно больших запаздываний такая стабилизация невозможна. С другой стороны, мотивацией к исследованию стабилизации путем введения запаздывания в обратную связь явились компьютерные эксперименты К. Пирагоса по стабилизации хаоса – стабилизации неустойчивых периодических орбит, погруженных в странный аттрактор той или иной хаотической системы. В диссертационной работе решена поставленная выше Г.А. Леоновым задача для случая двумерных и трехмерных управляемых систем. В частности, сделан вывод о том, что для обеспечения эффекта Пирагоса необходимо ввести зависящий от времени коэффициент усиления в обратной связи.

Далее, в нелинейной теории колебаний и в теории управления автогенераторами актуальными являются проблемы захватывания частоты автоколебаний внешним гармоническим воздействием. Проблемам захвата различных автоколебательных режимов под частоту внешнего воздействия посвящено немалое количество работ. Одним из важнейших свойств вынужденных периодических процессов является их устойчивость «в целом», когда явление захватывания наблюдается в любом режиме работы автогенератора.

Кроме того имеется ряд экспериментальных результатов, показывающих, что различные нелинейные системы, допускающие хаотическое поведение, могут быть стабилизируемы гармоническим или другим периодическим внешним воздействием. В связи с этим весьма актуальной является задача о стабилизации и исследование устойчивости дифференциальных систем с гистерезисными функциями. Последние являются математическими моделями систем автоматического управления, в которых возникает люфт, сухое трение, некоторые приборы имеют зоны нечувствительности, происходит пространственное запаздывание управляющего сигнала. Первые теоретические исследования систем с гистерезисными нелинейностями появились в 40-е годы в работах А.А. Андронова, А.А. Фельдбаума и Ф. Краутвига. Ряд результатов по исследованию двумерных систем с гистерезисом был получен также в более поздних работах Р.А. Нелепина, В.В. Казакевича, В.В. Петрова, Г.М. Уланова и др. В этих работах использовались методы фазовой плоскости и точечных отображений, а гистерезисные функции имели «стандартный» вид. Дальнейшее развитие техники требовало создания математических методов глобального исследования многомерных систем с гистерезисом. Такие методы появились в начале 60-х годов вслед за выходом работ В.М. Попова. В последующие годы в работах В.А. Якубовича, Н.Е. Барабанова, А.Х. Гелига, В.А. Брусина, Я.З. Цыпкина и др. были получены эффективно проверяемые частотные критерии различных типов поведения решений систем с гистерезисными функциями. Новые методы, развитые в теории абсолютной устойчивости, позволили исследовать системы с гистерезисными нелинейностями, удовлетворяющими достаточно общим ограничениям. В работах В.А. Якубовича, Н.Е. Барабанова, А.Х. Гелига, М.А. Красносельского, А.В. Покровского были даны строгие определения понятия гистерезисной функции, позволяющие обсуждать в общей постановке проблемы существования, единственности, продолжимости решений систем дифференциальных (и интегральных) уравнений с гистерезисом. В диссертационной работе, развивая идеи и методы разработанные В.А. Якубовичем и его учениками для исследования абсолютной устойчивости нелинейных систем, проводится исследование вопросов устойчивости и стабилизации дифференциальных систем с гистерезисными нелинейностями. Далее, в инженерной практике при проектировании различных энергоустановок (например, турбогенераторов) и технических систем, а также анализе и синтезе систем управления такими объектами, возникает задача о проведении нелокального анализа переходного процесса — от запуска машины в начальный момент времени до её выхода в рабочий режим.

Математическая постановка этой задачи дана в одной из работ Г.А. Леонова. Формализация указанной выше задачи укладывается при некоторой идеализации в математическую модель работы динамической системы «машина–регулятор Уатта» (точнее, модифицированный регулятор Уатта). В работе Г.А. Леонова проведен нелокальный анализ системы «машина– регулятор Уатта» в предположении, когда в уравнениях движения угловое ускорение и коэффициент трения являются постоянными, восстанавливающая сила — линейной, а центробежная сила имеет приближенный вид. В диссертационной работе поставленная Леоновым задача решена в общем виде без указанных выше предположений.

Цель диссертационной работы. 1) Разработка нового, простого алгоритма стационарной стабилизации по состоянию линейной стационарной управляемой системы.

2) Выяснение возможностей линейной обратной связи с запаздыванием для стабилизации неустойчивых линейных стационарных управляемых систем.

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

4) Получение эффективных частотных критериев устойчивости и стабилизации систем с гистерезисом.

5) Проведение нелокального анализа переходного процесса в динамической системе «машина–регулятор Уатта».

Задачи исследования. В соответствии с целью в работе поставлены и решены следующее задачи:

– построение нового, простого, алгоритма для стабилизации по состоянию линейных стационарных управляемых систем;

– получение необходимых и достаточных условий стабилизируемости двумерных и трехмерных линейных систем обратной связью с запаздыванием;

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

– проведение нелокального анализа переходного процесса при включении динамической системе «машина–регулятор Уатта».

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

Научные результаты. На защиту выносятся следующие научные результаты:

1) Новый алгоритм стабилизации по состоянию линейных управляемых систем.

2) Критерии стабилизируемости двумерных и трехмерных линейных управляемых систем обратной связью с запаздыванием.

3) Эффективный частотный критерий стабилизации внешним гармоническим воздействием автономных нелинейных систем, допускающих, в частности, хаотическое поведение.

4) Новые частотные критерии устойчивости дифференциальных систем с гистерезисными нелинейностями.

5) Достаточные условия асимптотической устойчивости динамической системы «машина–регулятор Уатта» с оценкой снизу области притяжения в фазовом пространстве системы.

Научная новизна диссертационной работы заключается в том, что: 1) Разработан новый, простой, алгоритм стабилизации по состоянию линейных стационарных управляемых систем.

2) Получены новые необходимые и достаточные условия стабилизируемости двумерных и трехмерных линейных управляемых систем обратной связью с запаздыванием двух видов: обычной и по Пирагосу.

3) Получен новый частотный критерий стабилизации при помощи внешнего гармонического воздействия автономных нелинейных систем, допускающих, в частности, хаотическое поведение.

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

5) Получено новое достаточное условие асимптотической устойчивости динамической системы «машина–регулятор Уатта».

Достоверность и обоснованность научных результатов. Результаты, полученные в диссертационной работе достоверны и обоснованны, они математически строго доказаны.

Практическая ценность работы. Диссертация носит теоретический характер. Ее результаты могут быть применены:

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

2) В различных технических системах с гистерезисом, при проектировании энергоустановок, в частности, турбогенераторов.

3) В учебном процессе при изучении курсов «Дифференциальные уравнения», «Математическое моделирование», а также при написании курсовых и дипломных работ.

Реализация результатов работы носит теоретический характер. Результаты исследований нашли свое применение в работах по выполнению Программы фундаментальных исследований Президиума РАН №19 «Управление механическими системами (проект № 1.4), по гранту РФФИ (проект № 04-01-00-250А), по программе «Университеты России», по гранту НШ-2257.2003.1 Совета по грантам Президента РФ для поддержки молодых российских ученых и ведущих научных школ, в лаборатории «Интел» при Санкт-Петербургском государственном университете и в учебном процессе Адыгейского государственного университета.

Вклад автора в разработку проблемы. Все основные результаты диссертационной работы получены автором самостоятельно.

Апробация работы. Основные положения диссертации докладывались на научных конференциях, в том числе международных, и обсуждались на научно-исследовательских семинарах. В частности, – на научно-исследовательском семинаре по теории управления при кафедре теоретической кибернетики в Санкт-Петербургском государственном университете (руководитель чл.-корр. РАН, профессор Г.А. Леонов);

– на научно-исследовательском семинаре по динамическим системам и теории управления в Адыгейском государственном университете (г. Майкоп; руководители:

профессор К.С. Мамий, доцент М.М. Шумафов);

– на первой научно-практической конференции «Дифференциальные уравнения и их применения» (Санкт-Петербург, декабрь, 1996);

– на международной конференции «Нелокальные краевые задачи и родственные проблемы математической биологии, информатики и физики» (Нальчик, декабрь, 1996);

– на международной конференции «Control of Oscillations and Chaos» (St. Petersburg, August, 1997);

– на третьей международной конференции «Дифференциальные уравнения и их применения» (Саранск, май, 1998);

– на международной конференции по прикладной и вычислительной математике («International Conference on Applied and Computational Mathematics»);

– на XII международной конференции «Устойчивость и колебания нелинейных систем управления» (конференция Пятницкого) (Москва, июнь, 2012).

Публикации. По теме диссертационной работы опубликовано 24 работы, в том числе тринадцать статей в журналах, рекомендованных ВАК РФ, четыре монографии.

Структура и объем работы. Диссертация состоит из введения, семи глав, заключения, библиографического списка и содержит 260 станиц текста, в том числе рисунков.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Якубовича, Н.Е. Барабанова и др. Гистерезисные функции, входящие в рассматриваемые системы, могут содержать в отличие от рассматривавшихся в литературе случаев несколько петель, которые обходятся в различных направлениях. На конкретном примере проведен сравнительный анализ полученных результатов с ранее известными. Показано, что в пространстве параметров полученное нами частотное условие выделяет большую область устойчивости, чем известное ранее частотное условие. Сформулируем соответствующие результаты.

Ниже всюду будем придерживаться следующего формального определения гистерезисной функции (В.А. Якубович). Пусть для каждого значения (-, + ) задано множество E[ ] «начальных значений» гистерезисной функции. Пусть далее любому t 0, каждой непрерывной на отрезке [0,t] функции ( ) и любому значению0 E[ (0)] сопоставлено число [,0], причем [,0] = 0. В этом случае говорят, что при t t задана гистерезисная функция [,0].

t Рассмотрим систему x = Ax + b, = c*x, (1) (t) = [,0 ]t, где A - постоянная (n n) - матрица, b и c - постоянные n -векторы, [,0]t - гистерезисная функция.

Под решением системы (1) понимается пара функций (x(t), (t)), t [t0,T ], где x(t) - абсолютно непрерывна, а (t) - суммируема по Лебегу функция удовлетворяющая уравнению (1) при почти всех t [t0,T ].

Введем передаточную функцию линейной части системы (1) () = с*(A - E)-1b.

Предположим, что гистерезисная функция удовлетворяет соотношению («условие секториальности») 0 (t)[,0 ]t (t)2 (0 < ). (2) Определение 1. Будем говорить, что система (1) обладает глобальной асимптотикой (глобально асимптотически устойчива), если любое ее решение стремится при t + к некоторому состоянию равновесия.

Определение 2. Система (1) называется абсолютно устойчивой, если 1) для любого решения x(t) системы (1) lim x(t) = 0, и 2) существует непрерывная возрастающая t+ функция (r), зависящая только от A, b, c и , (0) = 0, такая, что выполнено неравенство x(t) (| x(0) |) для любого 0 E( ), = c*x(0).

0 Определение 3. Стационарным множеством системы (1) назовем множество = {(x0,0 ) : Ax0 + b0 = 0, [c*x0,0 ]t 0}.

Определение 4. Система (1) дихотомична, если любое ограниченное при t решение стремится при t + к стационарному множеству.

Определение 5. Гистерезисная функция [,0]t обладает свойством предельной непрерывности, если из соотношений (t) и [,0 ]t * при t + следует, * что * E[ ] и [*,*]t *. (Здесь E[* ] R - множество, фигурирующее в * определении гистерезисной функции) Теорема 1. (О глобальной асимптотике.) Предположим, что в системе (1) A - гурвицева матрица, гистерезисная функция [,0]t удовлетворяет соотношению (18) и при = ограничена, а передаточная функция () невырождена. Пусть далее выполнены следующие условия:

1) система (1) диссипативна, 2) существуют числа > 0, > 0, 0 и такие, что выполнено неравенство / + Re( +i)(i) - i(i) 0, (3) 3) существует непрерывная функция F( ) и число , для которых выполнено неравенство [ (t),0 ]t - F[ (t)] [,0 ]t, (4) 4) 4 > ( )2.

Тогда для любого решения x(t) системы (1) lim x(t) = 0. Если, кроме того, гистерезисная t+ функция удовлетворяет свойству предельной непрерывности, то x(t) 0 является стационарным решением системы (1), и тогда система (1) обладает глобальной асимптотикой.

Нижеследующая теорема не предполагает диссипативности системы (1) и усиливает теорему 1.

Теорема 2. (Об абсолютной устойчивости.) Предположим, что выполнены все условия теоремы 1 кроме условия 1), а условие 3) формулируется так: существует непрерывная функция F( ), график которой принадлежит сектору {(, ) :0 / }, удовлетворяющая неравенству (4). Тогда утверждение теоремы 1 остается в силе. Более того, система (1) абсолютно устойчива.

Замечание. Утверждения теорем 1 и 2 остаются справедливыми, если частотное неравенство (3) заменить на частотное условие:

/ + Re( + i)(i) )( + ) (i) ( 0), которое несколько усиливает теоремы 1 и 2.

Сформулируем теперь теорему дихотомичности системы (1) с гистерезисной функцией [,0]t в другом классе функций. А именно, пусть функция [,0]t удовлетворяет следующим условиям:

1) если функция (t) абсолютно непрерывна на [0, + ), то и [,0 ]t – тоже;

2) при почти всех t [0, + ) выполнено соотношение d[,0 ]t d (t) d (t) 0 ( > 0), (5) dt dt dt причем если (t) 0, то d[0,0 ]t dt 0 ;

3) функция [,0 ]t обладает свойством предельной непрерывности;

4) если (x(t), (t)) - решение системы (1), то пара функций (x(t + c), (t + c)) - тоже (c = const).

Примером гистерезисной функции удовлетворяющей условиям 1)-4) является, например, люфт.

Теорема 3. (О дихотомичности.) Пусть в системе (1) A гурвицева матрица, гистерезисная функция [,0 ]t удовлетворяет вышеприведенным условиям 1)-4). Пусть далее выполнены условия 1)-4) теоремы 1. Тогда система (1) дихотомична.

Пример 1. Рассмотрим уравнение второго порядка с гистерезисной нелинейностью x + x + = 0, x + (6) = [x,0 ]t, где > 0, > 0 - постоянные, а = [x,0 ]t - гистерезисная функция. (Здесь (t) x(t).) Предложение 1. Если выполнено неравенство 2 0 < <, 2( + ) - то система (6) обладает глобальной асимптотикой в классе всех гистерезисных нелинейностей [x,0 ]t, удовлетворяющих «условию секториальности» (2) и обладающих свойством предельной непрерывности. (Здесь > 0 - число, фигурирующее в неравенстве (4).) На примере системы (6) проведен сравнительный анализ результатов, даваемых частотным условием (3) теоремы 1 с одной стороны, и известным частотным условием. Область глобальной асимптотики, выделяемая условием (3) в пространстве параметров шире, чем область, выделяемая известным частотным условием.

Пример 2. Рассмотрим систему второго порядка с гистерезисом x = -x - y, (7) y = x - , = [ y,0 ]t, где > 0, > 0 - постоянные, [y,0]t - гистерезисная функция.

Предложение 2. Пусть гистерезисная функция [y,0]t удовлетворяет «условию секториальности» (2) и обладает свойством предельной непрерывности. Тогда система (7) обладает свойством глобальной асимптотики, если 2 - 3 2 <, >.

( -1)2 ( - число, фигурирующее в неравенстве (4).) В третьей главе рассматриваются те же самые вопросы, что и во второй главе, в критическом случае, когда одно из собственных чисел матрицы линейной части системы рано нулю. Здесь получены соответствующие частотные критерии устойчивости в критическом случае. Рассматриваются два примера, соответствующие критическому случаю.

Невырожденным преобразованием x = Sz, z = (y,)* (det S 0) система (1) приводится к виду y = Py - q, = r* y + , (8) = -, = [,0 ]t, n-1 n-где y R, R, P - гурвицева (n -1) (n -1) - матрица, q, r R, и - числа.

Поскольку передаточные функции () и 1() систем (1) и (8) совпадают, то () = c*(A - E)-1b + / .

Положим = r*q - . Имеет место следующая теорема.

Теорема 4. (О глобальной асимптотике.) Предположим, что матрица A в системе (1) имеет одно простое нулевое собственное значение, а остальные собственные значения имеют отрицательные вещественные части. Пусть гистерезисная функция [,0 ]t удовлетворяет «условию секториальности» (2) и обладает свойством предельной непрерывности.

Пусть далее передаточная функция () системы (1) невырождена и выполнены следующие условия:

1) существуют числа > 0, > 0, 0 и такие, что выполнено частотное неравенство / + + 2 Re[ / 2 + ( + / 2)i](i) > + + i(i) 0 ;

2) существует непрерывная функция F( ), для которой при некотором > выполнено неравенство [,0 ]t - F( (t)) [,0 ]t, ± причем F( )d < + ;

3) 4 > ( )2.

Тогда система (1) обладает глобальной асимптотикой.

Теорема 5. (О дихотомичности.) Пусть выполнены все условия теоремы 4, кроме ± условия сходимости интеграла F( )d. Тогда система (1) дихотомична.

Рассмотрим теперь случай, когда гистерезисная функция [,0 ]t не удовлетворяет «условию секториальности» (2), но удовлетворяет условию (5).

Теорема 6. (О глобальной асимптотике.) Пусть в системе (8) матрица P - гурвицева, гистерезисная функция [,0 ]t ограничена, удовлетворяет условию (5), обладает свойством предельной непрерывности и такова, что существование предела lim [,0 ]t = 0 влечет существование конечного предела lim (t) = . Пусть далее пара t+ t+ (P, q) управляема и выполнены следующие условия:

a) существуют числа > 0, > 0, 0 и такие, что справедливо частотное неравенство -1 Re K(i) - K(i) - ReiK(i) + R, где K() = r*(P - E)-1q - , = r*q -.

b) выполнены условия 2) и 3) теоремы 4.

Тогда система (1) обладает глобальной асимптотикой.

Примеры. 1.Рассмотрим уравнение x + = 0, = [x,0 ]t, x + (9) где > 0 - постоянное число, а = [x,0 ]t - гистерезисная функция. Применим к системе (9) теорему 4.

Предложение 3. Система (9) обладает глобальной асимптотикой в классе всех гистерезисных функций [x,0 ]t, удовлетворяющих «условию секториальности» (2), если выполнено неравенство > , где (0,1).

Применяя к системе (9) теорему 6, получаем следующее Предложение 4. Система (9) обладает глобальной асимптотикой при всех > 0, если гистерезисная функция [,0 ]t удовлетворяет всем условиям, приведенным в теореме 6.

2. Рассмотрим систему y = - y /T + q, = ry + , (10) = -, = [,0 ]t, где y, , R ; q, r, и - числа такие, что r * q = -1, = T - (1- ), T > 0, 0 < < 1.

Применяя к системе (10) теорему 4 и 6, получаем следующие утверждения.

Предложение 5. Пусть гистерезисная функция [,0 ]t удовлетворяет «условию секториальности» (2). Тогда система (10) обладает глобальной асимптотикой, если выполнены неравенства 0 0 T > 1- , 0 < < 2 /, где = max{1, 2, 3}, 1 = (1-1/ - T ) T, 2 2 T (1- )2 + 2(1+ T ) + T (1- ) - 2 / T (1- ) + T -1/ +2 =, 3 =.

1+ T T Предложение 6. Пусть гистерезисная функция [,0 ]t удовлетворяет всем условиям, приведенным в теореме 6. Тогда система (10) обладает глобальной асимптотикой, если выполнено неравенство T > 1- .

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

Андронова, А.А. Витта, В.А. Плисса, В. Линдсея и др.). Одним из важнейших свойств вынужденных периодических процессов является их устойчивость «в целом», когда явления захватывания наблюдается в любом режиме работы автогенератора. Кроме того имеется ряд экспериментальных результатов, показывающих, что различные нелинейные системы, допускающие хаотическое поведение, могут быть стабилизируемы гармоническим внешним воздействием (работы Е.Н. Дудника, Ю.И. Кузнецова, И.И. Минакова, Ю.М. Романовского, П.С. Ланда, А.Ф. Ольховой, С.М. Перминова, В.В. Мигулина, Б.А. Сильнова и др.).

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

Рассматриваемая система имеет вид x(t) = Ax(t) + b( (t) + sin t), (t) = c*x(t), (11) (t) [,0], t n где x(t) R - вектор состояния, (t) R - вход, (t) R -выход, [,0] - гистерезисная t n функция, A - вещественная постоянная (n n) - матрица, b,c R - постоянные векторы, , - положительные числа. (Знак * означает операцию транспонирования.) Пусть гистерезисная функция [,0] удовлетворяет следующим условиям:

t 1) существуют константы 1, R такие, что d[,0 ]t d d 0 (12) dt dt dt для почти всех t, для которых (t) [1, ];

2) справедливо неравенство d[,0 ]t d d k (13) dt dt dt для почти всех t ;

3) функция [,0] ограничена t [,0 ]t l (14) для всех t и (t). Здесь , k,l - некоторые положительные числа.

Введем передаточную функцию системы (11) ( p) = c * (A - pE)-1b и следующие обозначения:

+ (t) = c * eAtb, = (t) dt, 1 + l 1 - l = lim p( p) = -c *b ; T =.

arcsin 2 - arcsin p (i ) (i ) Здесь предполагается, что 1 - l < (i ), + l < (i ) (15) Будем предполагать, что матрица A гурвицева. Тогда ясно, что < +.

При любых начальных данных t0, x0,0 решение системы (11) существует и продолжимо на полупрямую [t0, + ]. Число T является оценкой сверху наибольшего возможного времени пребывания решения x(t) системы (11) в полосе {x |1 < c * x < } при достаточно больших t. Сформулируем основную теорему.

Теорема 7. Пусть в системе (11) гистерезисная функция [,0] удовлетворяет t условиям (12)-(14). Если при некоторых числах 1 > 0,2 > 0 и > 0 выполнены следующие условия:

1) все полюсы передаточной функции ( p - ) имеют отрицательные вещественные части;

2) при всех 0 выполнено неравенство 2 1/ + Re (i - ) - 1 (i) - - 2 (i - )(i - ) - 0 ;

3) выполнены неравенства:

* | * - l | | + l | < 1, < 1, (i ) > l , (i ) (i ) * * где и соответственно левый и правый концы интервала (*, ) оси , на которую * проецируется петля графика гистерезисной функции [,0] ;

t 4) выполнено неравенство -k(1+ k ) - 2 > T, 1то для любых двух решений x1(t) и x2 (t) системы (11) справедливо предельное соотношение lim x1(t) - x2 (t) = 0. (16) t+ 2 Более того, существует - периодическое решение x0 (t) системы (11), к которой стремятся при t + все другие решения системы (11).

Последнее утверждение теоремы 7 адекватно захватыванию под частоту внешнего гармонического воздействия. Очевидно, что при достаточно малых и достаточно больших для некоторых 1 > 0,2 > 0, > 0 выполнены все условия теоремы 7.

Следовательно, при малых внешнее гармоническое воздействие с достаточно большой амплитудой стабилизирует систему (11).

Для случая = 0 соотношение (16) выполнено, если 0 / > kT B(0 ) k A(0 ) -1, где 0 = Re p* / 2 ( p* - ближайший к мнимой оси полюс передаточной функции ( p) ) и 2 A() = max (i - ), B() = max (i - )(i - ) - . (17) 0 Пример 1. Рассмотрим уравнение автогенератора с гистерезисом, на который действует внешняя периодическая сила:

d[,0 ]t (t) + a (t) + (t) + = sin t, (18) dt где [,0] - гистерезисная функция, a - константа. Из теоремы 7 следует следующее t Предложение 7. Предположим, что a > 0 и гистерезисная функция [,0] t удовлетворяет условиям (12)-(14), где = 0. Если выполнено неравенство a + q 1 - q 2 > karcsin - arcsin B(0 ) k A(0 ) -1, 8 (i ) (i ) -q = 8l 1+ a2(a 4 - a2 ), 0 = a / 4, в случае 0 < a < 2, и неравенство 0 + r 1 - r > karcsin - arcsin B(0 ) k A(0 ) -1, 4 (i ) (i ) a - a2 - r = l(a + a2 - 4) (1+ a2 )(a2 - 4)-1, 0 =, в случае a > 2, то в любом режиме работы автогенератора, описываемого уравнением (18), будет наблюдаться явления захватывания под частоту внешнего периодического воздействия. ( Здесь A() и B() - числа, определяемые (17).) Пример 2. Рассмотрим систему уравнений, описывающую динамику автогенератора радиодиапазона с гистерезисом - 2h + + g = 0, (19) = - f [, f0 ]t - sin t, где f, g, h,, - положительные константы, - малый положительный параметр, гистерезисная функция f [, f0 ]t имеет вид f [, f0 ]t = m + [,0 ]t, причем m > 0, g > 2hm, гистерезисная функция [,0 ]t удовлетворяет условиям (12)-(14) с = 0. Введем обозначения: ( p) = p2 - 2hp +1, - (g - 2mh) / 4m, если g < 2m(h +1) 0 = - (g - 2mh) / 4m + (g - 2mh)2 /16m -1/ 4, если g > 2m(h +1) 2 (i - 0 ) m(i - 0 ) - g(i - 0 ) M = max, N = max.

0 g(i - 0 ) + m(i - 0 ) [ (i - 0 ) + m](i - 0 ) + g(i - 0 ) Предложение 8. Если выполнено неравенство 0 / > kT N Mk -1, (20) то для системы (19) при малых > 0 имеет место явление захватывания в классе всех гистерезисных функций [,0 ]t, удовлетворяющих условиям (12)-(14), где = 0.

Рассмотрим случай достаточно большой амплитуды внешнего гармонического воздействия. Тогда 1 - 1 + 2l T , (i ) и вместо условия (17) будем иметь следующее условие стабилизации системы (19) k( - 1 + 2l) N Mk - >.

0 (i ) Предложения 7 и 8 являются обобщениями на класс гистерезисных функций соответствующих утверждений, доказанных Г.А. Леоновым для однозначных непрерывно дифференцируемых функций ( ).

В пятой главе разрабатывается новый алгоритм стабилизации линейных систем управления. Здесь дано новое, элементарное, доказательство теоремы о стабилизации линейного объекта управления по полному выходу. В ходе доказательства теоремы дан конструктивный метод построения стабилизирующего регулятора u = Kx, дающего к тому же заданное расположение корней характеристического полинома (задача о размещении полюсов — «pole assignment problem»). Предложенный алгоритм построения стабилизирующей матрицы K является наиболее простым и «экономным» из существующих, он удобен и эффективен с вычислительной точки зрения и предполагает лишь повторение по существу одной единственной операции – приведения матрицы к диагональному виду (например, к жордановой нормальной форме). Вместе с элементарным доказательством теоремы о приведении матрицы к жордановой нормальной форме, данным А.Ф. Филлиповым, предложенное нами доказательство теоремы о стабилизации по состоянию становится элементарным в полном смысле слова.

Рассматривается задача о статическом регуляторе по полному выходу для линейной системы n m x = Ax + Bu, x R, u R (21) где A и B – вещественные постоянные матрицы размеров n n и n m соответственно.

Требуется найти вещественную (m n) – матрицу K такую, чтобы система (21), замкнутая обратной связью u = Kx, оказалась бы асимптотически устойчивой, т.е. чтобы матрица A + BK была бы гурвицевой.

Решение выше сформулированной задачи дает следующая Теорема 8 (об управлении спектром матрицы). Пусть A и B – любые вещественные n (n n) –, (n m) – матрицы такие, что пара (A, B) полностью управляема, и { } C – j j=произвольный набор комплексных чисел, замкнутых относительно операции комплексного сопряжения. Тогда существует вещественная (m n) – матрица K такая, что спектр n n (A + BK) матрицы A + BK совпадает с набором { } : (A + BK) ={ }.

j j j =1 j =Теорема 8 была впервые доказана для частного случая m = 1 Дж.Е. Бертрамом в 19году, используя метод корневого годографа. Этот же случай рассматривали также Р. Калман, Дж. Риссанен и Х. Розенброк. Для случая m > 1 частные результаты были получены В.М. Поповым, С.Е. Лангенхопом, Дж.О. Симоном и С.К. Миттером, П.А. Бруновским.

Впервые полное доказательство вышеприведенной теоремы об управлении спектром матрицы было дано В.И. Зубовым (1966). (В зарубежной англоязычной литературе первенство незаслуженно приписывается У.М. Уонэму (1967)). Позже в своей работе 19Уонэм сформулировал проблему управления спектром матрицы (pole assignment problem) и в этой же работе дал другое, отличное от доказательства В.И. Зубова, доказательство сформулированной теоремы. После выхода работ Зубова и Уонэма было написано немалое число работ, в которых предлагались другие, альтернативные, доказательства теоремы 8.

Однако следует отметить, что все известные доказательства теоремы 8 достаточно громоздки для вещественного векторного случая ( m > 1). Предложенное нами доказательство является полным и наиболее простым. Проблема управления спектром матрицы систем с неполной обратной связью также рассматривалась многими авторами, обзоры работ которых можно найти, например, в статье Б.Т. Поляка и П.С. Щербакова. Задача об управлении решениями линейного дифференциального уравнения рассматривалась в работах И.Н. Сергеева.

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

Мотивацией к исследованию стабилизации путем введения запаздывания явились компьютерные эксперименты К. Пирагоса (K. Pyragas) по стабилизации хаоса – стабилизации неустойчивых периодических орбит, погруженных в странный аттрактор той или иной хаотической системы. Здесь, в третьей главе, получены необходимые и достаточные условия стабилизируемости неустойчивых двумерных и трехмерных систем с постоянными коэффициентами путем введения обратной связи с запаздыванием. При этом рассматриваются два вида запаздывания в обратной связи: обычная и по Пирагосу.

Доказанные теоремы в целом хорошо иллюстрируют эффективность введения запаздывания в обратной связи для стабилизации линейных управляемых систем. Результаты получены на основе D - разбиения пространства параметров рассматриваемых систем. Они позволяют сделать вывод о возможностях линейной стационарной обратной связи с запаздыванием для стабилизации линейных неустойчивых стационарных систем. Оказывается, что как для обычной обратной связи с запаздыванием, так и для обратной связи по Пирагосу линейная система в седловом случае не является стабилизируемой ни при каком постоянном коэффициенте усиления и ни при каком запаздывании. Поэтому для обеспечения эффекта Пирагоса необходимо ввести зависящий от времени коэффициент усиления в обратной связи. Приведем основные результаты.

Рассмотрим линейную стационарную систему x(t) = A x(t) + bu(t), y(t) = c x(t) (22) n n где x(t) R - вектор состояния, u(t) R - управление (вход), y(t) R - выход, A - вещественная постоянная матрица размера n n, b и c - соответственно вектор-столбец и вектор-строка размерностей n.

Рассмотрим два способа введения обратной связи с запаздыванием: обычная u(t) = ky(t - ) (23) и по Пирагосу u(t) = k[y(t - ) - y(t)], (24) где k 0 и > 0 - варьируемые параметры.

Основная задача: требуется найти значения параметров k и такие, чтобы замкнутая система (22),(23)/(22),(24) оказалась асимптотически устойчивой.

Здесь основная задача решается для случая двумерных и трехмерных управляемых систем.

Предложение 9. Пусть в (22) n = 1 и A > 0. Тогда систему (22) невозможно стабилизировать обратной связью Пирагоса (24) ни при каких значениях параметров k и 0, а обратной связью (23) — можно: областью стабилизации является интервал 0 < A < 1 .

Двумерные системы. Пусть в (22) n = 2 и систему (22) можно привести к следующему виду x1(t) = x2(t), x2(t) = - a1 x1(t) - a2 x2(t) - u(t), y(t) = c1 x1(t) + c2 x2(t), (25) где a1, a2 ; c1, c2 - вещественные параметры.

Возможны три случая: 1) c1 0, c2 = 0, 2) c1 = 0, c2 0, 3) c1 0, c2 0.

Теорема 9. Пусть в системе (25) c1 0 (c1 := 1), c2 = 0. Для того чтобы система (25) была стабилизируема обратной связью (23) необходимо и достаточно, чтобы a1 0, a2 > 0 или a1 > 0, a2 > - 2a1.

Теорема 10. Пусть в системе (25) c1 0 (c1 := 1), c2 = 0. Для того чтобы система (25) была стабилизируема обратной связью (24) необходимо и достаточно, чтобы a1 > 0.

Теорема 11. Пусть в системе (25) c1 = 0, c2 0 (c2 := 1). Тогда необходимым и достаточным условием стабилизируемости системы (25) обратной связью (23) является выполнение неравенства a1 > 0.

Теорема 12. Пусть в системе (25) c1 = 0, c2 0 (c2 :=1). Тогда для стабилизируемости системы (25) обратной связью (24) необходимо и достаточно, чтобы 2 2 a1 > 0, a2 > 0 или a2 0, a1 > a2 /16, где = min (cos + sin / ) ( -1,0419).

[0,2 ] Для системы (25), где a1 =1, a2 = -d, c1 = 0, c2 =1, интервал стабилизации для параметра d : 0 d < 4 / | |, даваемый теоремой 12, хорошо согласуется с интервалом стабилизации 0 < d <1,216, полученным Пирагосом с помощью компьютерного эксперимента.

Теорема 13. Пусть в системе (25) c1 = 0, с2 0 (c2 := 1). Тогда для стабилизируемости системы (25) обратной связью (23), необходимо и достаточно, чтобы выполнялось хотя бы одно из условий:

a) c1 > 0, b) c1 < 0, c1a2 < a1 0, c) c1 < 0, a1 > 0, c1a2 < a1(a1 + 2c1 ).

Нижней границей области стабилизации в плоскости параметров (a1, a2 ), определяемой условием c) теоремы 13, является ветвь параболы 2 4 2 (a1 + c1 )2 / c1 - a2 / c1 = с асимптотой a1 = c1a2 - c1.

Теорема 14. Пусть в (25) c1 0, c2 0 (c2 :=1). Тогда для стабилизируемости системы (25) обратной связью (24), необходимо и достаточно выполнение хотя бы одного из условий:

a) c1 > 0, a1 > 0 или b) c1 < 0, a1 > 0, a2 > c1.

Рассмотрим вопрос о стабилизации системы (22), где n = m = 2, B = E ( E - единичная (2 2) - матрица), посредством полной обратной связи по Пирагосу u(t) = K[x(t - ) - x(t)]. (26) Здесь K - варьируемая (2 2) - матрица, а > 0 - скалярный параметр (запаздывание). Пусть det(E - A) = 2 + a2 + a1 (a1, a2 R ) - характеристический полином матрицы A.

Теорема 15. Для стабилизируемости системы (22) (n = m = 2, B = E) обратной связью (26) необходимо и достаточно, чтобы a1 > 0.

Трехмерные системы. Пусть в (22) n = 3 и систему (22) можно привести к следующему виду x1(t) = x2 (t), x2 (t) = x3 (t), x3 (t) = -a1x1(t) - a2 x2 (t) - a3 x3 (t) - u(t), (27) y(t) = c1x1(t) + c2 x2 (t) + c3 x3 (t), где ai, ci (i =1,2,3) - вещественные параметры. Справедливы следующие теоремы.

Теорема 16. Пусть в системе (27) c1 0 (c1 :=1), c2 = c3 = 0. Тогда для стабилизируемости системы (27) обратной связью (23) необходимо и достаточно, чтобы a2 > 0, a3 > 0 или a2 < 0, a3 > 0, a2 < 2a1a3.

Теорема 17. Пусть в системе (27) c1 0 (c1 :=1), c2 = c3 = 0. Тогда система (27) стабилизируема обратной связью (24) в том и только том случае, если a1 > 0, a3 > 0.

Теорема 18. Пусть в системе (27). Тогда для c2 0 (c2 :=1), c2 = c3 = стабилизируемости системы (27) обратной связью (3) необходимо и достаточно, чтобы a1 > 0, a3 > 0.

Теорема 19. Пусть в системе (27) c2 0 (c2 := 1), c1 = c3 = 0. Тогда для стабилизируемости системы (27) обратной связью (24) достаточно, чтобы выполнялось хотя бы одно из условий:

8 a) 0 < a1 < a2a3 + a2 a2, a2 > 0, a3 > 0, 2 - ( - 8) - 4 b) a1 > a3 (-a2 ) + (-a2 ) - a2, a2 < 0, a3 > 0.

4 - (4 - ) 4 - Предложение 10. Необходимым условием стабилизируемости системы (27) обратной связью (24) является неравенство a1 > 0.

Теорема 20. Пусть в системе (27) c3 0 (c3 :=1), c1 = c2 = 0. Тогда для стабилизируемости системы (27) обратной связью (23) достаточно, чтобы выполнялось хотя бы одно из следующих условий:

a) a1 > 0, a2 > 0, + 3 3 a2 , b) a2 < 0, a3 > 0, a1 > a2 a3 + - - 3 - 27 + 3 3 a2 .

c) a3 < 0, a2 < a3, a1 > a2 a3 + - (6 - 3)2 - 3 Теорема 21. Пусть в системе (27) c3 0 (c3 :=1), c1 = c2 = 0. Тогда для стабилизируемости системы (27) обратной связью (24) достаточно выполнения хотя бы одного из следующих условий:

4 a2 , a) a2 > 0, a3 > 0, 0 < a1 < a2 a3 + (- ) 12 a1a +12(2 - 3) ab) a1 > 0, a3 > 0, a2 > -, a3 [ +12(2 - 3)] + 3 3 a2 , 0 > a2 > - 27 a3.

c) a1 > a2a3 - - - 3 3 (6 + 3) Здесь то же самое число, что и в теореме 12.

Теорема 22. Пусть в системе (22) n = m = 3, B = E ( E - единичная (3 3) - матрица).

Пусть далее det(E - A) = 3 + a32 + a2 + a1 - характеристический полином матрицы A.

Тогда система (22) стабилизируема обратной связью (26), где K - (3 3) - матрица, если a1 > 0, a3 > 0.

Условия a1 > 0, a3 > 0 являются также и необходимыми для стабилизируемости системы (22) (n = m = 3, B = E) обратной связью вида (26).

Сформулированные выше теоремы показывают возможности стабилизации двумерных и трехмерных стационарных линейных систем обратной связью с запаздыванием двух видов: обычной (23) и по Пирагосу (24). Теоремы 9,10,13 и 14 (двумерные системы) и 16, 17, 20 и 21 (трехмерные системы) хорошо показывают как введение запаздывания в обратной связи расширяет возможности обычной стационарной стабилизации u(t) = ky(t) (k = const) без запаздывания. Исключения составляют теоремы 11 и 12 ( n = 2 ;

случай c1 = 0, c2 0 ) и 18, 19 ( n = 3 ; случай c1 = c3, c2 0 ). Из теорем 10 и 17 следует, что в рассматриваемых случаях области стабилизации обратной связью с запаздыванием вида (23) и без запаздывания совпадают. А из теорем 12 и 19 следует, что при стабилизации обратной связью по Пирагосу (24) не происходит расширения области стабилизации, получаемой обычной стационарной обратной связью без запаздывания. Далее предложение 1 (n =1) и теоремы 15 (n = 2) и 22 (n=3) показывают, что если матрица разомкнутой системы имеет одно положительное или нулевое собственное число, то стабилизация систем обратной связью по Пирагосу невозможно ни при каких значениях варьируемых параметров:

коэффициента усиления и запаздывания. Поэтому для обеспечения эффекта Пирагоса следует взять коэффициент усиления зависящим от времени.

В седьмой главе рассматривается задача о стабилизации динамической системы «машина – регулятор Уатта». Рассматривается модифицированный регулятор Уатта.

Доказаны теоремы об асимптотической устойчивости положения равновесия системы, дающие также оценку снизу области притяжения в фазовом пространстве. Полученные достаточные условия устойчивости близки к необходимым. Здесь также проведен нелокальный анализ переходного процесса. Далее, рассматривается случай, когда коэффициент трения возмущен гауссовским «белым» шумом. Доказаны соответствующие теоремы об устойчивости по вероятности. Нелинейности, входящие в уравнения, описывающие динамику работы системы, принадлежат классу функций, не рассматривавшиеся ранее в литературе. А именно, в работах многих авторов изучались частные случаи системы дифференциальных уравнений, к которой сводятся рассматриваемые нами уравнения движения на нелинейные функции. В этих работах на нелинейные функции накладывалось условие «секториального типа» - график той или иной нелинейной функции должен лежать в некотором секторе, содержащемся в первой и второй координатных четвертях. Помимо этого, указанные функции у упомянутых работах подчинялись обобщенным в какой-либо форме условиям Гаусса – Гурвица. Рассматриваемые же нами нелинейности таковы (квадратичного типа), что они не удовлетворяют указанным выше условиям. Для исследования системы применятся специальный прием, использующий модифицированный метод функций Ляпунова.

Работа системы «машина – регулятор Уатта» при некоторой идеализации описывается следующей нелинейной системой дифференциальных уравнений:

J = F(x), (28) x + x + x = m m(r + x), где , , , J, m,r - некоторые положительные числа, F(x) непрерывно дифференцируемая на прямой (-,+ ) функция, = (t), x = x(t) - неизвестные функции.

Рассматривается также обобщенная система (28), когда коэффициент трения является функцией от x, а восстанавливающая сила x заменена на нелинейную функцию от x. Отметим, что в частном случае, когда 1) F(x) - линейная функция, 2) правая часть (центробежная сила) второго уравнения имеет приближенный вид mr(t)2 + m0 2 x, система (28) рассматривалась в одной из работ Г.А. Леонова. Там же была поставлена задача о проведении исследования системы «машина – регулятор Уатта» без предположений 1) и 2).

Предположим, что существует число x0 > 0 такое, что F(x0 ) = 0, F = F (x0 ) 0 и F(x) 0 при x x0. Тогда (0, x0 ), где 0 2 = x0 m(r + x0 ), есть единственное положение равновесия системы (28).

Задача заключается в том, чтобы перевести при включении систему «машина – регулятор Уатта» из неподвижного состояния ((t) 0, x(t) 0, x(t) 0 ) в стационарное рабочее состояние ((t) 0, x(t) x0, x(t) 0 ) – стационарный режим.

Вводя обозначения 0 y = F (x - x0 ) J, z = y, (y) = F(x0 + Jy / F ) J, 0 a = / m, b = / m, c = (-F )x0 J m, k = (-F ) (r + x0 ) J, систему (28) можно привести к виду 2 = ( y), y = z, z = -az - by + y - k + c (29) Здесь a > 0, b > 0, (0) = 0. Введем функцию ляпуновского типа 2 2 2V (, y, z) = 2(k - 3c) / 3 + 2y(k - c) + y2 (b - ) + (z + ay)2 (30) и множество ( = 30 ) K = {(, y, z) :V (, y, z) 0, 0 }.

Теорема 23. Пусть в системе (29) функция (y) непрерывно дифференцируема для всех y, (0) = 0. Пусть далее выполнены следующие условия:

1) F < 0, 2x0 < r ;

2) существует положительная константа > 0 такая, что 0 y((y) y -1) y0 для всех y [ y1, y2 ], где y1 < 1, y1 = 4F x0 (r + x0 ) J (r - 2x0 ), y2 = 2(-F )x0 J ;

J[J (r - 2x0 )2 + 3F f0mr(r + x0 )] 3) R 0 2(-F )x0 (r - 2x0 )[J (r + x0 ) - 3F f0mr]> , где f0 = 2m0 (r + x0 ). (31) Тогда все траектории системы (29), начинающиеся в области K, остаются в ней, и положение равновесия (0,0,0) асимптотически устойчиво по Ляпунову с областью притяжения K.

В случае, когда (y) y ( = 0), условие 2) теоремы 23 отсутствует, а условие 3) принимает вид J ( - 2x0 ) > 3.

(-F ) f0mr(r + x0 ) Сравнив условия 1) и 3) теоремы 23 в случае (y) y с необходимыми и достаточными условиями Рауса – Гурвица 0 F < 0, J ( - m0 2 ) > (-F ) f0m асимптотической устойчивости линеаризованной в окрестности положения равновесия (0,0,0) системы (29), можно убедиться, что эти условия близки друг друга.

Непосредственным следствием теоремы 23 является следующая Теорема 24. (О переходном процессе.) Пусть в системе (28) функция F(x) непрерывно дифференцируема для всех x, причем F(x0 ) = 0, F = F (x0 ) < 0, F(x) при x x0 для некоторого x0 такого, что 2x0 < r. Пусть далее выполняются следующие условия:

1) существует положительная константа k > 0 такая, что 0 (1- F(x) / F (x - x0 ))(x - x0 ) k(x - x0 )для всех x [x1, x2 ], x x0, где x1 = -x0, x2 = x0 (5r + 2x0 ) (r - 2x0 ) и 2k(r + x0 ) < r 2x0 -1;

2) имеет место неравенство (-F )R / J > k, где R - число, определяемое из (31);

3) < m.

Тогда для любого решения ((t), x(t)) системы (28) с начальными условиями (0) = 0, x(0) = 0, x(0) = 0 справедливы следующие соотношения:

(t) [0, ], = 30, lim (t) = 0, lim x(t) = x0, lim x(t) = 0.

t+ t+ t+ Теоремы 23 и 24 допускают обобщение и на случай, когда коэффициент = (x) функция от x, а x заменена на нелинейную функцию (x). Формулировки соответствующих теорем аналогичны теоремам 23 и 24 соответственно. Рассматривается также случай, когда коэффициент трения в системе (28) возмущен «белым» шумом (t) с коэффициентом диффузии = (, x).

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

1) Разработан новый, простой алгоритм стабилизации по состоянию линейных стационарных управляемых систем.

2) Разработаны новые алгоритмы стабилизации двумерных и трехмерных линейных управляемых систем с помощью обратной связи с запаздыванием. Получены коэффициентные необходимые и достаточные условия стабилизируемости систем, расширяющие область стабилизации, даваемой обратной связью без запаздывания.

3) Выявлены возможности стабилизации линейных систем обратной связью с запаздыванием: обычной и по Пирагосу. Показано, что в седловом случае линейная система не является стабилизируемой. Сделан вывод о необходимости принятия коэффициента усиления в обратной связи с запаздыванием зависящим от времени, чтобы обеспечить эффект стабилизации хаоса в компьютерных экспериментах К. Пирагоса.

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

5) Получены новые эффективные частотные критерии глобальной асимптотики абсолютной устойчивости и дихотомичности дифференциальных систем с гистерезисными нелинейностями. Рассмотрены примеры: уравнение второго порядка и система уравнений второго порядка с гистерезисом.

6) Получено новое достаточное условие (близкое к необходимому) для стабилизации динамической системы «машина-регулятор Уатта». Доказаны теоремы об асимптотической устойчивости положения равновесия системы, получена оценка снизу области притяжения в фазовом пространстве.

7) Проведен нелокальный анализ переходного процесса в динамической системе «машина-регулятор Уатта».

Автор выражает глубокую признательность своему научному консультанту членукорреспонденту РАН, профессору Геннадию Алексеевичу Леонову за постоянный интерес к работе автора, поддержку и обсуждение результатов.

Список публикаций по теме диссертации Монографии 1. Леонов Г.А., Шумафов М.М. Проблемы стабилизации линейных управляемых систем.

СПб. Изд-во С.-Петербург. ун-та. 2002. 307 с.

2. Леонов Г.А., Шумафов М.М. Методы стабилизации линейных управляемых систем.

СПб. : Изд-во С.-Петербург. ун-та. 2005. 420 с.

3. Leonov G.A., Shumafov M.M. Stabilization of Linear Systems. Cambridge: Cambridge Scientific Publishers, 2012. 408 p.

4. Леонов Г.А., Шумафов М.М., Тешев В.А. Устойчивость систем с гистерезисом.

Майкоп: Изд-во Адыгейского государственного университета.2012. 182с.

Публикации в изданиях, рекомендованных ВАК РФ 1. Шумафов М.М. О переходном процессе в динамической системе с регулятором Уатта // Вестник С.- Петербург. ун-та. Сер.мат., мех. астр. 2001. Вып. 1 (№1). С. 53-59.

2. Шумафов М.М. Об асимптотической устойчивости динамической системы с регулятором Уатта // Дифф. ур. 2002. Т. 39. №1. С. 57-62.

3. Леонов Г.А., Шумафов М.М. Элементарное доказательство теоремы о стабилизируемости линейных управляемых систем // Вестник С.- Петербург. ун-та.

Сер. матем. мех., астрономия. 2003. Вып. 3 (№17). С. 56-68.

4. Леонов Г.А., Шумафов М.М. Алгоритм пошаговой стабилизации линейного объекта управления // Изв. вузов. Северо-Кавказ. регион. Естеств.. науки. 2005. №2. С. 14-19.

5. Шумафов М.М. К задаче стабилизации двумерной линейной дискретной системы // Изв. вузов. Северо-Кавказский регион. Естественные науки. 2009. №5. С. 71-74.

6. Шумафов М.М. О стабилизации двумерных линейных управляемых систем обратной связью с запаздыванием // Вестник Адыгейского государственного университета.

2010. Вып. 2. С. 40-52.

7. Шумафов М.М. Стабилизация линейных стационарных управляемых систем второго порядка обратной связью с запаздыванием // Изв. вузов. Математика. 2010. №12. С.

87-90.

8. Leonov G.A., Shumafov M.M. Stabilization of Controllable Linear Systems // Nonlinear Dynamics and System Theory. 2010. V.10. №3. P. 235-268.

9. Leonov G.A., Kuznetsov N.V., Seledzhi S.M., Shumafov M.M. Stabilization of unstable control system via design of delayed feedback. International Conference on Applied and Computational Mathematics, 2011, pp. 18-25 [ISSN 1792-4235].

10. Leonov G.A., Shumafov M.M. Vibrational Stabilization and the Brockett problem // Differential Equations. 2011. V.47. №13. P.1853-1915.

11. Леонов Г.А., Шумафов М.М. Стабилизация линейных управляемых систем с запаздывающим аргументом в цепи с обратной связью // Докл. РАН. 2012 (в печати).

12. Шумафов М.М. Стабилизация систем с гистерезисными нелинейностями гармоническим внешним воздействием // Вестник Адыгейского государственного университета. 2012. Вып.3. С.11-19.

13. Шумафов М.М. Устойчивость систем дифференциальных уравнений с гистерезисными нелинейностями// Вестник Адыгейского государственного университета. 2012. Вып.3. С.20-31.

Статьи, материалы конференций 1. Шумафов М.М. О диссипативности и устойчивости по вероятности случайных процессов, определяемых некоторыми нелинейными дифференциальными уравнениями второго порядка / Третья региональная теоретическая конференция молодых ученых Северного Кавказа. (Тезисы докладов и сообщений.) Майкоп. 1990.

С.261-263.

2. Шумафов М.М. О диссипативности случайных процессов, определяемых некоторыми нелинейными дифференциальными уравнениями второго порядка // Дифф. уравнения.

– 1993. –Т.29. -№1. – С.175-176.

3. Тешев В.А., Шумафов М.М. Частотный критерий устойчивости систем дифференциальных уравнений с гистерезисными нелинейностями в критическом случае // Труды Физ. общ. Респ. Адыгея. - 1996. - №1. -С. 48-50.

4. Пономаренко Д.В., Тешев В.А., Шепелявый А.И., Шумафов М.М. Частотный критерий синхронизации гистерезисных систем / Первая научно-практическая конференция «Дифференциальные уравнения и их применения». СПб. 1996. С. 177.

5. Шумафов М.М. Частотный критерий устойчивости систем с гистерезисными нелинейностями в критическом случае / Международная конференция «Нелокальные краевые задачи и родственные проблемы математической биологии, информатики и физики». (Тезисы докладов.) Нальчик. 1996. С.106.

6. Шумафов М.М. О построении функций Ляпунова для двумерных стохастических систем / Тезисы докладов второй международной конференции «Дифференциальные уравнения и их применения». Саранск. 1996. С. 127.

7. Ponomarenko D.V. Teshev V.A. Shepeljavyi A.I. Shumafov M.M. Stabilization of Systems with Hysteresis by Periodic External Force / In book: International Conference on Control of Oscillations and Chaos. S. -Petersburg. 1997.

8. Тешев В.А., Шепелявый А.И., Шумафов М.М. Частотный критерий стабилизации гармоническим внешним воздействием систем с гистерезисными нелинейностями / Нелинейные динамические системы. Выпуск 1. С.Петербург: СПбГУ. 1997. С. 261-280.

9. Шумафов М.М. Функции Ляпунова для двумерных линейных стационарных стохастических систем // Труды Физ. Общ. Респ. Адыгея. 1997, №2. С.1-26.

10. Шумафов М.М. Об одной системе дифференциальных уравнений, возникающей в теории фильтрации / Тезисы докладов четвертой Северо-Кавказской региональной конференции «Функционально-дифференциальные уравнения и их приложения». - Махачкала. 1997. С. 102-103.

11. Тешев В.А., Шумафов М.М. О глобальной асимптотической устойчивости решений уравнений второго порядка с гистерезисными нелинейностями // Труды Физ. Общ.

Респ. Адыгея. 1997, №2. С.61-69.

12. Леонов Г.А., Тешев В.А., Шумафов М.М. О глобальной асимптотической устойчивости систем с гистерезисными нелинейностями / Труды третьей международной конференции «Дифференциальные уравнения и их приложения».- Саранск. - 1998. -С. 36-37.

13. Тешев В.А., Шепелявый А.И., Шумафов М.М. Явление захватывания в системах с гистерезисным элементом // Труды Физ. Общ. Респ. Адыгея. 1998, №3. С.49-57.

14. Тешев В.А., Шумафов М.М. Частотный критерий дихотомичности регулируемых нелинейных систем с гистерезисным элементом // Труды Физ. Общ. Респ. Адыгея.

1999, №4 С.34-39.

15. Shumafov M.M. On the stochastic stability of a nonlinear system perturbed by a “white” noise random process // Trudy Fiz. Obsch. Resp. Adygeya (FORA). 1999. №4. P. 118-124.

16. Шумафов М.М. Об асимптотическом поведении решений одной нелинейной трехмерной системы //Труды Физ. Общ. Респ. Адыгея. 2001, №6 С.48-50.

17. Тешев В.А., Шумафов М.М. Об абсолютной устойчивости и дихотомичности систем дифференциальных уравнений с гистерезисными функциями // Труды Физ. Общ.

Респ. Адыгея. 2001. №6. С. 51.

18. Шумафов М.М. О существовании гомоклинической орбиты в обобщенной системе Лоренца // Труды Физ. Общ. Респ. Адыгея. 2001. №6. С. 57-66.

19. Shumafov M.M. On the stability of a second-order nonlinear stochastic system // Trudy Fiz.

Obsch. Resp. Adygeya. 2002. №7. Р. 98-102.

20. Шумафов М.М. Об асимптотической устойчивости двумерной нелинейной динамической системы, возмущенной “белым” шумом // Труды Физ. Общ. Респ.

Адыгея. 2004. №9. С. 106-109.

21. Leonov G.A., Shumafov M.M., Stabilization of Linear Systems / Proc. of Fourth European Conference on Structural Control (4ECSC). S. –Petersburg. 2008. №2. P. 461-494.

22. Шумафов М.М. О диссипативности решений стохастических дифференциальных уравнений второго порядка // Вестник Адыг. гос. университета. Серия «Естественноматематические и технические науки». 2008. Вып.4 (32).С.11-17.

23. Сташ А.Х., Шумафов М.М. Устойчивость и ограниченность решений линейных дифференциальных систем // Труды Физ. Общ. Респ. Адыгея. 2008. №13. С. 12-21.

24. Шумафов М.М., Цей Р.М. Математическое моделирование и обратные задачи // Вестник Адыгейского государственного университета. Серия «Естественноматематические и технические науки». 2008. Вып.4 (32). С.18-24.

25. Шумафов М.М., Цей Р.М. Метод модулирующих функций и его применение при решении обратных задач // Вестник Адыг. гос. университета. Серия «Естественноматематические и технические науки». 2008. Вып.9 (37). С.9-22.

26. Шумафов М.М., Цей Р.М. Алгоритм решения задачи определения фильтрационноемкостных параметров газоносного пласта методом модулирующих функций // Вестник Адыгейского государственного университета. Серия «Естественноматематические и технические науки». 2008. Вып.9 (37). С. 23-26.

27. Шумафов М.М., Цей Р.М. Решение обратных коэффициентных задач методом модулирующих функций / /Наука- 2008. Ежегодный сборник научных статей ученых и аспирантов Адыгейского госуниверситета. Майкоп. Изд-во АГУ. 2008. С.8-16.

28. Шумафов М.М., Цей Р.М. Идентификация параметров газоносного пласта на основе решения обратной задачи теории фильтрации // Вестник Адыгейского госуниверситета. Серия «Естественно-математические и технические науки». 2009.

Вып.1 (43). С. 33-42.

29. Шумафов М.М., Цей Р.М. Разработка алгоритма для численного решения обратной задачи теории фильтрации методом модулирующих функций // Вестник Адыгейского государственного университета. Серия «Естественно-математические и технические науки». 2009. Вып.1 (43). С. 43-49.

30. Шумафов М.М., Цей Р.М. К вопросу об определении фильтрационно-емкостных параметров газоносного пласта // СамДиф: Конференция «Дифференциальные уравнения и их приложения». Самара, 29 июля – 2 июня 2009г, Тезисы докладов.

Самара: Изд-во «Универс групп». 2009. С. 69-70.

31. Шумафов М.М. О стохастической устойчивости некоторых двумерных динамических систем // Дифф.уравн. 2010. Т.46. №6. С.1-5.

32. Шумафов М.М. Стохастическая устойчивость двумерных линейных стационарных систем // Вестник Адыгейского государственного университета. 2010. Вып.1. С.9-21.






© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.