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

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

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

.лизуемые оценки Лю' и W будут иметь вид

 

 

К*22 1

V 1-- Г*

(— Я2 <2) +

С

<1,2>).

X' 20

_2L n*<2) 1 ^2

 

 

 

 

(4.54)

 

/

 

 

 

 

1 \ [

К*„

 

 

Г*

X ] {l) +

X ] ’

( 1. 2) ,

 

2 Г

К * и

1

(Я 1 — г * г

" - Z>),

 

 

 

 

где

 

 

 

 

 

1,2

 

 

 

 

 

 

 

1

 

 

 

 

 

( 1, 2)

V* п(1.2

 

 

 

 

 

УѴ*,,2

)

 

 

 

 

Л‘

 

 

 

 

 

 

 

 

 

1=1

 

(4.55)

 

 

 

 

 

 

"*1,2

 

 

 

 

-)** (1.2) __ 1

 

 

 

 

 

у \ рС.2)

 

 

 

 

 

Я2

 

 

'

 

 

 

 

 

 

 

 

/=1

 

 

В

соответствии

с (4.50)

числа экспериментов N і* и

Д^2* не зависят

от результатов экспериментов:

 

 

 

 

N \ = N \ = 4 ^ ,

 

(4.56)

а число Л7і,2* определяется из соотношения

 

 

 

 

 

 

=

“ г* а>

(4-57)

где

 

 

 

 

 

 

 

 

к .

(4.58)

к*„к*

 

Будем считать, что Ки* и К2 2 *, входящие в (4.54),

вычисляются по экспериментам типа (1) и (2), а г*, вхо­ дящее в (4.54) и (4.57), — по экспериментам типа (1, 2).

Организация экспериментов и обработка их резуль­ татов в данном случае могут быть следующими.

Сначала проводится по Мо/2 экспериментов типа (1)

и (2), т. е. эксперименты в каждой из двух точек с неза­ висимыми последовательностями воздействий Wt и W2 соответственно. Затем с последовательностью Wі>2 про­

водятся эксперименты типа (1, 2) в первой и второй точ­ ках, ■§ ходе которых все время вычисляется статистич^-

160

ское значение г* коэффициента корреляции. Эти экспе­

рименты прекращаются, как только фактическое число экспериментов станет больше расчетного. После оконча­

ния

экспериментов производится определение

оценок

W

и W по формулам (4.54) и (4.55).

 

Покажем, что даже при сравнительно небольшом чис­

ле

экспериментов оценки Аю' и Аао' близки по

точности

к оценкам Аю и Аго-

 

В качестве мер точности оценок Аю' и Аа/ примем ма­ тематические ожидания квадратов отклонений от их точ­ ных значений Хі и А2:

 

 

/ и ^ м к ^ о - я . ) 2],

(4.59)

 

 

К ' 2= Л4[(А'20- А 2) г].

 

 

 

Если К і

и К г

будут

близки к К і и Кг,

то это озна­

чает, что замена

К н , Кгг

и г их статистическими значе­

ниями при

вычислении

оценок допустима,

так как не

приводит практически к ухудшению точности оценок. Определим относительное увеличение дисперсии оценок:

 

К’і - К ,

'* -к .

(4.60)

 

 

Кг

Кг

 

 

Для рассматриваемого случая

 

 

Лі

м„

Л ',=

м 0 ■

(4.61)

Вычислим значения

дисперсий

К і

и К г , входящих

в выражения для е. С этой целью воспользуемся сов­ местно аналитическим методом и методом статистиче­ ских испытаний.

Будем рассматривать случай, когда значения R i и R 2 величины R , математические ожидания которых равны

искомым значениям Аі и Аг вероятностной характеристи­ ки, подчиняются нормальному закону распределения. В этом случае [3] статистические значения математиче­ ских ожиданий и корреляционных моментов независимы.

Следовательно, независимы между собой А*(1) и К*и ,

