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

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

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

 

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

 

 

 

 

 

 

 

 

Kq q =

°

при к ф і;

 

(2.29)

 

 

 

К

 

- —

п2 при k = і.

 

 

 

 

 

 

 

 

 

 

 

 

'

N

Чі

 

 

 

 

 

Учитывая

(2.29) и то, что

 

 

 

 

 

 

 

 

 

 

 

П

 

 

 

 

 

\_Щ_

 

 

 

/(«;<?)= П

 

1

exp

(2.30)

 

 

 

2 Яг

 

 

 

 

 

 

V2nqt

 

 

 

 

из

(2.24)

найдем

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

nj

 

 

 

 

(2.31)

 

К,° - 2 N

Ц

I

 

 

 

 

1

I ^

*)du

 

 

 

S

1!

 

 

 

 

 

 

 

 

Разность

дисперсий

К

и Ко

 

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

значе­

ния X* и оценки Хо в рассматриваемом случае будет

К

- К л

T{Jr

(и) {.(и;

q)du

 

| ^ >

) / > ; q) du

 

 

 

 

 

 

 

 

 

 

 

 

 

Qi

 

1 ^f(«;

q)du]

(2.32)

 

 

 

i- 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем функцию

 

 

 

 

 

 

 

 

 

 

F (и) = ¥ (и) — J ¥ (u') f (u'; q) du'

(2.33)

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

 

 

 

 

 

 

 

 

 

 

 

 

di

 

F

 

~ 1 Л (м’ 9) dn

(2.34)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Яг

 

 

f (и; q) du

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя функцию Л(и), формулу (2.32) перепи­

шем в следующем виде:

 

 

 

 

 

 

 

 

 

к

- к п

N

 

F2(u)f(u; q)du

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

я?

 

 

\

 

 

 

 

-4-21UM'*

 

 

f(w;

q)du

(2.35)

 

 

 

(=1

 

 

 

 

 

 

 

 

 

70

Непосредственной подстановкой коэффициентов ан можно убедиться в справедливости следующего соотно­ шения:

JР2(и)f («; q)'du — ~

УF (и) Яі

)f(u; q)du

 

 

1

 

 

О

 

(2.36)

 

f (и;

q) du.

 

Qi

 

 

При доказательстве данного соотношения необходи-

мо учитывать, что

2

 

 

и,

 

(2.37)

— -

1 / (м; q)du = 2.

 

Яі

Из (2.35) и (2.36) с учетом того, что f(u\ q)>Q, полу­

чаем неравенство (2.27), которое необходимо было до­ казать.

Найдем, при каких условиях в (2.27) имеет место знак равенства, т. е. оценка Ко и статистическое значе­

ние К* одинаковы по точности и выигрыш

в точности

от применения аналитического метода

отсутствует.

Из (2.35) и (2.36) нетрудно получить, что это имеет место, если

і=і

т. е. когда функция W(u) удовлетворяет условию

W (u )= J W(u')f(u'; q)du! > ~ 1}

Функция (2.39) является частным случаем функции сле­ дующего вида:

П

(2.40)

'F (и) = ЬГІ + У] Ь;аГ.

<=і

Нетрудно убедиться, что только при определенных зна­

чениях коэффициентов Ьо, Ьи ... ,

Ьп функция Ч- (и)

будет удовлетворять условию (2.39).

(2.40), вероятност­

В случае, если

Ч(и) имеет вид

ная характеристика

системы К линейно зависит от ком-

7J

понент вектора q вероятностных характеристик воздей­

ствия. Действительно, из (2.18) и (2.40) имеем

П

* = fy> + Е &г<7г.

(2.41)

і= 1

 

Таким образом, оценка Хо и статистическое значе­ ние X* могут быть одинаковыми по точности и выигрыш

в точности от применения аналитического метода к части системы может отсутствовать только в частном случае, когда X является линейной функцией <7, (г= 1, 2, ... , п ) .

Примеры, иллюстрирующие эффективность примене­ ния аналитических методов к части системы, будут рас­ смотрены в § 2.8.

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

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

системы, либо к системе, отличающейся

по структуре

от рассматриваемой, либо к исследуемой

системе, но

при других, не окончательно выбранных значениях ее параметров.

