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

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

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

Пример 2. Пусть в исходной системе осуществляется усреднение с весовой функцией р{т)>0. Тогда импульсная переходная функция

системы будет

t 1

W (t , т) = р (т) J P M dll

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

Wy (t, т) = .

Воздействия на исходную и упрощенную системы — нестационар­ ный нормальный белый шум с корреляционной функцией

K ( t u /s)=C(<l)c(/2)Ö(<l-M.

где 6 (tit2)—ö-функция; 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,

 

 

 

 

 

 

 

 

 

 

* 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

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