MATHEMATICAL MODEL OF OSCILLATION FORMATION WITH THE USE OF PASSIVE DIGITAL SYNTHESIS
52 Ю. А. Никитин
УДК 621.391
Ю. А. НИКИТИН
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ФОРМИРОВАНИЯ КОЛЕБАНИЙ С ИСПОЛЬЗОВАНИЕМ МЕТОДОВ ПАССИВНОГО ЦИФРОВОГО СИНТЕЗА
Рассмотрена математическая модель, предназначенная для анализа и расчета спектров двух- и многоуровневых колебаний с равномерной и неравномерной дискретизацией по фазе, а также многоуровневых колебаний, характеризуемых конечной точностью выбора узлов аппроксимации. Получены аналитические выражения для спектров колебаний с различными огибающими.
Ключевые слова: пассивный цифровой синтез, накапливающий сумматор, счетчик импульсов, делитель с дробно-переменным коэффициентом деления.
Методы двух- и многоуровневого пассивного цифрового синтеза (ПЦС) частот обеспечивают выполнение целого ряда жестких и взаимно противоречивых требований, предъявляемых к современной аппаратуре синтеза частот и обработки сигналов. К таким требованиям относят высокую долговременную стабильность фазы выходного колебания и ее преемственность при смене частот, малое время переключения с частоты на частоту, низкий относительный уровень дискретных побочных спектральных составляющих, малый уровень фазовых шумов.
Однако на спектр формируемого колебания существенным образом влияет алгоритм работы системы ПЦС как (цифрового) конечного автомата (КА). Наилучший спектр обеспечивается посредством методов оптимального синтеза, использование которых в случае двухуровневого ПЦС [1] позволяет аппаратно реализовать конечный автомат в виде накапливающего сумматора (НС) и/или делителя с дробно-переменным коэффициентом деления (ДДПКД), а в случае многоуровневого ПЦС — только НС, дополненный цифроаналоговым преобразователем или управляемым устройством задержки [2].
Целью настоящей статьи является описание работы конечного автомата применительно к теории синтеза частот и рассмотрение математической модели формирования периодического колебания с огибающей произвольной формы.
Рассмотрим структурные схемы некоторых простейших КА, применяемых в технике синтеза частот, и описывающие их работу алгоритмы [1—5].
Наиболее широко в системах активного и пассивного цифрового синтеза частот используются счетчики импульсов (СИ) (могут содержать поглотители импульсов) и накапливающие сумматоры (рис. 1, а, б соответственно).
а) б)
Fвых 1
Sk K fвх
Fвых
Q
Sk Р fвх
1 СИ
Q НС
Рис. 1
Функция выхода (сигнал переполнения) ρk СИ определяется как
ρk = ⎣k/P⎦ – ⎣(k–1)/P⎦; ρk ∈ (0, 1),
(1)
где ⎣*⎦ — оператор выделения целой части числа, меньшей или равной ему; при этом дли-
тельность такта на входе КА равна Т0.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Математическая модель формирования колебаний с использованием методов ПЦС
53
Функция переходов (текущая сумма) Sk для СИ определяется следующим образом:
Sk = P { k/P },
(2)
где {∗} — оператор выделения дробной части числа, 0 ≤ { k/P } < 1 [5]. Изменение коэффициента пересчета СИ с дробно-переменным коэффициентом деления (N)
осуществляется с помощью поглотителя импульсов и накапливающего сумматора (рис. 2). Такое решение КА широко применяют в системах активного цифрового синтеза частот на основе колец импульсно-фазовой автоподстройки частоты. В последние годы в зарубежной литературе накапливающие сумматоры также называют аккумуляторами фазы и ∆Σ-модуляторами: первое название пришло из теории и техники пассивного (прямого) цифрового синтеза (Direct Digital Synthesis — DDS), а второе — из теории и техники аналого-цифрового преобразования.
Из двух последовательно включенных НС первого порядка (рис. 3) состоит накапливающий сумматор второго порядка — НС2 (и т.д.) (рис. 4).
⎣(kt)/T0⎦
ДДПКД
P 1/⎣Ν⎦
NT0
Q ΣP
HC1 Sk
ρk {kQ/P}
D Такты
Q Sk
ρk =0; 1
Sk+1
Рис. 2
Рис. 3
Q ΣP
SL D
HC2
ΣP Sk
D
ρk =–1; 0; +1; +2
Y1 Σ
+ +
Сигналы управления –
на СИ (ДДПКД) Такты от ДДПКД Y2D
SL+1 Y2
D
D Sk+1
Рис. 4
Если принять емкость НС первого порядка (n =1) равной Р единиц (см. рис. 3), то на его входе можно записать число Q ∈ (0, 1, 2, … Р–1).
Функцию ρk для НС1 можно представить как
ρk = ⎣kQ/P⎦ – ⎣(k–1) Q/P⎦; ρk ∈ (0, 1),
(3)
а функцию Sk — как
Sk = P{kQ/P}.
(4)
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
54 Ю. А. Никитин
Сравнивая выражения (1) и (3), можно сделать вывод, что СИ является частным случаем
НС1 при Q = 1. В табл.1 приведены примеры записи функций переходов (текущей суммы Sk) и функции
выхода (сигнала переполнения ρk) для вариантов НС1 различной емкости Р.
Таблица 1
Функция 0
1
Номер такта 23456
7
8
Р = 7; Q = 2
Sk 0 2 4 6 1 3 5 0 2
ρk 1 0 0 0 1 0 0 1 0 Р=7; Q =3
Sk 0 3 6 2 5 1 4 0 3 ρk 1 0 0 1 0 1 0 1 0
Р=8; Q =3
Sk 0 3 6 1 4 7 2 5 0 ρk 1 0 0 1 0 0 1 0 1
В табл. 2—5 приведены примеры записи функций переходов и выходов для вариантов
НС2 различной емкости Р; в таблицах приняты следующие обозначения: SL — текущая сумма 1; SL+1 — текущая сумма 1 при задержке (D) на один такт; Sk — текущая сумма 2; Sk+1 — текущая сумма 2 при задержке на один такт; Y1, Y2 — сигналы управления; Y2D — сигнал управления Y2 при задержке на один такт; ρk — сигнал управления СИ.
Р = 5; Q = 2
Таблица 2
Функция 0
1
2
3
4
Номер такта 5678
9
10 11 12
SL 2 4 1 3 0 2 4 1 3 0 2 4 1 SL+1 0 2 4 1 3 0 2 4 1 3 0 2 4 Sk 2 1 2 0 0 2 1 2 0 0 2 1 2 Sk+1 0 2 1 2 0 0 2 1 2 0 0 2 1 Y1 1 0 0 1 0 1 0 0 1 0 1 0 0 Y2 0 1 0 1 0 0 1 0 1 0 0 1 0 Y2D 0 0 –1 0 –1 0 0 –1 0 –1 0 0 –1 ρk 1 1 –1 2 –1 1 1 –1 2 –1 1 1 –1
Р = 7; Q = 3
Таблица 3
Функция 0 1 2
Номер такта 3456 7
8
9
10
SL 3 6 2 5 1 4 0 3 6 2 5 SL+1 0 3 6 2 5 1 4 0 3 6 2 Sk 3 2 4 2 3 0 0 3 2 4 2 Sk+1 0 3 2 4 2 3 0 0 3 2 4 Y1 1 0 0 1 0 1 0 1 0 0 1 Y2 0 1 0 1 0 1 0 0 1 0 1 Y2D 0 0 –1 0 –1 0 –1 0 0 –1 0 ρk 1 1 –1 2 –1 2 –1 1 1 –1 2
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Математическая модель формирования колебаний с использованием методов ПЦС
Р = 8; Q = 3
Таблица 4
Функция 0
1
2
3
4
5
6
Номер такта 7 8 9 10 11 12 13 14 15 16 17
SL 3 6 1 4 7 2 5 0 3 6 1 4 7 2 5 0 3 6
SL+1 0 3 6 1 4 7 2 5 0 3 6 1 4 7 2 5 0 3
Sk 3 1 2 6 5 7 4 4 7 5 6 2 1 3 0 0 3 1
Sk+1 0 3 1 2 6 5 7 4 4 7 5 6 2 1 3 0 0 3
Y1 1 0 0 1 0 0 1 0 1 0 0 1 0 0 1 0 1 0
Y2 0 1 0 0 1 0 1 0 0 1 0 1 1 0 1 0 0 1
Y2D 0 0 –1 0 0 –1 0 1 0 0 –1 0 –1 –1 0 –1 0 0
ρk 1 1 –1 1 1 –1 2 –1 1 1 –1 2 0 –1 2 –1 1 1
55
Р = 10; Q = 3
Таблица 5
Функция 0
1
2
3
4
5
6
7
8
Номер такта 9 10 11 12 13 14 15 16 17 18 19 20
SL 3 6 9 2 5 8 1 4 7 0 3 6 9 2 5 8 1 4 7 0 3
SL+1 0 3 6 9 2 5 8 1 4 7 0 3 6 9 2 5 8 1 4 7 0 Sk 3 9 8 0 5 3 4 8 5 5 8 4 3 5 0 8 9 3 0 0 3
Sk+1 0 3 9 8 0 5 3 4 8 5 5 8 4 3 5 0 8 9 3 0 0
Y1 1 0 0 0 1 0 0 1 0 0 1 0 0 0 1 0 0 1 0 0 1
Y2 0 0 1 1 0 1 0 0 1 0 0 1 1 0 1 0 0 1 1 0 0
Y2D 0 0 0 –1 –1 0 –1 0 0 –1 0 0 –1 –1 0 –1 0 0 –1 –1 0
ρk 1 0 1 0 0 1 –1 1 1 –1 1 1 0 –1 2 –1 0 2 0 –1 1
В табл. 1—5 фоном выделен период повторения цикла счета.
Для анализа спектра выходных колебаний КА, дискретных во времени с шагом T0 = 1/f0
и/или квантованных по амплитуде с шагом s, в работе [6] рассмотрена математическая мо-
дель (рис. 5), пригодная для формирования цифрового колебания (Ав(t)) с периодической оги-
бающей любой формы. Модель позволяет учесть конечную точность воспроизведения оги-
бающей при использовании ограниченного набора строго периодических колебаний разной
амплитуды и частоты следования, причем на периоде неравномерности структуры выходного
потока импульсов Tс = 1/Fс укладывается целое число периодов вспомогательных колебаний
модели; здесь Tс =QTв = РT0, Tв = 1/Fвых. Аналогичная структура, пригодная для компьютерно-
го моделирования и анализа спектров цифровых двухуровневых колебаний, рассмотрена в
работе [7].
A0(t)
Идеальный
Aв(t)
формирующий
A2(t) фильтр
A1(t)
A=1
ИФФ
A3(t)
T0
Рис. 5
Например, для модели формирования цифровой моногармоники (ЦМ) с равномерной
дискретизацией по времени T0 целесообразно выбрать следующие исходные строго периодические колебания:
A0(t) — периодическую последовательность δ-импульсов с периодом следования T0; A1(t) = sin(2πQt) — моногармоническую огибающую требуемой (синтезируемой) частоты;
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
56 Ю. А. Никитин
A2(t) = A3(t) = 0; ИФФ — идеальный формирующий фильтр с шириной окна, равной T0. Для спектра ЦМ можно записать
AЦМ (n)
=
(−1)n P
sin(πQ / P) π(nP ± Q)
,
где n = 0,1, ...
Аналогично, если принять
A
ЦТ
(t)
=
⎧1 + ⎨⎩1 −
X X
/π /π
для для
− π ≤ X ≤ 0; 0 ≤ X ≤ π;
A2(t) = A3(t) = 0, X = t/T0, PT0/Q = π, 0 ≤ t ≤ T0, получим модель формирования цифрового треугольного колебания (треугольника — ЦТ) с
равномерной дискретизацией по времени T0, для спектра которого можно записать
AЦТ (k)
=
2λµ
sin[π(Q ± k) / πP(Q ± k) sin2[π(R
P] + 1)
/
, P]
где µ = 1 ∀ P ≡ 0 (mod 2); µ = cos[π(R +1) / P] ∀ P ≡ 1 (mod 2) ; λ = 1cos[π(R +1)] ; k = 0, 1, 2, ..., Q–1;
R = (–1)r–1 (±k)Pr–1; здесь Pr–1 — числитель предпоследнего (r–1)-го члена разложения коэффициента деления N = 2P/Q в цепную дробь по алгоритму Евклида [5].
Одна из важных характеристик выходных колебаний КА, применительно к задачам синтеза частот, — уровень E полезной составляющей спектра, отнесенный к амплитуде A моногармонического колебания требуемой частоты: G = 20lg(E/A). Параметр G позволяет оценить степень приближения формируемого конечным автоматом цифрового колебания к своему аналоговому прототипу. На рис. 6 показан график изменения параметра G в зависимости от коэффициента N = P/Q деления (частоты) КА для разных огибающих цифрового колебания: идеальной моногармоники с равномерной дискретизацией по времени, псевдо- и квазимеандра, цифрового треугольника с равномерной дискретизацией по времени [6, 7]. Здесь величины P и Q выражают соответственно номинальные значения частот опорного f0 и синтезируемого (выходного) Fвых колебаний в единицах общей меры Fс = f0/P = Fвых/Q — шага сетки синтезируемых частот.
G, дБ 0
Идеальная моногармоника
–3,9 Идеальный меандр
–7,8 Идеальный треугольник
Цифровой треугольник
Псевдо- и квазимеандр
Цифровая моногармоника
Рис. 6
1/N = Q/P
Предложенная математическая модель позволяет производить численный расчет спектров формируемых конечным автоматом двух- и многоуровневых колебаний с равномерной и неравномерной дискретизацией по фазе, а также многоуровневых колебаний, характеризуемых конечной точностью выбора узлов аппроксимации, что, в свою очередь, позволяет получать аналитические выражения в свернутом виде для спектров колебаний с различной огибающей.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Математическая модель формирования колебаний с использованием методов ПЦС
57
СПИСОК ЛИТЕРАТУРЫ
1. Никитин Ю. А. Широкополосный синтез частот с помощью цифровых структур // Изв. вузов СССР. Приборостроение. 1990. Т. 33, № 9. С. 39—47.
2. Гуревич И. Н., Никитин Ю. А. Синтез сетки частот с произвольным шагом // Радиотехника. 1992. № 4. С. 53—58.
3. Брауэр В. Введение в теорию конечных автоматов: Пер. с нем. М.: Радио и связь, 1987. 392 с.
4. Апериодические автоматы / Под ред. В. И. Варшавского. М.: Наука, 1976. 424 с.
5. Виноградов И. М. Основы теории чисел. М.: Наука, 1972.
6. Никитин Ю. А. Математическая модель формирования колебаний с периодической огибающей // Тр. НИИР: Сб. статей. М., 1999. 124 с.
7. Никитин Ю. А. Конечный автомат как элемент цифровой системы синтеза частот // Докл. X Междунар. науч.-техн. конф. „Радиолокация, навигация, связь“. Воронеж, НПФ „САКВОЕЕ“ ООО, 2004. Т. 1. С. 526—533.
Юрий Александрович Никитин
Сведения об авторе — канд. техн. наук; Филиал ФГУП НИИ радио — Ленинградский отрас-
левой НИИР, Санкт-Петербург; ст. науч. сотрудник; E-mail: yuriyan@list.ru
Рекомендована Институтом
Поступила в редакцию 21.05.10 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
УДК 621.391
Ю. А. НИКИТИН
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ФОРМИРОВАНИЯ КОЛЕБАНИЙ С ИСПОЛЬЗОВАНИЕМ МЕТОДОВ ПАССИВНОГО ЦИФРОВОГО СИНТЕЗА
Рассмотрена математическая модель, предназначенная для анализа и расчета спектров двух- и многоуровневых колебаний с равномерной и неравномерной дискретизацией по фазе, а также многоуровневых колебаний, характеризуемых конечной точностью выбора узлов аппроксимации. Получены аналитические выражения для спектров колебаний с различными огибающими.
Ключевые слова: пассивный цифровой синтез, накапливающий сумматор, счетчик импульсов, делитель с дробно-переменным коэффициентом деления.
Методы двух- и многоуровневого пассивного цифрового синтеза (ПЦС) частот обеспечивают выполнение целого ряда жестких и взаимно противоречивых требований, предъявляемых к современной аппаратуре синтеза частот и обработки сигналов. К таким требованиям относят высокую долговременную стабильность фазы выходного колебания и ее преемственность при смене частот, малое время переключения с частоты на частоту, низкий относительный уровень дискретных побочных спектральных составляющих, малый уровень фазовых шумов.
Однако на спектр формируемого колебания существенным образом влияет алгоритм работы системы ПЦС как (цифрового) конечного автомата (КА). Наилучший спектр обеспечивается посредством методов оптимального синтеза, использование которых в случае двухуровневого ПЦС [1] позволяет аппаратно реализовать конечный автомат в виде накапливающего сумматора (НС) и/или делителя с дробно-переменным коэффициентом деления (ДДПКД), а в случае многоуровневого ПЦС — только НС, дополненный цифроаналоговым преобразователем или управляемым устройством задержки [2].
Целью настоящей статьи является описание работы конечного автомата применительно к теории синтеза частот и рассмотрение математической модели формирования периодического колебания с огибающей произвольной формы.
Рассмотрим структурные схемы некоторых простейших КА, применяемых в технике синтеза частот, и описывающие их работу алгоритмы [1—5].
Наиболее широко в системах активного и пассивного цифрового синтеза частот используются счетчики импульсов (СИ) (могут содержать поглотители импульсов) и накапливающие сумматоры (рис. 1, а, б соответственно).
а) б)
Fвых 1
Sk K fвх
Fвых
Q
Sk Р fвх
1 СИ
Q НС
Рис. 1
Функция выхода (сигнал переполнения) ρk СИ определяется как
ρk = ⎣k/P⎦ – ⎣(k–1)/P⎦; ρk ∈ (0, 1),
(1)
где ⎣*⎦ — оператор выделения целой части числа, меньшей или равной ему; при этом дли-
тельность такта на входе КА равна Т0.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Математическая модель формирования колебаний с использованием методов ПЦС
53
Функция переходов (текущая сумма) Sk для СИ определяется следующим образом:
Sk = P { k/P },
(2)
где {∗} — оператор выделения дробной части числа, 0 ≤ { k/P } < 1 [5]. Изменение коэффициента пересчета СИ с дробно-переменным коэффициентом деления (N)
осуществляется с помощью поглотителя импульсов и накапливающего сумматора (рис. 2). Такое решение КА широко применяют в системах активного цифрового синтеза частот на основе колец импульсно-фазовой автоподстройки частоты. В последние годы в зарубежной литературе накапливающие сумматоры также называют аккумуляторами фазы и ∆Σ-модуляторами: первое название пришло из теории и техники пассивного (прямого) цифрового синтеза (Direct Digital Synthesis — DDS), а второе — из теории и техники аналого-цифрового преобразования.
Из двух последовательно включенных НС первого порядка (рис. 3) состоит накапливающий сумматор второго порядка — НС2 (и т.д.) (рис. 4).
⎣(kt)/T0⎦
ДДПКД
P 1/⎣Ν⎦
NT0
Q ΣP
HC1 Sk
ρk {kQ/P}
D Такты
Q Sk
ρk =0; 1
Sk+1
Рис. 2
Рис. 3
Q ΣP
SL D
HC2
ΣP Sk
D
ρk =–1; 0; +1; +2
Y1 Σ
+ +
Сигналы управления –
на СИ (ДДПКД) Такты от ДДПКД Y2D
SL+1 Y2
D
D Sk+1
Рис. 4
Если принять емкость НС первого порядка (n =1) равной Р единиц (см. рис. 3), то на его входе можно записать число Q ∈ (0, 1, 2, … Р–1).
Функцию ρk для НС1 можно представить как
ρk = ⎣kQ/P⎦ – ⎣(k–1) Q/P⎦; ρk ∈ (0, 1),
(3)
а функцию Sk — как
Sk = P{kQ/P}.
(4)
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
54 Ю. А. Никитин
Сравнивая выражения (1) и (3), можно сделать вывод, что СИ является частным случаем
НС1 при Q = 1. В табл.1 приведены примеры записи функций переходов (текущей суммы Sk) и функции
выхода (сигнала переполнения ρk) для вариантов НС1 различной емкости Р.
Таблица 1
Функция 0
1
Номер такта 23456
7
8
Р = 7; Q = 2
Sk 0 2 4 6 1 3 5 0 2
ρk 1 0 0 0 1 0 0 1 0 Р=7; Q =3
Sk 0 3 6 2 5 1 4 0 3 ρk 1 0 0 1 0 1 0 1 0
Р=8; Q =3
Sk 0 3 6 1 4 7 2 5 0 ρk 1 0 0 1 0 0 1 0 1
В табл. 2—5 приведены примеры записи функций переходов и выходов для вариантов
НС2 различной емкости Р; в таблицах приняты следующие обозначения: SL — текущая сумма 1; SL+1 — текущая сумма 1 при задержке (D) на один такт; Sk — текущая сумма 2; Sk+1 — текущая сумма 2 при задержке на один такт; Y1, Y2 — сигналы управления; Y2D — сигнал управления Y2 при задержке на один такт; ρk — сигнал управления СИ.
Р = 5; Q = 2
Таблица 2
Функция 0
1
2
3
4
Номер такта 5678
9
10 11 12
SL 2 4 1 3 0 2 4 1 3 0 2 4 1 SL+1 0 2 4 1 3 0 2 4 1 3 0 2 4 Sk 2 1 2 0 0 2 1 2 0 0 2 1 2 Sk+1 0 2 1 2 0 0 2 1 2 0 0 2 1 Y1 1 0 0 1 0 1 0 0 1 0 1 0 0 Y2 0 1 0 1 0 0 1 0 1 0 0 1 0 Y2D 0 0 –1 0 –1 0 0 –1 0 –1 0 0 –1 ρk 1 1 –1 2 –1 1 1 –1 2 –1 1 1 –1
Р = 7; Q = 3
Таблица 3
Функция 0 1 2
Номер такта 3456 7
8
9
10
SL 3 6 2 5 1 4 0 3 6 2 5 SL+1 0 3 6 2 5 1 4 0 3 6 2 Sk 3 2 4 2 3 0 0 3 2 4 2 Sk+1 0 3 2 4 2 3 0 0 3 2 4 Y1 1 0 0 1 0 1 0 1 0 0 1 Y2 0 1 0 1 0 1 0 0 1 0 1 Y2D 0 0 –1 0 –1 0 –1 0 0 –1 0 ρk 1 1 –1 2 –1 2 –1 1 1 –1 2
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Математическая модель формирования колебаний с использованием методов ПЦС
Р = 8; Q = 3
Таблица 4
Функция 0
1
2
3
4
5
6
Номер такта 7 8 9 10 11 12 13 14 15 16 17
SL 3 6 1 4 7 2 5 0 3 6 1 4 7 2 5 0 3 6
SL+1 0 3 6 1 4 7 2 5 0 3 6 1 4 7 2 5 0 3
Sk 3 1 2 6 5 7 4 4 7 5 6 2 1 3 0 0 3 1
Sk+1 0 3 1 2 6 5 7 4 4 7 5 6 2 1 3 0 0 3
Y1 1 0 0 1 0 0 1 0 1 0 0 1 0 0 1 0 1 0
Y2 0 1 0 0 1 0 1 0 0 1 0 1 1 0 1 0 0 1
Y2D 0 0 –1 0 0 –1 0 1 0 0 –1 0 –1 –1 0 –1 0 0
ρk 1 1 –1 1 1 –1 2 –1 1 1 –1 2 0 –1 2 –1 1 1
55
Р = 10; Q = 3
Таблица 5
Функция 0
1
2
3
4
5
6
7
8
Номер такта 9 10 11 12 13 14 15 16 17 18 19 20
SL 3 6 9 2 5 8 1 4 7 0 3 6 9 2 5 8 1 4 7 0 3
SL+1 0 3 6 9 2 5 8 1 4 7 0 3 6 9 2 5 8 1 4 7 0 Sk 3 9 8 0 5 3 4 8 5 5 8 4 3 5 0 8 9 3 0 0 3
Sk+1 0 3 9 8 0 5 3 4 8 5 5 8 4 3 5 0 8 9 3 0 0
Y1 1 0 0 0 1 0 0 1 0 0 1 0 0 0 1 0 0 1 0 0 1
Y2 0 0 1 1 0 1 0 0 1 0 0 1 1 0 1 0 0 1 1 0 0
Y2D 0 0 0 –1 –1 0 –1 0 0 –1 0 0 –1 –1 0 –1 0 0 –1 –1 0
ρk 1 0 1 0 0 1 –1 1 1 –1 1 1 0 –1 2 –1 0 2 0 –1 1
В табл. 1—5 фоном выделен период повторения цикла счета.
Для анализа спектра выходных колебаний КА, дискретных во времени с шагом T0 = 1/f0
и/или квантованных по амплитуде с шагом s, в работе [6] рассмотрена математическая мо-
дель (рис. 5), пригодная для формирования цифрового колебания (Ав(t)) с периодической оги-
бающей любой формы. Модель позволяет учесть конечную точность воспроизведения оги-
бающей при использовании ограниченного набора строго периодических колебаний разной
амплитуды и частоты следования, причем на периоде неравномерности структуры выходного
потока импульсов Tс = 1/Fс укладывается целое число периодов вспомогательных колебаний
модели; здесь Tс =QTв = РT0, Tв = 1/Fвых. Аналогичная структура, пригодная для компьютерно-
го моделирования и анализа спектров цифровых двухуровневых колебаний, рассмотрена в
работе [7].
A0(t)
Идеальный
Aв(t)
формирующий
A2(t) фильтр
A1(t)
A=1
ИФФ
A3(t)
T0
Рис. 5
Например, для модели формирования цифровой моногармоники (ЦМ) с равномерной
дискретизацией по времени T0 целесообразно выбрать следующие исходные строго периодические колебания:
A0(t) — периодическую последовательность δ-импульсов с периодом следования T0; A1(t) = sin(2πQt) — моногармоническую огибающую требуемой (синтезируемой) частоты;
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
56 Ю. А. Никитин
A2(t) = A3(t) = 0; ИФФ — идеальный формирующий фильтр с шириной окна, равной T0. Для спектра ЦМ можно записать
AЦМ (n)
=
(−1)n P
sin(πQ / P) π(nP ± Q)
,
где n = 0,1, ...
Аналогично, если принять
A
ЦТ
(t)
=
⎧1 + ⎨⎩1 −
X X
/π /π
для для
− π ≤ X ≤ 0; 0 ≤ X ≤ π;
A2(t) = A3(t) = 0, X = t/T0, PT0/Q = π, 0 ≤ t ≤ T0, получим модель формирования цифрового треугольного колебания (треугольника — ЦТ) с
равномерной дискретизацией по времени T0, для спектра которого можно записать
AЦТ (k)
=
2λµ
sin[π(Q ± k) / πP(Q ± k) sin2[π(R
P] + 1)
/
, P]
где µ = 1 ∀ P ≡ 0 (mod 2); µ = cos[π(R +1) / P] ∀ P ≡ 1 (mod 2) ; λ = 1cos[π(R +1)] ; k = 0, 1, 2, ..., Q–1;
R = (–1)r–1 (±k)Pr–1; здесь Pr–1 — числитель предпоследнего (r–1)-го члена разложения коэффициента деления N = 2P/Q в цепную дробь по алгоритму Евклида [5].
Одна из важных характеристик выходных колебаний КА, применительно к задачам синтеза частот, — уровень E полезной составляющей спектра, отнесенный к амплитуде A моногармонического колебания требуемой частоты: G = 20lg(E/A). Параметр G позволяет оценить степень приближения формируемого конечным автоматом цифрового колебания к своему аналоговому прототипу. На рис. 6 показан график изменения параметра G в зависимости от коэффициента N = P/Q деления (частоты) КА для разных огибающих цифрового колебания: идеальной моногармоники с равномерной дискретизацией по времени, псевдо- и квазимеандра, цифрового треугольника с равномерной дискретизацией по времени [6, 7]. Здесь величины P и Q выражают соответственно номинальные значения частот опорного f0 и синтезируемого (выходного) Fвых колебаний в единицах общей меры Fс = f0/P = Fвых/Q — шага сетки синтезируемых частот.
G, дБ 0
Идеальная моногармоника
–3,9 Идеальный меандр
–7,8 Идеальный треугольник
Цифровой треугольник
Псевдо- и квазимеандр
Цифровая моногармоника
Рис. 6
1/N = Q/P
Предложенная математическая модель позволяет производить численный расчет спектров формируемых конечным автоматом двух- и многоуровневых колебаний с равномерной и неравномерной дискретизацией по фазе, а также многоуровневых колебаний, характеризуемых конечной точностью выбора узлов аппроксимации, что, в свою очередь, позволяет получать аналитические выражения в свернутом виде для спектров колебаний с различной огибающей.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Математическая модель формирования колебаний с использованием методов ПЦС
57
СПИСОК ЛИТЕРАТУРЫ
1. Никитин Ю. А. Широкополосный синтез частот с помощью цифровых структур // Изв. вузов СССР. Приборостроение. 1990. Т. 33, № 9. С. 39—47.
2. Гуревич И. Н., Никитин Ю. А. Синтез сетки частот с произвольным шагом // Радиотехника. 1992. № 4. С. 53—58.
3. Брауэр В. Введение в теорию конечных автоматов: Пер. с нем. М.: Радио и связь, 1987. 392 с.
4. Апериодические автоматы / Под ред. В. И. Варшавского. М.: Наука, 1976. 424 с.
5. Виноградов И. М. Основы теории чисел. М.: Наука, 1972.
6. Никитин Ю. А. Математическая модель формирования колебаний с периодической огибающей // Тр. НИИР: Сб. статей. М., 1999. 124 с.
7. Никитин Ю. А. Конечный автомат как элемент цифровой системы синтеза частот // Докл. X Междунар. науч.-техн. конф. „Радиолокация, навигация, связь“. Воронеж, НПФ „САКВОЕЕ“ ООО, 2004. Т. 1. С. 526—533.
Юрий Александрович Никитин
Сведения об авторе — канд. техн. наук; Филиал ФГУП НИИ радио — Ленинградский отрас-
левой НИИР, Санкт-Петербург; ст. науч. сотрудник; E-mail: yuriyan@list.ru
Рекомендована Институтом
Поступила в редакцию 21.05.10 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9