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