CC BY
43
Бифуркации ^^^^^^^^^^^^^^^^^
^ш в динамических системах
Изв. вузов «ПНД», т. 18, № 3, 2010 УДК 621.391.01
БИФУРКАЦИИ И КОЛЕБАТЕЛЬНЫЕ РЕЖИМЫ В СЛОЖНОЙ СИСТЕМЕ С ФАЗОВЫМ УПРАВЛЕНИЕМ
В.П. Пономаренко
Приведены результаты исследования режимов динамического поведения автогенераторной системы с фазовой автоподстройкой частоты и автоматическим регулированием усиления, проведенного на основе математической модели с двумя степенями свободы в цилиндрическом фазовом пространстве. Установлено расположение областей параметров с различными динамическими состояниями системы, выделена область устойчивости синхронного режима, изучены процессы, развивающиеся в области неустойчивости этого режима. Выяснено, что изменением параметров инерционности цепей управления, начальной частотной расстройки и степени связи через управляющие сигналы можно стимулировать в системе генерацию разнообразных периодических и хаотических колебаний.
Ключевые слова: Системы с фазовым управлением, динамические состояния, устойчивость, бифуркации, аттракторы, синхронный и несинхронные режимы, динамический хаос.
Введение
В последнее время проявляется заметный интерес к исследованию свойств динамического поведения автогенераторных систем, использующих принцип фазового управления частотой генерируемых колебаний. Область научно-технических приложений таких систем достаточно широка и включает в себя, в частности, решение таких задач, как стабилизация частоты колебаний, оптимальный прием, синхронизация и обработка сигналов, создание новых коммуникационных технологий и др. [1-5]. С позиции теории колебаний и нелинейной динамики модели систем с фазовым управлением представляют собой особый класс нелинейных динамических систем с цилиндрическим фазовым пространством. Эта особенность систем с фазовым управлением объединяет их с такими объектами, как физические маятники, электрические машины, вибрационные механизмы, системы со сверхпроводящими переходами и многими другими объектами и системами из различных областей естествознания и техники. Все эти объекты имеют сходные математические модели,
следовательно, свойства динамического поведения, которые демонстрируют системы с фазовым управлением, будут наблюдаться и в динамике других объектов, динамические модели которых определены в цилиндрическом фазовом пространстве. По этой причине математические модели систем с фазовым управлением во многих случаях рассматривают как базовые модели объектов, состояние которых определяется циклическими координатами. Привлекательным свойством систем с фазовым управлением является легко реализуемая возможность эффективного воздействия на свойства и области существования возникающих динамических состояний с помощью изменения параметров цепи фазового управления.
В связи с большим научным интересом и прикладной значимостью систем с фазовым управлением изучение процессов и явлений нелинейной динамики в разнообразных вариантах таких систем продолжает оставаться актуальной задачей. Данная работа посвящена исследованию динамических режимов и бифуркаций в модели автогенераторной системы с фазовой автоподстройкой частоты (ФАП) и автоматическим регулированием усиления (АРУ) в цепи управления частотой колебаний генератора. Происхождение и традиционное применение систем с такой структурой связано с решением задач передачи и обработки информации в радиотехнических системах со сложными сигналами, которые получают все большее распространение в радиоавтоматике, радионавигации и радиосвязи. В решении этих задач системы ФАП-АРУ применяются как схемные реализации оптимальных алгоритмов следящей оценки изменяющихся параметров радиосигналов с амплитудной модуляцией [4,6,7]. При этом система ФАП обеспечивает совпадение частот входного сигнала и колебаний управляемого генератора, а петля АРУ производит оценку непостоянной амплитуды входного сигнала и регулирует в соответствии с этой оценкой коэффициент усиления в цепи автоподстройки частоты, уменьшая тем самым влияние изменения амплитуды сигнала на качество работы системы ФАП.
Систему ФАП-АРУ можно представлять как вариант адаптивной системы фазовой автоподстройки, в которой адаптация осуществляется следящей подстройкой коэффициента усиления петли фазового управления. Интересно отметить, что в [8] вариант системы фазовой автоподстройки с автоматической регулировкой усиления в цепи обратной связи рассматривается в качестве радиофизической интерпретации обобщенной модели модифицированного генератора с инерционной нелинейностью, для которой характерна неединственность состояний равновесия.
Динамика системы ФАП-АРУ в простейшем случае использования в цепях управления фильтров низких частот первого порядка исследована в [9], где выяснено, что система демонстрирует регулярное поведение и имеет два стационарных режима: синхронный режим - режим слежения за параметрами входного сигнала и периодический несинхронный режим. В фазовом пространстве соответствующей динамической модели режиму слежения отвечает устойчивое состояние равновесия, а несинхронному режиму - устойчивый предельный цикл вращательного типа (охватывающий фазовый цилиндр). В пространстве параметров выделены области: глобальной устойчивости синхронного режима, одновременного существования синхронного и несинхронного режимов, существования несинхронного режима. В данной работе мы рассмотрим особенности динамики системы в случае, когда в цепях
управления подсистемы ФАП и петли АРУ используются фильтры соответственно второго и первого порядка. Такие фильтры интересны тем, что они соответствуют моделям динамики оцениваемых параметров принимаемого сигнала, используемым в ряде актуальных прикладных задач синтеза оптимальных систем приема и обработки [4,6]. Обращение к исследованию динамики системы ФАП-АРУ вызвано, с одной стороны, стремлением выяснить новые явления, возникающие в системе при усложнении фильтра в цепи фазового управления и обязанные своим происхождением потере устойчивости режима слежения, а с другой стороны, наблюдаемой сейчас тенденцией к расширению функциональных возможностей систем с фазовым управлением за счет использования различных несинхронных режимов в качестве рабочих процессов систем (в частности, попытками создания на этой основе новых коммуникационных систем с хаотическими сигналами [10-12]). Организация статьи следующая. В первом разделе обсуждается математическая модель исследуемой системы. Второй раздел посвящен анализу устойчивости и условий реализации синхронного режима. В третьем разделе рассмотрены особенности сложной динамики системы в области неустойчивости режима слежения. Четвертый раздел содержит результаты исследования двухчастотных несинхронных режимов. В Заключении сформулированы полученные результаты.
1. Модель системы
Будем рассматривать один из наиболее интересных вариантов системы ФАП-АРУ, в которой имеется связь между подсистемами через управляющие сигналы [7,9]. Уравнения, описывающие динамику такой системы, получаем из выведенных в [7] уравнений для оценок A* и ф* амплитуды A и фазы ф входного сигнала. При этом в связи со сложностью рассматриваемой задачи будем считать амплитуду A постоянным параметром, принимающим произвольные положительные значения, и не будем учитывать влияние помех. В этом случае удобно ввести фазовое рассогласование ф=ф—ф* и отношение амплитуд x=A*/А, тогда уравнения динамики рассматриваемой системы в операторной форме (p=d/dt) можно представить в виде [9]
рц>/к\=у — K1(p)x sin ф, x=S0 + к2K2(p)(cos ф — x + ax sinф). (1)
В уравнениях (1) у=рф/к1 - начальная расстройка частот, к\ - полоса удержания системы ФАП, So=Ao/A, Ao - среднее значение амплитуды входного сигнала, к2 - коэффициент преобразования петли АРУ, a - степень связи через управляющие сигналы, K1(p) и K2(p) - операторные коэффициенты передачи фильтров низких частот в цепях управления подсистем. Типы этих фильтров определяются используемыми моделями изменения оцениваемых параметров ф и A [4,6].
Для рассматриваемых в данной работе фильтров с коэффициентами передачи Ki(p) = 1/[(1 + Tip)(l + T2P)] и K2(p) = 1/(1 + Тзр), где Ti, T2 и T3 - постоянные времени, из уравнений (1) получаем следующую математическую модель динамики
системы ФАП-АРУ:
dф
dT = y,
dy = z,
dT
dz
ц— = у — x sin ф — £iz — y,
dT
dx
£2—t = 0 — px + cos ф + ax sin ф,
dT
где t = kit, 0 = So/fc2, в = 1 + 1/k2, £1 = (Ti + Тг)кь £2 = Тз^/к2, ц = TiT2k1. Система уравнений (2) в силу периодичности правых частей по фазовой переменной ф с периодом 2п является динамической системой с четырехмерным цилиндрическим фазовым пространством U = {ф(шоё 2п), y, z, x} и семимерным пространством параметров Л = {у, 0, в, a, £i, £2, [}. Параметры 0, в, [х, £i, £2 положительны по физическому смыслу. В силу инвариантности системы (2) относительно замены (у, a, ф) ^ (—у, —a, — ф) достаточно рассматривать ее при значениях a > 0.
При исследовании динамики системы ФАП-АРУ на основе математической модели (2) в данной работе рассмотрены следующие две задачи. Первая задача связана с определением условий существования и реализации режима слежения, в котором система осуществляет автоматическую подстройку частоты и амплитуды колебаний управляемого генератора к аналогичным параметрам входного сигнала. Эта задача представляет наибольший интерес для традиционного применения систем ФАП-АРУ как устройств, реализующих оптимальные алгоритмы приема сигналов с непостоянной амплитудой. Вторая задача имеет целью получить представление о возможных типах несинхронных режимов и сценариях их эволюции при изменении параметров модели (2). Эта задача представляет интерес для проблемы расширения функциональных возможностей систем с фазовым управлением. В фазовом пространстве U режиму слежения отвечает устойчивое состояние равновесия модели (2), несинхронным (квазисинхронным и асинхронным) режимам соответствуют устойчивые предельные циклы, притягивающие торы и хаотические аттракторы колебательного, вращательного и колебательно-вращательного типа.
Квазисинхронные режимы, определяемые аттракторами колебательного типа модели (2), характеризуются непостоянством разности фаз ф, которая изменяется в ограниченном диапазоне значений. В асинхронных режимах, определяемых аттракторами вращательного и колебательно-вращательного типа модели (2), наблюдается либо возрастание (убывание) разности фаз ф, то есть имеет место вращение ф (режимы вращательного типа), либо чередование стадий неограниченного изменения (вращения) и колебательного поведения ф в ограниченном диапазоне значений (режимы колебательно-вращательного типа). В несинхронных режимах имеют место периодические или хаотические автомодуляционные колебания управляемого генератора системы ФАП-АРУ. В [7, 8] обоснована предпочтительность колебательно-вращательных хаотических асинхронных режимов для получения в системе с фазовым управлением широкополосных колебаний с хаотической фазой и равномерной в полосе спектральной плотностью.
Важно отметить, что и в решении первой задачи исследование несинхронных режимов имеет принципиальное значение - бифуркации, в результате которых появляются несинхронные режимы, определяют границы областей параметров модели (2), в которых обеспечивается существование и глобальная устойчивость режима слежения в исследуемой системе. Таким образом, знание особенностей поведения системы в несинхронных режимах и способов воздействия на развитие колебательных процессов имеет большое значение для правильного выбора режимов и значений параметров системы в различных конкретных приложениях. Решение поставленных задач проведено при помощи качественно-численных методов нелинейной динамики [3, 13, 14] и компьютерного моделирования с использованием программного комплекса «ДНС» [15].
2. Устойчивость режима слежения
Полагая в системе (2) dф/dт = dy/dx = dz/dT = dx/dT = 0, получаем следующие уравнения для определения координат состояний равновесия
Y — x sin ф = 0, 6 — ex + cos ф + ax sin ф = 0, y = 0, z = 0. (3)
На плоскости (ф, x) в полосе —п < ф < п первое уравнение в (3) определяет кривую Xi (ф) = у/ sin ф, второе уравнение - кривую Х2(ф) = (6 + cos ф)/(в — a sin ф); точки пересечения кривых Х^ф) и Х2(ф) определяют значения ф и x, соответствующие состояниям равновесия. Анализ взаимного расположения кривых Х^ф) и Х2(ф) показывает, что в зависимости от соотношения между параметрами у, 6, в и a может быть либо четыре, либо две, либо ни одной точки их пересечения. Границы областей существования состояний равновесия, определяемые из условий касания кривых Х^ф) и Х2(ф), находим из равенств
Yi(6, в, a) = sin3 фтг/(в cos фт), (4)
в которых ф^(г = 1, 2,..., 6) - координаты точек касания кривых Х^ф) и Х2(ф) -определяются из уравнения
в cos ф(6 + cos ф) — sin2 ф(в — a sin ф) = 0. (5)
Уравнение (5) при значениях a > в имеет два корня ф = фт1 и ф = фт2 (—п/2 < фт1 < 0, —п < фт2 < —п/2), когда 0 < 6 < \/1 — в2^2; четыре корня ф = фm1, ф = фт2, ф = фт3 и ф = фт4 ( —п/2 < фт1 < 0, —П < фт2 < —п/2, п/2 < фтз < фт4 < п), когда л/1 — в2^2 < 6 < 1; два корня ф = фт1 и ф = фт3 (—п/2 < фт1 < 0, п/2 < фт3 < п), когда 6 > 1. При значениях a < в уравнение (5) имеет четыре корня ф = фт1, ф = фт2, ф = фт5 и ф = фт6 ( — п/2 < фт1 < 0, —п < фт2 < —п/2, 0 < фт5 < п/2, п/2 < фт6 < п), когда 6 < 1, и два корня ф = фт1 и ф = фт5, когда 6 > 1. При этом фт1 соответствует касанию кривых Х1(ф) и Х2(ф) в области x > 0 при значениях y < 0; фт2 - в области x < 0 при значениях y > 0; фт3, фт4 и фт6 - в области x < 0 при значениях y < 0; фт5 - в области x > 0 при значениях Y > 0.
На рис. 1 представлен качественный вид кривых Yi, определяемых из (5) и выделяющих на плоскости (y, 6) области с различным числом состояний равновесия
Рис. 1. Разбиение плоскости параметров (у, 6) на области с различным числом состояний равновесия модели (2) при значениях 0 < а < в (а) и а > в (б)
при значениях 0 < а < в (рис. 1, а) и а > в (рис. 1, б). Для значений у и 6 из области Ci модель (2) имеет четыре состояния равновесия: x\, 0, 0), А2(ф2, x2, 0, 0),
Аз(фз,жз, 0, 0) и А4(ф4,Х4, 0, 0), причем Ai и А2 расположены в области x > 0, а Аз и А4 - в области x < 0 фазового пространства U. При значениях у и 6 из области Co модель (2) имеет два состояния равновесия А1 (ф1, xi, 0, 0) и А2(ф2, x2, 0, 0), при значениях у и 6 из области C2 - два состояния равновесия Аз(фз^з, 0, 0) и А4(ф4, x4, 0, 0).
Характер состояний равновесия устанавливаем из анализа коэффициентов характеристического уравнения
X4 + [(в^ - а^ sin ф- + £1£2)/(£2^)]Хз + [(в£1 - ав1 sin ф- + £2)/(£2^)]^2+
2 (6) + [(в — а sin ф- + £2xj cos ф-)/(£2^)]X + ^x- cos ф- — sin2 ф-)/(£2^) = 0,
где j = 1, 2, 3, 4. Состояния равновесия А1 и Аз могут быть как устойчивыми, так и неустойчивыми, состояния равновесия А2 и А4 - неустойчивые седлового типа. Условия устойчивости состояний равновесия А1 и Аз находим, применяя к уравнению (6) критерий Рауса-Гурвица, и записываем в виде
в^ — а^ sin ф- + £1£2 > 0, в£1 — а£1 sin ф- + £2 > 0,
в — а sin ф- + £2x- cos ф- > 0, вxj cos ф- — sin2 ф- > 0,
(в — а sin ф- + £2x- cos ф-)[(вц — а^ sin ф- + £1£2) х (7)
х (в£1 — а£1 sin ф- + £2) — (в — а sin ф- + £2x- cos ф-)£2^] —
— (в^ — а^ sin ф- + £1£2)2(вxj cos ф- — sin2 ф-) > 0, j = 1,3.
При выполнении условий (7) в системе ФАП-АРУ существуют стационарные режимы слежения, определяемые состояниями равновесия А1 и Аз. В режиме состояния равновесия А1 колебания управляемого генератора системы ФАП с точностью до величины ф1 синфазны, а в режиме состояния равновесия Аз - с точностью до величины фз противофазны с колебаниями на входе системы. Основным режимом слежения является режим состояния равновесия А1, так как в этом режиме выше точность оценки параметров входного сигнала, характеризуемая величинами ф1 и x1, которые определяются из уравнений (3). Из уравнений (3) следует, что при значениях 6 = 6s(y, в, а), где 6s = в — ау — л/1 — у2, в основном режиме слежения
Рис. 2. Бифуркационная диаграмма для модели (2), разделяющая плоскость параметров (ц, у) на области, соответствующие различным динамическим режимам, при значениях параметров а = 2, в = 1.1, б = 1.25, £1 = 1, £2 =2
осуществляется точная оценка амплитуды входного сигнала (х1 = 1) с величиной фазового рассогласования ф1 = агевт у; при б > бз величина XI > 1, при б < бз величина х1 < 1.
Уравнения (3) и (5), равенство (4) и неравенства (7) позволяют определять границы области Со существования основного режима слежения. На рис. 2 приведены кривые у = у1 и у = уз на плоскости параметров (ц, у), выделяющие область Со, построенные при значениях параметров а = 2, в = 1.1, б = 1.25, £1 = 1, £2 = 2. Линия у1, определяемая из (3)-(5), является границей области существования состояний равновесия А1 и А2. Линия уз - это граница области устойчивости состояния равновесия А1, определяемой условиями (7). Область С0 : {у1 < у < уз} является областью удержания режима слежения. Область Со разделяется кривой ур, соответствующей бифуркации петли сепаратрисы Пф седло-фокуса А2 (И,еА,1)2<0, 1шХ1>2 = 0, А,3 < 0, Х4 > 0, где ^1)2,з,4 -корни характеристического уравнения (6) для состояния равновесия А2), на области Сз : {шах(у1; ур) < у < уз} и Сз1 : {у1 < у < шш(ур, уз)}. Для значений ц и у из области Сз состояние равновесия А1 является единственным (по результатам численного исследования модели (2)) аттрактором в фазовом пространстве и. Следовательно, область Сз - это область захвата в основной режим слежения, для значений параметров из этой области режим слежения устанавливается в системе при любых начальных условиях.
При переходе с увеличеним ц (или с уменьшением у) через границу ур в область С31 из петли Пф появляется устойчивый (так как седловая величина о = А.4 + И,еА,г < 0, г = 1, 2, 3) 2п-периодический по ф вращательный предельный цикл Ьо. Численное исследование модели (2) показывает, что при движении в области С31 с изменением ц и у возможно исчезновение цикла Ьо в результате седло-узловой бифуркации и возникновение хаотического аттрактора вращательного типа через перемежаемость. На рис. 3, а, б приведены (ф, у)-проекции фазовых портретов и фрагменты временных реализаций колебательных процессов у(т), соответствующие режимам цикла Ьо и хаотического аттрактора. Установлено также существование сложных многооборотных вращательных предельных циклов. Таким образом, в области Сз1 модель (2) демонстрирует бистабильное поведение - в фазовом пространстве и одновременно существуют устойчивое состояние равновесия А1 и аттракторы вращательного типа. Следовательно, в зависимости от начального состояния системы в ней реализуется либо синхронный, либо периодический или хаотический асинхронный режим. Для реализации синхронного режима начальные значения фазовых переменных должны принадлежать области притяжения состоя-
Рис. 3. Несинхронные режимы системы при значениях у = —0.54, ц = 1.6 (а), 2.63 (б); у = 0.86, И, = 0.727 (в), 1.56 (г), 0.74 (д), 1.24 (е) (а = 2, в = 1.1, б = 1.25, £г = 1, £2 = 2)
ния равновесия А\. Области Со и С3 служат основными динамическими характеристиками системы ФАП-АРУ при решении задачи следящей оценки изменяющихся параметров входного сигнала. Подбором величины и знака параметра связи а можно изменять величины этих областей.
3. Особенности динамики системы в области неустойчивости режима слежения
Для значений параметров модели (2) вне области Со режим слежения не существует, стационарными режимами системы являются несинхронные (квазисинхронные и асинхронные) режимы, которые могут быть как регулярными, так и хаотическими. Вначале рассмотрим качественные изменения в поведении системы ФАП-АРУ, связанные с потерей устойчивости режима слежения, наблюдаемой при выходе из области Со. Для этого снова обратимся к двухпараметрической бифуркационной диаграмме (ц, у), приведенной на рис. 2. Будем использовать следую-
щие обозначения аттракторов модели (2): Si (г = 1, 2, 3, 4, 5, 6) - колебательные предельные циклы; Р) (] = 1, 2, 3, 4) - колебательные хаотические аттракторы; Ьк (к = 1, 2, 3, 4) - вращательные и колебательно-вращательные предельные циклы; Шг (I = 1, 2, 3, 4, 5, 6, 7) - хаотические аттракторы асинхронных режимов; Тт (т = 1, 2, 3, 4, 5) - двумерные торы.
При пересечении с увеличением ц границы уз области устойчивости Со в модели (2) наблюдается суперкритическая бифуркация Андронова-Хопфа, когда в решении характеристического уравнения (6) появляется пара комплексно-сопряженных корней с положительной реальной частью. В результате этой бифуркации в фазовом пространстве и мягко рождается устойчивый предельный цикл Sl, охватывающий состояние равновесия А1, то есть происходит мягкий переход системы ФАП-АРУ от режима слежения к периодическому квазисинхронному режиму. На рис. 3, в изображены (ф, у)-проекция фазового портрета и зависимость у(т), соответствующие режиму колебательного предельного цикла Sl.
Режим цикла Sl существует в области Сро, заключенной между кривой уз и частями кривых ус1, у^ь ус2, ус3 и у^2 (см. рис. 2). Кривые ус1, ус2 и ус3 соответствуют седло-узловым бифуркациям, а кривые у^1 и у^2 - бифуркациям удвоения периода колебательных предельных циклов модели (2). Бифуркационная кривая ур разделяет область Сро на области Ср1 и Ср2 : {уз < у < шт(ур, уа2)}. В области Ср1 предельный цикл Sl глобально устойчивый, то есть квазисинхронный режим, соответствующий этому циклу, реализуется в системе при любых начальных условиях. Область Ср2 - это область с бистабильным поведением модели (2). При значениях параметров из этой области в фазовом пространстве и одновременно существуют колебательный предельный цикл и аттракторы вращательного типа, которые могут быть как регулярными, так и хаотическими. Хаотизация колебаний осуществляется через бифуркации удвоения периода вращательных предельных циклов, а также через перемежаемость, как и в области Сз1.
С седло-узловыми бифуркациями колебательных предельных циклов, происходящих на линиях уС2 и усз, связано существование области Са на плоскости (ц, у) (см. рис. 2), ограниченной частями кривых ус2, ус3, у^1 и у^2, для точек которой в фазовом пространстве и, кроме предельного цикла S1, существует еще один предельный цикл колебательного типа S2. Соответствующие этому циклу (ф, у)-проекция фазового портрета и зависимость у(т) приведены на рис. 3, г. Квазисинхронные режимы системы ФАП-АРУ, определяемые циклами Sl и S2, различаются амплитудами колебаний переменных х и у. В области Са также имеет место бистабильное поведение модели (2), какой из этих режимов реализуется в системе - это зависит от начальных условий. При переходе из области Са в область Ср1 через кривые ус2 и ус3 происходит исчезновение одного из предельных циклов S1 и S2. Таким образом, область С* = Ср1 и Са, граничащая с областью Сз, - это область периодических квазисинхронных режимов системы ФАП-АРУ.
Теперь рассмотрим перестройки поведения системы, находящейся в квазисинхронном режиме, которые происходят при пересечении границ области С* с изменением параметра ц. Характер этих перестроек оказывается различным в зависимости от того, через какую из бифуркационных границ ус1, у^1 и у^2 происходит выход из области С*. При переходе через кривую уС1 (см. рис. 2) квазисинхронный режим исчезает, в системе наблюдается жесткий переход к периодическому асинхронному режиму колебательно-вращательного предельного цикла Ь1. На рис. 3, д даны
(ф, у)-проекция фазового портрета и зависимость у(т), соответствующие этому режиму. При пересечении кривой у^ вначале система мягко переходит к квазисинхронному режиму предельного цикла который появляется в результате бифуркации удвоения периода цикла 51, а затем следует жесткий переход от режима цикла 52 к периодическому асинхронному режиму колебательно-вращательного типа. При переходе через кривую у^2 и последующем удалении от нее наблюдается иное поведение системы - реализуется последовательность бифуркаций удвоения периода колебательного предельного цикла, которая завершается возникновением хаотического аттрактора Р1 в фазовом пространстве и. Таким образом, в этом случае имеет место мягкий переход системы к хаотическому квазисинхронному режиму аттрактора Рь Примеры соответствующих такому режиму (ф, у)-проекции фазового портрета и зависимости у(т) приведены на рис. 3, е.
Рассмотрим динамическое поведение модели (2) в области Сс, расположенной на диаграмме (ц, у) (см. рис. 2) справа от области Ср0 (Ср0 = С* и Ср2). Область Сс - это область со сложной динамикой модели (2). Для значений параметров из этой области в системе реализуются квазисинхронные и асинхронные режимы различной сложности. При этом в областях С^1 и С^2 (С^,С^2 С Сс) существуют и периодические квазисинхронные режимы предельных циклов 51 и 52 соответственно. Отметим, что правая граница области С^2 имеет сложный характер, она образована бифуркационными кривыми, соответствующими бифуркациям удвоения периода и седло-узловым бифуркациям предельного цикла Б2.
При исследовании модели (2) в области Сс в качестве варьируемого параметра был выбран параметр инерционности ц, который характеризует степень влияния фильтра второго порядка в цепи фазового управления. С целью получения представления о возможных несинхронных режимах системы ФАП-АРУ, развивающихся в области Сс, обратимся к однопараметрическим бифуркационным диаграммам (ц, у) (рис. 4, а) и (ц, х) (рис. 5, а), а также к (ф,у)- и (ф, х)-проекциям фазовых портретов и зависимостям у(т) и х(т), соответствующим аттракторам модели (2) при различных значениях параметра ц (рис. 4, б-н, рис. 5, б-г).
Диаграмма (ц, у) на рис. 4, а построена при значении расстройки у = 0.625; она характеризует режимы поведения и бифуркационные переходы, наблюдаемые при выходе из области С^ через границу у^1 и последующем удалении от этой границы. Эта диаграмма отображает эволюцию одновременно существующих квазисинхронных режимов предельных циклов 51 и 52 при изменении параметра ц от 1.16 до 2.1. Для динамики системы характерны переходы между периодическими и хаотическими несинхронными режимами. При этом для большей части исследуемого интервала значений ц в системе наблюдаются хаотические асинхронные режимы. Можно выделить следующие три характерные области значений параметра ц. В области М1 : {1.16 < ц < 1.423} чередующиеся асинхронные режимы сложных предельных циклов и хаотических аттракторов колебательно-вращательного типа существуют одновременно с квазисинхронным периодическим или хаотическим режимом, то есть модель (2) демонстрирует бистабильное поведение. При переходе в область М2 : {1.423 < ц < 1.676} в системе наблюдается жесткий переход от асинхронного режима хаотического аттрактора (рис. 4, б) к квазисинхронному режиму хаотического аттрактора Р1 (рис. 4, в).
В области М2 поведение системы характеризуется преобразованием режима аттрактора Р1 через периодический квазисинхронный режим предельного цикла 53
Рис. 4. Диаграмма развития квазисинхронных режимов предельных циклов 51 и Я2 при изменении параметра ц (а), проекции фазовых портретов и временные реализации у(т) для различных несинхронных режимов (б-н) при значениях ц=1.422 (б); 1.47 (в); 1.5 (г); 1.6 (д); 1.63 (е); 1.99 (ж); 2.24 (з); 2.74 (и); 3.07 (к); 4.15 (л); 4.6 (м,н) (у = 0.625, а = 2, в = 1.1, б = 1.25, £1 = 1, £2 = 2)
(рис. 4, г) в асинхронный режим колебательно-вращательного хаотического аттрактора (рис. 4, д). Примечательно, что с увеличением ^ происходит уменьшение числа витков вращательных траекторий в структуре аттрактора далее наблюдается явление перемежаемости типа «хаос-хаос». Это явление, иллюстрируемое рис. 4, е, характеризуется нерегулярным чередованием на хаотическом аттракторе Ш2 коротких стадий вращательных движений и длинных стадий колебательных движений. После исчезновения вращательных движений на аттракторе ^2 система возвращается к квазисинхронному режиму хаотического аттрактора Р\, который затем мягко преобразуется в периодический квазисинхронный режим предельного цикла 52.
При значениях ^ £ Мз : > 1.676} система, наряду с асинхронными режимами, демонстрирует квазисинхронные режимы иного типа, определяемые колеба-
Рис. 5. Развитие асинхронного режима предельного цикла Ь2 при изменении параметра (а) и несинхронные режимы системы при значениях |Л=1.37 (б); 1.405 (в); 1.54 (г) (у = 0.86, а = 2, в = 1-1, б = 1-25, £х =1, £2 = 2)
тельными предельными циклами и хаотическими аттракторами модели (2) с ограниченным диапазоном изменения разности фаз ф, превышающим 2п. Некоторые из таких аттракторов, соответствующих так называемым квазисинхронным режимам с проворотом фазы, изображены на рис. 4, ж,и,к. На рис. 4, з,л,м,н даны примеры аттракторов асинхронных режимов, которые сменяют квазисинхронные режимы при изменении ц. Переходы к хаотическим квазисинхронным режимам происходят через бифуркации удвоения периода колебательных предельных циклов с проворотом фазы. Хаотические асинхронные режимы возникают на базе хаотических квазисинхронных режимов через перемежаемость типа «хаос-хаос».
Диаграмма (ц, х) на рис. 5, а соответствует значению у = 0.86 и отражает развитие периодического асинхронного режима, на который система переходит в результате пересечения границы уС1 области С*. В соответствии с этой диаграммой, где стартовым состоянием системы является периодический асинхронный режим колебательно-вращательного предельного цикла Ь2 (рис. 5, б), динамика системы при изменении параметра ц от 1.37 до 1.56 характеризуется следующими явлениями: переходом к асинхронному режиму колебательно-вращательного хаотического аттрактора Ш3 (рис. 5, в) через бифуркации удвоения периода цикла Ь2; нерегулярным чередованием периодических и хаотических режимов колебательно-вращательного типа; жестким переходом от режима колебательно-вращательного хаотического аттрактора (рис. 5, г) к квазисинхронному режиму предельного цикла Б2 (см. рис. 3, г). При дальнейшем увеличении ц режим цикла 52 жестко сменяется режимом колебательно-вращательного хаотического аттрактора, область существования которого прерывается узкими «окнами» сложных предельных циклов. Отметим, что и в рассматриваемом случае хаотические асинхронные режимы колебательно-вращательного типа наблюдаются для большей части исследуемого диапазона изменения параметра ц.
Результаты, представленные на рис. 4 и 5, свидетельствуют о сложной картине перестроек динамического поведения модели (2) при движении внутри области Сс. Численное исследование модели (2) позволило установить, что закономерности ее динамического поведения, представленные на рис. 2-5 (потеря устойчивости синхронного режима, бифуркации периодических движений, переходы к хаотическому поведению по различным сценариям, нерегулярное чередование квазисинхронных и асинхронных режимов различной сложности, существование разнообразных колебательно-вращательных предельных циклов и хаотических аттракторов в фазовом пространстве и, преобразование регулярных колебательно-вращательных аттракторов в хаотические и обратно, изменение соотношения колебательных и вращательных движений в структуре аттракторов, существование колебательных аттракторов с проворотом фазы), характерны и для других значений параметров а, в, б, е1 и £2. На основании этого можно полагать, что результаты, приведенные на рис. 2-5, являются типичными для динамического поведения системы ФАП-АРУ.
Кроме этих результатов, при исследовании модели (2) обнаружено бистабиль-ное поведение системы при изменении параметра связи а, выяснено, что при больших значениях параметра инерционности £2 наблюдаются новые явления сложной динамики модели (2), связанные с возможностью существования притягивающих двумерных торов в фазовом пространстве и, которым соответствуют двухчастот-ные несинхронные режимы системы ФАП-АРУ. Обсудим качественные особенности процессов возникновения и развития таких режимов.
4. Двухчастотные несинхронные режимы
При численном исследовании модели (2) выявлена существенная роль параметра связи а в возникновении и развитии двухчастотных квазисинхронных режимов. Рассмотрим случай, когда значения параметров в = 1.1, 8 = 1.25, у = 0.1, ц = 5, е1 = 1.97, е2 = 150, а параметр а изменяется от 1.6 до 2.6. На рис. 6 и 7 приведены построенные в этом случае бифуркационная диаграмма {а, х} и (ф, х)-проекции фазовых портретов, (у, х)-проекции сечения Пуанкаре и временные реализации х(т), соответствующие аттракторам модели (2). Начальным состоянием системы при а = 1.6 является квазисинхронный режим предельного цикла Б4 (рис. 7, а). При увеличении а цикл Б4 теряет устойчивость, и в фазовом пространстве рождается притягивающий двумерный тор Т колебательного типа (рис. 7, б), которому в отображении Пуанкаре соответствует замкнутая инвариантная кривая Г1 (рис. 7, в). Тору Т1 соответствует двухчастотный квазисинхронный режим системы ФАП-АРУ.
1 .1 7
2.09 а 2.6
Рис. 6. Эволюция квазисинхроного режима предельного цикла £4 при изменении параметра связи а (у = 0.1, ц = 5, в = 1.1, 8 = 1.25, ех = 1.97, £2 = 150)
1.88 ф 3.55 -1.9 ф 3.53 1.0 у 1.04 -2.0 ф 3.65
■3.14 ф 9.42 -2.0 ф 3.67 0.39 у 0.49
Рис. 7. Аттракторы несинхронных режимов при значениях параметра а = 1.6 (а); 1.65 (б,в); 1.9 (г); 1.946 (д); 2.08 (е,ж,з); 2.097 (и,к); 2.35 (л); 2.436 (м,н) (у = 0.1, ц = 5, в = 1.1, б = 1.25, £1 = 1.97, £2 = 150)
Далее в интервале 1.623 < а < 1.92 имеет место чередование режима тора Т1 и режимов резонансных многооборотных колебательных предельных циклов, один из которых - девятиоборотный цикл S5 приведен на рис. 7, г. Начиная со значения а = 1.92, наблюдается искажение формы замкнутой инвариантной кривой Г1, что свидетельствует о постепенном преобразовании тора Т1 в хаотический аттрактор Р2 (рис. 7, д) типа тор-хаос [13]. Режим аттрактора Р2 прерывается режимами многооборотных колебательных предельных циклов, из которых наибольшую область существования имеет восьмиоборотный цикл S6. При дальнейшем увеличении а режим цикла S6 через бифуркации удвоения периода преобразуется в режим колебательного хаотического аттрактора Р3 (рис. 7, е,ж,з). При а > 2.088 в структуре аттрактора Р3 появляются витки вращательных фазовых траекторий, что свидетельствует о преобразовании его в колебательно-вращательный хаотический аттрактор Ш5 (рис. 7, и,к).
При изменении а в интервале 2.09 < а < 2.432 динамика системы характеризуется чередованием хаотических аттракторов и сложных предельных циклов колебательно-вращательного типа, которые отличаются числом витков колебательных траекторий. На рис. 7, л приведен один из таких циклов - цикл Ь3 с девятью витками колебательных движений. Хаотизация колебаний в режимах этих циклов осуществляется через бифуркации удвоения периода. Когда а превышает значение 2.432, система через перемежаемость «хаос-хаос» переходит от режима колебательно-вращательного хаотического аттрактора к квазисинхронному режиму хаотического аттрактора Р4 (рис. 7, м,н). С дальнейшим увеличением а режим аттрактора Р4 преобразуется через режим предельного цикла, подобного циклу S5 (рис. 7, г), в квазисинхронный двухчастотный режим тора Т2, подобного тору Т1
(рис. 7, б,в). Затем режим тора Т2 мягко преобразуется в квазисинхронный режим предельного цикла, подобного циклу Б4 (рис. 7, а) - система возвращается к исходному квазисинхронному режиму.
Существование двумерных вращательных торов в фазовом пространстве модели (2) установлено в случае, когда значения параметров в = 1.1, 8 = 1.25, а = 5, ц = 5, £\ = 1, е2 = 150, а расстройка у является варьируемым параметром. На рис. 8 и 9 приведены примеры (ф, у)-проекций фазовых портретов, (у, х)-проекций сечений Пуанкаре и реализаций колебательных процессов у(т), отвечающих различным регулярным и хаотическим режимам системы ФАП-АРУ при изменении у от 0.85 до 0.32. При у = 0.85 в фазовом пространстве и существует двухоборотный (4п-периодический по ф) вращательный предельный цикл Ь4 (рис. 8, а). При уменьшении у происходит потеря устойчивости цикла Ь4 через бифуркацию рождения устойчивого двумерного вращательного тора Т3 (рис. 8, б), которому соответствуют замкнутые инвариантные кривые Г 2 и Г 2 отображения Пуанкаре (рис. 8, в). Тору Т3 отвечает двухчастотный асинхронный режим системы.
При дальнейшем уменьшении у наблюдаются следующие явления: чередование режима тора Т3 и периодических асинхронных режимов на многооборотных вращательных предельных циклах; постепенное вырождение тора Т3 в хаотический аттрактор типа тор-хаос (рис. 8, г,д); чередование хаотических режимов (рис. 8, е,ж,з) и периодических режимов двухоборотного, четырехоборотного (8п-периодического по ф) (рис. 8, и) и трехоборотного (6п-периодического по ф) (рис. 8, к) вращательных предельных циклов. Затем наблюдается эффект потери
Рис. 8. Аттракторы асинхронных режимов при значениях у = 0.853 (а); 0.845 (б,в); 0.819 (г,д); 0.818 (е,ж,з); 0.8038 (и); 0.65 (к) (ц = 5, а = 5, в = 1.1, 8 = 1.25, ех = 1.97, £2 = 150)
Рис. 9. Асинхронные режимы, развивающиеся на базе трехоборотного вращательного предельного цикла при уменьшении расстройки у: у = 0.646 (а); 0.643 (б); 0.636 (в); 0.625 (г); 0.45 (д); 0.362 (е); 0.323 (ж); 0.32 (з,и,к) (ц = 5, а = 5, в = 1.1, 8 = 1.25, £х = 1, £2 = 150)
устойчивости трехоборотного цикла с рождением в фазовом пространстве и притягивающего двумерного тора Т4, которому соответствуют замкнутые инвариантные кривые Г3, Г2 и отображения Пуанкаре (рис. 9, а). Дальнейшее движение по параметру у приводит к разрушению тора Т4 и переходу к режиму хаотического аттрактора (рис. 9, б,в). Область существования аттрактора прерывается областями, в которых существуют периодические асинхронные режимы вращательных предельных циклов (рис. 9, г,д,е). При этом переходы между периодическими и хаотическими режимами осуществляются через прямые и обратные бифуркации удвоения периода, через перемежаемость, а также в результате седло-узловых бифуркаций предельных циклов. Далее на базе вращательного предельного цикла, (ф, у)-проекция фазового портрета которого приведена на рис. 9, е, снова возникает двухчастотный асинхронный режим тора Т5 (рис. 9, ж), который затем также преобразуется в хаотический асинхронный режим (рис. 9, з,и,к).
Заключение
В данной работе приведены результаты изучения режимов поведения и явлений нелинейной динамики в автогенераторной системе с фазовой автоподстройкой частоты и автоматическим регулированием усиления, полученные с помощью качественно-численного исследования соответствующей математической модели (2) -динамической системы с двумя степенями свободы в цилиндрическом фазовом пространстве. Установлено, что модель (2) может иметь аттракторы различной сложности - от состояния равновесия до хаотических аттракторов колебательного, вращательного и колебательно-вращательного типа, которым соответствуют синхронный,
квазисинхронные и асинхронные режимы исследуемой системы. Получены условия существования синхронного режима, определяющие границу области Со удержания этого режима.
Выяснено, что потеря устойчивости синхронного режима сопровождается мягким возникновением периодического квазисинхронного режима, в то время как исчезновение синхронного режима приводит к жесткому переходу системы к периодическому или хаотическому асинхронному режиму. Внутри области удержания Со выделены область Cs глобальной устойчивости синхронного режима (область захвата) и область Csi, соответствующая сосуществованию синхронного и асинхронных режимов различной сложности. Границы области захвата определяются бифуркациями модели (2), связанными с потерей устойчивости состояния равновесия, образованием петли сепаратрисы седло-фокуса и слиянием устойчивого и седлового состояний равновесия. Проанализировано изменение поведения системы при переходе через границы областей захвата и удержания. В области неустойчивости синхронного режима выделены область периодических квазисинхронных режимов, примыкающая к области захвата, и область со сложной динамикой модели (2), в которой у системы существуют периодические и хаотические несинхронные режимы колебательного, вращательного и колебательно-вращательного типа. Обнаружена сложная картина перестройки динамического поведения системы при движении с изменением параметров внутри области со сложной динамикой, свидетельствующая о неоднородной структуре этой области.
Сравнительный анализ результатов, полученных при исследовании модели (2), с результатами исследования динамики отдельной системы фазовой автоподстройки частоты с фильтром второго порядка, полученными в [12], показывает, что многие динамические режимы и явления, обнаруженные в такой системе (потеря устойчивости синхронного режима, бифуркационные механизмы возникновения несинхронных режимов, квазисинхронные режимы с проворотом фазы), свойственны и системе ФАП-АРУ. Это свидетельствует о сильном влиянии свойств поведения цепи фазового управления на динамику системы ФАП-АРУ. Наряду с этим установлена возможность существования в системе ФАП-АРУ двухчастотных квазисинхронных и асинхронных режимов, обусловленная наличием цепи автоматического регулирования амплитуды. Выяснено, что в процессах возникновения и развития таких режимов значительную роль играют параметр связи и начальная частотная расстройка. Примечательно, что, в отличие от отдельной системы ФАП, система ФАП-АРУ демонстрирует различное поведение в зависимости от знака начальной частотной расстройки у. Это связано с отсутствием у модели (2) свойства симметрии структуры разбиения пространства параметров относительно замены у ^ —у.
Для применения системы ФАП-АРУ в задачах следящей оценки параметров принимаемого сигнала наибольший интерес представляют результаты о границах области захвата Cs, которая является характеристикой устойчивости реализуемого системой ФАП-АРУ алгоритма оптимальной обработки по отношению к вариациям параметров системы, и результаты об изменении характера динамики системы при переходе границ областей захвата и удержания. Эти результаты также дают представление об особенностях поведения системы в процессе вхождения в режим слежения и после нарушения этого режима. В плане расширения круга возможного применения рассматриваемой системы представляют интерес приведенные в работе результаты о развитии и преобразовании несинхронных режимов в области со
сложной динамикой модели (2) при изменении параметров инерционности, степени связи и начальной частотной расстройки. Эти результаты могут быть использованы при решении задач управления свойствами и областями существования различных несинхронных режимов с помощью изменения указанных параметров системы с целью ее перевода в желаемый стационарный режим.
Библиографический список
1. Шахгильдян В.В., Ляховкин А.А. Системы фазовой автоподстройки частоты. М.: Связь, 1972.
2. Линдсей В. Системы синхронизации в связи и управлении. М.: Сов. радио, 1978.
3. Системы фазовой синхронизации / Под ред. В.В. Шахгильдяна, Л.Н. Белюсти-ной. М.: Радио и связь, 1982.
4. Тихонов В.И., Кульман Н.К. Нелинейная фильтрация и квазикогерентный прием сигналов. М.: Сов. радио, 1975.
5. Дмитриев А.С., Панас А.И.Динамический хаос. Новые носители информации для систем связи. М.: Физматлит, 2002.
6. Ярлыков М.С. Применение марковской теории нелинейной фильтрации в радиотехнике. М.: Сов. радио, 1980.
7. Кульман Н.К., Жеронкина Н.Н.Помехоустойчивость оптимального приема квазигармонического процесса со взаимно коррелированными амплитудой и фазой // Радиотехника и электроника. 1969. Т. 14, № 11. С. 2050.
8. Постнов Д.Э., Баланов А.Г. Синхронизация в хаотических системах со счетным числом состояний равновесия // Изв. вузов. Прикладная нелинейная динамика. 1997. Т. 5, № 1. С. 69.
9. Пономаренко В.П. О режимах работы и области захвата системы фазовой синхронизации с цепью автоматического регулирования усиления (АРУ) // Радиотехника и электроника. 1986. Т. 21, № 10. С. 2023.
10. Дмитриев А.С., Широков М.Е. Выбор генератора для прямохаотической системы связи // Радиотехника и электроника. 2004. Т. 49, № 7. С. 840.
11. Дмитриев А.С., Клецов А.В., Кузьмин Л.В. Генерация сверхширокополосного фазового хаоса в дециметровом диапазоне // Радиотехника и электроника. 2009. Т. 54, № 7. С. 709.
12. Матросов В.В., Шалфеев В.Д.Динамический хаос в фазовых системах: Учебное пособие. Нижний Новгород: Изд-во Нижегородского университета, 2007.
13. Шильников Л.П., Шильников А.Л., Тураев Д.В., Чуа Л. Методы качественной теории в нелинейной динамике. Часть 2. Москва; Ижевск: НИЦ «Регулярная и хаотическая динамика», Институт компьютерных исследований, 2009.
14. Анищенко В.С. Сложные колебания в простых системах. М.: Наука, 1990.
15. Динамика нелинейных систем. Программный комплекс для исследования нелинейных динамических систем с непрерывным временем: Учебно-методическая разработка / Сост. В.В. Матросов. Н.Новгород: ННГУ, 2002.
НИИ Прикладной математики Поступила в редакцию 28.09.2009
и кибернетики ННГУ После доработки 23.03.2010
BIFURCATIONS AND OSCILLATORY MODES IN COMPLEX SYSTEM WITH PHASE CONTROL
V.P. Ponomarenko
The results are produced of research of dynamical modes and bifurcation in a complex system with phase control, based on mathematical model with two degrees of freedom in the cylindrical phase space. The location of domains corresponding to different dynamical states of the system is established. The processes developing in the system as a result of loss stability of the synchronous mode, and scenarios of evolution of nonsynchronous modes under variation of system parameters are investigated. The possibility to stimulate the generation of various periodic and chaotic oscillations changing the parameters of control loops inertia, coupling power and the initial frequency detuning is revealed.
Keywords: Systems with phase control, dynamic states, stability, bifurcation, attractors, synchronous and nonsynchronous modes, dynamical chaos.
Пономаренко Валерий Павлович - родился в Мурманске (1941), окончил Горьковский университет (1963). С 1968 года работает в Научно-исследовательском институте прикладной математики и кибернетики при Нижегородском университете (НИИ ПМК при ННГУ). Защитил диссертации на соискание ученой степени кандидата технических наук (1972) и доктора физико-математических наук (1993) в области теории синхронизации, теории колебаний и динамики систем. С 1991 года - заместитель директора НИИ ПМК при ННГУ. Область научных интересов - теория колебаний и волн, применение нелинейной динамики в радиоавтоматике синхронизация автоколебательных систем, математическое моделирование. Один из авторов монографий «Фазовая синхронизация» и «Системы фазовой синхронизации». Опубликовал много научных статей по направлениям, указанным выше. 603005, Н. Новгород, ул. Ульянова, 10 НИИ Прикладной математики и кибернетики при ННГУ E-mail: povp@uic.nnov.ru
