Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы и устройства обработки сигналов в радиотехнических системах

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
4.93 Mб
Скачать

- 50 ~

УДК 621.396,96

Б.А.Федоров, И.Д.Шестаков М.А.Первенцев

ПОСТРОЕНИЕ И АНАЛИЗ ЭМЯШВПОСТИ ДВУШНАЛЬНЫХ PEMKTOPHffi ГРЕБЕНЧАТЫХ ФИЛЬТРОВ

Рассмотрены метод построения и анализ эффективности двухкана­ льного режекторного фильтра (ДРФ) при воздействии пассивных помех *

Ссвременное состояние микроэлектроники позволяет реализовать

цифровые режокторные фильтры (В&) рекурсивной и нерекурсивной стру­ ктуры большей эффективности, чем устройства череспериодной компен­ сации (ЧПК) и рециркуляторы, Целесообразность применения сложных ГФ

определяется в первую очередь эффективностью.

Наибольший вклад в эффективность обнаружения сигналов от дви­ жущихся целей (СДД) на фоне пассивных помех в системах с некогере-

ятиьм накоплением вносит режекторный фильтр. В теории

оптимальной

фмл1тоьции доказано, что

оптимальный фильтр СДД должен

иметь коэф­

фициент передачи К (со)

вида:

 

К (со) = [ S* (и>)/рп ( с о ) ] е ^ с

£

r.ju Sc (со) -

функция, комплексно-сопряженная со спектром обраба­

тываемого

сигнала; Рп (**>)

- энергетический спектр

пассивной

поме­

хи. .I-Hoe

строгое требование в приложениях теории оптимальной

фи­

льтрации

к задачам радиолокации определяет величину

К (со)

кait

функшпе,

имеющую коэффициент передачи, близкий к

единице

и нулю

соответственно

в

областях

средних частот спектра

сигнала

Sc {со) и

помехи Рп (со)

.

На основе

этого положения при построении

ДР£

будем

исходить из критерия максимального отношения сигнал/помеха на выхо­ де РК

Ставится задача повысить эффективность ЙБ типа ЧПК путем вве­

дении дополнительного

Р£, использующего те >ге блоки памяти, как и

и Г 1.1 Таким образом,

имеется два параллельных канала рпжекции

нас-ингой помехи. Предположим, что амплитудно-частотная характери-

•'I HKO ( ЛТ1Х) и

фаз очастотная

характеристика (ФЧХ)

первого

канала

пи'-иг«:<угея известными выражениями:

 

 

^ (со) -

2 п js/n

11

^ (со)- n .(3g-

?

( I )

г ;.:

n. ~ кратность ЧПК; coT мзчшприодный набег фазы помехи за

- 51 -

период Т ЛЧХ и ФЧХ второго дополнительного канала в рампах /1^1

будут синтезированы из условия максимума отношения сигнал/noivy а ни выходе когерентного вычитателя, который объединяет оба параллельных

канала. В этом случае дисперсии сигнала

и помехи е*п на выго­

де вь'читателя описываются соотношениями

[? .]

 

<„ +

?,.г (Т )

а отношение

сигнал/помэха q.

 

 

 

я ~ ( Ы

 

~ д * А ? < 2

<*)»

( 2 )

 

%2

 

гДе <*12 '

е2с '

' %п '

P<2 (th

9<?(г)п - дисперсии сигнала

и помехи на выходах ремекторных каналов и коэффициенты взаимной кор­ реляции между каналами по сигналу и помехе.

Для обеспечения Цнакс необходимо и достаточно, чтобы коэф­ фициенты корреляции каналов по сигналу и помехи были близки к значе­ ниям

 

 

 

 

 

 

 

 

 

( 3)

Условие

(3 )

относительно помехи можно выполнить,

если

АЧХ и ФЧУ

ка­

налов ДРФ будут

вьбрани,как в

[ з ]

 

 

 

 

 

 

 

 

 

* М рпТ)=% Ш ~„т )г

 

(4)

где

Ay

, К2 ,

щ

- АЧХ и ИХ каналов

в спектральной

пло­

тности

мощности помехи

AF^ Т , а относительно

сигнала - соответс­

твенно

 

 

 

 

 

 

 

 

 

 

 

К<Ш СТ )> К

FJA FCT ) * P ,

 

 

(5 )

 

 

 

 

 

 

 

 

 

где

Af^T

-

энергетическая

ширина полосы спектра

сигнала.

 

'При выполнении (4 ) и (5 ) максимальное отношение сигнал/помеха ЫГ». 1*1,1. .1- чР»: ‘ 1*

 

Унакс

