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

книги из ГПНТБ / Марченко Б.Г. Метод стохастических интегральных представлений и его приложения в радиотехнике

.pdf
Скачиваний:
12
Добавлен:
24.10.2023
Размер:
9 Mб
Скачать

дал ему исчерпывающее объяснение, которое с несущественными доработками принято и теперь, В Шоттки [ 100]. Подобно дробовому эффекту могут быть построены модели других шумов, описывающих­ ся линейными случайными процессами. Некоторые из них будут рассмотрены ниже.

§ 2. Тепловые и полупроводниковые шумы. Реверберационные процессы и шумы торошения льдов

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

Наряду с дробовым эффектом большое распространение в ра­ диотехнике и физике получила модель тепловых шумов, которая описывает флуктуации тока и напряжения в электрических цепях, обусловленные хаотическим движением электронов в проводнике. Хотя изучением этого вида шумов физики занимаются начиная с конца прошлого века, основополагающие результаты, удобные для применения в теории связи, были получены только в 1927—1928 го­ дах Джонсоном [91] и Найквистом [99]. Эти результаты, начиная с 1928 г., неоднократно подтверждались в экспериментальных рабо­ тах многих авторов.

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

а) модель, которая исходит из общих принципов статистической механики, и

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

Так как целью настоящей главы является установление воз­ можности представления рассматриваемых здесь процессов с помо­ щью стохастического интеграла вида (1), то нет смысла подробно останавливаться на рассмотрении указанных выше моделей теплово­ го шума. Исходя из теоремы Найквиста 199] и близости теплового шума к броуновскому движению, что было установлено многими авторами [35] с использованием модели б), сразу приведем некото­ рые доводы в пользу возможности представления теплового шума в виде (1).

Согласно теореме Найквиста для электрических цепей [11, 73,

120

99], спектральная плотность шумовой э. д. с. в малом частотном

интервале df

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

выражением *

 

 

 

 

 

 

S(J)df = AkTRp (I/ 1) df,

 

(409)

где f — частота, гц; S (f)

спектральная плотность теплового шу-

ма; k — постоянная Больцмана ^1,38-10 23

 

j ; Т — абсолютная

температура,

°К;

R — со­

 

 

 

 

 

противление,

ом;

р (f)

 

 

 

 

 

фактор Планка:

 

 

 

 

 

 

= Л _

Р(/) =

—1

 

 

 

 

 

exp

hf_

 

 

 

 

 

kT

kT

— 1

 

 

 

 

 

где h =

6,62

10

(410)

 

 

 

 

 

34 дж X

 

 

 

 

 

X сек — постоянная План­

 

 

 

 

 

ка . Обоснование выражения

 

 

 

 

 

(409) содержится в [11].

 

 

р (/), определенной со­

На рис. 6 приводится график функции

гласно (410) для значений /, удовлетворяющих условию 0 <

< 4 .

При

нормальной комнатной температуре

в диапазоне

частот,

используемых в радиотехнике [73], имеем

А|П

1, поэтому в ука-

 

kT

занном диапазоне р

и выражение (409) в этом случае упро-

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

 

 

 

 

 

 

 

 

S(f)df ^AkTRdf.

 

 

(411)

В работе [35] при помощи метода моментов в предположении

линейной внешней нагрузки и с использованием модели б)

доказы­

вается, что тепловой шум имеет нормальное распределение, не за­ висящее от времени.

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

Рассмотрим следующие три случая.

Случай 1. Если, исходя из приближенной формулы (411), пред­ положить, что тепловой шум содержит в полосе частот (—оо, оо) постоянную спектральную плотность, т. е. распространить эту формулу на бесконечную область частот, то получим, что тепловой шум является «непрерывным гауссовым процессом с независимы­ ми значениями» или, более точно, представляет собой обобщенную производную [15, 29] от винеровского процесса:

Z(t) = o ^ P - ,

(412)

