Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Пугачев В.Н. Комбинированные методы определения вероятностных характеристик

.pdf
Скачиваний:
21
Добавлен:
24.10.2023
Размер:
11.2 Mб
Скачать

корреляционной связи Хі * и

ц*. Д ля получения оцен­

ки Хіо это пересчитанное значение ошибки

вычитается

из статистического значения

Хі* величины X,.

В резуль­

тате оценки Хіо получается более точной, чем статисти­ ческое значение Хі*.

Как следует из (1.22), оценка каждой из компонент вектора X может быть произведена независимо от оцен­

ки других компонент. Поэтому в дальнейшем будем рассматривать только одну из компонент и для просто­ ты записи отбрасывать индекс і.

При этом оценка А*) вероятностной характеристики X

в соответствии с ( 1.22) будет

 

Я0 = А* - V

J sV - f * ) ,

(1.23)

где X* — статистическое значение искомой вероятностной

характеристики исходной

системы; ц и

ц* — точное и

статистическое значения вектора вероятностных харак­ теристик упрощенной системы; R и S — случайные вели­

чина и вектор.

Для практического использования формулы (1.23) необходимо знать корреляционные матрицу-строку K R S и матрицу Kss-

Так как упрощенная система может исследоваться аналитически (по крайней мере для определения век­ тора ц), то в принципе .возможно вычисление аналити­ ческим путем и точного значения корреляционной мат­ рицы Kss■Однако в некоторых случаях определение Kss

аналитическим путем может оказаться затруднитель­ ным. При этом вместо Kss можнб найти ее статистиче­

ское значение K ss,

вычисленное по тем же N

экспери­

ментам, по которым были найдены X* и ц*, т. е. по

формуле

N

 

 

 

Kls =

i r

(1.24)

/=і Исходная система не поддается аналитическому

исследованию, поэтому вместо корреляционной матрицы-

строки

K R S м о ж н о

найти

только

ее статистическое

зна­

чение

Kss по формуле

 

 

 

 

 

N

 

 

 

 

=

і г Е

(Ъ -

Ф - ^)т-

(1 -25)

 

 

;=1

 

 

20

В следующем параграфе будет показана возможность замены K R S и KSS их статистическими значениями KRS

и Kss при вычислении оценки вероятностной характе­

ристики. В соответствии с изложенными выше возмож­ ностями вычисления K R S и Kss помимо оценки Яо будем рассматривать оценки Яоі и Хоъ которые находятся по

соотношениям

яв1= я

* - к ; ^ о**-!*),

(1.26)

я01= я*

- / с;5(^;5)-*(і‘* - і*).

0-27)

Ниже будет показано, что оценки Ао и Яоі, А02 близки

по точности даже при сравнительно небольшом числе экспериментов.

Из изложенного ясно, что для практического исполь­ зования годятся только оценки Аоі и Xoz- Однако Хо

удобна для анализа точности и выигрыша в точности оценок Яоі и Хо2 -

1.3. Точность оценки АоВыигрыш в точности или в числе экспериментов

Найдем точность оценки Ло вероятностной характери­ стики X. Дисперсия Ко оценки Хо в силу (1.23) будет

К0 = М [(Я0 - Я)2];= М [{(Я* -

Я) -

KRSK ^ (* * - *)}*] (1.28)

или

 

 

 

 

К0 = М [(Я* - Я)2]-

2М [(Я* -

Я) KRSK~^ (р.* -

!!)]+

+ М [KRSK ^

(г* -

ц) KRSK ^ (К* - Н*)].

(1.29)

Так как KRSKSJ (р.* — р.) и (Я* — Я) — скалярные вели­

чины, а не матрицы и матрица KJ$ симметричная, то

KRSKss (V-* -

V)-=

-

(ц* -

Ц)1Т =

= (р.*

р.) K ss K RS=

(р.*

р.)

(1.30)

(я* - я) K RSK - <

(р.* -

ц) = V

“ 1(ц *-ц )

(я* - Я).

21

Далее, так как

 