/(* « ■ % • )"

 

(6 )

 

*

 

 

 

Технически1реализовать аналитические выражения

(З М 5 ) сложно.

Рассмотрим пример построения ДВЬ, Пусть в качестве первого

канала ДРФ используется

устройство ЧПК-1, АЧХ и ФЧХ которого

опи­

сываются выражением ( I ) .

 

 

 

Второй

канал ДВВпо

условию (4 ) выбираем с

АЧХ вида

A^*J-

* с? / sin 7 }

, а ФЧХ соответственно имеет вид

 

 

- 52

-

Уг (ы) - § - ы Т

(7 )

Структурная схема ДРФ и графические зависимости АЧХ и ФЧХ обоих ка­

цапов

приведены на рис. 1 ,а ,б

соответственно.

Анализ

АЧХ каналов

(рис.

1 ,6 )

в

полосе рекекции

Al'n Т

показал, что

они совпадал»? с

точностью

до

постоянного множителя

Ь = 0 ,5 ,

а ФЧХ

каналов -

не

удовлетворяют условию (4 ) . Для

выполнения условия (4) в первый

ка­

нал режекции последовательно включено нерекурсивное звено первого

порядка на двух блоках памяти с АЧХ и ФЧХ

вида

- {h a 1* 2* м Л *

tpf(uharc tg

а результирующие АЧХ и ФЧХ первого канала описываются выражениями

Кн(со) в 2 s i

n а2* 2а cos2tdT /

 

Ц{(и)

Л -CJ Т

-arctg

а

 

- 2

/ +Qcas2a>T

(8)

а

Рис Л Структурная схема ДРФ и графические зависимости АЧХ и ФЧХ ка­

налов, рассчитанные по выражениям (7 ) , ( 8 ) , приведены на рис. 2 , где показано, что АЧХ и ФЧХ удовлетворяют условиям (4 ) , ( 5 ) .

Структурная схема ДРФ, приведенная на рис. 2 , а , обладает из­ быточностью блоков памяти. Проведя преобразования коэффициентов импульсной характеристики первого канала, получим их значения и структуру ДР^ на трех блоках памяти. Коэффициенты импульсных ха­ рактеристик первого и второго каналов соответственно равны:

V ' -

A js a > 4 а ~а *

V о,

&0 ~ о .

- 53 -

Рис.2

Бсоответствии с коэффициентами импульсных характеристик АЧХ

иФЧХ каналов описываются выражениями

К^ш) « У [ 1-сд$соТ+ aco$2cSf-a cos5uiT]*+[s/nьуТ~аs i n a s i n J c j T ] 3 t

(pj(cu)*arcty[(smct}T-Q2o>T* as(n3b)T)/(i-coscoT + ...)] f

K2(a>)= 3L

j s i n

C J T I f

Ф2(со) = 3 ? - o > T t

(9)

Анализ

ФЧХ

каналов

по выражению (9) показал,

что фазовой ха­

рактеристикой первого канала ДВЕ можно управлять с помощью весово­

го

коэффициента а

, а

АЧХ второго канала ДРФ - с помощью измене­

ния

весового

коэффициента

L . Оба канола ДМ объединяются через

когерентный

вычитатель.

 

 

Эффективность режекции ДМ по коэффициенту прохождения помехи

К„п в

зависимости от ее

корреляционных

свойств проведена аналити­

ческим

путем по известным выражениям

[ 4

] . Анализ результатов

Кпп

, приведенных на рис. 3

при а

= 0,3 6 3 3 , показал, что

при

фиксированной величине

AFH T

имеется

оптимальное значение

L »

при котором величина

имеет минимальное значение, ^являющееся

предельной величиной нарёкзфСй®йог6 фильтра третьего порццка. Ре­

зультаты расчетов

величины Кпп

в зависимости

от JL

и пара­

метра

AFn T

сведены в

таблицу.

 

 

 

 

A Fn Т

0 .0 5

'

0 ,1 .

0 .1 5

0 ,2

0 .25

0 .3

« п т * *

-5 6 ,7 3

-3 9 ,9 1

-2^,51

-2 2 ,2 1

-1 6 ,8 5

-12*04

L

*

0,67125

0,6463

0,61

0 ,5 7

0 ,5 3

0 ,5

В

 

I

 

I

I

I

I

I

 

 

 

-.54

-

 

 

 

При использо­

 

 

 

вании ДРФ для

це­

 

 

 

пей защиты РЛС от

 

 

 

пассивных помех

 

 

 

необходимо перехо­

 

 

 

дить

к

адаптивным

 

 

 

методам

обработки.

 

 

 

Параметры L

и

 

 

 

AFn Т

, при

кото­

 

 

 

рых наблюдается

 

 

 

максимальное подав­

 

 

 

ление

помехи,

стро­

 

 

 

го

согласованы.

 

 

 

При этом отношение Мощностей остатков режекции

и v

на выхо-

дах каналов ДЭД практически равны, а их отношения

 

Таким образом, разработанная

структура

позволяет

предельно

редактировать пассивную помеху В широком диапазоне изменения норми­

рованной ширины спектра флюктуаций помехи

AFn T с

помощью управ­

ления единственного коэффициента

L

двухканального РФ, а коэффи­

циент

а « 0,36328125 является

постоянной

величиной-.

 

 

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

 

1.

Л .с. 1?628450. Устройство

компенсации пассивных помех / Д.И.

Попов,

В.А.Федоров,

Н.Д.Шестаков.

-

Опубл.

1 5 .1 0 ,7 8 .

- Бюл. №38.

2 .

Венцель Е.С.

Теория вероятности. -

М.: Фиэматгиз, 1962. -

766

с .

 

 

 

3. Тихонов В.И. Статистическая радиотехника. - М.: Сов. радио,

1966. - 677

с .

 

 

4 . Кузьмин С .З. Основы проектирования

цифровых систем обработ­

ки

сигналов.

- М.: Радио и связь, 1986. -

352 с .

- 55 -

УДК 6 2 1 .3 9 6 ,9 8

А.К.Бернюков, В,А,Егоров

МОДЕЛИРОВАНИЕ НА ЭВМ РАДЮНАВИГАЦИОННОГО ОРИШИКА С РАНГОВЫМ ОШЛ1ОДИТШМ

Процедуры принятия решений на борту времяимпульсных радиосис­

тем навигации

и посадки

самолетов

(ВИРСНГ1)

[ I ] - распознавание

помех, обнаружение-разрешение

сигнала Uc ( t )

и отражений от мест­

ных предметов

I/

({)

, оценивание навигационных параметров на

фоне шумовых помех

Ц

( i)

-

связаны с обработкой процесса

 

 

 

 

а

 

 

т= ис ш + £ % а ) + иш а ) .

Априорная неопределенность помеховых ситуаций, вызванная не­

известным набором местных предметов, изменением их свойств при об­ л ете,. влиянием на навигационные системы других электронных средств

в зоне управления воздушным движением, вызывает необходимость ис­

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

и объединенных в систему функционально-адаптивной обработки (ФАДО)

[ 2 ] .

 

 

 

 

 

 

( йВ<(1~2)

) лупа ан­

Специфика ВИРСНП -

сканирование узкого

тенны радиомаяка с

постоянной

скоростью Si

позволяет провести

квантование

рабочего сектора

<рс

на М=%/Дф микросекторов

по

углу Дцр л Si LН

соизмеримых с

шириной диаграммы направленности

антенны (ДНА), и организовать цифровую обработку реализаций (I )

в

квазиреальном времени

по

сегментам данных X- =? £x0f

 

,

где

f - f / j c

= 0 , 1 , 2 , . . .

-

номер периода сканирования.

В системе

мно-

роэтапноЦ ФАДО рре^смотрек этап

обнаружения сигнала

Uc

( t )

на

фо^ё^сигнадрцодрбнък отражений и шумовых негауссовых помех с нсиэвестнцм распределением и интенсивностью. Такую роль успешно выпол­ няет непараметричэркий обраруяште^ь, реализованный в виде модифи-

1{ированного

рангового

теста’ Вилко|сс6на [3 -4 ]

 

 

 

 

ГДе

ф1 -

элемент исследуемой

выборки X s

 

ху г - 9хл}

за, п

сканирований

антенны;

у,.

- элемент

скользящей

опорной помеховой

выборки У=

 

 

 

 

вэятой

из

т

смежных

с ана­

лизируемым дискретбв

по углу;

VJ -

сэлемент

ранговой

выборки

 

 

 

-

56 -

 

£ = [ г „ ...,

г*,

, г п | .

с _ порог.

 

Новостные

оценки непараметрического обнаружителя [3 - 4 J

не

отражают в полной мере специфики и тактики использования хМГСНИ: изменения азимутального направления на радиомаяк и соответственно перемещения сигнала внутри дискрета по углу* наличия в опорной выборке. вмелучевых сигнолоподобных отражений, разнообразия форм ДНА радиомаяка и схем построения бортовых приемников, а также парамет­

ров рангового обнаружителя гп, гг

Решение такой многофакторной задачи аналитическим методом

крайне затруднительно, методом полунатуркого моделирования нецеле­ сообразно, так как требуется сложный имитатор. В данном случае предложено решить поставленную задачу методом статистического мо­ делирования радионавигационного приемника с ранговым обнаружителем

на ЭВМ. Структура модели (рис. I ) состоит из блока I формирования сигналов Б1>С, отражений и шумов, преобразующей линейной и нелиней­ ной частей радионавигационного приемника (РП) 2 , рангового обнару­

жители

(1-Х)) 3 и

блока 4 статистической

обработки

реализаций (БСО).

 

В качестве

источника

 

 

 

р/.«

воздействий

восп-

 

 

 

Р«

/.-игея процесс

( I ) ,

РП

БСО

 

четвугощий радибнави-

 

 

 

 

 

.-донной обстановке, сос­

 

 

 

тоящей из радиомаяка, са­

 

 

 

молета

и нескольких мест­

 

 

 

ных предметов-отражателей. На выходе блока I формируются дискрет­ ные квадратурные составляющие, необходимые для определения комп­

лексных

огибающих сигнала

Uc (i)

, отражений 0 ц ( i)

и шумов

йш (i)

в

к -м дискрете

времени

J -го микросектора

(д<р*Ш%г)1-го

обзора

 

F Q . к - И - i c ),

 

 

 

 

 

 

ffJA J-

 

1 ^ ),

