Файл: Уриг, Р. Статистические методы в физике ядерных реакторов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 115

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

оо

Следовательно, J | х (t) | dt не сходится и фурье-преобразование

—ОО

неприменимо.

Эту трудность можно преодолеть, введя новую функцию, назы­ ваемую спектральной плотностью мощности, обозначаемую симво­ лом Ф (со), как фурье-преобразование автокорреляционной функ­

ции, т. е.

 

ОО

(4.15)

Ф*ЯИ = F [Ф** М]= 5 Фл* (т) &-'iaxdT.

— ОО

 

оо

| фая (т) | dx

Поскольку фа;* (я) СТреМИТСЯ К нулю при Т ± оо, ТО j

— оо

 

существует и определение (4.15) справедливо.

 

Сущность спектральной плотности мощности может быть пока­

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

в

уравнение (4.15)

определения

ф ^ (т):

 

 

 

оо

 

 

Ф * д » =

5 9**(T)e-Jet =

 

 

 

 

— оо

 

 

 

 

 

т

 

 

Т оо

Л

J

f f{t)f(t + x)dt

e~ i WT dx,

(4.16)

Г

 

 

 

 

 

 

где / (/) определяется как случайная функция, стационарная во времени, равная x(i) в интервале — Т ^ t ^ Т и нулю для 11\ > > Т, как показано на рис. 4.1 и б), т. е.

/(/) =

*(/),

( _ Г < г < Г);

 

/00 =

0,

(\t \> T ) .

(4-17)

Поскольку j |/ (t) | dt существует, то существует его

фурье-

— оо

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

 

 

оо

т

 

