книги из ГПНТБ / Солодовников В.В. Статистическая динамика линейных систем автоматического управления
.pdf476 |
ТОЧНОСТЬ ЛИНЕЙНЫХ |
СИСТЕМ С ПЕРЕМЕННЫМИ |
ПАРАМЕТРАМИ |
[ГЛ. X |
||
то можно формально записать, что |
|
|
|
|
||
|
f(t) = |
L~1(p, t)K {р, t) m (t), |
|
|
||
где |
L - 1 (p,t) — оператор, обратный L(p,t) . |
|
|
|||
|
При этом исходное уравнение (9.2) превращается в |
|
||||
|
D(p, t )x(t) = |
М (р , t)L~* (р, t)K(p, |
t)m(t), |
(10.73) |
||
где m(t) — «белый» шум, |
или в |
|
|
|
|
|
где |
D(p, t) х (() = M l (р, t) т ((), |
|
(10.74) |
|||
|
|
|
|
|
|
|
|
M , { p J ) = M { p , t ) L - 1{p, |
t ) K ( p , t). |
|
|||
|
Это уравнение формально совпадает с (10.2) и его корреляцион |
|||||
ная |
функция может быть |
исследована аналогично изложенному. П-ри |
||||
|
|
этом, очевидно, необходимо ввести |
||||
|
|
новые |
начальные условия. |
|
||
|
|
Фактически речь идет о последо |
||||
|
|
вательном соединении заданной си |
||||
|
|
стемы |
и |
формирующего |
фильтра |
|
|
Рис. 10.6. |
(рис. 10.6) и |
исследовании всей |
|||
|
системы |
в целом, как имеющей на |
||||
|
|
входе |
«белый» |
шум. |
|
|
Взаключение необходимо отметить, что изложенные здесь методы целиком применимы к линейным системам с постоянными парамет рами и, в частности, к исследованию динамической точности этих систем непосредственно после поступления на вход системы случай ного процесса.
7.Методы определения дисперсии и корреляционной функции случайных процессов на выходе линейных
систем с переменными параметрами
Внастоящее время распространены два метода решения задачи
анализа динамической точности систем автоматического |
управления, |
т. е. методы определения дисперсии и корреляционной |
функции на |
выходе заданной системы при случайном входном сигнале с известной корреляционной функцией.
Первый метод — аналитический. Он заключается в определении моментов распределения вероятностей ошибок системы методом ре шения дифференциального уравнения, описывающего систему, или вычисления интегралов типа рассмотренных в § 10 гл. IV по извест ной спектральной плотности входного сигнала и передаточной функ
ции |
системы. |
Серьезным недостатком этого метода является |
то, что |
|
он |
применим |
к узкому |
кругу задач, включающих в себя в основном |
|
линейные постоянные |
системы при стационарных случайных |
воздей |
||
7] |
ОПРЕДЕЛЕНИЕ ДИСПЕРСИИ И КОРРЕЛЯЦИОННОЙ ФУНКЦИИ |
4 7 7 |
|
ствиях. |
Но даже при этом в случае систем высокого порядка он при |
||
водит к сложным вычислениям и методам численного |
решения. |
|
|
Второй метод — экспериментальный. Он состоит |
в возбуждении |
||
исследуемой системы случайными воздействиями заданного типа, полу чаемыми при помощи различных генераторов шума и эксперименталь ной обработки случайного сигнала на выходе. В случае постоянных систем и стационарных воздействий эта задача решается при помощи коррелятора, вычисляющего корреляционную функцию стационарного процесса.
Если же воздействие нестационарно или система имеет перемен ные параметры, необходимо проводить запись достаточно большого числа реализаций выходного процесса и применять осреднение по совокупности. Недостатками этого метода являются необходимость большого количества данных для получения надежных результатов, сложность применяемой аппаратуры и громоздкость обработки.
Однако возможен и несколько иной подход к этой задаче. Он со стоит в том, чтобы, зная дифференциальное уравнение системы и характеристики случайного сигнала на входе, свести задачу к случаю «белого» шума на входе и получить соотношения, определяющие характеристики выходного случайного процесса через передаточную функцию системы. Мы применим этот метод к задаче определения корреляционной функции и дисперсии на выходе линейных систем с переменными параметрами при случайных воздействиях рассмотрен ного выше класса. Из предыдущего параграфа очевидно, что принци пиально достаточно рассмотреть только случай «белого» шума на входе системы, так как метод формирующих фильтров позволяет обобщить его на широкий класс стационарных и нестационарных случайных воздействий.
Итак, предполагая, что система описывается дифференциальным уравнением (10.1) и на вход ее подан «белый» шум, поставим задачу определения корреляционной функции Rx {t, т) нестационарного про цесса на выходе этой системы.
Выше было показано, что Rx (t, т) удовлетворяет дифференциаль ным уравнениям (10.51) — (10.54). Следовательно, фактически задача сводится к решению одной пары этих уравнений (вследствие сим метричности Rx (t, т) относительно плоскости t = x), порядок которых
равен порядку исследуемой системы и которые не включают в себя случайных функций. Эти уравнения могут быть решены приближен ными методами, рассмотренными в §§ 3 и 4 настоящей главы. Так как эти методы уже были рассмотрены выше, мы изложим здесь метод решения этих уравнений с применением моделирующих устройств. Целесообразно рассмотреть моделирование уравнений (10.51) и (10.52) или, что совершенно эквивалентно, (10.53) и (10.54), так как реше
ние каждого |
из уравнений, входящих в эти пары, имеет свои осо |
бенности, |
. . |
478 |
ТОЧНОСТЬ ЛИНЕЙНЫХ |
с и с т е м с п е р е м е н н ы м и п а р а м е т р а м и [ г л . X |
Рассмотрим вначале метод решения однородного дифференциаль |
||
ного |
уравнения (10.51). |
|
Основным вопросом, который необходимо решить при моделиро |
||
вании |
этого уравнения, |
является вопрос об определении начальных |
условий. Он сводится к |
определению дисперсии и производных от |
|
корреляционной функции |
в плоскости t — x. |
|
Как было показано в § 5, для рассматриваемого класса случайных
процессов эти производные существуют |
и непрерывны до ( 2 я — 1)-го |
|||
порядка |
во |
всей |
области / > х и t < |
х и имеют скачки при t — т, |
которые |
могут |
быть зафиксированы. |
Следовательно, производные |
|
от R t ( t , x ) |
существуют при подходе справа к плоскости t — x. Так |
|||
как уравнение (10.51) определяет корреляционную функцию в области t > х, то именно эти значения корреляционной функции и ее произ
водных выбираются |
в |
качестве |
начальных |
условий |
для |
его |
ре |
||||||
шения. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Один |
из возможных |
методов:) определения этих |
величин с |
при |
||||||||
менением |
моделей |
состоит в следующем. |
|
|
|
|
|||||||
|
Рассмотрим систему |
уравнений |
|
|
|
|
|
|
|||||
d. |
' d‘x (t) dJx (t) |
|
|
dl+ix(t) |
dhr (x) |
|
|
|
|
||||
dt |
dt1 |
d-J |
|
|
|
|
dti+ l |
dxi |
|
|
|
|
|
|
|
dlx (t) |
|
du |
'x(x) |
|
|
l, j — 0, |
1, . . . . |
n — 1, |
(10.75) |
||
|
|
dtl |
|
dxl+l |
-= t- s |
|
|||||||
|
|
|
|
|
|
|
|
|
|||||
гак |
как |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
d‘x (t) |
|
|
dl+ iX(x) |
|
|
|
|
|
|
|
|
|
dt |
dxl |
-= /-8 |
|
dxi+l |
|
|
|
|
|
где |
s > 0 — малая величина. |
|
|
|
|
|
|
||||||
|
Получающиеся при этом члены с /t-ми производными могут быть |
||||||||||||
выражены |
через |
|
производные от |
х |
более низких порядков |
из урав |
|||||||
нения (10.1): |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
dxn |
|
1 |
|
|
|
|
|
|
|
(10.76) |
|
|
|
dt'1 |
An (9 |
-fc =0 |
|
|
(1 =0 |
|
|
||||
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|||||
Для получения дисперсии и производных от корреляционной функции Rx (t, х) при t — x+ необходимо усреднить уравнение (10.75) по совокупности. При этом возникает задача определения средних
J) См., например, Du n c a n D., Response of Linear Time — Dependent Systems to Random Inputs, Journal of Applied Physics, t . 24, № 5, 1953 и книгу, цитированную на стр. 464.
480 ТОЧНОСТЬ ЛИНЕЙНЫХ СИСТЕМ С ПЕРЕМЕННЫМИ ПАРАМЕТРАМИ [ГЛ. X
производных корреляционной |
функции |
при |
т — t : |
|
|
|||||||||||
|
|
|
|
= |
R l+1,l-\-RJx l+1. |
i. J = |
о, |
1.......... Я — |
2, |
|
||||||
|
d ^ и—i, / _^ л—i, /+1 |
|
1 |
n- 1 |
|
|
|
|
|
|
||||||
|
|
a*(t)Rx |
, |
|
|
|
||||||||||
dt |
|
|
|
|
|
|
(0 |
|
|
|
||||||
|
|
|
|
|
|
A’ =0 |
1= 0, 1....... ft — 2, |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||
— |
pi |
|
|
|
|
|
|
|
|
|
diL+ J'k(t, t) |
|
|
|||
dt Hx |
|
|
|
|
|
|
la«0 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(10.81) |
||
|
|
|
|
|
|
//-1 |
|
|
\ |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
• |
У = |
0 ,1 .........ft- 2 |
, |
|||
|
|
|
|
|
|
*=o |
|
|
J |
|
|
|
|
|
|
|
d |
|
|
я-1 __ |
1 |
{ т |
|
|
d7+V-1 |
(/, |
t ) |
|
|
||||
dt |
Kx |
|
~ |
e„ (0 |
|
|
|
|
dtn~' |
|
k |
|
|
|||
|
l|i = 0 |
|
|
|
|
|
- , r |
|
|
|||||||
|
|
|
|
|
|
|
|
- Ц « * ( о [ я |
|
|
|
|
||||
где |
|
обозначено /?у |
= |
— :— r R r (t,t) |
, = r |
|
|
|
|
|||||||
|
|
|
|
|
|
|
dVd-d |
|
-vV |
|
|
|
|
|||
|
|
Если (t - |- у)-ая производная от корреляционной функции непрерывна |
||||||||||||||
при |
t = |
z, |
то |
R x‘ j = |
Rx |
и |
число уравнений уменьшается. Отметим, |
|||||||||
что |
если |
т = |
0, т. е. оператор уравнения |
(10.1) |
М(р, t) = b0(t), |
|||||||||||
то |
все производные до |
(2п — 1)-го порядка |
от Rx (t, т) |
непрерывны |
||||||||||||
при |
t = |
т, |
а производные |
от k{t, т) |
при t = |
т до (и — 2)-го порядка |
||||||||||
обращаются в нуль |
согласно |
(10.29). |
В этом |
случае система состоит |
||||||||||||
из |
|
л (п-|-1) |
|
„ |
вида |
|
|
|
|
|
|
|
|
|||
|
— |
—- уравнений |
|
|
|
|
|
|
|
|
||||||
|
|
J L |
r J‘ = |
|
+ |
|
|
|
t , j = 0, |
1...........ft — |
2, |
|
||||
JLdt R * i 1= R*+1,n |
|
Л |
S VСОЯ+* |
|
|
(10.82) |
||||||||||
|
|
|
|
|
|
|
|
ft= 0 |
i — 0 , * |
1.............n — |
||||||
|
|
|
|
|
|
|
|
|
|
|
2 , |
|||||
d |
p/i-i, л—l |
(0 |
|
|
^ |
|
|
|
|
|
|
|
|
|||
|
|
* |
|
|
4 ( 0 |
" |
an (0 |
fc=o |
|
|
|
|
|
|
||
Преимущество этого метода заключается в том, что он позволяет
определить R lJ в виде непрерывных функций времени и приспособ лен для использования моделирующих устройств.
9] |
ЛИНЕЙНЫЕ СИСТЕМЫ СО СЛУЧАЙНЫМИ ПАРАМЕТРАМИ |
485 |
||||
|
|
|||||
Таким |
образом, k{t, |
t — 0) как |
функция 0 может быть получена |
|||
путем последовательного |
соединения |
инверсных |
исследуемой системы |
|||
и формирующего фильтра (см. рис. |
10.7). |
|
|
|
||
Отметим, что для определения |
kn{t, t — 0) |
более |
целесообразно |
|||
решать на |
моделирующей установке не уравнение |
с воздействием |
||||
в виде дельта-функции, а однородное уравнение с начальными усло
виями |
вида |
(10.27). |
|
|
3. |
Заключительный этап решения задачи |
состоит в |
возведени |
|
в квадрат |
k { t , t — 0) и вычислении дисперсии |
по формуле |
(10.83). |
|
Повторяя процесс решения для различных фиксированных значений t, можно исследовать изменение дисперсии со временем.
9. Некоторые вопросы анализа линейных систем со случайными параметрами
Ранее мы рассматривали системы, параметры которых изменялись со временем вполне определенным образом, а входные сигналы явля лись случайными процессами. Однако существенный практический интерес представляет более общая постановка задачи. Дело в том, что большинство систем автоматического управления имеет параметры, которые меняются со временем в той или иной мере случайно. Дей ствительно, все элементы систем создаются с техническими допусками и поэтому должны рассматриваться как имеющие случайные пара метры в этих пределах. Кроме того, случайные параметры имеют системы, характеристики которых изменяются в зависимости от мгно венных значений входных сигналов. В радиотехнике, например, моду лирующие системы, а также систему автоматической регулировки
усиления иногда удобно рассматривать как системы |
со |
случайными |
|||||
параметрами. |
|
|
|
|
|
|
|
Отсюда следует, что общая трактовка вопросов |
анализа |
систем |
|||||
автоматического |
управления при |
случайных воздействиях |
состоит |
||||
в предположении, |
что не только входной сигнал, |
но и параметры |
|||||
самой системы изменяются случайным образом. |
|
|
|
|
|||
Отметим, что |
математически |
изучение таких |
систем |
заключается |
|||
в исследовании |
статистических свойств решений |
дифференциальных |
|||||
уравнений с коэффициентами, которые являются случайными процес
сами. Поэтому |
методы, применяемые при анализе |
линейных |
систем |
с переменными |
параметрами, не всегда пригодны |
для этих |
систем, |
так как они предполагают точное знание параметров как функций времени.
Рассмотрим некоторые вопросы анализа выходных процессов систем этого класса.
Итак, предположим, что параметры исследуемой линейной системы изменяются во времени случайным образом и к системе приложено случайное воздействие. Мы ограничимся в дальнейшем для простоты
