Радиоинтерферометр с автономными гетеродинами в разнесенных приемных каналах

Размещено на http: //www. allbest. ru/

Технологический институт Южного Федерального Университета

Радиоинтерферометр с автономными гетеродинами в разнесенных приемных каналах

А.П. Дятлов, Д.О. Троц



Аннотация

Получена 22 июля 2013 г.

Исследуются принципы построения радиоинтерферометра и стенда для его моделирования с целью обеспечения оптимизации основных характеристик при различных исходных условиях.

Ключевые слова: радиоинтерферометр, линейный тракт приемника, пеленг, автономный гетеродин, шумовой сигнал, база пеленгования, отношение сигнал/помеха, погрешность пеленгования.



Abstract

Principles of constructing of radio interferometer and test bench for his modeling to ensure optimization of the basic characteristics for different initial conditions are considered.

Keywords: radiointerferometer, linear receiver chain, bearing, autonomousheterodyne, the noise signal, the base direction finding, signal-to-noise ratio, accuracy of direction finding.



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

В данной статье с целью дальнейшего снижения погрешности пеленгования исследуются принципы построения РИ с автономными гетеродинами в каналах линейного тракта приемного устройства.

С целью уменьшения погрешности пеленгования в РИ с большими базами в работе [2] было предложено использование независимых атомных стандартов частоты с фиксацией выходных эффектов на выходе каналов в запоминающих устройствах и последующей постобработкой на основе совмещения записей для определения пеленга.

Однако такой подход имеет такие недостатки, как высокая стоимость РИ, низкое быстродействие и возможность пеленгования в узком частотном диапазоне. радиоинтерферометр пеленгование гетеродин погрешность

При необходимости функционирования РИ в широком частотном диапазоне предлагается использовать в качестве гетеродинов в приемных каналах двух независимых синтезаторов частоты, перестраиваемых по частоте от общего пункта управления (ПУ). С учетом современного уровня элементной базы, используемой при проектировании перестраиваемых широкодиапазонных синтезаторов частот имеется возможность реализации в каналах ЛТП гетеродинных напряжений с разносом частот до 500 кГц при перепаде температур от -20° С до +30° С [3].

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

С целью обеспечения функционирования в реальном масштабе времени между каналами приемника следует предусмотреть каналы связи по промежуточной частоте, а для уменьшения аппаратурной погрешности, обусловленной нестабильностью группового времени запаздывания в каналах связи предлагается использовать компенсационный алгоритм [3].

Рис. 1

С учетом вышеизложенного структура исследуемого РИ приведена на рис. 1, где А1, А2 -- антенны; ЛТП1, ЛТП2 -- линейные тракты приемника; Г1, Г2 -- гетеродины; М1, М2 -- модуляторы; УСВК1, УСВК2 -- устройства совмещенного встроенного контроля; Обн1, Обн2 -- обнаружители; РЛЗ1, РЛЗ2 -- регулируемые линии задержки; П1, П2 -- перемножители; ПСА1, ПСА2 -- параллельные спектроанализаторы; ПУ -- пункт управления; Упр1, Упр2 -- управители; ЛС1, ЛС2 -- линии связи; Ш1, Ш2,Ш3 -- шины управления; К1, К2 -- каналы РИ.

Принцип действия РИ состоит в следующем. При пеленговании по азимуту сектор пеленгования определяется шириной диаграммы направленности антенн А1 и А2 в азимутальной плоскости a. В случае пеленгования источников радиоизлучений в дальней зоне, когда фронт приходящих электромагнитных волн можно полагать плоским, на входы каналов РИ поступают процессы

y1(t) = S(t) + n1(t); y2(t) = S(t - s) + n2(t); s = ,

где S(t) -- шумовой сигнал (ШС), соответствующий квазибелому шуму (КШ); n1(t), n2(t) -- внутренние шумы; s -- временной сдвиг ШС, пропорциональный пеленгу ; d -- разнос антенн А1 и А2 (база РИ); c -- скорость распространения радиоволн.

Автокорреляционная функция ШС имеет вид:

Rs() =sinc(fs) cos(2fs),F = fв - fн,

где  -- дисперсия сигнала S(t); fs, fs -- несущая частота и ширина спектра ШС; fн, fв, F -- нижняя, верхняя границы и ширина рабочего частотного диапазона РИ.

Поиск ШС по частоте в ЛТП1 и ЛТП2 осуществляется синхронно и завершается обнаружением ШС после перестройки частоты Г1 и Г2, когда

