Файл: Уриг, Р. Статистические методы в физике ядерных реакторов.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