* В работе [ 11 ] в (409) вместо р (// |) приведено р (/), что не позволяет распро­ странить эту формулу на область / < 0.

121

где о2 = AkTR и w (t) — винеровский процесс. Так как g (t) яв­ ляется обобщенным случайным процессом, то выписать его плот­ ность распределения в явном виде в классе обычных функций не представляется возможным, поэтому под гауссовым процессом

подразумевается тот факт, что процесс

[w (t -f- Д^) — w(f)]

для

At > 0,

который с

уменьшением At

приближается к обобщенной

производной процесса ow (t),

имеет

плотность

распределения

 

Раі {х) =

 

ехр

 

 

Д *> 0.

 

Корреляционная

функция

этого

процесса

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

так:

М !

\w (t -f At) w (*)] [w(t +

s + A ) — w (t + s)J [ =

 

 

 

AP (A* — Is I) ПРИ

И < А /,

 

(413)

 

 

О

при \s\^> At.

 

Следовательно, его спектральная плотность существует и опре­ деляется следующим образом:

A t

 

 

,, ,

 

1 — cos 2nfAt

—/гл/s аа . . .

.

 

« = !

AP (A t - \s \) d s

=

----

-

°

— A t

 

 

 

 

 

 

 

При малых At эта

спектральная

плотность

 

близка к а 2, т. е.

при малых At функция 5д, (/) приближается к 5 (/), определенной согласно (411) при /£ (—оо, оо). Таким образом, хотя процесс £ (t), определенный согласно (412), и не существует в обычном смысле,

соответствующая ему

спектральная плотность, полученная

из

по­

следнего выражения

при

уменьшении

At, ведет

себя так,

как

будто бы

I (t) существует

и является

стационарным гауссовым

процессом

с постоянной спектральной плотностью о2.

(412),

В радиотехнике процесс g (t), определенный

согласно

обычно именуют «белым шумом». Для этого процесса Р (|£ {t)\~>A} = = 1 при любом положительном А . Физически осуществить процесс g (/) невозможно хотя бы потому, что он обладает бесконечной мощ­ ностью .

Таким образом, предположение о постоянстве спектральной плотности S if) на всем интервале частот/ £ (—оо, оо) для реального

физического теплового шума оказывается неприемлемым.

еди­

Случай 2. Предположим, что в (409) функция р |/| равна

нице в интервале(— /0, /0), f0> 0, равна нулю при |/ |l > f0, а

зна­

чения в точках 1/1 = /0

могут быть любыми, но конечными.

 

Не останавливаясь

на вопросе единственности линейного

ин-

122

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

V (f) — о \

sin 2я/о (t— т)

dw (т),

(414)

я (t— т)

где а 2 = 4kTR. Для этого достаточно рассмотреть спектральную плотность, соответствующую процессу (414) и сравнить ее с исход­ ной. Корреляционная функция процесса (414) согласно (298) опре­ деляется так:

со

sin 2nf0t

sin 2я/0 (t+ s)

 

sin 2n,f„s (415)

• (S) = ° 2 J

dt

nt

я {t-j- s)

 

 

Выражение для спектральной плотности процесса (414) с учетом [19], 3.741, (2), получаем в следующем виде:

 

а2

при

|/ |< / о ,

S (/) = а2 j sin^ f o e~i2nfs ds

-т р

при

|/1 = /о,

 

О

при

| / | > / 0.

Таким образом, центрированный гауссов процесс (414) имеет спектральную плотность, постоянную в полосе (—/0, /„) и равную нулю вне этой полосы, а следовательно, может быть принят в ка­ честве интегрального представления теплового шума. Однако такое представление все же остается очень идеализированным, так как его ядро является физически нереализуемым.

Итак, тепловой шум в рассматриваемом случае описывается слу­ чайным процессом, допускающим представление (1). Заметим, что процесс I (t), определенный согласно (412), также является в неко­ тором смысле линейным, так как формально может быть выражен в виде предела процесса (414) следующим образом:

СО

 

5(0 = а $ Пгп sin* f V - .$ -dw(T).

(416)

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

Случай 3. В наиболее общем случае, когда спектральная плот­ ность теплового шума определяется выражением (409), он может быть описан линейным процессом (1) по винеровскому:

CQ

 

о0(0 = ц ] ф0(г — x)dw(x),

(417)

123

где функция ф0 (і) определяется из уравнения (65), которое в рас­ сматриваемом случае принимает вид

 

 

и

С

p—tZKfsI

f I

 

 

Фо(*)Фо(* +

«)Я =

j

(h \f \ \

 

~ dt-

(418)

 

 

 

 

exp \~kT~) ~

1

 

Если дополнительно предположить, что функция ф0 (t)

является

физически реализуемой, то уравнение (418)

превращается

в урав­

нение типа Винера — Хопфа [62]:

 

( 2nkT

\

 

оо

 

со

 

 

 

 

 

Xcos (

т

XSI

(419)

J Фо (0 Фо (t +

S) dt =

j

-------- г г - ----- dx,

О

 

о

 

 

 

 

 

 

которое с учетом [17], 3.951, (5), можно записать еще и так:

оо

кТ /

h2

 

л2

,

2 2n2kTs

Г ... .. , . ,,

 

J *Ро (0 Фо (t + S) dt —

(

8n2fe2T2s2

2 cosech

h

 

о

 

 

 

 

 

 

 

 

Так как согласно (419) и с учетом [19],

3.411,

(1)

 

 

О

 

 

 

 

 

 

(421)

 

 

 

 

 

 

 

 

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

Выражение (421) с учетом представления (417) позволяет опре­ делить дисперсию теплового шума (аналогично [73])

h(0) = Mvl]{t) = » 1,8928 • 1 0 (дж ■сек ■ом), (422)

а выражение (298) с учетом (417) и (420) определяет его корреляци­ онную функцию

h(s) =

fiR

2n2k2T2R

cosech2

2n2kTs

(423)

2JI2S2

h

h

На рис. 7 приводится график нормированной корреляционной функции теплового шума, которая с учетом (422) и (423) определя­ ется следующим выражением:

ft(s)

3ft.2

2n2kTs

 

 

h

 

 

А(0)

4n*T2k2s2

2n4Ts

 

 

7,2907 • IQ"24

0,4114 • 10127’s

(424)

 

T2s2

sh (0,4114 • 1012rs)

 

 

При этом диапазон

изменения

s выбран из

условия

0 < s <

20 • ІО-12

 

 

интервале

корреля­

•< -----j,-----. Как видно из графика, на этом

ционная функция теплового шума убывает достаточно быстро-

124

Рассмотренные выше три случая позволяют сделать заключе­

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

процессов

для достаточно полного

описания реальных тепловых шумов.

Остановимся кратко

на простейшей линейной модели

полупро­

водниковых шумов. При анализе этих шумов обычно проводят раз­ личие между тепловым шумом, дробовым шумом и фликкер-шумом [11 ] Тепловой шум в полупроводниках обусловлен случайным дви­ жением носителей тока Он аналогичен тепловому шуму, возника­ ющему в проводниках, который, как было установлено выше, мож­ но описать интегральным представлением (1). Дробо­ вой шум в полупроводниках наиболее характерен, его спектр в области нижних частот равен постоянной ве­ личине независимо от ча­ стоты. По этому признаку он подобен дробовому шуму в электронных лампах, кото­ рый рассматривался ранее и, как было выяснено, мо­ жет быть описан интеграль­

ным представлением (1). Фликкер-шум характеризуется своеобраз­

ной

спектральной

плотностью

[11,

12,

59]

 

 

 

 

 

 

5 (со) = с2со~а, со 6(0,

оо),

 

(425)

где

а = 1

+

е, 8 — некоторая

постоянная

величина, 0

<; s <С 1,

со =

2яД

с — некоторая постоянная.

По

спектральной

 

характе­

ристике этот

шум

подобен фликкер-шуму

в электронных

лампах.

Фликкер-шум в полупроводниках значительно превышает уровень теплового и дробового шумов. При его анализе с помощью линей­ ной * теории обычно делают предположение о дискретности электрон­ ного заряда подобно тому, как это делалось в § 1 настоящей главы, т. е. предполагают, что группы зарядов (в § 1 рассматривались не группы, а отдельные электроны) появляются в случайные моменты времени и удовлетворяют предположениям а) и б) из § 1 настоящей главы. Следовательно, моменты их появления распределены по за­ кону (121) Так как эти заряды носят случайный характер и имеют одинаковые функции распределения, то, следуя рассуждениям, приведенным в § I, придем к представлению фликкер-шума в виде

(403)Однако теперь функция aNe. (t) будет представлять собой

импульс, связанный с зарядом группы, состоящей из N синхрон­

но

возникающих носителей, каждый из которых имеет заряд q =

 

со

=

j е(х) dx, aN — Nq, а не заряд, равный заряду электрона, как это

 

о

* Существует и нелинейная теория, развитая Моррисоном. Ее описание и ссылки на литературу даются в работе [11].

125

было в дробовом эффекте (т) — форма отдельного импульса). Следовательно, полупроводниковый фликкер-шум можно предста­ вить в виде

оо

 

£(0 = J e (f--T ) cfatj (т).

(426)

о

 

Теперь остается только определить е (т) таким

образом, если это

возможно, чтобы спектральная плотность процесса £ (t) описыва­ лась бы согласно (425).

Корреляционная функция процесса (426) с учетом (405) имеет следующий вид:

оо

 

h (s) —

J е (т) е (т + s) dxf

(427)

со

 

о

 

 

 

 

 

 

где а% = J хЧЗ~(х), 3 (х) — функция распределения зарядов

груп-

о

а е (т) удовлетворяет

условию физической

реали­

пы носителей,

зуемости .

 

 

 

 

Согласно теореме Винера — Хинчина [76] имеем

 

 

оо

оо

оо

 

S (со) =

j h (s)e~“°sds — Xa% j

J e (r) e(r -f s) elasdsdx.

(428)

— co

— oo 0

 

Функция

 

00

 

 

 

 

 

 

 

Sx(too) =

J e (s)е~ш ds

(429)

 

 

о

 

 

аналитическая в правой полуплоскости. С учетом (429) выражение (428) принимает вид

S (со) = ка% Sx (іа) Sj (— tco),

(430)

т. e. последнее выражение представляет собой факторизацию спект­ ральной плотности 5 (со) на прямой. і

Из (429)

 

ОО

■'*

e(0 = - s r

S (ісо)еш <ѣ.

(431)

Следовательно,

—00

 

00

 

 

 

5і(0) =

J е(т) dx = q,

(432)

 

о

 

поэтому должно выполняться условие

 

2

5(0)

(433)

ан

 

Если исходить из спектральной плотности, определенной вы­ ражением (425), то согласно (433) средняя интенсивность импуль­

126

сов носителей равна бесконечности. Кроме того,

функция

(425)

не интегрируема на интервале [0, оо).

Все это

не согласуется

с физикой явления. Значит, спектральную

плотность нельзя

опре­

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

аппроксимации. Исходя из

того, что область определения

(425)

не включает в себя точку

to = 0 и ее некоторую б-окрестность,

совершаем в (425)

замену переменной согласно соотношению

со2 =

= сох + б2. гДе 0

< б <С 1.

В результате получаем еще одну ап­

проксимацию спектральной плотности фликкер-шума в полупро­ водниках:

а

 

S (со) = с2 (б2 + со2)- ~ , (0 б [0, оо),

(434)

которая при малых, соответственно выбранных значениях б с на­ перед заданной степенью точности на любом конечном интервале из области со £ (0, оо) может аппроксимировать значения спектраль­ ной плотности, убывающие по гиперболическому закону, при этом постоянные с и п определены как в (425). Величину б при известных с и а можно определить экспериментально, путем измерения зна­ чения 5 (0) и вычислений по следующей формуле:

9

__1_

 

б = с а S

а (0),

(435)

полученной из (434). В случае спектра, определенного выражением (434), имеем

S ( < o ) = t a £ ------

 

 

 

 

(б +

 

«а) 2

(б — «а) 2

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с, =

 

 

— cqS

г (0)

 

 

 

 

 

 

 

V'-Ы{

 

 

 

 

 

 

 

 

При

0

с

учетом [6],

 

3.2, (3),

из последнего выражения

находим

 

 

 

 

 

 

 

 

 

 

 

 

 

e(t) =

-

се~ ы

1(0

 

 

gc exp { -[5 (0 ) с~2]

а

,

(436)

 

 

1—-

 

 

 

 

 

1—-

 

 

 

 

 

 

 

 

 

 

 

 

Ѵ Ч г ( у )

 

 

 

 

(0)Г

т)

 

 

 

где 1 (/) — единичная функция

 

Хевисайда:

 

 

 

 

 

 

 

 

1 (t) =

[0

при

^ с

0,

 

 

 

(437)

 

 

 

 

 

,

при

_

п

 

 

 

 

 

 

 

 

(1

I >

0,

 

 

 

 

которая

терпит

разрыв при t — 0 и в этой

точке

содержит

пре­

дел справа.

 

 

 

 

 

 

 

 

 

 

 

 

127

Нетрудно видеть,

что функция е (t), определенная согласно

(436), отличается от

функции Лагерра нулевого порядка только

постоянным множителем, поэтому для нее выполняется условие (432).

Следовательно,

в рассматриваемом случае форма импульса

е (t) определяется

спектральной плотностью S (со) с точностью до

постоянного множителя независимо от вида cf (x).B свою очередь спектральная плотность не накладывает никаких существенных ограничений на распределение адт-суммарного заряда группы, за исключением условия (433), благодаря которому нам пришлось пе­ ресмотреть аппроксимацию спектральной плотности и придти к ее новому виду (434). Заметим, что вопрос, какая из аппроксимаций (425) или (434) лучше, подлежит экспериментальной проверке. Таким образом, исходя из физической линейной модели шума в по­ лупроводниках, приходим к заключению, что он также может быть описан интегральным представлением (1).

Рассмотрим кратко линейную модель шумов, обусловленных рассеянием звука на неоднородностях морской среды, именуе­ мых реверберационными шумами. Теоретические результаты пер­ вой главы позволяют описать при помощи характеристических функционалов поле реверберационной помехи. Рассмотрим вначале модель и свойства реверберационной помехи в некоторой точке шумового поля. Следуя работе [63], будем различать объемную ре­ верберацию, реверберацию от слоя и граничную реверберацию, а основное внимание сосредоточим на рассмотрении модели объем­ ной реверберации, так как остальные виды не имеют принципиаль­ ных отличий. На практике часто приходится иметь дело с линейной комбинацией сразу двух или трех видов реверберации [63].

Будем считать, что рассеиватели в морской среде расположе­ ны дискретно. Эффективные сечения рассеивателей практически оказываются намного меньше расстояний, занимаемых излучае­ мыми сигналами в пространстве. Если предположить, что рас­ сеиватели в среде расположены статистически независимо, для рассматриваемой достаточно большой области среды средняя плот­ ность рассеивателей постоянна, то количество элементарных рассеян­ ных сигналов, приходящих в точку приема за время Ат, т. е. попав­ ших в интервал (т, т -f- Дт), подчиняется распределению Пуассона (121). Это можно обосновать аналогично обоснованию выражения (121), проведенному для случая дробового эффекта.

Обозначим через ak случайные амплитуды, а через tk — слу­ чайные моменты возникновения рассеянных сигналов. Если из­ лучаемый сигнал описывается функцией s (t) и каждый элементар­ ный рассеянный сигнал s (tk, t), пришедший в момент времени tk, зависит от состояния рассеивателя в момент прохождения через него

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

фикси­

рованного момента времени t может быть записан в виде

 

£ ( 0 = Ü a*s(**,g.

(438)

— ОС

 

128

При помощи определенного в § 2 гл. I процесса щ (т), прира­ щения которого связаны с ak соотношением (402), реверберационный шум можно записать в виде

t

со

 

 

Ш = Ш = i s (т, 0 djti (т) =

j

s(x, ОЛтДт),

(439)

— со

— со

 

 

где предполагается, что s (т, f) 0 при т >

t.

 

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

Как уже упоминалось в § 1 настоящей главы, при к ->■ оо про­ цесс, полученный в результате центрирования и нормирования про­ цесса (т), стремится к винеровскому процессу, поэтому и процесс (439) после центрирования и нормирования также стремит­

ся к гауссовому в общем

нестационарному процессу,

который в

этом случае можно записать так:

 

 

СО

 

V (t) =

^ s (т, t) dw (т).

(440)

 

— со

 

Обоснование этого выражения можно построить аналогично обос­ нованию соотношения (408).

Ядро

процесса (439) может

быть найдено из уравнения

 

ОО

 

 

 

[ s (т, 4 )s(r,

t%)dx — h(tv 4),

(441)

 

— СО

 

 

где h (4,

4) — корреляционная

функция процесса (439),

которая

может быть определена из теоретических предпосылок или экспе­ риментально.

Если в (406) и (407) функции е (4 — т) и е (і — т) везде заме­ нить соответственно на s (т, tk) и s (т, і), то получим многомерную характеристическую функцию и характеристический функционал реверберационного шума для некоторой точки шумового поля.

Для описания «-компонентного вектора реверберационной по­

мехи

(^ (t),

t,2

(t),

..., £„(4 ), воспользовавшись теоремой 6, можно

сразу

привести

выражение характеристического функционала,

который с учетом (127) можно записать в виде

 

 

 

 

 

(gif Sä,

,

gn) =

 

 

 

ОО

ОО

п

 

 

 

ехр

j

J

[ехр(гл: 2

S * ( T ,

f)dgk {t)) 1 ]d3r (x)dx). (442)

 

 

 

— оо — oo

f t = l

 

 

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

9 3 — 6 6 2

129

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