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

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

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

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

Добавлен: 16.10.2024

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

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

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

Поскольку т — переменная интегрирования, то, введя замену в пер­ вом интеграле

т = — т'

 

(4.113)

и опустив штрих, получим

 

 

О

со

 

аху (со) = 2 ^ — cpx;/ (—Т) е ^ г dx +

2 $ срэд (т) е~ iит с/т.

(4.114)

о

 

Используя равенство, даваемое уравнением (4.56)

 

Ф.л:У(—т) = (рух(т),

(4.115)

и меняя пределы интегрирования в первом интеграле, можно по­ лучить, что

• °о

оо

 

 

 

Gxy М = 2 5 фух (т) ejQTdx +

2 5

ФЛ.„ (т)

с/т.

(4.116)

о

о

 

 

 

Подстановка соотношения Эйлера для комплексной экспоненты

(4.103) в уравнение (4.116)

приводит к уравнению:

 

 

 

ОО

 

°ху Н = Сху И — jQrcy И

= 2 5 ф;у.-с (т) (cos (ОТ -f j si п сот) dx +

 

оо

b

 

 

 

 

+

2 ^ cpxy (т) (cos сот—j sin сот) dx =

 

 

о

 

 

 

оо

 

 

=

2 5 [фук (т) + фху (т)] cos сотdx—

 

 

о

 

 

 

ОО

 

 

— j2 5 [ф.г-у (т)— фу^^л Sin сотс/т,

(4.117)

 

о

 

 

из которого могут быть выделены действительная и мнимая части:

 

 

оо

 

 

СхУИ

= 2 I

[фад (т) + фух (т)] cos сотdx

(4.118)

Qxy И

= 2

$

[фзд (т)— Фуя. (т)] sin сотdx.

(4.119)

 

 

о

 

 

Из этих выражений видно, что Сху (со) — действительная четная функция частоты, в то время как Qxy (со) — действительная нечет­ ная функция частоты. Следовательно,

Сху (“ ) =

Сху (-со) =

Сух (со) = Сух (-со)

(4.120)

и

 

 

 

Qxy (®) --

Qxy ( М) --

Qyx (®)--Qyx( ®)-

(4.121)

88


Поскольку функция GX1J (со) комплексна, она может быть записана в полярных координатах:

Gxy (®) = I @ху (m) I eJ0*y (м) ( 0 < с о < «> ),

(4.122)

где | Gxy (со) | — амплитуда и 0ЗД (со) — фазовый угол

могут быть

выражены через мнимую и вещественную составляющие, а именно:

| Gxy (со) | = [С* (со) + Qly (со)] 1/2

(4.123)

и

 

Qxy (со)

(4.124)

0ад (“ )= arctg

CXIJ (со)

 

В практике используются два приема для измерения взаимной спектральной плотности. Первый заключается в измерении взаим­ ных корреляционных функций срЛ.у (т) и (рух (т), расчете действитель­ ной и мнимой составляющих по уравнениям (4.118), (4.119) и опре­ делении взаимного спектра и его фазового угла по уравнениям (4.123) и (4.124). Другой метод заключается в непосредственном из­ мерении вещественной и мнимой составляющих с помощью частотно­ полосного оборудования и определении амплитуды и фазы взаим­ ного спектра по уравнениям (4.123) и (4.124).

Взаимные корреляционные функции могут быть рассчитаны по

С (со)-

и Q (со)-спектрам с

помощью следующих

соотношений:

 

 

_ 1_

ej“* da = —

I

ФХу(со) ejt0T da +

 

Фзд(Т) = 2я” j Ф адИ

 

 

 

оо

 

оо

 

 

СО

_ /

-1—- Г Ф (со)Ыитс/со=— Гих 1X1 e - j “Tdco+ — Г

ху Ш) ei™ da=

2л J

ху '

 

2я J

2

 

2я J

2

0

 

 

0

 

 

о

 

 

 

оо

 

 

 

 

 

 

= ~Ь. I Т

[С“-Г ^ —

И ] tcos ®т—j sin ИТ] da +

 

 

о

 

 

 

 

 

 

1

оо

 

 

 

 

 

 

Г 1

 

 

 

 

 

+ —

[Сжг, (со)—jQxy (со)] [cos сот+ j sin сот] da =

 

 

О

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

="У И-[С*»^ +Сух ^

 

coswt+

 

 

 

 

о

 

 

 

 

 

 

+ - у №хи (“ )—Qyx (®)] S iп сотJdco+

 

 

 

 

oo

 

 

 

 

 

 

+] ~L I {т ^

^ +

 

cos“T+

 

 

 

 

0

 

 

 

 

 

 

+

у [Cyx (со) —Cxg (со)] sin c o t | da.

(4.125)

89



Используя равенства (4.120) и (4.121), соотношение (4.125) можно упростить и привести к виду:

оо

Фч, (т) = ^ J fc -^/ И cos MT+Q*» И sin ИТ1da- (4.126)

о

Связь между входным и выходным сигналами. Связь между вход­ ным и выходным сигналами, описанная в § 4.9, может быть выражена через односторонние спектральные плотности. Уравнение (4.93), включающее спектральную плотность мощности, записывается при этом следующим образом:

Gyy Н = Я* (со)Я (со) Gxx (со) = | Я (со) |» Gxx (со),

(4.127)

а уравнение (4.100), включающее взаимную спектральную плотность, как

Gxy (<*>) = Я (со) Gxx (со).

(4.128)

§4.12. Зависимость корреляционных функций

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

от среднего значения

Автокорреляционная функция. Рассмотрим случайную перемен­ ную, имеющую среднее значение р*. Автокорреляция определяется в соответствии с выражением

Фх* (т) = Е [х (t) х (t + т)]

(4.129)

и показана на рис. 4.12. График автокорреляционной функции сим­ метричен относительно вертикальной оси, а максимальное значение, равное среднеквадратическому значению находится при т = = О..При стремлении т к ± оо срг;с (т) стремится к р;, т. е. автокорре­ ляционная функция смещается по вертикали на величину р | из-за того, что х (t) имеет среднее значение рж, что и показано на рис. 4.12. Поскольку случайная переменная x(t), имеющая среднее значение р, может рассматриваться как сумма

х (t) = х' (t) + рж,

(4.130)

где х' (t) — та же переменная с нулевым средним значением, авто­ корреляционная функция х (t) может быть выражена как сумма ав­

токорреляционных функций х' (() и рд., т.

е.

 

Ф**(т) = Ф*'л;'(тО +

!1х-

(4.131)

Легко показать, что дисперсия а® равна разности между значением автокорреляционной функции при т = 0 и смещением горизон­ тальной оси.

90


Спектральная плотность мощности. Двусторонняя спектраль­ ная плотность мощности может быть получена с помощью преобра­ зования Фурье автокорреляционной функции. Применив преобра­ зование Фурье к уравнению (4.131), получим:

ФКЛ. (со) = Фх-х>(со) +2щ% б (со),

(4.132)

т. е. наличие среднего значения, отличного от нуля, добавляет дель­ та-функцию Дирака в начале координат, как показано на рис. 4.13.

Рис. 4.12. Автокорреляционная функция случайной переменной, имеющей среднее значение

Iх* * 0.

Рассмотрим связь между односторонней и двусторонней спек­ тральной плотностью мощности, определяемую уравнением (4.101):

ОО00

Фл- = ~ J Ф** N dc° =

j G*xИ dw’

(4-133)

-ОО

0 ”

 

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

Gxx (<й) = 2Фхх(а>),

(со >

0);

 

Gxx (со) = Ф ,я (0) =

2яр.5,

(со =

0);

 

<5**(«») = 0,

 

(со <

0)

(4.134)

и показано на рис. 4.13.

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

91

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

со

 

 

^ = 1 л j ° хх^

= + ^ •

(4 .1 3 5 )

о~

 

 

Если в это интегрирование не включать начало координат или если среднее значение равно нулю, площадь под кривой равна диспер­ сии Ох.

(и>0)

(и=01 1й]<01

Рис. 4.13. Спектральная плотность мощности пере­ менной, имеющей среднее значение

Щ2 Ф0.

Рис. 4.14. Взаимная корреляционная функция двух переменных, каждая из которыхимеет свое среднее значение.

Взаимная корреляционная функция. В случае взаимной корре­ ляции наличие постоянных составляющих рж и р ц в переменных х (t) и у (t) приводит к кривой взаимной корреляции, показанной на рис. 4.14, т. е. взаимная корреляционная функция смещается по вертикали на величину р^ру. Математически это выражается как

Ф*и (т) = Фх'в- W + р* Ру.

(4.136)

где х' и у' —■переменные со средними значениями, равными нулю. Взаимная спектральная плотность. Двусторонняя спектраль­

ная плотность определяется с помощью преобразования Фурье урав­ нения (4.136):

ф ху М = ф *'у М + 2яр с р„ 6 (со),

(4.137)

92