Цифровой широкополосный коррелятор Сибирского солнечного радиотелескопа

УДК 520.8.056

ЦИФРОВОЙ ШИРОКОПОЛОСНЫЙ КОРРЕЛЯТОР СИБИРСКОГО СОЛНЕЧНОГО РАДИОТЕЛЕСКОПА

А.В.Губин1, С.В.Лесовой2

Институт солнечно-земной физики СО РАН, 664033, г. Иркутск, ул. Лермонтова,126.

Приведено описание цифрового корреляционного приемника Сибирского солнечного радиотелескопа (ССРТ). Изложены основные принципы измерения спектра мощности при помощи коррелятора задержек, реализованного на программируемых логических интегральных схемах (ПЛИС). Показана возможность измерения с помощью ПЛИС спектра мощности сигналов с полосой частот шире, чем полоса частот самой ПЛИС. Это достигается путем демультиплексирования входных данных и параллельной обработки полученных цифровых потоков. Рассмотрены характеристики описываемого коррелятора как анализатора спектра мощности и как составной части приемной системы ССРТ. Приведены примеры наблюдений Солнца на ССРТ с помощью корреляционного приемника.

Ил. 6. Библиогр. 7 назв.

Ключевые слова: радиотелескоп; коррелятор; спектр мощности; демультиплексирование.

DIGITAL WIDEBAND CORRELATOR OF THE SIBERIAN SOLAR RADIO TELESCOPE A. V. Gubin, S.V. Lesovoi

Institute of Solar-Terrestrial Physics SB RAS, 126 Lermontov St., Irkutsk, 664033.

The article describes a digital correlation receiver of the Siberian Solar Radio Telescope (SSRT). It provides basic principles of measuring the power spectrum by means of the delay correlator, implemented on programmable logic integrated circuits (PLIC). The authors demonstrate the possibility to use PLIC-spectrum for measuring signal power with a frequency band wider than the frequency band of the PLIC itself. This is achieved by demultiplexing of the input data and the parallel processing of the obtained digital streams. The characteristics of the described correlator as a power spectrum analyzer and as a component of the SSRT receiving system are examined. Examples of sun observations by the SSRT with the help of the correlation receiver are given. 6 figures. 7 sources.

Key words: radio telescope; correlator; power spectrum; demultiplexing.

Задача по измерению спектра мощности сигнала или взаимного спектра мощности двух сигналов является одной из основных для приемных устройств, используемых в радиоастрономии. С развитием цифровых методов обработки сигналов измерение спектра мощности корреляционным способом становится всё более привлекательным, т.к. достигается максимальное перекрытие выборок анализируемого сигнала, а значит, и максимальная чувствительность. Наиболее подходящим для этой задачи является реализация на программируемых логических интегральных схемах, которые сочетают в себе легкую конфигурируемость -гибкость и высокую скорость обработки данных.

В общем случае кросскорреляция двух сигналов и есть функция их относительного сдвига

по времени: с (г) = ^//^^(О 52( t + г) ей, где г - временной сдвиг между сигналами, а время выборки Т = т г определяется числом отсчетов т для дискретных сигналов [1]. Фурье-преобразование с (г) и есть взаимный спектр мощности сигналов и 52(0. Если на оба входа коррелятора подать один и тот же

сигнал, то получим автокоррелятор. В этом случае Фурье-преобразование выходных данных коррелятора, согласно теореме Винера-Хинчина, дает спектр мощности или энергетический спектр входного сигнала. Измерение корреляции в ограниченном диапазоне задержек эквивалентно умножению корреля-

ционной функции на прямоугольную функцию ширины , что в частотной области равносильно свертке спектра сигнала с Фурье-образом прямоугольной