М [(Я *-ЯЛ = 4 - ;Ѵ

 

М [fr*- ? ) ( ? . * - ѵ ) т\ = - ± - K ss,

(1.31)

м[(ѵ* - V)а* - л ] = 4 - * s« = 4

Kl s ’

то из (1.29) с учетом (1.30) и (1.31) следует:

Ко ~ "лГ

^ R R ~ ^RS^' SS Ks/j] ■

(1.32)

 

Введем понятие коэффициента корреляции rRS вели­ чины R и вектора 5, определив его по соотношению

о -зз)

Коэффициент корреляции гя8 представляет собой матрицу-строку с т столбцами. В частном случае при т = 1 rRS представляет собой обычный коэффициент

корреляции двух случайных величин. Из (1.33) и (1.32) имеем

(1.34)

Величина rRSrSR представляет собой квадрат мно­

жественного коэффициента корреляции случайной вели­ чины R и случайного вектора 5 [3].

Дисперсию Ко оценки Ао удобно сравнить с диспер­ сией К статистического значения к* вероятностной ха­

рактеристики. В качестве сравнительной меры точности ?»о и к* будем использовать отношение

тіо = W o .

(1.35)

Отношение т|о представляет собой выигрыш в точно­ сти по дисперсии оценки вероятностной характеристики по сравнению с ее статистическим значением. Следова­ тельно, значение т)о определяет выигрыш от применения предлагаемого метода по сравнению с обычным мето­ дом статистических испытаний.

Так как

K = M [ ( V — X)*\ = KRR/N,

(1.36)

то в силу (1.34) и (1.35)

(1.37)

22

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

Если определяются п вероятностных характеристик

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

^ = 1

^ - 7-

0 = 1 . 2 .......п).

(1.38)

1

rHtS rSR(

 

 

где rR srSR — квадрат множественного коэффициента

корреляции случайной величины Ri и вектора S.

В ряде случаев необходимо находить значения веро­ ятностных характеристик е заданной точностью. При этом за счет применения изложенного метода можно получить выигрыш в числе экспериментов. Найдем зна­ чение этого выигрыша.

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

ских испытаний для определения г-й компоненты век­

тора вероятностных характеристик

потребовалось бы

Ni экспериментов, где

(1.39)

Ni = KR R IDi.

ча чі

 

Поскольку все п компонент находятся из одних и тех

же экспериментов, то необходимое число эксперимен­ тов N при обычном методе статистических испытаний

будет

іѴ= max (KR R ID i).

(1.40)

Естественно, что при числе экспериментов N диспер­ сии Кі статистических значений вероятностных харак­

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

(1.41)

При применении метода коррелированных процессов для определения Ей компоненты потребовалось бы N i0

экспериментов, где

Ni0:

KKRi ~ X s

(1.42)

 

А

 

23

Необходимое же число экспериментов для получения всех компонент равно

( 1 ~

r Rt S r S R )

(1.43)

NB— max

Dl

 

 

Выигрыш в числе экспериментов при применении метода коррелированных процессов по сравнению с чи­ стым методом статистических испытаний составит

шах ( К р в

/ D i )

 

_ N ______________ i

_________________

(1.44)

AT- - max [/С^л (1 —

г ^ 3 r S R ) / D i ] '

 

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

£=•1/(1—rRSrSR). (1-45)

В данном случае выигрыш § в числе экспериментов равен просто выигрышу rjo в точности по дисперсии, т. е.

£=гПо.

(1-46)

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

Di — K R R jMo.

(1.47)

При этом в силу (1.38) и (1.44)

е = шіп-ч,-0.

(1.48)

Таким образом, в обоих рассмотренных более или менее общих случаях величины выигрышей по точности определяют и выигрыш в числе экспериментов.

Поскольку квадрат множественного коэффициента корреляции определяет выигрыш по точности и, следо­ вательно, в частных случаях — выигрыш в числе экспе­ риментов, то ниже приведены значения квадрата мно­ жественного коэффициента корреляции r ^ s R для т = 1, 2, 3.

1. m = U

Г

f

Т

RS'

(1.49)

'

RS1S R ------

 

 

24

Зависимость выигрыша г|о от

|rws|

приведена

на рис. 1.1.

2. т 2 :

 

 

2rRS, rRSj rs,sa + гяэа

 

r RSrSR :

rRS,

(1.50)

 

 

 

rs,sa

 

 

 

 

3. m = 3:

 

 

 

 

 

 

 

 

 

 

 

 

r RSrSR:

''З Д

''s.S a ~~

r SaS3 +

 

X

 

2 r S,Sa ''x.Sa r SaS3

X [ry?.s'10

r s , s )

(1

 

r s , s )

”b r Rs 30

r s vs.)

2 r RSir RS1( r S,S.

rS1s/s,S3>)

^ rRS. r RS, (rS,S, r s . s / s . s )

 

2 r RS,r RS3(rV' STSs,s,

rS1S2r. SIS3.)]•

(1.51)

Покажем, что с увеличением числа вероятностных характеристик, рассматриваемых в упрощенной системе, значение выигрыша т)о растет. Для ,

этого достаточно показать, что если £ к т вероятностным характеристи- 8 кам прибавить (т + 1)-ю, то квад­

рат множественного коэффициента rRSrSR увеличивается, т. е. разность 6

Д=

