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