книги из ГПНТБ / Пугачев В.Н. Комбинированные методы определения вероятностных характеристик
.pdf.лизуемые оценки Лю' и 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) |
|
||||||
у |
J£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 + 2а (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
