
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfп—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