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

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

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

п—I

Р\г (k)

 

0 (k)

 

 

 

Сумма

 

входящая

в

— Рог (k)

— Pi о (k)

k=0

 

 

 

 

 

 

 

 

 

(3.89), не зависит

от выходных

переменных

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

da p(k) = 0, 1 и в

силу симметрии диаграммы направлен­

ности антенны (2.41) равна

нулю.

 

 

 

В соответствии' с (3.89) измеритель азимута должен

выделять в k-й и (k-—1)-й

позициях лишь

две комбина­

ции из четырех, т. е. 01 и 11, придавая каждой

из

них

соответствующий вес: V0i{k)

или Vn{k).

Тем

самым

упрощаются блоки Va^(k),

т,

е. аппаратурные

затраты

сокращаются вдвое.

 

 

 

 

 

Для оценки точности измерения азимута .воспользуем­ ся формулой Крамера — Рао (3.28), имея в виду ограни­ чения в связи с ее 'применением, отмеченные в п. 3.4.2.

Перепишем

(3.28) в эквивалентной

форме

 

а 2

[Я]

 

. = •

 

1

 

 

 

 

 

 

 

••{

д In wn

| X)

 

 

 

 

 

 

Ж

 

 

 

 

 

 

 

1

 

 

 

(3.90)

 

 

 

 

 

 

 

 

" M X ) '

 

 

 

 

д2 1пм>і(и|Х)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дХг

 

 

 

 

где

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

а Г

 

 

 

/ п ( Я )

= да,

д In wn

| X)

(3.90а)

 

 

 

 

 

 

— информационная

мера

Фишера

(информация, содер-

жащаяся в выборке и относительно

параметра X).

В соответствии с (3.90) дважды дифференцируем

1пЯ(й|Р) (3.87)

 

по р:

 

 

 

 

 

д* In я |р)

 

п—I

 

р"оо

{k) Poo {k) •

 

•Si

 

 

 

 

 

 

d00(k)

 

 

P\o (k)

 

 

 

 

 

ft=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Р"ог

(А) Рог (k) •

 

 

 

 

 

 

 

Рої

(*)

 

 

 

 

+

 

 

Р"го (k) P l t

(k) -

р>ї0 (k)

 

 

 

d„(k)

 

 

 

 

 

270

гдо

da4 (к) -- индикатор

перекода

(1-11);

р " ^ (/«)==

 

Производя статистическое усреднение (3.91) в соот­

ветствии с (3.90), получаем

 

 

 

 

 

 

 

 

-•>

 

 

/ я (Р) = -

/ » ,