(rRSrSR )m + t( Т RS?“R ) т

 

(1.52)

5

 

 

 

 

 

 

становится больше нуля. Здесь ин-

3

 

 

дексы

т + 1

и т при rRSrsR указы­

 

 

 

вают,

что

значения квадрата

мно­

 

 

 

жественного

коэффициента

 

взяты

1О 0,2 0,0- 0,6 0,8\t=RSI

при т + 1 и т

вероятностных ^харак­

 

 

 

теристиках

в упрощенной системе.

Рис. 1.1. Зависимость

Обозначая

через

Кт и Кт+и А т и

выигрыша

% от rRS.

Ат+1

соответствующие

значения

 

с

учетом

корреляционных

матриц

K s s

и K R S , можем

(1.33)

записать,

что

 

 

 

 

 

 

 

 

 

 

■\Ат+1КГ 1

АЛт

АтК~'АТ ].

(1.53)

Очевидны

КRR

т + 1

т + 1

m m *

 

'

следующие

связи между матрицами Ат+\

И А т> Кт+1 И Кт

Ат+1 — \\Ат К RS

||.

 

(1.54)

 

 

 

 

 

 

 

 

 

Кт +і

Кт В п1

 

 

(1.55)

 

 

 

 

 

 

 

 

 

В т К

Sm+ l5m+ l

25

где

Si" . K s

s

\\.

(1.56)

Bm = \\Ks

 

m + l 1

m- fl m

 

 

Применяя к матрице

Km+1 формулу

обращения [40]

блочной матрицы, найдем

 

 

 

 

 

Х - І ^ ±

К - ' В Г BmK~ 1

— L K ~'B T

 

 

m 1 И

m m

m

 

Cl

m m

(1.57)

ІГ\

 

 

 

 

 

1

m + l

------- BmK-'

 

 

 

 

 

 

а

rn

 

 

 

 

ll

 

где а — величина, определяемая выражением

 

 

а= К<' с

с

в т

к ~

х в Т

 

(1.58)

 

 

m +1

m + l

 

т

т

 

 

Величина а положительна, так как 1представляет

собой диагональный элемент матрицы, обратной кор­

реляционной матрице К т +

1.

 

 

в (1.53), после

поблоч­

Подставляя (1.54)

 

и

(1:57)

ного перемножения

матриц

получим,

что

 

 

 

аКRR

[.АтК~'ВТ ВтК~]Ат

 

 

 

1

 

m

m

 

 

m

m

 

 

 

B m!\-'A T К

 

AmK~'BTl

К

 

 

+

4 , ,

].

(1.59)

m m iw

 

m

m

m

RSm +1 ‘

"'®ш+1

 

 

 

 

 

Учитывая, что Лт /(

-1 пТ

 

и Дт /<

 

 

Л

— величины, а

 

ß

 

 

 

 

 

 

m m

 

 

 

m m

 

 

 

не матрицы и матрица

/С~ симметричная, можно

напи­

сать следующее соотношение:

 

 

 

 

 

 

 

 

АтК~'Вт= (АтК~'Вт )г= В тІ Г ' А т .

 

 

m m

 

 

 

m m

 

 

 

 

m m

 

 

Отсюда с учетом (1.59) имеем

 

 

 

 

 

 

 

 

 

,(ЛТ к~'в„

« s

 

 

) 2.

 

(1.60)

 

'

 

m

m

 

 

 

 

 

" аКRR

 

 

 

 

 

 

 

m + l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как а положительна, то

 

 

 

 

 

 

 

 

 

 

 

Д $г0.

 

 

 

 

 

 

(1.61)

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

26

1.4.Точность оценок Яоі и Я02

Впредыдущем параграфе была определена точ­ ность оценки Я0 и найден выигрыш в точности этой

оценки по сравнению с точностью статистического зна­ чения Я*. Однако оценкой Я0, как указывалось в § 1.3, практически пользоваться нельзя и вместо нее прихо­ дится применять оценки Я0і или Яог-

В данном параграфе будет показано, что при доста­ точно большом N точность оценок Яоі и Я02 практически

равна точности оценки Яо. Поэтому все соотношения для точности и выигрыша в точности или в числе экспери­ ментов, приведенные в предыдущем параграфе по отно­ шению к оценке Яо, можно будет считать справедливыми и по отношению к оценкам Яоі и Я02 при достаточно большом N.

Не нарушая общности результата, в настоящем па­

раграфе

рассмотрим случай

т = 1, т. е. когда S

и р —

не векторы, а величины.

 

 

 

 

Перейдем к анализу точности оценки Яоі. В соответ­

ствии с

(1.3), (1.4),

(1.25)

