Здесь |
|
|
|
|
|
|
|
U* У) = — |
$ |
/?„ (s) 11а (f, |
s) ds, |
£/р (*) = |
— |
j R 0(s) ilß (г*, |
s) fife; |
|
—oo |
|
|
|
—cc |
|
в частности, |
для |
четвертой |
модели в силу |
(13) |
имеем |
|
|
|
|
^ ds |
То / 2 |
|
|
|
и а { і ) = - ^ - Ъ |
^ 7?o(s) р/Г(^—s) X |
|
|
|
h l п >\t I |
0 |
|
|
|
|
|
|
' |
- с о |
— Г о / 2 |
|
|
|
X Рл(/ ■—s) U0 (t—и —А) и * (s—Д—ф) 0а dft. |
(20) |
Представим искомые функции R 0 (/), /?а (/), /?р (/) как интегралы Фурье от спектральных плотностей Gro if), Gro, if), GRß if). Тогда пу тем Фурье-преобразования интегральные уравнения (19) с разност ным ядром сведем к алгебраическим для спектральных плотностей
G j?o ( Ш 1 + * ( f ) ] |
= |
G u a (f), |
|
GRa( f n i + K ( m = Ga (f), |
(21) |
G«ß ( f ) l l + * ( D l |
= |
G p (f) . |
|
Здесь x if) — нормированное относительно уровня шума значение спектральной плотности мощности помехи, заданной (в нулевом при ближении четвертой модели) выражением вида
оо |
|
х(/) = § ті(т)е-/'2л^ т , |
(22) |
—с» |
|
Ga if) и Gß if) — спектральные плотности функций Ua it) и Щ (/). Решение для модели 3 не содержит членов с множителями или ин
дексами ß. Оно отличается от предыдущего заменой г| на ф в выраже ниях для к if) и Ua it). Для модели 3
|
|
со |
|
|
|
|
|
|
|
7t(f)= J ф ( т ) е - '2я^ т , |
|
(23) |
|
|
— со |
|
|
|
|
|
|
Ua it) —~~ |
Ra is) 9p i t s ) |
и A i t - fl) U% is - |
fl) fl2 ds dtt. |
(24) |
|
Nn JO |
|
|
|
|
|
|
—oo |
|
|
|
|
|
|
Выражения (21) |
для Grk if) и |
|
if) допускают дальнейшие пре |
|
образования. Применительно к модели 3 и у (F) = |
ö (Т) можно полу |
|
чить |
|
|
|
|
|
|
|
aGRaf ) |
aGa if) |
a |
2 |
J L U+x(f)l- |
(25) |
|
g ro (/) |
JUA if) |
|
3TJ |
T, |
df2 |
|
|
где Ttf = 2 ]/a — протяженность отражений по уровню 1/е.
§ 2.3.7. ВЫВОДЫ ИЗ АНАЛИЗА ОБНАРУЖЕНИЯ ДВИЖУЩЕЙСЯ ТОЧЕЧНОЙ ЦЕЛИ В ОБЛАКЕ ОТРАЖАТЕЛЕЙ И ИХ ОБСУЖДЕНИЕ
1.Имеет место определенная зависимость алгоритма оптимальной обработки от структуры облака.
2.Для сосредоточенного распределения отражателей по дальности (модели 1, 2) оптимальная обработка в целом (15) строго распадается
на |
межпериодную, |
характеризуемую набором коэффициентов Н„,, |
и |
внутрипериодиую, |
характеризуемую функцией £70 (/). Последняя |
остается такой же, как и в отсутствие мешающих отражателей. Межпе риодная обработка зависит от корреляционных моментов напряжений
D mß после внутриперподной обработки, |
несколько различающихся |
для моделей 1 п 2. |
рассматривать модель с 0 - |
Может |
возникнуть вопрос, стоило ли |
с р е д о т |
о ч е II и о й пассивной помехи? Действительно под таким |
наименованием она обычно не фигурирует. Однако часто заранее при нимают внутрипериодиую обработку согласованной, т. е. такой же, как II в отсутствие пассивной помехи. Это не нарушает оптимальности обработки в целом, когда отражения от помехи и цели имеют одинако вую структуру, т. е. помеха столь же сосредоточена по дальности, как п сама цель.
Наряду с различием в распределении отражателей цели и помехи по дальности, при обзоре пространства сказываются различия в рас пределении отражателей по угловым координатам. В литературе, от казываясь от рассмотрения углового обзора, часто исходят из форми рования пачечных сигналов в передатчике. Учесть различия в угловом распределении отражателей цели и помехи при этом также невоз можно.
Совместное использование предположений о согласованной внутрипериодной обработке и формировании пачки в передатчике [431 эквивалентно оптимизации обработки в целом с использованием моде ли сосредоточенной по дальности и угловой координате помехи.
2. Третья модель соответствует помехе, сосредоточенной только по угловой координате. Ей также соответствует замена обзора по угло вой координате излучением пачки импульсов ограниченной длитель ности. Нулевое приближение модели [70] может эффективно исполь зоваться, если протяженность облака в единицах временного запазды вания превышает длительность сигнала. Пределы применимости при ближения количественно определяются из [(25), § 2.3.6].
Пусть излучаемый сигнал является одиночным экспоненциальным
импульсом Uл (t) = Uü it) = е~і/хя (/ >- 0), принимающим нулевые зна чения для t < 0. В отсутствие разброса скоростей отражателей (т) = = 1 функция корреляции отраженных от облака колебаний согласно
[(4), (8), §2.3.6] будет ф (х) = (D0/2Af0)e_l х l/T“, нормированная спек тральная плотность их мощности % (/) = (О0/уѴ0)ти2 (1 -|- 4я/2г2п)_1 (далее полагаем D0xn2/N0 = 1). В случае 2я/т„ « 1 относительная поправка [(25), § 2.3.6)1 наиболее заметна и будет
a GHa (f)/GR0(f) ~ 2хЦх%. |
(1) |
310 |
§ 2.3.7. |
Значение (1) сравнительно мало, если протяженность отражений т* превышает в 3—4 раза длительность единственного облучающего им пульса т„. Когда же тп > Тф, кривая х (/)/т$ сужается, роль второй производной в [(25), § 2.3.6) резко возрастает и относительная по правка (1) может стать заметно больше единицы.
Аналогичное резкое возрастание поправки [(25), §2.3.6] имеет мес то, когда зондирующий сигнал представляет собой пачку радиоим пульсов, период которых заметно превышает протяженность отраже ний Тф. Действительно, спектр пачки состоит из гребней, ширина ко торых обратно пропорциональна длительности пачки. Чем меньше эта ширина, тем больше вторая производная в [(25), § 2.3.6]. Соответ ственно резко возрастает относительная поправка первого приближе ния. Использование третьей модели Е этом приближении не впол не подходит поэтому для тех реальных ситуаций, когда не только дли тельность пачки, но и период посылки значительно превышают времен
ную |
протяженность отражений Тф. |
4. |
Четвертая модель, как это следовало из постановки задачи, |
более правильно, чем третья, описывает такие реальные ситуации, ког да длительность отдельных отражений от облака существенно меньше длительности пачки.
Это же обстоятельство можно проиллюстрировать сравнением по правочных членов [(20), (24), § 2.3.6]. Примем для упрощения, что им пульс U0 (() и пачка Ua (t) имеют ограниченную длительность своих ненулевых значений. Функции /?0 (s) [(20), (24), § 2.3.6] приближают ся к последовательностям импульсов, имеющих длительность порядка сигнального, моменты действия которых соответствуют моментам дей ствия импульсов полезного сигнала. Тогда произведения /?0 (s) Uo (s —
— ft — •&) О2 и R0(s) U a (s — О-) О2 в [(20), (24), § 2.3.6) принимают не одинаковые максимальные значения. В первом из этих произведений соответствующие моменты времени s близки к моментам действия како го-то полезного импульса, иначе будет мала величина первого сомно жителя. Этот полезный импульс должен иметь номер /, в противном слу чае (при ограниченных пределах интегрирования по Ф) будет мала ве личина второго сомножителя. Значение ■&не может существенно пре взойти длительности импульса, иначе максимуму одного из сомно жителей будет соответствовать нулевое значение другого. При этом ве личина О2 заметно не превосходит квадрата длительности радиоимпуль са U0 (/). Из аналогичных соображений величина й-2 в [(24), § 2.3.6] ограничивается квадратом длительности пачки Ua (()■ Таким образом, асимптотические поправки четвертой и третьей моделей должны иметь совершенно различный характер.
5. Выражения т] (т) и х(/) для четвертой модели, кроме разброса ско ростей отражателей и ограниченного времени облучения, учитывают обновление состава отражателей от зондирования к зондированию вследствие обзора по угловой координате.
На рис. 2.3.10 иллюстрируется построение нормированной функ ции автокорреляции т} (т)/ц (0) в виде произведения трех нормирован
ных автокорреляционных функций: |
флюктуаций помехи pF (т), моду |
ляции диаграммой направленности |
рл (т), модуляции передатчиком |
§ 2.3.7. |
311 |
ру (т). Последняя функция соответствует сумме по т в выражении 1(10 а), § 2.3.6], взятой после соответствующей нормировки.
Спектральная плотность к (/) является преобразованием Фурье функции г| (т). Непосредственно ее можно найти как свертку функций y'(F), I Ga (/) |2, I Gu (/) I2, где GA (f) и Gy (/) — спектральные плотности
A{t) и U(t).
6.Обработка для четвертой модели в нулевом приближении реали
зуется оптимальным фильтром с частотной характеристикой /<опт (/) =
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= G |о (/) е-і^Но, |
где |
GR0 (/) |
соответствует [(21), |
§ 2.3.6], |
а |
кон |
|
|
|
|
|
|
станта і0 подбирается из ус |
|
|
|
|
|
|
ловия |
осуществимости. |
Это |
|
|
|
|
|
|
соответствует |
рассмотренно |
|
|
|
|
|
|
му уже |
в § |
1.2.3 обобщению |
|
|
|
|
|
Л/т) |
теории обнаружения |
на фоне |
|
|
|
|
|
белого |
шума |
применительно |
т |
к |
і |
ш |
м |
|
к |
стационарному |
небелому |
|
шуму. Нормированная спект |
|
|
|
|
|
|
ральная |
плотность |
мощности |
|
|
|
|
S>FM |
|
этого |
шума |
определяется |
|
|
|
|
|
|
выражением |
|
а3 = |
1 |
+ |
к (/).87* |
ßA(T)j>a(T)pF(z) |
|
|
|
|
Оптимальный |
фильтр |
можно |
|
|
|
|
при |
этом составить из звена |
|
£ £ І |
|
|
с |
характеристикой |
1/сг (/), |
|
|
|
«обеляющего» |
шум, |
и звена с |
|
|
|
|
|
|
характеристикой G'ua |
(/)/о (/), |
Рис. 2.3.10. К построению функции авто |
оптимальной |
для |
|
сигнала |
корреляции |
пассивном помехи. |
на |
выходе |
обеляющего |
зве |
|
|
|
|
|
|
на. На рис. |
2.3.11 |
иллю |
|
|
|
|
|
|
стрируется |
построение опти |
мальной |
амплитудно-частотной |
характеристики фильтра |
KonT(f) = |
= [ Gua (f) |/[1 + |
|
и (/)] |
при обнаружении пачки периодически следую |
щих когерентных радиоимпульсов. Эта характеристика (рис. 2.3.11, в) имеет гребенчатую структуру. В ней можно выделить гребни подавле ния помехи и гребни накопления сигнала.
7. Характеристику рис. 2.3.11, в можно практически реализовать только при сопровождении по скорости или обзоре по скорости целей. Обычное же послекомпенсацнонное некогерентное накопление опти мально только при независимых (быстрых) флюктуациях импульсов сигнала от цели после компенсации. Приближенную оптимизацию мож но провести, заменяя в 1(1), § 2.3.5] (Ja {t) на U0 (t ) и сохраняя корре ляционную функцию помехи. Все это соответствует оптимальной филь трации одиночного радиоимпульса на фоне когерентной помехи с меж
периодной |
корреляцией. Построение характеристики | /<опт (/) | = |
= |G0 (/)(/[ 1 |
+ и (/)1 с гребнями подавления, реализуемой, например, |
системами многократного череспериодного вычитания, иллюстрирует ся на рис. 2.3.12.
8. Рассмотренные разновидности фильтрации обеспечивают не толь ко селекцию по скорости, ио и дополнительную селекцию по дально сти (см. § 1.2.3). Повышения разрешающей способности не наблюдает-
ся, если спектр сигнала имеет прямоугольную огибающую. Именно этот случай имелся в виду при построении рис. 2.3.11 и 2.3.12.
9. В отсутствие скоростных различий может наблюдаться эффект расширения отдельных гребней спектра (когда их форма отличается от прямоугольной), это означает сокращение длительности пачки в линей
ном фильтре, а для модели 4 — повышение разрешающей |
способности |
по угловой координате. |
|
10. Если зондирующий импульс U0 (/) имеет прямоугольный ампли |
тудно-частотный спектр, то а л г о р и т м о б р а б о т к и |
д л я и о |
д е л и 2 с сю т в е т с т в у е т н у л е в о м у п р и б л и ж е н и ю
JLJLJl A J t M
L U t z lf)
'0 f
JLJLJüL tlSf1
Рис. 2.3.11. Амплитудно-частотный спектр сигнала в виде когерентной пачки радиоимпульсов (а); распре деление спектральной плотности мощности результирующей помехи (б); амплитудно-частотная характе ристика оптимального фильтра (в).
Рис. 2.3.12. Амплитудно-частотный спектр одиночного радиоимпульса в пределах выбранной области частот (а); распределение спектральной плотности мощности результирующей помехи (б); амплитудно-частотная характеристика
фильтра (в).
м о д е л и |
4. Действительно, помножим произвольное |
р-е уравнение |
из уравнений [(16), |
§ 2 .3 .6] для, модели 2 |
на G*0 (/) еі2я^ |
и просумми |
руем все уравнения по р. Замечая, что величины Ь тр, и tm — |
зави |
сят только от разности т — р, сводя двойную сумму к произведению, |
получаем |
|
|
|
|
|
|
[ 2 |
Н т Go (!) е ~ /2я/Ч * V д , |Де - /2я/ |
= Gg (/) 2 |
е/ 2 |
я |
(2) |
L т |
|
|
J т —(д, |
jx |
|
|
|
Первый сомножитель в левой части равенства (2) — это |
К0ПТ (/), вто |
рой |
сомножитель |
представляет собой |
величину, совпадающую |
с |
І1 + |
X (/)], |
правая |
часть равенства (2) — это Gua (/)• Таким образом, |
мы действительно пришли к алгоритму для модели 4. |
|
|
|
11. |
Если повышение разрешающей способности по дальности и угло |
вой координате недостаточно эффективно, |
можно, таким образом, с рав |
ным успехом пользоваться моделями точечного и протяженного рас |
пределений помехи. |
При недостаточной эффективности скоростной се |
лекции, неравномерном спектре сигнала и интервалах менаду импуль сами, существенно превышающих временную протяженность помехи, предпочтительна модель 4.
§ 2.3.8. НЕКОТОРЫЕ ДОПОЛНЕНИЯ К § 2.3.5 — 2.3.7
Остался незатронутым ряд теоретических вопросов, на которых кратко оста новимся в данном параграфе, а именно:
— влияние изменения параметров помехи по дальности на алгоритм обна
ружения;
— влияние изменения параметров помехи по угловой координате на алго ритм обнаружения; ,
—качественные показатели обнаружения при известных статистических параметрах помехи;
—влияние некогерентностн сигнала на алгоритм и качественные показа
тели обнаружения.
1. В л и я н и е и з м е н е н и я п а р а м е т р о в п о м е х и п о д а л ь н о с т и рассмотрим применительно к случаю, когда оно происходит сравни тельно медленно, причем радиолокатор осуществляет периодическое импульсное зондирование. Используем метод приближенной замены параметров помехи (§2.2.12) и примем марковский закон изменения этих параметров по дальности. Все это вновь приведет к схемам рис. 2.2.16, с тон лишь разницей, что выходы элементов решетки (подключенных к ним усилителей) заменяются соответствую щими выходами линий задержки на время Т, 2Т и т. д. Сумматор 2 0 на схеме рнс. 2.2.16, а образует цепь когерентного накопления сигнала. Угловая само настройка схем рис. 2.2.16 по помехе заменяется при этом скоростной fill). Если когерентное накопление сигнала в схемах рис. 2.2.16 применить не удает ся, следует предусмотреть послекомпенсациониое некогереитное накопление.
2. В л и я н и е и з м е н е н и я п а р а м е т р о в п о м е х и по у г л о в о й к о о р д и н а т е может сказываться от элемента к элементу углового разрешения и в пределах каждого элемента. Изменение параметров от элемента к элементу учитывается в процессе самонастройки схем обработки, последова тельно во времени, как и при изменении по дальности. Изменение в пределах каждого элемента поддается учету путем многоканального приема, например с помощью фазированной решетки. Пусть, например, вместо одной из моделей облака [(3), § 2.3.5] используется модель
|
D (fi, F, 0 = |
da 5 (ф)ѵ (F - aQ e " K \ |
(1) |
Используя (1), решающую функцию найдем из уравнения |
|
|
ОО |
|
|
|
R |
(t, х) + JJ ф (t, |
s, |
X, $)R (s, £) dsdl = UA(t), |
(2) |
|
-----СО |
|
|
|
где X—координата |
раскрыва, ф (П s, |
х, | ) —нормированная |
автокорреляцион |
ная функция помехи. Принимая далее коэффициент а малым |
и используя метод |
последовательных приближений, решение (2) можно искать |
в виде |
R (/, х) = Ro(t, X) + о#! (/, А-). |
(3) |
Слагаемое R0 (t , х) соответствует независимым угловой и скоростной обработке, слагаемое (t, х) является добавочным и учитывает особенности пространствен но временной обработки: в каждом угловом направлении целесообразно преи мущественное подавление колебаний своей допплеровской частоты [136].
3. К а ч е с т в е н н ы е п о к а з а т е л и о б н а р у ж е н и я когерент
ного сигнала зависят от коэффициента использования энергии к. |
Например, для |
когерентного сигнала в случае модели |
т о ч е ч н о й помехи [(15), (16), § 2.3.6], |
энергетическое отношение сигнал/помеха [43] будет |
|
I S ' V |
J 2/ |
2 |
2 ^ ( D - i ) m(iV |
(4) |
ІИ |
I / |
m, n |
m, n |
|
Здесь Ф ~ 1)Пщ—элемент матрицы, обратной матрице ||Dmfi||. В отсутствие пассивной помехи Dnn= N оАт , а (0 -1)ш)Л= Д ЛѴ. Поэтому коэффициент использо вания энергии
|
|
( 5 ) |
При |
выводе (5) использовалось обозначение |
Amß для элемента единичной |
матрицы |
(Дш^ = 0 при т Ф р., Дш^ = 1 при /и = |
ц). |
4.Некогерентность полезного сигнала при наличии пассивной помехи
(также, как и при ее отсутствии, см. §2.2.4—2.2.10) ограничивает оптималь ное время когерентного накопления. Становится целесообразной замена коге рентного накопления некогерентным. Соответственно изменяются качественные показатели обнаружения [149]. Заметим, что при учете влияния малой некоге рентности выкладки существенно упрощаются при использовании методов, опи санных в § 2.2.6—2.2.7 (можно избежать обращения матриц высокого порядка).
Ч А С Т Ь Ш
В О П Р О С Ы О Б Щ Е Й Т Е О Р И И Р А З Р Е Ш Е Н И Я
Естественный переход от полного разрешения к квазиполному раз решению— обнаружению (см. § 1.1.1) упростил его анализ, позволяя расширить область применения теории. Необходимо, однако, углуб ление понятий самой теории, тем более, что, оставаясь на позициях об наружения, мы не касались вопросов измерения. Понятие полного разрешения позволяет единообразно за счет выбора матриц стоимости и априорных вероятностей перейти к различным случаям разреше ния — обнаружения — измерения, определения числа разрешаемых объектов и т. д., в том числе к случаю квазиполного разрешения —■об наружения. Этим вопросам посвящена гл. 3.1 [51, 61, 66, 115, 116J. В гл. 3.2 даются простейшие примеры разрешения — измерения [70, 75, 85а, 106, 115, 116, 144, 145]. Приступая к чтению § 3.1.1, жела тельно повторно просмотреть § 1.1.1.
Глава 3.1
РАЗНОВИДНОСТИ ПОЛНОГО РАЗРЕШЕНИЯ
§ 3.1.1. КВАЗИПОЛНОЕ РАЗРЕШЕНИЕ - ОБНАРУЖЕНИЕ КАК РАЗНОВИДНОСТЬ ПОЛНОГО РАЗРЕШЕНИЯ
Уже отмечалось (см. § 1.1.1), что квазиполное разрешение — об наружение — это частный случай полного разрешения, определяемый надлежащим выбором стоимостей решений.
Поясним это, интересуясь заранее каким-то одним процессом А ѵ из всей их совокупности (ѵ = 1..., т), например сигналом от цели, нахо дящейся на ѵ-м участке дальности. Среди 2т возможных с о б ы т и й наличия и отсутствия процессов выделим 2т~1 событий, в которых про цесс А ѵ отсутствует (нулевая группа событий, 1 — 0). Им можно отве сти, например, первые номера / = 0, 1, ..., 2т_1— 1. Последующие 2т_1 номеров отведем событиям, соответствующим наличию А ѵ (первая группа событий, / = 1).
Аналогичную нумерацию і и |
аналогичные группы с номерами k = 0 |
и k = 1 введем для р е ш е н и й |
о наличии и отсутствии процесса Аѵ. |
Беспокоясь о правильности |
только этих решений, для их стоимо |
стей а и (і, / = 0, |
1, ..., 2т — |
1) зададим всего четыре фиксирован |
ные градации a kl, |
зависящие от k и / (k = 0,1; I = 0,1). Последние на- |
значаются в зависимости от того, к какой группе, нулевой или первой относятся числа і и /; в соответствии с і выбирается значение А, а в со ответствии с /' выбирается /. Это как раз и означает, что стоимость решения выбирается в зависимости от правильности решения об одном только ѵ-м сигнале, что определяет квазиполное разрешение.
Выражение [(3), § 1.1.1] для среднего риска Q=QV можно предста вить тогда в виде суммы четырех слагаемых,
Q= У> |
“kiPiPki- |
(1) |
к, 1 = |
0 |
|
При а 00 = а и = 0 выражение (1) можно привести к виду [(5), § 1.1.2].
§ 3.1.2. КВАЗИПОЛНОЕ РАЗРЕШЕНИЕ — РАЗЛИЧЕНИЕ
Рассмотрим различение каких-либо т1 значений параметра интере
|
|
|
|
|
|
|
|
сующего нас сигнала п р и н а л и ч |
и и |
/п„ |
м е ш а ю щ и х с и г - |
н а л о в. Общее число |
разрешаемых |
в смысле |
полного разрешения |
сигналов составляет т = |
тг + т0, причем из т1различаемых сигналов |
с неодинаковыми значениями |
параметра |
одновременно имеется |
один |
и только один сигнал. Стоимость решений а и |
устанавливается в за |
висимости от качества оценки |
параметра |
этого |
интересующего |
нас |
сигнала. Разрешение, |
которое в данном случае имеет место, назовем |
к в а з и п о л н ы м |
р а з р е ш е н и е м — р а з л и ч е н и е м . |
Матрица стоимости квазиполного разрешения — различения содержит
ml |
2т ‘+"’о |
элементов. Выбор элементов матрицы стоимости может |
быть произведен по крайней мере несколькими способами. |
|
Так называемая п р о с т а я |
м а т р и ц а |
с т о и м о с т и пре |
дусматривает |
о д и н а к о в у ю |
стоимость ущерба a hl для всех не |
правильных |
решений А Ф I, не детализируя |
его степени. Например, |
для всех неправильных решений можно принять одинаковую положи тельную сс > 0, а для правильных — нулевую стоимость ущерба. Алго ритм решения н е и з м е н и т с я, если для неправильных реше
ний задать нулевую, а |
для правильных — о т р и ц а т е л ь н у ю |
стоимость ущерба а < 0. |
|
К в а д р а т и ч н а я |
м а т р и ц а с т о и м о с т и учитывает |
увеличение ущерба, пропорциональное квадрату ошибки в оценке па раметра. Если индексы А, I соответствуют последовательным равноот стоящим значениям параметра, значения а ц пропорциональны (А —
-If-
К в а д р а т и ч н а я м а т р и ц а с т о и м о с т и с н а с ы - щ е и и е м сводится к квадратичной при малой величине ошибок и приближается к простой, когда ошибки превзойдут некоторый порог (например, порог срыва автосопровождения) [115].
Остановимся на некоторых (в целом искусственных) способах введе ния квадратичных матриц стоимости с насыщением, позволяющих охватить с единых позиций широкий класс используемых критериев. Замечая, что тела неопределенности ограниченных по спектру сигна лов, имея квадратичные вершины, спадают по обе стороны от них
практически до нуля, элементы матрицы стоимости можно задать соот ношением
где С — константа, а рй, — коэффициент автокорреляции пли взаим ной корреляции каких-то сигналов.
В качестве таких сигналов могут быть взяты, например, сами раз личаемые сигналы. Полагая при этом
|
Ры — |
^ uk(t)ut(t)dt |
У |
|
I |
1Uh(/) I2 dt $ Ui(()\2 dt, |
(2) |
|
|
|
|
|
— со |
— о |
|
|
где |
Uk (/), |
Ut (/) — комплексные |
амплитуды, приходим |
к стоимости |
a kl, |
которая [115]: |
k |
|
|
|
|
|
1) обращается в нуль при |
= |
/, |
когда совпадают |
оцениваемое |
иистинное значения параметра;
2)изменяется пропорционально (/г — /)2, если разность значений параметра не превышает ширины пика автокорреляционной или взаимокорреляционной функции;
3)приближается к насыщению, когда рк1 приближается к нулю;
4)не зависит от начальных фаз сигналов.
Остановимся на одном видоизменении элементов (2) матрицы стои мости, переходя от комплексных амплитуд к мгновенным значениям, так что phl заменяется на
Pcofef= J |
uh(0 wi(0 dt / у |
jj |
u\{t)dt |
щ (/) dt . |
(3) |
—00 |
/ |
—20 |
—OO |
|
|
Величина akl приобретает при этом зависимость от разности фаз раз личаемых сигналов, резко возрастая, когда последняя приближается к я, Зя, ... Значения аы определяются при этом уже не квадратичной функцией (с насыщением), а некоторой синусоидальной, имеющей эту квадратичную в качестве огибающей своих минимумов.
Близкие свойства имеет матрица стоимости с элементами
ahl = C y u k (t) — Ui(t)]2dt, |
(4) |
— со
оо
которые не отличаются от элементов (1), (3), если ^ и\ (I) dt =
—00
со
= § и] (f)dt = 1/2. Матрица подобного вида учитывает несовпадение
—со
каждого из параметров различаемых сигналов, в том числе начальной фазы и амплитуды. По отношению к ошибкам амплитуды матрица с эле ментами (4) является квадратичной (без насыщения). Матрицы с эле ментами (1), (2) или (1), (3) несовпадения амплитуд не учитывали.