CC BY
62
УДК 621.396.662
ПЕРЕХОДНЫЕ ПРОЦЕССЫ В СИНТЕЗАТОРАХ
С ФАЗОВОЙ АВТОПОДСТРОЙКОЙ ЧАСТОТЫ
ПРИ АДАПТИВНОЙ КОМПЕНСАЦИИ ПОМЕХ ДРОБНОСТИ
С.К. Романов1, Н.М. Тихомиров1, А.В. Леньшин2
1ОАО "Концерн "Созвездие", г. Воронеж
e-mail: skromanov@rambler.ru; tikhomir@sozvezdie.su;
2Военный авиационный инженерный университет, г. Воронеж e-mail: andrey-lenshin@yandex.ru
Предложена методика определения длительности переходного процесса при адаптивной компенсации помех дробности, генерируемых дельта-сигма-модуляторами в дробных делителях частоты с переменным коэффициентом деления систем импульсной фазовой автоподстройки частоты. Адаптивная компенсация помех дробности осуществляется с помощью специального дополнительного устройства, подключаемого на вход фильтра нижних частот.
Ключевые слова: адаптация, синтезатор частот, система импульсной ФАПЧ, дробное деление частоты, импульсный частотно-фазовый детектор, управляемый генератор, дельта-сигма-модулятор, адаптивная компенсация помех дробности.
TRANSIENTS IN PHASE-LOCKED LOOP SYNTHESIZERS WITH ADAPTIVE COMPENSATION OF FRACTIONAL SPURS
S.K. Romanov1, N.M. Tikhomirov1, A.V. Len'shin2
1OAO "Concern "Sozvezdie", Voronezh, Russia e-mail: skromanov@rambler.ru; tikhomir@sozvezdie.su
2Military Aviation Engineering University, Voronezh, Russia e-mail: andrey-lenshin@yandex.ru
A technique is proposed for determination of transient's duration in adaptive compensation of fractional spurs generated by delta-sigma modulators in fractional frequency dividers of phase-locked loop with a variable division coefficient. The adaptive compensation of fractional spurs is implemented using a dedicated additional device connected to the low-pass filter input.
Keywords: adaptation, frequency synthesizer, pulse PLL, fractional frequency division, pulse phase-frequency detector, controllable generator, delta-sigma modulator, adaptive compensation of fractional spurs.
В синтезаторах частот (СЧ), построенных на основе импульсно-фазовой автоподстройки частоты (ИФАПЧ), нашли широкое применение дробные делители частоты с переменным коэффициентом деления (ДДПКД). Преимущества подобных синтезаторов общеизвестны [1, 2]. Однако наличие ДДПКД в ИФАПЧ приводит к появлению помех дробности (ПД) в выходном сигнале СЧ. Для уменьшения уровня ПД в низкочастотной части спектра выходного сигнала СЧ в составе ДДПКД используется схема дельта-сигма-модулятора (ДСМ) в разных модификациях. Расчеты уровня ПД в ИФАПЧ с ДДПКД и ДСМ по линейной модели [3] показывают высокую эффективность применения
Рис. 1. Структурная схема системы ИФАПЧ с адаптивной компенсацией помех дробности
ДСМ. Однако в широкополосных системах ИФАПЧ в области высоких частот уровень ПД оказывается недопустимо большим значение. Для устранения этого нежелательного явления предлагается использовать специальные устройства компенсации: ненастраиваемые [4-10] и автоматически настраиваемые [11-14]. В настоящей работе рассматриваются автоматически настраиваемые устройства компенсации и основное внимание уделено определению времени переходного процесса компенсации как недостаточно представленному в цитируемой литературе. При моделировании адаптивной системы ИФАПЧ не учитываются и не рассматриваются ряд нелинейностей, присущих реальной адаптивной системе ИФАПЧ.
На рис. 1 приведена структурная схема исследуемой адаптивной системы ИФАПЧ. С выхода импульсного частотно-фазового детектора (ИЧФД) (на вход которого приходит опорный сигнал е0 (£) и сигнал обратной связи ес (£), где £ — текущее время), на вход фильтра нижних частот (ФНЧ) с передаточной функцией О (в) подается через сумматор СУМ1 сигнал гд (£). Сигнал гд (£) имеет импульсный характер и его вид зависит от типа используемого ИЧФД (рассматриваем трехуровневый ИЧФД с зарядовой накачкой) и очередности поступления фронтов импульсных сигналов е0 (£) в моменты времени пТ0 и ес (£) в моменты времени £п [15], где Т0 — период е0 (£), п — номер импульса е0 (£).
Напряжение еф (£) с выхода ФНЧ подается на устройство, которое моделирует управляемый генератор (УГ) и ДДПКД с коэффициентом деления частоты N. Это устройство образовано амплитудно-импульсной модуляцией на интервале времени от пТ0 до (п + 1)Т0, т.е.
(п+1)То (п+1)То
¿ТУА В этом выра-
при достижении условия
f гд(t)dt « [
uTn
uTn
жении ¿д(4) = гмsign(^дcм [пТо])[1(4 - пТо) - 1(4 - пТо - тдсм [пТо])], ¿м — амплитуда импульсов тока накачки ИЧФД;
^ДСМ(£) = ^дсм[пТо][1(4 - пТо) - 1(4 - (п + 1)То)]
— сигнал на выходе МУ;
2п п а
^дсм [пТо] = ^ ^ - ;
о к=1
Т
тдсм [пТо] = abs(^дcм [пТо]) —;
2п
¿туа(4) = ¿ТУА[п][1(4 - пТо) - 1(4 - пТо - ттуап)];
¿ТУА [п] — амплитуда импульсов тока ТУА; тТУАп — длительность импульсов тока ТУА на интервале времени от пТо до (п + 1)То.
На рис.2 условно показаны сигналы ¿д(4), ¿ТУА(4) и ^дСМ(4) в адаптивной системе ИФАПЧ, по оси X отложено время 4 в секундах, по оси У — токи ¿д (4), ¿ТУА (4) в амперах и ^дСМ (4) в радианах. С учетом работы [4], схемы на рис. 1 и графиков на рис. 2 условие компенсации ПД можно записать в виде
¿ТУА [п] = ^дсм [пТо ]. (1)
2п ТТУАп
Обычно выбирают ттуап = (4 ... 8)(4п+1 -)/^о, где (¿п+1 -¿п)/^о
— период колебаний сигнала УГ. Импульсы тока ¿ТУА (4) по длительности формируют из сигнала УГ, а по уровню — из сигнала ТУА с помощью выходного напряжения интегратора:
Рис.2. Эпюры сигналов 1Д(Ь) ¿тоа^дсмв системе ИФАПЧ с адаптивной компенсацией помех дробности
¿ТУА= ¿ТУА[п][1(£ - пТо) - 1(£ - пТо - Гтуап)] =
= Цм(^ТУА^ДСМ[пТо][1(* - пТо) - 1(* - пТо - ТТУАП)],
где Uint (£) — напряжение на выходе интегратора И, а $ТУА — крутизна характеристики управления аттенюатора ТУА.
Постановка задачи. Необходимо отметить, что точное выполнение соотношения (1) возможно только при использовании автоматической (адаптивной) системы настройки на минимум уровня помех дробности на выходе ФНЧ. Поскольку в работах [11-14] не даны рекомендации по нахождению параметра интегратора и зависимостей времени адаптации от и характеристик собственно системы ИФАПЧ, поэтому задачей настоящей работы является устранение пробелов в данном вопросе. Для решения поставленной задачи введем допущения:
Допущение 1 — частота опорного сигнала 1/То выше частоты среза ИФАПЧ, т.е. рассматриваем систему ИФАПЧ как непрерывную и описываемую с помощью системы обыкновенных дифференциальных уравнений в пространстве состояний.
Допущение 2 — тдСМ [пТо] ^ То — это соотношение обеспечивается в реальном устройстве при исследовании воздействий помех ДСМ и позволяет перейти от нелинейной системы ИФАПЧ к линейной системе ИФАПЧ.
Допущение 3 — тдСМ [пТо] ^ тТУАп — эта ситуация также, как правило, обеспечивается в реальных адаптивных системах ИФАПЧ.
С учетом этих допущений и [15] перейдем к анализу импульсно-непрерывной адаптивной системы ИФАПЧ, представленной на рис. 3, на котором в соответствии со схемой рис. 1 как аналоги блоков ТУА и интегратора И введены устройство перемножения П1 и аналоговый интегратор с передаточной функцией Кп/в, а также формирующие устройства ФУ1 и ФУ2. На выходах ФУ1 и ФУ2 сформированы соот-
Рис. 3. Линейная схема системы ИФАПЧ с адаптацией
Рис.4. Фильтр нижних частот цепи управления системы ИФАПЧ
ветствующие сигналы:
(t) = 2п • sign(^cm[nTo])[1(t - nTo) - 1(t - nTo - тдсм[nTo])]; To
^ (t) =
тТУАп
-^ДСМ[nTo][1(t - nTo) - 1(t - nTo - ТТУАП)].
Из сравнения двух адаптивных схем ИФАПЧ (на рис. 3 и 1) следует,
ЪТ ПС 2птТУАп
что Л^ = Gint¿туа . т—.
¿м То
При анализе адаптивной системы ИФАПЧ наряду с общими соотношениями (в качестве примера) будем ориентироваться на использование в цепи управления УГ простейшей схемы ФНЧ — ФНЧ4, представленной на рис.4. В отличие от подходов к исследованию в цитируемой литературе в настоящей работе сигнал управления вфа(4), снимаемый с ФНЧ4 для управления процессом адаптации, определяется как разность напряжений на конденсаторах С1 и С2, т.е. вфа(4) = ис1 (4) - ис2(4). Предлагаемый подход позволяет избежать проблем, связанных с нежелательным влиянием постоянной составляющей выходного сигнала ФНЧ4 на процесс адаптации.
Математическая модель адаптивной системы ИФАПЧ. Для получения математической модели адаптивной системы ИФАПЧ будем использовать дифференциальные уравнения, описывающие ее в пространстве состояний, вида
XX = АХ + Б{^1 (4) + ^ (4)};
(4) = Л* F (^дсм [пТо ])ефа(0, (2)
ефа (4) = СХ,
где X — вектор состояния размерностью к (напряжения на конденсаторах, токи в индуктивностях ФНЧ и фаза Фу (4)) УГ непосредственно системы ИФАПЧ; А — квадратная матрица состояния; Б — вектор управления; С — вектор-строка для вычисления выходной координаты вфа(4) ФНЧ4.
Для линейной схемы системы ИФАПЧ с адаптацией на рис. 3 имеем
-1/Т3 1/Тз -гм/(2пЯоС1) 1/Т1 -1/Т1 0
2П£уг 0 0
X = [Uci (t); UC2(t); Фу(t)] , A =
где Ti = Л1С2, T3 = RiCi, B =
iM/(2nCi) 0 0
C = [ 1 -1 0 ].
Для проведения дальнейших выкладок удобно перейти от системы уравнений (2) к уравнениям состояния адаптивной ИФАПЧ, расширенный вектор состояния Х^ которой имеет вид Х^ = [X; (£)]. Тогда систему уравнений (2) запишем в виде
ХЕ = ЛЕХЕ + А Б^2[пТо]
Кт^ (^дсм[ПТо])С 0
размером (к + 1) х (к + 1), БЕ = [Б; 0] — блочный вектор управления размером (к + 1), для ФНЧ4
где As[nTo] =
(3)
— блочная матрица
AsM=
-1/Тз 1/Ti 2nSyp
1/T3 -1/Ti
-iM/(2nNoCi) ¿M^2[t]/(2nCi)
KintF(^дсм [t]) -KintF(^дсм [t])
0 0 0
0 0
BE(t) =
iM(t)/(2nCi) 0 0 0
Уравнение (3) представляет собой неоднородное линейное дифференциальное уравнение с периодическими коэффициентами, так как ЛЕ[пТо] = ЛЕ[пТо + Тдсм], ЫпТо] = ^[пТо + Тдсм], где Тдсм — период ДСМ-последовательности и может в зависимости от типа ДСМ, порядка, емкости накопителей и числа а, подаваемого на вход первого накопителя ДСМ, составлять значение до 232 и более. Решение подобных уравнений не известно. Данную задачу можно упростить,
если интервал времени от пТо до (п + 1)То разбить на три подын-
Т
тервала п = тдсм[пТо] = аЪв(^дсм[пТо]) , Т2 = ттуап - тдсм [пТо]
2п
и т3 = То - тТУАп (см. рис. 2), на которых матрица АЕ постоянна. Обозначим матрицу на этих трех подынтервалах времени соответствующими индексами ЛЕ1, ЛЕ2, ЛЕ3, где
" Л Б ^2 [пТо] "
Kint Р (^дсм[ПТо])С 0
— блочная матрица размера (к + 1) х (к + 1), Л£2[пТо] = Л^пТо],
Л0
Кт^ (^дсм[ПТо])С 0
ASi[nTo] =
Аез [nTo ] =
0
0
— блочная матрица размера (к + 1) х (к + 1); для системы с ФНЧ4
A Si[nT о]
-I/T3 1/Ti 2nS
I/T3 -1/Ti
-iM/(2nNoCi) гм^2 [nTo]/(2nCi)'
УГ
KintF (^дсм [nT 0]) -Kint F (^дсм [nT 0])
0
0
Аеэ[пТ0] =
~ -I/T3 I/T3 -¿M/(2nNoCi) 0'
1/Ti -1/Ti 0 0
2nSyr 0 0 0
KintF(^дсм[ПТo]) -KintF(^дсм[ПТс]) 0 0
Поскольку матрица Asi[nTo] особенная (одно из собственных значений равно нулю), то с учетом введенных матриц и векторов решение (3) запишем в виде рекуррентного соотношения
Xs[(n +1)To] =
= e AS3 (T0 -тТУАп ) e AS2 [тТУАп тДСМ [nT0]] { gAS1 ТДСМ [nT0] Xs [nT0] +
Г Т1ГТ i Л ТДСМ[пТо1 .2 тДсм[пТо]
+ тдсм[пТо][1 + Asi-2--+ ASi-3--+ .. .]Bs^i[nTo]},
' (4)
где eAs"Ti — переходные матрицы на соответствующих подынтервалах времени t¿; Xs[nTo] — вектор состояния адаптивной ИФАПЧ в момент времени nTo; I — единичная матрица. Следует отметить, что в ряде разложения eAsiT«CM[nTo] при выполнении условия тдсм[пТ0] ^ To можно ограничиться несколькими первыми членами.
На рис. 5 приведен пример расчета по выражению (4) переходного процесса установления координаты AUint(t) = 1 — Uint(t) (кривая 1) от времени под действием ^дсм [nTo] (кривая 2), зависимость AUint(t) при представлении ^i(í) и (t) в виде £(t)-функций Дирака (кривая 3) (об этом подробно будет сказано далее). Начальное значение вектора X принималось нулевым, а AUint(0) = 0,1. ДСМ-последовательность формировалась как ДСМ третьего порядка типа MASH 1-1-1, значения емкостей накопителей принимались равными m = 2i6, число на входе первого накопителя принималось равным a = 1,тТУАп = 4ТУГ, период опорного сигнала выбирался равным To = 0,1 мкс, а F(ж) = ж, схема ФНЧ4.
Из анализа кривых на рис. 5 следует, что AUiní(t) в стационарном состоянии стремится к некоторому "малому" постоянному, отличному от нуля, значению с небольшими флуктуациями, вокруг которого наблюдается неидеальная компенсация помех ДСМ. Дополнительные ис-
■0,04 -1-1-1-1-1-1-1-
0 0,5 1,0 1,5 2,0 2,5 3,0 3,5 Х10"5
Рис. 5. Переходный процесс адаптации в системе ИФАПЧ
следования с использованием (4) показали, что с уменьшениемтТУАп [п] уменьшается и постоянная составляющаяДЦ"^(Ь) при Ь ^ то, а при замене ^(Ь) и (Ь) на 8(Ь)-функции ДЦ^ (Ь) ^ 0 при Ь ^ то. Необходимо отметить что выражение (4) громоздко и требует довольно больших вычислительных ресурсов для своего решения при больших
ТДСМ.
Таким образом, при полностью известных параметрах адаптивной ИФАПЧ выражение (4) позволяет определить время переходного процесса, при котором отклонения параметров состояния системы от стационарного значения достигают некоторого "малого" значения, задаваемого при проектировании системы.
Упрощенное решение задачи. Упростим решение задачи определения времени адаптации, представив идеализированный случай, когда на входы СУМ1 подаются сигналы в виде амплитудно-модулированных 8 (Ь)-функций вида
го
<Р1б (Ь) = То ^дсм [пТо] 8 (Ь - пТо)
п=0
и
го
¥25(Ь) = -То^ДСМ[пТо}иш(Ь) 8(Ь - пТо),
п=0
где То — период 8(Ь)-функций, а ¥дСМ[пТо] — амплитуда сигнала ДСМ.
Подставим в (2) <^1(5 (4) и (4) вместо и <£2(£) соответственно, а также введем расширенный вектор состояний Х^ = [X; ДЦ^(4)]. Тогда систему уравнений (2) запишем в виде
ХЕ5 = АЕ5 (*)ХЕ5, (5)
где
Asä (t) =
A -Ври (t)
KintF (рдсм (t))C 0
— блочная матрица размера (k +1) x (k + 1), для системы со схемой ФНЧ4
Asä (t) =
-1/Тз 1/Ti 2nSyr
1/Т3 -¿м/(2nNoCi) -¿M^ii (t)/(2nCi) -1/Ti 0 0
0 0 0
KintF(одсм(£)) -KintF(^дcм(í))
0
0
Уравнение (5) представляет собой однородное линейное дифференциальное уравнение с периодическими коэффициентами, так как As5 [nTo] = As5 [nTo + Тдсм], где ТДсм — период ДСМ-последователь-ности. Особенность дифференциального уравнения (5) является наличие ^(¿)-функций в составе матрицы A^(t). Используя методику решения дифференциальных уравнений, в состав которых входят обобщенные функции (к ним относятся £(t)-функции и ее производные [16]), решение (5) представим в виде рекуррентного соотношения
Xsä [(n + 1)To] = eA^nT°]T0 Aix [nTo ]XS5 [nTo],
(6)
где Aäx [nT0] =
является блочной единичной
1 -БТо^дсм [пТо] 0 1
диагональной матрицей размера (к + 1) х (к + 1), в которой (к + 1)-й столбец заменен на [-То^дСМ[пТо]Б; 1]; для схемы ФНЧ4
Aäx[nTo] =
1 0 0 - ¿M7o^cM[nTo]/(2nCi)
010 0 0 1 000
0 0 1
определяет скачок вектора состояния в момент времени пТо.
Выражение (6) значительно проще (4) и по нему определена кривая 3 переходного процесса, приведенная на рис. 5. Сравнение кривых 1 и 3 показывает, что даже при довольно малых отклонениях ДЦ^ (4) от положения равновесия они очень близки. Поэтому в первом приближении оценка времени адаптации может быть проведена с использованием выражения (6).
Определение времени адаптации. Перейдем к детальному определению времени адаптации с применением выражения (6). Если по
синтезу непосредственно ИФАПЧ есть большое число работ, то по определению влияния на процесс адаптации авторам такие работы неизвестны. Для исследования этого вопроса предлагается следующее:
— полагаем, что матрица А, вектор В, вектор-строка С известны, т.е. непосредственно система ИФАПЧ спроектирована с использованием какого-либо критерия качества;
— считаем, что на входы СУМ1 адаптивной системы ИФАПЧ подаются тестовые сигналы "меандр из -функций" вида
те те
= Горм{^ - 2пТо) - ^ - 2пТо - То)
n=0
n=1
и
те те
^(t) = -Т^м{^ 5(t - 2nTo) - ^ 5(t - 2nTo - To)}Umt(i) (7)
n=0
n=1
с периодом 2Т0 и амплитудой <^м, равной с отклонению ^дСМ[пТ0] на периоде ТдСМ =
реднеквадратическому
\
Едсм/ 0 ^Дсм[«то]
т^-^ ,т-е-
г=0
(ТДСМ/Т0 )
известен тип ДСМ, используемый в системе.
Используя (6), для сигнала (7) запишем решение на интервале времени от пТ0 до (п + 2)Т0:
ХЕ5 [(п + 2)Т0] =
= еАЫ(п+1)То]То АА 5ж[(п + ВД]еА -[пТо]То А 5я[пТ)]ХЕ5 [п^] =
= А 5жХЕ5 [пТ0], (8)
где _
А Еэ [пТ0 ] =
A 0
KmtF (^M)C 0
A 5ж[пТ0] = A Ea[(n +1)T)] =
1 -BT0^M 0 1
A
0
KmtF (-рм)С 0
А5ж [(n + 1)T0] =
1 BT0^M 01
А 5х = еА- К-+1)То]То А 5ж[(п + ВД]еА^То]То А 5ж[пТ0]
эквивалентная переходная матрица. Для системы ИФАПЧ со схемой
ФНЧ4 эти матрицы имеют вид
A £з[пТ о] =
-1/T3 1/T1 2nS
УГ
1/T3 -1/T1 0
-v/(2nN0C1) 0 00
KmtF (^м ) -KmtF (^м)
0 0
0 0
A Ез[(п +1)T 0] =
-1/T3 1/T1 2nS
УГ
1/T3 -1/T1 0
-¿м/(2п^0С1) 0
0 0
KmtF (-^м) -KmtF(-^M) 0
" 10 0 -¿мТ0^м/(2nCi) "
0 0 0
A äx[nT0 ] =
010 0 0 1 000
0 0 1
A sx[(n + 1)T0] =
1 0 0 гмТ0Рм/(2пС1)
0 1 0 0
0 0 1 0
0 0 0 1
Адаптивную систему ИФАПЧ с матрицей А¿ж будем считать эквивалентной адаптивной системой ИФАПЧ (ЭАИФАПЧ), являющейся стационарной линейной импульсной системой, переходный процесс в которой определяется в виде однородного разностного уравнения (8), его решение имеет вид
ХЕ(*) = (А5ж)2пХЕ(0), (9)
где £ = п2Т0, п = 1, 2, 3 ....
Решение (9) требует значительно меньше затрат машинного времени нежели (6), а также позволяет детерминировать задачу адаптации, заменив практически случайный процесс ^дСМ[пТ0] на детерминированный процесс (/31(^).
Результаты исследований. Далее приведем результаты исследований ЭАИФАПЧ четвертого порядка со схемой ФНЧ4. Для параметрического синтеза элементов непосредственно ИФАПЧ используются ряд критериев качества [14]. Применим достаточно популярный частотный критерий качества, гарантирующий устойчивость системы с определенным запасом. Это — задание запаса устойчивости по фазе ^зап ФЧХ на частоте среза шс амплитудно-частотной характеристики разомкнутой системы.
Выражения для параметрического синтеза элементов ФАПЧ в этом случае (задаемся значениями частоты ш0 и запаса по фазе ^зап) имеют вид [17]
^с = ^о2 Т, Т /Т2 = 1 + 2tg2 (^зап) + 1,
T =
1 (f
-о V T2
0,25
¿м S
где -о =
м°УГ , T2=RCiC2/(Ci + C2), Ci + C2 = ¿мSyr
(С + С2' -1 ' ш2
С = (С + С2) Т2/Г1, С2 = (С + С2) - С1, Я1 = Т1 /С2.
Введем также еще некоторые нормированные параметры: КщЬ связывающий К^ с другими параметрами ЭАИФАПЧ, Кп = = К1п1гмТ02^м/(2пС1), ¿¡0 = Т0 — нормированная базовая частота ИФАПЧ, ¡а = Ьа/Т0 — нормированное время адаптации (в качестве времени адаптации £а примем время, за которое координата (£) уменьшается в 100 раз при всех остальных нулевых начальных координатах).
Для примера на рис. 6 приведены результаты расчета переходного процесса адаптации в ЭАИФАПЧ (кривые 1, 2, 3) и ИФАПЧ (кривая 4) для случая К^ = 0,21, ¿¡0 = 0,0628. В этом примере система ИФАПЧ синтезирована на запас по фазе ^зап = 45°, Т0 = 0,1 мкс, АЦ^(0) = 0,1, Т(х) = х. На рис.6 кривая 1 — результаты расче-
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 Х10" Рис. 6. Переходный процесс адаптации в системах ЭАИФАПЧ и ИФАПЧ
тов AUint(t) по (9), кривая 2 — результаты расчетов AUint(t) по (6), кривая 3 — зависимость [nT0] (показана условно, как решетчатая функция) для = 0,02, кривая 4 — результаты расчетов AUint(t) по (4) для тестового сигнала удсм(t) = yMsign[sin(nt/T0)].
Из анализа кривых на рис. 6 следует, что время ta для этого случая составляет 4,8 мкс. Кривая 1 представляет собой сложного вида форму огибающей для кривой 2. Кривая 4, как и кривая 1 на рис. 5, стремится к некоторому постоянному значению, которое для предлагаемого тестового сигнала удсм (t) можно получить из (4) как
Xs [nTo ] = ASc eAE3 (To-Ттуап } eA^2 [ттуаП-ТДсм [nTo И {2птдсм [nTo ][I+
n—YOQ
+ Asi ТдсмМ0. + A2
тдсм
[nTo ]
2 ' 3!
где Л^с определяется инверсной матрицей от суммы матрицы I и матрицы еАЕ3 (То-ТтуАп)еАЕ2ТтуАп . В этой инверсной матрице Л^с знаки элементов последнего столбца изменены на противоположные. В кодах МЛТЬЛБ эти операции можно представить как Л^32 = = еАЕ3(То-ТтуАп)еАЕ2ТтуАп, ЛЕ32 = -ЛЕ32= 1), ЛЕС = шУ(1 + + ЛЕз2).
Результаты обработки расчетов времени адаптации в ЭАИФАПЧ с использованием (9), введенных нормированных параметров, для запаса по фазе ^зап = 45° и ^(х) = х приведены на рис. 7.
На этом рисунке по осям X и У отложены соответственно КК^ и ¡а, а в качестве дополнительного параметра используется значение частоты а0. Из анализа графиков на рис. 7 следует, что с увеличением значения КК^ и уменьшением а0 время адаптации уменьшается.
+ ...]BE },
Рис. 7. Зависимости длительностей процесса адаптации в системе ЭАИФАПЧ
На рис. 7 приведены результаты использования значений параметров Кы ^ 0,0345 и (0 ^ 0,125. Дополнительные исследования показали, что при уменьшении К^ относительно значения 0,0345 при частоте (0 ^ 0,125 в к раз время адаптации увеличивается примерно в к раз. Интервалы времени адаптации в системе при частотах (0 = 0,125 и (Ш0 = 0,25 мало отличаются друг от друга, т.е. изменения значений параметра а0 при выборе а0 ^ 0,125 слабо влияют на время переходного процесса.
Приведем на рис. 8 некоторые результаты по расчету спектров фазовых шумов с использованием (4) в адаптивной системе ИФАПЧ с параметрами: частота среза разомкнутой системы ИФАПЧ /с = 1,25 МГц, ^зап = 45°, Кы = 0,168, ДСМ типа МЛ8И111, т = 216, а =1, (¡¡0 = 0,75, ттуап = 4Туг.
Время адаптации в системе с такими параметрами при начальном векторе Х^[0] = [0; 0; 0; 1,1] составляет примерно 35мкс. После завершения переходного процесса с помощью //¿-преобразования в системе МЛТЬЛБ вычислялись спектры фазовых помех сигнала УГ на частоте 2,210000152587891 ГГц для двух случаев: отсутствие компенсации (при АЦп1 = -1) — кривая 1 и адаптивная компенсация по модели системы по схеме рис. 3 — кривая 2.
На рис. 8 по оси X отложено значение частот помех дробности, по оси У отложено значение уровня фазовых помех дробности в децибелах. Сравнение кривой 1 и 2 показывает, что на частотах в районе частоты среза системы ИФАПЧ уровень компенсации достигает зна-
Рис. 8. Спектры фазовых помех сигнала УГ в адаптивной системе ИФАПЧ с ДСМ
чения 49,7 дБ. К этому показателю уровня компенсации необходимо отнестись критически и считать его идеальным, потому что модель адаптивной системы ИФАПЧ (см. рис. 3) не учитывает ряд нелиней-ностей, присущих реальной адаптивной системе ИФАПЧ и не рассматриваемых в настоящей работе.
Выводы. Предложена методика определения длительности переходного процесса адаптации при компенсации помех дробности. Методика предполагает разбиение на три этапа:
первый — синтез непосредственно системы ИФАПЧ с использованием каких-либо критериев качества;
второй — с использованием введенных нормированных параметров адаптивной системы ИФАПЧ нахождение времени адаптации в эквивалентной адаптивной системе ИФАПЧ, когда сигнал с ДСМ заменяется "эквивалентным" сигналом и при этом определяется значение коэффициента усиления интегратора Kint;
третий — уточнение времени адаптации при использовании реальных сигналов с выхода ДСМ и устройства компенсации по предложенным формулам.
СПИСОК ЛИТЕРАТУРЫ
1. Варфоломеев Г. Ф. Спектр помех дробности в системе фазовой АПЧ с дробным делителем частоты // Техника средств связи. Сер. ТРС. - 1978. -Вып. 10 (21). - С. 66-71.
2. К о з л о в В. И., П а л е н к о в А. В., Р я п о л о в А. А. Синтезатор частот с модуляцией дробных коэффициентов деления в петле ФАПЧ // Электросвязь. -1988. -№№ 9. - С. 48-50.
3. Р о м а н о в С. К., Марков И. А. Определение помех дробности в синтезаторах частот с системами ФАПЧ, использующих дельта-сигма модуляторы в дробных делителях частоты // Теория и техника радиосвязи: Науч.-техн. сб. / ОАО "Концерн "Созвездие". - Воронеж, 2006. - Вып. 1. - С. 97-102.
4. Романов С. К., Марков И. А. Тихомиров Н. М. Пути уменьшения помех дробности в синтезаторах с системами ИФАПЧ, использующих дельта-сигма модуляторы в дробных делителях частоты // Теория и техника радиосвязи: Науч.-техн. сб. / ОАО "Концерн "Созвездие". - Воронеж, 2007. - Вып. 1. - С. 7077.
5. Meninger S. E., and Perrott M. H. A fractional-N frequency synthesizer architecture utilizing a mismatch compensated PFD/DAC structure for reduced quatization-induced phase noise, IEEE Trans. Circuits Syst // Analog Digit. Signal Process. - Vol. 50, no. 11. -P. 839-849, Nov. 2003.
6. M e n i n g e r S. E., and Perrott M. H. A 1 MHz bandwidth 3,6 GHz 0,18 um CMOS fractional-N synthesizer utilizing a hybrid PFD/DAC structure for reduced broadband phase noise // IEEE Journal of Solid-State Circuits. - Vol. 41, no. 4. -P. 966-980, April 2006.
7. P a m a r t i S., and G a 11 o n I. Phase-noise cancellation design tradeoffs in delta-sigma fractional-N PLLs // IEEE Trans. Circuits Syst // Analog Digit. Signal Process. - Vol. 50, no. 11. - P. 829-838, Nov. 2003.
8. PamartiS., Jansson L. and G a 11 o n I. A wideband 2,4 GHz delta-sigma fractional-N PLL with 1 Mb/s in-loop modulation // IEEE J. Solid-State Circuits. -Jan. 2004. - Vol. 39, no. 1. - P. 49-62.
9. Wang K., Swaminathan A., Galton I. Spurious-tone suppression techniques applied to a wide-bandwidth 2,4 GHz fractional-N PLL // IEEE International Solid-State Circuits Conference. Digest of Technical Papers. February 2008.
10. W a n g K., S w a m i n a t h a n A. and Galton I. Spurious tone suppression techniques applied to a wide-bandwidth 2,4 GHz fractional-N PLL solid-state circuits // IEEE Journal. - 2008. - Vol. 43. - P. 2787-2797.
11. Gupta M. and Song B. -S. A 1,8-GHz spur-cancelled fractional-frequency synthesizer with LMS-based DAC gain calibration // IEEE J. Solid-State Circuits. -2006. - Vol. 41, no. 12. - P. 2842-2851.
12. Swaminathan A., Wang K., Galton I. A wide-bandwidth 2,4 GHz ISMband fractional-N PLL with adaptive phase-noise cancellation // In IEEE ISSCC Dig. Tech. Papers. - 2007. - P. 302-303.
13. Swaminathan A., Wang K., Galton I. A wide-bandwidth 2,4 GHz ISM-band fractional-N PLL with adaptive phase-noise cancellation // IEEE Journal of Solid-State Circuits. - 2007. - Vol. 42, no. 12. - P. 2639-2650.
14. H s u C. -M., Straayer M. ., Perrot M. H. A Low-noise wide-BW 3,6GHz digital fractional-N frequency synthesizer with a noise-shaping time-to-digital converter and quantization noise cancellation // IEEE Journal of Solid-State Circuits. - 2008. - Vol. 43, no. 12. - P. 2776-2786.
15. Левин В. А., Малиновский В. Н., Романов С. К. Синтезаторы частот с системой импульсно-фазовой автоподстройки. - М.: Радио и связь, 1989.-232 с.
16. РозенфельдА. С. Скачки и дельта-импульсы вектора состояния линейных электрических цепей // Изв. вузов СССР. Радиоэлектроника. - 1978. - Т. XXI. -№ 11. - С. 67-73.
17. Романов С. К., Тихомиров Н. М., Рахманин Д. Н. Методика определения быстродействия синтезаторов частот с коммутацией токов накачки и постоянных времени фильтра нижних частот // Вестник МГТУ им. Н.Э. Баумана, Сер. Приборостроение. - 2010. - № 3 (80). - С. 79-93.
Статья поступила в редакцию 30.01.2012
Станислав Константинович Романов — заместитель начальника научно-технического центра в ОАО "Концерн "Созвездие". Автор более 110 научных работ в области синтеза частот.
S.K. Romanov — deputy head of scientific and technical center of OAO "Concern "Sozvezdie". Author of more than 110 publications in the field of frequency synthesis.
Николай Михайлович Тихомиров — д-р техн. наук, начальник научно-технического центра в ОАО "Концерн "Созвездие". Автор более 100 научных работ в области синтеза частот.
N.M. Tikhomirov — D. Sc. (Eng.), head of scientific and technical center of OAO "Concern "Sozvezdie". Author of more than 100 publications in the field of frequency synthesis.
Андрей Валентинович Леньшин — д-р техн. наук, профессор Военного авиационного инженерного университета (г.Воронеж). Автор более 130 научных работ в области синтеза частот и радиопередающих устройств.
A.V. Lenshin — D. Sc. (Eng.), professor of Military Aviation Engineering University (Voronezh). Author of more than 130 publications in the field of frequency synthesis and radiotransmitting devices.