[ -

dp2

 

 

 

 

 

 

 

 

 

и—I

 

 

 

 

 

 

 

 

 

/ Л х ? ( ^

(3.92)

 

 

 

 

 

(/г)

 

 

 

 

 

 

 

Й=0 а ,

р

 

 

 

 

где

р ( а ? ) абсолютная

вероятность

состояний

и к _ 1 = х а ,

и-к — х^ (1.90). Преобразуя (3.92) с учетом (1.90а), а так­

же принимая во внимание, что

Роо (/ г ) =

1 — Рої

 

Р'оо (k) =

Р'ОІ ( / Г ) .

 

 

Р ' о о ( ^ ) = Р ш (*)•

 

 

 

получаем

 

 

 

 

 

 

 

 

 

^ [ Й и , „ =

1 Д п ( Р ) =

'

 

 

Г«—1

1

 

 

 

 

 

 

 

= £ Е1" <?Р J

Pal (*)

/>.оИ*)

 

(3.93)

 

 

ft=0

а = 0

 

 

 

 

 

 

 

где р а (fe) — абсолютная

[вероятность

(1.89).

Легко

убе­

диться, что для

нульсвязной

цепи

Маркова

выражение

(3.93) соответствует формуле (3.64).

 

 

 

 

Соотношение

(3.93)

можно

обобщить на случаи

мно­

госвязной цепи

Маркова,

многоуровневого

квантования

при измерении любого

скалярного параметра X. Так, на­

пример, для v = 2 и произвольного

числа уровней

кван­

тования, равного Ги потенциальная точность измерения параметра X равна

•"Й-НЕЕ^^-Г-

(3'94)

4=0 apf

 

где а, р, у = 1, г, г = г, ——{ 1, р а р т (k) = у

271

J • .

При совместной

оценке компонент вектора параме-

->

 

тров Я=(А.і, %2, ... ,

h)T элементы информационной ма­

трицы Фишера (3.28а) для двухсвязной цепи Маркова соответствуют выражению

&=0

а, 3, т

г д е р в ? т ( / г ) = р а 3 7

І*(*)]; І., і = " Г Т ; а - Р- т = ^

Совместная оценка компонент вектора Я производится многомерными измерителями [69].

3.5.Ц И Ф Р О В Ы Е И З М Е Р И Т Е Л И ,

С И Н Т Е З И Р О В А Н Н Ы Е Э В Р И С Т И Ч Е С К И М МЕТОДОМ

3.5.1: Цифровые многоканальные частотомеры. Проб­ лемы цифрового спектроанализа в настоящее время при­ влекают все большее внимание широкого круга спе­ циалистов по автоматической обработке информации. Автоматическое управление, радиоастрономия, гидроло­ кация, радиолокация — вот далеко не полный перечень областей науки и техники, связанных со спектральным анализом случайных процессов.

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

При эвристическом синтезе цифровых частотомеров используются некоторые характерные приемы. В филь­ тровых частотомерах используются фазовые отсчеты входного процесса (В . 1) . В корреляционных частотоме­ рах при определении взаимнокорреляцнонной функции входного процесса и опорного колебания чаще всего ис­ пользуется полярный метод с предварительным центри­ рованием [4]. В обоих случаях применяется бинарное квантование. Малоразрядное преобразование в обоих

272

случаях значительно упрощает техническую реализацию устройства, а потерн во времени наблюдения, оценивае­ мые по относительной эффективности (3.37), не превы­ шают 3-4-6 дБ.

Однако оба класса частотомеров имеют один недоста­ ток, который в зависимости от области применения су­ ществен в большей или меньшей степени. Если анализи­ руемый спектр является достаточно гладким и содержит не более одной ярко выраженной дискретной компонен­ ты, то кроме упомянутых потерь во времени наблюдения никаких дополнительных эффектов не возникает. Если же в спектре имеются по крайней мере две дискретные ком­ поненты, причем одна из них значительно слабее другой, то происходит подавление слабой компоненты и появ­ ляются ложные спектральные составляющие {88, 89]. По­ добные системы являются квазилинейными при любом числе дискретных компонент, если все они примерно оди­ наковы по мощности и для каждой из них выполняется соотношение q<0,7 [86].

•В автокорреляционных анализаторах спектра вначале вычисляется коэффициент корреляции R(x) анализируе­ мого процесса и(1). Технически наиболее просто реали­ зуется полярный метод [4], при котором центрированный процесс u(t) дискретизируется по времени и бинарно квантуется по нулевому уровню. Для нормального слу­ чайного процесса u(t) вероятность совпадения его зна­ ков, т. е. появления комбинаций 11 или 00, равна [46]

со оо

р (т) = j * | ш)2 (и,, u„) du^du^ -\-

0 о

оо

- f - j "

j " доа (и,, M,)'rfu,cfu2 = - ^ - a r c c o s [ — R (t)].

— C O

— C O

Вероятность совпадения знаков одной полярности равна (1.100)

/><и)(т) =РФО)(Г) = (1/2я) arccos[—і?(т)].

В соответствии с этими соотношениями легко получить

выражение для коэффициента

корреляции

 

R

(т) = —cos [лр (т) ] =

^-cos [2яр(И) (т)] =

 

 

=—cos {2л/?(оо) (т)].

 

Как показано в [4], это соотношение для R(x)

справед­

ливо и для

синусоидального колебания u(t)

sin(coo^+

18—1410

273

+ ф) с постоянной

амплитудой Е, частотой соо и равномер­

но распределенной

в интервале (0,2я) фазой ср. Корреля­

торы совпадения полярностей хорошо известны и описа­ ны в литературе [4, 39].

Энергетический спектр анализируемого процесса определяется далее с помощью алгоритма быстрого пре­ образования Фурье [122] в универсальной или специали­ зированной ЭЦВМ . При этом одновременно осуществля­ ется коррекция (уменьшение боковых лепестков) путем введения весовой функции [39, 122], правда за счет неко­ торого снижения разрешающей способности. Наряду с полярным методом вычисления коэффициента корреля­ ции находят применение и другие методы: множитель­ ный, релейный, вычисление через условное среднее И Т . д . [4].

Перейдем теперь к фильтровым частотомерам. При­ мером подобного устройства является фазовый корреля­

тор

(см.

п. 2.5.3). Такая функциональная

схема

(рис.

2.33)

лишь поясняет принцип действия и не являет­

ся экономичной при ее технической реализации. С целью упрощения весовые коэффициенты принимаются равными

1 или_0, что соответствует

съему

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

Q и Q триггеров

регистра

сдвига. Для

суммирования

кодов, снимаемых

с регистров

сдвига,

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

использовать параллельные счетчики, 'принципы построе­ ния которых описаны в книге [93]. И, наконец, число раз­ рядов регистра сдвига может быть взято / = 7- m l n /7'A , а не п^>1. Здесь f m i n = l/Timnминимальная частота в спектре анализируемого процесса, 7д<;1/2/тах, /шах максималь­ ная частота анализируемого спектра, n = Tu/TR, Тн — вре­ мя наблюдения (накопления). Иными словами, входной процесс разбивается на п/l участков, каждый из которых обрабатывается в двух регистрах сдвига (двух квадра­ турных каналах) и параллельных счетчиках. Дальней­ шая обработка (накопление) в пределах интервала Тв зависит от метода сканирования антенного луча РЛС . При сканировании с постоянной скоростью накопление проще всего осуществляется с помощью оператора экспо­ ненциального сглаживания, который может быть выпол­ нен на основе накапливающего сумматора (3.73). В РЛ С с электронным сканированием, когда антенный луч фик­ сируется в заданном направлении в течение времени Г ш накопление производится в цифровом интеграторе (на­ капливающий сумматор). Число параллельных счетчиков,

274

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

каждом квад­

ратурном канале равно М= (fm a x—/mm)/Af,

где Af по­

лоса пропускания парциального фильтра.

 

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

всего осуществляется

с помощью

схемы, изображенной

на рис. 3.23,а. Входное

колебание

u(t) после гетеродини-

рования и двухстороннего амплитудного ограничения по­ ступает на первые входы схем Иі, И2 , а на вторые их входы подаются сдвинутые на 1/4 периода следования

5/*

гл

а)

б)

Рис. 3.23. Квантование фазы на четыре уровня.

импульсов дискретизации Г д последовательности импуль­ сов их и ы2- С выходов г/і и у2 снимается двоичный код фазы колебания у {t) относительно опорной частоты /= = l/TR. Соответствие между фазой и двоичным кодом та­ ково:

«5

О-Иі/2

и/2-=-ти

И-4-3-/2

Зт!/2Ч-2тс

Ф

 

 

 

 

код yty2

11

ГО

о о '

01

Эпюры напряжений, поясняющие работу схемы, изо­ бражены на рис. 3.23,6.

3.5.2. Цифровой частотный дискриминатор. Построен­

ный эвристически цифровой

частотный

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

с аналоговыми сигнальным и опорным входами

изобра­

жен на рис. 3.24 {123]. Принятое колебание u(t)

после

двухстороннего ограничения

совместно

с опорной

часто-

18*

275

топ Won(0 подается -на фазовые детекторы двух квадра­ турных каналов. Выходное напряжение фазового детек­ тора каждого квадратурного канала подается на счетные входы триггеров. Выходные напряжения триггеров и фор­ мирователей коротких импульсов, соответствующих по­ ложительным перепадам потенциала, поступают и а логи­

ческую схему, состоящую из схем И и ИЛИ.

Если

частота входного колебания превышает опорную

a(t)>

> ( 0 о п ( ' 0 . то на вход реверсивного счетчика по суммирую­ щему входу будут поступать счетные импульсы, если же со (0<со 0 ц (0. т о счетные импульсы будут поступать по

1 I 1 - І

Рис. 3.24. Цифровой частотный дискриминатор.

276

вычитающему входу. Двоичный код разностной частоты снимается с выходов разрядов счетчика. Реверсивный счегчик выполняет функции интегратора. Данный частот­ ный дискриминатор является астатической системой, так как выходной код пропорционален f [w(t)—a)on(t)]dt.

Если же со(^) =со 0 п ( 0 . но входное и опорное

напряжения

не совпадают по фазе, то двоичный код

реверсивного

счетчика не изменяется.

 

 

 

Вместо реверсивного

счетчика

можно

использовать

устройство,

реализующее

оператор

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

сглаживания

(3.73). В этом случае

выходной код будет

пропорционален разностной частоте, усредненной за вре­ мя, определяемое параметром сглаживания а (3.73).

3.5.3. Цифровой фазовый дискриминатор. Временные интервалы между моментами пересечения нулевого уров­

ня опорным напряжением u0n(t)

и напряжением сигнала

(B.I) u(t)

приближенно можно

использовать для

опре­

деления фазы

напряжения

u(t)

относительно

опор­

ного напряжения Uon(t).

Рассмотрим

один из возможных

вариантов

построения

цифрового фазового

дискримина­

т о р а — дискретного

фазометра

(ДФ) [123].

Структурная

схема ДФ и эпюры напряжений, 'поясняющие его рабо­

ту, приведены на рис. 3.25. На вход умножителя

частоты

подается

опорное

колебание u0n(t),

частота

которого

умножается в

г4

раз число уровней квантования

фазы, на рис. 3.25

ГІ = 8). С помощью

детектора

«нулей»

процесса

«оп('0

(Дет. «О») выделяются импульсы, соот­

ветствующие моментам времени, когда колебание часто­

ты 'Гій)оп(0 проходит через нулевой уровень

сверху вниз

(с отрицательной производной). В качестве

детектора

«нулей», может быть использован, например, триггер Шмитта. Импульсы с выхода формирователя Fi подают­ ся на вход счетчика с малым временем переходного про­ цесса, работающего в циклическом режиме счета. На­ пряжение частоты a(t) подается на вход детектора «ну­ лей» (Дет. «Ог»), выход которого соединен со входом формирователя F2. Назначение и работа Дет. «Ог» и F% аналогичны рассмотренным выше.

На вход схемы И подаются импульсы F2, а также за­ держанные на время Хсч переходного процесса (установ­ ления кода) в счетчике импульсы Fi. Таким образом, импульс считывания фазы г|з,- в двоичном коде поступает с выхода схемы И только после окончания переходного

277

процесса в счетчике. Несинхронность «нулей» опорной

частоты

и колебания u(t) учитываются в

формирователе

Fi импульсов длительностью

т і ^ Г і , где

Ту — период

ко­

лебания

гісо о п (0 - Двоичный

код фазы

ярі ( і = І , 2,

... )

снимается за каждый период частоты a>(t) с выхода ре­

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

(PC). Если фазы

колебаний u(t) и

f

Aem.„Qt

 

Счетчик

um(t) '»,f

 

 

t—і

 

 

 

 

0

 

 

 

5*

 

 

 

Fz

 

 

 

¥ A

a(t)

 

 

 

 

aj

 

Двоичный под

Рис. 3.25. Цифровой фазовый дискриминатор.

278

Uon(t)

совпадают (кривые /, 4(a),

рис. 3.25), то

с

выхо­

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

го. При сдвиге фазы частоты u(t)

как в

сторону

опере­

жения

'[кривая 4(6)],

так и запаздывания

[кривая

4(e)]

происходит

изменение

величины

считываемого

кода

фазы [рис. 3.25,6 эпюра а|п(/)].

 

 

 

 

Отсчеты

фазы

tyi(t)

обычно производятся

из

смеси

полезного сигнала и помехи, поэтому часто выполняют усреднение на интервале времени Тн , на котором ty(0 = = const. Однако непосредственное усреднение может при­ вести к грубой ошибке [124], так как фаза смеси сигнала и помехи представляет собой циклическую случайную ве­ личину, определенную на интервале (0—360°). В самом деле, пусть, например, отсчеты фазы группируются отно­

сительно точки

я|зо = 0°: і|н = 300°, г[)2=330°, 1рз=30о , і|34 =

= 60°. Усредняя,

находим

 

4

 

; = i

тогда как истинное значение средней фазы в данном слу­ чае равно i|}o=0° (360°). По этой причине такое усредне­

ние возможно

лишь

для

гармонических

функций

ty(t).

 

Дискретные

же

отсчеты фазы с выхода ДФ необходи­

мо

подать в

блоки

вычисления

эшФг, cosФх, которые

за­

тем

должны

усредняться

в

накопителях

£ ,,

2 2 .

Далее

следует

образовать

 

отношениг

S,/S2 = t g ^ , '

в

соответ­

ствии с

которым 41

=

Arctg

(2,/Е,).

 

 

 

 

 

Основное ограничение в применении подобного фазо­

вого дискриминатора связано

с требованием

высокого

быстродействия счетчика, регистра совпадений, а также

устройств

схемы усреднения.

Так,

например,

при

fou=

=

10 МГц,

г і = 1 6 , r i X f i = 1 6 0

МГц.

При этом

Дет. «0i»,

Fu

счетчик, регистр совпадений (PC) должны

работать

на

частоте

160 МГц, а ДеТ. «0г» и блоки усреднения

на частоте

10 МГц.

 

 

 

 

 

3.5.4. Синтез квазиоптимальных

измерителей. Синтез

широкого класса измерителей параметров связан с реше­ нием уравнения максимального правдоподобия (3.2). Для марковской аппроксимации последовательности дис­ кретных переменных {и} на выходе квантизатора сущест-

279

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