(3 )

где i i

-

шаг дискретизации в пределах дискрета Aty=Ktm c &i ,

F (j

K 'A i) - дискретные значения ДНА радиомаяка/ f

и t

временные положения сигнала и отражения соответственноf

°*

Предусмотрено моделирование четырех типовых ДНА:

 

Ш ч т (К ,8 )/(К А %(в)=!Ш*0(г 8)/0<г 8)а, ${*)*е*р[-(К3в)г] ,

^ (в )ф £ п [я е ^ 8 )] т (1 С А t * f A t Ki/ c o 6 » ( n в )],

(4)

 

 

 

 

 

 

 

- 57

-

 

 

 

 

 

 

где

Kj

2

ь 4

- константы,

характеризующие

ширину ДНА;

в

-

угол

азимутальной плоскости относительно направления на максимум;

60

^0 ,6 9 2 6 ;

б , » 0 ,3 2 1 9 ;

- -0 ,0 1 4 5

[ 5 ]

 

 

 

 

 

 

Независимые квадратурные составляющие

сигнала и отражения на

I -м обзоре формируются по правилам:

 

 

 

 

 

 

 

UCf f t ] =

Uc co s(2 JtZ [L l) ,

UC2 [L] - Uc s in (2 ftl t L) ) ,

 

 

 

 

