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

книги из ГПНТБ / Марков В.Н. Теория управления (устойчивость, стабилизация, оценки) учебное пособие

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

 

 

 

 

-<^1X1

-W ti

 

 

 

'

 

Q-1

'

•(-%&

■+

)

f

 

где oi^

 

 

 

 

^

 

и

 

" абсолютные значения

мнимых и вещественных

частей

корней полинома

Q z x (w ) .

 

 

 

 

Наиболее

часто используются для

аппроксимации

к

 

выражения вида

 

 

 

 

 

 

Х-ЧХ1 Л1

 

 

 

 

 

•*|tl

<3«.

. &

4

с

 

e ' Cwyxt

iih.p.itij^O'

Cbonmj

 

, -

-

 

Учитывая условие физической осуществимости и применяя к

(I0 .I2)

преоразование Фурье получим:

 

 

 

 

 

 

 

 

 

 

m

 

 

(10.17)

Для одномерного случая

будем иметь

 

 

 

 

 

 

 

S x ^ ^ v o ) -

 

 

 

. (ЮЛ8)

Подобно

тому,

как было получено (10.15), (10.17)

можно предс­

тавить в виде:

./

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

2>х ХЦ 4 )

 

 

(ЮЛ9)

Проиллюстрируем теперь использование данных'выводов и

соотношений на примерз анализа простой одномерной системы.

Пример.

Пусть

система представляв^ собой апериодическое зве­

но (например \С,~ цепочка) с постоянной времени <"V

й общим

коэффициентом усилений

(см. рис.12).

 

 

 

X (t)

а д .

к

 

TS +

i .

Рис. 12

- ТЗО -

На входа этого звена начиная с момента 'Ь* действует возмущение » представляющее собой стационарную слу­ чайную функцию с математическим ожиданием М. х и корреля­ ционной функцией

*СзЛ** ~^ь) *•

^

До момента ~^о

система находилась в-покое. Найти

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

Yu,.

В соответствии

c-(lu .2),

(10.3) и учитывая,^что

 

U -U .-y -

4

о,~

можно записать

 

 

 

tov^U) «■

]

- *9

d t =•

ч m

 

т

Чг —

 

 

Т

£

т Д-С * V m » U - Y TJ .

Выражение для дисперсии непосредственно получим из выражения

(ЮЛ) для корреляционной функции: * *

>9UJ * Ц ,и9 = 1]

*j

g(.+- X) g-U- «J К*. u - »j <U * t

 

 

4%

4

^

 

 

V

^

Ac-4

_оцг-©1

 

-k.- о-г

 

_ _ _

~-=f~

1)5,19

 

Q.

d

£

ДV olts

 

4%•**

 

 

 

 

 

А»

 

 

 

 

 

lc

г -1*

 

V

-<*'**-*i cAo екх.

* Г

°W ^

 

- ТЗ Г -

Проинтегрируем сначала по *9

*

для чего разобьем

интервал

интегрирования по переменной

9

на две

части:

{-to

)

и ( ^ I . ). Тогда будем иметь

 

 

 

 

г*

 

Ч ? -Alr-e)

{

 

г

Jl е>

 

d (« tW r

 

$ * ] h

V*) - -|-V WX

 

 

 

 

 

*fc.

 

 

 

 

 

Долее процесс интегрирования не представляет труда и оконча­ тельно получим:

I

•+^iX-v~ "Ьэ)

 

 

 

 

-t

 

 

Как нами отмечалось ранее, здесь

мы имеем дело

с неуста^

повившимся режимом, что и выразилось

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

 

и

от времени. Устремляя

в этих выражениях -Ь

ко©

,

получим:

 

 

 

 

to'j (.^ * lonx =e/*vt^

 

 

 

Таким образом, в .установившемся режиме на, выходе будем

 

иметь стационарный процесс.

Установившееся значение

диспер­

 

сии можно получить другим путем. На основании (10.14) и (9.20),

учитывая, что

З ^ К ^ О о )

,

для установившегося

значения

 

получим:

 

 

с*

J

 

 

 

 

 

 

\

^

 

 

 

 

 

 

 

 

I

 

i

1 v i

р

В соответствии

о

(9.19) при

к * о 9 =

v _-<Ч*1 „

г сг^ы

 

 

 

тп:

 

В рассматриваемом примере

 

 

*оМ Ыг .

 

 

 

 

■г

\Кс.

_

kv*

( . .

2.''С O'*

 

№ Цч)\-1-VT.iUJ

'

л+Vu' >

ъ аы 3 '

h*x\d')&x^*\)

 

Тогда

 

 

. OQ

4

 

 

 

 

W^cri oi

 

 

dv>

 

 

3) м -

тг

 

 

trV + \)

 

 

 

 

 

 

 

 

 

- тз? _