Наличие указанных значений вероятностных харак­ теристик, которые в дальнейшем будем называть пред­ варительными вероятностными характеристиками, обус­ ловлено процессом проектирования системы, в ходе которого исследователи производят выбор и уточнение структуры и параметров системы. Кроме того, происхо­ дит непрерывное уточнение представления о системе. Предварительные вероятностные характеристики могут быть получены или в результате аналитических иссле­ дований или с привлечением метода статистических испытаний.

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

П

При этом под вероятностными характеристиками «упро­ щенной системы» следует понимать предварительные вероятностные характеристики, а под «исходной систе­ мой»— исследуемую. В данном случае для получения оценки помимо аналитических значений предваритель­ ных вероятностных характеристик необходимо найти их статистические значения, для чего потребуются допол­ нительные статистические испытания «упрощенной си­ стемы». Поэтому использование предварительных веро­ ятностных характеристик, как следует из гл. 1, будет целесообразным тогда, когда «упрощенная система» действительно является простой и определение статисти­ ческих вероятностных харктеристик не требует дополни­ тельного большого объема вычислений.

Если же предварительные вероятностные характе­ ристики найдены с привлечением метода статистических испытаний, то использование их для оценки вероятност­ ных характеристик системы может производиться с по­ мощью метода, изложенного в настоящем параграфе.

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

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

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

Пусть вероятностная характеристика Я представляет собой математическое ожидание некоторой случайной величины R. Отметим, что в случае применения метода, изложенного в гл. 1, за R в этой главе следует прини­ мать величину R —/C*BS^C_1SS(5 —ц).

Обозначим через и Х2 значения вероятностных

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

73

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

 

Лг-Л*[Яі] (1 = 1, 2),

(2.42)

где Ri — значения

случайной величины R для г-й точки.

Предположим,

что для определения

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

значения Аі* вероятностной характеристики Хі в первой

точке ранее было

произведено

 

независимых экспери­

ментов. Из этих

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

выделим Nо

первых

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

Значение No

будет

определено

ниже.

По первым Nо

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

и

по

остальным

N і— Nо

могут быть найдены статистические значения Хю* и 1ц*

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

,

-

1

No

 

(2.43)

R

 

 

 

'No

1 3

 

 

;=1

 

 

 

 

 

N,

 

2

* •

 

1

(2.44)

 

7j R,

І^Ѵ.-Л'о

/=ЛГ„+1

где Ra — значение Ri в j-u эксперименте.

В силу независимости экспериментов величины Аю* и Ац* являются независимыми.

Статистическое значение вероятностной характери­ стики в первой точке, найденное по всем Л'і экспери­ ментам, будет

(2-45)

і=1 Величина Хі* является искомым статистическим зна­

чением вероятностной характеристики Аі при перво­ начальных исследованиях, т. е. в первой точке. Для целей же оценки вероятностных характеристик во вто­ рой точке с использованием результатов, полученных в первой точке, будут необходимы величины Ію* и Ап*, которые зависят от числа N0. Оптимальное значение

числа іѴ0, как будет показано ниже, зависит от резуль­ татов статистических испытаний во второй точке, однако оно не превышает N t/2. Отсюда следует, что для при­

влечения при последующих возможных расчетах в дру­ гой точке результатов, полученных в первой точке, не­ обходимо запоминание значений Я ц для js^Ni/2 и Аі*.

74

Действительно, после того, как по экспериментам во второй точке будет найдено число N0, имея значения /Д; для js^'Ni/2, можно по формуле (2.43) найти Яю*, а за­

тем, зная Яі*, вычислить и Яц* по соотношению

Данная формула вытекает из (2.43) — (2.45). Статистические испытания для получения оценки

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

осуществляются при тех же самых реализациях слу­ чайных -воздействий, что и N0 экспериментов в первой точке. В принципе это требует запоминания No реали­

заций случайных воздействий, однако при расчетах на

вычислительных

машинах и образовании воздействий

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

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

чисел {10, 17] можно ограничиться запоминанием только опорных чисел, из которых получались эти последова­ тельности. Остальные эксперименты производятся с не­ зависимой последовательностью случайных воздействий. Обозначая общее число экспериментов во второй точке через N2, можем аналогично (2.43) и (2.44) написать,