Фже((0) = Н

г п

Г

\

f(t)f(t+x)& -i™dtdx.

(4.18)

г - > о о

2.1

J

J

*

 

 

 

— оо

Т

 

Подставив

 

т' =

 

 

 

 

 

t

+ х

(4.19)

66


в уравнение (4.18) и разбив его на два интеграла, получим

ф^ И =

Г

оо

J /(T')e-J'^'dT' =

7!‘|^-^г

 

 

— Г

— оо

= lim4 r [F^ * F ^

]= lim -5F

1 ^ И 1 2.(4-

оо

2 7

Г ^ о о 2

/

где звездочкой помечено сопряженное значение. Равенство (4.20) дает другое определение функции спектральной плотности мощности.

Рис. 4.1. Случайная переменная x(t) (а) и усеченная слу­ чайная переменная f(t) (б).

Величины Фжа; (со) и cpxx (т) фактически содержат одну и ту же информацию относительно статистических свойств х (t). Кроме того, преобразование уравнения (4.15) является обратимым*, так как Фхх (со) -> 0 при со ± 00 , и автокорреляционная функция и спектральная плотность мощности представляют парное фурьепреобразование, т. е.

 

со

 

 

ф.та(®)=

$

q*kc('c) e_i<BTdT,

(4.21)

Ф** (т) = - ^ -

I

фхх (®) eJtDT dco.

(4.22)

* Как будет рассмотрено ниже, исключение составляет физически не­ реализуемый случай белого шума.

3*

67

Эти уравнения известны как соотношения Винера — Хинчина, названные так в честь двух математиков, которые первыми указали на эквивалентность автокорреляционной функции и спектральной плотности. Поскольку ср.™ (т) и Фхх (со) симметричны, отношения Винера — Хинчина можно выразить, используя косинусоидальное преобразование. В этом случае уравнения (4.21) и (4.22) примут вид:

Фхх (со) =

jj срхх (т) cos сотdx,

(4.23)

qi (т) =

1

С Ф . (co)coscoTrfco.

(4.24)

----

 

2п

J

 

 

со

 

Удачной аналогией, которая действительно приводит к понятию спектральной плотности мощности применительно к Фхх (со), яв­ ляется представление х (t) в качестве напряжения, прилагаемого к сопротивлению в 1 ом. Тогда ф$. является средней рассеиваемой мощностью, а Фм (со) — средней мощностью на единичную часто­ ту, вычисленную для данной частоты. Такое представление может

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

(4.22) для

частного случая, когда т = 0, т. е.

 

 

 

1

 

00

(4.25)

Ф.г-д,- (0) = Ф* = —

 

jj Фх-х И сЬ.

 

00

 

Для линейной частоты / уравнение (4.25) переходит в

 

Ф**(0) = ,Ф5=

$

фхх (/) df-

(4.26)

оо

 

 

Это значит, что среднеквадратическое значение х((), которое, как было показано равенством (4.8), эквивалентно cpxx (0), определяется площадью под кривой спектральной плотности мощности Фхх (/).

§4.5. Специальные автокорреляционные функции

испектральные пло'рности мощности

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

68


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

Ф** (т) = Фи (т) + Фа» (т) +

...+ срй (т) +

 

+ • • • + <PNN(т) =

S

Фг г (т),

(4.27)

 

i =

1

 

где x { t ) = x 1 (t) + X'i (t) + . . . + x i (t) +

. . . + x N (t) = ^ x t {t)

(4.28)

 

 

i= i

 

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

ф** И = Фц И

+ ф .а N - I - ... + Фш (со) = f

ф „ (со). (4.29)

 

 

г= 1

 

 

,

_— 2Жаг

 

С

 

О

-<J

О

+Ь)

Рис. 4.2. Графики функций ср1Х (т) и Ф** (со) для постоянной

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

Постоянная. Постоянная может рассматриваться как переменная, которая не имеет частотной составляющей кроме со = 0 и всюду кор­ релирует сама с собой. Если константа имеет значение (амплитуду) а, то

ф** (т) = Е [а• а] = а2,

(4.30)

оо

 

Фа;а.(со)= ^ а2е— dx 2яа2б(со).

(4.31)

—ОО

 

На рис. 4.2. изображены функции ср** (т) и Фжж (со) для постоян­ ной.

Синусоида. Для синусоиды

х (t) = A sin со„ t

(4.32)

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

 

Л2

(4.33)

Фх.г (Т) = — cos СО0 Т.

69



Поскольку рассматривается процесс стационарный во времени, положение начала произвольно. Следовательно, наличие фазового угла в уравнении (4.32), даже равного 90° (что может изменить си­ нус на косинус), не повлияет на результаты, представляемые урав­ нением (4.33). Важность этого процесса состоит в том, что амплитуда автокорреляции не стремится к нулю с возрастанием т, как в случае непериодических изменений. Следовательно, используя большие значения т, можно установить наличие периодических составляю­ щих в шуме.

Спектральная плотность мощности для синусоиды должна вы­ ражаться, как это интуитивно чувствуется, дельта-функциями Ди­ рака, существующими при ± ® 0. Это легко показать:

ф = «

= т I cos ш0 т cos сотс/т =

 

= —

[б(со-(-со0) + 6(со—со0)].

(4.34)

Автокорреляционная функция и спектральная плотность мощности синусоиды приведены на рис. 4.3.

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

Ф-г-Дю) — а.

(4.35)

70

Автокорреляционная функция при этом имеет вид

оо

ср. (т) = С cos on cl(£>= ад (т), (4.36) 2я J

— оо

что означает, что функция х (t) полностью некоррелирована, кроме точки т = 0, т. е. функция полностью случайна. Графики функций Фд-а; (т) и Фж* (ю) приведены на рис. 4.4.

 

<Р,М

а

ФиМ

 

^

а

 

zj

 

 

О

+Т -и

0

Рис. 4.4. Графики функций фХ1 (т) и Ф** (со) для бе­ лого шума.

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

Ф « И = а. (М <Ю о);

Ф . г - * (со) =

о, ( I со I > со0),

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

 

оо

 

<й0

Фсс.-c (т) = — | Фп И

eiB td(0 = - ^ - j cos coxdco :

 

 

—©о

a

flCOp

SinШо Т

ят

Я

(ОрТ

(4.37)

(4.38)

Графики фжа; (т) и Фхл; (со) для низкочастотного белого шума пред­ ставлены на рис. 4.5.

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

71