 
        
        книги из ГПНТБ / Пугачев В.Н. Комбинированные методы определения вероятностных характеристик
.pdf| корреляционной связи Хі * и | ц*. Д ля получения оцен | |
| ки Хіо это пересчитанное значение ошибки | вычитается | |
| из статистического значения | Хі* величины 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