fs - fos1 = fs1; fs - fos2 = fs2; fos1 - fos2 ? F; fn ?fs; F<fn;

где fos1, fos2 -- частоты гетеродинов Г1 и Г2; fn -- полоса пропускания ЛТП1 и ЛТП2; F -- величина расстройки частот гетеродинов Г1 и Г2 при их автономной перестройке по целеуказаниям от ПУ; fs1, fs2 -- частота ШС после преобразования частоты в ЛТП1 и ЛТП2.

После прохождения процессов y1(t) и y2(t) через ЛТП1 и ЛТП2 с учетом аппаратурных погрешностей, обусловленных неидеальностью их фазочастотных характеристик, получаем:

y11(t) = K1[S1(t - 1) + n1(t - 1)]; y21(t) = K2[S2(t - s- 2) + n2(t - s- 2)], Д0=1-2;

где K1, K2 -- коэффициенты передачи ЛТП1 и ЛТП2; 1, 2 -- групповое запаздывание, вносимое ЛТП1 и ЛТП2; Д0-неидентичностьЛТП1 и ЛТП2 по групповому запаздыванию.

При использовании в РИ совмещенного встроенного контроля, как показано в работе [1], обеспечивается благодаря УСВК1 и УСВК2 корректировка группового запаздывания в ЛТП1 и ЛТП2, таким образом, чтобы 1 = 2.

В Обн1 и Обн2 осуществляется энергетическое обнаружение, если

U01(T0) = >Uп; U02(T0) = >Uп,

где U01(T0), U01(T0) -- напряжение на выходе энергетических обнаружителей Обн1 и Обн2; T0 -- постоянная интегрирования в обнаружителе; Uп -- пороговое напряжение.

При fпT0>> 1 характеристики помехоустойчивости рассчитываются из следующих соотношений [5]:



Pпо = Ф; Ф(x) = ; Pлт = 1 - Ф;

 = ;  = ,

где Pпо, Pлт -- вероятность правильного обнаружения и ложной тревоги; Ф(x), arc Ф(x) -- функция Лапласа и обратная ей функция;  -- входное отношение сигнал/помеха по мощности;  -- отношение сигнал/помеха по напряжению на выходе Обн1(2); gп -- нормированный порог; -- дисперсия шума в полосе fп.

Максимальное время поиска сигнала S(t) в ЛТП1(2) по частоте и его обнаружения  равно

 = .

После обнаружения ШС начинается этап подстройки РЛЗ1 и РЛЗ2 с целью обеспечения оценки временного сдвига s. При этом отличительной особенностью исследуемого РИ является необходимость перехода от корреляционной к корреляционно-фильтровой обработке. Учитывая, что средние частоты сигналов S1(t) и S2(t) отличаются на априорно неизвестную величину F, то для выделения свернутого по спектру выходного эффекта используется ПСА, частотный диапазон которого перекрывает интервал изменения F.

После завершения подстройки РЛЗ1 и РЛЗ2 на выходе перемножителя П1(2) одного из каналов ПСА1(2), центральная частота которого совпадает с разностной частотой F, выделяется компонент «сигнал-сигнал»

S12(t) = ;

S22(t) = ;

hi(t) = fksinc(t) cos(2t);  ? F,

где l, 2 -- временные сдвиги, вносимые РЛЗ1 и РЛЗ2; 11, 12 -- временные сдвиги, вносимые линиями связи ЛС1 и ЛС2; hi(t) -- импульсная реакции i-го канала ПСА со средней частотой  и полосой пропускания ?

При

1 = s - l2 и 2 = -s + 11

на выходе i-го канала ПСА1(2) компоненты «сигнал-сигнал» представляет собой квазигармонические процессы

S13(t) = Kпrs(1) cos(2Ft); S14(t) = Kпrs(2) cos(2Ft); rs() = sinc(рfs),

где Kп -- коэффициент передачи П1(2); rs(1), rs(2) -- коэффициенты автокорреляции сигнала в 1-м и 2-м каналах РИ; 1, 2 -- временные рассогласования за счет неточности подстройки РЛЗ1 и РЛЗ2, а также неидеальности функционирования УСВК1 и УСВК2.

При использовании в качестве линий связи комбинированных многожильных кабелей с компенсирующими усилителями, например, кабель марки SAT100n, можно обеспечить:

1) равенство задержек, вносимых ЛС1 и ЛС2 (11 = 12);

2) передачу в одной оболочке как радиосигналов (ЛС1 и ЛС2), так и управляющих сигналов (Ш1, Ш2, Ш3);

