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

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

.pdf
Скачиваний:
92
Добавлен:
30.10.2023
Размер:
14.48 Mб
Скачать

где Т — длина интервала времени,

в

котором исследуется

появление сигнала.

 

 

формы корреляционная

В случае сигналов произвольной

функция

 

 

 

 

 

!х(^—.To)=v('c)v('co)=y] QghmgM

(56.18)

стремится к нулю при

g. h

 

 

 

| т — т01 —> оо.

Если при

условии

 

 

| ч; — т01 > Ди

 

 

(56.19)

функция р. (т — т0)

принимает малые

значения,

котовыми

можно пренебречь,

то

при дискретных значениях (56.14)

опять можно считать, что выполняется условие ортогональ­ ности (56.10).

Задача об М ортогональных сигналах имеет и другие

интерпретации и применения. Например, сигналы могут

быть разнесены не по времени прихода,

а по частоте и

вследствие этого удовлетворять условиям

(56.10).

Рассмотрим Л4-канальный приемник,

принимающий ре­

шение о

наличии сигнала m(t, k.)

при условии

А )>

[при всех k = 1,..., /И,

отличных от k], (56.20)

k

 

что такое измерение дает

Вычислим вероятность того,

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

Afe (или zk) параметр k

равен истинному значению k0.

Обозначим эту вероятность чеоез Р(г); она,

очевидно, за­

висит от значения г =

Пусть F(z) есть

вероятность

ложной тревоги в каждом канале при пороге z (если этот

канал использовать для обнаружения каждого из сигналов m(t, k) в отдельности). Тогда вероятность Р(г) в силу не­

зависимости случайных величин \ в М — 1

канале (k

^/г0) равна

 

Р(г) = [1 — F(z)]M~1.

(56.21)

Если плотность вероятности величины z в &0-м канале, где

действительно имеется сигнал

m(t, kB), обозначить

через

p(z), то вероятность

Рт — вероятность произвести правиль­

ное измерение k при

любом значении z— будет равна

ОО

 

— F(z)]M~xp(z)dz

(56.22)