а также Я*(2) и К*г2. Напомним, что значения К*п и

К*.2 находятся соответственно по экспериментам типа(1)

и (2). С учетом независимости экспериментов различного

типа получаем, что все

четыре указанные

величины не­

зависимы между собой

и являются независимыми по от­

ношению к значениям А**(1,2), а” (і,2) и г*,

определяемым

11—288

 

161

по экспериментам типа (1,2). Статистический коэффициент корреляции г* находится через статистические корреля­

ционные моменты, определяемые по экспериментам типа (1,2). Поэтому, если бы число экспериментов типа (1,2) не определялось через статистический коэффициент корре­

ляции г*, то значения Я**(1,2)

и Я**(І,2)

не

зависели

бы

от г*. Зависимость Я**(|,2) и Я**(1,2) от

г* вызвана только

лишь зависимостью числа экспериментов JV*12 типа (1,2) от

4*

*#/10)

**/10)

 

 

 

г*, по которым найдены Я,

1' '

и 12 1 ’ .

 

 

Найдем

сначала

Кі-

Обозначим

через

фі(/Сц*,

К2 2 *,

г*) условное математическое ожидание величины (Яю'— —Лі)2, т. е. ее математическое ожидание при определен­

ном значении Ки*,

К2 2 * и г*. Так как

 

 

 

Я',

Я, — - у - (я, ( 1 - -

/1,) -|- ■ (Я. • ( 1. 2) я,)

 

2 /

к * „

 

[(я 2- (,,2). _

я £) - ( я ;

( 2 )

- Я,)],

К*ы 1+

Y'i — r'

 

(4.62)

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

Л 4 [ ( я ;( , ) - Я , ) Ч

м [(я ;(|)- я2)ч

Л4[(Я, '(1,2 )

ж [ ( я ; * ^ - я 2) ^

К

и

II

 

/V*,

М 0

 

К

г г

2

К и

*2

 

М0 ’

 

 

 

2К„

 

^*1.2

‘ /И0К і- /- * 2 ’

]

Кц

2

2.Кц

^

2И„J/ 1—/-*2

м [(я**(|,2) -

 

 

 

К,2

 

_

2г Ѵ к л й і

Я,) (Я**(І,2) — Я2)1 = JV*I

 

М

01/' г —Г* 2

получаем

 

 

 

 

 

 

 

 

(/< *„,

/С*22,

г*)

=

К и

2

V 1— г* 2

 

 

 

 

 

 

 

 

 

I 1

К * и

к і2

 

 

 

 

 

 

2

К і , Я * 2 2 | Л — Г * 2 (1 + Ѵ\ — Г * 2)

 

у

U K22_ ________ ___________

(4.63)

/ К*"

К’*

 

Г*Г

 

 

 

 

К и

К * 2 2

Ѵ \ — г* 2( 1 + У" 1— Г* г)

 

162

Отметим,

что

если

в

(4.62)

Ки*,

Kzz* и г*

заменить

 

 

на Кц,

Kzz и г, то оценку Яю' необходимо заменить на Ям-

 

 

Поэтому в силу (4.61)

*

и

г. )= ^*. и .

 

( 4 . 6 4 )

 

 

 

 

 

«

М

 

Имея условное математическое ожидание -фі(/Си*,

 

 

Kzz*,

г*),

можем написать выражение для Кі':

 

 

 

 

=

 

J

5 5

ф І.

С

Ѵ

 

к\ =

 

К

*

мdK*udK*S)) / ( г 2dr*,

 

 

 

(К*и))/ ■ // ,2 (

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 4 . 6 5 )

где /і(/Сц*), fz(Kzz*), H r* ) — дифференциальные

законы

 

 

распределения величин Ки*, Kzz* и г*.

 

 

 

 

 

 

 

Предполагая число экспериментов Л40/2 достаточно

 

 

большим, можем в пределах реально возможных значе­

 

 

ний

Ки*— Ки

и

Kzz*— Kzz разложить

 

функцию

 

фі(/Си*,

 

 

Kzz \

г*)

в ряд по степеням этих отклонений, ограничив­

 

 

шись членами второго порядка. Тогда, учитывая, что

 

 

 

 

 

м

 

 

М \ К * ц - К ц \ = О,

 

 

 

 

 

 

 

 

 

 

К*и Ки

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К

и

 

 

 

 

 

 

л и

 

 

 

 

 

получим

 

 

 

 

[ 1

+

і

/

(

г

dr**],)

/

(

г( 4* . )6

6

)

где

 

 

К’г-

 

КЛ4„и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X (г*)

Kl—Г*2 (1 + K l—г*2)

 

(г* ~ г) + Ж

(2г* ~ г)]-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.67)

 

 

Из

(4.60) — (4.62)

находим

выражение

для

 

относи­

 

 

тельного увеличения дисперсии оценок:

 

 

 

( 4 . 6

8

)

 

 

 

 

 

 

*=

ti(r*)f(r*)dr*.

 

 

 

 

Сложность определения значения е аналитическим путем по данной формуле обусловлена тем, что практи­ чески невозможно найти выражение для закона распре­ деления f(r*) статистического значения коэффициента

корреляции. При заданном числе экспериментов в {3] имеется выражение для дифференциального закона рас­ пределения величины г*, однако для вычисления е он не

годится, так как в нашем случае число экспериментов

11*

163

N i,2* не задано, а определяется через случайное значе­ ние самого же г* по формуле (4.57). Поэтому дальше

для определения е был применен метод статистических испытаний.

При статистических испытаниях вырабатывались по­

следовательности двух случайных

нормальных коррели­

 

 

 

 

рованных величин, вычис­

 

 

 

 

лялась

 

последователь­

 

 

 

 

ность статистических зна­

 

 

 

 

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

кор­

 

 

 

 

реляции

и находилось то

 

 

 

 

значение г*, для которого

 

 

 

 

выполняется

 

условие

 

 

 

 

(4.57). Далее по формуле

 

 

 

 

(4.67)

определялось

зна­

 

 

 

 

чение

 

функции

%(/■*)•

 

 

 

 

Среднее

значение

функ­

 

 

 

 

ции

%(г*)

по

п

таким

 

 

 

 

экспериментам

дает

ста­

 

 

 

 

тистическое

значение ве­

 

 

 

 

личины е. Поскольку при

 

 

 

 

г*,

близком

к

единице,

 

 

 

 

число

N і,2*

может

 

быть

 

 

 

 

очень малым, то значение

 

 

 

 

УѴі,2* выбиралось из усло­

 

 

 

 

вия (4.57), но не менее 5.

Рис.

4.2.

Относительное

увеличе­

На

рис. 4.2

приведена

полученная зависимость е

ние е дисперсии оценки в зависи­

от г

для

Мо, равного 50,

мости

от

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

корреля­

ции г

для

различных М 0.

 

200 и 500. Полученные за­

 

 

 

 

висимости

показывают,

что значение е является небольшим. Таким образом, за­ мена дисперсий и коэффициента корреляции их стати­ стическими значениями в формулах для оценок вполне

допустима.

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

основания.

В заключение данного параграфа покажем на про­ стейшем примере, как зависит эффективность предла­

164

гаемого метода от взаимного расположения взятых то­ чек (значений параметров), для которых находятся искомые вероятностные характеристики. Пусть искомая вероятностная характеристика модели системы пред­ ставляет собой дисперсию установившегося процесса на выходе модели системы с передаточной функцией К(р) = =>k/(l +рТ) при воздействии на нее стационарного бело­

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

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

значение квадрата процесса У

на

выходе системы:

 

R = Y\

 

 

 

Поскольку У нормален,то

 

 

 

где гу у

— коэффициент корреляции значений

установив

щегося

процесса для двух значений параметров при одном

и том же воздействии.

 

 

 

Легко получить, что

 

 

 

 

_ 2 ] П \ Т \

 

 

 

