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

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

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

Результаты решения уравнения (3.69)

представлены

на рис. 3.18 (кривые в, г) . Они достаточно

близки к ха­

рактеристикам, полученным в соответствии с (3.68). Та­

ким

образом,

при 7 2 ^ 1 , PN opi=0,203,

с

увеличением

Ф /Jjvopt

уменьшается.

Как показывают

расчеты на

Э Ц В М ,

величина

а[(3]

слабо

зависит

от

UQ. Поэтому

с

практической

точки

зрения

целесообразно

выбрать

Рк = 0,1.

 

 

 

 

 

 

 

 

Обобщенные

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

потенциальной

точности

измерения азимута при бинарном квантовании

приведе­

ны на рис. 3.19. Они вычислены

по формуле

(3.67) (кри­

вая а соответствует нефлуктуирующей ЭОП цели, кри­

вая б — шумоподобным

флуктуациям ЭОП) . Расчет кри-

 

 

 

 

0,5

 

 

0,203

 

 

 

 

 

 

0,f

 

 

 

10

 

 

0,01

 

V

 

 

 

 

 

г'

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\

 

 

 

 

 

 

\

 

 

 

0,001

to

£г;т

o,i

Ю

f,dS

-го

-20 -10

Рис. 3.18.

Зависимости

рКат>\.=

Рис. 3.19. Обобщенные характери-

 

 

 

 

стики потенциальной точности

из­

 

 

 

 

мерения

азимута.

 

вой б проводился по формуле, приведенной в- работе [118]. Как показывают расчеты, проведенные методом Монте-Карло на ЭЦВМ, погрешность характеристик, представленных на рис. 3.19, не превышает 10%.

3.4.3. Синтез цифровых измерителей и анализ эффек­ тивности при зависимых отсчетах [174]. В п. 2.6.3 обсуж­ дались особенности статистического синтеза цифровых обнаружителей при зависимых отсчетах дискретных пе­ ременных на выходе квантизатора. Здесь мы рассмотрим статистический синтез цифровых измерителей: однокаиальных, многоканальных и дискриминаторных. Во всех трех случаях в рамка» метода максимума правдоподо­ бия достаточной статистикой является логарифм функ260

цни правдоподобия или в соответствии с (3.9) условный коэффициент правдоподобия. Для простоты изложения будем считать, что связность цепи Маркова v = l . Обоб­ щения на случай \>>1 не представляет затруднений. Не­ обходимость рассмотрения статистически связанных от­ счетов выборочных переменных обсуждалась в п. 2.6.3.

 

При

построении

многоканальных неследящих

измерите­

лей

необходимо

решить

систему

 

(3.1)

или

определить

max In Л (и | Я),

т. е. в

каждом

канале

должна

формиро-

—•

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

ваться

статистика

[см. (2.204)]

 

 

 

 

 

 

In Л (и | * , ) = t

daWa

(0)з + £

S

da?

(k) Wa?

(k)it (3.70)

 

 

 

ec = l

 

ft=l

a,

3

 

 

 

где

/ =

1, M; k =

1, n — 1;

 

 

 

 

 

W

 

= Ш

 

 

; W ( 0 ) j = = l n t o ( ° , X

i )

• 73.71)

В качестве оценки Я скалярного параметра Я выбирается

номер

канала

х,

в котором In Л (и\ Ях ) = та х {1пЛ(и| Я)}.

 

 

 

 

 

х

 

 

Структурная

схема

измерителя (3.70)

приведена

на

рис. 3.20. Пергходныз вероятности paV

входящиэ

в

выражения

для весовых коэффициентов

(3.71), опреде­

ляются по методике п. 1.7.2. При произвольной

связности

цепи Маркова v > l число входов в блоки

^ ^ ( ^ Ь равно

v + 1 ;

каждый

вход

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

т. Схема

(рис. 3.20) является универсальной для оценки параме­

тров широкото класса сигналов

и помех. Дополнив ре-

 

 

 

-+

шающее устройство (РУ) блоком сравнения

1пЛ(ы т ах)

с поротом с, получим систему

совместного

обнаруже­

н и я — измерения.

 

 

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

метры

нескольких целей, то следует изменить блок РУ.

В этом

случае необходимо сравнивать с порогом с вели-

—>

чину \nA(u\Xj), /=-1, М, в каждом канале. В качестве оценок параметров решающим устройством должны вы-

261

uft)