Интеграл вычисляется стандартными приемами разложения подин­

 

тегральной функции на простые дроби.

Интегралы от дробно-ра­

 

циональных функций часто встречаются в практике решения за­

 

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

специаль­

 

ные таблицы формул,Естественно, результат совпадает с получен

 

иым

выше.

 

 

 

 

 

 

 

 

 

Отметим, что если в системе, заданной дифференциальным

 

уравнением требуется учесть случайные начальные условия то

 

поступают

следующим образом:

 

 

 

 

 

 

а) Определяют математическое ожидание выхода, решая не­

 

однородную систему с правой частью,

соответствующей матема­

 

тическому

ожиданию входной случайной функции, при начальных

 

условиях,

определяемых их математическими ожиданиями.

 

г

б) Получают решение

однородной системы, задаваясь началь

 

ными условиями в общем виде.

При этом эти решения

являются

 

линейной функцией начальных условий.

Полагая начальные ус­

 

ловия случайными величинами, вычисляют дисперсии и корреля-

 

цион-ные функции для решений однородной системы.

 

 

в) Определяют дисперсию и корреляционную функцию выхода,

 

•обусловленные действием 'входного случайного сигнала при

 

нулевых начальных условиях.

 

 

 

 

 

 

г) Для определения дисперсий и.корреляционных функций пе

 

ременных,

соответствующих суммарному движению системы,следу­

 

ет воспользоваться формулами -лла дисперсии и корреляционной

 

пункции суммы случайных

,уикций.

 

 

 

 

 

Пусть в нашем примере имеют место

начальные условия в -

 

м о м е н т

“ t “

' t o

( с и с т е м а

I - г о

п о р я д к а ,

з а д

о т к л о н е н и е ) ,

п р н ч . ы ,

N 1

 

и

з а д а н ы .

- 133 -

х

S

-имеет решение

*

л

/I

V

*-

Л

,

d

*L

^ где Л

 

 

 

 

__ -fc-v

 

 

 

 

.с £

"т’’ , $.>•

-«-

 

Тогда

 

 

 

 

 

 

 

 

-

Y *

 

 

 

Т

* Ч 4^ « № » Л

 

Т > , 3 > ^ * 1 > * Л

 

и ч*

Для определения

результирующих значений

 

 

в случае некоррелированности входного сигнала со случайным начальным состоянием системы следует сложить flv и 3) слу­

чайного сигнала выхода от каждого из рассмотренных факторов:

3L+J * «tjjL-t) +

,

+ \ 4t y

- т34 -

3. О преобразовании закона распределения. Если входио$нУме& гауссовское распределение, то, как известно, век­ тор на выходе линейной динамической системы также нормален.

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

Остановимся далее на качественной оценке влияния линейных

систем на трасформацию распределения входного вектора в ка­ нале вход-выход. Ограничимся стационарными системами и бу­

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

системы в момент^* О

нулевым. Рассмотрим систему,

представленную на рис. II. В

соответствии

с

(9.5)

-t

 

 

 

 

 

 

 

 

^f(>)

J S'Gfc-'XjXcy A t .

(10.20)

Разобьем

интервал интегрирования

 

 

на ^ интервалов

длиной

Д-

и обозначим

 

 

 

 

 

 

 

 

 

W

v y X t x M x

 

 

(Ю.21)

При достаточно

большом ^ и

малом

Ь

^ Хц,на интервале

 

 

 

можноМ О Г :

считать постоянным и равным

Х .иь) .

Тогда можно написать:

 

 

 

 

 

 

 

 

 

 

 

V t )

=

s-M

 

,

( 10. 22)

 

 

 

 

 

 

 

 

V

-

Х(^)

 

,

-

v r 'U ,

 

 

 

 

 

 

 

 

 

 

Если бы слагаемые здесь были независимы,

то на основании цент­

ральной предельной теоремы можно было бы

считать,

что закон

распределения ординат

случайной

функции

при росте ^

- Т35 -

стремится к нормальному*

Однако величины Ч являются зависимыми (они выражаются

через ординаты одной и той же случайной функции). Если бы ,

стремление

не приводило к увеличению зависимости

между

Сх \\х ,-^к!) * закон распределения "Чй^также

стремился

бы к нормальному. На самом деле в этом случае свя­

зь между слагаемыми только увеличивается и предельная тео­ рема для суммы (10.22) неприменима. Тем не менее никоторые выводы о свойствах закона распределения сделать можно.

Чем меньше зависимость между *ЧС С л, , т.0, чем быстрее убывают корреляционные моменты между составляю­

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

корреляционная

функция

, тем ближе будет закон pao- I

пределения

к нормальному.

 

В центральной предельной теореме кроме большого числа • слагаемых требуется еще однородность слагаемых. Наиболее благоприятным в этом отношении является случай одинакового