% с п * р9 ц. со6<2*у9 [11),

ц л т - Р 9 ц * * л @ * $ т ) ,

 

 

а »

где

Z

CL]

и

Yy [ i ]

- независимые псевдослучайные числа

с

равно­

мерным законом распределения на интервале

(О Д ),

определяющие фнп

ктуации

фазового фронта сигнала.

 

 

 

 

 

 

 

 

Для формирования дискретных квадратурных составляющих шума

известными корреляционными функциями использован алгоритм

[6 ]

 

4lfy2

=

^0

 

1,2

 

 

 

 

 

(6)

 

 

'

 

U-L

 

 

 

 

 

 

 

 

 

где

 

i2

-

независимые

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

ных псевдослучайных

чисел

с

параметрами

(О Д ); £0 [ t i

-

весовые

коэффициенты,

определяемые, спектральной

плотностью шума;

 

L

 

параметр,

выбираемый

исходя из

точности формирования корреляцион­

ной функции. Для гауссовой корреляционной функции RQ[i]a]fSfte

где ft «(/STAt -ye /(T c ■

 

 

для корреляционной функции вида

 

 

 

 

 

 

 

где

2 x A t(p c /( T e

Ава ) .

 

 

 

 

Пару квадратурных последовательностей независимых нормальных

псевдослучайных чисел можно получить преобразованием пары последо­

вательностей независимых,

равномерных на интервале (О Д)

псевдо­

случайных чисел

Ж91 [ с , j

, к ]

:

 

z w J ll j , K]

=/-

i

СЦ , <7•"«(2 я Ъ9 Щ , * ] ) ,

 

2~ш2 ft>ji7} * )/-2 £пif