л

ufk)

 

 

 

 

 

 

 

 

 

 

 

Рг.с.

 

 

 

 

 

 

п-1 n-г

і

о

 

 

 

 

• в

 

 

 

 

 

 

 

• • и /

 

/V) I

 

 

 

 

*

Л=<*>

Wc,(0)M

-ігЧҐ1 -1

РИС: 3.20. Структурная схема измерителя ( 3 . 7 0 ) .

даваться номера каналов %, у, р, ... , в которых произош­ ло превышение порога обнаружения. Одноканальный неследящий измеритель времени запаздывания имеет при v = l вид, соответствующий рис. 3.17 (с учетом штрихо­ вых линий).

Если последовательность отсчетов и(1г) на выходе квантизатора образует однородную односвязную цепь Маркова, то в соответствии с (2.205) в многоканальном измерителе должны формироваться статистики

-> ->

In Л* (и | h)

=

In Af t . , І Я,-) + S da9 (k) Wa9J,

(3.72)

 

 

 

 

a. p

 

 

где In Л0 (її I Xj) =

2

dWa

(0), а коэффициенты Waj,

Wa9j

 

 

a

 

 

 

 

 

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

из (3.71). Входящие в

(3.71)

переходные

вероятности

pe pS A r(*,Aj) =

> e p s A r(Aj) . (a .

Р =

1, г) не

зави­

сят от k, поскольку процесс u(t) на входе квантизатора стационарный.

2 6 2

Структурная схема

алгоритма

(3.72) изображена на

рис. 3.21. Для упрощения

реализации из схемы исключе­

ны начальные веса Waj,

j

= [,M,

что можно сделать при

большом объеме выборки. Схема содержит М накопите­ лей _27з, каждый из которых может быть выполнен в двух видах. Для РЛС с электронным сканированием антенны и последовательным анализом Вальда накопи­ тель £ есть просто дискретный интегратор, т. е. замкну-

uft)

П

й(к)

Рг.с

 

 

 

 

1 О

РУ

Х = д€

 

 

тах%

 

•Рис. 3.21. Структурная схема алгоритма (3.72).

гая система, состоящая из сумматора и элемента за­ держки накопленной информации на один такт:

 

y(k)=y{k—l)+x(k),

 

где x(k)—входная

переменная; y(k)—выходная

пере­

менная. Подобный накопитель можно выполнять на осно­

ве накапливающего

сумматора. Для РЛ С с механиче­

ским сканированием

антенны при использовании крите-

рия Неймана — Пирсона накопитель

^ есть дискретное

апериодическое звено, т. е. оператор

экспоненциального

сглаживания

 

 

y(k) = (l-a)y(k-l)+ax(k)

=y(k-l)

+

+ a[x(k)—y(k—\)],

 

(3.73)

где а = 2 / ( л + 1 ) .

За основу построения такого накопителя можно так­ же взять накапливающий сумматор. Решающее устрой­ ство (РУ) так же, как и в схеме рис. 3.20, может быть выполнено в двух видах.

263

При увеличении связности цепи Маркова боле І і воз­ растает число входов в блоки W B ? / - Разрядность регист­ ра адреса постоянного запоминающего устройства, храня­

щего веса Wa?, равна m{v-\-l),

где т — разрядность пе

ременных на выходе квантизатора. Объем оборудования схемы рис. 3.21 примерно в п раз меньше, чем у схемы рис. 3.20.

Рассмотрим теперь построение оптимальных цифро­ вых дискриминаторов. В соответствии с уравнением оптимального дискриминатора (3.8), формулой (3.9) и алгоритмом (2.204) имеем выражение для оценки

 

 

 

 

п - 1

 

 

 

 

 

 

 

 

 

 

 

£

2

d«p (ft) У„р (ft)

 

 

 

 

 

 

 

ft=l aS

 

 

 

 

 

(3.74)

 

 

 

0

и—1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

S rf«p

(*)

Z„p (ft)

 

 

 

 

 

 

 

fc=l

чР

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

v a

(k)—w

.(A) =

І - m

 