что статистические значения Яго* и Ягі* вероятностной характеристики Яг, найденные соответственно по N0

первым и

по остальным

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

будут

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.47)

 

 

 

 

 

/=1

N2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.48)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где R 2j —■значение R2 в ;-м эксперименте.

 

 

 

 

 

 

В

 

силу

независимости

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

 

величины

Яго*

и

 

 

независимы. Так как первые N0 экспериментов

в Я2і*

 

 

 

 

 

 

 

 

 

 

 

 

 

первой и второй точках произведены с одной и той же

последовательностью воздействий, то

Яю*

и'

Яго*

явля­

ются зависимыми. Очевидно, что

Ян*

 

 

 

Яю*,

 

и

 

Аналогично,

 

не

 

не зависит от

Яго*

Я2і*.

Я2і*

 

 

Яю*,

 

Я2о*

 

 

 

 

 

зависит от

 

и Ян*.

Оценку Я2о вероятностной характеристики Яг во вто­

рой

 

точке

будем строить

через

Яю*, Ян*,

Яго*

и

Ягі*

75

в классе линейных несмещенных оценок

оптимальной

в смысле минимума ее дисперсии. Таким образом,

Я2о = яЯго* + b% 2 1* + сЯю*+ dXu*,

(2.49)

где а, Ь, с и с/ — коэффициенты, выбираемые из условия

несмещенности

оценки Я2о и минимума ее

дисперсии.

 

Так как

 

^ І Я І0*] =

М [Я *11] =

Я1,

 

 

 

 

 

 

 

\

 

 

 

 

 

 

ЛГІЯМ*1 =

Л І[Я*„] =

А3

1

 

 

и величины

Яі

и

Яг

могут быть различными,

то условие

несмещенности оценки

 

 

 

 

(2.51)

принимает вид

 

 

 

44 [Яго] = Яг

 

 

 

 

 

 

 

 

 

 

d-\-b —■1, I

 

 

 

(2.52)

 

 

 

 

 

 

 

c + d = 0. I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда и из (2.49)

 

получаем

 

 

 

(2.53)

 

 

Яго= (1&)Яго* + ^Ягі*d{ Яю*— Я ц *).

 

но

Учитывая высказанные выше замечания относитель­

независимости

и зависимости величин

 

(2.48)

и

 

а также соотношения (2.43),

 

(2.47),

 

(2.44), Яю*, Яи*, Яго*

и Я2і*,

 

 

 

 

 

 

 

 

 

 

 

 

(2.53), получим выражение для дисперсии D2 оценки:

 

+ <Г-Г=д ? ----- 2(1 - V d j

f .

(2.54)

 

JV1 JVо

1

о

 

В данной

формуле Ки

и К2 2 — дисперсии

случайных

величин

и /?2, а /Сіг — корреляционный

момент их

в первых N0 экспериментах.

 

 

 

Значения b и d, обеспечивающие минимум диспер­

сии D2, будут определяться из уравнений

 

О

Ь ) ^ - + Ь

к,

 

I Кіг

= 0,

N*— N„

Ул

 

 

(2.55)

 

/fl.

 

6) ~

 

 

 

 

 

d ^ r + d У , - N ' — (1

:0.

 

' ЛГ.

 

76

Отсюда получим

 

Nt

 

 

 

 

Ъ

N t - N t

 

 

 

N>

N,

-r2

 

 

 

 

 

 

N i - N ' N2- N 0

\

(2.56)

 

 

N,

-N,

d =

.-ш/К»_____

 

 

N,

 

 

 

V Ku

Nt

r2 )

 

 

 

N ^ - N , N2- N 0

 

где r — коэффициент корреляции -случайных величии R\ и Ri в первых No экспериментах,

г = * ? !. .

(2.57)

УКиК»

'

Подставляя (2.56) в (2.54), найдем значение диспер­ сии Di при оптимальных коэффициентах 6 и rf:

 

Nt

 

r2

К*г

~ N t- N t

D , N * - N n

Nt

N,

(2.58)

 

N t - N ,

N2- N 0

Будем предполагать, что желаемая точность получе­ ния вероятностных характеристик в первой и второй точках определяется одним и тем же числом М экспери­

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