и (1.26)

 

 

 

 

 

 

N

 

 

 

ДЯ01=

Я01- Я

=

- ^ ^ № - Я ) -

 

 

 

 

 

1 = 1

 

 

 

1 N

 

 

 

N

 

 

— лЙг-SiRj~Я) (Sj~^Е (Sl” ^“

 

 

/■=1

 

 

 

/=1

 

 

1 N

 

N

 

N

 

 

- ж - 2

(Ri - Я)

(S, -

р ),£ (Si - р).

(1.62)

 

1 = 1

 

/ =

1

1 = 1

 

Точность оценки Яоі будем определять математическим ожиданием квадрата іЛЯоі, т. е. величиной

/Соі= Л4і[(АЛ.0і) 2].

(1.63)

Из (1.62) и (1.63) видно, что в выражение

для /Соі

будут входить центральные моменты до шестого поряд­ ка включительно от случайных величин R и S. Здесь

и в дальнейшем будем предполагать, что все необходи­ мые нам моменты величин R и S существуют. При вы­ числении Коі следует учитывать, что величины R и S

с различными индексами независимы в силу независи­ мости различных экспериментов.

27

После ряда преобразований можно получить, что

К

N RR

КRS^SS ) +

di. I

I

A

I

(1.64)

N2“ Г

hiа " Г

N*

'N5

Здесь

величины Лгвыражаются через моменты /(

9 :

 

 

 

 

 

 

 

 

 

 

0-65)

Выигрыш в точности оценки Яоі по сравнению с точ-

ностью статистического значения Я* будет

 

 

 

 

 

 

 

 

К

 

 

 

 

 

 

 

 

 

 

^ 01

 

 

 

 

 

_

 

 

 

1

 

 

 

 

(1.66)

 

 

 

9

В,

, Вг , В3

Ri

 

 

 

 

 

 

 

 

(* — rRS) +

N +ДГ2 + jV3 + N*

 

 

 

где

 

 

 

 

 

 

 

 

 

(1.67)

 

 

 

Вг =

А-іІКщ?

 

 

 

Отсюда получаем, что

 

 

 

 

 

 

 

 

 

lim т|01 =

1/(1

- r \ s)

 

 

 

(1.68)

 

 

УѴ-»со

 

 

 

 

 

 

 

или с учетом

(1.37)

lim т]01 =

V

 

 

 

(1.69)

 

 

 

 

 

 

 

 

 

N-*оо

 

 

 

 

 

 

Из сравнения

(1.37)

и (1.66), а также из

(1.69) следует,

что оценки Яоі и Яо

при

достаточно

большом

N близки

по точности.

Найдем ориентировочно те значения N, начиная с ко­

торых оценки Яо и Яоі практически совпадают по точно­ сти. Сравнение оценок Яо и Яоі по точности удобнее всего вести в относительных величинах, например путем срав­ нения выигрышей г|о и т)оі-

Рассмотрим два частных случая. В первом случае будем считать, что вероятностные характеристики Я и ц представляют собой математические ожидания, а вели­ чины К и 5 подчинены нормальному закону распреде­ ления. Во втором случае в качестве Я и ц будем рас­ сматривать дисперсии некоторых случайных величин X

и У, подчиняющихся нормальному закону распреде­ ления.

28

Для первого случая из (1.37) и (1.66) имеем

 

71 о = 1/ ( 1 -r2RS),

(1.70)

 

1

 

 

1

(1.71)

1

„2

4

 

г RS

 

.

1

 

 

1 +

2 J

(гЯ.у1 ж

 

 

 

 

1=1

 

 

где

 

 

 

 

(1.72)

C

/ W

= W

- 4

) -

 

Так как rjo и г|оі зависят только от rRS и іѴ, то, исклю­ чая из (1.70) и (1.71) величину rRS, можно написать, что

Т)01= 'ф(тіо; 77).

(1-73)

На рис. 1.2 приведена зависимость т]оі от тіо для раз­ личных N, полученная методом статистических испыта­ ний. Как видно из рисунка, уже при N ^ 3 0 значения rjo

и rjoi близки, т. е. оценки Ло и Яоі близки по точности.

Рис. 1.2. Зависимость Поіот т)0

Рис.

1.3. Зависимость т)0і от Цо

для

различных N в случае,

для

различных N в случае,

і когда

А, и р — математические

когда А, и р — дисперсии,

ожидания.

Во втором случае цо и rjoi будут зависеть только от коэффициента корреляции гх у случайных величин X и У

и числа экспериментов N. Причем, как будет

показано

в § 1.8,

 

7і01= 1 / ( 1 - 4 ) .

(1.74)

29

Соседние файлы в папке книги из ГПНТБ