
книги из ГПНТБ / Климов, В. А. Некоторые прикладные методы анализа и синтеза сложных автоматических систем с использованием ЦВМ
.pdf364
ставляюцей. |
В качестве |
примера изложим алгоритмы для |
^ -й со |
|
ставляющей. |
|
|
|
|
Предварительно фиксируется порядок предыдущей j |
~ J -й со |
|||
ставляющей. |
_ |
|
_ |
|
Вели |
то |
I, и если |
то^_,= |
2. |
Затем по соотношению (6.51) получаем |
|
|
|
|
i = |
' |
(6*54) |
|
После этого для j- |
составляющей вычисляется параметр p(j^) |
||||
по формуле (1.84) с учетом (6.54). |
|
|
|
||
Если р- -с I, то |
j составляющая имеет первый порядок и |
||||
ее передаточная функция соответствует |
(4.138). |
|
|
||
Если |
$» I, то j- |
составляющая имеет второй порядок и ее |
|||
передаточная функция соответствует (4.139). |
|
L |
|||
При использовании функций (4.138) |
и (4.139) |
в качестве |
|||
должно приниматься число (6.54). |
|
|
|
||
|
Последняя А составляющая |
|
|
||
Для того чтобы правильно осуществить окончание процесса |
|||||
формирования передаточных функций составляющих процессов, |
в |
ходе процедур для всех составляющих выполняется сравнение чи
сел |
L |
с величиной |
п . |
|
|
|
|
Если окажется для |
J -й составляющей |
|
|||
|
|
|
|
i - |
= /7» |
(6.55) |
то |
это означает, |
что данная составляющая не существует, и по |
||||
казывает, что предыдущая j |
~ I -я составляющая является послед |
|||||
ней ( |
А - составляющей).. |
Условие (6.53) используется как |
||||
признак окончания процессов. |
|
|||||
|
В использовании условия (6.55) может возникнуть затрудне |
|||||
ние. Действительно, если окажется для очередной J - 1 |
состав |
|||||
ляющей |
. |
. |
. |
|
|
|
|
|
|
|
I =■ i-j-i - п - 11 |
|
|
то это показывает,что |
эта составляющая является последней |
|||||
( |
Л - |
составляющей) |
и имеет первый порядок. Однако опреде |
лить затем величину (6.55) не представляется возможным, так как для J— 1 составляющей нельзя определить величину р (р ^ ) (один из коэффициентов в формуле будет иметь отрицательный ин
декс). Уто свидетельствует о том, что указанный коэффициент не существует. Поэтому принимается условно здесь
365
? j -i = 0,90, |
(6.56) |
Величине p^_f может быть присвоено и любое другое число, мень шее единицы.
х х
х
Таким образом, в результате расчетов по описанным выше
алгоритмам оказываются известными для каждой j -й составляю щей параметры и общий порядок i = Lj уравнений уже выделен ных составляющих. Это позволяет, как уже указывалось, выбрать формулы для передаточных функций составляющих и определить их
коэффициенты. Оказывается известным также общее число состав ляющих Л .
§ 3. МГОРИТШ ОПРЕДЕЛЕНИЯ ЗАПАСОВ УСТОЙЧИВОСТИ
Алгоритмы определения запасов устойчивости должны состоять из процедур вычисления параметров т 3 и процедур определения
параметров |
р |
для колебательных составляющих. |
Параметры |
должны вычисляться по соотношениям (4.140), |
|
а параметры |
р |
по формуле, которую мы запишем для промежуточ |
ной составляющей. Эта формула имеет вид
(6.57)
Формула (6.57) легко записывается по характеристическому уравнению для j -й колебательной составляющей, которая соот
ветствует знаменателю функции (4.139). При использовании фор
мулы (6.57) необходимо учитывать, что оценивать по этой фор муле колебательности для составляющих процессов можно при вы полнении исходной предпосылки метода, т.е. при условии, когда все параметры т^ равны или меньше единицы (4.140). Это усло вие будем записывать
(6.58)
Только при выполнении условия (6.58) могут применяться и все
другие излагаемые ниже алгоритмы.
371
Запишем также выражение для производной выходной коорди
наты. Из (6.79)
^вь/х,^ ~ 2е°*4 £(Z7oiCcojcos to£-(Cd+J?co)slncoi]. (6.84)
Для оценки качества процессов по колебательным составляю щим необходимо знать максимальные значения выходных координат
и максимальные скорости их изменения. Очевидно, что для опре
деления максимального значения координаты необходимо^найти
значения этой' координаты в трех точках
t = 0 ; |
t |
- t■zэ., |
|
* |
|
|
|
|
|
(ё-85) |
|
и выбрать из них наибольшее. На рис.6.3 |
показаны примеры,ког |
||||||||||
да максимальные значения соответствуют одной из трех точек |
|||||||||||
(6.85). |
|
|
|
|
|
|
|
|
|
|
|
Очевидно также, |
что для определения максимальной скорости |
||||||||||
изменения координаты х д |
-необходимо найти значения этой ско |
||||||||||
рости для трех точек |
/ |
- |
. |
|
. |
|
|
_ |
‘ |
|
|
|
5Г |
|
|
|
л |
( 6 . 8 6 ) |
|||||
t - 0 •, t - t 3r! |
2 |
'со |
И |
£ = £ |
|
_ . + - г - • г г • |
|||||
^ |
^э,1 + |
2 |
со |
|
|||||||
и выбрать из них наибольшее. На рис.6.4 |
показаны кривые,когда |
указанные максимальные скорости соответствуют одной из трех
точек (6.86). Вводится также условие, |
что для второй точки |
|||
цри |
|
|
|
|
^ ~ ^Э,1 |
1 - ^ - |
о |
(6.87) |
|
|
г сю |
|
|
|
значение х Вь1Х ^ае определяется, |
так как значение скорости в |
этой точке не имеет Физического смысла.
Из описания излагаемых алгоритмов следует,что для определе ния искомых характеристик качества процессов необходимо предва рительно определять точки t=t3 , и t = t3Z.Для этого воспользуем ся материалами, имеющимися в [9J.
Введем в рассмотрение комплексную функцию |
|
Ф ( t) = 2 е* t R* ( P i) j a t . |
(6 .88) |
W iW
Нетрудно видеть, что ф (f) цредставляет собой логарифмиче
скую спираль, для которой угол поворота радиуса-вектора равен oot • Причем перемещения проекции точек этой спирали по веще
ственной оси соответствуют процессу (6.79), точнее его пере-
373
y{t) = te1* (D + jc)(cosu)t +<fsincot). |
(6.89) |
Выполняя перемножения, находим
ф(£) = 2е [(Pcoscot - Csin cot)+j- (Ccosoot-<-l)si.ncot)j. (6.90)
Таким образом, на самом деле вещественная часть (6.88) совпа
дает с выражением для переменной части (6.79).
Из полученного результата вытекает как следствие, что точ ки с вертикальными касательными для кривой (6.88) соответству
ют экстремумам процесса а'й •. Из рис.6.5 видно, что угол по ворота радиуса-вектора щжвои (6.88) до точек экстремумов, а
следовательно, и моменты t и t3 г зависят от положения началь