= ^ Ц ^ ;

 

(3.75)

"Pv

;

"p V

'

d\

 

i°ap(ft)w

P„p(ft)sw

 

v

 

 

ZAK)—

 

 

pa?(k)

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

по

мето

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

дике п.

 

1.7.2.

Начальные

значения

2 d a ( 0 ) V a ( 0 ) и

 

 

 

 

 

 

 

 

 

a

 

 

2 ^ „ ( 0 ) г а ( 0 ) в

выражение

(3.74)

не включены,

что

спра-

a

 

 

 

 

 

 

 

 

 

 

 

вгдливо

в

установившемся

режиме.

 

 

 

При обработке сигналов с неизвестным временем при­ хода и усреднении по k на интервале постоянства пара­ метра X, суммы в числителе и знаменателе (3.74) следует заменить на операторы сглаживания: скользящего сред­ него (движущееся окно) или экспоненциального сглажи­ вания (3.73). Последний оператор предпочтительнее, так как он значительно проще в реализации, а характери­ стики накопления для обоих операторов практически одинаковы. Еще одно упрощение возможно, если неизме-

ряемые параметры a(t)—неэнергетические,

а энергети-

54

ческий параметр — амплитуда — нормируется с помощью схемы АРУ или логарифмического УПЧ. Тогда, как ука­ зывалось в п. 3.2.1, знаменатель в (3.74), т. е.

д21п Р(и\Х)/дХ20, может быть вычислен заранее и введен .в оптимальный дискриминатор. Теперь уравнение дискриминатора принимает вид

X(k)

= X0-e(k)

= l 0 ~

{e(k-

1) +

 

+

a[V{u(k-

 

\),Z(k))-7(k-

1)]},

(3.76)

где

 

 

 

 

 

 

V Си (k -

1), u(k)) = 2

(k) V a p (k)

=

 

= 2 ^ p K P W ] .

 

(3-77)

Выражение для сигнала рассогласования є{k) на /г-м шаге соответствует рекурсивному алгоритму оператора экспоненциального сглаживания (3.73).

На основ? уравнения (3.76) в принципе можно постро­ ить цифровой следящий измеритель, формируя оценку из­ меряемого параметра X{k) с помощью сглаживающих це­ пей. Однако при этом возникает серьезноэ затруднение, связанное с" большой разрядностью ПЗУ, хранящего веса

