Файл: Уриг, Р. Статистические методы в физике ядерных реакторов.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 16.10.2024
Просмотров: 152
Скачиваний: 0
В работе [10] найдено выражение для вероятности Р того, что по грешности определения амплитуды и фазы передаточной функции будут меньше некоторого значения, т. е.
Р — Prob |# 1 - |Я | |
< s in e |
и |
|
\ н \ |
|
|
|
|
I — У% |
|
к[Г |
|
|
(6.84) |
|
1 |
2 |
2 |
|
1 — уху cos |
|
е |
Можно преобразовать выражение (6.84), чтобы получить число сте пеней свободы
/е = -— |
2 In (1 —Р) |
(6.85) |
|
1 2 |
|||
In |
1 — Уху |
|
|
1 — Уху C O S 2 |
6 |
||
|
которое очень полезно для планирования эксперимента. К сожале нию, необходимо предположить известным значение функции коге рентности у1у (со), так как она непосредственно связана с качест вом измерения. Обычно выбирается низкая оценка yly, чтобы до стигнуть требуемой точности измерения.
Можно найти погрешность е подстановкой тригонометрического равенства
|
|
cos2 е = 1 — sin2 е |
(6.86) |
в выражение (6.34) и решением его относительно е, т. е. |
|
||
е л* sine = |
1 |
1 |
(6.87) |
Уху (о>) |
И — УдУ®)] ( 1 _ р ) 2 / А |
||
|
|
|
Это означает, что с вероятностью Р справедливы три неравенства (смысл всех неравенств один и тот же):
||Я | — |Я | | < |
| Д | sine, |
(6.88) |
|
Я | (1 —sin е) < |
| Я | < |
| Я | (1 + sine); |
(6.89) |
| й | < |
| Я | < |
if ii |
(6.90) |
1 +sm 8 |
|
1—Sin 8 |
|
10— е | < |
10 | < |
| 0 -be [. |
(6.91) |
Эти соотношения графически изображены на рис. 6.12, где с вероят
ностью Р оценки Я (со) попадают в заштрихованную область. Для значений е < 0,2 в выражениях (6.88), (6.89) и (6.90) можно исполь зовать
sin е « е . |
(6.92) |
168
В работах [12] и [13] также исследовались погрешности как откло нения от экспериментальных значений. Опять предполагалось, что измеренная величина несмещенная. Для общности рассматрива лась система с многими независимыми и коррелированными входами. Мы рассмотрим общее выражение для погрешности и ограничимся случаями одного и двух входов.
Рассмотрим систему, изображенную на рис. 4.18 с q входами, где передаточные функции # г (со) даны взаимнокорреляционным соот ношением из уравнения
(4.211)
G | » = 2 |
Hh (<o)Gtk(a) |
£ |
||
|
k = 1 |
|
(6.93) |
° -г |
|
|
|
g |
|
|
|
|
|
5 |
для |
случая, |
когда |
входы |
сз |
t |
||||
коррелированы. Для случая, |
|
|||
когда |
входы |
некогерентны, |
|
|
соотношение |
(6.93) |
можно |
|
упростить и привести к виду |
0 |
<Действительная ось |
|
|||||||
Сг!/М = |
ЯгМ <ЗгДсй). |
(6.94) |
|
|
||||||
|
|
|
|
|||||||
|
гг |
^ |
случаев |
(1 |
|
Рис. 6.12. Оценка статистической погреш- |
||||
Для |
ОООИХ |
|
|
ности в измерениях передаточ- |
||||||
—а)—процентные доверитель- |
|
ной функции [11]. |
|
|||||||
ные |
пределы |
для |
истинной |
|
|
|
|
|||
амплитуды | Hi (со) | |
и истинного фазового |
угла 0; (со) — для каж |
||||||||
дой |
частоты f |
и каждой системы (i |
= 1, |
2, ... q) даются в виде |
||||||
|
|
[| Hi (со) | - |
Д, (со)] < |
IН, (со) | < |
[| Hi (со) | + Д,- (со)], |
(6.95) |
||||
|
|
[ 0 £ |
(со) — |
6 f (со)] < |
е |
г ( с о ) < [ 0 г (со) |
+ 6 г (со)], |
( 6 . 9 6 ) |
где выражения для погрешности определения амплитуды (радиаль
ной) Дг (со) и фазы 6г (со) даются для каждого рассматриваемого слу чая. Эти соотношения показаны схематично на рис. 6.13, где с (1 —
— а) — доверительной вероятностью значения Ht (со) нахо дятся внутри круга.
Множество некогерентных входов. Для случая, когда все входы некоррелированные, Дг- (со) и 6г (со) даются в виде
Дг И = |
k— 2 |
1/2 |
(6.97) |
Ои (СО) |
|||
|
|
( Р п „ л., а) [ 1 - Ъ с И ] | ^ М |
|
|
$i (со) = arc sin [Дг (со)/1 Я г (со) I], |
(6.98) |
|
где k — число статистических степеней свободы для каждой |
оцен |
||
ки спектральной |
плотности; (F„ti „2>а) — 100%-ные точки ^-распре |
||
деления с пх = 2 и /г2 = |
k — 2 степенями свободы; Ъи (со), Gyy (со)— |
169
измеренные спектральные плотности г'-го входного и выходного сиг
налов соответственно; yh — измеренное значение функции, коге рентности между выходом у (t) и г'-м входом xt (г1), т.е.
У & И = I G ,ji И \2/ G u И ° и и И - |
( 6 - 9 9 ) |
Для этого случая точность увеличивается с ростом числа степе ней свободы для каждой спектральной оценки, увеличением (при ближением к единице) фактора когерентности у}я (со) и увеличением
отношения Gyy (cd)/Gh (со), которое представляет собой квадрат ам плитуды передаточной функции Ht (со).
Рис. 6.13. Статистические ошибки в измерениях пере даточной функции с многими входными сигналами.
Один вход. Для случая, когда имеется только один входной сиг нал, все соотношения (6.94)—(6.99) применимы с подстановкой ин декса i = l .
Множество когерентных входов. Когда входные сигналы коррелированы, оценка статистических ошибок усложняется. Выраже ния (6.95) и (6.96) остаются справедливыми, за исключением того, что величина Д* (со) имеет вид
_ 2q_ |
По.а) 1 -Ту.««Р) |
| 1/2 |
(6. 100) |
(Fn |
|||
п — 2<7 |
1 -Т?.*(со) |
Gil (со) I |
’ |
где q — число входов. Степени свободы ^-распределения даны при пх = 2q и /г2 — п — 2q\ yl.x (со) — измеряемая величина много канальной функции когерентности между выходом у (t) и всеми
измеренными входами; у].х — измеряемая величина многоканаль ной функции когерентности между входом л: (/) и другими измеряе мыми входами, исключая хг (/). Множество функций когерентности дается выражением
yf-x ( ® ) = 1 — [ G Hсо)( G * J ~ ~ 1 , |
(6. 101) |
170
где GxX обозначает i-ю диагональ обратной матрицы |
[G** (со)]-1 |
||
и матрица GAA(со) дается выражением |
|
|
|
' G n N |
G J 2 (co) . . . |
G lg (со )' |
|
^21 (®) |
^22 (®) • ■• |
Gnq (ft)) |
(6. 102) |
G q i H |
G q 2 (a>) . . . |
G 5 ? ( co) . |
|
Множество функций |
когерентности yj).A.(со) между выходом у (t) |
|
и всеми входами дается в виде |
|
|
y l x |
(со) = 1 - [Gyy (со) G'ijxx(со)] ~ \ |
(6.103) |
где Gyxx (со) есть i-ый диагональный элемент обратной расширенной матрицы [Gj,** (со)]-1. Расширенная матрица G ухх (со) имеет вид
’ G y y (со) |
G y l (со) |
G y2 (со) |
. . . G yq (со) |
|
|
G i y И |
О ц (со) |
G12 (со) |
. . . |
G lq (со) |
|
Gyxx (®) = G oy (со) |
С 21 (со) |
G22 (со) |
. . . |
G 2q (со) |
(6.104) |
■Gg„(С О ) |
G g l (со )G 9 2 ( co) |
. . . |
G qq (■ ) _ |
|
|
Выражение для б; (со) то же, что и в соотношении (6.98). |
|
||||
Два когерентных входа. |
Если входные сигналы системы с двумя |
входами коррелированы, то матрица (6.102) превращается в матрицу
|
Gii И |
G12 ( со) |
|
G x x (со) = |
G22 (со) |
||
|
G 2i ( со) |
||
и обратная матрица есть |
|
|
|
[оя* и ] - ' = |
G22‘ |
(^) |
G (со) |
4 |
|
G u ( со) |
|
|
— G 2i ( co) |
||
где А — детерминант матрицы Gxx (со) |
|
||
Gu (со) G12 (со) |
|
G22 ( со) — | G12 (со) |2]. |
|
А = |
[G1X(со) |
G%i (со) G22 (со)
(6.105)
(6.106)
(6.107)
Следовательно, первый и второй диагональные элементы обратной матрицы [G** (со)]-1 суть
Gxxl( u ) - ° ^ |
G 22 (со) |
(6.108) |
|
G u (co) G 22(co) - | G 12(co) P ’ |
|||
Д |
|
||
G n (со) ________________ G u (со)_______________ |
(6.109) |
||
Gxx, (со) ■ |
G n (со) G 22 (со) — | G 12 (со) р |
||
A |
|
171