rYlv — т1+

Т2

 

где Ті и Т2— значения параметра

Т в двух точках.

Отметим, что в данном случае коэффициент корреля­

ции гѵ ѵ , а следовательно, и г

не

зависит

от взятых

IІІ2

 

 

 

значений ki и k2, т.е. выигрыш не зависит от того, какие значения параметра k рассматриваются.

Пусть для определенности значения параметров Уі и

У2 таковы, что

Т2>Тр, тогда, используя формулу (4.50),

получим

 

-

Ло= 1 + Уі/Уг-

Отсюда, например, видно, что если даже Т2 в два раза больше Ті, то выигрыш т)о= 1,5.

165

4.4. Подоптимальный метод определения

 

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

в трех точках

 

с требуемой точностью

 

 

При произвольном выборе чисел

N i , N i th , . . . , N ti2

i

экспериментов различного типа, но заданном общем чис­ ле экспериментов N, только лишь за счет оптимальной

обработки результатов экспериментов можно получить выигрыш в точности или в крайнем случае ту же точ­ ность, что и при обычном определении вероятностных характеристик, когда проводятся эксперименты либо типа (і), либо типа (1, 2 ,...,/). Отсюда вытекает, что,

задаваясь желаемой точностью определения вероятност­ ных характеристик и более или менее произвольным пе­ рераспределением чисел N i , Ni'h, ■., N i p ......i различных