3) возможность увеличения базы пеленгования d до десятков километров, поскольку в данном варианте РИ радиосигналы передаются на промежуточной частоте и кабели вносят затухание не более нескольких дБ на 1 км.

С целью уменьшения аппаратурной погрешности пеленгования, вносимой нестабильностью l1 и l2, оценки временных сдвигов  и , фиксируемые в ПУ, используются в компенсационном алгоритме для получения оценки временного сдвига  и оценки пеленга

= ;  =.

С учетом вышеизложенного среднеквадратичная погрешность пеленгования зависит от таких параметров как значение пеленга , база пеленгования d и среднеквадратичная погрешность оценивания временного сдвига сигналов

= ; d = ; ? ,

где d -- время прохождения ШС расстояния, соответствующего базе пеленгования d;  - сектор пеленгования.

Погрешность зависит от флюктуационной погрешности 1, аппаратурной погрешности 2, методической погрешности 3



= .

Среднеквадратичная флюктуационная временная погрешность 1 определяется из соотношений [2]

1= ; fkTб = 1; fk = ; g =  при < 1; > 1,

где -- отношение сигнал/помеха по напряжению на выходе i-го полосового фильтра канала ПСА1(2); g -- отношение сигнал/шум по напряжению на выходе i-го канала ПСА1(2); fa -- рабочий частотный диапазон ПСА1(2); fk -- полоса пропускания каждого канала ПСА1(2); nk -- количество каналов в ПСА1(2); T -- постоянная интегрирования на выходе каналов ПСА1(2).

Среднеквадратичная аппаратурная погрешность 2 определяется неидеальностью функционирования УСВК1(2)

(р12 = м1 - 2)

и нестабильностью временного сдвига в ЛС1(2)

(ЛС1(2) = l1 - l2)

2= .

Среднеквадратичная методическая временная погрешность 3 при дискретном законе перестройки РЛЗ1(2) зависит от шага временного поиска ш и равна



3=  при ш<ks,

где ks -- интервал корреляции ШС. На этапе пеленгования осуществляется дискретная перестройка РЛЗ1(2) до тех пор пока на частоте F в одном из каналов не появится гармоническое колебание, соответствующее компоненту «сигнал-сигнал». Окончательная подстройка РЛЗ1(2) обеспечивается Упр1(2), куда вводится дифференцирующая цепь для придания дискриминационной характеристике нечетного характера.

Длительность этапа пеленгования равна

Tэ= .

Реальная чувствительность РИ рассчитывается следующим образом

Ps1 = kT0fnNn; Ps2 = kT0fnNn; kT0 = 4·10-21 Вт/Гц,

где Nn -- коэффициент шума ЛТП1(2); Ps1, Ps2 -- реальные чувствительности РИ в режимах обнаружения и пеленгования,,- входное отношение сигнал/помеха по мощности в режимах обнаружения и пеленгования. Длительность одного сеанса функционирования РИ равна

TРИ = Tп + Tэ.

С целью проверки теоретических расчетов при оптимизации основных параметров функциональных узлов и РИ в целом при различных исходных данных описываются принципы построения моделирующего стенда, реализованного на основе продуктов программы системотехнического моделирования SystemView.

Функциональная схема стенда приведена на рис.2

Рис. 2

На входе стенда используется три генератора шума: ФУ 0 для формирования шумового сигнала (ШС) и ФУ 30, 32 для формирования внутреннегошума n(t) РИ. Приемная часть РИ моделируется двумя каналами линейных трактов (ЛТП1 и ЛТП2). ЛТП 1,2 состоит из входных ПФ (ФУ 10, ФУ11), смесителей (ФУ 14, ФУ 15), гетеродинов (ФУ 16, ФУ 17), и полосовых фильтров на промежуточных частотах (ФУ 18, ФУ 19). При этом для имитации аппаратурных погрешностей РИ используется возможность задания на основе меню для моделей ПФ, соответственных форм АЧХ и ФЧХ, а изменением значений частот гетеродинов можно устанавливать различную величину частотного разноса F. Фиксированная линия задержки(ЛЗ) (ФУ 7) позволяет моделировать изменения значений пеленгов.

Коррелятор имеет три выхода и состоит из перемножителя (ФУ 21), регулируемой линии задержки (РЛЗ) (ФУ 20), трех параллельно включенных узкополосных фильтров УФ (ФУ 22, ФУ 23, ФУ 24), трех параллельно включенных детекторов огибающей (ФУ 25, ФУ 26, ФУ 27) и фильтров нижних частот (ФУ 28, ФУ 29, ФУ 30).

