CC BY
50
УДК 621.373.1+517.938
В.В. Зайцев, И.В. Стулов, А.Н. Шилин1
СУБГАРМОНИЧЕСКАЯ СИНХРОНИЗАЦИЯ АВТОКОЛЕБАНИЙ В ДИСКРЕТНОМ ВРЕМЕНИ
Исследованы автоколебания в системе, осциллирующей в дискретном времени и находящейся под действием внешнего гармонического сигнала. В качестве объекта выбрано предложенное ранее авторами дискретное отображение осциллятора Ван дер Поля. Режимы синхронизации автоколебаний исследованы методами численного эксперимента, гармонического баланса и медленно меняющихся амплитуд. Установлено, что эффект субгармонической синхронизации в дискретном времени можно реализовать при произвольном отношении частот свободных автоколебаний и внешнего сигнала.
Ключевые слова: автоколебательная система, дискретное время, подмена частот, синхронизация, захват частоты.
1. Предварительные сведения
В физике колебаний известно явление субгармонической синхронизации (субгармонического захвата) [1]. Оно возникает при воздействии на автоколебательную систему с частотой свободных автоколебаний По внешнего гармонического сигнала с частотой Пя , лежащей в окрестности МПо, где М — целое число. При выполнении определенных условий автоколебания меняют свою частоту По на частоту Па = П8/М, т. е. происходит захват частоты автоколебаний. Механизм захвата — это взаимодействие между вынужденными колебаниями с частотой внешнего воздействия и гармониками автоколебаний [1].
В настоящем сообщении показано, что это же явление в системах, осциллирующих в дискретном времени, имеет свою специфику, обусловленную эффектом подмены частот (наложения). Эффект имеет место для гармоник автоколебаний, не удовлетворяющих теореме Котельникова. Они воспринимаются автоколебательной системой как подмененные [2]. Частота ш' подмененной гармоники с номером М определяется выражением
П' = П (МПа - ¡¡ааг(МПа)), ВД = { 1 -П° ^^ , (1)
х© Зайцев В.В., Стулов И.В., Шилин А.Н., 2015
Зайцев Валерий Васильевич (zaitsev@ssu.samara.ru), Стулов Игорь Валерьевич (stulov@mail.ru), кафедра радиофизики и полупроводниковой микро- и наноэлектрони-ки, Самарский государственный университет, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
Шилин Александр Николаевич (Shilax@yandex.ru), АО "РКЦ"Прогресс", 443009, Российская Федерация, г. Самара, ул. Земеца, 18.
где О' = ш'/ша, Оа = ша/ша, ша —частота автоколебаний, — частота дискретизации. В результате такой подмены субгармонический захват может наблюдаться и для частот, не находящихся в целочисленном отношении с частотой свободных автоколебаний. Эффект с участием подмененной М-й гармоники имеет место для автоколебаний с частотами Оа > 1/2М. На более низких частотах имеет место классический вариант субгармонической синхронизации.
2. Численный эксперимент с дискретным осциллятором Ван дер Поля
Для демонстрации субгармонической синхронизации в дискретном времени используем дискретное отображение осциллятора Ван дер Поля вида [3; 4]
y[n] - 2а cos(2nQ0)y[n - 1] + a2y[n - 2]= yF (y[n - 1],y[n - 2]) + x[n - 1]. (2)
Здесь Qo — нормированная на частоту дискретизации собственная частота линейного контура аналогового прототипа, y > 0 — параметр глубины положительной обратной связи, а — параметр диссипации, связанный добротностью Q прототипа: а = exp (-nQo/Q)- Внешнее гармоническое воздействие на осциллятор
x[n] = E cos(2nQsn)
характеризуется амплитудой и нормированной частотой Qs. Функция кубической нелинейности в правой части отображения (2) имеет вид
F (y[n - 1],y[n - 2]) = cos(2nQo)g(y[n - 1]) - g(y[n - 2]), g(y) = y - ys/3.
При превышении параметром y порога самовозбуждения Yst = 1-а2 ДВ-осцилля-
Рис. 1. Амплитудный спектр автоколебаний
тор (2) генерирует свободные автоколебания у[п] с амплитудным спектром А(О), изображенным на рис. 1 (для = 0.21, Q = 20 и Y = 0.15). При частоте автоколебаний Оа = 0.21 спектральные линии подмененных третьей (30) и пятой
(50) гармоник в полном соответствии с формулой (1) расположены на частотах П'3 = 0.37 и П5 = 0.05.
Приведенные ниже рисунки отображают результаты численного эксперимента в дискретной автоколебательной системе с параметрами По = 0.21, Q = 2° и 7 = 0.15. Подмененная частота третьей гармоники в этом случае П' = 0.37. Значение амплитуды сигнала синхронизации фиксировано Е = 0.25, а его частота Пя медленно возрастает в диапазоне П^ < Пя < Пн, содержащем П'. Генерируемый сигнал подвергнут полосовой фильтрации на центральной частоте По, и затем методом аналитического сигнала выделены огибающая и фаза основной гармоники автоколебаний. Полученные таким образом частото-частотная и амплитудно-
Рис. 2. Частото-частотная характеристика
частотная характеристики показаны на рис. 2 и 3 непрерывными линиями. В центре характеристик расположена область захвата частоты, а по краям — области биений. В целом характеристики типичны для явления захвата частоты. Но они зеркальны по отношению к характеристикам захвата на основном тоне.
3. Статические характеристики субгармонического захвата частоты
Детали механизма захвата частоты в исследуемой системе можно выяснить, записав установившееся решение разностного уравнения (2) в виде суммы
У[п] = Уа[п]+ ys[n], (3)
Здесь первое слагаемое
У a [n] = 2 Aexp (j2nüan) + k.c. = 1 AZn + k.c.
Рис. 3. Амплитудно-частотная характеристика
— это автоколебания с комплексной амплитудой А и частотой Оа, удовлетворяющей уравнению (см. (1))
а второе
^ (3Оа - ¡¡ООт(ЗОа)) = О3
Уэ[пП\ = ^Вехр (]2пОвп) + к.с. = ^BZn + к.с.
(4)
— вынужденные колебания с частотой Оя и комплексной амплитудой В в линейном контуре с частотной характеристикой Н(О):
В = Н (О3)Е =
Zя
Z2S - 2а cos(2пО0)Zs + а2
Е.
С учетом суммы (3) уравнение (2) принимает вид
Уа[п] - 2а сов(2пОо)уа[п - 1] + а2уа[п - 2} = 1Е (у[п - 1},у[п - 2]).
(5)
Методом гармонического баланса разностное уравнение (5) сведено к комплексному алгебраическому уравнению
(1 - 4А2) А = Н(Оа)(1
Н(ОоИ 1 -7АО ) А = Н(Оа) ( 1 - 7АГ - 7
^ ..2_ 1 В|2_ ЧА1В
4'"' 21 1 4 А
(6)
где Ао — амплитуда свободных автоколебаний с частотой Оо. Помимо очевидного нулевого решения А = 0, уравнение (6) имеет два решения вида
3
\А(Оа)\2 =4 - вг(Оа) (4 - АО) - 2\В\2 ±у/Б(Оа, В.
(7)
Здесь
ДП„, В) = (4 - вг(Па) (4 - АО) - 4\В\2 - (А(Па) (4 - АО))0
0) - 7ВП \В\2 - (вг(Па) (4 - АО))0 (8)
и введены обозначения
в(Па) = н(По)/н(Па), вг = Де(в), вг = Мв).
При малой расстройке частот — 2^\Па - По\ С 1 функция в(Па) приобретает простой вид:
в(Па) = 1+ 32Ц Па- По =1+ зп(Па). (9)
По
Для этого приближения решения (7) записываются как
\А(Па,Е)\0 = ао - 3\В\2 ± ^(Ао-^^^^^еТ-П^-А^. (10)
Анализ динамики процесса субгармонического захвата (см. ниже) показывает, что устойчивому решению соответствует знак плюс в формуле (10). График амплитудно-частотной характеристики, построенный по (10), показан на рис. 3 точками. Теоретическая зависимость \А(Па)\ количественно хорошо согласуется с экспериментальной. При этом точечный график зависимости Па(Пя), построенный по уравнению (4), на рис. 2 совпадает с прямолинейным участком экспериментального графика. Тем самым подтверждается адекватность механизма с участием подмененных гармоник автоколебаний исследуемому явлению захвата частоты в дискретном времени.
Режим субгармонического захвата может реализоваться при положительных значениях выражения (8). Граница области захвата определяется решением нелинейного уравнения Е(Па,В) = 0. При использовании приближения (9) для нее нетрудно получить простое выражение
Па(\В\) = По ^ 1 ± ЗД(4В- А0)УАО - 7\В\^ , П = ^ (3Па - ¡¡оат(ЗПа)). (11)
На рис. 4 граница (11) показана непрерывной линией. Зона режима субгармонического захвата — это область I.
4. Динамика субгармонического захвата частоты в приближении ММА
Для исследования динамики процесса субгармонической синхронизации воспользуемся методом медленно меняющихся амплитуд (ММА), адаптированным для применения к нелинейным системам с дискретным временем [4]. Комплексную амплитуду автоколебаний уа[п] в суперпозиции (3) будем считать медленной по сравнению с Zn функцией дискретного времени:
У а [п] = 1 А[п^ + к.С.
При этом комплексная амплитуда вынужденных колебаний у8 [п] по-прежнему считается постоянной. С физической точки зрения это означает, что обратная связь в системе включается после установления вынужденных колебаний в "холодном" контуре.
Методом ММА уравнение (5) преобразуется в укороченное уравнение для комплексной амплитуды:
A[n] = A[n - 1] - пQ(1 + jn)A[n - 1]+
2 i В 1 ,г и (12)
+ 2 | 1 - 4 A[n - 1]|2 - 1 |В|2 - 1 V ^ | A[n - 1]
Линеаризация уравнения (12) в окрестности одного из решений уравнения (6) приводит к следующему уравнению для отклонений комплексных амплитуд от их стационарных значений:
АА[п] = АА[п - 1] - пО0(1 + ]щ)АЛ[п - 1]+ + 2 (1 - 2|А|2 - 1 |Е|2) АЛ[п - 1] - 4 (Л*Б + 2А2) АА*[п - 1]. (13)
Основное назначение этого линейного уравнения — исследование устойчивости различных режимов автоколебаний.
Устойчивому решению Л = 0 соответствует режим захвата на основном тоне. При этом автоколебания с частотой По подавляются, и в системе происходят лишь вынужденные колебания с частотой Уравнение (13) в этом случае имеет наиболее простой вид: ( )
АА[п] = Р(Оа, \Е\)АЛ[п - 1],
Р(па, \Е\) = 1 - пО-(1+ „) + 2(1 - 1 |Е|2) . Отсюда следует, что на плоскости параметров (О.а, |Е|) граница области захвата
на основном тоне определяется решением уравнения P(Па, |—|) = 1:
-( Па)
^2 + 4 (l - пQ 1 - 4п2(П - П,^,
(14)
П = Ft(3Па - ¡1ааг(3Па)).
На рис. 4 граница (14) показана пунктирной линией. Из-за слабой частотной зависимости граница практически прямолинейна. Захват на основном тоне ns происходит в области III, лежащей выше линии (14). В областях II в системе наблюдаются биения между автоколебаниями и вынужденными колебаниями.
Рис. 5. Временные зависимости амплитуд
Эти выводы подтверждаются результатами моделирования переходных процессов с помощью укороченного уравнения (13). На рис. 5 представлены временные зависимости амплитуд автоколебаний при выборе параметров внешних воздействий из трех областей: I— —| = 0.75, Hs = 0.37, II — —| = 0.75, Hs = 0.379 и III— —I = 1.25, ns = 0.376. Линия I — это выход амплитуды на стационарный уровень в режиме захвата. В процессе установления наблюдается "ступенька" в окрестности неустойчивого значения (знак "минус" в формуле (10). Линия II — огибающая биений между сигналами автоколебаний и внешнего воздействия, линия III — релаксация отклонения амплитуды к устойчивому значению .
На диаграмме режимов автоколебаний (рис. 4) часть области I занимает область IV. Режим, реализующийся в этой области, зависит от предыстории системы. Если процесс изменения параметров воздействия соответствует переходу I ^ IV, то в области IV реализуется режим субгармонического захвата. Если же имеет место переход III ^ IV, то в области IV наблюдается захват на основном тоне ns.
Рисунок 6 иллюстрирует поведение системы в описанной ситуации. На нем изображены графики зависимости амплитуды автоколебаний |А| от амплитуды внешнего воздействия — |, построенные по формуле (10) при фиксированной частоте внешнего воздействия ns = 0.37. Непрерывная линия соответствует знаку
А|
1.67
1.33 1
0.67 0.33
°0 0.5 1* " |В|
Рис. 6. Зависимость амплитуды автоколебаний от амплитуды внешнего сигнала
плюс в формуле, пунктирная — знаку минус. На отрезке линии между точками Ы\ и М2 параметры системы находятся в области IV. При плавном изменении амплитуды \Б\ наблюдаются скачкообразные изменения амплитуды |А|. Скачок Т происходит при переходе между областями IV ^ III, скачок Т2 — при переходе III ^ I. Переходы между областями I ^ IV и IV ^ I в точке Ы\ не сопровождаются скачком амплитуды автоколебаний.
Заключение
При анализе воздействия внешнего гармонического сигнала на осциллятор Ван дер Поля установлено, что явление субгармонической синхронизации в дискретном времени можно реализовать при произвольном отношении частот синхронизированных автоколебаний и синхронизирующего сигнала. В этом варианте субгармоническую синхронизацию предлагается использовать для дробного деления и умножения частоты дискретных сигналов.
Литература
[1] Хаяси Т. Нелинейные колебания в физических системах. М.: Мир, 1968. 432 с.
[2] Оппенгейм А., Шафер Р. Цифровая обработка сигналов. М.: Техносфера, 2006. 856 с.
[3] Зайцев В.В., Линдт С.В., Шилин А.Н. Осцилляторы Ван дер Поля и Рэлея в дискретном времени // Вестник Самарского государственного университета. 2014. № 7(118). С. 104-114.
[4] Зайцев В.В. О дискретных отображениях осциллятора Ван дер Поля // Физика волновых процессов и радиотехнические системы. 2014. Т. 17. № 1. С. 35-40.
References
[1] Hayash C. Nonlinear oscillations in physical systems. M., Mir, 1968, 423 p. [in Russian].
[2] Oppenheim A., Schafer R. Digital signal processing. M., Tekhnosfera, 2006, 856 p. [in Russian].
[3] Zaitsev V.V., Lindt S.V., Shilin A.N. Van der Pol and Rayleigh oscillators in discrete time. Vestnik Samarskogo gosudarstvennogo universiteta [Vestnik of Samara State University], 2014, no. 7(118), pp. 104-114 [in Russian].
[4] Zaitsev V.V. About discrete mapping of the Van der Pol oscillator. Fizika volnovykh protsessov i radiotekhnicheskie sistemy [Physics of Wave Processes and Radio Systems], 2014, Vol. 17, no. 1, pp. 35—40 [in Russian].
V.V. Zaitsev, I.V. Stulov, A.N. Shilin2
SUBHARMONIC SYNCHRONIZATION OF SELF-OSCILLATIONS IN DISCRETE TIME
Self-oscillations in the system oscillating in discrete time and being under the influence of an external harmonious signal are investigated. As an object of researches the option of discrete map of the oscillator of Van der Pol offered earlier by authors is brought to attention. The modes of synchronous self-oscillations are analysed by means of methods of numerical experiment and harmonic balance. It is established that the effect of subharmonic synchronization in discrete time can be realized at any relation of frequencies of synchronized self-oscillations and the synchronizing signal.
Key words: self-oscillatory system, discrete time, substitution of frequencies, synchronization, capture of frequency.
Статья поступила в редакцию 28/VJ//2015. The article received 28/VT7/2015.
2Zaitsev Valeriy Vassilievich (zaitsev@ssu.samara.ru), Stulov Igor Valerevich (stulov@mail.ru), Department of Radiophysics and Semiconductor Micro- and Nanoelectronics, Samara State University, 1, Acad. Pavlov Street, Samara, 443011, Russian Federation.
Shilin Alexandr Nikolaevich (Shilax@yandex.ru), JSC SRC Progress, Samara,18, Zemetsa Street, Samara, 443009, Russian Federation.