типов экспериментов, можно получить выигрыш в об­ щем числе N экспериментов или в крайнем случае не

проигрыш. Этот выигрыш будет тем значительнее, чем ближе будет перераспределение указанных чисел к опти­ мальному.

Сказанное, а также результаты, полученные при оптимальном решении задачи в случае двух точек, позво­ ляют построить подоптимальное решение для случая трех точек.

Пусть требуемая точность определения оценок Лоі, ког и коз вероятностной характеристики к в трех точках за­ дается необходимым числом экспериментов Мй в каждой

из трех точек при обычной организации экспериментов и независимом определении вероятностных характери­ стик.

Из анализа случая двух точек вытекает, что при этих предположениях рационально числа экспериментов типа (1), (2) и (3) брать одинаковыми, независимыми от ко­ эффициентов корреляции и равными числу М0, поделен­

ному на число точек.

N 2 и N s

выберем рав­

Таким образом, значения N t,

ными

 

 

Ml==N2 = N

^ ^

(4.69)

Так же как и в случае двух точек, числа эксперимен­ тов типа (1,2), (1,3) и (2,3) выберем такими, чтобы при равенстве одного из коэффициентов корреляции единице соответствующее число экспериментов было равно нулю.

166

П р и м е м

 

 

N.i,k :

М„

(4.70)

т Ч і

где а — некоторая величина, значение которой выбирает­

ся исходя из требуемой точности получения оценок Мо, Л20 и ^зо.

Число экспериментов типа (1, 2, 3) возьмем равным

нѵлю, т. е.

 

Ni,2,3= 0.

(4.71)

Общее число экспериментов N будет

