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

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

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

U], т. е. от вектора значений внешних воздействий в /-м эксперименте, то Я,- 'коррелирована с Rj и S jt но не зави­ сит от Ru и Sh при І ф к .

С учетом (3.5), (3.10), (3.28) и (3.29) дисперсии и корреляционные моменты, входящие в формулу (3.18) для коэффициента с, будут равны

Ң K s s \

= -дг %рр 1

V

___Ф JS

(3.30)

'» «

N ^ Р Р '

 

К

JL к

 

 

іч*о

N ^ s p ’

 

где Kss, Kpp, Kns, К п р , K S P —■дисперсии и корреляцион­

ные моменты случайных величин R, S

и Р\

Р=-\і + а[Ѳ(07)—#„].

(3.31)

Из (3.18) и (3.31) получаем формулу для коэффици­

ента с:

 

Kps~~ Крр

(3.32)

К . « - 2К е п + К\рр

 

которую можно записать в более компактном виде:

__KRQ

(3.33)

R QQ

 

Здесь Кпд и KQQ — корреляционный момент случайных величин R и Q и дисперсия случайной величины Q, рав­

ной

(3.34)

Q = S Р.

С учетом (3.33) формула для оценки

Ло принимает

вид

 

ав = а * - ^ ( и * - * в).

(3.35)

A QQ

 

Полученная формула для оценки по своей структуре совпадает с формулой для оценки, полученной в первой главе. Здесь, как и в гл. 1, вместо оценки Яо можно рас­ сматривать оценку Яоь которая получается путем замены KnQ и K QQ на их статистические значения K * PQ и K * QQ-

И0

Таким образом, практически реализуемая оценка Яоі бу­ дет

я„, = я*

V

ѵъУ

(3.36)

(и-* -

Причем, как указывалось в гл. 1, точность оценок Яо и Ящ уже при сравнительно небольшом числе эксперимен­ тов примерно одинаковая, что позволяет анализировать точность оценки Яоі по величине дисперсии оценки Яо.

Для К*RQ и K*QQ можно использовать формулы, вы­

текающие из (3.39) и (3.34):

1 Л

1

^— **«) ~ а Іѳ(у з)

/= 1

(3.37)

N

 

^QQ'- N 5 ] { ( S j - f t ,) - « [8 М ) -?*]}*•

i = 1

В формулы (3.37) входит матрица-строка а. Компо­

ненты ее могут быть получены следующим путем. Так как при достаточно большом числе экспериментов q* близко к q, то

__ Ф

_

ф 0

(3.38)

dqt

~~

d q \ ‘

 

Значение же производной öpo/dqt* можно приближенно определить как отношение изменения Лио вероятностной характеристики модели щ за счет изменения A q f і-й ве­ роятностной характеристики q f к величине этого изме­

нения.

Таким образом, для определения и необходимо про­ вести дополнительные исследования модели системы ли­ бо аналитическим методом, либо методом статистических испытаний. Причем, если применяется метод статистиче­ ских испытаний, то для определения Дцо/Дqi* значение вероятностной характеристики модели для исходного вектора вероятностных характеристик и измененного на Л<7і* значения его і-й компоненты целесообразно полу­

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

Описанный выше способ получения а с учетом фор­ мул (3.36) и (3.37) позволяет вычислить оценку Яоі.

I l l

3.3. Точность оценки вероятностной характеристики. Выигрыш в точности оценки или в числе натурных экспериментов за счет использования результатов теоретических исследований

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

шом числе экспериментов можно считать с достаточной степенью точности равной дисперсии оценки ?„оі-

Дисперсию Л'о найдем, подставляя с из (3.18) в (3.17):

(3.39)

Учитывая (3.30) и то, что

 

 

К ц

N h RR,

(3.40)

получаем

 

 

 

 

К п

К RR

1 -

(KRS ~ KRP)2

(3.41)

 

N

 

(^SS 2Ksp + Kpp)

 

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

Выигрыш в точности оценки К о ( К о і ) от использования

результатов

теоретических исследований будем опреде­

лять отношением дисперсии статистическаго значения К *

к дисперсии

оценки ;Хо:

т]0= Кп /Ко.

(3.42)

В силу (3.40) и (3.41)

 

(3.43)

( K R S K R P ) 1_________

КRR (Кss 2^SP+ Kpp)

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

112

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

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

Очевидно, что наибольший выигрыш в точности или

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

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

Предельное значение rjoo выигрыша rjo будем опреде­ лять в предположении, что модель и система полностью идентичны и ошибки измерительной системы нулевые. Значения R, S и Р при этих условиях обозначим через Ro, S a и Ра. Очевидно, что

Ra = S 0.

(3.44)

Отсюда и и з '(3.43) получаем предельное значение вы­ игрыша:

L- ‘ .,KJ S A (Rs.

2КS„PÖ

(3.45)

й

RJoJoS S КгPntP n— Кs»p»

 

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

ности системы

 

т

 

(3.46)

V

І=1

где Ui — значения воздействий.

Значения

воздействий U-t (/ = 1, 2, ..., т) будем счи­

тать независимыми между собой, подчиняющимися нор-

8 — 288

113

мальному закону распределения с нулевыми математиче­ скими ожиданиями и средними квадратическими откло­ нениями Оі ( /=1 ,2, ... , т).

Так как искомой вероятностной характеристикой си­

стемы является дисперсия

,и = 7)у

значения Y, то

(3-47)

Очевидно, что

s . =

