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