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

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

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

Следовательно,

£> = М [ ( 6 Л о о ) * ) ,

( 5 . 8 1 )

где

т

 

 

 

 

 

М.. =

т г £

а< . '

<5-82)

 

9 = 1

 

 

Отметим, что в силу

(5.79)

для бЯоо можно написать

другую формулу, которая будет необходима нам в даль­

нейшем:

т

0 0

£

[н-*я («go) -

(0)].+ ft**g ( - °9о) -

Н* (0)1

(5.83)

 

2

 

 

п^=1

 

 

 

 

В настоящем параграфе будем рассматривать такую организацию статистических испытаний упрощенной си­ стемы для определения оценки бло, при которой будет обеспечена заданная величина D0 дисперсии D значения

6Яоо этой оценки.

Исследуем сначала случай т = 1, т. е. когда опреде­

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

этого

случая

 

 

 

 

 

8Я = ~ а з 2,

 

 

 

1

2

(5.84)

 

 

оЯо= - Y

а0<з\

 

а 0 =

~ 2

(3о) — Iх* (0)] 4 -

К ( — О0) — Iх* (°)]}-

 

а0

 

 

 

8Я00 = 4 " (I»1* Ю — Iх* (°)] +

[Р* ( — ао) — Iх* (0)1}.

Возможные

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

определения

бЯо (ао) будем

обозначать через « + », «0» и «—», указы­

вая тем самым, что эксперименты проводятся при откло­ нении параметра от номинального значения соответст­ венно на +оо, 0 и —ао- Все возможные типы экспери­ ментов разобьем на три группы:

I группа — эксперименты типа ( + , —, 0);

II группа — эксперименты типа ( + , 0), (—, 0),

( + . —);

III группа — эксперименты типа ( + ), (0), (—).

230

Эксперименты типа ( + , —, 0) проводятся таким об­ разом, что при отклонении параметров на -ТсГо, —0о_и 0

используется одна и та же последовательность воздейст­

вий на

систему. Аналогично,

эксперименты типа ( + ,

0) — с

использованием одной

последовательности воз­

действий и т. д. Эксперименты различных типов прово­ дятся с независимыми последовательностями воздейст­

вий.

 

 

 

 

 

 

 

 

 

 

 

рациональна

 

Покажем, что для получения оценки

 

 

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

первой группы. С этой

 

 

біо

 

 

 

целью определим дисперсии D0 оценки

біо

для указанных

групп экспериментов.

 

 

 

 

 

 

 

 

Вероятностная характеристика ц является математи­

ческим ожиданием случайной величины S. Поэтому

 

Aßo

 

 

 

 

