На правах рукописи
Перкин Алексей Александрович
Устойчивость математических моделей систем фазовой синхронизации
05.13.18 - математическое моделирование, численные методы и комплексы программ
Автореферат диссертации на соискание учёной степени кандидата физико-математических наук
Великий Новгород - 2012
Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Санкт-Петербургский государственный архитектурно-строительный университет на кафедре прикладной математики и информатики.
Научный консультант: доктор физико-математических наук, доцент Смирнова Вера Борисовна
Официальные оппоненты: доктор физико-математических наук, доцент Чурилов Александр Николаевич (Санкт-Петербургский государственный морской технический университет, профессор) кандидат физико-математических наук, доцент Кондратьева Наталья Васильевна (Санкт-Петербургский государственный политехнический университет, доцент)
Ведущая организация: Федеральное государственное бюджетное учреждение науки Институт проблем машиноведения Российской академии наук (ИПМаш РАН)
Защита состоится 28 мая 2012 г. в 16:00 на заседании диссертационного совета Д.212.168.04 при Новгородском государственном университете имени Ярослава Мудрого по адресу 173003, г. Великий Новгород, ул. Большая Санкт-Петербургская, д.41.
С диссертацией можно ознакомиться в библиотеке Новгородского государственного университета имени Ярослава Мудрого.
Автореферат разослан л____ апреля 2012 г.
Ученый секретарь диссертационного совета Д.212.168.к.ф.-м.н., доцент Токмачев М.С.
Общая характеристика работы
Актуальность темы. Системы, основанные на принципе фазовой синхронизации, широко применяются в различных областях механики, радиотехники, электроники, связи. К ним можно отнести системы фазовой автоподстройки частоты (ФАПЧ), системы электропривода, вибрационные машины, системы дальней связи.
Задачи проектирования и анализа работы систем фазовой синхронизации требуют подробного изучения влияния значений конструктивных параметров системы на установление в ней того или иного режима.
Особенно важно определить область значений параметров, при которых для любых начальных состояний системы в ней устанавливается режим синхронизма. Такие системы называются глобально асимптотически устойчивыми. К задаче о глобальной устойчивости тесно примыкает задача об определении частот периодических режимов и задача об оценке наибольшего отклонения выхода системы от его начального состояния.
Задачам асимптотического поведения систем фазовой синхронизации различной природы посвящено много работ. Фазовые системы второго и третьего порядков, в том числе и с запаздыванием в петле обратной связи, рассматривались в исследованиях В.В. Шахгильдяна, А.А. Ляховкина, Ю.Н. Бакаева, Л.Н. Белюстиной, В.Н. Белыха, Л.З. Фишмана, В.Н. Сафонова, С.И. Евтянова, Б. Бисваса, П. Банерджи, А. Бхаташариа, И.И. Блехмана, Л. Шперлинга, А.С. Сомолиноса, О.Р. Рыкина, Е.Л. Урмана, Е.И. Юревича и многих других авторов.
К решению задач асимптотического поведения фазовых систем применялись различные качественные, качественно-численные, численные методы.
Поскольку системы фазовой синхронизации являются нелинейными системами автоматического управления, большое значение для их исследования получил второй (прямой) метод А.М. Ляпунова, широко применяемый и развиваемый с 30-х годов XX века. При этом специфика систем фазовой синхронизации: неединственность состояний равновесия и периодический характер нелинейности, потребовала построения новых классов функций Ляпунова.
Началом исследований в этом направлении явилось введение Ю.Н. Бакаевым и А.А. Гужем периодической функции Ляпунова для фазовых систем третьего порядка. В 70-е годы XX века Г.А. Леоновым в рамках второго метода Ляпунова была предложена серия методов, позволяющих конструировать функции Ляпунова специального вида для многомерных систем фазовой синхронизации: метод периодических функций Ляпунова, метод нелокального сведения, метод инвариантных конусов. Г.А. Леоновым, Ю.А. Корякиным, Ф. Райтманном, В.Б. Смирновой, О.Б. Ершовой, Л.С. Сперанской эти методы были распространены на бесконечномерные и дискретные системы. Результаты формулировались в терминах передаточной функции линейной части системы в виде эффективно проверяемых частотных неравенств.
Второй метод А.М. Ляпунова является методом математически обоснованным и, следовательно, не может приводить к ошибочным результатам, как приближенные методы. С другой стороны, оценки областей устойчивости, получаемые с его помощью в пространстве параметров системы, достаточно близки к истинным областям устойчивости.
Разработанные Г.А. Леоновым периодические функции Ляпунова успешно использовались затем другими исследователями (Н.В. Утина, А.И. Шепелявый, Я. Янг, Л. Хуанг, П. Лу) для изучения проблем устойчивости и проскальзывания циклов систем фазовой автоподстройки.
Результаты, доставляемые вторым методом А.М. Ляпунова, носят характер достаточных условий устойчивости. Они могут быть улучшены путем конструирования новых функций Ляпунова.
Задача построения новых функций Ляпунова и, как следствие, улучшения оценок областей глобальной асимптотической устойчивости, является весьма актуальной задачей теории систем фазовой синхронизации.
Цель работы состоит в получении новых эффективно проверяемых критериев глобальной асимптотической устойчивости и отсутствия периодических режимов второго рода определенной частоты для математических моделей сосредоточенных и распределенных систем фазовой синхронизации, а также реализации этих критериев для построения оценок областей устойчивости в пространстве параметров конкретных систем и оценок частоты колебаний второго рода для этих систем.
Для достижения поставленной цели необходимо было решить следующие задачи:
1. Построить модификации периодических функций Ляпунова, применяемых при исследовании устойчивости сосредоточенных систем управления с периодическими нелинейностями.
2. Разработать новый класс функционалов Попова для исследования асимптотического поведения распределенных систем с периодическими нелинейностями.
3. Получить многопараметрические частотные критерии глобальной асимптотической устойчивости фазовых системы, описываемых дифференциальными и интегро-дифференциальными уравнениями.
4. Преобразовать критерии глобальной устойчивости в критерии отсутствия у фазовой системы периодических режимов второго рода определенной частоты.
5. На основе частотных условий устойчивости установить оценки числа проскальзываний циклов для сосредоточенных и распределенных фазовых систем.
6. Разработать численные методы и соответствующие им алгоритмы для проверки частотных критериев устойчивости для систем фазовой автоподстройки частоты (ФАПЧ) с запаздыванием в петле обратной связи.
7. Разработать комплекс программ для реализации алгоритмов проверки частотных критериев для систем ФАПЧ с пропорционально-интегрирующим фильтром.
Методы исследования. В работе использованы второй метод Ляпунова, метод априорных интегральных оценок Попова, теорема Якубовича-Калмана о разрешимости специальных матричных неравенств, методика построения периодических функций Ляпунова, теория тригонометрических рядов.
Научная новизна диссертационного исследования заключается в следующем:
1. Для систем фазовой синхронизации, описываемых системами обыкновенных дифференциальных уравнений и интегро-дифференциальных уравнений Вольтерра с дифференцируемыми нелинейными функциями, получены Х новые многопараметрические частотные критерии глобальной асимптотической устойчивости;
Х новые частотные оценки числа проскальзываний циклов;
Х частотные условия, гарантирующие отсутствие циклов второго рода определенной частоты.
2. Разработаны численные методы и соответствующие им алгоритмы, обеспечивающие проверку частотных критериев для систем фазовой автоподстройки частоты с запаздыванием в петле обратной связи.
3. Разработан комплекс программ для оценки полос захвата и частоты биений для систем фазовой автоподстройки частоты с пропорционально-интегрирующим фильтром.
Практическая значимость диссертационного исследования заключается в следующем:
Х для различных математических моделей систем фазовой синхронизации установлены как достаточные условия глобальной устойчивости, так и условия отсутствия периодических решений второго рода, и получены новые оценки числа проскальзываний циклов;
Х полученные критерии устойчивости, а также реализующие их алгоритмы и программы могут быть эффективно использованы для построения оценок областей глобальной устойчивости конкретных систем фазовой синхронизации в пространстве их параметров;
Х полученные оценки числа проскальзываний циклов, а также реализующие их алгоритмы и программы позволяют выделять внутри областей устойчивости подмножества, характеризующиеся заданным числом проскальзываний циклов.
Комплекс программ, разработанный в диссертации, зарегистрирован в отделе регистрации программ для ЭВМ ФГУ ФИПС.
Достоверность изложенных в работе результатов. Области устойчивости в пространстве параметров конкретных систем синхронизации, полученные с помощью результатов диссертации, близки к областям устойчивости, полученным путем решения уравнений численными методами.
Основные положения, выносимые на защиту:
1. Многопараметрические частотные критерии глобальной асимптотической устойчивости для математических моделей сосредоточенных и распределенных систем фазовой синхронизации.
2. Частотные критерии отсутствия периодических решений второго рода для математических моделей систем фазовой синхронизации с сосредоточенными и распределенными параметрами.
3. Частотные оценки числа проскальзываний циклов для многомерных и бесконечномерных фазовых систем.
4. Численные методы и алгоритмы проверки частотных критериев для систем с дробно-рациональной передаточной функцией и запаздыванием в петле обратной связи.
5. Комплекс программ для построения оценок полос захвата и частоты биений для систем фазовой автоподстройки частоты с пропорционально-интегрирующим фильтром при наличии и отсутствии запаздывания в петле обратной связи.
Апробация работы. Основные результаты диссертационной работы докладывались на 57-й, 58-й, 59-й, 60-й, 61-й международных научно-технических конференциях молодых ученых СПбГАСУ (2004, 2005, 2006, 2007, 2008 гг.); 66-й, 67-й, 68-й научных конференциях профессоров, преподавателей, научных работников, инженеров и аспирантов СПбГАСУ (2009, 2010, 2011 гг.); международном конгрессе "6th International Congress on Industrial and Applied Mathematics" (Zurich, Switzerland, 16-July, 2007); X международном семинаре им. Е.С. Пятницкого, "Устойчивость и колебания нелинейных систем управления" (Москва, 3-6 июня 2008г); международной конференции "Sixth EUROMECH Nonlinear Dynamics Conference" (Saint-Petersburg, Russia, June 30 - July 4, 2008); VII всероссийской научной конференции "Нелинейные колебания механических систем" (Нижний Новгород, 22-26 сентября, 2008г.);
международной конференции "3-rd IEEE Multi-Conference on System and Control" (SaintPetersburg, Russia, July 8-10, 2009); международных конференциях "4-rd International Scientific Conference on Physics and Control" (Catania, Italy, 1-4 September, 2009), ''5-th International Scientific Conference on Physics and Control" (Leon, Spain, 5-8 September 2011); XI международной конференции "Устойчивость и колебания нелинейных систем управления (конференция Пятницкого)" (Москва, 1-4 июня, 2010г.); XIII конференции молодых ученых "Навигация и управление движением" (С.-Петербург, ЦНИИ "Электроприбор", 15-17 марта 2011г.); 14-й международной научно-технической конференции ''Моделирование, иденитификация, синтез систем управления'' (Канака, Украина, 11-18 сентября 2011г.) Публикации. Все результаты диссертации опубликованы в работах [1-24]; работы [1-4] опубликованы в изданиях, рекомендованных перечнем ВАК. В научных работах соискателя, написанных в соавторстве, В.Б. Смирновой принадлежит общее руководство работой, соавторам А.И. Шепелявому и Н.В. Утиной принадлежит исследование устойчивости дискретных систем фазовой синхронизации.
Объем и структура работы. Диссертация объемом 134 страницы состоит из четырех глав, заключения и списка литературы. Библиография содержит наименований.
Содержание работы Первая глава носит вводный характер. В ней, прежде всего, обосновываются исследуемые в диссертации математические модели систем фазовой синхронизации (СФС). В ней описываются основные задачи, возникающие при исследовании асимптотического поведения СФС, и методы, применяемые для их решения. Особое внимание уделяется прямому методу А.М. Ляпунова и методу априорных интегральных оценок В.-М. Попова. В этой главе сформулированы также основные результаты диссертации.
Вторая глава посвящена исследованию асимптотики сосредоточенных многомерных фазовых систем. Математической моделью таких систем является система обыкновенных дифференциальных уравнений.
& z = Az + Bf ( ), (1) & = C*z + Rf ( ).
Здесь A, B, C, R - действительные матрицы размерностей m m, ml, ml и l l соответственно, z(t) -- m -векторная функция ( z(t) = zk (t) ), (t) -- l -векторная k =1,Km функция ( (t) = (t) ), f ( ) - l -векторная функция с компонентами ( ) j j j j=1,Kl ( j = 1,Kl). Символом * здесь и далее обозначено эрмитово сопряжение.
Предполагается, что пара (A, B) наблюдаема, пара (A,C) управляема, а матрица A гурвицева. Предполагается, что ( ) ( j = 1,K,l) - -периодическая непрерывно j j j дифференцируемая функция, причем j ( )d < 0 ( j = 1,K,l).
j Предполагается также, что функция ( ) имеет на промежутке [0, ) два нуля j j 1 j < , причем 2 j ( ) +( ( )) 0 ( j = 1,K,l;k = 1,2).
kj kj Система (1) является системой непрямого управления. Ее основной характеристикой является передаточная матрица линейной части системы от входа f к & выходу (- ) K( p) = -R + C*(A - pEm)-1B ( p C), где через Em обозначена единичная m m -матрица. Все результаты этой главы сформулированы в терминах K( p).
Система (1) обладает счетным множеством положений равновесия как устойчивых в малом, так и неустойчивых.
Определение 1. Говорят, что система (1) глобально асимптотически устойчива, если любое ее решение стремится к какому-либо ее положению равновесия.
Определение 2. Решение {z = z(t), = (t)} системы (1) называется предельным циклом второго рода, если можно указать такое положительное число T и такие целые числа I 0 ( j = 1,K,l), что для всех t 0 выполнены равенства j z(t + T ) = z(t), (t +T ) = (t) + I ( j = 1,K,l).
j j j j 2 Величина = называется частотой предельного цикла второго рода.
T Во второй главе устанавливаются достаточные условия глобальной асимптотической устойчивости системы (1), а также достаточные условия отсутствия у системы (1) предельных циклов второго рода определенной частоты.
Введем в рассмотрение такие числа kj ( j = 1,K,l;k = 1,2), что d ( ) j j 1 j 2 j, R, d j и сформируем диагональные матрицы Ak = diag{k1,K,kl} (k = 1,2).
Введем далее числа j ( )d j = ( j = 1,K,l), j j | ( ) | d j ()d 0 j j 0 j = () | (1 - -1 'j ())(1- 11'j ())d 0 j | j 2 j j Параграфы 2.2 и 2.3 второй главы посвящены вопросу глобальной асимптотической устойчивости. В параграфе 2.2 доказано следующее утверждение.
Теорема 1. Пусть существует такая диагональная матрица = diag{1,K,l}, положительно определенные диагональные матрицы = diag{1,K,l}, = diag{1,K,l}, = diag{1,K,l} и такие числа ak [0,1] (k = 1,K,l), что выполняются следующие требования:
1) для всех 0 справедливо неравенство * * Re{K(i) - K (i)K(i) -(K(i) + iA1-1) (2) -(K(i) + iA2 )}- 0 (i2 = -1);
2) квадратичные формы 2 2 Gk (,, ) k + k + + kakk + ka0k , k 0k где a0k = 1- ak (k = 1,K,l), являются положительно определенными.
Тогда справедливы предельные соотношения z(t) 0 при t +, (t) c при t +, j j где (cj ) = 0 ( j = 1,K,l).
j Доказательство теоремы 1 основано на применении прямого метода Ляпунова и частотной теоремы Якубовича-Калмана. При построении периодической функции Ляпунова реализовано обобщение процедуры Бакаева-Гужа. Как и в случае стандартной процедуры Бакаева-Гужа, функция Ляпунова конструируется в виде квадратичная форма плюс интеграл от нелинейности, причем компоненты векторной нелинейной функции под знаком интеграла имеют те же периоды и те же нули, что функции , но j обладают нулевым средним на периоде. Однако используемая в диссертации нелинейная функция имеет более сложную структуру, чем в стандартной процедуре.
Теорема 1 является обобщением известных теорем о глобальной асимптотической устойчивости фазовых систем, доказанных с использованием периодических функций Ляпунова. Она позволяет использовать при проверке частотного неравенства большее число варьируемых параметров.
В параграфе 2.3 второй главы описано применение теоремы 1 для построения оценок полос захвата систем фазовой автоподстройки частоты (ФАПЧ) с пропорционально-интегрирующим фильтром (ПИФ) и синусоидальной характеристикой фазового детектора.
Система (1) при m = l =1 является распространенной математической моделью ФАПЧ с ПИФ. Здесь (t) описывает разность фаз эталонного и подстраиваемого генераторов.
Нелинейная функция в данном случае является скалярной. Она задается формулой f ( ) = sin - ( (0,1)), где - нормированная начальная расстройка.
Передаточная функция рассматриваемой системы имеет вид Tp + K( p) = T (0 < 1).
Tp +Области устойчивости ФАПЧ с ПИФ в пространстве параметров, доставляемые теоремой 1, сравниваются как с областями устойчивости, полученными численным решением системы (1), так и с областями, доставляемыми периодическими функциями Ляпунова, использованными Г.А. Леоновым [Гелиг А.Х., Леонов Г.А., Якубович В.А.
Устойчивость нелинейных систем с неединственным положением равновесия. М.:
Наука, 1978]. (Последние являются частным случаем функций Ляпунова, построенных при доказательстве теоремы 1, когда ak = 1 (k = 1,K,l) ). Для значения = 0.2 результаты проведенного сравнения представлены на рис. 1. Область устойчивости, полученная с помощью теоремы 1, лежит ниже кривой (3). Ниже кривой (1) лежит истинная область глобальной устойчивости. Ниже кривой (2) лежит область, полученная Леоновым Г.А.
при a1 = 1. Использование обобщенной процедуры Бакаева-Гужа позволило расширить область глобальной асимптотической устойчивости, доставляемую прямым методом Ляпунова.
Рис.1. Рис.2.
В параграфе 2.4 второй главы обобщение процедуры Бакаева-Гужа использовано с целью определения гарантированной верхней оценки для возможной частоты циклов второго рода.
Теорема 2. Пусть для какого-либо набора диагональных матриц , > 0, > 0, > и чисел ak [0,1], удовлетворяющих требованию 2) теоремы 1, частотное неравенство (2) выполнено для = 0 и для всех 0. Тогда система (1) не имеет предельных циклов второго рода частоты 0.
Доказательство теоремы 2 основано на разложении T -периодических фунций & z(t), (t) и f ( (t)) в ряды Фурье.
В параграфе 2.5 второй главы теорема 2 применена для оценки областей отсутствия режима биений в системе ФАПЧ с ПИФ и синусоидальной характеристикой фазового детектора, описанной в параграфе 2.3. Полученные области сравнивались с областями, доставляемыми процедурой Бакаева-Гужа с фиксированным параметром a1 = 1 [Сперанская Л.С. Автореферат канд. дисс. ''Частотные условия существования периодических решений нелинейных дифференциальных уравнений.'' Нижний Новгород, 1994].
Для случая = 0.2, T = 10 результаты приведены на рис.2. Кривая (3) соответствует теореме 2, кривая (2) - результатам Л.С. Сперанской. Области, где отсутствуют циклы второго рода, лежат левее соответствующих кривых. Теорема выделяет для каждого значения начальной расстройки более обширную гарантированную область отсутствия циклов II рода. На рис. 2. показана также кривая зависимости частоты биений от значений начальной расстройки, полученная численным методом - кривая (1).
В третьей главе результаты второй главы распространяются на фазовые системы с распределенными параметрами, математической моделью которых является система интегродифференциальных уравнений Вольтерра t & (t) = 0(t) + Rf ( (t - h)) - (t - ) f ( ( ))d (t 0;h 0). (3) Здесь (t) = (t), (t) = (t), f ( ) = ( ) - вектор-функции, j 0 0 j j j j=1,K,l j=1,K,l j=1,K,l R - l l -матрица, (t) - l l -матричная функция.
kj k, j=1,K,l Предполагается, что функции 0 j (t) непрерывны и стремятся к 0 при t +, а функции (t) измеримы. Предполагается также, что существует такое c > 0, что kj 0 ject,ij (t)ect L2[0,+) ( j,i = 1,K,l).
Относительно функции f ( ) остаются в силе все предположения главы 2.
Передаточная матрица линейной части рассматриваемой системы от входа f к выходу & (- ) имеет вид K( p) = -Re- ph + (t)e- ptdt ( p C). (4) Система (3) также, как и система (1), обладает счетным множеством положений равновесия.
В параграфе 3.2 главы 3 устанавливаются два частотных критерия глобальной асимптотической устойчивости системы (3).
Теорема 3. Пусть для некоторых положительных диагональных матриц = diag{1,K,l}, = diag{1,K,l}, = diag{1,K,l}, = diag{1,K,l} и чисел ak [0,1] (k = 1,K,l) выполнены требования 1) и 2) теоремы 1, где K( p) вычисляется по формуле (4). Тогда для решения (t) системы (3) справедливы предельные соотношения &k (t) 0 при t +, (t) ck при t +, k где k (ck ) = 0 (k = 1,K,l).
Теорема 4. Пусть для некоторых положительных диагональных матриц = diag{1,K,l}, = diag{1,K,l}, = diag{1,K,l}, = diag{1,K,l} и матричной функции (4) выполнено требование 1) теоремы 1. Пусть далее компоненты матриц , , удовлетворяют неравенствам 2 > ( j = 1,K,l), j j 2 j j где j ( )d j = 2 j j -1 -1+ j (1-1 j' ( ))(1-2 j' ( )) d.
| ( ) | j j j 0 j Тогда справедливо заключение теоремы 3.
В основе доказательства теорем 3 и 4 лежит метод априорных интегральных оценок Попова. Второй критерий доказан с помощью модификации процедуры БакаеваГужа. К функционалам Попова применяется обобщенная процедура Бакаева-Гужа, разработанная во второй главе.
В параграфе 3.3 приведены результаты применения теоремы 3 к исследованию системы ФАПЧ с ПИФ, синусоидальной характеристикой фазового детектора и запаздыванием в петле обратной связи. Ее передаточная функция имеет вид Tp + K( p) = T e- phT.
Tp +Нелинейная функция имеет вид f ( ) = sin - ( (0,1)).
Полученная граница области глобальной асимптотической устойчивости для h = 0,1 приведена на рис. 3. Это - кривая (2). Кривая (1) - граница, полученная в работе [Л.Н. Белюстина, М.С. Киняпина, Л.З. Фишман. Динамика систем фазовой синхронизации с запаздыванием. В кн. Теоретическая электротехника, № 48, Львов, 1990] качественно-численными методами. Области устойчивости лежат ниже своих границ.
Рис.3.
В параграфе 3.4 теорема 3 применяется для исследования задачи о самосинхронизации неуравновешенных роторов, находящихся на общей колеблющейся платформе. Динамика медленно изменяющихся фаз неуравновешенных роторов описывается системой двух дифференциальных уравнений второго порядка с синусоидальной нелинейной функцией, зависящей от их разности. Система дифференциальных уравнений сводится к интегральному уравнению вида (3), где (t) - разность медленно меняющихся составляющих фаз роторов. С помощью теоремы сформированы алгебраические требования на параметры механической системы, обеспечивающие синхронную работу двух роторов. Установлено, что условия самосинхронизации, доставляемые теоремой 3, являются менее жесткими, чем условия самосинхронизации, полученные с помощью стандартной процедуры Бакаева-Гужа.
Параграф 3.5 третьей главы посвящен изучению периодических решений второго рода системы (3).
Определение 3. Решение (t) системы (3) называется периодическим решением второго рода, если существует такое число T > 0 и такие целые числа Ik 0 (k = 1,K,l), что (t +T ) = (t) + Ikk (k = 1,K,l).
k k 2 Величина называется частотой периодического решения второго рода.
T Теорема 5 Пусть для каких-либо положительных диагональных матриц = diag{1,K,l}, = diag{1,K,l}, = diag{1,K,l}, = diag{1,K,l} и чисел ak [0,1] (k = 1,K,l), удовлетворяющих условию 2) теоремы 1, частотное неравенство (2), где матрица K( p) определяется формулой (4), выполнено для = 0 и всех 0 > 0. Тогда система (3) не имеет периодических решений второго рода частоты 0.
Теорема 5 применяется в параграфе 3.5 к определению областей отсутствия периодических решений второго рода для систем ФАПЧ с ПИФ и запаздыванием в петле обратной связи.
Четвертая глава диссертации целиком посвящена задаче о числе проскальзываний циклов в фазовой системе. Рассматриваются как сосредоточенные, так и распределенные фазовые системы. В параграфе 4.1 рассматривается сосредоточенная фазовая система со скалярной нелинейностью. Ее математическое описание имеет вид & z = Az + b(), (c,b ) (5) = c*z + (). Rm, R & Здесь A - m m гурвицева матрица, пара (A,b) - управляема, пара (A,c) -- наблюдаема.
Функция ( ) -- -периодическая непрерывно дифференцируемая функция с отрицательным средним и двумя нулями на периоде.
Определение 4. Говорят, что решение системы (5) проскальзывает k циклов, если существует такое число T > 0, что | (T ) - (0) |= k, и для всех t справедливо неравенство.
| (t) - (0) |< (k +1).
В параграфе 4.1 периодические функции Ляпунова, построенные в параграфе 2.2, используются для построения оценок числа проскальзываний циклов для решений системы (5).
Пусть 1 = ( ), 2 = sup[0,) ( ). Введем в рассмотрение функции inf [0,) -1 - ( ) = (1-1 ( ))(1-2 ( )), x j ( )d + (-1) k rj (k,, x) = ( j = 1,2), | ( ) | d x j ( )d + (-1) k r0 j (k,, x) = ( j = 1,2) ( ) |( ) | d и матричные функции arj (k,, x) 1 Tj (k,, x) = arj (k,, x) a0r0 j (k,, x), 2 0 a0r0 j (k,, x) где , , , , a0, a - параметры.
Построим матрицы A b O c Q =, L =, D =, O 0 1 z(t) d y = и функцию (t) = ( (t)).
( (t)) dt Рассмотрим квадратичную форму переменных y Rm+1, R G(y, ) = 2y*H (Qy + L) + y*DD*y + y*LD*y -(D*y -1-1 )* (2 - D*y) + y*LL*y, * где H = H - (m +1) (m +1) -матрица, , , , - параметры.
Выполнение частотного неравенства * -Re{K(i) - | K(i) |2 -(K(i) + 1 i) (6) -(K(i) + i )}- 0, где - K( p) = c*(A - pE) b - , гарантирует существование такой матрицы H, что выполнено неравенство G(y, ) 0, y Rm+1, R. (7) Теорема 6 Пусть существует такое число a [0,1], положительные , , , число 0 и натуральное k, что выполняются условия:
1) для всех 0 справедливо частотное неравенство (6);
2) матрицы Tj(k,, y*(0)Hy(0) - I) ( j = 1,2), где a0 = 1- a, I = y*(t)Hy(t), inf tR+ * являются положительно определенными для некоторой H = H, удовлетворяющей (7).
Тогда для решения (5) при всех t 0 справедлива оценка | (t) - (0) |< k. (8) Параграф 4.1 содержит также несколько модификаций теоремы 6.
В параграфе 4.2 результаты параграфа 4.1 распространяются на решения интегродифференциального уравнения t & (t) = (t) + ( (t - h)) - (t - )( ( ))d (t 0;h 0), (9) у которого функции 0(t) и (t) обладают теми же свойствами, что функции 0 j и в jk главе 3, а функция ( ) описана в параграфе 4.1.
& Передаточная функция линейной части уравнения (9) от входа к выходу (- ) имеет вид:
K( p) = -e- ph + (t)e- ptdt ( p C). (10) Теорема 7 Пусть существуют такие положительные числа , , , , число a [0;1] и натуральное число k, что выполняются условия:
1) для всех 0 справедливо частотное неравенство (6), где K( p) определяется формулой (10);
2) матрицы Tj (k,,Q), где a0 = 1- a, а величина Q вычисляется с помощью специальной процедуры, являются положительно определенными.
Тогда для любого решения уравнения (9) при всех t 0 выполнена оценка (8).
В параграфе 4.3 описана процедура применения теоремы 7 к системе ФАПЧ с ПИФ и синусоидальной характеристикой фазового детектора. Для числа проскальзываний циклов приводятся конкретные оценки, выраженные через параметры системы.
В параграфах 4.4 и 4.5 результаты параграфов 4.1 и 4.2 распространяются соответственно на многомерные и бесконечномерные фазовые системы с векторной нелинейностью.
Заключение.
В диссертации решены следующие задачи:
1. Построено обобщение периодических функций Ляпунова для фазовых систем дифференциальных уравнений. Разработан новый класс функционалов Попова для исследования асимптотического поведения интегродифференциальных систем с периодическими нелинейными функциями.
2. С помощью новых периодических функций Ляпунова и функционалов Попова получены новые многопараметрические частотные критерии глобальной асимптотической устойчивости математических моделей сосредоточенных и распределенных систем фазовой синхронизации.
3. Получены новые многопараметрические частотные критерии отсутствия у математических моделей фазовых систем периодических решений второго рода.
4. На основе частотных условий глобальной асимптотической устойчивости установлен ряд частотных оценок числа проскальзываний циклов для систем фазовой синхронизации.
5. Разработаны численные методы и алгоритмы проверки частотных критериев для систем фазовой автоподстройки частоты (ФАПЧ), содержащих запаздывание в петле обратной связи.
6. Разработан комплекс программ для реализации проверки частотных критериев для систем ФАПЧ с пропорционально-интегрирующим фильтром.
Публикации автора по теме диссертации Работы, опубликованные в изданиях из перечня ВАК:
1. Смирнова В.Б., Утина Н.В., Шепелявый А.И., Перкин А.А. Частотные оценки числа проскальзываний циклов в фазовой системе с векторной нелинейностью // Вестник СПбГУ сер.1, вып.1, 2009, C.33-43.
2. Смирнова В.Б., Утина Н.В., Шепелявый А.И., Перкин А.А. Покоординатные оценки векторного выхода многомерных систем с фазовым управлением. // Вестник СПбГУ, сер.1, вып.3, 2009, C.70-78.
3. Перкин А.А., Смирнова В.Б., Утина Н.В., Шепелявый А.И. О применении метода периодических функций Ляпунова // Вестник СПбГУ, сер.1, вып.3, 2011, C.31-40.
4. Перкин А.А. Устойчивость фазовых систем управления с дифференцируемыми нелинейностями // Гироскопия и навигация, №2, 2011, C.110.
Другие публикации по теме диссертации:
5. Перкин А.А., Смирнова В.Б. Частотные критерии устойчивости интегродифференциальных уравнений с периодическими нелинейными функциями. // ''Актуальные проблемы современного строительства.'' 57-я международная научнотехническая конференция молодых ученых. Сборник докладов Ч.I;, С.-Петербург.
СПбГАСУ. 2004, C.43-6. Перкин А.А., Смирнова В.Б. Глобальная асимптотика систем фазовой синхронизации // ''Актуальные проблемы архитектуры, строительства и транспорта. 58-я международная научно-техническая конференция молодых ученых. Сборник докладов'' Ч.II; С.-Петербург. СПбГАСУ. 2005, C.206-210.
7. Перкин А.А., Смирнова В.Б. Оценка областей устойчивости систем с фазовым управлением // ''Сборник докладов победителей конкурса грантов 2005г.'' С.-Петербург.
СПбГАСУ. 2006, C.75-103.
8. Смирнова В.Б., Перкин А.А. Оценка частоты периодических режимов II рода систем с фазовым управлением // ''Актуальные проблемы современного строительства.'' 59-я международная научно-техническая конференция молодых ученых. Сборник докладов. Ч.III; С.-Петербург. СПбГАСУ. 2006, C.121-131.
9. Смирнова В.Б., Утина Н.В., Шепелявый А.И., Перкин А.А. Об асимптотических свойствах многомерных фазовых систем с векторным управлением // X международный семинар им. Е.С. Пятницкого ''Устойчивость и колебания нелинейных систем управления.'' Тезисы докладов. Москва. ИПУ РАН. C. 271-210. Smirnova V.B., Shepeliavyi A.I., Utina N.V., Perkin A.A. The problem of cycleslipping for multidimensional phase control systems. // Proceedings of Sixth Euromech Nonlinear Dymanics Conference. 2008, St.Petersburg, Russia Смирнова В.Б., Утина Н.В., Шепелявый А.И., Перкин А.А. О задаче проскальзывания циклов для многомерных систем с векторным фазовым управлением. // Труды VIII Всероссийской научной конференции ''Нелинейные колебания механических систем''. Издательский дом ''Диалог культур'' т. I, Нижний Новгород, 2008, C. 314-318.
12. Smirnova V.B., Perkin A.A. Multiparametric frequency-domain criterion for stability of distributed systems with multiple equilibria. // PAMM, volume 7, Issue 1, p.20300452030046; 13. Перкин А.А., Смирнова В.Б. Многопараметрические оценки числа проскальзываний циклов в системах фазовой синхронизации. // ''Актуальные проблемы современного строительства.'' 61-я международная научно-техническая конференция молодых ученых. Сборник материалов. Ч.IV; С.-Петербург. СПбГАСУ. 2008, с.100-106.
14. Smirnova Vera, Perkin Alexey, Shepeliavyi Alexander. Frequency-domain criteria for gradient-like behavior of phase control system with vector nonlinearities. // 3-rd IEEE Multi-Conference on Systems and Control. Saint-Petersburg, July 2009. Papers. P.142-115. Перкин А.А., Смирнова В.Б. Частотные оценки решений фазовых систем с дифференцируемыми векторными нелинейностями // Доклады 66-й научной конференции профессоров, преподавателей, научных работников, инженеров и аспирантов СПбГАСУ. Ч. II, С.-Петербург, СПбГАСУ, 2009, C.147-152.
16. Alexey A. Perkin, Vera B. Smirnova, Alexander I. Shepeliavyi, Natalia V. Utina The employment of periodic Lyapunov function for asymptotic analysis of multidimensional phase control systems // Proceedings of 4-th International Scientific Conference of Physics and Control. Catania, 2009. 17. Перкин А.А., Смирнова В.Б., Утина Н.В., Шепелявый А.И.
Многопараметрические частотные критерии в теории устойчивости многомерных фазовых систем. // XI международная конференция ''Устойчивость и колебания нелинейных систем управления (конференция Пятницкого)''. Тезисы докладов. Москва.
ИПУ РАН. 2010 C.315-318. Перкин А.А., Смирнова В.Б. Устойчивость сингулярно-возмущенных уравнений фазовой синхронизации. // Доклады 67-й научной конференции профессоров, преподавателей, научных работников, инженеров и аспирантов СПбГАСУ. Ч. IV, С.Петербург, СПбГАСУ, 2010, C.188-192.
19. Перкин А.А., Смирнова В.Б., Шепелявый А.И. Прямой метод Ляпунова в исследовании асимптотики фазовых систем. // Труды ИПММ НАН Украины, т.21, 2010, C.177-120. Перкин А.А. Частотные условия отсутствия у многомерной фазовой системы предельных циклов II рода. // Доклады 68-й научной конференции профессоров, преподавателей, научных работников, инженеров и аспирантов СПбГАСУ. С.Петербург, СПбГАСУ, Ч. 4, 2011, C.108-112.
21. Alexey A. Perkin, Vera B. Smirnova, Alexander I. Shepeliavyi. Gradient-like properties of distributed and discrete phase systems. // Proceedings of 5-th International Scientific Conference of Physics and Control. Leon, 2011.
Перкин А.А. Устойчивость фазовых систем управления с дифференцируемыми нелинейностями.// Навигация и управление движением. Материалы XIII конференции молодых ученых, С.-Петербург, 2011, C.305-323. Перкин А.А., Перьева Е.Л., Смирнова В.Б., Шепелявый А.И.
Многопараметрические частотные оценки числа проскальзываний циклов для фазовых систем с дифференцируемыми нелинейностями. // Сб. тезисов Четырнадцатой международной научно-технической конференции ''Моделирование, идентификация, синтез систем управления''. Москва-Донецк, 2011. С. 16-17.
24. Перкин А.А. Комплекс программ. ''Частотный анализ сосредоточенных и распределенных систем фазовой автоподстройки частоты с пропорциональноинтегрирующим фильтром'' // Свидетельство о государственной регистрации программы для ЭВМ № 20116192 Авторефераты по всем темам >> Авторефераты по техническим специальностям