функции или в/лс-функцией ) , где Л[ = 1/т.

Поэтому спектральные данные, полученные с помощью коррелятора, должны быть проверены на влияние боковых лепестков в/лс-функции на конечный результат.

Спектральный коррелятор, или коррелятор с задержками, реализован по схеме, представленной на рис. 1. На один вход коррелятора подается сигнал через линию задержки, которая представляет собой сдвиговый регистр. Задержка в данном случае определяется тактовой частотой г = и количеством

1Губин Алексей Владимирович, младший научный сотрудник; тел.: 89086473888, е-mail: gubin@iszf.irk.ru

Gubin Aleksei, radio astronomy, Junior Research Worker, tel.: 89086473888, e-mail: gubin@iszf.irk.ru

2Лесовой Сергей Владимирович, кандидат физико-математических наук, заведующий лабораторией, тел.: (3952) 564507,

е-mail: svlesovoi@gmail.com

Lesovoi Sergey, Candidate of Physico- Mathematical sciences, Head of the Laboratory, tel.: (3952) 564507, e-mail: svlesovoi@gmail.com

Рис. 1. п-канальный коррелятор с задержками

ячеек регистра. Выход линии задержки подключен непосредственно к перемножителям, число которых определяет число каналов коррелятора п. Длина данного сдвигового регистра определяет положение, соответствующее г = 0, т.е. при количестве ячеек в регистре, равном половине каналов п/2, значение, соответствующее г = 0 , будет находиться в середине выходного массива данных коррелятора.

Второй вход коррелятора представляет собой сдвиговый регистр с параллельным выходом. Каждый выход регистра соответствует определенной задержке и соединен с соответствующим перемножителем, т.е. длина сдвигового регистра определяется числом перемножителей или числом каналов п. В нашем случае на оба входа подается один и тот же сигнал, т.е. измеряется автокорреляция. Так как функция автокорреляции симметрична относительно г = 0, то для упрощения реализации (уменьшения объема используемой логики) можно получать автокорреляцию в диапазоне задержек от 0 до г п/2 , а недостающую половину дополнять программным способом. Как будет показано ниже, такое уменьшение объема логики позволяет повысить скорость обработки данных при реализации автокоррелятора на ПЛИС. С каждым тактом выход каждого перемножителя добавляется к соответствующему аккумулятору до тех пор, пока не будет выработан сигнал сброса, определяемый временем накопления. По этому сигналу данные фиксируются на выходе аккумулятора, затем аккумулятор обнуляется и весь цикл начинается заново.

Максимальная скорость битовых потоков в интегральных микросхемах, используемых в спектральных

корреляторах, обычно составляет десятки Мбит/с, что по крайней мере на порядок меньше, чем у устройств цифровой выборки. Преобразование выхода устройства выборки из последовательной формы в параллельную, т.е. демультиплексирование во временной области, позволяет получать оптимальные скорости поступления данных на коррелятор. Рассмотрим коррелятор, в котором выход каждого устройства выборки демультиплексируется на N субпотоков, а параллельная архитектура коррелятора обеспечивает обработку N субпотоков одновременно. Любые N смежных отсчетов выборки попадают в разные субпотоки, т.е. в каждый субпоток попадут только ^ые отсчеты. Чтобы получить все произведения для пары входных сигналов, необходимых в корреляторе с задержками, требуется вычислить корреляцию каждого из субпотоков одного сигнала со всеми субпотоками другого (рис. 2). Определим кросскорреляционную дискретную функцию сигналов длиной п как ссу = (дк+у/к+ь)к, где д и / — дискретные сигналы, Ь = п/2. В случае разбиения входных данных на N субпотоков, логично разбить и значения кросскорреляционной функции: ссу = сс^ , где у = N(5 + а. Очевидно, что для N = 2 можно записать:

сс° сс1

,а1+р

1

где сср и сср - четные и нечетные отсчеты кросскорреляционной функции, а д° , д^ и /0, - четные и нечетные субпотоки входных данных. Далее, если поток разбит на три части, мы можем записать:

Рис. 2. п-канальный коррелятор с задержками с демультиплексированием входных данных во временной

области на три потока субпотока

а\+р

+3 + 1

2 0 1 д<;+р 1 9<;+р+1/

Наконец в случае, когда входные данные разбиты на N субпотоков:

/ сс'° \

ССр =

1 r-r-N-i ,

+Р

а\+р

ЛГ-1

1

%+Ч3

\/2+W

ые р;

> £<+//+1 ' ' ' 9<;+р+1)

При данной схеме демультиплексирования скорость поступления потоков в каждый субкоррелятор уменьшится в N раз относительно исходной скорости цифровой выборки, при этом общее количество ячеек корреляторов возрастет в раз, но длина каждого субкоррелятора уменьшится в N раз по сравнению с исходным. Таким образом, количество логических элементов ПЛИС, необходимых для вычисления всех значений функций корреляции (для всех задержек), увеличивается в раз при демультиплексировании на N субпотоков и соответствующем понижении необходимой скорости работы логических ячеек в раз. Данное свойство коррелятора задержек можно использовать не только для измерения спектра сигналов с полосой частот шире максимальной полосы частот ПЛИС. Появляется возможность «выходить» за пределы одного чипа, реализуя субкоррелятор для каждой пары из потоков в отдельной ПЛИС и увеличивая тем самым число спектральных каналов в раз.

Волноводный тракт ССРТ смонтирован по схеме двоичного дерева [2], сигналы синфазно суммируются со всех 128-ми антенн плеча Восток-Запад (ВЗ), как и плеча Север-Юг (СЮ). Суммарные сигналы от всех антенн плеча ВЗ и плеча СЮ подаются на корреляционный приемник, измеряющий автокорреляции суммы и разности сигналов ВЗ и СЮ. Аналоговая часть выполнена по схеме супергетеродинного приемника с подавлением зеркального канала полосовым волно-водным фильтром, после которого сигнал подается на СВЧ-усилитель. АЧХ СВЧ-фильтра и СВЧ-усилителя, а также положение частоты гетеродина приведены на рис. 3. На входе АЦП формируется сигнал с шириной полосы частот Д/ = 1 3 5 М Гц и центральной частотой / = 6 7 , 5 М Гц.

Для оцифровки сигнала используются АЦП от Analog Device AD5463 на отладочной плате ADS5463 Texas Instruments [3]. Данный АЦП имеет разрядность выходных данных 12 бит и позволяет оцифровывать сигналы на частотах до 500 МГц, среднее значение джиттера 150 фсек. Требуемое питание: +5В для аналоговой части и +3.3 для цифровой. Ввод аналогового сигнала осуществляется через коаксиальный разъем SMA-типа. Выходные данные поступают на высокоскоростной разъем QTH-30X2D-A. Уровни напряжений выходных данных соответствуют стандарту LVDS. Максимальная скорость передачи данных до 800 Мбит/с.

Разрядность входной шины данных составляет 9 бит и обусловлена шумовым поведением младших разрядов АЦП, которые ухудшают качество данных,

Рис. 3. Блок-схема корреляционного приемника ССРТ и амплитудно-частотные характеристики фильтров

с указанием положения частоты гетеродина

Основные характеристики анализатора спектра мощности

Эффективное число каналов 1 m 1 п' МГц 1 max' МГц Af, МГц fn fn - 1, МГц Динамич. диапазон, Дб Временное разрешение, сек

210 5643 5779 136 0,53 26 0,5

понижая соотношение сигнал-шум [4]. Разрядность выхода умножителя понижена с 18 до 12 бит для дальнейшей возможности многократного накопления результата корреляции. При данной разрядности одним 18-битным DSP-блоком формируется восемь перемножителей [5]. Выходные данные аккумулятора имеют разрядность 32 бита.

Дальнейшая обработка сигналов производится на базе отладочного набора Altera Stratix IV GX FPGA, в котором использован чип EP4SGX230, содержащий 228 000 логических элементов и 161 18-битный DSP-блок и являющийся на сегодняшний день одним из самых производительных из доступных ПЛИС [6]. Наличие двух LDVS совместимых разъемов HSMC (High-Speed Mezzanine Connector) делает сбор данных с АЦП наиболее удобным как по скорости, так и по помехоустойчивости, что позволяет вести обработку данных в реальном времени.

Реальная скорость работы ПЛИС зависит от объема используемой логики. При общей загруженности чипа более 67% максимальная скорость умножителя, реализованного на DSP-блоке, составляет около

, т.е. обрабатываемая полоса частот 4 5 М Гц, что в три раза ниже требуемой. Таким образом, необходимо демультиплексирование входных данных как минимум на три потока (рис. 3), что позволит увеличить требуемую полосу частот до . При заданном числе каналов 210 и разбиении данных на три потока получаем, что на один коррелятор требуется 630 перемножителей. В нашем случае необходимо два автокоррелятора и соответственно 1260 перемножителей, что занимает 98% DSP-блоков в используемой ПЛИС EP4SGX230. В ПЛИС имеются логические элементы общего типа, с помощью которых можно реализовать гораздо большее число перемножителей. Однако реализация умножителей на логических элементах характеризуется меньшей примерно на 30% скоростью работы и не позволяет обрабаты-

плексировании входных данных.

Корреляционный приемник работает под управлением встраиваемого процессора NIOS II фирмы Altera, реализованного в этом же чипе [7]. Собственно коррелятор был разработан с помощью пакета Quartus фирмы Altera. Для разработки программного обеспечения использовался пакет для встроенных систем NIOS II IDE. Система с процессором NIOS II - это эквивалент "компьютера в чипе", который содержит процессор, комбинацию периферии и памяти в одной ПЛИС. При этом обеспечивается полный набор функций 32-битной архитектуры ПК, включая обработку 32 прерываний, доступ к множеству внутричиповой периферии по шине Avalon и интерфейсам с внешней памятью. Коррелятор выполнен в виде периферийного устройства, подключенного к процессору через шину Avalon. Данные, накопленные коррелятором, посредством механизма прямого доступа к памяти вводятся в память процессора. Удаленный компьютер, управляющий процессом наблюдений на ССРТ, имеет доступ к режимам работы коррелятора и к его данным посредством сетевого протокола TCP/IP.

Каждый автокоррелятор содержит 210 каналов и измеряет половину значений функции автокорреляции, соответствующую положительным г. Значения, соответствующие отрицательным формируются программно, путем зеркального отображения, затем данные дополняются до 512 точек для выполнения быстрого преобразования Фурье. При частоте оцифровки входных данных расстояние между каналами по частоте в результирующем спектре должно составлять , что подтверждается измерениями, проводимыми в СВЧ диапазоне, когда тестовый СВЧ сигнал подается непосредственно на антенную решетку ССРТ. Измерения остальных характеристик канала проводились в диапазоне промежуточных частот (ПЧ). Получено, что перекрытие частотных каналов имеет место на уровне 0 . 7 5 , уровень боковых ле-

вать требуемую полосу частот даже при демульти- пестков составляет около 20% (рис. 4).

I  с hanne 1,214( 111.28),215( 1 U .8) jjflj-xjl

SO 00 .....

6000 lf\

4000 Ii -

2000 J IL

0 WW ш

07.640 110.240 Частота, МГц I 12.840 15.

Рис. 4. Характеристики 2-х соседних частотных каналов анализатора спектра (сплошная - 214, пунктир - 215 каналов). Каналы перекрываются на уровне 0.75, уровень боковых лепестков составляет примерно 20%. По оси У - выходная мощность в произвольных единицах

Ш ро1чег,сМ58(82МНг),159

^QA -

103 г -

.о5 г/;

101

-20 -10 0 Мощность, Дбм 10 20 |

Рис. 5. Динамическая характеристика 2-х соседних каналов анализатора спектра. Линейный участок составляет 26 дБ. Ограничение сверху обусловлено реализацией цифровой обработки. По оси Y - выходная

мощность в произвольных единицах

На рис. 5 представлена динамическая характеристика одного из частотных каналов автокоррелятора. Мощность тестового сигнала генератора на частоте 8 2 М Гц увеличивалась от -2 0 дБм до 2 0 дБ м. Измеренный на ПЧ динамический диапазон данного приемного устройства составляет порядка 2 6 дБ и ограничивается в данном случае реализацией алгоритма цифровых вычислений.

В цифровом корреляционном приемнике сохранен проектный метод ССРТ получения двумерных изображений - использование синфазно-противофазной модуляции. Для реализации этого метода одновременно измеряется автокорреляция суммы и разности плечей ВЗ и СЮ; соответствующие спектры мощности, полученные программно, записываются в выходной файл в том же формате, что и для хранения данных ССРТ. Это позволяет использовать существующее программное обеспечение ССРТ для обработки данных цифрового корреляционного приемника. Одной из главных причин разработки цифрового корреляционного приемника для ССРТ, наряду с высокой стабильностью цифровых устройств, является отсутствие у такого приемника паразитной амплитудной модуляции, связанной с реализацией синфазно-противофазной модуляции. Традиционные фазовые модуляторы, используемые на ССРТ, представляют

собой волноводные устройства. Точность изготовления таких устройств недостаточна для того, чтобы полностью скомпенсировать амплитудные искажения сигналов при изменении их фазы. При условии, что полезный сигнал от диска спокойного Солнца составляет лишь около 1% от полных сигналов на выходе волноводного тракта ССРТ, уровень паразитной амплитудной модуляции сравним с откликом от диска Солнца. Цифровой корреляционный приемник свободен от паразитной амплитудной модуляции, так как изменение фазы происходит после оцифровки сигналов. Еще одно преимущество цифрового корреляционного приемника - линейность динамической характеристики, которая в акусто-оптическом приемнике не превышает 20 Дб. Но, пожалуй, основным преимуществом в данном случае является идентичность и стабильность характеристик спектральных каналов, т.к. большим недостатком акусто-оптических систем является нестабильность и неравномерность данных характеристик и соответственно потребность в постоянной юстировке характеристик приемника. Пример изображения Солнца, полученного на ССРТ (длина волны 5.2 см), показан на рис. 6. Хорошо видно, что аппаратная функция автокоррелятора не влияет на полученное изображение. В противном случае наблюдались бы вертикальные полосы вблизи ярких источ-

Рис. 6. Изображение Солнца в интенсивности микроволнового излучения, полученное на ССРТ (длина волны 5.2 см) с помощью цифрового корреляционного приемника

ников. С другой стороны, боковые лепестки диаграммы направленности ССРТ отчетливо видны на полученном изображении. Это показывает возможность применения описываемого автокоррелятора для использования в приемной системе ССРТ.

Таким образом, предложенная методика демультиплексирования входных данных позволяет существенно увеличить полосу частот обрабатываемого сигнала и легко выходить за пределы одного чипа, наращивая тем самым число спектральных каналов.

Разработанный цифровой корреляционный приемник предназначен для ССРТ, однако может использоваться в любых задачах для измерения спектра мощности в реальном времени в полосе частот до с максимальным числом частотных каналов 420. Данные параметры могут варьироваться в зависимости от предъявляемых требований изменением прошивки модуля коррелятора.

Библиографический список

1. Томпсон А.Р., Моран Д.М., Свенсон Д.У. Интерферометрия и синтез в радиоастрономии / пер. с англ. под ред. Л.И. Матвеенко. 2-е изд. М.:ФИЗМАТЛИТ, 2003. 258 с.

2. V.V.Grechnev, S.V.Lesovoi [et al.] The Siberian Solar Radio Telescope: the current state of the instrument, observations and data // Sol. physics. 2003. Vol.216, №1-2. Р.239-272 2004

3. ADS5463 - Texas Instruments http://www.ti.com/lit/ds/symlink/ads5463.pdf

4. Thompson A. R. Quantization Efficiency for Eight or More

Sampling Levels. MMA Memo.220, National Radio Astronomy Observatory, Socorro, NM, 1998.

5. DSP Blocks in Stratix IV Devices www.altera.com/literature/hb/stratix-iv/stx4-siv51004.pdf

6. Stratix IV GX FPGA Development Board. Reference Manual www.altera.com/literature/manual/rm_sivgx_fpga_dev_board.pdf

7. Nios II Software Developer's Handbook CTp.40 http://www.altera.com/literature/hb/nios2/n2sw_nii52010.pdf

УДК 54.084;53.083.8

МЕТОДИКА СРАВНИТЕЛЬНОГО ИЗУЧЕНИЯ ОБРАЗЦОВ НАНОМОДИФИЦИРОВАННОГО ЧУГУНА МЕТОДАМИ СКАНИРУЮЩЕЙ ЭЛЕКТРОННОЙ МИКРОСКОПИИ

Н.Н.Иванчик1, С.С.Колесников2, В.В.Кондратьев3

Национальный исследовательский Иркутский государственный технический университет, Физико-технический институт, 664074, г. Иркутск, ул. Лермонтова, 83.

Представлены результаты научных исследований с отработкой основ методологии сравнительного изучения проб черных металлов методами сканирующей электронной микроскопии. Изучаемые пробы представлены двумя типами чугуна - наномодифицированным и обычным. При исследовании использованы возможности систем регистрации упругорассеянных электронов, энергодисперсионного анализа и системы картрирования примесей на поверхности проб. Ил. 5. Библиогр. 8 назв.

Ключевые слова: электронная сканирующая микроскопия; методика изучения материалов; наноструктуриро-ванные материалы; наночастицы; характеристическое рентгеновское излучение.

METHODOLOGY AND PRACTICE OF THE COMPARATIVE STUDY OF CAST IRON-WITNESS AND MODIFIED CAST IRON SAMPLES BY THE METHODS OF SCANNING ELECTRONIC MICROSCOPY N.N. Ivanchik, S.S. Kolesnikov, V.V. Kondratyev

National Research Irkutsk State Technical University, Physico-Technical Institute, 83 Lermontov St., Irkutsk, 664074

The article presents the results of researches where the basics of the methodology of comparative study of ferrous metals samples by the methods of scanning electronic microscopy are worked off. The studied samples are represented by two types of cast iron - nanomodified and normal. The research used the possibilities of the systems for recording elas-tically scattered electrons, energy and dispersive analysis, and the system for mapping impurities on sample surfaces. 5 figures. 8 sources.

Key words: electronic scanning microscopy; methods to study materials; nanostructured materials; nanoparticles; characteristic X-rays.

1Иванчик Николай Николаевич, студент, техник отдела инновационных технологий, тел.: 89500952880, e-mail: kolan2004@yandex.ru

Ivanchik Nikolai, Student, Technician of the Department of Innovative Technologies, tel.: 89500952880, e-mail: kolan2004@yandex.ru

2Колесников Сергей Сергеевич, ведущий инженер отдела синтеза наноструктур, тел.: 89501265890, e-mail: kvv@istu.edu Kolesnikov Sergey, Leading Engineer of the Department of Synthesis of Nanostructures, tel.: 89501265890, e-mail: kvv@istu.edu

3Кондратьев Виктор Викторович, кандидат технических наук, начальник отдела инновационных технологий, доцент кафедры квантовой физики и нанотехнологий, тел.: 89025687702, e-mail: kvv@istu.edu

Kondratyev Victor, Candidate of technical sciences, Head of the Department of Innovative Technologies, Associate Professor of the Department of Quantum Physics and Nanotechnologies, tel.: 89025687702, e-mail: kvv@istu.edu