книги из ГПНТБ / Пугачев В.Н. Комбинированные методы определения вероятностных характеристик
.pdfПример 2. Пусть в исходной системе осуществляется усреднение с весовой функцией р{т)>0. Тогда импульсная переходная функция
системы будет
t 1
W (t , т) = р (т) J P M dll
.о
В качестве упрощенней примем систему, которая производит про стейшее усреднение. Импульсная переходная функция такой упрощен ной системы равна
Wy (t, т) = — .
Воздействия на исходную и упрощенную системы — нестационар ный нормальный белый шум с корреляционной функцией
K ( t u /s)=C(<l)c(/2)Ö(<l-M.
где 6 (ti—t2)—ö-функция; c( t ) — функция, определяющая интенсив ность белого шума (c(t)>0).
Искомые вероятностные характеристики исходной и упрощенной систем — дисперсии значений У и Z соответствующих процессов. В данном случае коэффициент корреляции rYz будет определяться
по соотношению |
t |
|
|
- |
|
2 |
|
|
^ с2 (т) р (т) d-c |
|
|
.0 |
“ г t |
' |
|
Г t |
dv |
||
f с2 (т) |
Г С2 (т) |
(т) Дс |
|
.0 |
|
. _0 |
|
а выигрыш г)о — по формуле
1
Для оценки значения выигрыша воспользуемся неравенством, приведенным в [5]:
пп
0 < / я 1< х і < Л І 1, 0 < т2 Рі A!s.
В соответствии с этим неравенством получаем, что
r h > 4(аІ/2+ а - > /у г,
50
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
* |
|
^ ___ Стах |
(cp) max |
|
|
|
||
|
|
|
|
|
|
|
Стгп |
(cp)mi n |
|
|
|
|
|
|
|
|
^тах = max с (x) |
|
cmin = |
min c (x), |
|
||||
|
|
|
|
|
|
X |
|
|
|
|
X |
|
|
(cp)max = |
max c (x) p (x); |
(cp)min = |
min c (x) p (x). |
|
|||||||
|
|
|
|
|
X |
|
|
|
|
X |
|
|
Отсюда для выигрыша rjo получаем следующее неравенство: |
||||||||||||
|
|
|
^Зо ^ |
*0отіп — |
1 — |
16 |
( а 1/ 2 + |
а . - ' 12) ' * |
|
|||
|
|
|
|
|
|
|
|
|||||
Зависимость |
отіп от а представлена на рис. 1.8. |
|
||||||||||
Для того чтобы яснее представить возможные значения выигры |
||||||||||||
ша тіо, определим величины а для |
следующих |
конкретных |
случаев: |
|||||||||
а) |
с(т)=Со, |
т. е. воздействие— |
?отіп |
’ |
|
|
||||||
стационарный |
белый |
шум. |
Тогда |
|
|
|||||||
|
|
Ct —Pmaxlpmin, |
|
|
|
|
|
|
||||
|
где |
ртах = |
max р (х) |
|
|
|
|
|
|
|||
|
|
|
|
|
X |
|
|
|
|
|
|
|
|
и Pmin = |
т І П р ( х ) . |
|
|
|
|
|
|
||||
|
|
|
|
|
t |
|
|
|
|
|
|
|
Как |
видно |
из рис. 1.8, чем |
|
|
|
|
||||||
МеНЬШе Ct—рmaxIРтп іп* тем больше |
|
|
|
|
||||||||
выигрыш Г)Оmin- |
Однако даже при |
|
|
|
|
|||||||
а = 3 |
выигрыш |
|
T)o^Tiomtn = 2,3. |
|
|
|
|
|||||
б) |
Усреднение |
является |
опти |
|
|
|
|
|||||
мальным. В данном случае р(т) = |
|
|
|
|
||||||||
= 1/с2(т); |
|
|
|
|
|
|
|
|
|
|
|
|
|
(ср) Viах—1/сmin» |
|
|
|
|
|
|
|||||
|
(CP) min ” |
1ІСтах\ |
|
|
Рис. 1.8. |
Зависимость |
мини |
|||||
|
tt“ (Стая/Cmift)2* |
|
|
|||||||||
|
|
|
мального |
ВЫИГрЫШа |
TJomifi |
|||||||
Пример 3. |
Пусть |
вероятност |
от а. |
|
|
|
||||||
ные характеристики |
исходной |
и |
|
|
|
|
||||||
упрощенной систем представляют собой математические ожидания
соответственно р-й и q-й степени |
случайной |
величины X (рФя), |
т. е. |
|
S = X«. |
£ = Хр ; |
!X= M[S]; |
Обозначая через Мь, момент k-n степени случайной величины X, можем написать, что
2 |
{Mp+q — MpMq)2 |
> s = |
(M2P- M 2p )(M,q - M q)2 |
Выигрыш 7j0 определяется через r^s по формуле (1.37). Значение
выигрыша т]о при q= 2 и р=4 и 6 и различных законов |
распределе |
ния случайной величины А' приведены в табл. 1.1. |
|
4* |
51 |
|
|
|
|
|
|
Т а б л и ц а 1.1 |
|
|
|
|
|
|
Закон распределения |
||
|
Значения |
|
|
|
Равновероят |
|
|
|
Р И |
|
Нормальный |
Треугольный |
|||
|
|
|
ный |
||||
|
|
|
|
|
|
|
|
р — 4 , |
<7=2 |
|
|
4 |
|
1 2 ,2 |
іб |
р — 6 , |
<7=2 |
|
|
1 ,6 |
|
5 ,1 |
6 , 2 |
Рассматривались |
следующие |
законы |
распределения: |
|
|||
нормальный |
|
|
|
|
|
|
|
|
|
fx (Х) - ѴТп |
2 |
о> . |
|
||
|
|
|
|
|
|||
ра вновероятный |
|
J _ |
|
|
|
|
|
|
|
|
\х\ < а , |
|
|
||
|
fx ( * ) = |
2а |
|
|
|||
|
|
|
|
|
|||
|
О |
х\ > а; |
|
|
|||
|
|
|
|
|
|||
треугольный |
|
0 |
X < |
0, |
|
|
|
|
|
|
|
|
|||
|
|
|
2х |
* 0 < |
X < |
а |
|
|
|
|
а |
|
|||
О, X > а.
Из данной таблицы следует, что даже при значительном разли
чии степени нелинейности исходной |
( р= 4 и 6) и упрощенной (q=2) |
|
систем выигрыш тіо является весьма существенным. |
характеристики |
|
Пример 4. Рассматриваемые |
вероятностные |
|
исходной и упрощенной систем при применении метода коррелирован ных процессов могут различаться по физическому смыслу. Представ ляет определенный интерес оценить выигрыш для подобных случаев. Отметим, что пример 3 относится именно к такому случаю, гак как в нем вероятностные характеристики равны математическим ожида ниям различных степеней случайной величины.
В настоящем примере вероятностная характеристика исходной системы равна вероятности р того, что случайная величина X будет находиться в определенных пределах, а для упрощенной системы —
дисперсии о2 этой величины, т. е. |
|
|
|
||
\ = |
р = М [/?]; |
/1 |
при |
I 2 Г |< а , |
|
\о |
при |
I X I |
|||
|
|
||||
|
р. = ч2 = Д4 [SJ; S = X |
2. |
|||
Для вычисления выигрыша предположим, что X подчиняется нормальному закону распределения и имеет нулевое математическое
52
ожидание. Нетрудно показать, что тогда
|
|
|
|
|
2 |
|
2 (а/и)* ?*(«/») |
|
|
(1.140) |
||||
|
|
|
|
0 ? s - |
Ф (й/а) [1 — Ф (а/а)] |
’ |
||||||||
|
|
|
|
|
||||||||||
где Ф(а/<т) — интеграл вероятностей, а |
|
|
|
|
||||||||||
|
|
|
|
?(а/о) = |
- р |
= - е_1/2(а/0)а |
• |
|
(1.141) |
|||||
|
Так как |
|
|
|
|
р = |
Ф (а/а), |
|
|
|
(1.142) |
|||
|
|
|
|
|
|
|
|
|
|
|||||
то |
2 |
|
|
|
|
|
можно в неявном |
виде |
выразить через р. |
|||||
rRS И| следовательно, т)0 |
||||||||||||||
|
Зависимость |
выигрыша |
т]0 от вероятности р, полученная из |
|||||||||||
(1.37), (1.140)-т-(1.142), приведена на рис. 1.9. |
Как следует из ри |
|||||||||||||
сунка выипрыш г)о является суще |
|
|
|
|
|
|||||||||
ственным для достаточно широко |
|
|
|
|
|
|||||||||
го |
диапазона |
значений |
вероятно |
|
|
|
|
|
||||||
сти р. |
|
|
|
|
|
|
|
|
|
|
|
|
||
Г. Р асчетны й |
п р и м е р |
|
|
|
|
|
|
|
|
|
||||
ра |
В качестве |
расчетного |
приме |
|
|
|
|
|
||||||
рассмотрим |
анализ |
точности |
|
|
|
|
|
|||||||
системы, |
схема |
которой |
представ |
|
|
|
|
|
||||||
лена на |
рис. |
1.10. |
из |
безынер |
|
|
|
|
|
|||||
|
Система |
состоят |
|
|
|
|
|
|||||||
ционного |
измерителя |
дискретного |
|
|
|
|
|
|||||||
действия, |
линейного |
нестационар |
|
|
|
|
|
|||||||
ного дискретного фильтра, |
линей |
Рис. |
1.9. Зависимость выиг |
|||||||||||
ного экстраіполятора |
непрерывного |
|||||||||||||
рыша |
тіо от |
вероятности р. |
||||||||||||
действия |
и нелинейного |
стацио |
||||||||||||
нарного объекта |
управления. |
|
|
|
|
|
|
|||||||
|
Безынерционный измеритель дискретного действия производит |
|||||||||||||
измерение полезного воздействия на систему U(t) |
в дискретные мо |
|||||||||||||
менты времени th—kT, где |
k — целое число, |
а |
Т — период дискрет |
|||||||||||
ности. |
|
|
|
|
|
|
|
|
|
|
|
|
||
Рис. 1.10. Схема исследуемой системы.
Полезное воздействие U(t) является линейной функцией времени со случайной скоростью изменения, имеющей равновероятный закон распределения, т. е.
U(t) = Üi. |
(1.143) |
где Ü — случайная величина, равновероятно распределенная |
в преде |
лах от — Ümax до Ü„ , a X 9 |
|
53
Измерение полезного воздействия |
производится с |
ошибками |
V(t],) (вредное воздействие). Значения |
V(tk) подчиняются |
нормаль |
ному закону распределения, в различные моменты времени незави симы, имеют нулевые математические ожидания и одинаковые сред ние квадратические отклонения Gv-
Таким образом, измеренные значения полезного воздействия
будут |
|
W(th) = U(tk) + V(tk). |
(1.144) |
Линейный нестационарный фильтр осуществляет определение скорости U полезного воздействия по значениям W(h) методом наи
меньших квадратов. Найденное значение скорости U(th) на выходе такого фильтра будет
k
б £ iW{tt)
= k ( k + \ ) (2 è + I)- ' |
(1.145) |
Линейный экстраполятор непрерывного действия на основании значе
ния Ü (th) выбрасывает |
воздействие |
X |
(t) на систему управления |
|||||||
объектом |
|
|
|
|
|
|
|
|
|
|
|
X ( t ) = Ü ( t h)t |
при tk< |
t < t h+1. |
|
(1.146) |
|||||
e{t) |
Управление объектом производится с момента t0 |
по |
разности |
|||||||
воздействия X(t) |
и координаты объекта X0(t): |
|
|
|||||||
|
|
z{t)=X[t) —XB(t). |
|
|
(1.147) |
|||||
|
Уравнения нелинейного стационарного объекта имеют следую |
|||||||||
щий вид: |
|
|
|
|
|
|
|
|
|
|
|
T0X 0.(t) + |
X 0 (t) = |
F[e(t)], |
|
(1.148) |
|||||
|
|
— В |
|
при е (t) |
|
— b, |
|
|
|
|
|
|
В |
е (t) |
|
|
|
|
|
|
|
|
F [•(<)] = |
- у |
при — й < |
е (/) |
&, |
|
(1.149) |
|||
|
|
В |
|
|
при е (t) |
|
b. |
|
|
|
где |
7о — постоянная |
времени |
объекта; |
В и Ь— параметры, |
характе |
|||||
ризующие нелинейность объекта. |
|
|
системы является дис |
|||||||
|
Искомой вероятностной |
характеристикой |
||||||||
персия ошибки системы, |
т. |
|
е. разности |
У |
полезного |
воздействия |
||||
U(t) |
и координаты объекта |
А'о(0 на момент |
времени |
/V |
|
|||||
|
|
|
Y=U(t ' o) - X0(t'o). |
|
|
(1.150) |
||||
Из анализа воздействий и системы следует, что математическое ожидание У равно нулю. Поэтому искомая дисперсия X=D у величи-
54
Ны У будет определяться соотношениями
Л = Dy = M[R],
(1.151)
Я = У2.
Приведенная система является нелинейной и нестационарной. Она включает в себя как элементы непрерывного, так и дискретного дей ствия. При таких условиях единственно целесообразным методом исследования является метод статистических испытаний с примене нием по мере возможности аналитических методов.
В связи с этим рассмотрим применение метода коррелированных процессов.
Упрощенную систему построим из элементарных систем, имею щих одинаковую структуру, но различные параметры. При выборе структуры элементарной системы будем руководствоваться следую щими соображениями. Ошибка системы У зависит как от полезного, так и от вреднего воздействия. Если не учитывать вредного воздей ствия и пренебречь нелинейностью системы, то ошибка будет про
порциональна скорости U. Дискретный нестационарный фильтр про изводит сглаживание вредного воздействия, и, по-видимому, это сглаживание является одним из существенных преобразований дан ного воздействия в системе. Исходя из этих общих сведений об исходной системе, элементарные системы, входящие в упрощенную, будем строить таким образом, что их ошибки будут равны
Zt = c tÜ - X v (t't), |
(1.152) |
где Xy(t'o) — часть воздействия X(t) на систему, обусловленная вредным воздействием V(/); с4— параметр элементарной системы.
Учитывая линейность уравнений (1.144), (1.145) и (1.146) исход ной системы и свойство нестационарного дискретного фильтра, эти уравнения можно заменить на следующие эквивалентные уравнения:
|
k |
ІѴ (h) |
|
\ |
|
|
6 £ |
|
|
|
|
Ѵ {th)^ ~ k ( k + U { 2 k + \ ) |
’ |
I |
(1.153) |
||
x v {t) = |
[V (th) t |
при th< |
/ < / Wx, |
|
|
X{ t) = |
U(t) + X v (t). |
|
|
|
|
Данная замена удобна тем, |
что величина |
Іѵ (І'о), |
входящая |
||
в упрощенную систему, вычисляется в исходной системе. Таким обра зом, для моделирования упрощенной системы, по существу, не тре буется дополнительных расчетов. Следовательно, построенная ука занным способом упрощенная система, состоящая из ряда элемен
тарных систем, практически являетея |
частью исходной системы. |
||
В качестве вероятностных характеристик упрощенной системы |
|||
взяты дисперсии случайных величин Zi, т. е. |
|
||
— М (6Д, 1 |
(1.154) |
||
Si = z f . |
I |
||
|
|||
55
Из (1.143) и |
(1.153) нетрудно получить, чго если і'п = ІТ, где |
I — целое число, то |
, |
|
|
|
|
6/<j2 |
|
|
|
|
Ü2 |
1 |
|
Р-г |
|
|
V |
|
|
+ |
< |
|
|
|
(* + І)(2/+1) |
3 |
|
|||||||
|
|
|
||||||||
Ksгs г - 2 |
|
6h2v |
|
|
|
|
|
6 / o |
||
(/ + |
1) (2/+ |
1, |
|
+ |
4c‘ |
(/ + |
1 ) ( 2 / + I) X |
|||
|
|
|
02 |
I |
A -----[A |
|
|
|||
|
|
X |
^max |
|
(І.155) |
|||||
|
|
|
3 |
ь |
<( |
4 5 |
u max’ |
|||
|
^S.S. — 2 |
|
|
6/02 |
|
+ 4etfj X |
||||
|
(/+1) (2/+ 1) |
|||||||||
|
|
I 3 |
|
|
|
|||||
|
|
G/o |
|
|
Ü |
|
|
|
|
|
X |
|
|
|
|
’-'max |
„2 „ 2 . |
Ul |
|||
(/ + |
1) (2/ + 1) ~ |
|
3 |
|
||||||
|
|
+ Ci ci |
45 umax- |
|||||||
Расчеты производились при |
следующих значениях параметров: |
|||
Ün.ax = - f > |
0Р = 1 |
; |
Т'0 = 1,57'; 6 = 2 , 5 ; |
|
В |
1,5 |
■ t0 = |
27'; t \ = \ 7 T . |
|
|
||||
Поскольку для определения статистических значений вероятност ных характеристик упрощенной системы практически не требуется дополнительных расчетов, то было взято несколько значений пара метра Сі упрощенной системы от с,=0,2Г до с{=107’.
Рис. 1.11. Зависимость выигрыша т]0 от отношения с/Т.
Было проведено N=100 экспериментов с исходной и упрощенной системами. По результатам этих экспериментов определено статисти ческое значение квадрата коэффициента корреляции rRS в зависимо сти от коэффициента с при заданном значении Т. Значение rRS в со ответствии с (1.37) позволяет найти выигрыш г|0 от применения ме тода коррелированных процессов. Полученная зависимость выигрыша г)о от коэффициента с представлена на рис. 1.11. Из этого рисунка
56
видно, что в довольно широких пределах изменения коэффициента с выигрыш г)о является существенным.
Оценка Лоі вероятностной характеристики системы вычислялась для двух упрощенных систем. В первом случае упрощенная система состояла из одной элементарной (т—1), во втором — из двух ( т = —2). При т = 1 параметр с, упрощенной системы принимался равным
2Т, а при т = 2 Сі = 2Г и |
= 0,27’. |
|
в табл. |
1.2. |
|
|||
Результаты |
моделирования |
и расчетов сведены |
|
|||||
|
|
|
|
|
|
Таблица |
1.2 |
|
т |
1* |
У |
f*a |
1V |
N |
^OI |
Т,о |
|
|
|
|
|
|
|
|
|
|
1 |
1,139 |
1,617 |
1,495 |
___ |
— |
1,038 |
7,6 |
|
2 |
1,139 |
1,617 |
1,495 |
0,156 |
0,175 |
1,044 |
9,2 |
|
Из |
таблицы |
видно, |
что, несмотря на |
существенное различие |
ве |
|||
роятностных характеристик исходной и упрощенной систем, выигрыш в точности или в числе экспериментов является большим.
Так при т = 2 выигрыш г)0=9,2, т. е. при 100 проведенных экспе риментах точность оценки А<и соответствует примерно 920 экспери ментам. В рассматриваемом примере выигрыш г|0 практически равен выигрышу в стоимости исследования, так как упрощенная система, по существу, является составной частью исходной.
1.9.О возможных путях развития метода
Взаключение настоящей главы рассмотрим возмож ные пути развития изложенного выше метода опреде
ления вероятностных характеристик, использующего совместно результаты статистических испытаний и ана литического упрощенного исследования.
Полученная оценка была оптимальной в классе ли нейных оценок. Очевидно, что возможности этого метода расширились бы, а его эффективность возросла, если рассмотреть и нелинейные оценки. Однако здесь не сле дует увлекаться, так как более сложные нелинейные оценки будут требовать для их реализации определения из экспериментов более сложных взаимных вероятност ных характеристик исходной и упрощенной систем, чем корреляционная матрица KRS*, которая требовалась для
изложенного метода. Так, например, могут понадобиться статистические значения моментов более высокого по рядка, чем второго. Точность лее статистических значе ний моментов ухудшается с ростом порядка момента. Поэтому оптимальная нелинейная оценка после замены
57
в ней точных значений необходимых вероятностных характеристик на статистические может настолько ухуд шиться по точности, что практическое использование ее станет нецелесообразным.
В качестве первого шага на пути получения нелиней ных оценок может быть, например, построение нелиней ной оценки, зависящей линейно от К* и полинома второй
степени от компонент вектора ц*—р.
Что может дать нелинейная оценка? Во-первых, бу дучи более общей, чем линейная, и включая в себя в частном случае линейную, нелинейная оценка в прин ципе может быть более точной. Во-вторых, для нелиней ных оценок, по-видимому, в меньшей степени будет критичным не совсем удачное построение упрощенной системы. Тем самым при получении одинаковой точ ности линейной и нелинейной оценок можно в послед нем случае облегчить построение упрощенной системы.
Другой возможный путь развития предложенного метода заключается в совершенствовании методики по строения упрощенных систем.
По-видимому, целесообразно для определенных клас сов исходных систем создать набор типовых элементар ных систем, из которых рационально строить упрощен ные системы, и создать более совершенную методику отбора этих элементарных систем в процессе статисти ческих испытаний.
Г л а в а 2
ОПРЕДЕЛЕНИЕ ВЕРОЯТНОСТНЫХ ХАРАКТЕРИСТИК МОДЕЛИ СИСТЕМЫ ПРИ ЧАСТИЧНОМ ПРИМЕНЕНИИ АНАЛИТИЧЕСКИХ МЕТОДОВ. ИСПОЛЬЗОВАНИЕ АПРИОРНЫХ СВЕДЕНИЙ
2.1. Возможные пути частичного применения аналитических методов и использование априорных сведений
В предыдущей главе был рассмотрен один из спосо бов применения аналитических методов для оценки вероятностных характеристик модели системы, заклю чающийся в использовании результатов аналитических исследований не самой исходной системы (модели),
58
а некоторой упрощенной системы. В данной главе рас сматриваются возможные пути применения аналитиче ских методов непосредственно к исследованию исходной системы.
Будем предполагать, что вероятностные характери стики исходной системы не могут быть найдены чисто аналитическим путем. Поэтому по отношению к иссле дованию исходной системы можно говорить только лишь о частичном применении аналитических методов наряду с методом статистических испытаний.
Стремление хотя бы частично решить задачу опре деления вероятностных характеристик аналитическими методами проявляется практически при каждом иссле довании. Это объясняется преимуществами аналитиче ских методов, дающих более наглядные и общие резуль таты, чем метод статистических испытаний, и позволяю щих в ряде случаев существенно упростить задачу, исследуемую при статистических испытаниях.
Идея совместного использования аналитических ме тодов и метода статистических испытаний не является новой [16, 29]. Однако ввиду отсутствия соответствую щих обоснований совместно использование этих методов для исследования исходной системы обычно произво дится в значительно более узких пределах, чем это возможно. Кроме того, отсутствие представления о воз можном выигрыше хотя бы при частичном использова нии аналитических методов и необходимость применения в конечном счете метода статистических испытаний порой приводят исследователя к неправильному реше нию ограничиться только лишь этим методом с потерей всех преимуществ от совместного использования обоих методов.
Наиболее широкое развитие частичного использова ния аналитических методов заключается в их примене нии к части исследуемой системы. Обычно исходная система разбивается на ряд более простых подсистем. Это разбиение часто соответствует сложившемуся раз делению труда у группы исследователей, когда каждая из подгрупп занимается определенными элементами системы и физическими процессами. При этом каждая из подгрупп исследователей, глубже понимая и всесто ронне изучая процессы в части исследуемой систе мы, может значительно лучше найти возможные анали тические методы исследования.
59
