Цифровая следящая система с конечным временем затухания свободного процесса Текст научной статьи по специальности «Математика»

УДК 62-50

DOI: 10.17586/0021-3454-2019-62-12-1078-1086

ЦИФРОВАЯ СЛЕДЯЩАЯ СИСТЕМА С КОНЕЧНЫМ ВРЕМЕНЕМ ЗАТУХАНИЯ СВОБОДНОГО ПРОЦЕССА

А. И. Коршунов

Военно-морской политехнический институт ВУНЦВМФ „Военно-морская академия им. Н. Г. Кузнецова",

198514, Санкт-Петербург, Россия E-mail: [email protected]

Показана возможность с помощью линейного дискретного корректирующего устройства получить в линеаризованной модели цифровой следящей системы конечное время затухания свободного процесса. Это время не превышает целое число периодов дискретизации, равное порядку непрерывной части системы, таким образом достигается предельно возможное быстродействие системы. Получены условия затухания свободного процесса за конечное время: управляемость и наблюдаемость дискретной модели непрерывной части системы. Для соблюдения условий необходимо отсутствие общих множителей числителя и знаменателя передаточной функции непрерывной части. Предложена методика выбора простейшей передаточной функции дискретного корректирующего устройства, учитывающая требования „грубости" системы. Рассмотрен пример цифровой следящей системы с непрерывной частью 3-го порядка.

Ключевые слова: цифровая следящая система, свободный процесс, конечное время затухания

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

Э. Джури называл апериодическим такой характер протекания переходных процессов [1]. Отечественные ученые, например В. П. Шипилло [2], также уделяли внимание дискретным системам с конечным временем затухания переходных процессов.

В настоящее время расчету цифровых регуляторов, обеспечивающих конечное время затухания свободных процессов, уделяется недостаточное внимание. Более того, даже в учебниках для вузов этот вопрос освещен недостаточно или даже некорректно, как например, в [3]. Это может создать у учащихся неправильное представление о проблеме. В работе [3] рассмотрен расчет цифрового регулятора, обеспечивающего в системе стабилизации частоты вращения электропривода постоянного тока достижение заданной скорости за один период дискретизации. Это значение достигается в конце первого периода дискретизации, однако свободный процесс при этом не только не заканчивается, он оказывается медленно затухающим. Поскольку в моменты квантования частота вращения принимает заданное значение, цифровой регулятор перестает реагировать. Вследствие этого переходный процесс протекает, как в разомкнутой системе, т.е. очень медленно.

В отечественной литературе расчет цифровых автоматических систем рассмотрен с приемлемой полнотой [4—6], но системы с конечным временем затухания свободных процессов рассмотрены недостаточно в работах, посвященных расчету ЦСС.

Передаточная функция дискретного корректирующего устройства, обеспечивающего конечное время затухания свободных процессов. Рассмотрим ЦСС, линеаризованная модель которой представлена на рис. 1.

ИИЭ

а

Дф

ФНП

НЧ

вх а

л

Т &[У7] й(г) (1-е-РТ)/Р »нч(р)

авых(0

Рис. 1

Идеальный импульсный элемент ИИЭ учитывает квантование сигналов по времени с периодом Т. Дискретный фильтр ДФ с передаточной функцией О(г) преобразует дискретные значения рассогласования &[\Т] в а^уГ] по корректирующей программе, реализуемой обычно в микропроцессоре. Фиксатор нулевого порядка ФНП учитывает факт запоминания $![\Т] на период Т в выходном регистре микропроцессора. ФНП представлен передаточной функцией формирующего элемента

1 - ерТ »фэ (Р) =-•

(1)

В передаточной функции непрерывной части НЧ »Нч (Р) учтено и преобразование а^ vГ ] в напряжение в микропроцессоре. Интегрирующие свойства исполнительного двигателя обеспечивают первый порядок астатизма ЦСС, что учтено множителем р в знаменателе передаточной функции:

кЯп-1( р)

»нч (Р) =

р<2п-1( Р)

(2)

где

Оы(Р) = (Т1Р + 1)(Т2Р +1)••• (Тп-1Р +1), Т >0, г = 1, 2,

