УДК 621.391
В.В. Котов, С.И. Ширшин, В.Н. Бурканов ЦИФРОВОЙ ПРЕОБРАЗОВАТЕЛЬ СИГНАЛОВ
Рассматриваются способы реализации устройств обработки сигналов с применением цифрового квадратурного модулятора и формирователя комплексного сигнала на основе дискретного преобразования Гильберта.
V.V. Kotov, S.I. Shirshin, V.N. Burkanov SIGNALS DIGITAL CONVERTER
The article considers the ways of realization of devices of processing of signals with using the digital QM modulator and discrete shaper of code signal based on Gilbert discrete converter.
Устройства с квадратурными модуляторами (демодуляторами) являются универсальными устройствами, предназначенными для обработки сигналов независимо от вида модуляции. Существуют интегральные микросхемы аналоговых и цифровых квадратурных модуляторов, например фирмы Analog Device.
В данной работе решается задача смещения во времени и по частоте исходного квазигармонического широкополосного сигнала с финитным спектром. При этом необходимо минимизировать допущения при аналитическом формировании сигнала за один период частоты дискретизации.
Для определенности, задача решалась по следующим исходным параметрам: частота дискретизации f = 90 МГц, частота входного сигнала f0 = 30 МГц, полоса частот входного сигнала Af = 30 МГц (f ± 15 МГц), шаг смещения по частоте не хуже Af < 2f2c •10 6, смещение во времени до 4 мс. Исходя из этого определяются как аналитические, так и схемотехнические подходы к решению задачи.
Для реализации цифровых устройств с использованием свойств аналитического сигнала широкое распространение получили методы формирования последовательности отсчетов сигнала, сопряженного по Гильберту. В одном случае для этого может быть использовано конечное дискретное преобразование Фурье (КДПФ), а в другом -приближенные преобразования на базе дискретного преобразования Гильберта. На пути реализации первого метода необходимо осуществить прямое и обратное КДПФ для перехода из временной области в частотную и наоборот, которое может быть выполнено на основе применения алгоритмов БПФ. Однако использование этого метода для создания устройств с поточной обработкой сигналов встречает трудности реализации преобразователя, работающего в реальном времени, т.е. за время одного периода частоты дискретизации (At), что в нашем случае составляет At = 11 нс.
Второй метод к получению отсчетов входного сигнала, сопряженного по Гильберту, имеет простую реализацию, хорошо подходит для поточной обработки, но основан на использовании конечного числа членов бесконечного ряда. Весовые
коэффициенты ряда убывают медленно, поэтому для того, чтобы достичь хорошей фильтрации, число коэффициентов должно быть очень большим.
В то же время известно
2 3 к„-1 кп к0+1 . Нг конечное дискретное
преобразование Гильберта, представляющее свертку в аналитическом виде. Однако прямое использование такого преобразования для
формирования отсчетов сигнала, сопряженного по Гильберту, сопряжено с большими вычислительными затратами времени и памяти, кроме того ориентировано на блочную обработку сигнала.
В данном случае применим метод вычисления конечного дискретного
преобразования Гильберта, обеспечивающий экономию вычислительных ресурсов и возможность простой аппаратной реализации при поточном преобразовании [1]. На рис. 1
приведена структурная схема преобразователя сигнала.
Сигнал, дискретизованный с шагом &, поступает на сдвиговый регистр. Коммутатор осуществляет подключение к отводам сдвигового регистра в соответствии со смещением во времени. Если время запаздывания т0, то коммутатор считывает последовательность отсчетов с отвода к0, обеспечивая задержку на к0 тактов. Эта последовательность и[п-к0] поступает на формирователь дискретного комплексного сигнала (ДКС), который образует сигнал:
^[п] = и [п] + [п],
где и[п]=Ке(^[п]} - реальная составляющая, равная отсчетам входного сигнала; у[п]=1ш(^[п]} - мнимая составляющая; ] - мнимая единица.
Отсчеты сигнала у[п] определяются выражением [ 1]:
^ _1
у(п) = ^ Ь { и [п + (2/ +1)]_ и [п _ (2/ +1)]}, п = 0,1,... ,
Рис. 1
Ng
(2i +1)—
Ng J
где Ng - размах конечного дискретного преобразователя Гильберта.
Образованный таким образом ДКС имеет свойства аналитического сигнала. Далее, квадратурные составляющие u[n] и v[n] умножаются на cos (p n) и sin (p n) соответственно, где p=oAt (со - круговая частота смещения), а затем суммируются (при
о>0) или вычитаются (при о<0).
2=0
2
Поскольку погрешность квадратурной модуляции связана с искажениями, обусловленными цифровым формирователем квадратурных составляющих, рассмотрим результаты численного анализа формирователя ДКС с помощью методов спектрального анализа для входного сигнала в виде стационарного случайного процесса, которые позволяют оценить искажения АЧХ и вносимый формирователем фазовый шум.
Пусть спектры входного действительного сигнала и[п] и выходного комплексного сигнала м>\р\ в основной полосе частот, полученные по ц-й реализации случайного процесса, определяются ^-точечными конечными дискретными преобразованиями Фурье:
2 ы-- ( 2П Л
ич[к] = тт£и[п]ехР( -7 — .пк\,
N п=о V N )
1 ы~і ( 2п Л
^[к]=тг Ё [п]ехр(- 7^пк I,
N п=о V N )
к = 0,1,..., N/2 -1.
Тогда частотная характеристика формирователя может быть оценена с помощью функции дискретной частоты:
н [к ] = От[к ]/Ои [к ],
1 Ь
где си^[к ] = - Ё К[к ] [к ] -
Ь ц=1
взаимная спектральная плотность;
Ои [к] = — Ё|иц [к]2 - спектральная
Ь ц=1
плотность входного сигнала; Ь -количество реализаций случайного процесса.
Модуль |Н[к] | -
определяет амплитудночастотную характеристику
формирователя, аргумент Н[к], равный
^[к]=агс1в (1т [Н[к]}ГКе[Н[к]} ) -
фазочастотную характеристику.
Для определения функции дискретной частоты Н[к] было проведено статистическое
моделирование на ЭВМ
формирователя с входным сигналом в виде дискретного белого шума:
где Рп - средняя мощность (дисперсия) шума; ^[п] -
нормально распределенные
случайные числа со средним, равным нулю, и дисперсией, равной единице.
Результаты моделировались для размаха преобразователя Ng=32 (см.
Рис. 2. АЧХ а) и ФЧХ б) формирователя ДКС при Ыд=32
а
рис. 2). Можно отметить две особенности амплитудно-частотной характеристики. Во-первых, видно, что имеется область частот, где I Н[к] I близка к единице, как в идеальном случае. Во-вторых, существуют две области частот, в которых дискретная функция I Н[к] I нелинейная и падает значительно ниже единицы. В относительных частотах, нормированных относительно частоты Найквиста п, этими областями являются: //_/П<1/16 и //_/П>15/16. При моделировании с различными значениями область нелинейности уменьшается при увеличении Л§, а при N§=32 амплитудная характеристика линейна на большинстве частот (90% основной полосы частот) с усилением, близким к единице.
Рассмотрение фазочастотных характеристик показывает, что фазовый угол флуктуирует вблизи нуля, при этом величина фазовых флуктуаций в области линейности не превосходит 0.007 рад (и практически не зависит от размаха преобразователя Гильберта).
Структурная схема цифрового поточного формирователя отсчетов комплексного сигнала на основе преобразования Гильберта при N§=32 показана на рис. 3.
Рис. З
Структура цифрового преобразователя сигнала [1] приведена на рис. 4. Схема содержит АЦП, цифровую линию задержки (ЦЛЗ), блок управления, ЦАП.
Блок управления реализован с использованием программируемой логической интегральной схемы D1 (ПЛИС) фирмы Altera и микроконтроллера D2 связи с компьютером. Команды управления, обмен данными, мониторинг платы преобразователя осуществляется от компьютера (по RS232), для чего разработано соответствующее программное обеспечение. На рис. 5 приведено окно программы.
1 I
I ALJII I
I-------------------------
Uin
-e
*-&■
»—*
ADC
CLK
i1 Q » Q
O
i
IJAII l
Uout
-Q-
•»-e-
DAC
OUT Data In
r----------------r
UH3 I
RAM
RAM
I----------
! naBM
h
EPLD
EPLD:
EPLD: I/O
I/O
EPLD:
LOAD
DataDAC
EPLD:
I/O
PortB MC In/Out
Port A PortC
r-------------------------
I Bjiok ynpaBneHM I
. J_______________________
^Управление ЦФОС / DSP v2.0.4.1 / норма
liLJJJDJJj
Настройки
Задержка сигнала
Канал 1 Канал 2
□ программе Помощь
'Щ
3102 ^ нее К Канал 1 -1200 Л. Т- Гц
[2500014 нсек Канал Z 37007 л. ж Гц
✓ о
L13. г, п
Включить AM сигнала по дальности I-Работа
Тест
ЦФОС не отвечает
17 Включить обработку сигнала
Гатов Питание
501
301
10'
Т'С
Инициализация Windows NT, v 5.1.. build 2G00
Выбран порт СОМ1 ск. 2 U1 = 100 U2 = 100 Ч Модуляция откп.
Рис. 5
Экспериментальные исследования [2] показали работоспособность устройства в соответствии с заданными параметрами.
Необходимо отметить тот факт, что применялась и исследовалась доступная элементная база (ПЛИС) с быстродействием 2,8...3 нс, т.е. среднего быстродействия по современным меркам. Таким образом на данный момент можно синтезировать более быстродействующую конфигурацию преобразователя сигналов, работающую с широкополосными сигналами (с полосой не менее 200 МГц), в том числе с применением DSP фирмы Texas Instruments.
ЛИТЕРАТУРА
1. Бурканов В.Н. Цифровая имитация отраженного от цели сигнала / В.Н. Бурканов, В.В. Котов, С.И. Ширшин // Радиотехника и связь: материалы Междунар. науч.-техн. конф. Саратов: СГТУ, 2004. С. 71-81.
2. Коломейцев В.А. Цифровая обработка сигналов в радиотехнике /
В. А. Коломейцев, В. В. Котов // Радиотехника и связь: материалы Междунар. науч.-техн. конф. Саратов: СГТУ, 2005. С. 123-127.
Котов Владимир Васильевич -
ассистент кафедры «Радиотехника»
Саратовского государственного технического университета
Ширшин Сергей Иванович -
доктор физико-математических наук, доцент кафедры «Радиотехника»
Саратовского государственного технического университета
Бурканов Владимир Николаевич -
кандидат технических наук, доцент кафедры «Радиотехника»
Саратовского государственного технического университета