закона

распределения. У нас это будет

тогда, когда коэффи­

циенты Оч постоянны и не зависят от индекса,

т.е. когда

Q,(T) ^

c>n,s+ . В этом случае М\1Цлл)

определяется

равенством

 

 

 

 

 

 

Ш . Й - Д 'Х с - '& . ч )

 

 

UJ.2J

н.-аледисе означает, что рисе, атрпваеыая

система

имеет оеско-

печно узкую полосу.пропускания ь районе

нудеьс ;

частоты.

Таким осразом

получаем, что. чем уже полоса .пропускания, тек

с л.а

к нормальному закону бу;.ет закон

р.аспр-.у ели:ия аг

*36 -

нат функции на выхода. Выбрав для такой системы достаточно большим промежуток (о} -Ь ) можно получить сколь угодно близ

кий к нормальному закон распределения.

Итак, можно сделать следующий качественный вывод;

Проходи через линейную систему случайный процесс преобре-

тает свойства, сближающие его с нормально распределенным процессом, причем эффективность этой ’‘нормализации” будет

тем оольше, чем уже полоса пропускания; чем больше время

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

ного процесса.

С учетом начальных условий выход системы может быть

записан

в виде:

 

.

Vv.

 

 

 

. /V

 

 

 

1 U J

 

(10.2<t)

\с$

” независимые

интегралы однородного

дифференци­

ального

уравнения,

£4

- постоянные, линейно

связанные с

начальными условиями. Если начальные условия случайны то

случайны и Cjj. Однако не являются независимым» и сла­

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

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

распределен нормально.

Но даже в том случае, когда

он не

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

из сла­

гаемых не превалирует

над дисперсией "Х ш , то плотность

как правило,

оказывается ближе

к плотное**»

нор-

пильного распределен;},,

чем плотность

\

 

рС$) .

 

Следует отметить,

что приведенные

качественные

сооб­

ражения дают иногда основания для проведения в рамках vjp-

реляционной теории исследований езете.., czAC'y*:;Z2U

кейац э-зиенту.

§ П • Соотношения в пространстве состояний.

I . Дифференциальное .уравнение относительно ковари­ ационной матрицы. До смыслу переходной функции состояний (первое уравнение.(9.1)) и решения (9.2), очевидно, что если система (9.1) возбуждается чисто случайным процессом,

то вектор состояний

представляет собой марковский

случайный процесс.

 

 

--

Так как линейный

преобразования

гауссовского

вектора

сохраняет гауссовское

свойство (см.

например \i*t]

), то

гауссовский марковский случайный процесс всегда может бы­

ть представлен как вектор состояния непрерывной линейной

ы

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

 

 

 

 

X

<=■

 

. ( И Л )

где

XCt?

-

К

-

мерный вектор,

 

 

гН*)

-

Чх

-

мерный вектор белого шума,

 

 

 

 

М -1 Xlio}}

^ ftv-xl'bj

(П*2)

-№i[X(4.j-"ixt+UlXib)-

а 1C*.

(и.з)

 

 

 

 

 

 

Vn w l-9

( I I .-'О

.M jW b i) -

 

 

 

- »nw<4)Y - T lU O ^ V tjjc ii.V )

Блок

схема

такого

представления изображена на рис.

13.

- Т38 -

Применял

к ( I I . I )

операцию математического

ожидания,

получим дифференциальное

уравнение относительно

 

№ *(*) - Л(Л}тпх1^+

 

Ш *5)

 

 

Ж*. О»)

- задано.

 

Оощее решение

( I I .5) дает

зависимость

V(\x -

. iloo

мотрим, каким образом определится ковариационная матрица

такого процесса

К .

 

 

(

-п

 

К **

K,«>»V) -

MtXtylJU)] J .

(И.6)

Вычитая (11.3) из

(1 I.I)

и умножая

затем результат на

LXU) - tnx^

T- txtvjj

получим

 

 

 

Транспонируем полученное

уравнение

'Хи^СХ^З ■*Х^Ш"лL9 ч Х<-9 Mdlfв \t>• о■I•а)

Просуммируем (11,7)

и (11.8)

i l x i y t b j } =

t {X t+X i^kl* )

“VЪ U;

LX

11 Xoj Ife}1 u;

и применим операцию математического ожидания

К х = Ж ^ К х О К х/ и ? + S(.v)l!it) t J x 9 & W II,9 )

I

где

 

Ь-9 -- N l|X L*)[Jjrt^T1j

- uai[j:u:a ( ft. Xp)

взаимно"-; кор'.елыдо вектора состояил!- и вектора

воздействия. Леии.л>ауя .уа,..амслталнну^

матрицу с,:сь-.,ыы

, решение относите;, т о -

центр дрог з*;.ю V:, сое т

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