Яп-1(р) = Ъп-1 Рп-1 +... + ЪР +1, Т > 0, Ьг > 0 г = 1, 2, . Передаточная функция разомкнутой ЦСС имеет вид [4—6]:

»Ц

., п -1, п -1.

где

НЧ(г) = 2-

"Цсс( г ) = 0( г )»НЧ ( г ),

Кп-1( Р) 1 = ^п-1( г)

Р2бп-1 (Р) I (г -1)(г - 4) ••• (г - ёп-1)

-Т/Т

е 1 г.

(3)

(4)

Яп-1 (г) = Рп-1гп -1 +... + Рг гг +... + р1г + Ро.

Коэффициенты вг определяются в результате разложения »Нч (Р)/ Р на простейшие дроби, их 2-преобразования и сложения результатов. Аналитические выражения для коэффициентов вг очень сложны и громоздки, численные значения вг можно получить, используя прикладные программы, например Ма1ЬаЬ.

Заметим, что дискретная передаточная функция НЧ помимо полюса г=1, лежащего на границе области устойчивости

|г| < 1, (5)

может иметь нули вне этой области даже при Яп-1(Р) = 1.

Для затухания свободного процесса в линеаризованной модели ЦСС за конечное время необходимо, чтобы существовало ограниченное управляющее (скалярное) воздействие, переводящее ее из любого исходного состояния в нулевое за конечное время [7]. В рассматриваемом в настоящей работе случае непрерывная часть управляема и наблюдаема при отсутствии общих множителей числителя и знаменателя ее передаточной функции [7].

Управляющее воздействие на непрерывную часть ЦСС является ступенчатой функцией, длительность ступенек равна шагу квантования по времени Т.

Для выявления существования ступенчатого управляющего воздействия, обеспечивающего конечное время полного затухания свободного процесса, то есть возвращения непрерывной части в нулевое состояние из любого исходного состояния, рассмотрим процесс в естественном фазовом пространстве НЧ [8]:

й 0

• = А0 + Ни, Л

Т

"вых = с 0,

(6)

где

0 =

" х1" " 1"

х2 , Н = Н2 , с = 0 , А =

_ хп _ _ Нп _ 0

0 1 0

0 0 0 0 «1 «2

1

ап-1

, Х1 "вых,

п-к -1

Нп-к = Ьк - X а]+кН], Н = Ьп-ъ к = п - 2, п - 3, 1, 0, } =1

I

