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

книги / Твёрдотельная фотоэлектроника. Физические основы-1

.pdf
Скачиваний:
15
Добавлен:
12.11.2023
Размер:
15.42 Mб
Скачать

При £2 = h двумерная случайная величина х\Х2 вырождается в одномерную величину х\ = х 2 — представляет собой просто средний квадрат:

5 * (ii,ti)= jОО x\p{xi,t{]dx = {x\) = {xl) = a2x {ti) + { x i f

(5.1.10)

— ОО

 

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

Исследование случайного процесса, а также реакции на случайный процесс различных устройств существенно упрощается при его стационарности. Слу­ чайный процесс называется строго стационарным (или «стационарным в узком смысле»), если его многомерная плотность вероятности p (x i,ti\хъМ', ■■■xn,tn)

зависит только от интервалов <2 - t i ; *з - h; ... tn - t\ и не зависит от поло­ жения этих интервалов в области изменения аргумента t.

Во многих приложениях теории случайных процессов условие стационарно­ сти ограничивается требованием независимости от времени только одномерной и двумерной плотностей вероятности (такой случайный процесс называется «стационарным в широком смысле»). В этом случае автокорреляционная функ­ ция зависит опять не от самых моментов времени t\ и t2, а только от интерва­ ла между ними T = t2 ~ t\. Для стационарных процессов (в узком и широком смысле) выражения (5.1.4)—(5.1.6) и (5.1.9) можно записать без обозначения фиксированных моментов времени:

ОО

(5.1.11)

ОО

ОО

(5.1.12)

— ОО

(5.1.13)

ОООО

(5.1.14)

где х\ = x(t) и Х2 = x(t + т). При т = t^ — t\ = 0 автокорреляционная функция равна среднему квадрату случайного процесса в любом сечении t — ti.

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

Практически при усреднении по ансамблю оперируют с реализациями N (N » 1) идентичных систем на некотором интервале времени Та - Д ля эргодического процесса достаточно взять одну систему и получить ее представитель­ ную реализацию (несущую в себе все основные черты исследуемого процесса) за гораздо большее время Т ~ N T A . За «экономию» в числе систем приходится расплачиваться временем.

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

Для эргодического процесса соотношения (5.1.11)—(5.1.14) эквивалентны следующим выражениям, в которых операция усреднения по интервалу времени обозначена прямой чертой сверху:

Т

2

Х(4) = гЙЗо? /

X^

dt'

(5.1.15)

 

 

 

 

т

 

 

 

х2{t)=

lim

^

f

x2(t)dt,

(5.1.16)

V'

Т —>ооТ J

W ’

 

 

 

 

_ т

 

 

 

<rx = x 2(t) -

(s(i))

 

(5.1.17)

 

 

 

 

Т

 

 

 

 

 

 

2

 

 

Вх (т) = x(t)x(t + т) =

lim

— /

x(t)x (t + T)dt.

(5.1.18)

 

 

Т ►ooT J

 

 

T

2

Из определения (5.1.18) вытекает четность функции Вх (т) = Вх (—т) для эргодических стационарных процессов.

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

При анализе стационарных шумовых процессов применяют автокорреляци­ онную функцию (иногда ее называют ковариационной функцией) в несколько ином виде, пригодном для описания шумов, то есть центрированных величин — x (t) — x(t):

Определенная таким образом функция автокорреляции (ковариационная

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

вательно, Вх (т = 0) = сг%. Если т достаточно велико, то значения |x(t) —:r(f)J

и |ж(t + т) —x(t)] делаются независимыми и Вх (т) при т —>оо стремится к

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

Вводится также нормированная автокорреляционная функция

Ях (т) = ^ Р - .

(5.1.20)

°т

 

Очевидно, что Дх(0) = 1.

В качестве характеристики связи между значениями двух стационарных и стационарно-связанных эргодических случайных процессов X (t) и Y (t) вводят взаимно-корреляционные функции, определяемые выражениями

Т

2

Вху(т) = (x(t)y(t + т)) = x {t)y (t + т) =

lim ^

[

x(t)y(t + r)dt,

 

Т —>оо 1

J

 

 

 

~1

(5.1.21)

 

 

2

 

Вух(т) = (y(t)x(t + т)) = y (t)x (t + т) =

lim i

[

y(t)x(t + r)dt.

 

T —>ocl

J

 

 

 

_ T

 

 

 

2

 

Как и для детерминированных колебаний, взаимно-корреляционные функ­ ции стационарных случайных процессов не изменяются, если сдвиг на т одной из функций x(t) или y(t) заменить на сдвиг в обратном направлении другой функции. Поэтому

Вху (т) = x(t)y{t + т) = x(t - r)y(t),

Вух (т) = y(t)x(t + r) = y{t - r)x(f).

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

Вху (т) —Вух ( т ) ,

ВуХ (т) = Вху (—т ) .

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

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

Определим корреляционную функцию алгебраической суммы двух случай­ ных процессов X (t) и F (i) с нулевыми средними (ж = у = 0): s(t) = x ( t ) ±

=*=2/(*)• Из соотношения (5.1.18)

В3(т) = s(t)s(t + т) = [x(t) ±y{t)\ [x(t + т) ± y(t + т)] =

= x(t)x (t + т) ± x ( t ) y ( t + r) ± y ( t ) x ( t + т) + y(t)y(t + T) =

= Вх (т) ± Вху (т) ±

Вух (г) + By (т).

При т = 0 имеем Вх (0) = <7%, Ву (0) = crj), Вху (0) = Вух (0)

и дисперсия сум­

марного и разностного процессов составляет а3 = В3(0) =

+ (Ту ± 2Вху (0).

Если процессы X (t) и Y (t) независимы, то дисперсия (средняя мощность)

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

 

<?l = al + crl-

(5.1.22)

В противном случае мощность процесса s(t) в зависимости

от знака Вху(0)

может быть больше или меньше суммы дисперсий. Так для двух совпадающих (полностью коррелированных) процессов X (t) = Y (t) флуктуации в каждый момент складываются и а3 = 4ах = 4сгу.

5.1.2. Спектральная плотность мощности (энергетический спектр) случайного процесса. Обозначим через Gxw(u>) спектральную функцию и- реализации x„(t) с продолжительностью Тр стационарного центрированного случайного процесса X(t).

Для нахождения спектральных характеристик случайного процесса X (t) осуществим операцию, обратную переходу от спектра периодического колеба­ ния к спектру импульсного сигнала. Нарастим ^-реализацию большим числом N —►оо других реализаций с такой же длительностью. При этом спектральная функция Gx (ш) случайного процесса X (t) усредняется по реализациям:

^_____

1

^

х N

Tf

Gx (u>) = GXU(OJ) =

lim — 'S^Gxviui) =

lim — V ]

/ x u(t)exp{-jut)dt.

 

N-+oo N

*—'

N —>oo N

I

 

 

U=1

u=1- i / 2

Достаточно поменять порядок суммирования и интегрирования в этом со­ отношении, чтобы убедиться, что среднее значение или спектральная функция Gx (ш) такого процесса равна неинформативному нулю:

 

Л

N

\

°?____

Gx (u)=

\ u ^ i — '^ 2 x l/( t ) \ e x p ( - j u t ) d t =

/ x u(t)exp{-jujt)dt = 0.

- о о

'

"=1

'

- о о

Пределы интегрирования можно сделать бесконечными, так как вне интервала [-Гр/2, Гр/2] х „(f) = 0; кроме того, для центрированного процесса xu (t) = 0.

Если бы ^-реализация наращивалась во времени ее периодическим повто­ рением, то в результате получился бы линейчатый спектр с амплитудами

Gx (ч>) = Gx (ч^п) = ЧтгСр­ одна ко в случае усреднения разных реализаций

Gx (ч>) Gxu(j-^) dnTp —I^n.|бхр (jipnjTp —|an| exp (j(pn)Tp —0.

Так как амплитуды и фазы в разных реализациях не коррелированы, то |dn|exp(j</?) = |an| - exp(j<pn) — среднее от произведения равно произведению средних.

Гармоники с одинаковой частотой в спектрах разных реализаций стационар­ ного случайного процесса хаотически сдвинуты во времени друг относительно друга, поэтому их фазы равномерно распределены в интервале [0,27г]. Это об­ стоятельство и обуславливает нулевой результат усреднения: exp[jV„] = 0.

Тем не менее, спектральный анализ широко используется для исследования флуктуирующих величин. Шумы можно охарактеризовать спектральной плот­ ностью среднего квадрата их величины (средней мощностью, выделяемой на нагрузке 1 Ом), поскольку величина среднего квадрата не зависит от соотноше­ ния фаз суммируемых гармоник. Спектральную плотность средней мощности (средняя мощность на заданной частоте, приходящаяся на полосу частот 1 Гц), имеющую размерность энергии [Вт/Гц = Вт • с = Дж], называют также энер­ гетическим спектром случайного процесса.

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

Энергия любой из реализаций x v {t) с конечной длительностью Тр стацио­ нарного процесса X (t) выражается через модуль спектральной функции этой реализации Gxv{u) — уравнение (4.1.33):

 

Тр/2

оо

 

 

 

& =

J x v2 ( t ) d t = ^

J

GXv М

dui.

(5.1.23)

-Г р/2

 

 

 

 

Тогда средняя мощность реализации xu(t) на интервале Тр

 

 

7р/2

 

о° . .

 

 

 

1

f |6*„М | du.

(5.1.24)

% = % ] x* ® d t = h l

 

 

-Tp/2

Здесь GXu (ч>) /Г р — спектральная плотность средней мощности ^-реализации с продолжительностью Тр.

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

Усреднив (5.1.24) по ансамблю реализаций, получают спектральную плот­ ность мощности [Вт/Гц] или энергетический спектр случайного процесса

W{w)

&xv (w) I

(5.1.25)

 

При этом средняя мощность случайного процесса X (t)

 

Оо

 

 

 

 

 

 

JL

- fJ

GXV (&)

 

 

2.!—

|СЖ1Ди;п)|

^ | а п|2 > 0.

ТР

Т7.

Трг

"27

7^2

 

п

Р

 

 

2тг

 

 

 

 

 

 

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

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

% о И = ( ^ ) ) 227гй(а;) + W{u>).

(5.1.26)

При интегрировании по f=u>/2n первое слагаемое дает (а;(£)^ — мощ­

ность постоянной составляющей, а второе — мощность флуктуационной со­ ставляющей, то есть дисперсию а%:

ОО

<?l = ^ J w M dw .

(5.1.27)

—ОО

 

5.1.3. Соотношение между энергетическим спектром и автокорреляци­ онной функцией случайного процесса.Установим связь между энергетиче­ ским спектром случайного процесса и его автокорреляционной функцией. В соответствии с (4.1.16) спектральная функция его ^-реализации xu(t) с про­ должительностью Тр

 

Гр/2

GXu (ч^)

/ х и (t) exp (—jwt) dt.

 

-ГР/2

При этом из (4.1.29) следует, что

гР/ 2

/x u(t - r)exp(-ju>t)dt = Gxu(u>)exp(—j u r ) .

-r P/2

Подставим эти соотношения в (4.1.29):

Г р /2

оо

J х и (t) x„(t - т) dt =

J GxvG*xuexp ( J U T ) <ко.

—Т р /2

- о о

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

Гр/ 2 Т р /2

П т —

/

x v(t)xv( t - T ) d t = —

/

xu(t)x„(t

- r)dt = x v{t)xu{t -

т) =

Гр—*oo J p

J

 

 

 

Ip

J

 

 

 

 

 

 

 

- Г р / 2

 

 

 

—Г р /2

 

 

 

 

 

 

 

=

lim

7 l - a^ h,)L-

eXp(jfa;t)dh;= —

/ 1 ^ 1 exp(jut)du

 

 

Tp—>oo 2?r J

T p

 

 

'

2ir

J

T

KJ

 

 

 

—oo

 

 

 

 

 

—oo

 

 

 

ИЛИ

 

 

 

 

 

oo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5 (r )

= —

 

[

W(u)exp(jwT) dw.

 

 

(5.1.28)

 

 

 

 

27Г

у

 

 

 

 

 

 

 

—OO

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

Прямое преобразование Фурье автокорреляционной функции приводит оче­ видно, к спектральной плотности мощности

ОО

 

W (u )= J В(т)ехр(—jwT)dT.

(5.1.29)

— ОО

Уравнения (5.1.28) и (5.1.29) представляют собой математическую запись теоремы Винера-Хинчина.

Так как функция автокорреляции стационарного процесса четная, то

 

ОО

ОО

 

W (ш) =

J В (г) exp (—jwr) dr = 2 J В(т) costотdr,

 

—оо

О

(5.1.30)

 

ОО

 

В(т) = -

I W (u>)cosu>Tduj.

 

 

тг J

 

 

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

п

т

0

деленных интервалом т, и спек-

 

Wi (to)

h*-—|

трально. Здесь имеется аналогия с

 

У ////У ////Л

S888Sa

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

 

 

 

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

-(DJ-WQ

 

 

процессов. Однако при детермини­

 

 

 

рованных процессах связь между

 

 

 

мгновенными значениями функци­

Р и с . 5.1.2.

Широкополосный и

узкополосные

ональная («жесткая»), в то время

энергетические спектры

 

как при случайных процессах эта

 

 

 

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

Таким образом, автокорреляционная функция является достаточно полной характеристикой стационарного случайного процесса, она позволяет анализи­ ровать поведение во времени этого процесса. При т = 0 она равна среднему квадрату случайной функции (или дисперсии — для центрированной случай­ ной величины), при т —►оо — квадрату ее средней величины. Преобразование Фурье, примененное к автокорреляционной функции, представляет исчерпыва­ ющие сведения о спектральной плотности мощности случайного процесса.

Определим автокорреляционную функцию для нормального стационарного процесса с нулевым средним и энергетическим спектром, равномерным на всех

частотах -оо < и>< оо (белый

шум). Если в (5.1.28) подставить W (и) = Wo,

то получим

ОО

 

 

 

£ ( г ) = ^

[ exp(jur)du> = WO5(T ).

(5.1.31)

2тг

J

 

Корреляционная функция равна нулю при всех значениях т, кроме т = 0, при котором В(т) обращается в бесконечность — дисперсия (мощность) бе­ лого шума бесконечно велика. Поэтому, когда говорят, что корреляционная функция реального процесса есть дельта-функция, то это, конечно, является приближением.

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

лосные фильтры (рис. 5.1.2).

 

 

 

После фильтра с равномерным

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

h = wi/27r корреляционная функция становится

 

 

U>1

 

LJ1

 

В\ (т) = —

[ Woexp (JUJT) du = —-

[ cosurdu —

27г

J

2тг

J

 

-UJl

 

-UJi

 

 

= Wo2f1S^ ^ = a

f -S ^ ^ = a ^ R 1(T). (5.1.32)

 

 

 

UJ\T

0J\T

P и c. 5.1.3. Нормированные корреляционные функции случайных процессов с энергети­

ческим

спектром равномерным

в полосе: |u>|<u>i и \и\ <

= ш\/5

(а); UQ -

£1/2

^

^ CJQ +

£ 1 / 2

( б )

 

 

 

 

 

Здесь

а\ — дисперсия

или средний квадрат

шума на

выходе

фильтра,

R \{T ) — нормированная корреляционная функция.

 

 

 

 

На рис. 5.1.3с приведены нормированные корреляционные функции i?i(T), а также R 2 {T ), полученная для узкополосного фильтра с равномерным пропус­ канием в полосе частот от нуля до Fi = fii/27r — в несколько раз меньшей f\ (эта полоса обозначена на рис. 5.1.2 штриховкой). Сужение спектра приводит к соответствующему растягиванию графика R 2 (T ) по оси т .

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

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

U J Q -f~^i /2

и>о Н-Г21 /2

Wo

J

Вз(т) = 2тт

 

—(«-«-’о -h^i/2)

exp(jujr)duj+

J

exp(jcor)d(

(CJQ $"2I/2)

LIL

sin-niT

2

= (W0 2Fi) I —g^-cosw oT I =<7зД3(т). (5.1.33)

Огибающая функции R3(т) — штриховая линия на рис. 5.1.36 имеет вдвое большую протяженность, чем /?2 (т), а частота высокочастотного заполнения R3(т) равняется центральной частоте спектра шума и>о.

Примерный вид реализации узкополосного шума приведен на рис. 5.1.4. Шумовые колебания происходят в среднем с частотой а>о> ПРИ этом огибающая

Р и с . 5.1.4. Реализация случайного процесса, корреляционная функция которого показана на рис. 5.1.36 (масштабы по осям t и т разные) [15]

шума изменяется со временем сравнительно медленно — подобно шумовому процессу с наивысшей частотой Пь

5.2.Примеры случайных процессов

5.2.1.Биномиальный процесс. К биномиальному распределению приводит решение задачи Бернулли. Постановка этой задачи излагается следующим образом. Проводится m независимых испытаний, при каждом из которых происходит одно из противоположных событий: А (вероятность его появления