Первый выход коррелятора  (ФУ 22) предназначен для фиксации частотного разноса гетеродинов F. Второй выход коррелятора с ( ФУ 23) предназначен для грубой оценки задержки , вносимой ЛЗ (ФУ 7), и соответственно пеленга в режиме поиска на основе последовательной перестройки РЛЗ (ФУ 20) линейно-ступенчатым законом. Третий выход коррелятора с УФ 3 (ФУ 24) и дополнительно включенной дифференцирующей цепью (ФУ 31) предназначено для точного оценивания задержки  и соответственно пеленга  в режиме слежения при использовании дискриминационной характеристики.

Для обеспечения настройки и проверки нормального функционирования всех функциональных узлов и РИ в целом используется набор универсальных анализаторов (ФУ 32 - ФУ 42), обеспечивающих контроль временного и спектрального представления выходных процессов.

На рис 3 а, б, в приведены спектральные распределения ШС на выходе ФУ 10(11), ФУ 18, ФУ 19.

Рис. 3 а,б,в



Рассмотрим пример использования стенда при следующих исходных данных:

1. ШС имеет среднюю частоту  и ширину спектра ;

2. Разнос между частотами гетеродинов (ФУ16, ФУ17) может находится в интервале ;

3. Диапазон поиска по пеленгу ;

4. База РИ ;

5. Диапазон поиска по задержке равен

;

;

6. Частота дискретизации процессов .

В процессе моделирования необходимо осуществить:

1) предварительную оценку пеленга при последовательном поиске путем оценки времени запаздывания , вносимого РЛЗ;

2) оценку разноса частот между гетеродинами ;

3) точную оценку пеленга в режиме слежения путем оценки времени запаздывания , вносимого ЛЗ.

При проведении первого этапа моделирования полагаем ; , а также закон перестройки РЛЗ (ФУ 20) линейно-ступенчатым с параметрами

; ; ; ; Tп=NTш ,



Где  - шаг перестройки РЛЗ;  - длительность шага перестройки РЛЗ; -номер и количество шагов перестройки РЛЗ;  - полоса пропускания ;  - длительность поиска по задержке (пеленгу).

На рис 4 а,б,в приводятся результаты моделирования при установке в стенде , c, N=12; Гц, .

Рис.4 а,б,в

Из рис.4в следует, что при отсутствии помехи n(t) погрешность оценивания временного сдвига  и пеленга  в режиме поиска при б=равна



; ;

При  и d=3 км имеем

.

При проведении второго этапа моделирования полагаем , а границы полосы пропускания  соответствует диапазону изменения разноса F.

Для оценки частот гетеродиновFможно использовать спектральное распределение входного процесса на выходе , которое приведено на рис.5, из которого следует, что свертка спектра сигнала происходит на частоте , соответствующей частному разносу гетеродинов.

Рис.5

При проведении третьего этапа моделирования с выхода ДЦ при подаче а)

y(t)=S(t); б)y(t)=S(t)+n(t)



снимается дискриминационная зависимость

,

; ; ;

где  - полоса пропускания УФ3;  - шаг перестройки РЛЗ в режиме слежения. На рис. 6 приведена дискриминационная характеристика коррелятора при подаче на вход РИ только сигнала.

Рис. 6

На рис. 7 приведена зависимость выходного отношения сигнал/шум по напряжению к входному отношению сигнал/помеха по напряжению , а на рис. 8 дискриминационная характеристика  при .



Рис. 7

Рис.8

При подаче на вход РИ аддитивной смеси y(t)=y(t)+n, как следует из рис.8 , а , при  и , где  - среднеквадратичные погрешности оценивания задержки и пеленга в режиме слежения. Результаты моделирования хорошо согласуются с теоретическими расчетами.

Разработанные в данной работе принципы построения РИ могут быть использованы при построении многоэлементных РИ различного назначения.



Литература

1. Дятлов А.П., Дятлов П.А. Совмещенный встроенный контроль линейных трактов широкодиапазонных супергетеродинных приемников.-М.: Специальная техника, №1, 2010, с.29

2. Губанов В.С, Финкельштейн А.М.. Фридман П.А. Введение в радиоастрономию.-М:Наука, 1983. -280 с.

3. Патент РФ №218874. Корреляционный пеленгатор. Авторы: Дятлов А.П., Евдокимов Ю.Ф.

Размещено на Allbest.ru