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

книги из ГПНТБ / Ширман, Я. Д. Разрешение и сжатие сигналов

.pdf
Скачиваний:
82
Добавлен:
22.10.2023
Размер:
13.3 Mб
Скачать

Здесь

 

 

 

 

 

 

 

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.6.

309

§ 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). Повышения разрешающей способности не наблюдает-

312

§ 2.3.7.

ся, если спектр сигнала имеет прямоугольную огибающую. Именно этот случай имелся в виду при построении рис. 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

0

,

F

 

в)

 

Рис. 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.

Если повышение разрешающей способности по дальности и угло­

вой координате недостаточно эффективно,

можно, таким образом, с рав­

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

пределений помехи.

При недостаточной эффективности скоростной се­

§ 2.3.7.

313

лекции, неравномерном спектре сигнала и интервалах менаду импуль­ сами, существенно превышающих временную протяженность помехи, предпочтительна модель 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

 

314

§ 2.3.8.

Здесь Ф ~ 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 (можно избежать обращения матриц высокого порядка).

§ 2.3.8.

Ч А С Т Ь Ш

В О П Р О С Ы О Б Щ Е Й Т Е О Р И И Р А З Р Е Ш Е Н И Я

Естественный переход от полного разрешения к квазиполному раз­ решению— обнаружению (см. § 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т~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). Последние на-

316

§ 3.1.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].

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

§ 3.1.2.

317

практически до нуля, элементы матрицы стоимости можно задать соот­ ношением

Чі = С( 1— P j ,

(1)

где С — константа, а рй, — коэффициент автокорреляции пли взаим­ ной корреляции каких-то сигналов.

В качестве таких сигналов могут быть взяты, например, сами раз­ личаемые сигналы. Полагая при этом

 

Ры —

^ 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) несовпадения амплитуд не учитывали.

318

'

§ 3.1.2.

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