n (A ß o )= M {S (A 'ß o )],

 

 

 

 

(5.86)

где

— неслучайное отклонение

параметра от

номи­

 

 

 

 

 

 

— случайная

величина,

зави­

нального значения;

 

 

 

сящая от

неслучайного отклонения

 

 

 

 

 

 

 

 

 

S ( A ß o )

 

 

Aßo.

 

 

 

 

 

 

Разложим S

(Aßo)

в ряд по

 

Aßo,

ограничив­

шись в

 

 

 

 

степеням

 

 

минимальным

 

 

 

даю­

 

 

разложении

 

числом членов,

щим отличное от нуля значение 51. Таким разложением будет разложение до второго порядка включительно, т. е.

S (Aßo) = Со + C i Aßo Т C 2(Aßo)2,

(5.87)

где С0, Ci, С2— случайные величины, зависящие от слу­

чайных воздействий на упрощенную систему.

 

и

В

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

(5.86)

и

(5.87)

для

р * (

+

ао), р *(0 )

 

ц*(—его) получим следующие выражения:

 

 

 

 

Р* ( +

 

 

Л"

 

 

С , P o +

С г Р о )’

 

 

 

 

3о) — дгГ /=1 (C 0J +

 

(5.88)

 

 

 

( ° ) =

N'

 

 

 

 

 

 

 

 

 

 

^ /=1

с -

 

 

 

 

 

 

 

 

 

 

 

£

 

 

 

 

 

 

 

 

 

р* ( -

*с) =

 

N’

<С »3 ~

С >Яо +

 

 

 

 

 

 

W

Л

С , р ,

 

 

 

 

 

 

 

 

 

 

 

 

 

/'=1

Здесь C0j, Cij, C2j — значения С0, Си С2 в /-м эксперимен­ те, а N' — число экспериментов.

231

При экспериментах первой группы для определения р*( + сГо), ц*(0) и р*(—(То) используется одна и та же последовательность воздействий, т. е. значения C0j, Cij и C2j во всех формулах (5.88) одинаковые. Обозначая число экспериментов N' для первой группы в каждой

точке (при каждом отклонении парамера) через АД_о, из (5.85) и (5.88) получим

 

9

Л У г - о

Csj-

(5-89)

 

?Я°° = Л ^ 7 7

S

 

/=•

 

 

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

будет

ІѴ1 = ЗУѴ+_0.

 

(5.90)

 

 

Дисперсия Di оценки бЯоо для первой группы экспе­

риментов будет

 

 

 

 

 

4

 

К 22,

 

А

N Ч— о К,

 

(5.91)

где Кгг — дисперсия коэффицента

С2.

из со­

Полагая А равным требуемому значению Dv,

отношения (5.91) найдем необходимое значение общего

числа экспериментов для первой группы:

 

^ = Д Г 0о-

(5-92)

-L/0

 

Можно показать, что из экспериментов второй груп­ пы имеет смысл рассматривать только эксперименты ти­ па ( + , 0) и (—, 0), исключив эксперименты типа ( + , —) .

Очевидность

данного

положения

вытекает из

того, что

в выражение

(5.85)

для 6Дю входят разности ц*( + его)—

—р*(0) и р*(—0о)—[х*(0).

 

 

Пусть N+ о и АОо — числа экспериментов типа '( + , 0)

и (—, 0). В соответствии е (5.88)

 

 

 

 

 

N+o

 

 

р* (+ ао) - ѵ>*( ° ) = ж д - S

(с '.яо + C'2j4

),

 

 

 

 

 

(5.93)

Р*(— °e) - I 1* (°) =

цг

J ] ( - c " , j3o + c ' v 02 ),

j

 

 

 

 

 

где C'ij, C ' 2j

и C" ij ,

C " 2j — значения С4 и С 2 в экспери­

ментах соответствующего типа.

 

 

232

В силу независимости экспериментов различного ти­ па C'ij, C'zj не зависят от С"і„ С"г;-

Из (5.85)' и (5.93) получим значение дисперсии Dz

оценки бЯоо для второй группы экспериментов:

Da =

- j - 17777(^пао +

2^ і23о+

Ä 22ао) +

 

+

" 2

K + K

} }’

(5-94)

где Ки, Kzz и /(i2— дисперсии и корреляционный момент случайных величин Сі и С2.

Общее число экспериментов для второй группы равно

N2 = 2(N+0 + N - O) ■

(5.95)

Полагая дисперсию Dz равной требуемому значению Do и выбирая числа N+0 и УѴ_0 оптимальными в смысле

минимума общего числа экспериментов, получаем

=+ 2/б12а0 -(- К22?1 +

 

+

Кп - 2 К 12о0 +

К2/ 0 У*1 .

(5.96)

Покажем, что, за редким исключением, Л/2>>ЛѴ Из

(5.92)

и (5.96)

имеем

2

 

 

 

АД

 

 

 

 

 

Лг23о

 

 

 

лд

^ l l

Н"" 2 iC j2 3 0 + ^ 2 2 30 T "

^ l l ■—

2/V ,2 3 O - f - К г 2 3 0 ^

 

Если /Сц=£0, то при малом а0

 

(5.97)

 

 

 

 

^

3

/С22

2

 

(5.98)

 

 

ЛД ~

2

/Сп

0

 

 

 

 

 

Отсюда видно,

что А^і/ТѴг

представляет

собой

малую ве­

личину порядка о2о.

В крайнем случае, когда Л'ц= 0, корреляционный мо­ мент тоже равен нулю и

JVI/JV2=1,5.

Таким образом, как правило, АД<СЛД и только в ис­ ключительном случае N I > N 2, но не более чем в полтора

раза. Следовательно, если не проверять при проведении экспериментов величину отношения (5.97), то безусловно целесообразно использовать первую группу эксперимен­

233

тов. При этом практически всегда получим выигрыш в числе экспериментов, а возможный проигрыш не будет

более чем 1,5.

 

 

третьей

группы. Пусть

Перейдем к экспериментам

N+, No и АП — числа экспериментов типа

( + ), (0) и (—)•

В силу (5.88)

м+

 

 

 

 

 

 

 

Н* (+ 3,) =

^ 5 ]

(Соз + С'іНо + С'іР о)>

 

/=Л

 

 

 

І‘* ( 0 ) = д г £ с " оі,

 

} (5.99)

 

/=1

 

 

 

 

 

 

 

N

 

 

 

Ѵ-*( - 3») =

лС

(С'"оі -

С " \ р 0+

С ' \ р І ),

 

/=1

 

 

 

где C'ij, С"и, С"'іj — независимые между собой значения Со, Сі и С2 в экспериментах соответствующего типа.

Из (5.85) и (5.99) получим значение дисперсии £>з

оценки для третьей группы экспериментов:

 

 

=

j щ (Коа+

Кц320 Ң- Кгр\ + 2 К0р 0+

2К0р~ +

+

2 / ( 123 Q

) + д

Р ( Л о о

+

+

^ 2 2 ° 0

~

 

2 ^ 0 1 3 0 +

2 ^ 0 2 3 0

2 ^ . 2 3 0 )

і+

д ^ - ^ О О »

(5.100)

где Кдг — корреляционный

момент

случайных

величин

Сч и Су.

 

 

 

 

 

 

 

 

 

Общее число экспериментов для третьей группы рав-

но

 

Ns — N+ + N - + N0.

 

 

 

(5.101)

 

 

 

 

 

Полагая Dz равной требуемому значению D0 и выби­ рая числа АП, А/_ и N0 оптимальными в смысле миниму­

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

N

3

1

1

1/ Х

~ 4-

 

 

D0

\ ѵ

4 о» 4 -

К оо + ^ 1 1 3 0 +

^ 2 2 3 0 +

 

2 ^ 0 1 3 0 + 2 ^ 0 2 3 0 Н ~ 2 ^ і а 3 0

+4П ^. + '<.,’о4 '<4 - 2 К , л + 2 К , / 0 - 2К „ 4 }.

(5.102)

234

Если КтФ®, то при малом а0

(5 .103)

Ne~ Km/Do

и

(5.104)

N 1/N3 3/С2204о//Соо-

Следовательно, при КооФО отношение Ni/N3 представ­

ляет собой малую величину порядка о40, т. е. JV3>7VI.

 

В крайнем случае, когда /Соо=0, то /<01= 0

и Кю = 0.

Если при этом КнфО, то из

(5.92) и (5.102) при малом

сто получим

 

 

Д9__ Жы 2

 

 

 

 

 

 

 

 

 

 

(5.105)

 

 

 

 

 

N,

~ Ки 3° ’

 

 

 

 

 

 

 

 

 

 

т. е. отношение N\IN3 является малой величиной, порядка

(Т2о, и, следовательно, Nз^>'АЦ

 

 

когда

 

= 0

и

И

 

только

в исключительном случае,

/Соо

 

(при

этом

К

оі

= К 2=К

і

2=0), ів

силу

 

и

/Си= 0

 

 

 

 

о

 

 

 

(5.92)

(5.102)

 

 

 

Ni/Nз = 3;

 

(5 .106)

 

т. е.

N i> N 3.

 

 

 

 

 

 

если

не

проверять,

имеют

ли место

Таким образом,

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

Перейдем к общему случаю. Будем предполагать, что для определения aqQпо формуле (5.79) значения \L*q(oqo),

ц*(0) и p%(—<т9о) находятся с использованием одной и той же последовательности воздействий. Целесообраз­ ность такого предположения была показана выше для случая т= 1.

Рассмотрим два способа определения ачо. Первый способ, когда все aqо определяются с использованием

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

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

Рассмотрим первый способ. Обозначим через S +q, S-q и So случайные величины, математические ожидания

которых равны рЦс^о), цд(—адо) и ц(0), т. е.

Ѵ-Ч(57о) =

М [“5+9ІІ 1

=

(5Л07)

**(0) = M [ S 0]. I

235

Пусть N — число экспериментов

для

определения

каждой из величин

ц.*</(а,іо)

и p*f/(—

ст^о)

(<7=1, 2, . . ,,т)

и р*(0). Тогда

 

 

 

 

 

 

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

ния всех aqо, т. е. оценки біо,

будет

 

 

 

 

 

N 1=(2m + l)N.

 

 

(5.108)

Значения р.% (ад0), р*д (—а до ), р*(0)

 

определяются

по формулам

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

I1*? (

39о)---

 

 

 

 

(5.109)

 

 

 

 

/= 1

 

 

 

 

 

“*<0> = T

£ S.J,

 

 

 

 

 

 

 

1= 1

)I

 

 

где S+gj, S —qj и S 0j — значения S+g, S_g и So в /-м экспе­

рименте.

Для каждого / все случайные величины S+gj, S - gj и

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

Из (5.83) и (5.109) получим, что

т

( , JL

1

 

ч ,о = у;<| - г Ё

Z,j} ’

(5.110)

(7=1

 

 

 

Zqj = ~2~ [(5+93

"~ S oj) +

(S_4j~ S oj)l

(5.111)

Дисперсия біоо при первом способе будет

тт

(5112)

 

 

<7 = 1

r = 1

 

 

где

K z zг — корреляционный

момент

случайных величин

Z g j

И Zrj.

 

 

 

 

Полагая

равным требуемоему

значению D0,

най­

дем

необходимые значения

чисел экспериментов

N и

236

общего числа экспериментов А\:

тт

М -— 9=і г='

*

(5.113)

D

 

и а

т

 

 

т

 

 

Ё

£ Kz 7

 

 

 

*<1лг

(5.114)

N , = (2т -f-1)

 

-------

Перейдем ко второму способу. Обозначим через Nq

число экспериментов для определения каждой из величин цМсТдо), p*g(—ergo) и р*(0), входящих в ад0. Отметим,

что при втором способе значение ц*(0) для получения каждого ctqo находится независимо по соответствующему числу Nq экспериментов.

Общее число экспериментов для получения всех aqo, т. е. для вычисления оценки 6Ха, будет

т

 

ЛЛ = 3 £ Л Ѵ

(5.115)

< 7 = 1

Значения ,n%(üqo), ц*д(—Одо), ц*(0) для определения aqо будут находиться по формулам

 

 

У

 

Р*Ч (3Чо) =

Щ

/=1

 

 

 

 

М-%( ~ °?„)

Na

S _ q j ,

(5.116)

Чо'

/=і

 

 

 

 

 

 

/='

 

В последней формуле индекс (q) у величины

S {^ ука­

зывает на то, что это значение было получено при опре­ делении aqо.

Из (5.83) и (5.116)

находим,

что

 

т

 

Ч о =

Е

(5.117)

 

а=.\

 

237

где

 

 

 

 

 

- 4 -

KS J - s i? )+

(s - «

-

S!?)J-

(5-118)

Все величины

Zq j (q= \, 2, . .

m

и

/= 1 , 2,

. . Nq)

являются между собой независимыми. Поэтому диспер­ сия 6Лоо при втором способе равна

т

D-=S^v,- (5Л19)

9=1

Числа Nq выбираем оптимальными в смысле минимума общего числа экспериментов N2 при условии, что дис­ персия D2 равна заданной Do-

Нетрудно получить, что

Na

V

l v

N..

(5.120)

 

 

 

Г=1

 

 

 

N = 3

EYKh‘

(5.121)

 

Dn

 

 

 

 

 

Отношение чисел экспериментов при первом и втором

способах будет

 

 

 

 

 

 

т

т

 

 

 

Е Е NZq?r

 

N, __+

1 9- і

г=\

(5.122)

 

 

 

 

Если 6 ^ 1 , рационально

использовать первый

способ,

если б>1, то лучше применять второй способ. Перейдем к вопросам практической реализации рас­

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

небольшое число экспериментов, соответствующих пер­ вому способу. Организация экспериментов при первом способе может быть принята следующей. С одной и той же реализацией воздействий проводятся эксперименты: при номинальных значениях параметров, при отклонении только первого параметра на -f стю, при отклонении толь­ ко первого параметра на —аю, при отклонении только

238

второго параметра па + а 2о, при отклонении только вто­ рого параметра па —а2о и т. д. Затем реализация воз­ действий изменяется и снова проводятся эксперименты в указанном порядке.

После проведения небольшого числа экспериментов (число экспериментов N при каждом отклонении пара­

метров порядка 10—20) вычисляются статистические зна­

чения К* . величин

К 7 _

(а,

r =

1,

2,

т) и затем ана-

/д/.г

 

'Ѵѵ

величина

 

логично (5.122)

находится

 

 

 

 

 

т

т

 

 

 

 

 

 

2

2

к*.

 

 

_

+ 1

<7=1

г—Л

(5.123)

с

— 3 —

 

 

 

 

2 у к*

Ѵ<7=1

Если окажется, что ö*=^l, то определение ачо необ­

ходимо осуществлять первым способом, т. е. продолжать ранее начатые эксперименты. Эксперименты следует за­ кончить как только фактическое число их N станет удов­

летворять условию

2 2 к*

N

q — \ г — 1

 

(5.124)

 

о»

 

 

 

 

 

В этом неравенстве

 

K*z г

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

значе­

ние Kz z , найденное по

всем

проведенным на

данный

момент времени экспериментам.

Если 6*>1, то для определения ас,о рационально при­

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

при

отклонении только первого параметра на +0ю и

при

отклонении только первого параметра на —стю-

Затем реализации воздействий изменяются и проводятся эксперименты: при номинальных значениях параметров, при отклонении только второго параметра на + 020 и

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

239

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