а{, ЬI — коэффициенты передаточной функции Wнч (Р):

^НЧ(Р) = \+А1 +п^Ь.р + гь0 , а0 = 0. (7)

РП + апЛрп 1 +... + а1р

Обозначим и1 входной сигнал, действующий на НЧ на интервале 1Т<<(1+1)Т ,

0г- = 0(1Т), I = 0, 1, 2, 3, ... Проинтегрировав систему дифференциальных уравнений (6),

получим

01 = Н (Т )00 +1Н (Т - т)Ни0 Лт = Н (Т )00 + Ьи0,

0

Т

где Н (Т) = еАТ — матричная экспонента, Ь = | Н (Т - т)Нй т — вектор-столбец.

0

Аналогичным образом получим

Т

02 = Н (Т )01 +1Н (Т - т)Ни1йт = Н 2(Т )00 + Н (Т)Ьи0 + Ьи1,

(8)

0т = Н'"(Т )00 + Нт-1(Т )Ьи0 + Нш_ z (Т )Ьи1 +... + Н (Т )Ьит -2 + Ьит-1. Введем в рассмотрение вектор-столбец

и = [Uo, ul,..., ит-1]

и ( п х т )-матрицу

V = [Нт-1(Т)Ь Нт-2(Т)Ь|...|Н(Т)Ь|Ь]. Это позволяет записать выражение (9) в виде

гт-2.

VU = 0т - Нт (Т)0С

0т =0

= 00 00 =-Нт (Т)00.

(9)

(10) (11) (12)

Полагая 0т = 0, получим систему п уравнений с т неизвестными управляющими воздействиями и0,иь...,ит-1, переводящими непрерывную часть из состояния 00 в нулевое 0 = 0.

При т < п за т периодов управления можно перевести непрерывную часть в нулевое состояние 0т = 0 только из исходных состояний 00, принадлежащих области значений

матрицы (Нт (Т))-1V . Множество этих состояний является только частью п-мерного естественного фазового пространства непрерывной части, совпадающей с множеством всех линейных комбинаций т столбцов этой матрицы [9]. Оно имеет размерность т при линейной независимости столбцов матрицы.

При т = п и линейной независимости столбцов матрицы (Нт (Т))-1V возможен перевод непрерывной части из любого исходного состояния в нулевое. Следовательно, невырожденность квадратной ( п х п )-матрицы

12 ..... ' " (13)

(Нп(Т))-У = [Н-1(Т)Ь Н-2(Т)Ь

Н1-п (Т)Ь Н-п(Т)Ь]

является необходимым и достаточным условием существования и единственности управления (10), переводящего непрерывную часть за п периодов при любых начальных условиях в нулевое состояние.

Заметим, что для невырожденной матрицы Н (Т) полученное условие совпадает с условием управляемости линейной дискретной системы со скалярным управлением [10].

Из формул (12) и (13) при т = п получаем

-1

и = -[ Н-1(Т )ЬН -2(Т )Ь

Н1-п (Т )ЬН-п (Т )Ь] 00.

-п,

(14)

Поскольку ступенчатое управление переводит непрерывную часть ЦСС в нулевое состояние за конечное время, желательно определить передаточную функцию дискретного корректирующего устройства ^(2), решающего эту задачу в замкнутой системе. Наиболее просто определить ^(2), используя аппарат дискретных передаточных функций.

Обеспечить затухание свободного процесса за конечное время может передаточная функция замкнутой ЦСС вида

ф ( ) ^Цсс(2) ф цсс( г) = ■—--

ёпа2 1 +... + ё12 + £0

1 + Wц

(15)

"цсс(2) гп

позволяющая определить искомую передаточную функцию дискретного корректирующего устройства О(г)

Б( 2 ) =

^Цсс( 2) = ёп -1 гп-1 +... + £12 + ё0 Оп (г )

Wнч (2) 2п - £п-12п -1 - ... - £12 - £0 кЯт (2У

(16)

где

^сс(2) =

фЦСС(2 )

1 - фт

ёп-12 '

п-1

ё12 + ё0

'ЦСС(2) 2П-ёп-12п -...- ё12-ё0 (17)

Оп ( 2 ) = ( 2-1)( 2-й1)...( 2-йп-1).

Как известно [5], для обеспечения „грубости" замкнутой скорректированной системы передаточная функция разомкнутой системы (^Цсс( 2) (3)) должна содержать все нули и полюсы дискретной передаточной функции непрерывной части (Wнч (2) (4)), равные и большие по модулю 1, т.е. лежащие на границе области устойчивости (5) и вне нее.

Из полюсов ^Нч (2) необходим 2=1 не только по соображениям „грубости", но и для сохранения первого порядка астатизма системы. Следствием этого является необходимость наличия у Wцсс( 2) (17) полюса 2=1. Это требует выполнения условия

" 1 £12-ё0 2=1 = 0 или ёп-1 +... + ё1 + ё0 = 1 (18)

2п-ёп-12^1

Чтобы сохранить нули ^нч (2), лежащие вне области устойчивости (5), целесообразно оставить все нули, т.е. корни ^п-1(2). Такой подход позволяет упростить дискретное корректирующее устройство. Для этого достаточно, чтобы числитель желательной передаточной функции замкнутой ЦСС Ф цсс( 2) (15) отличался от числителя ^нч (7) (4) только постоянным множителем. Выбрав в качестве этого множителя ( к (Рп_1 +... + Р1 + Ро))1, получим

( и_1 V1 (п_1 V1 8 = кРг к 2 рг = рг 2 вг . (19)

V г =0 У V г=0 У

Использование формулы (19) гарантирует выполнение условия (18). Таким образом, получаем

Оп (7) = (7 _ 1)(2 _^1)...(Г _ П

Я( 7 ) = ■

(п_1 Л

( п_1

V

2 Р/ (( _8п_12>7 1 _ .. _812 _8о) к 2 Рг'^ _Рпч^ 1 _ ... _Р^ _Рс

V/=о у V/=о

Сократив общий множитель (г-1) в числителе и знаменателе полученного выражения, получим передаточную функцию дискретного корректирующего устройства

( г_^1)...( г_ёп _1)

Б( 2 ) = ■

к(4п_1 412_4о)

где 4п_/ = Рп_/ + .. + Р1 + Ро г = 1, 2, . ., п.

Таким образом, порядок О(г) равен п-1, т.е. на единицу меньше порядка ^нч(2). Подстановка передаточной функции П(г) (2о) в выражения передаточных функций

(2о)

" цсс (7) (3) и Ф цсс (7) дает соответственно

Р п_ ^1 + ... + Р1Г + Ро

Wцcc(7)=

(Рп_1 + ... + Р1 + Ро)^ _Рп_1*Г1 _... _Р1^ _Р

Рп- ^ 1+ + М + Ро (1 _1)(4п_ 1 +... + 417 + Чо)

.п _1

»цсс(2) = вп_1!п_1 +... + Р^ + Ро = ёп^!"1 +... + 81! + Во

Ф

цсс

(7) = ■

(21)

(22)

1 _ »ЦСС (г) (Рп_1 +... + Р1 + Ро) ^ 7п

Переходную характеристику ЦСС ЩуТ], необходимую для оценки качества управления в переходных режимах, получим с помощью обратного Z-преобразования произведения

2 - ._! , , „ 2_(п_1) , „ __п> 7

Ф ЦСС( 2 )--7 = (ёп _12 +... + В!

+ ёо ')"

h[vT ] = ё п _1 1[( V _1)Т ] +... + 81 -1[(у_ п+1)Т ] + ёо-1[(У_ п)Т ]. Выражение (24) можно записать в более удобном для восприятия виде

о, V = о,

(23)

(24)

(25)

Ь^Т]= «¡28п-г, V = 1, 2, ..., п_1, /=1

1, V > п.

Свободная составляющая выходной величины &выхсв (^) в моменты квантования /=vT при $вх [vT ] = 1[уТ ] принимает значения

[уТ ]=

1, V = 0,

V

1 -X ёп-1, V = 1 2, ..., п - 1, г=1

0, V > п.

(26)

Можно показать, что и Звых св ^) = 0 при t > пТ в случае наблюдаемости импульсной системы [10]. Достаточно установить равенство нулю свободного процесса при t = пТ, т.е. 0св (пТ) = 0.

При 0сВ(пТ) * 0, получаем 0св((п + 1)Т) = Н(Т)0св(пТ), 0сВ((п + 2)Т) = Н2(Т)0сВ((п + 1)Т), ..., 0св((2п - 1)Т) = Нп-1(Т)0св((2п - 2)Т), откуда ^ххв^Т] = сТ0св[пТ], Звыхсв [(п + 1)Т] = = сТ 0св [(п + 1)Т ],..., Звыхсв (2п + 1)Т ] = сТ 0св[(2п + 1)Т].

Поскольку Звых [vT] = 0 при V > п, для определения 0^)| =пТ получаем систему п линейных уравнений, которую можно записать в матричной форме

" сТ ' " 0"

сТН (Т) 0[пТ ] = 0

_ сТНп-1(Т) _ 0

или

сНТ (Т)

с( НТ (Т ))п-1

0[пТ ] = 0.

(27)

Невырожденность матрицы наблюдаемости

сТН (Т )

сТНп-1(Т )

[7] доказывает

равенство 0св (t) | г=пТ = 0 и, следовательно, полное затухание свободного процесса в непрерывной части за п периодов дискретности при любых начальный условиях.

Возможно полное затухание свободного процесса в непрерывной части за т<п периодов дискретности при начальный условиях, принадлежащих области значений (п х т )-матрицы

[Н-1(Т)Ь Н-2(Т)Ь

Н1-п (Т )ЪН~Г 1(Т )Ь].

-п,

Пример. Рассмотрим ЦСС с Wнч(р) = к /[р(Т1 р + 1)(Т2р +1)], при Т1=0,1 с, Т2=0,02 с, Т=0,01 с, к=10 с-1. По выражениям, полученным в [11], рассчитаем параметры дискретной передаточной функции непрерывной части

^нч (2) = к в222 + в12 + в0 , (2 -1)(2 - йх)(2 - й2)

где в2 = 7,2024 -10-5, в1 = 2,4904 -10-4

= 0,90484, й2 = 0,60653 . Заме-

в0 = 5,3367-10"

тим, что ^нч(2) имеет нули 21=-3,2282, 22=-0,2295, первый из которых лежит вне области

устойчивости (5). Это требует соблюдения условия „грубости". Вычисление по формуле (20) передаточной функции дискретного фильтра Б( 2) дает

Б( 2) = .

2

(2 + «12 + «0)

к (Ч222 + Ч12 + %)

где а1=-^1-^2=-1,5114, 00=^2=0,54881, д2=в2 + в1 + во = 3,7444-10

-3

(28)

91=в1 + в0 =

= 3,0241-10-3, д0=в0 = 5,3367 10-4

Т

На рис. 2 представлена математическая модель ЦСС с конечным временем затухания свободного процесса, построенная в системе Ма1:ЬаЬ.

Тглги^р Реп

Рис. 2

Переходная характеристика ЦСС И(г), дискретные значения рассогласования $[vT] и его преобразованные по корректирующей программе значения vT ] представлены на рис. 3, откуда видно затухание свободного процесса за три периода дискретности (3Т).

е 0,8 о,б о,4 о,2 о

-о,2 -о,4 -о,б -о,8

о 1 2 3 г/Т

Рис. 3

Процессы в ЦСС при линейно возрастающем воздействии с единичной скоростью нарастания представлены на рис. 4. Поскольку начальное значение линейно возрастающего воздействия равно нулю, первый период система, находящаяся при нулевых начальных условиях, не работает. Свободный процесс при этом затухает за два периода дискретности (2Т).

9

о,о3 о,о2 о,о1 о

-о,о1 -о,о2

о о,о1 о,о2 о,о3 г, с

Рис. 4

Если линейно возрастающее воздействие имеет небольшое начальное значение, процесс его отработки начинается в момент г=о, а свободный процесс затухает за три периода дискретности (3Т). При увеличении начального линейно возрастающего воздействия до значения, соответствующего приращению за период дискретности, процесс его отработки начинается при г=о, а свободный процесс затухает за два периода дискретности (2Т) (рис. 5).

е

0,04 0,03 0,02 0,01 0

-0,01 -0,02

0 0,01 0,02 0,03 t, c

Рис. 5

Выводы

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

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

3. Предложена методика выбора простейшей передаточной функции дискретного корректирующего устройства, учитывающая требования „грубости" системы.

СПИСОК ЛИТЕРАТУРЫ

1. Джури Э. Импульсные системы автоматического регулирования. М.: Физматгиз, 1963. 567 с.

2. Шипилло В. П. Операторно-рекуррентный метод расчета электрических цепей и систем. М.: Энергоатомиздат, 1991. 311 с.

3. Анучин А. С. Системы управления электроприводов. М.: Изд. дом МЭИ, 2015. 373 с.

4. Бесекерский В. А., Попов Е. П. Теория систем автоматического регулирования. М.: Наука, 1972. 768 с.

5. Бесекерский В. А. Цифровые автоматические системы. М.: Наука, 1976. 575 с.

6. Бесекерский В. А., Ефимов Н. Б., Зиатдинов С. И. и др. Микропроцессорные системы автоматического управления. Л.: Машиностроение, 1988. 365 с.

7. Заде Л., Дезоер Ч. Теория линейных систем. Метод пространства состояний. М.: Наука, 1970. 704 с.

8. Коршунов А. И. Определение коэффициентов векторно-матричного дифференциального уравнения линейной непрерывной части системы со статическим преобразователем // Изв. вузов. Приборостроение. 2016. Т. 59, № 2. С. 107—112.

9. Ланкастер П. Теория матриц. М.: Наука, 1969. 280 с.

10. Изерман Р. Цифровые системы управления. М.: Мир, 1984. 541 с.

11. Коршунов А. И. Основы теории управления. Руководство к курсовому проектированию. Петродворец: ВМУРЭ им. А. С. Попова, 1998. 136 с.

Сведения об авторе

Анатолий Иванович Коршунов — д-р техн. наук, профессор; Военно-морской политехнический институт ВУНЦ ВМФ „Военно-морская академия им. Н. Г. Кузнецова", кафедра радиоэлектроники; E-mail: [email protected]

Поступила в редакцию 16.09.19 г.

9[v7] 0 «(t)

91[v7]/500

Ссылка для цитирования: Коршунов А. И. Цифровая следящая система с конечным временем затухания свободного процесса // Изв. вузов. Приборостроение. 2о19. Т. 62, № 12. С. Ш78—Ш86.

DIGITAL TRACKING SYSTEM WITH FINITE DECAY TIME OF FREE PROCESS

A. I. Korshunov

Naval Polytechnic Institute "N. G. Kuznetsov Naval Academy", 198514, St. Petersburg, Russia

E-mail: [email protected]

The possibility to obtain a finite decay time of free process in a linearized model of a digital tracking system using a linear discrete correction device is demonstrated. The decay time does not exceed the integer number of sampling periods equal to the order of the system continuous part, and therefore the maximum possible speed of the system is achieved. The conditions of free process damping for finite time are derived: controllability and observability of the discrete model of the continuous part of the system. To comply with the conditions, the absence of common multipliers of the numerator and denominator of the transfer function of the continuous part is necessary. The method of choice of the simplest transfer function of the discrete correcting device considering requirements of system roughness is proposed. An example of a digital tracking system with a continuous part of the 3rd order is analyzed.

Keywords: digital tracking system, free process, finite decay time

REFERENCES

1. Jury E.J. Sampled-data control systems, NY, Wiley, London, Chapmen and Hall, 1958.

2. Shipillo V.P. Operatorno-rekurrentnyy metod rascheta elektricheskikh tsepey i sistem (Operator-Recurrence Method for Calculating Electrical Circuits and Systems), Moscow, 1991, 311 p. (in Russ.)

3. Anuchin A.S. Sistemy upravleniya elektroprivodov (Electric Drive Control Systems), Moscow, 2015, 373 p. (in Russ.)

4. Besekerskiy V.A., Popov E.P. Teoriya sistem avtomaticheskogo regulirovaniya (The Theory of Automatic Control Systems), Moscow, 1972, 768 p. (in Russ.)

5. Besekerskiy V.A. Tsifrovye avtomaticheskie sistemy (Digital Automatic Systems), 1976, 576 p. (in Russ.)

6. Besekerskiy V. A., Efimov N. B., Ziatdinov S. I. et al. Mikroprotsessornyye sistemy avtomaticheskogo upravleniya (Microprocessor-Based Automatic Control Systems), Leningrad, 1988, 365 p. (in Russ.)

7. Zadeh L., Desoer Ch. Linear system theory: the state space approach, 1963.

8. Korshunov A.I. Journal of Instrument Engineering, 2016, no. 2(59), pp. 107-112. (in Russ.)

9. Lankaster P. Theory of matrices. NY, London, 1969.

10. Isermann R. Digital control systems, Berlin etc., 1981.

11. Korshunov A.I. Osnovy teorii upravleniya. Osnovy teorii i sistem avtomaticheskogo upravleniya (Bases of the Theory of Management. Bases of the Theory and Systems of Automatic Control), Petrodvorets, 2017, pp. 136. (in Russ.)

Data on author

Anatoly I. Korshunov — Dr. Sci., Professor; Naval Polytechnic Institute "N. G. Kuznetsov Naval

Academy", Department of Radio Electronics; E-mail: [email protected]

For citation: Korshunov A. I. Digital tracking system with finite decay time of free process. Journal of Instrument Engineering. 2019. Vol. 62, N 12. P. 1078—1086 (in Russian).

DOI: 10.17586/0021-3454-2019-62-12-1078-1086