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

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

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

переход от компонент

вектора V (/, s) к

компонентам

вектора

V (t, s), можно записать в следующем виде:

 

 

V (t, s) =

__і_

s),

 

М

2 (/,

 

откуда

 

 

 

 

V (t, s) = А (t, s) M ~ (t, s) V (t, s),

s Ф 0,

(384)

где элементы матрицы* А (t,

s) и M (t, s) вычисляются по рекуррент­

ным формулам (382) и (383). Эти рекуррентные формулы,

а также

формулы (368), оказываются весьма удобными при использовании ЭВМ для вычисления численных значений элементов искомых мат­

риц преобразования

компонент.

 

 

 

 

_

 

 

 

Выпишем в явном виде элементы матриц А (t,

s) и М (/, s) при а =

= 2 и s Ф 0. Согласно (383) с учетом (382) имеем

 

 

 

 

 

 

^21(ß>®) —

(<>(+ S)

Ö 3 1 (*. S) =

fen(M)

 

 

 

 

fen (t, 0

feu(U)

 

 

 

 

 

 

 

 

 

 

 

 

 

Й41(/>

s) — feil (6 ^ ~f~ S)

 

 

 

 

 

 

 

 

 

 

 

 

feil ß

0

 

 

 

 

 

 

 

 

An (t ■о к

1 (6 < +

8) - feil t +

S) h12(t, t)

 

 

 

 

 

 

feil (*. 0 feil (* +

S, t +

s) — Af, (f,

t +

s)

 

 

^

S4

fen (t, t) hl2 (t + s , t + s ) — fen (/, t +

s) hl2 (t, I + s)

 

Ы

'

 

 

Au(<.0AuP + s.^ + s)-A ? ,ft^ + s)

 

(385)

Ö43 (t,

S) =

h-i (t,

S) fe22 {U t + s) ■

fel2 (t<0 feil (<■* + S)

 

 

 

feu (/» 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ß32(/> s) Ö42(^i s) h22 s)

 

 

 

 

 

 

Согласно (382) с учетом (383) имеем

 

 

 

 

 

 

 

fen (/> s) — Лц (t, f),

h2

2 (/> s) feu (Л 0 fen (* + *

f + s) - А?, (/, f +

s)

 

 

 

 

 

 

 

feu

t)

 

 

 

 

fe33 (*. s) =

{[/ІЦ (*. 0 fe22(*, 0

— fel2(/, /)] [fell (/, f) fell (/ +

S ./ + S ) -

— /і2іі(/, / + s)[ — [feu(/, 0 fe2i (/> / + s) — fen (^> / +

s) h12 (t,

t)\2) X

 

 

 

 

 

 

 

 

 

 

 

 

 

(386)

К fen' (t, t) [hn (t, t) hu (t —j—s, / -j- s) — fen (/,

t +

s)]

,

 

__

 

 

 

 

 

3

 

 

 

 

 

 

fe«4 (*> s) =

fe22 (* +

s. t +

s) — 2 ] of* (/, s) hkk (t,

s).

 

 

* Фактически некоторые элементы матриц А (/, s) и М (<, s) от s могут и не зависеть, однако для общности всегда будем пользоваться обозначениями

ац (t, s) и Кц (/, s).

110

Наряду с этими формулами

часто бывает удобно пользовать­

ся и такими соотношениями:

 

 

 

 

 

^ 2 2 (^» 0 “

аМ (ßl

/і]д (tj S) “f- Ö32 (/, S) /і22 (^l 5) -|—/igg(/,

S),

 

 

 

 

2

_

 

_

 

 

 

/i22 (* +

s, t + s) = 2

(t, s)

(/, s) + hgg(/, s),

 

(387)

 

 

2

k=\

 

 

 

 

 

 

 

 

 

 

 

 

 

h22 (t, t + s) = 2

a3k (t, s) aAk (t, s) hkk (t, s) +

a43 (г', s) /Г38 (/,

s).

 

 

fe=sl

 

 

 

 

 

у (/),

Если

существует

некоторый

нелинейный

функционал

определенный согласно (370) с

зависимыми

компонентами

(371),

то в результате замены компонент согласно (384) получим

 

У У) У (t + s) = Gt,s [ѵг(t, s).......... v2n(it, s)] =

Gt [xt (t), ...

 

 

........xn(01 GH-S[XJ.(t +

s), . . . ,

xn(t + s)].

 

(388)

Из последнего выражения согласно (247) получим

 

 

$ (S, t) =

М \у (t) g{t +

s)] = МGt,s [ѵг {t, s),

. . . ,

v2n (t, s)] =

c00...o (t, s),

 

 

 

 

 

 

 

 

(389)

где coo.,.0 (t, s) — коэффициент

Фурье — Эрмита для

функции

Gt,s(xb х2, ..., х2п), определяемый согласно (356).

_Таким образом, при помощи преобразования компонент вида (384) получены все необходимые выражения для вычисления сме­

щенной корреляционной функции

процесса у

(і), определенного

согласно (370). Применение этих

результатов

будет рассмотрено

в § 5 и 6 гл. III.

 

 

Г л а в а I I I ПРИЛОЖЕНИЯ МЕТОДА ИНТЕГРАЛЬНЫХ ПРЕДСТАВЛЕНИЙ В РАДИОТЕХНИКЕ И ТЕОРИИ СВЯЗИ

§1. Основные положения анализа шумов. Дробовой эффект

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

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

1) задание случайных процессов при помощи характеристи­ ческих функций или характеристических функционалов и, в част­

ности, при помощи соответствующих

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

2) задание случайных процессов

при помощи интегральных

идругих операторов;

3)задание случайных процессов при помощи дифференциальных уравнений.

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

1)свойства функций распределения или характеристических функций значений процесса в отдельные моменты времени;

2)спектральные и корреляционные свойства.

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

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

112

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

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

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

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

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

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

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

8 3 — 6 6 2

113

При рассмотрении многих шумовых процессов, обусловленных дискретностью электронного заряда, авторы [8, 21, 59, 60, 67]

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

1.За малый промежуток времени (т, т + Ат) вероятность по­

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

АДт +

о (Ат),

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

равна

о (Ат).

Этого предположения достаточно, чтобы сделать вывод о том, что вероятность Рк (т) появления ровно k импульсов за время т определяется выражением (121), в котором к теперь равна средней плотности числа импульсов за единицу времени.

В работах [21, 671 авторы, отправляясь от предположения а), дают более подробные выкладки, приводящие к распределению (121) для появления k импульсов за время т в случае дробового эф­ фекта, хотя в 167] при этом рассмотрен некоторый конечный ин­

тервал т £ 10, Т\, вместо бесконечного интервала т £ [0, оо).

 

Обозначим через л (т) случайную функцию такую,

что я (0) = 0

и я (т +

Ат) — я (т) равно числу появившихся импульсов за

про­

межуток

времени от т до т -f- Ат. Предполагаем,

что для

я (т)

имеет место сформулированное выше предположение 1.

 

2.

На

непересекающихся интервалах времени (т^ т2) и (т3, т4),

т, <

т2 <

т3 < т4, величины я (т2) — я (т4) и я (т4) — я (т3)

вза­

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

Из предположений 1

и 2 следует, что я (т) является

процессом

Пуассона с параметром

к [26], поэтому

 

 

М [я (т -f- Ат) — я (т)[ =

ААт,

 

М [я (т -}- Ат) — я (т)]2 = А,Ат

о (Ат).

(390)

Центральное место в теории шумов занимает теорема о наложе­ нии случайных возмущений [67], частный случай которой носит название теоремы Кемпбелла [26]. Остановимся кратко на содержа­ нии этой теоремы.

Пусть каждое событие (например, вылет электрона в дробовом эффекте) имеет некоторую интенсивность а и вызывает эффект, величина которого через время t после момента возникновения этого события равна е (t), причем

ооос

а = j e(t)dt< oo

и j e 2(0 < # < сю,

(391)

0

о

 

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

Если выходной контур рассматриваемой системы таков, что

llé

эффекты, вызываемые отдельными событиями (вылетом электронов), складываются линейно, и если предположить, следуя работе [59], что продолжительность отдельного события (время выхода электро­ на) является очень короткой по сравнению с длительностью е (/), т. е. перекрытием исходных импульсов можно пренебречь, то пол­ ный эффект в момент t от всех предшествующих событий равен

 

оо

 

W )

= 2 e ( t - t k),

(392)

 

k = £ ---- ОО

 

где tk — всевозможные

предшествующие t случайные

моменты

появления событий.

 

 

Если воспользоваться рассмотренным выше пуассоновским про­ цессом я (т), то сумму (392) можно записать в виде

t

Ш = 1 (0 = J e(t — т) dn (т).

(393)

— оо

 

При этом точки роста реализации процесса я (т) совпадают с точ­ ками tk.

Если дополнительно предположить,

что

 

e(t) — 0 при ( <

0,

(394)

то выражение (393) можно записать в виде

 

ОС

 

 

£ (/)= J e(t — x)dn(x),

(395)

— ОО

 

 

где £ (0 — стационарный и линейный в узком смысле процесс, по­ рожденный пуассоновским процессом и подробно рассмотренный в § 2 гл. I .

Сформулируем теперь в терминах математических ожиданий тео­ рему (в [67 ] эта теорема сформулирована в терминах временных сред­ них), которая впервые была сформулирована и доказана Кемп­ беллом в работах [87, 881.

Теорема о наложении случайных возмущений (теоре­ ма Кемпбелла). Среднее значение процесса £ (/), выраженного соглас­ но (395), определяется выражением

оо

 

Щ (t) = К j' е (т) dx,

(396)

— оо

 

средний квадрат флуктуаций относительно среднего значения оп­ ределяется согласно

оо

М К (0 - щ>(012 = а. J е2 (Т)dx

(397)

— ОС

 

идля процесса £(/) имеет место гипотеза об эргодичности.

До к а з а т е л ь с т в о . Эта теорема является следствием тео­

ремы 3, приведенной в гл. I. Выражения (396) и (397) сразу следуют из выражения (118) как частный случай.

8*

115.

В пользу подтверждения гипотезы об

эргодичности процесса

£ (t) можно сказать следующее. В работе

[67] доказательство тео­

ремы о наложении (последнее не претендует на математическую

строгость,

так как рассматриваются

только значения і £ [О, Т ])

проводится

в терминах временных

средних, а

так как (396) и

(397) совпадают с соответствующими

временными

средними, полу­

ченными в [671, то это свидетельствует о выполнении эргодической гипотезы, по крайней мере, в рамках первых двух моментов

В дополнение к теореме о наложении приведем еще выражение для корреляционной функции процесса £ ((), которое может быть

получено непосредственно из теоремы 3

или

из выражения (118):

 

 

 

 

оо

 

 

 

 

М {[£(/)- М £ (01 [£(* + s) - М £ (/)]}= Я J

е(х)е(х +

s)dx.

(398)

 

 

 

 

— оо

 

 

 

 

Одномерная

характеристическая

функция

процесса £ (/)

по­

лучается в рассматриваемом случае

из выражения (116) при п = 1

в виде

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

MeiuUt) = exp (Я J [eiueW— 1] dx).

 

(399)

 

—со

 

 

 

 

 

 

 

Одномерная

плотность распределения процесса £ (0

получается

в результате применения к (399) преобразования Фурье

 

 

 

00

оо

 

 

 

 

 

 

р(х) — -2^г j exp [ — шх + Я

j

[e'ue(r) IJdxjdu.

(400)

 

—оо

—-со

 

 

 

 

 

Выражение (400) полностью совпадает с

выражением,

при­

веденным в [671 под номером 1.4—7 и

полученным там

иным

спо­

собом

 

 

 

 

 

 

 

 

Частный, но важный для приложений случай получения чис­

ловых значений характеристической

функции

(399)

и

плотности

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

(400), когда е (т) = ßc~ßT

при

т > 0,

рассмотрен в

§ 3 гл I Графики этих функций приводятся на рис. 3 и 4. Из рис. 4

можно заключить, что с уменьшением

ß кривые плотности распре­

деления приближаются к гауссовым.

 

 

 

 

 

 

 

Пользуясь

результатами гл. I можно

без

особых

трудностей

записать любую конечномерную характеристическую функцию про­ цесса £ (t) или его характеристический функционал, а следователь­ но, путем преобразования Фурье найти конечномерные распределе­ ния процесса £ (t). Это значит, что процесс полностью задается сво­ им ядром е (х).

В [671 дается выражение для ряда Эджворта плотности (400), из которого можно оценить быстроту сходимости плотности (400) к нормальному закону при Я -> оо.

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

одинаковы и равны а все интенсивности рассматриваемых событий, т. е. рассмотрен дробовой эффект в режиме насыщения. Однако на практике возможен случай, когда в дробовом эффекте нельзя не учитывать различие начальных скоростей отдельных электронов. Рассмотрим этот случай подробнее. При этом сохраним все основ­ ные предпосылки, сделанные в начале настоящего параграфа, не­ изменными, за исключением того, что теперь для каждого события ak будут случайными независимыми и одинаково распределенными величинами. В таком случае суммарный эффект в некоторый момент времени t может быть записан следующим образом:

оо

5 ( 0 = 2

«**(*-'*)•

(401)

k = — оо

 

 

Для последующих выкладок

удобно воспользоваться рассмот­

ренным в §2 гл. I вероятностным процессом {ях (т), —о о < ;т <

оо},

предполагая, что точки роста реализаций процесса я х (т) совпадают с моментами начала импульсов tk, и в этих точках реализации ях (т) получают приращения, равные ак > 0. Промежутки времени меж­ ду событиями для процесса я х (т) удовлетворяют предположению а) с соответственно измененным значением X, а сам процесс я х (т) удовлетворяет условию 2. Однако я х (т) отличается от я (т) тем, что приращения процесса я х (т) уже не равны числу появившихся им­ пульсов за время от т до т + Ат, а равны суммарному приращению интенсивности за этот период:

Я! (т -ф Ат) — ях (т) = 2 ak- (402)

Пуассоновский процесс я (т) можно рассматривать как частный случай процесса я х (т), когда ах = а 2 = • • • = а.

Воспользовавшись процессом я х (т), выражение (401) можно за­ писать в виде

t

СО

 

5(0 = 5 ( 0 = J e(f — T) CÜTX(T) =

j e(t — x)dn1(x).

(403)

— оо

— со

 

Для этого процесса в [67] приводится аналог теоремы Кемпбелла, которая ниже сформулирована в принятых здесь обозначениях.

Обобщенная теорема о наложении случайных воз­

мущений. Пусть е (т) £ Ln (— с», оо). Тогда для процесса £ (/) се­ миинварианты п-го порядка определяются согласно выражению

 

оо

оо

 

 

х„ К (0] = я. J xnd& (х) J еп(т) dx,

(404)

 

— оо

— оо

 

 

где К среднее число событий,

а <Т (х) — функция распределения

интенсивностей ak.

Выражение

(404) является

частным

Д о к а з а т е л ь с т в о .

случаем выражения (128),

полученного

на основании теоремы 3.

117

Корреляционная функция процесса £ (t) может быть получена с использованием теоремы 3 согласно (128) в следующем виде:

М {[£(/) -

М£(0] [£(* +

s) - М£(/)]} =

 

= X

j

x2d&~(х) j е (т) е (т + s)dt.

(405)

 

—со

— оо

 

 

Выражение для

конечномерной

характеристической

функции

стационарного процесса £ (/), определенного согласно (403) с учетом (126), можно записать в виде

 

 

М ехр 2

ukl(tk)1=

 

 

6=1

 

оо

оо

п

 

= exp {X f

С [ехр (іх 2

— т))— \]d<F (x)dx], (406)

—оо — оо

k = \

 

а выражение для характеристического функционала, согласно (127), такое:

 

 

 

оо

 

 

 

М ехр [i

J

i(t)dg{t)] =

 

 

 

— оо

оо

оо

оо

 

= ехр (X I

J

[ехр (іх

^

e{t — x)dg(t))— 1] dtf (x)dx\. (407)

— оо — оо

— оо

 

Остановимся на случае, когда среднее число вылетов электронов в единицу времени становится очень большим, т. е. когда X -*■ оо. При этом дисперсия процесса я х (т), т > 0, определяемая согласно (123) формулой

оо

х2і (т)] = (Тя.т = Хт I x2dSr{x),

— оо

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

1 оо

(сія.т) 2 [я,(т) — Хт j xd&~(x)]

— оо

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

рирования и нормирования функции £ (/), при таком предельном переходе стремится к гауссовому. Это можно показать, пользуясь методом характеристических функций, следующим образом. С уче­ том (406) характеристическая функция процесса

оо оо

сгё-1 [1(0 — А. j

e(x)dx j xdS" (*)]

—-оо

— се

118

определяется выражением

/ (и) = exp X

іх

е (х)j — 11dS"(х) dx

сю

оо

 

 

■iX j

j" х -^ —е(т)dcT (x)dx\,

 

—оо — оо

 

 

где of — дисперсия процесса \ (t).

 

 

Раскладывая показательную функцию ехр іх

е (т)| в степен-

ной ряд, приводим последнее выражение к виду

 

 

 

 

 

(

Хи2

оо

 

оо

 

 

 

/ (и) =

 

 

 

jj

x2dS'(х) -f-

 

 

exp I-----^ е2(т) dx

 

 

 

 

 

 

\

ч

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

+

Я-

j

j

2 kl

іх - ^ е(хЦк d^(x)dxy

 

 

I

J 'I T

 

 

 

 

 

 

, * =

з

 

 

 

 

 

Учитывая,

что of =

Я, j

x2dSr (x), получаем

 

 

 

 

 

 

 

— 00

 

 

 

 

 

 

 

 

 

 

 

 

/(«) =

 

 

 

 

 

 

©O

 

 

 

OO

OO

QQ

 

 

 

= exp

~

j e2(x)dx;exp{ j

j

[ £

ek (

t )

( x )

 

 

-OO

 

 

—OO OO fc=3

£

 

Так как по условию о\ и X при Ä ,-> оо являются величинами

одного порядка, то of при k > 3 имеет более высокий

порядок ро-

ста, чем X. Следовательно, lim —g- =

0 при всех

k >

3, поэтому в

 

 

 

 

 

 

Я-*оо

CF]

 

 

 

 

пределе при А,-»- оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ди) =

ехр

4 “

j

e2(T)dx\,

 

 

т. е. при

X ->

оо

 

 

 

 

 

— ОО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

оо

 

 

оо

 

 

o f'fë fO — Я,

j

e(x)dx [

xd<F(x)]-*- j e(t — x)dw (x). (408)

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

В заключение отметим, что дробовой эффект изучал и впервые

119

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