Рт= J P(z}p(z)dz = J [1

б

oJ

 

 

349

или

 

М—1

 

 

со

 

 

 

,И—1

 

 

р.=S CL, (-1 )‘р' (ч Р w * = У 1 )‘с^7‘.

 

й=0

 

 

0

 

 

 

 

k= О

 

(56.23)

где

есть биномиальный

 

коэффициент

 

 

 

 

СМ-1 — k\ (/Ц — 1 _ k)l

 

 

 

k\

 

 

 

В случае сигнала с

 

неизвестной фазой,

когда

 

 

 

 

 

 

 

 

 

 

_|л+г*

 

 

 

 

 

/?(г) = е

2,

/?(г)==ге

2

/0(/и-г),

 

(56.25)

пользуясь тождеством (33.46), имеем

 

 

 

 

 

_

__ kz2

 

<х>

__ р- + (& + 1) г»

 

 

 

 

 

F*

= e

 

 

 

 

 

2 /0(f^)dz=

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

е

!«+»

 

-

------\.

*тт гГг=

 

 

 

=t+i-J

ге

 

 

 

 

 

 

 

О

 

 

 

 

 

Л|л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

е~ 2(i+Tf

 

(56.26)

 

 

 

 

 

 

= А + 1

 

 

 

 

 

 

 

 

 

 

откуда

из формулы (56.23) получаем следующее выражение:

 

 

 

 

 

 

 

 

е~2(*+1)

 

<56-27)

 

 

^^(-чЧ-.-т+т--

 

 

 

 

А=0

 

 

 

 

 

 

 

В случае флюктуирующего

сигнала

с неизвестной

фазой,

когда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(z) = e

_ —

 

 

 

 

z*

(56.28)

 

 

2,

р(г) = гА-е 2ll+^,

мы получаем более

простое выражение

 

 

 

 

 

М—1

 

 

 

 

 

 

 

 

 

 

р =У(—...(56.29)

v

 

 

Zj v

 

'

м-11 + k (1 + р.) ’

*=0

350

поскольку тогда

ОО /

—г 1 Г ~

F'=r+^jze '

о

I__ \

zs

1

Г+Г

Т,

dz=i+Mi+rt. (56.30)

При двух сигналах (7И = 2) формула (56.27) принимает вид

 

 

 

 

 

 

 

 

<56-31)

а формула (56.29)

дает

 

 

 

 

 

 

 

 

Р„=1-2-^-

 

 

(56.32)

При любом М в обоих случаях имеем

предельные соотно­

шения

 

 

 

 

 

 

 

 

 

lirnPm= 1,

ИгпРт = ~,

(56.33)

 

р.->оо

 

 

р.->0

 

 

 

которые, разумеется, можно

было ожидать заранее. Вто­

рое соотношение (56.33)

является следствием тождества

 

М—1

 

 

М—1

 

1

 

 

2 (-1)4-1 *Фт=Е<- 1)4с«-' j

 

 

*=0

 

 

 

k=0

 

о

 

1

М—1

 

 

 

1

 

 

 

= \dx J, (- l)ftC^_1Xfe

= j(l-x)M-1rfx=4- . (56.34)

6

Л=0

 

 

о

 

 

 

 

При больших М расчеты по формулам (56.27) и (56.29)

становятся громоздкими, поэтому целесообразно

вывести

асимптотические

формулы,

 

пригодные при М >оо. Для

этого можно аппроксимировать функцию

(56.21)

следую­

щим образом:

 

 

при з<г0, 1

 

 

 

Р(г) = 0

 

(56.35)

 

 

р (z)= 1

при Z>20, J

 

 

причем

г0 будем

определять из требования, чтобы при

г = г0 точная функция (56.21) была равна

 

 

 

Р(г0) = [1-^(г0)Г_1=4

(56-36)

или

 

 

 

 

 

 

 

 

 

F(z0) = e

2

 

= 1—2

м

’.

(56.37)

351

Отсюда

z0=l/-21n(l-2 м~‘)

(56.38)

Действительно, при увеличении Л1 переход функции Р(г) от значений P{z)^0 к значениям Р(г)я=,\ стано­

вится все более резким, поэтому аппроксимация (56.35)

дает все более точные результаты. Подставляя выражение

(56.35) в выражение (56.22), получаем

Рт=\ P{z}dz=D{zol

(56.39)

т. е. вероятность правильного измерения равна

вероят­

ности правильного обнаружения одного сигнала при использовании порога (56.38), зависящего от А/, для которого можно написать приближенное выражение

 

 

2о = 1/ 2 In(при/И > Г).

(56.40)

 

 

 

 

На рис.

54

и 55 приведена зависимость вероятности

Рт

от отношения

сигнал/помеха ц,

вычисленная для Л4 = 2,

4,

8 и 16 по

точным -формулам

и для Й4=16, 200, 2 00Э

и 20 000 — по приближенным. На рисунках видно, что уже при А4=16 расчеты, произведенные по приближенным фор­ мулам, дают достаточно удовлетворительные результаты Из предыдущих формул следует, что разность 1—Р при увеличении р уменьшается для сигналов с постоянной амплитудой гораздо быстрее (экспоненциально), чем для сигналов с флюктуирующей амплитудой. Таким образом, правильное измерение дискретного параметра нефлюктуи­ рующего сигнала с большой надежностью произвести го­ раздо легче, чем такое же измерение для флюктуирующе­ го сигнала, а для мало надежного различения сигналов с постоянной амплитудой требуется, наоборот, большее отношение сигнал/помеха. Это обстоятельство мы уже обсуждали раньше в теории обнаружения.

Из рис. 54 и 55 и формул (56.39) и (56.40) видно, что

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

дит под знаком логарифма.

В § 55 мы вывели для вероятности измерения непре­ рывного параметра х неравенство (55.02), в правой части

352

101^

Рис. 54. Зависимость вероят­

Рис. 55. Зависимость вероят­

ности правильного различе­

ности правильного различе­

ния Рт от отношения сиг-

ния Рт от отношения сиг-

нал/помеха в одном канале

нал/помеха в одном канале

при неизвестной фазе,

при неизвестной фазе и не­

сигнала.

известной амплитуде

Сплошные кривые — по точной

 

сигнала.

 

формуле (56.27), пунктирные — по

Сплошные кривые — по

формуле

приближенной формуле (56.39).

(56.29),

пунктирные — по

прибли­

 

 

женной формуле (56.39).

которого фигурирует полная

вероятность ошибки при раз­

личении двух сигналов tn (t, т -|- е)

и m(t, -с— е). Эти сиг­

налы при достаточно большом г можно считать ортого­

нальными, и тогда

по формулам (55.02)

и (55.08) мы по­

лучаем

 

 

 

/>(1*

-

 

(5б-41)

где

 

 

 

 

оо

г’

 

А = у=- J

 

(56.42)

23—483

353

В данном

параграфе

мы

показали,

что различить

М

ортогональных сигналов

при

весьма

большом 34 труднее,

чем при /14 = 2

или

34 = 4.

Поэтому

ясно, что при

боль­

ших значениях

М,

получающихся

 

по

формуле (56.17),

вероятность ошибки

Р (| х — х01 > е)

должна

значительно

превышать величину,

стоящую в правой

части

неравенства

(56.41).

 

 

 

 

 

 

 

 

 

 

 

В настоящее время лишь рассмотрение „системы М

ортогональных сигналов"

позволяет

получить какие-то ко­

личественные

(хотя

и приближенные)

результаты

о

воз­

можностях измерения параметра сигнала в широком интер­

вале значений этого параметра. Только при вероятности

Рдостаточно близкой к единице, можно считать, что

грубое измерение непрерывного параметра х достаточно надежно и что данный максимум коэффициента правдопо добия связан с сигналом, а не является случайным выбро­

сом. Если при разбиении априорного интервала 0<х<7'

на М частей в соответствии с формулами (56.14), (56.15) и (56.19) получается вероятность Рт, достаточно близкая к

единице, то можно ставить вопрос о дальнейшем уточнении значения х, т. е. более точном измерении.

Теоретические возможности, связанные с точным из мерением произвольного параметра т, были рассмотрены

в § 52 и конкретизированы в § 53—55 для измерения вре­

мени прихода сигнала. Если сигнал имеет

постоянную

амплитуду

и помехи

достаточно

слабы,

то

измеренное

значение

т

является

нормальной

случайной

величиной,

среднее

значение которой равно

То — инстинному

значе­

нию т, а дисперсия определяется формулой

(53.16).

Также

можно пользоваться

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

(55.02)

и

рис.

53, ко

торые должны давать приближенную оценку ошибки при «точном» измерении, проводимом вслед за грубым изме­ рением, поскольку тогда при «точном» измерении интервал измерения т достаточно узок.

На практике измерение параметра сигнала при наличии помех естественно разбивается на две различные опера­ ции: грубое измерение и точное, причем результаты точно

го измерения должны приниматься в расчет только при их согласии с результатами грубого. Заметим, что экспери­ ментальное исследование функций Л(т), <р(т) или$(т) для непрерывного параметра часто возможно лишь в резуль­

тате построения этих функций по отдельным точкам

(56.14), полученным в результате разбиения исходного ин-

354

тервала на большое число малых интервалов Ат, причем

помехи в разных интервалах являются приближенно независимыми. Таким образом, задача об М ортогональ­ ных сигналах отображает свойства реальных технических устройств.

Задача о системе М. ортогональных сигналов была впервые поставлена и решена (в более частных предполо­ жениях) В. А. Котельниковым.

§ 57. ОБНАРУЖЕНИЕ СИГНАЛА С ДИСКРЕТНЫМ ПАРАМЕТРОМ

ИИЗМЕРЕНИЕ ЭТОГО ПАРАМЕТРА

Впредыдущих параграфах мы рассматривали лишь

«чистое» измерение, т. е. предполагали заранее известным

присутствие полезного сигнала, причем только одного — с некоторым неизвестным значением параметра а, т или k,

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

ется в дополнительном исследовании.

Обнаружение сигналов mk (/) —m(t, k) с дискретным

параметром k (&='!,...,Л1) и измерение этого параметра можно производить с помощью различных AI-канальных приемников, которые мы будем условно называть прием­ никами I, II и III типов.

Приемник I типа производит обнаружение по коэффи­ циенту правдоподобия

 

 

 

 

м

 

 

 

 

 

Л=£аа,

 

(57.01)

а

именно,

 

Й=1

 

 

 

 

 

 

 

 

 

 

при

Л>Л решает,

что

сигналы

имеются,

(57.02)

 

при

А<А решает,

что

сигналов нет

 

 

и,

кроме того, при А> Л*

принимает

решение

о наличии

сигнала с

максимальным Aft. Этот

приемник нуждается в

модификации, если сигналы могут,

например,

появляться

попарно.

сигналы mk (t) других неизвестных параметров не

 

Если

имеют, то каждый канал

приемника

должен образовывать

величину

 

 

h

%

 

 

 

 

 

 

,

(57.03)

 

 

Л,= ^е

2

23*

355

где pk есть вероятность появления k-vo сигнала при усло­

вии, что какой-то полезный сигнал присутствует, а

определяется формулами (56.07). Таким образом, осуществ­

ление приемника I типа является с практической точки зрения довольно громоздкой задачей, поскольку после линейной обработки входных данных (при получении ?/г )

требуется образовывать величины (57.03) и затем их сум­ мировать по формуле (57.01). Однако прие^яик I типа

является оптимальным приемником обнаружения в смысле гл. V, так что характеристики всех других приемников

должны уступать характеристикам приемника I типа. По­ скольку вычисление характеристик приемника I типа свя­ зано с определенными трудностями (см. далее), детально

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

невелики.

Приемник II типа производит предварительное испы­

тание всех М каналов по правилу

 

 

при А,

>2.

сигнал m.(t) может быть,

i

(57.04)

k

к

"

7

J

при АД,<2Д;

сигнала mk (f) нет,

 

где 2^ — порог

в k-м.

канале

приемника. Если

 

во всех

каналах Aft<^2&, то выносится решение об отсутствии сиг­ налов, если хотя бы в одном канале Aft>2ft, то выносится

решение о

наличии сигналов. Если заранее известно, что

может быть

не

больше

одного

сигнала,

то

в

качестве

такового

предлагается

сигнал

с

максимальным Aft. Если

число сигналов не

ограничено,

то

испытание по

формуле

(57.04) является окончательным.

 

производит

измерение,

Приемник

III

типа сначала

 

т. е. выбирает из

всех Л4 величин Av . .. , Ам наибольшую,

а затем, испытывая выбранное значение Лй,

 

 

 

при Аа>2а решает, что есть

сигнал

 

 

i

при Аа<2а

решает, что сигналов нет,

 

 

)

где 2а—те

же пороги,

что и

в

формуле

(57.04). Этот

приемник можно применять только при условии, чтх> появ­

ление двух или более полезных сигналов mA(Z) a priori

исключено. Если же может появиться только один сигнал,

то приемники II и III типов, как легко видеть, эквивалент-

356

ны, поскольку они производят те же самые операции, но в другой последовательности.

Можно сказать, что в приемнике I типа сначала произ­ водится обнаружение, а затем измерение, в приемнике II

типа обнаружение и измерение проводятся параллельно,

а приемник III типа сначала измеряет параметр k, а потом проверяет наличие сигнала с измеренным параметром «на обнаружение». Приемники II и III типов осуществить на практике гораздо проще, чем приемник I типа, поскольку вместо отдельных Л А в каждом канале приемника можно об­

разовывать любую монотонно возрастающую функцию *А-

При исследовании характеристик введенных выше при­ емников мы ограничимся системой М ортогональных рав­ новероятных сигналов с одинаковым отношением сигнал/помеха (см. § 56), причем одновременное появле­ ние двух или более сигналов будем считать исключенным.

Ясно, что в силу равноправности всех сигналов [см. фор­ мулы (56.05) и (56.08)] все пороги в формулах (57.04)

и (57.05) следует брать одинаковыми, и тогда вероятности ложной тревоги и правильного обнаружения во всех кана­ лах приемника II типа будут одинаковы

F = F1 = . .. = FM, D - Dx = . . . = DM.

(57.06)

Вычислим

вероятность ложной тревоги в

приемнике

II типа. Случайные величины

в различных каналах не­

зависимы в

силу формулы (56.09), поэтому вероятность

правильного

необнаружения в Л4-канальном

приемнике

II типа по формуле умножения вероятностей независимых

событий равна

 

 

 

1-Р=(1— F)~,

(57.07)

и вероятность ложной тревоги

F' получается равной

 

F=l — (1- - F)M.

(57.08)

Вероятность пропуска сигнала приемником II типа, оче­

видно, равна

1 — D' = (l — D)(l — Л'И-1,

(57.09)

 

поскольку приемник II типа может обнаружить сигнал и за счет ложной тревоги в одном из М— 1 «пустых» кана­ лов. Вероятность правильного обнаружения равна

О' = 1 — (1 — D)(l --F)M~l.

(57.10)

'357

Формулы (57.08)

и (57.10) связывают вероятности F'

и D' в ЛТ-к'анальном

приемнике с вероятностями F и D

в отдельном канале, производящем обнаружение одного сигнала. Заметим, что F' и D' характеризуют лишь «обна­ руживающие» функции приемника II типа и не касаются его «измерительных» функций. При достаточно малом F мы можем пользоваться приближенными формулами

F’ = MF, D' = D.

(57.11)

Особенно простой результат эти формулы дают для

сигналов с неизвестной амплитудой и фазой

(§ 34). Дей­

ствительно, в силу формул (34.25) и (57.11)

мы имеем

/ р/ 1+р-

(57.12)

м)

где р есть отношение сигнал/помеха в каждом канале. Эта формула позволяет в явном виде найти значение р, необ­

ходимое для обеспечения заданных F' и D':

 

 

Inill -р77/- + Inin iviМ

 

 

(57.13)

 

 

р.= --------- i--------1.

 

 

 

 

In Д7-

 

 

 

 

 

В одноканальном

приемнике

обнаружения,

т.

е.

при

М—1, соответствующее значение р равно

[ср.

формулу

(34.24)]

 

 

 

 

 

 

 

 

 

In -рг

 

 

 

(57.14)

 

 

 

 

 

 

так что формулу (57.13) можно переписать в

виде

 

 

 

 

 

 

 

 

(57.15)

Мы видим, что

для

достижения

заданных

F'

и

D' в

AI-

канальном приемнике требуются тем большие значения

ц,

чем больше М,

причем необходимое приращение и.

пропор­

ционально In М.

358

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