N = N1+ Na + Nl + 2(Nl, t + N1' l + Nt, l) =

1 + (1 —

1 \

' 1 2 +'

г13 + г\23

(4.72)

 

 

Матрицы Qi и Qi,h в соответствии с § 4.2 имеют вид:

Qi.

Qi.s

 

 

0

0

 

 

0

Ü

0

Q ,=

0

0 0

 

Q2=

0

1

 

IT-0

 

 

0

0

 

 

0

Д22

 

0

 

 

0

0

 

 

 

0

0

0

 

 

 

 

Q.

 

0

0

0

 

 

 

 

 

0

0

1

 

 

 

 

 

 

К3г

 

 

 

1

1

 

 

 

1

 

Г12

1-

-Г2и

 

 

 

УМі^22

1 —

г 2

1

 

 

 

1

4 2

 

Г12

 

 

1

1

 

 

V'KuRti

,212

^22 1-

4

 

 

і-г ? з

 

 

 

 

 

 

l -

V К^К33

1„—

 

 

Ка

1— Г13

 

167

0

 

0

 

 

0

 

 

 

1

 

1

 

1

г 23

 

 

 

 

 

0

22

1 -

4

 

 

 

 

я

 

- ^ 22^33

1 ;—

4 > 3

 

 

 

 

 

0

1

 

Г23

 

1

1

 

 

 

 

 

К з з 1

- 4

 

 

 

 

1 —

г 23

 

 

 

 

 

 

 

где Гш— коэффициенты

корреляции случайных

величин

Ri и Rh в экспериментах типа

(г, k),

 

 

 

Оценка L0 вектора L значений вероятностных харак-

:еристик в трех точках в соответствии с (4.9) будет

L0 = P ;' {Q1L*1+ Q 2L*2 + Q3L*3 +

+

а [(1 - г 21 2 )

QU2L*„2 + (1 -

4 ) QutL\

з +

 

+

( 1 - 4 ) Q 2,3L*2, 3]},

(4.73)

где

 

 

 

 

 

 

Ро— Qi + Q2

+ Q3

+

 

+ 4(1

- 4 ) < Э ..2 + (1 - ^)Q..3 +

(1 - 4 ) Q

2. 3J. (4.74)

Корреляционная матрица Ко оценки L0 определяется

по соотношению (4.11) и равна

 

 

 

 

 

440

р

 

(4.75)

 

 

3

0

 

 

 

 

 

Дисперсии оценок Кі0 (і = 1, 2, 3), равные диагональ­ ным элементам матрицы Ко, будут

(1 +2<X)2- 4 « 3

д . = 4 ’ 3

Л40

кК.22 Q (1 + 2a)*-rf3a» і\2 М0

кК33 о (1+2

где

А = (1 + 2а)3 — 2ri;irI3rs3<x3--

~ ( 4 + 4 + 4 ) (I + 2а) а2

(4.76)

(4.77)

168

Так как в соответствии с требованиями по точности получения оценок

Kis^Kii/M0,

(4.78)

то величину а следует выбирать минимальной, но при выполнении условий, вытекающих из (4.76) и (4.78),

(1 + 2а)2Д ,а2

< 1 ,

 

3

д---- “—

 

д

 

(1 +2а)2-

 

(4.79)

3

Д

 

 

 

 

(1

+ 2 а ) 2 - 4 « 2

 

 

Если учесть, что в соответствии с постановкой задачи требование по точности получения оценок является ориентировочным, то можно найти значение а, напри­ мер, из уравнения

( 1 + 2 а ) « - а * (г?2 + 4 з + г |з ) /3

(4.80)

д

 

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

мых из уравнения (4.80), а не из системы неравенств (4.79), точность оценок по дисперсии будет отличаться от необходимой незначительно.

Действительно,

относительное отклонение

дисперсии

от допустимого,

например,

при

і — 1 и k = 2 будет

 

(1 + 2а)2-

 

 

(1 + 2а)*-а«(г?2 + г?3+ г|,)/3

_ _ 3

Д

 

~ 3

 

Д

 

_

Е~

 

 

 

(1 + 2 а)2- ^ « 2

~

 

 

 

 

3

Д

 

 

 

 

 

+ г 23

2

 

Г І2 + г 13 + г 23

 

 

 

3

 

 

— Г12

 

3

' 12

 

( л Г

+

2)

г?2

 

4 ~ г *2

 

 

 

 

 

 

 

 

 

 

1 +

Г 9,

 

 

 

 

 

2_

2

 

2_

 

(4.81)

 

 

3

 

4 — ' 12

9 ‘

 

 

 

 

 

 

169

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