книги из ГПНТБ / Марков В.Н. Теория управления (устойчивость, стабилизация, оценки) учебное пособие
.pdf
|
|
|
|
-<^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 |
й общим |
||||||||
коэффициентом усилений |
1с |
(см. рис.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 ^ |
|
t«
- ТЗ Г -
Проинтегрируем сначала по *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:, сое т |