Va? [u(k— 1), u(k)\. Дело в том, что последовательность дискретных переменных u(k) на выходе квантизатора, принимающих значения х , а = 1 , г, образует неоднородную

цепь Маркова, так как

каждая из переменных и (k)

зави­

сит от параметра X(k),

изменяющегося в большом

дина­

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

нального

неследящего

измерителя, изображенногона

рис. 3.20.

 

 

 

 

 

 

 

Для преодоления этого затруднения произведем

центрирование входного

процесса

u(t)

(3.10). Пусть нам

известна

оценка входного

сигнала

S(t,

K{t)^S(t,

l ( r ) ) ,

а сам входной сигнал S(t,

X(t))

нормирован. Оценка па­

раметра

 

 

 

 

 

 

 

 

X{t)^=XK(t) +

 

bX =

XK{t)-\-X-\-b%

(3.78)

265

состоитиз трех компонент: истинное значение параметра

%n(t);

флуктуационная компонента 6Я с малой дисперси­

ей и нулевым

средним;

смещение Я, вызванное

динами­

ческой ошибкой, погрешностью нормировки, аппаратур­

ной погрешностью.

 

 

 

 

 

 

Будем считать, что ДЯ мало, тогда

 

 

 

s

(t, я (/))= 5 (/, яи (0) +

и d S (

t ' i l t )

)

L + • • •

После

центрирования

 

 

 

 

 

 

v (t) =

u(t)-S

(t, £(t))

її (і) -

ДЯ й

? ( ^ ( 0

)

х = х

. (3.80)

При малых отклонениях ДЯ процесс v(t) можно считать стационарным, а получаемую при квантовании цепь Маркова — однородной.

Вычислим теперь весовую функцию (3.75), (3.77). Пе­ реходя к дискретному времени, с учетом (1.90а) получаем

 

 

 

 

 

 

 

 

|

|

w2

(у, ,о2 ) dVidvz

 

 

 

 

 

 

 

 

 

 

 

 

 

- ,

(3.81)

 

 

 

 

 

 

 

 

 

ton.

 

 

 

где y a

( k - l )

=

xa(k

-

l ) - 5 , (Q =

n(k~

 

;

 

уй (k) =

хъ

(k) -

Sa (£,) =

л(й) -

AASS'|.

(3.82)

Дифференцируя

(3.81)

с

учетом

(3.82),

имеем

 

 

д 1 п

Р^Ь/ь'У?)

 

 

Л,а (* . - !) . У (*)) X

 

 

 

дХ

 

 

х=>.„

 

 

X

дуа

 

 

дХ

(5г/р

 

aXj

Х=х„

(3.83)

 

 

 

 

 

 

 

 

 

 

 

 

 

В свою

очередь,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а (дя) as,

a^s

 

 

as,

а\

х=х„.

 

 

~ Ж " Ж +

А я а х 2

х=х„

— с ая х=х п

 

 

 

 

а(дя)

as2

 

 

az s

 

 

 

(3.84)

 

і

 

 

 

 

 

 

as,

ая

 

 

~ а Г

Ж

~т ~'Л Я дх2

х=х„

— с ая; х=х„

266

поскольку отклонение ЛА по предположению мало, а про­ изводная д(АХ)/дХ=с-~ Л из-за инерционности цепей сглаживания.

Обозначив

Ра$ У* (й-

1

(3.85)

Ра$ (Уа\У$)

с учетом (3.83) и (3.84) получим выражение для весовой

функции дискриминатора

(3.77):

 

 

 

 

 

 

 

 

 

 

 

dS,

 

 

 

V [Уа(k-l),

цЩ

=cW,

а(k-l),

ур (k)\ ^ х=х, [ +

 

 

+

 

cWa[yJk

 

 

 

 

(3.86)

Структурная схема цифрового следящего измерителя,

формирующего

оценку X{k) изменяющегося

параметра

X(t),

приведена

на рис. 3.22. Подобный измеритель

явля­

ется

структурно-инвариантным

как «по отношению

к ко­

дировке параметра

X в сигнале

S(t, X(t)),

так и по от-

u{t)=s(t,*.(t))+nrt)

 

 

 

 

 

 

 

 

 

 

п

шРг

 

e

t

aг

 

s(t,m))

 

 

Wt

 

- Ш)

 

 

 

ПЗУ

 

 

 

 

 

 

7-1

 

7f

m

j

t

 

 

 

 

 

 

m сц

 

 

 

 

 

 

 

 

ГОН

 

 

 

 

 

 

 

 

6

г

5

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.22. Структурная схема цифрового следящего измерителя.

ношению к статистическим свойствам помехи n(t), кото­ рая должна быть центрирована. В состав цифрового следящего измерителя, построенного в соответствии с из­ ложенными принципами, входят: I) квантизатор, 2) ком­ бинационная часть ( П З У ) , формирующая весовые коэф­ фициенты Wu W2, 3) арифметические блоки, 4) усред­ няющее устройство дискриминатора (3.73), 5) цифровые

267

сглаживающие цепи (СЦ) и 6) генератор опорного •на­ пряжения ('ГОН).

Комбинационная часть 2, формирующая веса Wu W2, выполняет функции ослабления мешающего действия по­ мехи. Сглаживающие цепи (СЦ) представляют собой цифровой фильтр (см. § 4.5, 4.6), формирующий как сглаженную оценку параметра в текущий момет време­

ни X(k),

так и

экстраполированную

на (/г + 1)-й

мо­

мент

%*'(k).

 

 

 

 

Кроме

параметра

Я (/г) с выходов сглаживающих

цепей

можно

получить

его

производные Я',

Я". Блоки 2,

3 и 4

и составляют, собственно, дискриминатор, который явля­ ется универсальным и пригоден практически для выделе­ ния различных параметров радиолокационных сигналов на фоне помех. Конкретные величины весовых коэф­ фициентов WyW2 должны быть рассчитаны заранее и храниться в памяти ПЗУ. Иными словами, данная структура является параметрической. 'Возможность запи­ си .в ПЗУ любой информации и делает ее структурноинвариантной по отношению к помехе с любым законом распределения в рамках выбранной дискретной аппрок­ симации, необходимо только, чтобы помеха была центри­ рована. Не представляет труда модициф.ировать схему при связности цепи Маркова v > l .

Генератор опорного напряжения (ГОН) формирует опорный сигнал как функцию экстраполированного зна­ чения параметра А*(&), а также и производную этого сигнала по параметру К, т. е. dS/d%. Наиболее громозд­ кими блоками являются сглаживающие цепи и арифме­ тические блоки. Благодаря центрированию и нормировке

процесса v(t)

квантование может быть

малоразрядным:

т = 2ч-6.

 

 

 

 

В

п. 2.6.4

рассматривался

вопрос о

связи двух мето­

дов

синтеза

цифровых обнаружителей:

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

и по аналоговому прототипу. Все сказанное в п. 2.6.4

от­

носится и к методам синтеза цифровых

измерителей;

так,

при

уменьшении интервалов

квантования Д«а —>-0,

а =

= 1, г, гУ-СО оптимальные алгоритмы, синтезированные обоими методами, совпадают. Однако техническая реали­ зация подобного цифрового измерителя затруднительна.

В качестве примера оценки эффективности цифрового измерителя рассчитаем потенциальную точность измере­ ния углового положения цели при бинарном квантова-

268

нии, когда последовательность бинарных случайных ве­

личин на

выходе

квантизатора

образует

односвязиую

цепь Маркова [121].

 

 

 

 

 

 

 

Структурная схема

одноканального

измерителя

ази­

мута при статистически связанных отсчетах

u(k)

на

вы­

ходе квантизатора

приведена на

рис.

3.17.

Индикатор

Я — это датчик текущего азимута

(см. рис. 3.13). При би­

нарном

квантовании

процесса

u(t)

блоки

Va?

(k),

k — Q, її—1 могут быть

упрощены.

В самом деле, запишем

логарифм функции правдоподобия для односвязной не­

однородной цепи

Маркова

с

двумя

состояниями

(1.10):

 

1 п Р ( а | р ) =

£

[d0 0 (fe) I n р 0 0

r f 0 l

(fe)Inp0 1

(A)

 

+

 

 

+ rf10(fe)lnp10(/e)

+

rf„(/e)lnpI1(/e)], .

 

 

(3.87)

где

rfoo(0)=rfi0(0)=rfo(0),

 

rfoi(0)

= d u ( 0 ) = r f i ( 0 ) — индика­

тор

начального

состояния

(1.8);

p0 o(0) = /?ю(0)

 

=po(0),

Рої (0) = рн (0) = рі(0) — вероятность

начального

состоя­

ния. Упростим .выражение

(3.87), использовав

очевидные

соотношения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Poi(k)=l

 

—poo (k),

ри (k)

=

1 — pw

(k),

 

 

 

 

. doi{k) = \—d0o(k),

du(k)

=

 

l—dw(k).

 

 

 

 

Тогда

формула (3.87) принимает вид

 

 

 

 

 

 

 

 

 

л - І

 

 

 

 

 

 

 

 

 

 

 

 

 

 

In P (и

I p) =

2

K ,

(k)

[In pn (k) -

In (1 -

 

p01

(k)]

+

 

 

+

d,, (k)

[In pu

(k) -

In (1 -

plt

(k)))

 

+

 

 

 

 

+

l n ( l _ P o i

( A ) )

+

l n ( l і0(Щ-

 

 

 

 

 

(3.88)

 

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

с

(3.2)

имеем

уравнение

максимума

функции

правдоподобия

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

л - 1

 

 

 

 

 

 

 

 

 

 

 

 

d l n P H P )

 

 

(d

(

h )

ffn

(k)

rf

 

„ .

w

 

<?p

 

 

2 j \ a ° , W

Л , ( А ) л . ( Л ) +

а

, ,

1 *

;

Л

 

 

 

 

 

fc=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л—1

 

 

 

 

 

 

 

 

 

 

 

 

p ' e p (k) =

 

 

fe=0

 

 

 

 

 

 

 

 

 

 

 

где

(?pep

(pft)/(?p;

a,

p ^

0, 1.

 

 

 

 

 

 

 

269

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