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