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

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

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

В этих переменных мы имеем также

J] uhxh.

У] ав^ёхи. — У] Vft’ ~

aA/ix/i’

h

d,h

k

h

(59.04)

и кратный интеграл (59.01) распадается на произведение Н простых интегралов

P(xlt..., х„) =

00

= (2^" JeXPUZS^--2- XPft^}dZ1’ ’ 'dtH —

—оо k k

(59.05)

«=1 —оо

каждый из которых легко вычисляется

(59.06)

(59.07)

Из теооии пэеобразований квадратичных форм известно, что определитель, составленный из их коэффициентов,

остается

при ортогональных преобразованиях (59.02) инва­

риантом,

поэтому

 

 

 

 

De,iM=Пр.= . •• • ?«■

<59-08)

 

 

k

 

 

Далее по формуле (59.04) получаем

 

S 7T=S

где

= S 7Г Vkh’

(5909)

й

g,h

 

Й

369

24—483

 

 

 

в то время как формула (59.03) дает

^=SP,V/A'

(59’10)

/

 

 

Используя соотношения, которым

должны

удовлетворять

коэффициенты а^Л ортогонального

преобразования (59.02)

 

2ал=2а^ал/==5г/1’

 

 

 

(59,П)

 

k

 

 

i

 

 

 

 

 

 

 

 

мы легко

получаем формулу

 

 

 

 

 

 

 

 

R'gl^Q-lh = V7j — аig. а.,11 а.,Ма..kh

=

рй 8ik

ig kh = 8gh'.,

 

SI

j.k,l

 

 

 

j,k

 

 

 

 

(59.12)

 

 

 

 

 

 

 

 

 

 

 

показывающую,

что

Q h — элементы

матрицы,

обратной

матрице

ЦТ? ft||.

Поэтому

формула (59.07)

в развернутом

виде гласит

 

 

 

 

 

 

 

 

 

 

 

р(хг,..

хн)—1......... — ехр /

5-V Qahxaxh\ .

 

И

/(2K)"Det||M|

П

 

4*^£^

 

(59.13)

 

 

 

 

 

 

 

 

 

 

 

При рассмотрении

комплексных

случайных

и

величин

нужно иметь в виду, что в силу

формул

 

(58.14)

(58.18)

мы имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

dxhdyK = -±- dzhdz* h , duhdvh =

dwh*dw h

,

(59.14)

поэтому,

определяя плотность

вероятности для комплекс­

ных величин так,

что

p{zr,

z* ,

. .., zH,

*z H)dzl*dz ... dzH*dz H

есть вероятность

того, что

эти

величины

попадают в

со­

ответствующий

элемент многомерного

пространства,

мы

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

р (гх, z\ ,. . ., zH, z* H) = Я,‘'ИеХР{~^Х(^г/г +

h

+ wA) ~ 4"J] Rghwg wh' } *dW'dw . . . dwHdwH.

(59.15)

370

С помощью унитарного преобразования

=

wh = Yi а'^8’

(59,16)

л

g

 

коэффициенты которого по определению удовлетворяют соотношениям

 

L

a, .а*=

V а,

а*

kh

=8.,

 

 

(59.17)

 

gk

hk

/ j kg

 

 

gh'

 

 

'

 

 

k

 

 

k

 

 

 

 

 

 

 

 

 

 

 

мы можем

привести,

эрмитову

матрицу

к

главным осям

 

J]

 

Ч

pJV

 

(рй>°)’

 

<59-18)

причем

g,h

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

k

 

h.

Zh=

k

 

 

 

 

 

 

h

 

 

h

 

 

 

 

 

 

 

 

