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

книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации

.pdf
Скачиваний:
47
Добавлен:
25.10.2023
Размер:
14.41 Mб
Скачать

Лд-Л/2

Л0 Л и Ао+А/2

Я

а.)

Л(*)

БВП

Рис. 3.26. Двухпороговый иепараметрическиГі дискриминатор.

Остальные

блоки выполняют операции в

соответствии

с уравнением (3.116).

 

 

Оценка (3.116) есть нелинейная функция от случай­

ных величин ku k2, k3.

Для приближенного

определения

среднего

и дисперсии

заменим нелинейную функцию

(3.116) ее разложением в двумерный ряд Тейлора в окре­

стности

средних значений случайных

величии

2 и k\ =

= ki + k2

[128]. Аналогичным

образом

поступим

и с квад­

ратом функции (3.116). По

найденным

математическим

ожиданиям двух разложений

определим

дисперсию

° а [ * І = л , [ ( * ) ] - « ? W -

Опуская громоздкие выкладки, записываем окончатель­ ный результат:

 

т, [Я]

і

л

1 рх

— рг

(3.117)

 

 

2

PZ —

PT '

 

 

1

 

а» [Я]

— А ) 3

•РХ—РІ

+ 2 р , (2Р2

1)], (3.118)

4 « ( р 2

 

 

 

 

 

где

Рі=р(Л.<(Ао—Д/2)]; Р2=р(Л,<'(Хо + Д/2)].

290

Всвою очередь, потенциальная точность оценки

(3.116), найденная по описанной выше методике, равна

0 м Л

. =

( Р 2

w X,

(3.119)

1-

lmni

/;

- А ) 2

 

Относительная эффективность оценки (3.116), опреде­ ляемая по формулам (3,37), (3.118), (3.119) для нор­ мального закона распределения параметра X при Яо = Яп, равна

 

• Л

 

А

, = х„)

 

w у Я

0 -

е(Х):

(3.120)

 

 

 

Аналогично, относительная эффективность (3.38) оценки (3.116) по сравнению с оценкой по среднему арифметическому в точке Хо—Хи будет равна

^ w - 2 G 2 ( / 7 ; , 7 / , ) 2 -

(3-121)

Зависимость e(X)=f(A), т. е. от расстояния между по­ рогами Д при точной настройке дискриминатора (Хо=Хя), приведена в табл. 3.1.

 

 

 

 

 

 

 

Т А Б Л И Ц А

3.1

Д/о

0

0,2

0,5

1

2

3

4

5

6

е{\) 0,637 0,685 0,746 0,804 0,740 0,491 0,256 0,111 0,028

Т А Б Л И Ц А 3 . 2

Д/о

0,2

0,5

1

2

3

4

5

6

 

0,043

0,119

0,262

0,496

0,655

0,54

0,269

0,13

—1

0,005

0,013

0,045

0,172

0,386

0,675

1,025

В соответствии с табл. 3.1 относительная эффектив­ ность в значительной степени зависит от отношения А/а. Для поддержания заданной эффективности необходимо одновременно с оценкой среднего значения получать

I Q*

90.1

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

Заметим, что высокая относительная эффективность оценки (3.116), определяемая по (3.21), возможна при Д/о>1 лишь в точке Ло=Ли , т. е. с вероятностью р = 0.

Перейдем теперь к задаче синтеза и анализа эффек­ тивности цифрового алгоритма оценки дисперсии. Эта задача имеет также и самостоятельное значение, напри­

мер, при построении

цифровых

систем

автоматической

регулировки усиления.

 

 

 

 

 

Возьмем вначале среднее арифметическое двух оце­

нок (3.107) с порогами в точках

Ко+А/2

и Ко—А/2,

полу­

чим алгоритм дискриминатора с двумя

порогами:

 

 

я = я 0

-

i - y - v / i

t

(

3 1 2 2 )

 

 

 

max

 

 

 

где k'i — число

превышений порога

Ко—А/2; кг — число

превышений порога Я.о + А/2.

 

 

 

 

Сравнивая

(3.122)

и

(3.116),

замечаем, что оценкой

для ш т а х может служить

выражение

 

 

 

 

w = = ^ - = Ь = - % - .

 

(3.123)

 

max

/;Д

„Д

 

\

'

Величина Wman для любого симметричного закона рас­ пределения связана с а зависимостью wmax—\/aa. Следо­ вательно, двухпороговое устройство позволяет дать па­ раметрическую оценку для ст в виде

л

лД

пА

ігм\

Приведем значения а

для

различных распределений: а —

= У 2-а. — нормальный

закон, а = |/"6— закон

Симпсона,

а = у\2 — равномерный закон, a = j/2—закон Лапласа.

Оценка (3.124), так же как и оценка (3.116), нели­ нейно зависит от случайных величин k'i и k2. Пользуясь методикой, о которой упоминалось выше при выводе со­ отношений (3.117) и (3.118), определяем математическое ожидание и дисперсию оценки (3.124):

я » і Й ~ — г - ^

г>

 

(3-125)

1 1 1

«(л> — Pi)

 

К

'

(3.126)

292

Для нахождения потенциальной точности при совместной

оценке % и о необходимо найти элементы информацион­ ной матрицы Фишера (3.28а). Можно показать, что информационная матрица имеет вид

;?/ог

О

О

2/г/ог

Тої да потенциальная точность оценки

(3.124) определит­

ся из выражения

 

 

 

 

 

 

 

д

 

1 J min

2на ( р 2

РіУ1

(3.127)

 

Относительная

эффективность

(3.38)

оценки (3.124)

с учетом (3.126)

и (3.127) равна

 

 

 

Д2

w (х0 ~Y Х0 = ХНj j

 

 

4 А ІРІ

— РІ)

(3.128)

 

 

 

Относительная эффективность оценки (3.124) по сравне­ нию с оценкой

для которой асимптотическое значение дисперсии состав­ ляет а2/2п, равна

(3.129)

2 Д » Л

Зависимость относительной эффективности е(а) от рас­ стояния между порогами для нормального распределе­ ния при А.о=А,л приведена в табл. 3.2.

Отметим, что относительная эффективность, опреде­

ляемая по (3.128)

и (3.129), характеризует лишь ско­

рость

сходимости

(отношение

дисперсий), поскольку

оценка

(3.124)

всегда

смещенная.

Смещение за­

висит

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

дискриминатора

относи­

тельно

медианы функции

распределения

(А,н—Аю). Для

определения этой зависимости заменим функцию

J F I ( A A I I )

ее разложением в ряд Бесселя в точке К,

ограничившись

2 9 3

первыми

тремя членами [13]:

 

 

 

 

 

 

 

Л ( л Л и ) ~0,5 + С[ (г—го)

+ (Сз/6) (г—го) X

 

 

 

 

 

 

Х [ ( г - г 0 ) 2 - 1 ] ,

 

 

 

(3.130)

где /•= (л—Ко)/А, fo=

п—Ко)/А

При К=,Ко—Л/2

и

Я =

= Яо + Д/2 из

(3.130)

получим

значения

pi

и

 

 

 

 

p i « 0,5— (0,5 + го) (сі—Сз/6)

 

 

 

 

 

 

(0,5 + г 0 ) 3 ( с з /6),

 

 

(3.131)

 

 

/ ; 2 « 0 , 5 + (0,5—г0 ) (СІсз/6)

+

 

 

 

 

 

 

+ ( 0 , 5 - г 0 ) 3 ( с з / 6 ) .

 

 

(3.132)

С учетом (3.125), (3.131), (3.132) получим нормирован­

ную величину смещения оценки (3.125):

 

 

 

 

 

щ [ а / а ] -

1 « А

 

 

 

(3-133)

 

 

 

 

0

а с 3 (0,25 +

3TQ Y )

 

 

 

где у2=

1 — 6( С і / С з ) .

 

 

 

 

 

 

 

 

Для

нормального

закона

при

Д/а^ 1

справедливы

следующие соотношения:

 

 

 

 

 

 

 

 

с 3 ~ — 0,334 (А/о) 3 ,

7 2 ^ 7 , 1 5 ( а / Д ) 2 .

 

 

Перепишем выражение (3.133) в виде

 

 

 

 

 

щ

-

1 »

 

^

 

 

5 -

1.

(3.134)

 

1

1 1

 

 

7,15с 2 г — 0,25 — 3rQ

 

 

'

Из формулы

(3.134)

следует, что минимум

смещения ле­

жит в окрестности точки Г о =0 . В табл. 3.2 приведена за­ висимость смещения от расстояния между порогами при

Хо = '%ц.

Определим теперь дискриминационные и флуктуационные характеристики дискриминаторов (3.107) и (3.116). Зависимость математического ожидания оценки параметра от расстройки дискриминатора (Яи—Ко) назы­ вается дискриминационной характеристикой. При извест­ ной функции распределения А(лДн) дискриминацион­ ные характеристики дискриминаторов (3.107) и (3.116) строятся по выражениям (3.109) и (3.117) соответствен­ но. Сравнение дискриминационных характеристик удобно

производить, используя

разложение функции

Fi(K/Ku)

в ряд (3.130). Запишем

аналитические выражения

дис-

294

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

нормированной

оценки /'= (А,—Яо)/А:

 

— дискриминатор Тейлора (3.107)

 

' « . М ^ . - ' о / П

(3-135)

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

Бесселя (3.116)

±

 

* 2 — 0,25 — Зг*

Названия дискриминаторов соответствуют названию ряда, который использовался при .синтезе алгоритма оценки. Параметры г0 и у были введены выше.

При оценке неслучайного параметра функцию распре­ деления Fi(X/Ku) можно задать следующим образом:

 

 

{ 1 при ЯЗэА,,,

 

что соответствует отношению

о/Л1—ИХ

 

Тогда дискриминационная

характеристика

вырож­

дается в зависимость:

 

 

 

 

 

( — 1 / 2 ш Д при

го<ґ0,

 

— дискриминатор Тейлора,

 

 

 

і

оо при /-„< 0,5,

 

» я , Й = (

0

при 0 , 5 < г 0

< 0 , 5 ,

(3-138)

[ - f -

оо

при г 0

^ 0 , 5

 

 

яепараметрический дискриминатор Бесселя.

В другом предельном случае, когда о/Д>-оо, у2 >-°о и дискриминационная характеристика вырож­ дается в прямую линию mi[r]=ro.

Дискриминационные характеристики при произволь­ ном отношении о/Д будут располагаться между граница­ ми ГПІ{Г] = Г0 и границами, определяемыми соотношениями (3.137) и (3.138).

2?5

На рис. 3.27 приведены дискриминационные

характе­

ристики

непараметрического

дискриминатора

Бесселя,

построенные

.по

выражениям

(3 . 117) — при

0/Д = О,6,

(3.136) —при

оШ=°°,

(3 . 138) — при о/Л = 0. При

пост­

роении предполагалось, что Fi{%l%a)

соответствует

нор­

мальному

закону.

 

 

 

 

 

 

Вследствие ограниченности разрядной сетки цифро­

вых функциональных

узлов, входящих в состав дискри-

 

 

 

 

6Vn

 

 

 

 

 

Флуктуаци-

*,г- -0,6

 

 

 

 

 

\онная

харах-

 

 

у

-

 

 

\mepacmuxa

 

 

 

 

 

 

 

 

 

0,8- -0,4

 

 

 

 

 

 

 

 

 

 

6/6=0

 

 

 

 

 

-0,2

У'б/А

=0,6

 

 

 

-0,4

 

-0,2

- 4 * -

0,2

о,ч

Й

 

 

 

 

 

 

 

 

 

 

 

 

 

--0,2

 

 

 

 

-0,8-

j 6/А -оо

*

-1,2 + • -0,6

Рис. 3.27. Дискриминационная и флуктуационная характеристики иепараметрического дискриминатора Бесселя.

минатора, дискриминационные характеристики при 0,5<

< г < 0,5 будут располагаться параллельно оси г, что показано на рис. 3.27 штриховыми линиями.

Флуктуационная характеристика определяется как ко­ рень квадратный из дисперсии оценки .в функции от рас­ стройки дискриминатора и—Я,о). Дисперсии оценок за­ даны соотношениями (3.110) и (3.118). Флуктуационная характеристика непараметрического дискриминатора Бесселя (3.116) при а/А=0,6 приведена на рис. 3.27, ,

Г Л А В А Ч Е Т В Е Р Т А Я

Ц И Ф Р О В ЫЕ ФИЛЬТРЫ

4.1.В В Е Д Е Н И Е

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

Основные затруднения при технической реализации ЦФ с заданными характеристиками связаны с обеспече­ нием требуемого быстродействия при большом объеме памяти. В этой связи в радиолокации можно выделить три характерные группы ЦФ.

1. Согласованные и квазиоптимальные ЦФ, исполь­ зуемые в устройствах обнаружения целей и измерения их координат, т. е. при первичной обработке. Такие ЦФ реализируются в виде специализированных ЭЦВУ, функ­ ционирующих в реальном масштабе времени (практиче­ ски на частотах, измеряемых единицами и десятками ме­ гагерц). Требуемое число разрядов аналого-цифрового преобразования для выделения сигналов при некоррели­ рованных помехах невелико: при обработке некогерент­ ных сигналов т = \, когерентных — т ^ > 2 . Подобные ма­ лоразрядные нерекурсивные ЦФ, в сущности, являются обнаружителями или измерителями. При выделении же сигналов движущихся целей на фоне пассивных (корре­ лированных) помех число разрядов преобразования со­ ставляет обычно 810.

Высокая тактовая частота затрудняет реализацию таких ЦФ в виде программ ЭЦВМ .

297

2.ЦФ, осуществляющие полосовую и гребенчатую фильтрацию. Основное требование к таким ЦФ — обеспе­ чение заданных амплитудно-частотных и фазо-частотных характеристик. Число разрядов преобразования, как пра­

вило, велико и составляет 8—12, разрядность арифмети­ ческих операций 812. В зависимости от ширины спектра входного процесса и несущей частоты ЦФ этой группы реалнзпруются либо в виде специализированных ЭЦВУ, либо в виде программ ЭЦВМ . Особое место занимают ЦФ — анализаторы «нулей» процесса.

3.В задачах сглаживания и экстраполяции параме­ тров и координат траекторий, т. е. при вторичной (траекторной) обработке, а также в Р Л С сопровождения целей находят применение ЦФ, реализуемые как программным

способом в универсальных или специализированных Э Ц В М , так и в виде специализированных ЭЦВУ . С по­ мощью этой группы ЦФ решаются задачи сглаживания и экстраполяции координат и параметров траекторий це­ лей в реальном масштабе времени. Число разрядов пре­ образования составляет 8—12, в ЭВЦУ или Э Ц В М про­ изводится обработка информации, практически непрерыв­ ной по уровню (с точностью до младшего разряда ман­ тиссы чисел, представимых в Э Ц В М или Э Ц В У ) .

Первая группа ЦФ рассматривалась во 2 и 3 главах, а вторая и третья группы будут рассматриваться в дан­ ной главе.

4.2.П Р И К Л А Д Н Ы Е ЗАДАЧИ Т Е О Р И И

ЛИ Н Е Й Н Ы Х Д И С К Р Е Т Н Ы Х СИСТЕМ

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

Т

I

где u(l), v(t) —соответственно входной и выходной про­ цессы; a,, hi — вещественные коэффициенты,

298

Передаточная функция системы, описываемой уравне­ нием (4.1), является дроб.но-рациональной:

у / 0

1 — gmS"1

+

flm-,Sm-'

+ .. . + a , S +

fl,

_A(S)

 

\^)

6,S« +

6 l _ 1 S i - '

+ ... + u I S + 6 e

 

~B(S)'

 

где S =

a-|-J(o;

i =

Y—1;

a

— вещественная

переменная;

со — ч а с т о т а . У

физически

 

реализуемой

системы

т<^1.

Техническая

реализация

уравнения (4.2) — т а к

назы­

ваемое прямое программирование '[8] — связана с исполь­ зованием интеграторов, усилителей и сумматоров.

Частотная

характеристика

системы получается из

(4.2)

подстановкой

5 = « о .

Система

будет

устойчивой,

если все корни

знаменателя

B(S)—полюсы

передаточ­

ной

функции K(S)

—-лежат

в

левой

части

плоскости S,

т. е. имеют отрицательные вещественные части. Уравне­ ние B(S)=0 называется характеристическим уравнени­ ем системы. Импульсная переходная характеристика си­

стемы

определяется

как обратное

преобразование

Лап­

ласа передаточной функции (4.2):

 

 

 

 

 

h{t) = -%a§K{S)esidS.

 

(4.2а)

Полиномы

с

вещественными коэффициентами Л (5) и

B(S)

имеют только

действительные

или комплексные по­

парно

сопряженные

корни — нули S° и полюсы Sp

. При

разложении этих

полиномов

на множители

 

 

4(S)

=

( S - S ? ) (S-S°2)...

( S - S ° J ,

 

 

B(S)

=

{S — Spl)

( S - S p ) . . . { S - S p )

 

каждый действительный корень даст линейный двучлен относительно 5, а каждая пара комплексно-сопряженных корней — квадратный трехчлен с вещественными коэффи­ циентами относительно 5. Каждый 'Нулевой корень даст множитель S. С учетом вышесказанного, передаточная функция (4.2) может быть представлена в виде следую­ щего разложения [24,-135]:

/к

Л (S - ос,) П (S - a s - t Y » ) (S - ah + < t h)

K{S) — C -E±

Щ

-, (4.3)

П

(5 - Pp) П (5 - P* -

^) (S - P*+іхл •

299

Соседние файлы в папке книги из ГПНТБ