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

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

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

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

Добавлен: 16.10.2024

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

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

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

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

(4.103):

0*в И = 2 Ф зд((о),

(со > 0 );

 

О ,> ) =

Фэд(0),

(© =

0);

 

@ху (®) =

0>

( и <

0).

(4.138)

Графики этих функций представлены на рис. 4.15.

Рис. 4.15. Взаимная спектральная плотность двух переменных для случая, когда обе име­ ют средние значения.

§ 4.13. Когерентные функции

Наличие постороннего шума или других входных переменных может отрицательно влиять на качество измерений. Другими факто­ рами, ухудшающими измерения, являются отклонения от предпола­ гаемой линейности системы и конечное время измерения. Обычным показателем качества измерения для отдельной системы вход — выход является когерентная функция, которая может рассматри­ ваться как показатель внутренней согласованности измерений. Рассмотрим динамическую систему без учета каких-либо практи­ ческих ограничений в частотной области, в которой могут проводить­ ся измерения. В § 4.11 были получены соотношения, связывающие сигналы вход — выход:

G»y = I

И I 2 °хХИ

(4.139)

и

 

 

Gxy (со) =

# е(<в) Gxx (со).

(4.140)

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

93

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

I Нс(со) |«

I Gxy (со) PlGxc (to)2

1Gxy (со) |2

(4.141)

у2(со)

Gуи(CO)lGxx (СО)

Gxx (СО) Gyy (со)

\НР(“ ) I2

 

Когерентная функция может быть также определена через дву­ сторонние спектральные плотности:

у2(со) = IФху (со) |»/ФетИ ФУ И-

(4-142)

Если все предположения, использованные при выводе уравнений (4.139) и (4.140), справедливы, то когерентная функция равна еди­ нице. С помощью неравенства (4.78) можно показать, что когерент­ ная функция удовлетворяет неравенству

у2(со) = I Ф х у ( « ) I*

(4.143)

Ф.Г.Х- (СО) Фуц (СО)

 

т. е. когерентная функция всегда меньше или равна единице. Когерентная функция может иногда рассматриваться как мера

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

§4.14. Корреляция при использовании двух детекторов

иизмерение спектральной плотности

Взаимная корреляция выходных сигналов двух детекторов. Рас­ смотрим систему, показанную на рис. 4.16, где Н (со) символизирует ядерный реактор, в который поступает возмущение i (/). Выход­ ной сигнал—нейтронная плотность—измеряется двумя регистриру­ ющими системами, расположенными рядом, с передаточными функ­ циями # ! (со) и # 2 (со). Выходные сигналы детекторов ух (t) и уг (t) искажаются детекторными шумами (/) и т 2 (() соответственно. Применив интеграл свертки, получим, что выходной сигнал реак­ тора должен описываться выражением

00

 

n(t)= $ h(k)i(t— k)dl.

(4.144)

94


Входные сигналы детекторных систем имеют вид

 

хг (t)

=

т1 (t)

+

п (t)

(4.145)

и

 

 

 

 

 

x2 (t)

=

т2 (t)

+

п (t).

(4.146)

С помощью интеграла свертки и уравнений (4.145) и (4.146) выход­ ные сигналы ух (t) и у 2 (t) описываются следующими выражениями:

Рис. 4.16. Структурная схема корреляционного экспе­ римента с двумя детекторами, имеющими собственные шумы.

 

 

ОО

 

 

 

00

 

 

0 1 (0 = S

Ai(ll)*i(i! —

 

S hiWlmxit — \i) + n(t — р )]ф

 

ОО

 

 

 

0 0

0 0

 

=

5

hx(ii)nix(t— р.)ф + 5

S h1([i)h(X)i(t— y,— K)d'kd\i,

(4.147)

 

— 00

 

 

 

— 00 — 00

 

 

 

 

 

 

ОО

 

 

 

 

 

 

Уг — S

h2(l)m2{t— l ) d l -{-

 

 

 

 

+ 5

$

h(l)h(r\)i(t— l — 4\)dT\d£.

(4.148)

по

Отметим, что p,, Я, £

и t) являются переменными интегрирования

времени.

Взаимная

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

ух (t) и

у2 (t)

определяется как

 

 

 

 

 

 

фи (Т) =

Ф » , У, ( Т ) = Е [ 0 1 (t) У2 (t + Т)] .

(4.149)

95


Если подставить соотношения (4.147) и (4.148) в уравнение (4.149), используя соответствующие определения корреляционной функции, то получим

 

СО

0 0

 

Ф12 (т)“

S

S ^1 (М1) ^ 2 (£) ф/и, т, (т -f- р 1) d\l dl -{-

 

 

— 00 — 00

 

00

 

 

 

+ 5 5 5 К О) К (s) h ft) ф„г, I (т + p — £—Ti) dp dl dr\ +

 

