Следовательно,
£> = М [ ( 6 Л о о ) * ) , |
( 5 . 8 1 ) |
где |
т |
|
|
|
|
|
М.. = |
т г £ |
а< . ' |
<5-82) |
|
9 = 1 |
|
|
Отметим, что в силу |
(5.79) |
для бЯоо можно написать |
другую формулу, которая будет необходима нам в даль
нейшем:
т
8Я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), (—).
Эксперименты типа ( + , —, 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' — число экспериментов.
При экспериментах первой группы для определения р*( + сГо), ц*(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 в экспери |
ментах соответствующего типа. |
|
|
В силу независимости экспериментов различного ти па 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), то безусловно целесообразно использовать первую группу эксперимен
тов. При этом практически всегда получим выигрыш в числе экспериментов, а возможный проигрыш не будет
более чем 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)
Если КтФ®, то при малом а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
Пусть 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 и |
общего числа экспериментов А\:
тт
М -— 9=і г=' |
* |
(5.113) |
— |
D |
|
и а |
т |
|
|
т |
|
|
Ё |
£ Kz 7 |
|
|
|
*<1лг |
(5.114) |
N , = (2т -f-1) |
|
------- |
Перейдем ко второму способу. Обозначим через Nq
число экспериментов для определения каждой из величин цМсТдо), p*g(—ergo) и р*(0), входящих в ад0. Отметим,
что при втором способе значение ц*(0) для получения каждого ctqo находится независимо по соответствующему числу Nq экспериментов.
Общее число экспериментов для получения всех aqo, т. е. для вычисления оценки 6Ха, будет
< 7 = 1
Значения ,n%(üqo), ц*д(—Одо), ц*(0) для определения aqо будут находиться по формулам
|
|
У |
|
Р*Ч (3Чо) = |
Щ |
/=1 |
|
|
|
|
М-%( ~ °?„) |
Na |
S _ q j , |
(5.116) |
Чо' |
/=і |
|
|
|
|
|
|
/=' |
|
В последней формуле индекс (q) у величины |
S {^ ука |
зывает на то, что это значение было получено при опре делении aqо.
Из (5.83) и (5.116) |
находим, |
что |
|
т |
|
Ч о = |
Е |
(5.117) |
|
а=.\ |
|
где |
|
|
|
|
|
- 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, __2т + |
1 9- і |
г=\ |
(5.122) |
|
|
|
|
Если 6 ^ 1 , рационально |
использовать первый |
способ, |
если б>1, то лучше применять второй способ. Перейдем к вопросам практической реализации рас
смотренного в настоящем параграфе метода получения оценок aq0. Для того чтобы выяснить, каким способом лучше определить aq0, необходимо сначала произвести
небольшое число экспериментов, соответствующих пер вому способу. Организация экспериментов при первом способе может быть принята следующей. С одной и той же реализацией воздействий проводятся эксперименты: при номинальных значениях параметров, при отклонении только первого параметра на -f стю, при отклонении толь ко первого параметра на —аю, при отклонении только
второго параметра па + а 2о, при отклонении только вто рого параметра па —а2о и т. д. Затем реализация воз действий изменяется и снова проводятся эксперименты в указанном порядке.
После проведения небольшого числа экспериментов (число экспериментов N при каждом отклонении пара
метров порядка 10—20) вычисляются статистические зна
чения К* . величин |
К 7 _ |
(а, |
r = |
1, |
2, |
т) и затем ана- |
/д/.г |
|
'Ѵѵ |
величина |
|
логично (5.122) |
находится |
|
|
|
|
|
т |
т |
|
|
|
|
|
|
2 |
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 и
при отклонении только второго параметра на —ого. Да лее реализации воздействий снова изменяются и прово дятся эксперименты: при номинальных значениях пара