[l,j, к] -sift (2л2^ Cl,J, к"]) •

^

Программа позволяет

моделировать линейный и логарифмический

приемник, к выходу которого подключены последовательно детектор и интегратор. Модель преобразующей части 2 получается пут1 комби­ нации алгоритмов типовых элементов г дискретного комплексного фи­ льтра (УДЧ), блока взятия модуля и нелинейного'элемента (детек­ тора), дискретного фильтра (интегратора) со сбросом в конце диск­ рета по углу.

Дискретную огибающую суммарного напряжения на выходе УПЧ мо­ жно выразить через квадратурные составляющие (5) и (6 ):

 

 

 

 

 

 

-

58 -

 

 

 

 

 

 

 

 

 

 

U fA j.lO

~ eU{(Qli 4 ^ h Kh

° ^ l U0j f Г д .'О + Ч и Д ^ Л *’^ *

 

 

 

 

 

 

 

P

 

чг j

>

 

 

 

+ (^^ca [g>J+в I

«'<,<» сд.Ф % fц *J) f

 

 

 

 

 

 

Ja#

*

J

( 8 )

где

G*

- дисперсия шума,

 

0 - отношение сигнал/шум.

 

Без учета инерционности реального детектора напряжение на его

выходе

можно записать

в

виде

 

 

 

 

 

 

 

W|

[

 

 

 

{ 9)

где

т?

-

I для

линейного

и

т) *

2 для квадратурного детекторов.

 

При наличии интегратора со сбросом на выхода детектора сигнал

интегрируется в пределах дискрета по углу

 

 

 

 

 

иа [L.H = £ * Ч о С ч ,

< » »

 

 

 

 

 

 

tc*1

*

 

Смоделирован также пиковый детектор, вырабатывающий сигнал

Ugg [ LJ] ^^рДсетц*

Г •*J*

 

Моделирование рангового обнаружителя

сигнала на выходе радионавигационного приемника в блоке 3 произво­

дится в виде (2 ) . Последующая статистическая обработка в блоке 4

необходима для определения характеристик обнаружения для разных моделируемых ситуаций и оценки точности результатов. Оценка вероя­ тности обнаружения рр критерию Неймана-Пирсона и оценка распреде­

ления плотности реррятррсти находятся* как

 

$ =

2 * K f s ) ,

 

W /S) = P iS^ -A S*s< s j

/ A S ,

( I I )

 

*

~4S 4 S <

 

 

 

 

'

 

ГДе

d\

 

 

- оценка вероятности

попадания

случайных значения ранговой

статистики

S в Д -й интервал,

Nmax

“ максимальное

число

моделируемых статистик,

 

 

максимальное значение

статистики

-гр

интервала

.

 

Адекватность статистической

модели •реальным процессам

qfjejiQ-

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

оценками

аналитического

решения из

' [ 4 J

дДя случая обнаружения

только на фоне шума при

m

« 8 ,

гг « 4

и Q 0 , J , 2 . Точность

моделирования оказалась не хуже 5 %.

 

В результату моделирования строились гистогдомш зависимое-

тей^оценок вероятностей

правильного

Р

и ложного обнаружения

( сС * D

при

Q S Q) оф

порога

С

(рис. 2> а), по которым полу­

чены зависимости

P - J (Q )

при фиксированных уровнях ложной

тревоги

оС

(рис.

2 ,6 ) .

В процессе исследования на модели полу­

чены следующие результаты.

 

 

 

 

- 59 -

Рис. 2

1 . Взаимное пространственное перемещение самолета, радиомаяка

и местных предметов вызывает изменение

соотношения

сигнал/шум

и

сдвиг сигнала внутри дискрета по углу..

 

 

 

2 . Виды амплитудных характеристик

приемника (линейная, квад­

ратичная)

не влияют на характеристики обнаружзния.

 

 

3 . Усечение объема выборки ( т- п

) рангового обнаружителя

 

оказывает

существенное

влияние

при Q

2 -3 , причем сильнее ска­

 

зывается

уменьшение п

(число

накоплений рангов

г'• ) , чем ft*

 

(число вычислений ранга

) .

 

 

 

 

4 . Нестационарность опорной шумовой выборки, вызванная нали­ чием отраженных сигналов, вызывает некоторое снижение эффективно­ сти обнаружения и по своему действию э!Гв|Гвалёнтна уменьшению соот­ ношения сигнал/шум.

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

Оптимизация выбора параметров"рангового обнаружителя ( м п ) должна быть связана с характером аэропортов и тактикой их исполь­ зования.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1 . Сосновский А .А ., Хаймович И.А. Авиационная радионавигация; Справочник. - iM.: Транспорт, - 256 с .

2 . Бернюков А.К. Микропроцессорная фушщионалъно-адаптипиая обработка сигналов систем радионавигации в бортовой подсистем */