(59.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вводя вещественные переменные по формулам

 

 

 

 

tk = rk + isk^ 4 = ^ + 4’

 

 

 

(59.20)

мы будем иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dtkdik — 2drttdsk, \t^ = rk2-\-s2k,

 

 

 

 

Ф(<?4+<.О = =Л + ’1Л.

 

 

 

(59.21)

откуда, пользуясь формулами (59.06), получаем

 

 

Jjexp{ 2

 

 

 

4-рЖ|2}^ dt: =

 

 

 

 

 

 

 

 

 

2

ГЯ 1

«1

(

R

R

 

 

= 2

J е

 

 

drk J

е ~trlksk

2"

 

PASA

,

 

 

 

—оо

 

 

—00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*2р

 

4

 

 

l^l"

 

 

 

(59.22)

 

 

е

 

-е 29k

,

 

 

— —

 

 

 

 

24*

 

 

 

 

 

 

 

 

 

 

 

 

 

371

так что

Р(2‘’21......= (4я)"П

еХр I

“2j "tr/ ’

<59-23)

 

 

 

 

 

 

 

 

k

 

 

 

k

 

 

 

Соотношение (59.08) справедливо

и

для

унитарных преоб­

разований, а

 

 

 

 

 

 

 

 

 

 

 

 

 

L ~7Г=2

 

 

где

 

= S TZ

(59'24)

к

 

8,Л

 

 

 

 

 

 

k

 

 

 

 

Поскольку

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2 Pja/g

 

и S

 

 

 

(59-25)

 

 

 

 

I

 

 

 

I

 

 

 

 

 

 

то матрица ||Qgft|| обратна

матрице

 

и

плотность

веро­

ятности

для Н

комплексных

нормальных

величин

 

равна

 

 

 

 

Р (^1 ’ ^1

>

■ • • 5

Zfl1

 

 

 

 

 

 

 

_ ________ _ exo J-----L-У

8h

!•

(59.26)

 

dW7 Det ||Р ;ц

 

)

 

2

Zj

8 h /

 

 

 

 

 

 

 

 

 

 

 

g,h

 

 

 

 

Если

по формулам

(58.14)

и

(59.14)

вернуться

 

к ве­

щественным

случайным

величинам

хг, .

.., хн и ylt ...,

ун и представить

эрмитову

матрицу

Qgh

в

виде

 

 

 

 

 

 

Q

— Q(D_LiQ(2>

 

 

/р-р

 

 

 

 

^■gh

 

 

i

^-gh

 

 

 

 

 

 

 

 

 

Q(2)=_Q(2)

 

 

(

}

 

 

 

^■gh

^hg'

^-gh

 

^-hg’

 

 

 

то плотность

вероятности будет

равна

 

 

 

 

 

 

 

 

 

 

 

У1, ..., уи)=

 

 

 

~ (2K)wDet ||7?

Д| еХР

 

 

 

 

(XgXfl +

 

 

 

 

 

 

 

 

 

 

S-h

 

 

у

 

'

 

 

 

 

+

 

 

 

 

 

 

 

 

(59.28)

Заметим,

что

если

мы имеем Н нормальных случайных .

величин

 

 

..., хн,

то

введением

Н дополнительных

случайных величин У1,

..., уи,

удовлетворяющих

усло­

виям (58.12),

можно

всегда

прийти к комплексным нор-

372

мальным величинам

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

выше типа.

При этом

коэффициенты R^h

можно

выбирать различным

образом.

Если, например, считать

=

то

величины (58.22)

бу­

дут вещественными.

В §

43 мы имели дело со

случаем,

когда выбор комплексных чисел

R h

подсказывался общей

формулой для корреляционной функции, полученной

ниже

в § 69.

 

 

 

 

 

 

Формула (59.13) дает //-мерное распределение Гаусса

для вещественных случайных величин, формула

(59.26) —

для комплексных.

 

 

 

 

 

 

С точки зрения теории вероятностей значение гауссовых случайных величин определяется тем, что величины, каж­ дая из которых является суммой достаточно большого чи­ сла независимых (или слабо зависимых) случайных вели­ чин, имеют распределение, близкое к гауссовому. В теории вероятностей исследуются условия, при которых гауссово распределение (59.13) является предельным в этом смысле.

§60. РАСПРЕДЕЛЕНИЕ ВЕРОЯТНОСТЕЙ ДЛЯ ОГИБАЮЩИХ

Вгл. VII мы пользовались для комплексных случайных величин несколько иными обозначениями. Мы исходили

из вещественных случайных величин

mx = Gxcosyx, t»x==Gxsin yx (х= 1, ..., L) (60.01)

и считали их гауссовыми. При этом естественно предполо­ жить, что одномерная плотность вероятности />(ух) и двух­

мерная

плотность

/’('fx, Yx) удовлетворяют соотношениям

/?('L) = /?(Yz + a)-

МГЛА) = Ж. + Мх + Н (60.02)

где & — произвольный

фазовый угол. Эти соотношения

означают, что при

смещении всех фаз ух на одинаковый

угол

мы приходим к

новым фазам, обладающим той же

вероятностью. Из соотношений (60.02) вытекают тож­

дества

cosyx = 0, sinyx = 0,

cos(yx4-yx) = 0, sin (Yx+Yx”) = °-

. (60.03)

cos yx cos yx = sin yx sin yx, cos ух sin Yx= — sin yx cos yx,

373

откуда получаем

и = 0, v = 0,

(60.04)

если сделать дополнительное предположение, что величины

Gx

и

Gx (амплитуды

или

огибающие) независимы

от

фаз

Yx

и ух. Если огибающие Gx нормировать по формуле

 

 

 

 

 

 

G2 = 2,

 

 

(60.05)

то

 

будут коэффициентами корреляции, поскольку

 

 

 

 

 

 

г‘х’=1.

 

 

(60.06)

 

Другие свойства коэффициентов гхХ} и гхХ

выражаются

формулами

 

 

 

 

 

 

(60.07)

 

 

 

 

Г^ = -Г^, г{2) = 0.

 

Если ввести комплексные случайные величины

 

 

 

 

 

wx = zzx— zwx —Gxe~!Tx,

ay’ = zzx -ф- ivx — GG'’-,

(60.08)

то

их моменты будут равны

 

 

 

 

 

 

 

 

 

w =*ау

=0,

 

 

 

 

 

ww,—w*

0,

 

(60.09)

 

 

 

w

х

= 2г™,

 

 

 

 

 

 

 

 

Л

хЛ’

 

 

 

 

где комплексный коэффициент

корреляции гхХ

равен

 

 

 

 

 

 

 

>■"=&.

 

(60.Ю)

 

При L—1 эрмитова матрица ||г”|| и обратная

ей

мат­

рица

||<7™|| сводятся

к единице

и плотность

вероятности

(59.13) принимает вид

 

 

 

л

 

 

 

 

 

 

 

 

 

(60.11)

 

 

P(«1^i) = iexp.{-4-(«?+^)|=^e

2.

При переходе от декартовых координат

к полярным

координатам Gn Yj

нужно учесть соотношение

 

 

 

 

 

 

dii.dv^G.dG^,,

 

(60.12)

374

откуда

 

 

 

 

 

 

=

e

2,

 

 

(60.13)

так что одномерные

плотности

вероятности

для

каждой

величины G\ и

ух

равны

 

 

 

 

 

 

 

 

 

/’(Yi)—2^

(О’Су1<2г),

|

 

 

 

 

 

 

 

 

0?

 

 

 

 

 

(60.14)

 

p(Gl)=Gle“y (0<G1<oo). )

 

 

При L = 2,

обозначая г'" = ре ге,

мы имеем

 

ре

 

 

ll^zxh

1 — р2

1

- pe~k

(60.15)

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

-peis

 

 

откуда

 

 

 

 

p(ultu2; v1( v8) =

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

Of +G2 —2pG‘ lG2 cos (7i —

+ ь)

(60.16)

/п V 2 z |

вхр

 

 

 

 

2(1 ~P2)

 

 

(2л)2(1 —р2)

 

Г

 

 

 

 

 

 

 

и при переходе к

 

полярным координатам получаем

 

 

 

 

 

/?(Gj, G2,

y2)

 

 

 

 

OiG2

 

 

 

 

О2 + О|

2pGiG2

cos (и — 7г + t)

(60.17)

(2я)г (1 —/>2)еХр

 

 

 

 

 

2(1-/О

 

 

 

 

 

 

 

 

 

 

 

 

 

Чтобы найти

двухмерную плотность

вероятности

p{Gx,G^y

необходимо

проинтегрировать

выражение (60.17) по ух

и у2. При интегрировании по Yi

нужно

использовать тож­

дество

 

 

 

 

 

 

 

 

 

 

 

 

2тс

 

 

 

 

 

 

 

 

 

 

 

(б0 J8)

A. J exp ^ c°s/b-72 + e)^Ti = Jo (pGjG^ ,

о

 

 

 

 

 

 

 

 

 

 

 

 

а последующее интегрирование

по

у2 дает

множитель 2х.

Окончательно получаем

 

 

 

 

 

 

 

/?(Gi, G2) =

 

O1O2

exp

 

Of + G2 К /pGiG2\

(60.19)

 

 

2(1—Р2))0 (J — p2J’

 

 

 

 

1 —р2

 

 

 

375

Если G;

и G2 нормированы иначе, например,

 

 

Gf = Gf=2a2,

 

(60.20)

так что

 

 

 

 

 

 

 

■w\w2~2^po~l\

 

(60.21)

то формула (60.19), как легко

показать,

приобретает вид

n (G

G) =__ _____ ехо /____ \f ( PG'G*

\

Р

 

с‘(1 - />2)ех9 \

2а2 (1 - *р)

/

°^2(1 - Р2})'

 

 

 

 

 

(60.22)

Заметим

в заключение, что

данное

выше обоснование

формул (60.04) является недостаточно строгим. Ввиду то­ го, что из этих формул вытекают многочисленные и важ­ ные следствия, следует отметить, что строгое изложение

вопроса должно основываться на определении огибающей

и фазы случайного процесса, данном В. И. Бунимовичем (см. его книгу). Обычно же эти величины вводятся с ма­ тематической точки зрения недостаточно корректно, что не

позволяет дать убедительного доказательства формул

(60.04).

§ 61. ТЕПЛОВЫЕ И ДРОБОВЫЕ ШУМЫ КАК ПРИМЕРЫ НОРМАЛЬНЫХ СЛУЧАЙНЫХ ПРОЦЕССОВ

Случайные величины возникают при рассмотрении «вы­ борок» случайных процессов, т. е. значений различных слу­ чайных функций времени в дискретные моменты. При этом

нормальными, или гауссовыми, случайными процессами на­

зываются такие процессы, выборки из которых имеют рас­ пределение Гаусса (см. § 59), или, что то же самое, ха­ рактеристическую функцию вида (58.07).

В теории оптимальных приемников обычно предполага­ ют, что помеха является стационарным случайным процес­ сом нормального типа. Часто делают еще одно дополни­

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

но случайный процесс»), может быть представлена, как беспорядочное наложение возмущений с нулевой длитель­ ностью (ср. § 12 и 17). Существенно отметить, что соб­ ственные шумы приемника — тепловые и дробовые, которые будут подробно рассмотрены ниже,—дают основа­ ние для такой идеализации. Однако представление о по-

■376

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

альности и, что

особенно важно,

не позволяет включить

в

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

интересный

вид

помех, встречающихся

в

радиолокации, а именно помех, обусловленных хаотиче­

скими отражениями (см. гл.

XI).

 

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

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

тате флюктуаций возникает случайная

электродвижущая

сила <§(/), спектральная интенсивность которой равна

5^(03)= 2W,

(61.01)

где k—постоянная Больцмана, R—сопротивление. Это со­

отношение применимо к любому элементу электрической цепи, причем R в общем случае определяется по формуле

/? = ReZ(«>),

(61.02)

где Z (®) — импеданс данного элемента

при частоте о>.

Мы видим, что для чисто омического сопротивления, не зависящего от частоты, флюктуационная электродвижущая сила $('/) является белым шумом, поскольку ее спектраль­

ная интенсивность является константой. Фактически, ра­ зумеется, при достаточно высоких частотах постоянство Sg (со) всегда нарушается вследствие остаточных емкостей

или индуктивностей, скин-эффекта или квантовых явлений.

Нужно еще учесть, что при росте частоты квазистационар-

ное рассмотрение цепи и вместе с ним формула (61.01) те­

ряют свой смысл.

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

Дробовой эффект в электронных лампах возникает как результат дискретной структуры электронного потока. Ес-

377’

ли, например, мы имеем диод, к которому приложено до­ статочно большое анодное напряжение, так что все испус­ каемые катодом электроны поступают на анод (рис. 56,а), то вследствие флюктуа­ ций катодной эмиссии на средний анодный ток J

накладывается случай­ ный «дробовой» ток /(О со спектральной интен­

сивностью

S^^eJ, (61.03)

Рис. 56. Электронные потоки в диоде:

а) обычный диод—электроны движутся от катода к аноду; б) .равновесный' диод—элек­ троны движутся навстречу друг другу и полный анодный ток равен нулю.

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

/(/), возникающий бла­ годаря дробовому эффек­ ту, также можно считать

белым шумом, поскольку его спектральная интенсивность не зависит от со. Фактически, если под /(/) понимать ток во внешней цепи диода, формула (61.03) применима при условии

«ГСО<1,

(61.04)

где То есть время пролета электрона через диод или дли­ тельность импульса во внешней цепи, появляющегося при прохождении одного электрона от катода к аноду. При ®То>1 происходит резкое спадание спектральной интен­ сивности дробовых шумов.

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

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

ный заряд уменьшает дробовой эффект, так что для диода

формула (61.03) видоизменяется

следующим образом:

S{(M) — FaeJ,

(61.05)

где F есть так называемый коэффициент депрессии дробо­ вого эффекта, a J—анодный ток.

Отметим, что электронная лампа является неравновес­

378

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