CC BY
21
Известия Института математики и информатики УдГУ
2015. Вып. 2 (46)
УДК 62-54.14, 681.511.4 © О. Г. Антоновская, В. И. Горюнов
МОДЕЛИРОВАНИЕ ПРОЦЕССА СКЛЕЙКИ ФАЗОВЫХ ТРАЕКТОРИЙ В СИСТЕМЕ С КОМБИНИРОВАННЫМ ЧАСТОТНО-ФАЗОВЫМ УПРАВЛЕНИЕМ
В работе излагаются основные моменты применения методики исследования динамики систем с широтно-импульсной модуляцией управляющего сигнала при учете его кусочного постоянства к исследованию динамики системы синхронизации с комбинированной широтно-импульсной модуляцией управляющего сигнала. Методика основана на применении метода точечных отображений и состоит в использовании для рассмотрения динамики системы фазовых подпространств, соответствующих постоянству управляющего сигнала. В частности, предложена процедура моделирования процесса перехода от режима фазового управления к режиму частотного управления в системе с комбинированным частотно-фазовым управлением.
Ключевые слова: управляемые системы, математическое моделирование, метод точечных отображений.
Введение
Известно, что в системе с комбинированным частотно-фазовым управлением (КИЧФД) [1] контур с фазовым управлением функционирует в случае, если импульсы синхронизирующего опорного сигнала (ОС) чередуются во времени с выходными импульсами счетчика (С). Частотное управлением в системе с КИЧФД включается, если процесс чередования импульсов ОС и С нарушается. При появлении во времени подряд двух и более импульсов С принудительно устанавливается минимально возможное значение сигнала управления, а при появлении подряд двух и более импульсов ОС — его максимально возможное значение.
Вследствие принципиальной нелинейности процессов в системе с КИЧФД в режиме фазового управления [2] разработка процедуры моделирования процесса перехода от режима фазового управления к режиму частотного управления, необходимая для построения математической модели системы управления в целом, является актуальной.
§ 1. Использование метода точечных отображений для построения математической модели системы
Для разработки процедуры моделирования процесса перехода от режима фазового управления к режиму частотного управления в системе с КИЧФД воспользуемся опытом математического моделирования процессов в системах с широтно-импульсной модуляцией и кусочно-постоянным характером сигнала управления [3]. С этой целью введем в рассмотрение в качестве безразмерной фазовой координаты относительную степень заполнения С так, что после заполнения фазовая координата в становится равной единице, после чего происходит мгновенный сброс на значение в = 0, и при этом на выходе С появляется управляющий импульс. Тогда динамика изменения координаты в в безразмерном времени т, нормированном на период следования импульсов ОС, описывается уравнением [4]
а(йв/йт) = д(п(т)) (0 < т < 1, 0 < в < 1). (1.1)
т
в
в = 1 в = 0
появляется управляющий импульс; а — безразмерный показатель С, причем значение а = 1 соответствует неуправляемому ПГ, когда и(т) = 0 д = 1 + 5и — линейная характеристика ПГ с крутизной 5 > 0.
Сущность использования метода точечных отображений для изучения динамики СЧ с КИЧФД состоит в использовании для рассмотрения динамики СЧ фазовых подпространств, соответствующих постоянству структуры дифференциального уравнения.
Согласно логике работы КИЧФД режим частотного детектирования (ЧД) с минимально возможным значением и(т) = и включается в режиме имиульсно-фазового детектирования (ИФД) при появлении повторного импульса С, а с максимальным значением и(т) = й — при появлении повторного импульса ОС. При чередовании импульсов С и ОС в режиме ФД ближайший приходящий во времени импульс ОС запускает генератор пилообразного напряжения так, что на его выходе появляется сигнал управления:
и(т) = (и - и)т + и (0<т<1). (1.2)
При использовании импульсного фазового детектора типа «выборка-запоминание» и(т) = = и(то) = сопв^ где то — момент прихода импульса С. Поскольку в режимах ЧД и ИФД структура уравнения (1.1) неизменна, можно считать, что при и(т) = и уравнение (1.1) определено в подпространстве Щ, при и(т) = Л — в подпространстве Пз, в режиме ИФД с и(т), определяемой выражением (1.2), — в подпространстве П2.
Рассмотрим свойства фазовых траекторий в каждом из подпространств.
Согласно (1.1) траектории движения в Щ на развертке тора 0 ^ т ^ 1, 0 ^ в ^ 1 суть прямые с наклоном
с1т/с16 = а/а, (1.3)
где а = д(и).
Необходимо отметить, что, поскольку КИЧФД переходит в режим ЧД с и(т) = и только при появлении повторного импульса С, анализ поведения траекторий уравнения (1.3) следует проводить только для начальных условий во = 0, 0 ^ то ^ 1. Рассмотрим свойства траекторий при различных значениях параметра а. Начнем со случая а < а, когда траектории в(т) на развертке тора расположены полого с йт/йв < 1. В момент достижения изображающей точкой
в=1 в=0
т = 1 т = 0
после чего ИТД продолжает движение.
Непосредственно из отмеченного следует, что при а < а в режиме ЧД в основном появляются импульсы С и одиночные импульсы ОС. Поскольку появление подряд двух импульсов ОС исключено, при а < а ИТД остаются в подпространстве Щ, то есть при любых начальных условиях СЧ не выходит из режима ЧД.
При а > а траектории движения в Щ по сравнению со случаем а < а на развертке тора идут круче биссектрисы координатного угла, так что становится возможным появление подряд двух импульсов ОС, необходимых для перехода КИЧФД в режим ФД. Согласно (1.3)
(0 , то )
в точке (в1,1), всегда соответствует импульсу ОС, причем
01 = (о/а)(1-то) (О<01<1). (1.4)
(в1, 0)
(в2, 1) в2 < 1
02 = (а/а)(2-го) (0 < 02 < 1). (1-5)
Из (1.5) следует, что в2 < 1 при т0 > т01 (т01 < 1). Это означает, что по координате т0 существует граница
Гц ты = 2-(а/а). (1.6)
В полуполосе а ^ а, 0 ^ то ^ 1 граница Г1 определяет две области. При начальном значении т0, лежащем в области, расположеиной левее Г1, КИЧФД остается в режиме ЧД, и при этом траектория движения порождает точечное отображение Тц с функцией исследования (ФИ):
Тп: г = то + (се/а) — 1 (0 ^ т0 ^ ш, а < а < 2а). (1.7)
При начальном значении то, лежащем в области, расположен ной правее Гх, КИЧФД переходит в точке (02) 0) в режим ИФД, и при этом до момента появления импульса С, то есть до достижения траекторией движения границы в = 1, осуществляется переходное движение с и(т) = и. Тем самым порождается переходное отображение Т\2 с ФП:
Т\2 : т = то + (а/а) — m (roi ^ то ^ 1, а < а < 2а; 0 ^ то < 1, а ^ 2а), (1.8)
где m = E[tq + (а/а)}, Е — знак взятия целой части, ат)2.
В подпространстве Пз, соответствующем режиму ЧД с и(т) = и, согласно (1.1) траектории движения на развертке тора 0 ^ т ^ 1, 0 ^ в ^ 1 суть прямые с наклоном
с1т/с1в = а/а, (1.9)
где а = g (и).
Необходимо отметить, что, поскольку КИЧФД переходит в режим ЧД с и(т) = и только при появлении повторного импульса ОС, анализ поведения траекторий уравнения (1.9) следует проводить только для начальных условий 0 ^ во ^ 1,то = 0.
Проводя анализ поведения траекторий аналогично тому, как это было сделано для случая и(т) = и, можно установить, что при а > а появление подряд двух импульсов С исключено, так что ИТД остаются в подпространстве Пз, то есть СЧ не выходит из режима ЧД.
При â/2 < а <а ИТД в зависимости от величины во либо остаются в Пз, либо переходят в подпространство П2, соответствующее режиму ИФД. При а < а/2 ИТД при любом во за
Пз
Т33 : в = в0- (а/а) (0 < в0 < в01, â/2 < а < а), (1.10)
где граница Гх задается уравнением
Г2: 0oi = 2 -(â/a) (â/2 < а <а) (1.11)
и отображением
Î32 : т = (а/а)( 1 - в0 + Е[в0 + â/a]) - 1 ЭХ < во < 1, а/2 < а < а; 0 < в0 < 1, а < а/2).
(1.12)
Переход из Пз в П2 реализуется при появлении подряд двух импульсов С.
В режиме ИФД детектор с «выборкой и запоминанием» в интервале следования во времени двух соседних импульсов С в моменты tq и г обеспечивает выполнение условия и(т) = u(tq) = = const. Поэтому, согласно (1.1), траектории движения в П2 определяются соотношением
в(т) = (д(хо)/а)(т - т0) (т0^т^т). (1.13)
На развертке фазового тора в режиме ИФД траектория (1.13) с учетом цилиндричности по координате т состоит из двух отрезков прямых с одним промежуточным переносом ИТД по координате т при в = вз. Поэтому условиями существования траектории (1.13) в режиме ИФД являются неравенства в(т = 1), вз < 1. Разрешая указанные неравенства относительно параметра а, неодим, что их можно представить в виде неравенств а > аз(то) и а < а4(то) соответственно, где
аз (то) = 5(и(то))(1 - то) (0 < то < 1), (1.14)
а4(то) = ^(и(то))(2 - то) (0 < то < 1). (1.15)
Границы а = а3, а = а4 в иолупол осе 0 ^ то ^ 1, а> 0 образуют три обл асти: Gi, G2, G3. Для точ ек (то, а) средней обл асти G2 режим ИФД сохраняется и после появления в момент т = т импульса С. Фазовая траектория (1.13) при этом порождает отображение
Т22: т = т0 + (а/д(и0)) - 1 (т0,т € П2, а3(г0) < а < а4(т0)). (1.16)
Если то на плоскости то, а лежит в области Gi, соответствующей неравенству а < аз(то), то реализуется случай двух импульсов ОС и КИЧФД переходит в режим ЧД с и(т) = и, порождая отображение
Т21: т = т0 + (а/д(и0)), (1.17)
и при этом точка во = 0, tq = т становится начальной точкой движения в Щ.
Если то на плоскости то, а лежит в области G3, соответствующей неравенству а > а^(то), то реализуется случай двух импульсов ОС и КИЧФД переходит в режим ЧД с и(т) = й, порождая отображение
Т23: в = (1/а)д(ио)(2-то) (т0 € П2, в € П3, а > а4(т0)). (1.18)
При этом точка во = в, tq = 0 становится начальной точкой движения в П3.
§2. О возможности качественного исследования поведения системы
Отображение T22, определяющее динамику СЧ в режиме ИФД, имеет единственную неподвижную точку. Подставляя в (1.16) условие tq = т = т*, находим, что для линейной характеристики (1.2) ИГ
а = а(т*) = д(и(т*)) = 1 - S[(û - и)т* (2.1)
График функции а = а(т*) (2.1) удовлетворяет при 0 ^ т* ^ 1 системе неравенств аз(т*) < < а < а4(т*), и, следовательно, неподвижная точка то = т* существует во всей полосе удержания, то есть при а < а < а.
Неравенства г* > 0, г* < 1 и g (и) > 0 определяют на плоскости a, S треугольную область существования синхронного режима и невырожденности (д(и) > 0) характеристики ПГ:
1 — al — а „ 1 , .
0 < ——,—— < S < --(и<0,й>0). 2.2
\и\ \и\ \и\
Из вида ФИ (1.16) следует, что условие устойчивости неподвижной точки
т0 = r*|dr/dro(ro =т*)| < 1
нарушается на границе
N-: S = 2а/(й-и), (2.3)
которая отсекает от треугольной области существования (2.2) ее левую верхнюю часть (в сторону уменьшения а и увеличения S). Тем самым определена область существования и устой-
то = т*
С целью реализации возможности качественного анализа устойчивости в «большом» непо-то = т* а, S > 0
существования и невырожденности границы 71 : а = 2а, 72 : а = таха:3(то), 73 : а = min а^(то).
то то
Тогда в каждой из трех образующихся подобластей области существования и невырожденности оказывается определена конкретная группа точечных отображений. Так, в части области существования над границей Yi (в сторону больших значений S) достаточно воспользоваться представлением графиков отображений T23, T33, T32 и T22. При этом отображепия T23, T33, T32 следует представлять на одной комбинированной диаграмме. В первом квадранте диаграммы следует отождествлять ось абсцисс с осью во. Тогда в интервале 0 ^ во ^ 0oi следует представить график отображения Т33, причем величины в следует указывать на оси ординат. В интервале doi ^ во ^ 1 следует представить график ФИ отображения Т32, причем величины т следует указывать на оси ординат. В четвертом квадранте следует представить график ФП T23. то
в отображения — указывать на оси абсцисс. При таком подходе предоставляется возможность качественного формирования графика ФП произведения T32T33T23, где k — степень, определяемая при построении.
Например, при a, a3max ^ а ^ а, 2а график ФП Т22 определен при всех 0 ^ то ^ 1 (рис. 1). На границе а = 2а ФП Т22 при tq = 0 равна единице, а при переходе через границу в сторону больших значений S становится определенной в интервале tq- ^ tq ^ 1, где г(го-) = 1.
При этом в интервале 0 ^ то ^ то_ определена ФП Т23, так что итоговый график ФП подпространства П2 в себя в этом интервале соответствует произведению Т32Т23. В силу этого для рассмотрения условий глобальной устойчивости неподвижной точки т* отображения Т22 в этом случае необходимо исследовать характер графика ФП отображения Т32Т23.
_ 1
а
_ 1
а
то-
то
_ А т Л к
/ —-V
/ /
/ ^ / 1
/ / Т33 —»-
/ 1 Л 4Т32
/ 17
А а 1 1 2-Шг - 1 а.
-\ —^01 1 -^^ о
\ Т23 1 4 \ 1 \ \ 1
\ \ - — \
\ \ \ \
Рис. 1. График ФП отображения П2 в себя при а, аътах ^ а ^ а, 2а
Рис. 2. Построение ФП произведения Т32Т23 при а/2 < а <а, бог > 2а/а
1
0
о
1
При одновременном выполнении неравенств а/2 < а <а, 9о1 > 2а/а способ качественного графического построения ФП произведения Т32Т23 приведен на рис. 2, а на рис. 3 — соответствующий вид графиков ФП отображения подпространства П2 в себя. На рис. 3 использовано обозначение т23 = Т—10о1- Количество периодов графика ФП в интервале 0 ^ то ^ то_ увеличивается по мере приближения графика ФП Т33 к биссектрисе и при а —> а — 0 неограниченно возрастает. При этом график ФП Т32Т23 поднимается вверх, уменьшаясь по амплитуде при т23 —> 1. При 001 = 2а/а Т23 = 0, и на рис. 3 пилообразность графика ФП Т32Т23 исчезает, так что остается только один ниспадающий участок.
При а —> а/2 + 0 001 0 график ФП стягивается в точку во = 0, в = 1, правый конец графика Т32 опускается до нуля, а единственный ниспадающий участок графика ФП Т32Т23 пересекает биссектрису г = то в нижнем левом углу квадранта то,т > 0, доказывая тем самым возможность существования неподвижной точки то = т ** цикл а Т32Т23, отличной от неподвижной точки то = т * отображен ия Т22-
То есть при указанном подходе представляется возможность убедиться, что отображение Т32 играет роль ловушки, из которой величина т ИТД поступает в интервал определения отображения Т22- Таким образом, можно качественно сформировать графики ФП в квадранте то,т > 0. При этом в интервале то- ^ то ^ 1 представлен график отображения Т22, а в интер-
Рис. 3. График ФП отображения П2 в себя при а/2 < а < а, в01 > 2а/а
вале 0 ^ то ^ то_ — сформированный график произведения T32T33T23, то есть предоставляется возможность качественного исследования процессов в СЧ в «большом». Отметим, что график произведения указанных отображений напоминает собой неравномерную пилу, располагающуюся по высоте в диапазоне (2а/а) — 1 ^ т ^ а/а. График ФП отображения Т22 располагается при т_о ^ то ^ 1 в верхней части квадрата.
Непосредственно из рисунка графиков ФП видно, что если точка (т_о, (2а/а) — 1) лежит выше биссектрисы координатного угла то, т > 0, то итерации с течением времени приводят в режим ИФД, обеспечивая глобальную устойчивость неподвижной точки отображения T22. Если указанная точка лежит ниже биссектрисы, что реализуется при больших S, то произведение отображений T32T23 имеет неподвижную точку. Последнее означает либо потерю глобальной устойчивости неподвижной точки то = т* режима ИФД, либо увеличение длительности переходных процессов вследствие прохождения ИТД в окрестности неустойчивой неподвижной точки седлового типа [5].
Список литературы
1. Левин В.А., Малиновский В.Н., Романов С.К. Синтезаторы частот с импульсно-фазовой автоподстройкой. М.: Радио и связь, 1989. 232 с.
2. Горюнов В.И. К теории систем импульсно-фазовой автоподстройки частоты (ИФАПЧ) // Изв. вузов. Сер. Приборостроение. 1974. № 10. С. 40-43.
3. Антоновская О.Г., Горюнов В.И. Метод точечных отображений и исследование динамики синтезаторов частот. Гамбург: Ламберт, 2014. 101 с.
4. Антоновская О.Г., Горюнов В.И. Об особенностях методики исследования динамики системы с широтно-импульсной модуляцией и запоминание сигнала управления / / Вестник Нижегородского университета. 2008. № 6. С. 135-140.
5. Горюнов В.И., Ерусланов В.Н., Лобашов II.II. Техническая полоса захвата одноконтурного синтезатора частоты // Техника средств связи. Сер. Техника радиосвязи. Вып. 2. 1990. С. 88-94.
Поступила в редакцию 24.08.2015
Антоновская Ольга Георгиевна, к. ф.-м. н., доцент, кафедра математики, Нижегородский государственный архитектурно-строительный университет, 603950, Россия, г. Нижний Новгород, ул. Ильинская, 65. E-mail: olga.antonovsckaja@yandex.ru
Горюнов Владимир Иванович, к. ф.-м. н., ведущий научный сотрудник, Нижегородский государственный университет, Институт информационных технологий, математики и механики, 603950, Россия, г. Нижний Новгород, проспект Гагарина, 23. E-mail: pmk@unn.ac.ru
O. G. Antonovskaya, V. I. Goryunov
Modeling the process of gluing phase trajectories in systems with a combined frequency-phase control
Keywords: control systems, mathematical modeling, point mappings method. MSC: 37E05, 37G15, 93C55, 93D20
The paper presents the use of the method of investigating the dynamics of systems with width-pulse modulation of control signal taking into account its piecewise constancy to studying the dynamics of the synchronization system with combined width-pulse modulation of control signal. The method is based on applying the point mappings method and consists of using the phase subspaces corresponding to constant control signal. In particular, the procedure of modeling the transition process from phase control regime to frequency control regime in the system with combined frequency-phase control is given.
REFERENCES
1. Levin V.A., Malinovskii V.V., Romanov S.K. Sintezatory chastot s impul'sno-fazovoi avtopodstroikoi (Frequency synthesizers with pulse-phase auto-tuning), Moscow: Radio i svyaz', 1989, 232 p.
2. Goryunov V.I. On the theory of systems with pulse-phase frequency auto-tuning, Izv. Vyssh. Uchebn. Zaved., Priborostroenie, 1974, no. 10, pp. 40-43 (in Russian).
3. Antonovskaya O.G., Goryunov V.I. Metod tochechnykh otobrazhenii i issledovanie dinamiki sintezatorov chastot (Point mappings method and frequency synthesizers dynamics study), Hamburg: Lambert Academic Publishing, 2014, 101 p.
4. Antonovskaya O.G., Goryunov V.I. On some peculiarities of the method for studying the dynamic of systems with pulse-width modulation and control signal storage, Vestn. Nizhegorod. Univ., 2008, no. 6, pp. 135-140 (in Russian).
5. Goryunov V.I., Eruslanov V.N., Lobashov N.I. Technical capture wave range of single-circuit frequency synthesizer, Tekhnika sredstv svyazy. Ser. Tekhnika radiosvyazy, 1990, no. 2, pp. 88-94 (in Russian).
Received 24.08.2015
Antonovskaya Olga Georgievna, Candidate of Physics and Mathematics, Associate Professor, Mathematics Department, Nizhni Novgorod State University of Architecture and Civil Engineering, ul. Il'inskaya, 65, Nizhni Novgorod, 603950, Russia. E-mail: olga.antonovsckaja@yandex.ru
Goryunov Vladimir Ivanovich, Candidate of Physics and Mathematics, Leading Researcher, Nizhni Novgorod State University, Institute on Information Technologies, Mathematics and Mechanics, pr. Gagarina, 23, Nizhni Novgorod, 603950, Russia. E-mail: pmk@unn.ac.ru
