УДК 62-50
DOI: 10.17586/0021-3454-2018-61-12-1036-1044
ПОВЫШЕНИЕ ДИНАМИЧЕСКОЙ ТОЧНОСТИ ЦИФРОВЫХ СЛЕДЯЩИХ СИСТЕМ АСУ ТП МЕТОДОМ КОМБИНИРОВАННОГО УПРАВЛЕНИЯ. Ч. II. ВЫСОКИЙ ТЕМП ВЫЧИСЛЕНИЯ РАССОГЛАСОВАНИЯ
А. И. Ермоленко, А. И. Коршунов
Военно-морской политехнический институт ВУНЦВМФ „Военно-морская академия им. Н. Г. Кузнецова",
198514, Санкт-Петербург, Россия E-mail: [email protected]
Повышение производительности автоматизированных систем управления технологическим процессом (АСУ ТП) ужесточает требования к точности цифровых следящих систем (ЦСС) опорно-поворотных устройств. Компенсации скоростной ошибки и ошибки по ускорению часто достаточно для удовлетворения требованиям динамической точности. Рассмотрена компенсация скоростной ошибки и ошибки по ускорению ЦСС современных АСУ ТП с достаточно высоким темпом вычисления рассогласования. При увеличении темпа вычисления рассогласования ЦСС комбинированного управления по свойствам приближается к непрерывной следящей системе комбинированного управления, параметры которой определены. Предложен метод расчета ЦСС комбинированного управления. Суть метода состоит в обеспечении в непрерывной следящей системе показателей качества управления, несколько лучших, чем заданы для ЦСС, и выборе периода вычисления рассогласования, обеспечивающего заданное качество управления ЦСС. Такой подход позволяет получать в ЦСС комбинированного управления высокое качество как свободных, так и переходных процессов, вызванных изменением задающего воздействия.
Ключевые слова: компенсация, скоростная ошибка, ошибка по ускорению, высокий темп вычисления рассогласования
Повышение производительности технологических процессов обусловливает ужесточение требований к управляющим ими автоматизированным системам управления (АСУ ТП) и, в частности, к электромеханическим цифровым следящим системам (ЦСС), осуществляющим механическое управляющее воздействие на объект управления. При повышении динамичности изменения управляющих воздействий, оцениваемой обычно их первой и второй производными по времени, возрастают требования к динамической точности ЦСС, оцениваемой скоростной ошибкой и ошибкой по ускорению.
Использование радикального средства повышения динамической точности ЦСС — исключения скоростной ошибки и ошибки по ускорению за счет повышения порядка астатизма ЦСС до третьего — делает следящую систему „условно устойчивой" [1], поскольку уменьшение коэффициента преобразования, например вследствие насыщения усилителя, может вызвать ее неустойчивость. Поэтому более целесообразна компенсация указанных ошибок за счет использования управления по задающему воздействию. При этом дискретный фильтр (ДФ), выполняющий его дополнительную обработку, имеет невысокий порядок и только нулевые полюсы передаточной функции [1, 2].
В ЦСС проектируемых АСУ ТП период вычисления рассогласования может быть выбран намного меньшим постоянных времени непрерывной части ЦСС. Это позволяет приблизить ЦСС по свойствам к непрерывным следящим системам, что дает им определенные преимущества. Задача настоящей статьи — рассмотрение особенностей расчета ЦСС комбинированного управления с высокой частотой вычисления рассогласования.
Будем полагать, что вследствие охвата элементов системы отрицательными обратными связями с постоянными параметрами можно считать, что передаточная функция непрерывной части (НЧ) ЦСС имеет вид:
Жнч (р) =-К-, (1)
р(тр + 1)(Г2р +1)' ^
где коэффициент преобразования К и эквивалентные постоянные времени Т1 , Т2 обладают достаточным постоянством. Увеличение скоростной ошибки, вследствие стабилизации К за счет обратной отрицательной связи по скорости, может быть компенсировано комбинированным управлением.
Дискретную передаточную функцию разомкнутой ЦСС можно представить следующим образом [3]:
' 1 1 .Яо(*)
Ж (7,0) = 7Жнч (р)[ = К—2<
р2(Т р + 1)(Т2 р +1) ] Qз( *У
(2)
где
= в272 + + р0, 63 = (^ - 1)(7 - ^1)(Г -= + ао^2 + а^ + СЦ), ^ = е"Т/Т1,
12 = е-Т/То, а2 = -(1 +1 + 12), а1 = 11 + 12 + а0 = -1^, 1 + а2 + а1 + а0 = 0, в2/ К = Т+(Т1 + 72х 1 + ^ -1)+т211 -112 62,
Т - Т2
в1/К = -(11 +12)(Т + Т1 +Т2)-(Т1 +Т2)(йхй2 -1)-2 Т1 Т11
Т1 - Т2
Р0/ К = (Т+Т1 + т2)йхй2 +
То2!. - Т1212 Т1 - Т2
При Т / Т << 1, г = 1, 2 с приемлемой точностью можно разложить коэффициенты знаменателя аг (г = 0,1,2) по степеням Т, ограничиваясь членами степенного ряда второго порядка малости:
а2 =
-3 + Т -1Т12 + Т22 Т2
Тэ 2 Т12Т22
а1 =
Т1 3 - 2— + -Тэ 2
^ 2
К + ТТ + _1 V Т2Т22 тэ2 ,
[1+0(Т)], Т" [1 + 0(Т)],
(3)
а=
Т - 1Т_
- + Тэ - 2 Тэ2
[1+0(Т)],
ээ
и разложениями коэффициентов числителя рг-, ограничиваясь членом разложения в степенной ряд третьего порядка малости:
во = Р0 = 6ТКТТ3[1 + 0(Т)], в1 = 4в0 = -ттт3[1 + 0(Т)], 6т1т2 3т1т2
(4)
где 0(Т) — величина порядка малости, большего Т, Тэ = (Т1Т 2)/(Т1+ Т 2).
Поскольку из теории [1] известно, что для компенсации скоростной ошибки и ошибки по ускорению ЦСС в закон управления при первом порядке астатизма системы необходимо вводить первую и вторую разности задающего воздействия, заранее не известного, передаточная функция ДФ должна иметь вид:
ф(z) = q2 + q1 z 1 + q0z 2 = (q2z2 + q1z + q0)/ z2, qt = const, i = 0,1,2. (5)
Согласно структурно-динамической схеме ЦСС комбинированного управления (рис. 1, а) получаем передаточную функцию замкнутой ЦСС комбинированного управления:
Фк( z, 0) =
(1+ф( z) )w (z, 0) = [(1+q2) z 2 + q\z+40 ] (p 2z 2 + p1z+p 0)
1 + W (z, 0)
z5 + (a 2 + p 2) z4 + (a1 + p1) z3 + (a0 + p 0) z2
(6)
и передаточную функцию эквивалентной системы с управлением по отклонению (рис. 1, б) в разомкнутом состоянии [4]:
W3 (z, 0) =
Фк (z,0) =(1 + ф( z) )W (z,0)
1- Фк (z, 0) 1- ф( z )W (z,0)
(7)
Рис. 1
Подстановка в формулу (7) выражений (2) и (6) с учетом необходимого для компенсации ошибок 3-го порядка астатизма Жэ(г,0) дает:
W3 (z, 0) =
Rt( z) = R4( z)
Q5(z) (z-1)3 5*2 (z)
(8)
где
Я4 (г) = [ г2 + г 2ф(г Щ (г )= [(1 + г 2 + 4 г + д0 ] (р 2г2 + ^ + р 0), 65 (г) = г 203 (г) --г2ф№2(2)=25 + (а2 -Р242)г4 + (а1 -Р241>3 + (а0 -Р042 -Р240>2 -
-(Р 041 + Р1^0) г-р 040 = (г-1)3 £2(г), г) = г2 + + %
Поскольку полиномы Q5(z), dQ5(z)/dz, С^5 (г)/Сг2 при подстановке г=1 обращаются в нуль, получим систему линейных уравнений для вычисления коэффициентов передаточной функции ДФ ф(г) :
42 + 41 + 40 = 0,
(4Р 2 + 3Р1 + 2Р0)42 + (3Р2 + 2Р1 + Р0)41 + (2Р2 + Рх)40 = 3 + а2 + а1, (9)
(6Р2 + 3Р1 + Р0)42 + (3Р2 + Р1)41 + в240 = 9 + 5а2 + 2а
(отсутствие в правых частях системы параметра а0 объясняется тем, что 1 + а 2 + а1 + а0 = 0), откуда:
42 = [3в 2 + 6Р1 + 9Р0+а2(Р2 + 3в1 + 5Р0) + а1(в1 + 2Р0>]/(Р 2 + Р1 + Р0)2,
41 = [-3Р2 -9Р1 -15Р0 + а2( -4Р1 -8Р0) + а1(Р2 -Р1 -3Р0)]/(Р2 + Р1 + Р0)2, (10)
40 = [3Р1 + 6Р0+а 2 (-Р 2 + Р1 + 3Р0) + а1(-Р 2 + Р 0)]/(Р2 + Р1 + Р0)2,
где Р2 + Р1 + Р0 = КТ(1- С1)(1 - С2) [3].
Подстановка в полученные выражения приближенных значений параметров а;-, Р/ (/=0, 1, 2) дает приближенные значения для 41 (/=0, 1, 2):
42 = 40 = ^2 [1+ 0(T)], 4 ="2^0 =-2 [1 + 0(T)]. (11)
KT KT
3
Деление полинома Qs(z) на (г-1) позволяет определить коэффициенты полинома £2( г) = г + s1z + 50 :
¿1 = 3 + а2 - 42Р2, ¿0 = 6 + 3а2 + а1 - 3Р242 - Р4 - Р241. (12)
В результате подстановки выражений (3) и (4) получим приближенные значения
¿1 = 55[1 + 0(Т)], ¿0 = 5[1 + 0(Т)], 5 = Т /(6ТЭ). (13)
2
Преобразование передаточной функции г ф(г):
г2ф(г)=4ог2 + 41г + 40 = 42(г -1)2 + (41 - 242)г + % -42 с учетом легко доказываемого равенства:
41 -24о =-(40 -42) =1/КТ
2
позволяет разложить г ф(г) по степеням (г-1):
г2ф(г) = 4о(г -1)2 + -17(г -1), (14)
Кт
что соответствует выделению во входном сигнале замкнутого контура ЦСС 0 вх1[п— первой и второй разности задающего воздействия:
0вх1[пТ] = 0вх [пТ] + К- А0вх [(п - 2)Т] + 42А\х [(И - 2)Т], (15)
Кт
где
А0вх [(п - 2)Т ] = 0Вх [пТ ] - 0Вх [(п - 1)Т ], А20вх [(п - 2)Т ] = 0вх [пТ ] - 20вх [(п - 1)Т ] + 0вх [(п - 2)Т ]. Подстановка в формулу (8) выражений (3), (4), (11), (13), (14) дает:
W (z,0) =
z2 + ^^ [1 + 0(T)](z -1)2 + -1; (z -1) KT2 KT
пЪ
KT3
[(1 + 0(T)) z 2 + 4(1 + 0(T) z +1 + 0(T)]
—Л-m- <16>
(z -1)3 [ z2 + 55(1 + 0(T)) z + 5(1 + 0(T))]
Частотная передаточная функция эквивалентной разомкнутой системы получается в результате подстановки в (8) [1]:
г = 1 + Ят/2 , г = Ркт , (17)
1 - рК Т/ 2 1 - рК Т/ 2
где
л 2 roT
X = — tg--(18)
T2
— абсолютная псевдочастота [1, 2], го — реальная круговая частота.
Перейдя к пределу в W3 (jX) при Т^0, с учетом Tim X = го, получим:
lim W3(jX,0) = W3(jro) = (T + T2Xjro)2 +3jro+K . (19)
TTMjro)3 v '
Нетрудно проверить, что Жэ (ую) — частотная передаточная функция разомкнутой непрерывной следящей системы эквивалентной следящей системе комбинированного управления, структурно-динамическая схема которой представлена на рис. 2, а.
Рис. 2
Таким образом, при достаточной малом значении Т свойства ЦСС комбинированного управления (рис. 1, а) близки к свойствам непрерывной следящей системы комбинированного управления (рис. 2, б).
Аналитическая оценка достаточной малости Т представляет сложную задачу и лишена практического смысла, поскольку проще оценить достаточную малость Т путем цифрового моделирования, например, в системе Ма1ЬаЬ 81шиНпк.
Расчет ЦСС комбинированного управления предполагает решение двух задач:
1) обеспечение быстрого затухания свободного процесса в замкнутом контуре управления,
2) обеспечение допустимой колебательности переходных процессов при комбинированном управлении.
Варьируемыми параметрами являются коэффициент преобразования К и шаг квантования по времени Т.
При решении первой задачи необходимо подбором К обеспечить с небольшим запасом в непрерывной модели ЦСС заданные показатели качества управления в переходных режимах (затухание свободного процесса).
Для второй задачи следует определить наибольшее значение Т, при котором ЦСС комбинированного управления удовлетворяет заданным требованиям. Если это значение Т оказывается слишком малым для реализации, можно увеличить запас в показателях качества непрерывной модели.
Пример. Рассчитаем ЦСС комбинированного управления с компенсацией скоростной ошибки и ошибки по ускорению при требованиях к качеству переходных процессов в замкнутом контуре управления ЦСС, заданных допустимым перерегулированием одоп = 30 % и
допустимым временем переходного процесса ^пдоп = 0,3 с. Для ЦСС комбинированного
управления одоп = 35 % , ¿пдоп = 0,3 с.
Для непрерывной модели замкнутого контура управления критический коэффициент преобразования [5]:
Ккр = — + — = +-
Т1 Т2 0,02 0,01
1 = 150 с"1.
Приняв средний запас устойчивости по амплитуде АЬ = 14 дБ, получаем значение К = Ккр/5=30 с1. Поскольку первая частота сопряжения логарифмической АХЧ больше К, можно принять частоту среза юс=К и определить запас устойчивости по фазе:
у=90° " агС^К^)" агС£(КТ2) = 90° - агС^, 6) - агС^,3) = 42,3°, тоже близкий к среднему значению.
Определив резонансную частоту замкнутого контура управления, с учетом (Т + Т2) - 1 = =0,8 > 0 по формуле:
ф|2 + Т22 + Т22)2 + 3/,2/22(2 К (Т + Т2)]/3 )Т1-1Т2-1 = 27,11 с-1
Юр =1
вычислим показатель колебательности замкнутого контура непрерывной модели:
М =
Щнч (УЧ)
1 + Щч С/Юр)
К
К - (Т + Г2)юР + /Т1Г2Юр(1 - ю2)
= 1,2268.
р .^^^р^ "р^
Полученные значения юс и М позволяют удовлетворить заданным требованиям к качеству переходных процессов в замкнутом контуре.
Результат моделирования переходной характеристики непрерывной модели замкнутого контура (рис. 3, кривая 1) позволяет определить 0,16 с, а = 20 %.
Ь, у.е.
1
0,8 0,6 0,4 0,2
4
..........Г ' \ 2 1
4~
/
7
0
А, у.е. 1,4
1,2
1,0
0,8
0,6
0,4
0,2
0
0,05
0,1 Рис. 3
0,15
0,2 г„ с
Л 2
Г \\
}
50
200
250 ю, с-
100 150 Рис. 4
Согласно амплитудно-частотной характеристике непрерывной модели ЦСС комбинированного управления (рис. 4, кривая 1) получим значение ее резонансной частоты юрк=41,242 с-1 и показателя колебательности Мк=1,3461. Переходная характеристика (рис. 3, кривая 2) дает лучшее значение (а=19,5 % и ¿п=0,125 с), чем у замкнутого контура, тогда как при компенсации только скоростной ошибки а=45,1 %, ¿п=0,2067 с. Объяснить это можно следующим. Компенсация (теоретическая) ошибки, пропорциональной третьей производной, позволяет получить передаточную функцию непрерывной модели ЦСС комбинированного управления в виде:
К
Фк( Р) = (1 + ф( Р))-
ТТ Р3 + (Т+Т2) Р2 + Р+к
= 1,
ф( р нтт2 р3 +(т1+т2) р2 + Р )/к
что обеспечивает а=0, ¿п=0.
Решая вторую часть задачи, примем Т= ТУ10=2-10-3 с. Вычисление параметров дискретной передаточной функции ДФ ф( г) по точным формулам (2) и (10) с пятью значащими цифрами дает: ^2=283,33, ^=-550,00, ^2=266,67. На математических моделях исходной ЦСС и ЦСС комбинированного управления, построенных в системе Ма1ЬаЬ 6.5 81шиНпк 5 (рис. 5), получены их переходные характеристики (см. рис. 3, кривые 3 и 4 соответственно). Для
исходной ЦСС получены а=21,84 %, ¿п=0,207 с, это несколько хуже, чем у непрерывной модели, но вполне удовлетворяет заданным требованиям.
С идо!*
Рис. 5
Компенсация скоростной ошибки и ошибки по ускорению в ЦСС комбинированного управления привела к изменению показателей качества: а=28,45 %, ¿п=0,145 с. Полученные показатели качества ЦСС комбинированного управления также вполне удовлетворяют заданным требованиям, но хуже, чем у ее непрерывной модели.
Амплитудно-фазовая частотная характеристика замкнутой ЦСС комбинированного управления, зависящая от реальной частоты (см. рис. 4, кривая 2), дает значения =41,154 си М к =1,5224, что больше, чем в ее непрерывной модели, это объясняет увеличение а.
Представляет практический интерес оценка влияния уменьшения шага квантования по времени Т на показатели качества управления ЦСС (а и tп) в обоих случаях. Уменьшив Т
вдвое (Т = 10-3 с), получим: д2=1066,7; ^=-2100,0; д2=1033,3. Изменив в математической модели (см. рис. 5) параметры ДФ и шаг квантования по времени Т, получим в результате моделирования переходные характеристики исходной ЦСС и ЦСС комбинированного управления (рис. 6, кривые 3 и 4). Их анализ показал соответственно а=21 %, ^=0,157 с; а=24,1 %, ^=0,139 с. Очевидна большая близость переходных характеристик ЦСС к характеристикам их непрерывных моделей (кривые 1 и 2 на рис. 3 и 6).
Рис. 6
Точность ЦСС комбинированного управления при Т=2 10-3 с оценена путем моделирования отработки гармонического воздействия 9вх (^ = 9вх т зт(ю^, 9вх [пТ] = 9вх т зт[юпТ]
при ю=2 с т = ю9
-1
о — 25° °вх т ■
Максимальная скорость изменения воздействия
2 2
— 50° / с, а максимальное ускорение ввх т = ю 9вх т —100° / с . Результаты отра-
вх т вх т ' ~ > " ..^..^тшу ^ ^^^^^ ^х т
ботки гармонического воздействия представлены на рис. 7, где 9вх — задающее воздействие, 91 и 92 — рассогласование ЦСС комбинированного управления и исходной ЦСС. Даже при
100-кратном увеличении из рис. 7 невозможно определить точное значение 91 т — 1,3°-10 .
Сравнение с амплитудой рассогласования исходной ЦСС 92 т —1,67 ° показывает повышение
в 1300 раз точности за счет комбинированного управления. Для сравнения рассчитана амплитуда рассогласования непрерывной модели комбинированной ЦСС:
9 т —
9в
1 + Щэ (/ю)
9в
25°
Щэ (/ю)| 18717
—1,3 °-10
—з
Щ (/ю) —
(Т + /2 )(/ю)2 + /ю+К Т1Т2С/Ю)3
хорошо согласующаяся с определенной выше ошибкой ЦСС комбинированного управления.
8, у.е. 30
20
10
0
-10
-20
-30
10081 0вх
1062
1-
I
0,5
1,5
Рис. 7
2,5
3 /, с
Ошибка слежения определена по линеаризованной модели, не учитывающей квантование по уровню и погрешности в определении параметров непрерывной части ЦСС и в установке параметров ДФ, вызывающих снижение точности.
Таким образом, по работе можно сделать следующие выводы:
1) компенсация скоростной ошибки и ошибки по ускорению в ЦСС комбинированного управления позволяет существенно повысить ее динамическую точность;
2) выбрав не слишком малый шаг квантования по времени (темпа вычисления рассогласования), возможно получить достаточное качество управления в переходных режимах как в замкнутом контуре, так и в ЦСС комбинированного управления.
список литературы
1. Бесекерский В. А. Цифровые автоматические системы. М.: Наука, 1976. 576 с.
2. Бесекерский В. А., Попов Е. П. Теория систем автоматического регулирования. М.: Наука, 1972. 768 с.
3. Коршунов А. И. Основы теории управления. Руководство к курсовому проектированию. Петродворец: ВМУРЭ им. А. С. Попова, 1998. 136 с.
4. Бесекерский В. А., Федоров С. М. Применение эквивалентной передаточной функции при расчете следящих систем комбинированного управления // Тр. I Междунар. конгр. Междунар. федерации по автоматическому управлению. М.: Изд-во АН СССР, 1961. Т. I. С. 154—165.
5. Коршунов А. И. Основы теории управления. Ч. II. Основы теории и систем автоматического управления. Петродворец: ВМПИ, 2017. 167 с.
0
1
2
Артем Игоревич Ермоленко Анатолий Иванович Коршунов
Сведения об авторах курсант; Военно-морской политехнический институт ВУНЦ ВМФ „Военно-морская академия им. Н. Г. Кузнецова", факультет САУ д-р техн. наук, профессор; Военно-морской политехнический институт ВУНЦ ВМФ „Военно-морская академия им. Н. Г. Кузнецова", кафедра радиоэлектроники; E-mail: [email protected]
Поступила в редакцию 10.10.17 г.
Ссылка для цитирования: Ермоленко А. И., Коршунов А. И. Повышение динамической точности цифровых следящих систем АСУ ТП методом комбинированного управления. Ч. II. Высокий темп вычисления рассогласования // Изв. вузов. Приборостроение. 2018. Т. 61, № 12. С. 1036—1044.
INCREASE IN DYNAMIC ACCURACY OF DIGITAL TRACKING SYSTEMS OF APCS BY COMBINED MANAGEMENT METHOD. PART II. THE HIGH RATE OF MISMATCH CALCULATION
A. I. Ermolenko, A. I. Korshunov
Naval Polytechnic Institute "N. G. Kuznetsov Naval Academy", 198514, St. Petersburg, Russia E-mail: [email protected]
Increase in productivity of industrial control system raises requirements to the accuracy of the digital tracking systems (DWS) of basic and rotary devices. Compensation of a high-speed mistake and mistake on acceleration often appears to requirements of dynamic accuracy is often sufficient for meeting the requirements. In part II, compensation of a high-speed mistake and mistake of acceleration of TsSS modern industrial control system in which rate of calculation of a mismatch can be chosen rather high is considered. At increase in rate of calculation of a mismatch of TsSS of the combined management on properties comes nearer to the continuous watching system of the combined management which parameters are determined. The method of calculation of TsSS of the combined management is proposed. The essence of the method consists in providing the continuous tracking system with management quality indicators which are slightly better than the set in TsSS, and the choice of mismatch calculation period providing the prescribed quality of TsSS management. As a result, it is possible to receive in TsSS of the combined management, a high quality of the free processes as well as of transition processes caused by change of the setting influence.
Keywords: compensation, high-speed mistake, acceleration mistake, high rate of mismatch calculation
REFERENCES
1. Besekerskiy V.A. Tsifrovye avtomaticheskie sistemy (Digital Automatic Systems), 1976, 576 p. (in Russ.)
2. Besekerskiy V.A., Popov E.P. Teoriya sistem avtomaticheskogo regulirovaniya (The Theory of Automatic Control Systems), Moscow, 1972, 768 p. (in Russ.)
3. Korshunov A.I. Osnovy teorii upravleniya. Rukovodstvo k kursovomu proektirovaniyu (Bases of the Theory of Management. The Management to Course Design), Petrodvorets, 1998, l36 р. (in Russ.)
4. Besekerskiy V.A., Fedorov S.M. Trudy I mezhdunarodnogo kongressa Mezhdunarodnoy federatsii po avtomaticheskomu upravleniyu (Proceedings of the I International ingress of the International Federation on Automatic Control), Moscow, 1961, vol. I, pp. 154-165. (in Russ.)
5. 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, 167 р. (in Russ.)
Data on authors
Artem I. Ermolenko — Cadet, Naval Polytechnic Institute "N. G. Kuznetsov Naval Academy",
Department of Automatic Control Systems Anatoly I. Korshunov — Dr. Sci., Professor; Naval Polytechnic Institute "N. G. Kuznetsov Naval
Academy", Department of Radio Electronics; E-mail: [email protected]
For citation: Ermolenko А. I., Korshunov А. I. Increase in dynamic accuracy of digital tracking systems of APCS by combined management method. Part II. The high rate of mismatch calculation. Journal of Instrument Engineering. 2018. Vol. 61, N 12. P. 1036—1044 (in Russian).
DOI: 10.17586/0021-3454-2018-61-12-1036-1044