—00

 

 

 

CO

 

 

 

"b 555

(p ) ^ 2 (£) ^ ft) Ф«, £ft ~Ь li — i ft- X) dp dl dX -ft

+ Ш hi ^)>h (l) h (X) h ft) Фгг ft + p - 1+ Я - T))dp dl dXd11.

(4.150)

— OO

 

 

 

Поскольку шумы тх (() и nil (0 по предположению не коррели­ руют друг с другом и с входным сигналом г (/), взаимные корре­ ляционные функции cp„,tm (т), ср„м- (т), (pnuJ (т) равны нулю и пер­ вые три члена в (4.150) исчезнут. При этом получим:

оо

Фи ft) =

 

= 5555 /11 (И-) h 2 (l) h (X) h ft) ф „ (т + р - | + X—Ti) dp dl dX ф . (4.151)

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

(/) = /?2 (/),

(4.152)

то, введя замену переменных

 

 

 

 

 

=

р,

 

(4.153)

и исключив 1]

и I,

 

 

р = т1—Я

 

(4.154)

можно получить:

 

 

 

 

 

со

оо

оо

 

 

 

Фи м = 5

5

5 M P)M <* + P)dp X

 

 

00

 

— оо — ос

-0 0

 

 

 

 

 

 

 

 

 

 

 

X

5

h(X)h(p + X)dX

Фг£ ft— а —Р) dp da

 

 

-00

 

 

 

 

 

 

 

 

00

00

 

 

 

 

 

 

=

I

5

Ф*1

Ai(o) Фай (р) Ф» (т —

P)dpda.

(4.155)

— 00 — со

96


Для входного сигнала, имеющего характер белого шума,

 

Фгг (т — ар) = К8 (т—а —р),

(4.156)

где К —константа, б—дельта-функция Дирака. Введя

новую за­

мену переменной

 

 

 

 

 

 

0

= р +

сг

 

(4.157)

и исключив р из уравнения .(4.155),

получим

 

 

со

 

 

 

 

Ф12 (т)=

й фя, kt (а) ФЛ,,. (0—а) Д6

(т—0) da dQ=

 

 

— 00

 

 

 

 

оо

 

 

со

 

 

'= К $

фя,я, ( O f ) Фяя(9—or)do 5

8 (т—0)d0 =

 

 

00

 

 

 

 

 

фя,я, (o)yllk(x— o)dc,

(4.158)

оо

 

 

 

 

 

поскольку § б(т— 0) dQ равен

1 и 8

(т— 0) отлично от нуля толь»

— ОО

 

 

 

 

 

ко при т = 0 . Предположим, что импульсные характеристики h(t) и h± (t) имеют вид:

h(t) — A exp (—at),

(t^O);

 

 

h(t) = 0,

 

( f < 0);

(4.159)

hid) = В exp ( —fit),

( t > 0);

 

 

K(t) = 0,

 

(t<0).

(4.160)

Тогда

 

 

 

 

 

Фяя(т) — 4 ~ e_a |T|>

(4-161)

 

 

2a

 

 

 

Фй.я1(т)= -^ е-Р 1'*1

(4.162)

и уравнение (4.158) перейдет в

уравнение

 

 

Ф12

И2

0 2

е —р lx—cr| da =

 

 

— e - a IffI л _

 

 

 

2(5

 

 

 

КА2В2

Г е - а

(orI—Р 1г—а|

(4.163)

4аР

J

 

 

 

 

Если интеграл разбить на три части при а = О и а =

т и

под­

ставить (о — т) вместо | т — a J и — о вместо | а | для областей

ин-

4 З а к . 5 7 6

97


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

 

КА2В2

о

т

 

Фгг(?)

аст — 3 ( т — а) da + ^ е ~ аст- Р <х- а>da +

4аР

 

 

о

 

 

 

 

 

+ t е-°°+ Р (г-°) da = —

— (йе-ах — ае~Рт) =

 

J

 

 

2сф(Р-—0 а 2)

 

т

 

 

 

 

 

 

= /С1 е—« —/Сае—Р \

(4.164)

где /Сх и К г— постоянные. Поэтому взаимная корреляционная функ­ ция выходных сигналов двух детекторов, когда импульсные харак­

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

Р » а ,

(4.165)

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

Если частотная характеристика аппаратуры постоянна во всем диапазоне частот, представляющем интерес, то во всех случаях для

практических целей можно предполагать, что

 

Лх (t) = В8 (0

(4.166)

Фл, /ii (?) = В26 (т).

, (4,167).

98