
книги из ГПНТБ / Лихарев В.А. Цифровые методы и устройства в радиолокации
.pdfРезультаты решения уравнения (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п Р(и\Х)/дХ2=В0, может быть вычислен заранее и введен .в оптимальный дискриминатор. Теперь уравнение дискриминатора принимает вид
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