=\i ( S

^

V

 

 

 

=1

 

 

 

 

=

D„

\1

2 2

(3.48)

 

2jao..

 

 

 

 

ém*

г ,

 

Отсюда следует,

что

 

 

 

 

 

 

 

 

т

 

 

 

 

Ро =

! >

’ (#* .

(3.49)

 

 

1=1

 

 

 

где (з2)* — статистическое

значение дисперсии

а2.

Из (3.49) нетрудно получить выражение для /V

 

 

 

т

 

 

 

 

Л =

2

а Ѵ \

(3.50)

г= 1

Учитывая, что Ui независимы и подчиняются нормально­

му закону распределения, из (3.45), (3.47) и (3.50) полу­ чаем

\ о

(3.51)

Если величина едет* одинакова для всех г, то

роо — t r i ,

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

Рассмотрим другой случай. Пусть Ui представляют

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

ряют условию

аі = Ае ІІС,

где А — некоторая постоянная величина; с — память си­

стемы дискретного действия (по числу значений воздей­ ствий) ; тогда

1- е - 4'с

Ъо ~ (1 - е - 2'е)*'

На рис. 3.1 представлена зависимость предельного выигрыша от памяти с. При 'больших с

цоо^с.

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

Рис. 3.1. Зависимость предельного выигры­ ша г|оо от памяти с.

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

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

Для рассмотренного выше частного случая выигрыш г)о при наличии ошибок измерений значений процессов

8*

115

и воздействий будет

Ъ = ------ --------------

 

т

1

\2

 

 

т

 

 

(/=1 А=1

/

 

{ т

\

Mfe

т

/

1

pm

 

S «М+. 2

/

ÜLi

 

ai al (3?+°д;) (aÄ+ °üfe)

 

Ѵ=1

i~1

 

 

 

 

i^=k

 

 

 

 

 

 

 

(3.52)

где ад и ад. — средние

 

квадратические отклонения оши­

бок измерений значений У и

процесса и воздействия.

Пусть выполняются соотношения

 

 

 

(»=1. 2,..., m),

(3.53)

 

 

 

где у — некоторая величина.

 

 

 

Неравенства (3.53)

 

означают, что отношения средних

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

Из (3.51) — (3.53) нетрудно получить, что при этих

условиях

 

1

(3.54)

 

(1 +Y2)4

 

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

ворять неравенству

 

 

3 — ^оо — ^

4У!гІо°

(3.55)

Чоо

1+ 4Y2^OO

 

Так, например, при -у=0,1 и г]о»=5 относительное умень­ шение выигрыша не будет превышать 0,08.

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

116

3.4. Оценка вероятностных характеристик системы по натурным испытаниям и теоретическим

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

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

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

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

117

и расчетные значения вероятностной характеристики мо­ дели системы. Значения вероятностных характеристик системы и модели, полученные по окончательным натур­ ным экспериментам, будем обозначать соответственно через %*, р* и ро, а значения вероятностных характери­

стик модели и системы, найденные по предварительным натурным экспериментам, — через р*п и ропФормулы для X* и р*, а также способ получения ро были приве­

дены в § 3.1 и 3.2. Формула для рп* аналогична форму­ ле для р*. Способ получения р0п в принципе такой же, как и способ получения ро, однако, как будет показано ниже, зная ро, можно роп вычислить значительно проще.

Пусть число окончательных и предварительных экс­ периментов равно соответственно N, и Ап. Тогда в соот­

ветствии с § 3.2 можно написать следующие выражения для X*, р* и рп*:

2* N

/= і

(3.56)

і=і

где Sin—’Значение S в і-м эксперименте с моделью си­ стемы, проведенном по результатам і-го предваритель­ ного натурного эксперимента.

Статистическое значение q* вектора q вероятностных

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

периментов, в соответствии с § 3.2 будет

N

(3-57)

/=1

Аналогично, статистическое значение qn*, найденное

по Ап предварительным натурным экспериментам, будет

Nu

q*a=^YiHUia)’ (3-58)

і=і

118

где Ui n — значение вектора Ü в і-м предварительном на­

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

Расчетное значение ро вероятностной характеристики модели системы находится по способу, описанному в § 3.1. После того как найдено р0, расчетное значение роп по предварительным натурным экспериментам мо­ жет быть получено более простым путем. Действительно, для ро справедливо соотношение (3.2G). Аналогично для роп можно записать, что

роп= р Т сі(Ѵп*_ <7и)-

(3.59)

Отсюда и из (3.26) получаем формулу, по которой удобно вычислить р0іІ:

роп = ро+ а ( < J n * Q * ) ■

(3.60)

Напомним, что необходимое для вычисления роп по данной формуле значение а находится для других целей

(см. § 3.2). Таким образом, если ро определено, то вы­ числить рок по (3.60) не представляет особого труда.

Оценку Яо вероятностной характеристики системы най­ дем через Я*, р*, ро, рп* и роп в классе несмещенных линейных оценок, оптимальных в смысле минимума дис­ персии.

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

 

Яо—оХ* + b\i* + cpo+ dpn* + (fpon,

(3.61)

где а, Ь, с, d,

е — коэффициенты, выбираемые, исходя из

несмещенности оценки Яо и минимума ее дисперсии.

Нетрудно

получить условие несмещенности:

 

 

а=1,

(3.62)

 

b + c + d + e = 0.

(3.63)

Предварительные и окончательные натурные экспери­ менты будем считать независимыми, поэтому Я*, р* и

ро не зависят от рп* и роп.

Ко оценки Яо бу­

С учетом условия (3.62) дисперсия

дет

 

 

 

к . = К „ + 2 І Ж , < + 2 с К ^ + Ь ‘ К „

+ 2

Ь с К +

 

+

f

.

(3.64)

r n r on

r o n r on

 

119

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