Ni = M,

(2.59)

а

 

D2 = ^ .

(2.60)

Подставляя (2.59) и (2.60) в (2.58), получим уравне­ ние для нахождения общего числа N2 экспериментов во

второй точке:

К22

М N0

г

М N2N0 м

Т 2

(2.61)

~NT - N T N 2- N 0

Решая это уравнение, найдем, что

Мг — (Мг Nl) г2

(2.62)

77

Число экспериментов N0 выберем таким, чтобы общее число экспериментов N2 для получения оценки вероятно­

стной характеристики во второй точке было минималь­ ным. Нетрудно получить, что No, обеспечивающее мини­ мум N2 , будет

N0

— M 1 + V 1 — т2

(2. 63)

Из этого соотношения видно, что

 

 

N0<M[ 2 = NJ2.

(2.64)

Данное неравенство

использовалось выше.

полу­

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

чения опенки %2 о вероятностной характеристики системы

во второй точке с требуемой точностью, в силу

(2.62) и

(2.63) будет

 

 

N a = M

2 r ';ZlL -

(2.65)

2

1 4 - Y \ — г2

v

Выигрыш г|2о вчисле экспериментов при оценке вероят­

ностной характеристики во второй точке, обусловленный использованием результатов статистических испытаний в первой точке, определим отношением числа М экспе­

риментов, необходимого для получения вероятностной ха­ рактеристики с требуемой точ­ ностью при обычном определе­ нии вероятностной характери­ стики, к числу экспериментов N2 , необходимому при рассма­

триваемом методе:

t)2 o = M /N 2.

(2.66)

Из (2.65) и (2.66) получаем

Рие. 2.2. Зависимость вы­

Ъо

1+ Ѵ Т ^ Т

(2.67)

2Ѵ)

 

игрыша т|2 о в числе экспери­

 

 

 

ментов от абсолютной вели­

Зависимость выигрыша цго

чины коэффициента корре­

ляции г.

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

 

от аб­

 

солютной величины

коэффи­

циента корреляции г представлена на рис. 2.2. Как вид­ но из этого рисунка, выигрыш г\2 о может быть весьма су­

щественным.

78

При оптимальном числе экспериментов М> в соответ­ ствии с (2.53), (2.56), (2.59) и (2.65) формула для оцен­ ки Аго принимает вид

ho* + Ѵ \ - г2 Лгі*

К

1 +

Ѵ \ "^Г 2

 

Kn (1 + К 1 — Г2)2(Я,.* - яп*).

(2.68)

 

Оценка А2о не -может быть практически реализована, так как для ее получения необходимо знать дисперсии К и и Кгг и коэффициент корреляции г. Однако, как бу­

дет показано ів гл. 4, в аналогичных случаях, даже при небольшом числе М, можно Кн, К2 2 и г заменить на ста­ тистические значения Кн*, К2 2 * и г*, найденные по соот­

ветствующим экспериментам при незначительном ухуд­ шении точности оценки. Таким образом, вместо оценки А20 будем рассматривать практически реализуемую оцен­

ку

„ _ _

ho* + V\ - г * * А „ *

 

20 ~

І +

 

 

 

г*2

(Я,„* А.*)-

(2.69)

(1 +

V 1 — А*2)2

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

ции находится по первым .Ѵ0 экспериментам в первой и второй точках, а статистические значения Кн* и К2 2 *

дисперсий— по всем экспериментам для соответствую­ щих точек.

Число экспериментов М0 в силу (2.63) зависит от ко­

эффициента корреляции г. Для практического опреде­ ления числа No вместо г будем пользоваться его стати­

стическим значением г* подобію тому, как это делается в гл. 4, где рассматривается одновременное определение вероятностных характеристик б нескольких точках. Та­ ким образом, вместо числа N0 будем пользоваться чис­ лом N*о, вычисляемым по аналогичной формуле:

N*о =

М -

V 1— г*2

(2.70)

1+

V 1— Г*2“ ’

 

 

 

а общее число экспериментов во второй точке будем определять числом N*2 вместо УѴ2, причем в соответствии

с (2.63) и (2.65)

N*2=2N*O. (2.71)

79

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