УДК 621.396:62-503.53 Дата подачи статьи: 24.12.19
DOI: 10.15827/0236-235X.132.661 -670 2020. Т. 33. № 4. С. 661-670
Метод синтеза адаптивной радиотехнической следящей системы
С.В. Буцев 1, д.т.н, ведущий научный сотрудник, butsev57@pk.ru 1 Научно-исследовательский центр ЦНИИ ВКС, г. Тверь, 170026, Россия
Радиотехнические следящие системы различного целевого назначения функционируют при наличии одновременно неопределенностей как параметров отслеживаемого процесса, так и обобщенного возмущения.
Автором статьи разработан метод синтеза алгоритмов функционирования адаптивной радиотехнической следящей системы, состоящей из адаптивного фильтра и адаптивной системы управления.
Синтез адаптивного фильтра, отслеживающего входной сигнал, включает разработку скользящего алгоритма адаптации на основе непосредственной оценки параметров фильтра, в частности, весовых коэффициентов усиления, формируемых блоком адаптации и используемых в основном фильтре следящей системы.
В статье предлагается использовать двухуровневую структуру, включающую два этапа синтеза: основного контура управления (оптимального регулятора) и контура адаптации (адаптивного регулятора). Оптимальный регулятор синтезируется на основе принципов теории оптимального управления стохастическими процессами при условии, что параметры объекта управления радиотехнической следящей системы постоянны и известны, а внешние возмущения не меняются (либо отсутствуют). Синтез адаптивного регулятора проводится для случая одновременного наличия неопределенности параметров объекта управления радиотехнической следящей системы и внешних возмущений, действующих на него (модель возмущений объекта управления радиотехнической следящей системы), на основе рекуррентной модификации метода идентификации по наименьшим квадратам.
Предложенный подход позволил формализовать задачу функционирования адаптивной радиотехнической следящей системы в условиях неопределенности параметров обобщенного возмущения. Получены новые соотношения для оценок коэффициентов усиления адаптивного фильтра и матрицы передачи адаптивного регулятора.
Разработанная система обеспечивает эффективное функционирование в рамках формализованного описания качества на основе квадратичного функционала, характеризующего точность работы системы и затраты на ее управление.
В статье приводятся некоторые результаты оценки эффективности функционирования синтезированной адаптивной радиотехнической следящей системы.
Ключевые слова: радиотехническая следящая система, объект управления, обобщенное возмущение, управляющее устройство, адаптивный фильтр, адаптивная система управления, адаптивный регулятор, алгоритмы функционирования.
Проектирование и разработка высокоточных радиоэлектронных систем, составной частью которых являются радиотехнические следящие системы (РТСС) различного целевого назначения (например, радиолокационные, радионавигационные, радиопеленгационные), размещаемые на различных носителях (стационарных, возимых, мобильных), - актуальная на сегодня проблема. В общем случае РТСС состоит из объекта управления (ОУ) (например, фазированной антенной решетки) и совокупности устройств (систем управления на основе автоматических регуляторов, следящих приводов и других подсистем), которые обеспечивают управление ОУ и являются для него управляющимустройством (УУ) [1].
Одно из направлений разработки таких систем - использование при синтезе УУ, состоящего, как правило, из оптимальных фильтра и регулятора, методов теорий стохастической фильтрации и оптимального управления [2]. Данные методы позволяют сформировать управление ОУ РТСС, оптимальное по минимуму интегрального функционала качества. Эффективность таких систем во многом зависит от степени соответствия реальных характеристик РТСС и условий их функционирования, таких как параметры отслеживаемого входного процесса, ОУ исполнительной части РТСС и внешних возмущений, действующих на ОУ системы (внешней среды), тем условиям и характеристикам, которые принимаются при их син-
тезе. В реальных условиях функционирования РТСС данное соответствие, как правило, не соблюдается, а классические методы теории автоматического управления [1, 2], как и современные точные методы [3], могут работать неоптимально из-за параметрических и структурных неточностей в описании УУ и ОУ. Для решения этой проблемы необходимо применить адаптивный подход. Однако в случае его применения возникают сложности при идентификации параметров нестационарных моделей входных полезных сигналов, ОУ и внешней среды. Это особенно актуально при возможной скоротечности протекания процессов изменения данных параметров.
В связи с этим возникает необходимость проектирования и разработки адаптивных РТСС, функционирующих в условиях неопределенности, когда отсутствует точная информация о полезном входном сигнале, ОУ и внешней среде и данная информация изменяется в процессе функционирования РТСС, то есть при наличии неопределенности параметров обобщенного возмущения, действующего на систему.
Синтез адаптивных следящих систем в настоящее время хорошо проработан (например [4, 5]). Однако в большинстве публикаций в данной предметной области рассматривается синтез либо алгоритмов адаптивной фильтрации полезного сообщения [3], либо адаптивной системы управления (адаптивного регулятора) при наличии только неопределенности внешних возмущений, действующих на ОУ системы [4, 6], или только неопределенности, вызванной изменением (нестационарностью) параметров ОУ исполнительной части системы в процессе ее функционирования [7, 8].
В реальных условиях функционирования РТСС возникают ситуации, когда неопределенными одновременно могут оказаться как параметры отслеживаемого входного процесса, так и ОУ РТСС и внешних возмущений, действующих на ОУ системы (обобщенное возмущение). На практике такие ситуации чаще всего возникают при размещении РТСС и их функционировании на мобильных носителях (наземных, воздушных и др.), когда определяющими являются нестационарность параметров входных полезных сигналов, скоротечность процессов, связанных с вибрационными, ветровыми, температурными, помеховыми и другими негативными воздействиями на ОУ. Поэтому возникает актуальная задача синтеза РТСС, функционирующей при наличии не-
определенности параметров обобщенного возмущения.
Цель - разработать метод синтеза алгоритмов функционирования адаптивной РТСС при наличии неопределенности (нестационарности) параметров обобщенного возмущения.
Постановка задачи
Пусть РТСС с ОУ, описываемым дискретными уравнениями вида
г(/+1) = А (/>(/) + В{г)и{г) + ю (/), (1)
а (г) = С(/')г(/),
предназначена для отслеживания входного полезного процесса (сообщения), например, координат объекта
а(г+1) = Ф(г)а(г) + Г(/)п(/), (2)
при наличии наблюдений (измерений) вида
></) = £д(/)[а(/)- ¿(0] +и(/), (3)
где г{г) - и-мерный вектор состояния; а (г) -п-мерный вектор выхода; А(г), Б(г), С(г) - матрицы размером пхп; и(г) - и-мерный вектор управления; а(г) - и-мерный вектор полезного сообщения; у(г) - и-мерный вектор наблюдения (при синтезе полагается, что модель дискриминатора представлена в виде линеаризованного динамического эквивалента); Ф(г) - матрица состояния размером пхп; Г(г) - матрица при шумах возбуждения размером пхп; ц(г) -п-мерный вектор шумов возбуждения, представляющий собой нормальный белый шум с нулевым математическим ожиданием (МО) и корреляционной матрицей (КМ) М{ц(г)цТ()} = = Бфу, М{...} - оператор МО, Б- - симметрическая неотрицательно определенная матрица дисперсий шумов возбуждения; га (г) - п-мер-ный вектор возмущений, представляющий собой нормальный белый шум с нулевым МО и
© симметрическая неотрицательно определенная матрица дисперсий шумов возмущения; кд - матрица коэффициентов усиления дискриминатора размером пхп; п(г) - п-мерный вектор шумов измерения, представляющий собой нормальный белый шум с нулевым МО и КМ М{п(г)пТ(]')} = QЪij, Q - симметрическая положительно определенная матрица дисперсий ошибок измерения.
Необходимо найти соответствующее допустимое управление иопт.(г), определенное на интервале дискретизации [0, т] и зависящее от текущих значений отслеживаемого процесса, такое, чтобы показатель качества вида
/ ( т ) =
т-1 I Т
= xi [7(О - к (ОМО] б"1 [у(1) - к (ОМО] +
1=0
+ ит (1) Я -и(1)} (4)
(Я - симметрическая положительно определенная весовая матрица) принимал минимальное значение на множестве допустимых управлений и (и(г) с и) при условии, что некоторые элементы матриц Ф, Г, Б, характеризующие отслеживаемый процесс, и некоторые элементы матриц А, В, А», характеризующие состояние ОУ РТСС, могут быть неизвестны, а область пространства возможных значений ошибки функционирования (рассогласования) РТСС 8а(/') = а(/) - а(г') должна быть существенно меньше некоторой заданной величины, например, полуширины дискриминационной характеристики Д (|8а(г) | << Д).
Так как рассматриваемая РТСС наблюдаема и управляема, при синтезе адаптивной РТСС будем использовать принцип разделения [5], согласно которому УУ РТСС должно состоять из адаптивного фильтра для оценки полезного сообщения при наличии неопределенности его параметров и адаптивной системы управления для управления ОУ РТСС при наличии неопределенности его параметров и внешних возмущений, действующих на него.
Известны несколько подходов к синтезу адаптивных фильтров [3-5]. Один из них состоит в том, что неизвестные параметры включаются в число оцениваемых и рассматривается задача фильтрации расширенного вектора состояния [3]. Существует подход к синтезу адаптивных фильтров на основе оценки неизвестных статистических характеристик сообщения и помех с использованием их для изменения параметров фильтра [5]. Данные подходы приводят к аналитическим усложнениям алгоритмов фильтрации и значительным вычислительным затратам. Еще один подход к синтезу адаптивных фильтров основан на том, что оценки неизвестных параметров формируются отдельным блоком адаптации и используются для корректировки параметров основного блока фильтрации [4]. Реализация его проще, чем в [3, 5], однако он также требует значительных вычислительных затрат, так как на каждом шаге адаптации необходимо пересчитывать матрицы дисперсий ошибок фильтрации и коэффициентов усиления фильтра. В статье предлагается обойти данную трудность за счет реа-
лизации скользящего алгоритма адаптации на основе непосредственной оценки параметров фильтра, в частности, весовых коэффициентов усиления, формируемых блоком адаптации и используемых в основном фильтре [9, 10].
Для синтеза адаптивной системы управления при одновременном наличии неопределенности параметров ОУ РТСС и внешних возмущений, действующих на него, предлагается использовать двухуровневую структуру [2], включающую два этапа синтеза: основного контура управления (оптимального регулятора) и контура адаптации (адаптивного регулятора).
Оптимальный регулятор синтезируется на основе принципов теории оптимального управления стохастическими процессами по критерию минимума показателя (4) при условии, что параметры ОУ РТСС постоянны и известны, а внешние возмущения не меняются (либо отсутствуют) (эталонная модель РТСС) [6, 11].
Синтез адаптивного регулятора проводится для случая одновременного наличия неопределенности параметров ОУ РТСС и внешних возмущений, действующих на него, по критерию минимума ошибки функционирования модели возмущения РТСС на основе рекуррентной модификации метода идентификации по наименьшим квадратам [6].
Синтез адаптивного фильтра
Для синтеза адаптивного фильтра воспользуемся подходом из работы [9]. Если матрицы Ф, Г, Б и Q полностью определены, то оптимальным по критерию минимума показателя (4) для оценивания сообщения (2) при наблюдении (3) является фильтр Калмана [12], описываемый соотношениями: а (г) = а(г / г -1) + К (г) х
х [Я0 " / г -1) + кл а(г)], а(1 /1 -1) = Ф(оа(г -1), (5)
Да (1 / 1 - 1) = Ф (1) Да (1 - 1)ФГ (1) + Г(1) ДГ (1), Dа (О = Да (1 / 1 - 1) Х
X 11 - [Па (1 / 1 - 1) + а(кякТд )-1 ]-1 } Да (1 / 1 - 1),
К (1) = Да (1) К б Л где а (1) - оптимальная оценка отслеживаемого сообщения; а,(г/г -1) - экстраполированная оценка отслеживаемого сообщения; К(г) - матрица коэффициентов усиления фильтра.
При неопределенности матриц Ф, Г, оптимальную матрицу К(г) определить по (5) не-
возможно, и в фильтре будет использоваться некоторая неоптимальная матрица коэффициентов усиления К (/). Тогда оценки, формируемые данным фильтром, также будут неоптимальными:
а* (г) = а* (г /г -1) + К * (г) х
х [^(г) -кда (г / г -1) + кд а(г)],
а* (г /г -1) = Ф*(г)а* (г -1) , (6)
где Ф*(/) - аналог матрицы состояния Ф(/) в оптимальном фильтре, но с другими параметрами.
Поставим задачу отыскания такой матрицы коэффициентов усиления К (/), при которой функционал (4) достигает своего минимума при неопределенности Ф, Г, Для этого воспользуемся методикой, предложенной в [13] и заключающейся в том, что оптимизация (4) рассматривается с точки зрения теории оптимального управления. В [9] показано, что (4) имеет минимум при тех значениях К (/), которые минимизируют матрицу дисперсий ошибок экстраполяции сообщения Оа(г/г - 1) из (5), зависящую от К*(/). При этом полагается, что Щ/// - 1) является состоянием системы, а К*(г) - управлением, которое нужно оптимизировать на интервале функционирования системы [0, т] так, чтобы функционал (4) имел минимум, причем матрица К рассматривается как функция вектора к(/) = [&п(0, ..., £„„(/)], компоненты которого ку(/) являются элементами данной матрицы.
Для нахождения оценки К * (г), минимизирующей функционал (4), вводится вспомогательная функция:
Н [К* (г), 1(г +1)] = [.у (г) - кя (г) 8а (г )]Т 2-1 х х [у (г) - к (г) 8а (г)] + иТ (г) К~1и (г) + (7)
+ХТ (г +1)К * (г),
где А,(/+1) - функционально сопряженные неопределенные множители Лагранжа [14].
Тогда канонические уравнения для К (/+1) и Х(г') можно записать в виде
дН [ К * (г), А,( г +1)]
K * ( i + 1 ) =
Ч i ) =
d^( i +1) дИ [ K * (i ), x(i +1)]
dK '(i )
. . dâ*T ( i +1/i ) . = X(i +1)-k -Ц—'-!-Qx
( ) д dK (i ) Q
x I~y (г +1) - £да* (г + l/г) + ka à (i + l) J,
с граничными условиями на концах [0, т] -К*(0) = ко, Х(т) = 0 [14].
Используя подход, предложенный в [3], получаем следующие уравнения для матрицы оценок коэффициентов усиления фильтра
К'(I) :
К * (г +1)= К * (г) +
да*Т (г +1/г)
+P ( i +1)-
■кд Q-I x
dK * (i)
x ^y (i +1) - кдa (i + l/г) + кд ä (г +1)] P-1 ( i + 1)= P-1 ( i ) + кд x
dâ*T (i + Vi) _ifdâ*T (i +1/i)'"
(8)
dK * ( i )
Q "
dK * ( i )
kT -
-[y (г +1) - кла (i + l/г) + kRä (г +1)] Q'1 x
xk„
d dâ*T (i + Vi)
дК* (г) дК* (г) ' где Р(/'+1) имеет смысл дисперсии ошибки оценки коэффициентов усиления фильтра.
Первая и вторая производные, входящие в (8), получаются путем дифференцирования (6) по К (/). Тогда имеем:
да*Т (г + 1/г)
dK * (i)
dâ*T (0 -1)
dK * (i) xO*T (i + 1),
d dâ*T (i +1/i) dK' (i ) dK' (i ) =
d dâ*T( i/i -1)
(9)
[1- кдИ'T (i)]-2k,
dot* (i/i -1)
dK * (i )
[дК * (г) дК * (г) хФ*Т (г + 1).
Объединив (6), (8) и (9), получаем уравнения адаптивного фильтра:
а* (г +1/г) = Ф* (г +1 )ос* (г) ,
а* (г +1) = а* (г + 1/1) + К * (г + 1)х
х (г +1) - кда* (г + 1/г) + А:да (г + 1)],
. , ,да*Т (г +1/г) .
К (г +1)= К (г)+^ (г +1) + кд 2-1 х
х ^(;' +1) - кда* (г + 1/г) + Ада (г +1)],
w ч ./s dâ*T ( i +1/i ) .
P - (i + 1) = P - ( i ) + kg j* + 1 ) Q -1 X
dK * ( i )
( dâ*T (i + Vi )
dK * (i )
(10)
kT -
T
-[у(1 + \)-кДа (г + 1/г) + £да(г + 1)] ¿Г1 кд х
9 9а *т (1 + V1)
х---- ,
9К (1) 9К (1)
9а*т (1 +1/1) 9К* (I) =
= [1 - кдК* (')] + [у (') - кда* У -1) + кда (')]} X
хФ*т (1 + 1) ,
9 еа *т (1+11)_ 9К * (I) 9К * (I) =
={ е еа*т(ф-1)[1 -кк-т(1)]-2к ес^)}х
|ек* (1) ек* (1) 1 д ()] д ек* (1) ]
хФ*т (1 +1).
Из (10) следует, что адаптивный фильтр состоит из двух блоков (рис. 1): основного блока фильтрации, описываемого первым и вторым соотношениями, и блока адаптации, описываемого третьим, четвертым, пятым и шестым соотношениями, формирующего оценки коэффициентов усиления К (1) , используемые в основном блоке фильтрации вместо значений К(г).
Предложенный подход к синтезу адаптивного фильтра оправдан, если неизвестными являются и Г(г), так как только они влияют на К(г). В случае неопределенности Ф(г) использованием К (1) в фильтре можно только частично компенсировать ее неопределенность, так как она определяет не только К(г), но и Ф*(г) в фильтре. Поскольку в реальных условиях функционирования РТСС формирование Ф(г) осуществляется на основе некоторой известной
статистической информации, этого достаточно для получения эффективной оценки полезного сообщения, формируемого данным фильтром.
Синтез адаптивной системы управления
Для синтеза адаптивной системы управления воспользуемся идентификационным подходом [5], позволяющим в условиях неопределенности параметров ОУ РТСС и внешних возмущений получать оценки возмущений, действующих на ОУ РТСС, и использовать их для формирования управления, компенсирующего данные неопределенности [13].
Представим эталонную модель функционирования ОУ РТСС в виде разностного уравнения:
гм(/) = Лм(/>м(/) + Вм(1)им(1),
¿м(0 = С(г>м(0, (П)
где матрицы АмЦ). ВмЦ). С{г) и векторы хмЦ). им{г), ам (г) имеют смысл и размерность, аналогичные в (1) матрицам /1(7). В(г), ('(¡) и векторам 2(1). м(/), а (г). При этом в соответствии с подходом [5] будем полагать, что матрицы Ам(г) = А (г) - Г(г) и ВмЦ) = В (г) - О(г), а векторы гмЦ) = г(г') - 8^(г) и ымЦ) = и(г) - 8м(г), где О(г) - матрицы возмущений размером пхп; Ьг(г), 8и(г) - п-мерные векторы возмущения векторов состояния и управления соответственно.
Вычитая из первого и второго соотношений (1) соответственно первое и второе соотношения (11), получаем модель возмущения ОУ РТСС:
8г(7+1) = ^(0 8г(г) + О(г) 8и(г) + ю (г), (12)
â( kT ) ä(kT)
( kT ) У ( kT )
ï(kT)
Блок адаптации
Основной блок фильтрации
1С
( kT )
* ( kT )
k
V i
â'(kT/kT - T)
Ф' (kT )• S-1
Оптимал ьный регулятор
* ( kT )
Модель объекта управления
Адаптивный фильтр
Л kT )
kT )
Объект управления
Цкт)
Su' ( kT )
Адаптивный регулятор
Se( kT )
,(кТ)
ЦкТ)
Рис. 1. Блок-схема адаптивной РТСС Fig. 1. Block diagram of the adaptive RTMS
u
u
Se(/) = ¿(0 - äM (0= C(i)<bz(i). (13)
Синтез адаптивной системы управления сводится к определению управления 8м(/), компенсирующего возмущение Ss(z'), которое может быть получено, например, на основе одношагового закона оптимального управления [11] в предположении, что скорость изменения параметров ОУ исполнительной части РТСС меньше скорости адаптации; переменные Ss(z') поддаются измерению, а погрешность их измерения мала.
Нахождение управления 8м(/) требует, чтобы матрицы F(i) и G(i) были известны в каждый момент времени управления. Воспользуемся идентификационным подходом [5]. Для этого, записывая i-ю строку неизвестных параметров в (12) в виде вектора
ef(i) = [/л(о,f(i), gfl(0,gn(0], (14)
i = 1, ..., n,
и представляя выходы и входы в (12) вектором xT(i) = [8z1(i), ..., 8zn(i), 8м1 (i), ..., ÖMn(i)], а состояние системы вектором SzT(i) = [8zi(i), ..., Szn(i)], трансформируем (12) к виду
Sz(i+1) = xT(i)ei(i) + ffli(i), i = 1, ..., n. (15) Минимизируя ошибки функционирования модели возмущения ОУ РТСС (15) относительно вектора ei(i) по критерию предсказания вида
Jпр. (e)={£ [ц (i+!)- x (oe, (о]} min, (16)
на основе метода идентификации по наименьшим квадратам [6] получаем следующие соотношения для его оценки:
e, (,+1)=e, (,)+w (i) x(i) [ (i+i) - xT (i)e, (i)],
W(i +1) = W(i) - у (i)W(i)x(i)xT (i)W(i), (17)
¥(i) = [xT (i)W(i)x(i) +1]-1,
W(0) = ц/, ц >> 1, 8(0) = 90. После определения по (17) оценок матриц F(i), G(i) может быть сформирован закон оптимального управления ОУ РТСС с учетом возмущающих факторов, удовлетворяющий требованиям уравнений (1), (12) и минимизирующий показатель качества (4). Используя уравнение, описывающее отслеживаемый процесс (2), и уравнение функционирования ОУ РТСС (11), введем расширенные уравнения состояния и ошибки слежения:
zp(i+\) = Ap(i)zp(i) + Bp(i)u(i),
8ам (г) = а (г) - äM (г) = Ср (г)zp (г),
где Ар (г) =
Ф(г) 0 0 AM (г)
> B (0 =
0
Bm (0.
Сp(i) = [I - С(i) ], zp(i) =
а(г)
ZM (i)
ниям Ss(/); z (i) =
Применяя принципы теории оптимального управления стохастическими процессами [5, 11], можно получить управление ОУ РТСС иопт.(г'), минимизирующее показатель качества (4) с учетом обобщенного возмущения и удовлетворяющее требованиям соотношений (1) и (12): и0111 (г) = и (г) + 8и (г) =
= - К' опт. (г) 2р (г) -8К * опт. (г) 8z (г) =
= -[К-1 + БТр (г)Б(г + 1)Бр (г)]-1 х
хВр (г)Б(г +1)Ар (г) Zp (г) -
-[К-1 + £Т (г)2(г)£(г)2^(18)
где 5 (г) = [ Ар (г) - Бр (г) К' ^ (г)] Б (г +1) + ; Б(т) = Q1; К*ош (г) - матрица передачи оптимального регулятора; 8 К*пт (г) - матрица передачи адаптивного регулятора; и*(г) - оптимальное управление при отсутствии возмущений; 8и*(г) - оптимальное управление, минимизирующее ошибки слежения системы по возмуще-а (г)
_ z (г) ]
Из (18) следует, что адаптивная система управления ОУ РТСС при обобщенном возмущении включает двухуровневую структуру: оптимальный регулятор при отсутствии возмущений (первое слагаемое в (18)) и адаптивный регулятор при наличии возмущений (второе слагаемое в (18)).
Блок-схема синтезированной адаптивной РТСС приведена на рисунке 1.
Пример
Для оценки эффективности разработанного метода воспользуемся подходами, предложенными в [9, 10]. В качестве показателя эффективности функционирования адаптивной РТСС принимается выигрыш в точности адаптации Р [10], который достигается использованием адаптивной РТСС вместо неадаптивной, синтезированной при известных характеристиках отслеживаемого процесса, ОУ РТСС и внешних возмущений, действующих на ОУ системы. Оценить его можно отношением дисперсий ошибок слежения в неадаптивной
D4 (i) и адаптивной (i) РТСС вида
I
ß=.
Dt, (i)
D, (i) '
(19)
где °ЪА (1) = А (1)| „ (0 = ^ , (г) = Оя(1) ^ ,
5Н w ^ / и (г) = и (г) 5
въ (0 = М {[а(0 - а(0][а(0 - а(/)Г} .
Выигрыш в точности адаптации Р рассчитывался по соотношению (19) для случаев наличия параметрической и структурной неопределенностей характеристик отслеживаемого сообщения (2) с учетом соотношений (10), (17), (18). При параметрической неопределенности полагалось, что Ф = е-, Ф* = е- (в
матрице состояния параметр у отличается от
*
его значения у , принятого при синтезе системы, только величиной), а при структурной -
Ф = е- , Ф* = е- (параметр у отличается и
ч >1!
структурой, и величиной от его значения у , принятого при синтезе системы). Результаты расчетов получены для установившегося режима фракционирования РТСС, задаваемого условием АО - АС - 1) = 0,01, при следующих исходных данных: А = В = Г = С = кд = 1, Аю = 0,1, 0 = 0,25, Я = 0,1, = 0,5, ц = 10. Семейство зависимостей выигрыша в точности адаптации Р = Ду), рассчитаное при нескольких фиксированных (эталонных) значениях параметра ушах = у* (0,01; 1; 100) для параметрической (сплошные линии) и структурной (штриховые линии) неопределенностей характеристик отслеживаемого сообщения приведены на рисунке 2. Зависимости рассчитывались для отклонения параметра у от его фиксированного значения у* = ушах. Видно, что выигрыш в точности адаптации Р зависит от того, насколько параметры реального отслеживаемого процесса отличаются от значений параметров, принятых при синтезе системы. Если у = ушах, то Р = 1 и выигрыша нет. По мере увеличения отклонения Ду = ушах - у увеличивается выигрыш в точности адаптации Р, причем Р при структурной неопределенности выше, чем при параметрической. Поэтому структурная неопределенность оказывает на точность функционирования РТСС более сильное негативное влияние, чем параметрическая. Вместе с тем анализ зависимостей, приведенных на рисунке 2, показывает, что введение адаптации в РТСС на основе разаботанного метода позво-
ляет существенно снизить влияние неопределенностей характеристик отслеживаемого сообщения на эффективность функционирования РТСС, в частности, практически полностью устранить влияние параметрической неопределенности дисперсии аддитивного шума измерений на работу РТСС при отклонении значения рассматриваемого (отслеживаемого) параметра у от его фиксированного значения у* на 25-35 %. Кроме того, исследование частных задач фильтрации параметров сообщения следящей системой, приведенное в [10], показало, что отказ от нормировки на матрицу дисперсий аддитивного шума измеренй 0 в разработанных в статье соотношениях приводит лишь к некоторому ухудшению процесса адаптации, заключающемуся в его затягивании по времени.
Исследования показали, что использование предлагаемого метода и реализующих его алгоритмов при функционировании высокоточных радиолокационных систем на основе фазированных антенных решеток позволяет повысить точность оценки (экстраполяции) параметров траектории движения летательных аппаратов (например, их координат), описываемой полиномом третьей степени, по сравнению с оцениванием на основе стандартных алгоритмов калмановской фильтрации [15, 16] при наличии неопределенности (отклонения) параметров отслеживаемого процесса от предполагаемых (эталонных, принятых при синтезе системы) и ограниченной выборке измерений данных параметров в среднем на 15-25 % [17],
ß3 2,5
2
1,5
1
0,5
\ 4
\ \ S
\ 4 \ \
> N N t >
0,001
0,01
0,1
10
100
Рис. 2. Семейство зависимостей выигрыша в точности функционирования адаптивной РТСС относительно неадаптивной
Fig. 2. Family ofpayoff dependencies in the accuracy of the functioning of the adaptive RTMS relatively non-adaptive
v
что, как следствие, приводит к повышению точности сопровождения данных объектов радиолокационными системами при недостаточной априорной информации о полезном входном сигнале и параметрах фазированной антенной решетки системы.
Выводы
В статье предложен подход к синтезу управляющего устройства адаптивной РТСС, включающий в себя синтез алгоритма адаптивной фильтрации и алгоритма адаптивного управления по обобщенному возмущению. Синтез адаптивного фильтра проведен на основе метода адаптивной фильтрации, предполагающего непосредственную регулировку коэффициентов усиления фильтра, а синтез алгоритма адаптивного управления - на основе параметрической идентификации по методу наименьших квадратов и методов теории оптимального стохастического управления. Предложенный подход позволил формализовать задачу функционирования адаптивной РТСС в условиях неопределенности параметров отслеживаемого процесса, состояния ОУ исполнительной части РТСС и внешних возмущений, действующих на ОУ данной системы (обобщенного возмущения). Получены новые соотношения для оценок коэффициентов усиления адаптивного фильтра и матрицы передачи адаптивного регулятора при наличии возмущений ОУ.
Разработанная на основе предлагаемого метода адаптивная РТСС обеспечивает эффек-
тивное функционирование в рамках формализованного описания качества на основе показателя, определенного квадратичным функционалом, который характеризует точность работы системы и затраты на ее управление, что подтверждается приведенными в статье результатами оценки эффективности функционирования синтезированной адаптивной РТСС в условиях воздействия на нее обобщенного возмущения.
Практическая реализация алгоритмов функционирования РТСС, реализующих разработанный метод, позволит повысить эффективность функционирования РТСС различного целевого назначения, в частности, радиолокационных систем на основе фазированных антенных решеток, размещаемых на мобильных носителях (наземных, воздушных и др.), когда определяющими негативными факторами, влияющими на их работу, являются нестационарность входных полезных сигналов, скоротечность процессов, связанных с негативными воздействиями на систему внешней среды, например, при сопровождении траекторий движения летательных аппаратов в условиях наличия неопределенности параметров отслеживаемых координат, вибрационных и помеховых воздействий на систему.
Предлагаемый подход также может быть использован при разработке оптико-электронных локационных систем, требующих получения и обеспечения более высокой точности функционирования по сравнению с радиотехническими системами.
Литература
1. Репин В.Г., Тартаковский Г.П. Статистический синтез при априорной неопределенности и адаптации информационных систем. М.: Сов. радио, 1977. 432 с.
2. Сейдж Э.П., Мелс Дж. Идентификация систем управления. М.: Наука, 1974. 246 с.
3. Первачев С.В., Перов А.И. Адаптивная фильтрация сообщений. М.: Радио и связь, 1991. 160 с.
4. Фрадков А.Л. Адаптивное управление в сложных системах. М.: Наука, 1990. 296 с.
5. Красовский А.А. Справочник по теории автоматического управления. М.: Наука, 1987. 712 с.
6. Цыпкин Л.З. Адаптивно инвариантные дискретные системы управления // Автоматика и телемеханика. 1991. № 5. С. 96-124.
7. Куржанский А.Б. Управление и наблюдение в условиях неопределенности. М.: Наука, 1977. 392 с.
8. Ци Э., Бар-Шолом Я. Активно адаптивное управление нелинейными стохастическими системами // ТИИЭР. 1976. № 8. С. 52-61.
9. Буцев С.В. К вопросу оценки чувствительности следящих систем к неопределенностям характеристик отслеживаемого процесса // Автоматика и телемеханика. 1992. № 8. С. 25-31.
10. Буцев С.В., Зайцев А.Г., Тикменов В.Н. Оценка чувствительности адаптивных следящих систем к неопределенностям характеристик входных сигналов // Радиотехника. 2004. № 5. С. 29-34.
11. Егоров А.И. Основы теории управления. М.: Физматлит, 2004. 504 с.
12. Синицын И.Н. Фильтры Калмана и Пугачева. М.: Логос, 2007. 772 с.
13. Буцев С.В. Алгоритмы функционирования адаптивной оптической системы фазового сопряжения // Квантовая электроника. 1996. Т. 23. № 8. С. 753-758.
14. Перов А.И. Статистическая теория радиотехнических систем. М.: Радиотехника, 2003. 400 с.
15. Мазурова М.А. Выбор коэффициентов диффузии для фильтра Калмана с компенсацией скоростной ошибки по дальности // Радиотехника. 2019. № 10. С. 32-42. Б01: 10.18127^00338486-201910(15)-06.
16. Перов А.И. Синтез и анализ оптимального алгоритма контроля слежения по задержке в навигационном приемнике // Радиотехника. 2019. № 9. С. 5-13. Б01: 10.18127/|00338486-201909(14)-01.
17. Буцев С.В., Горюнов Г.А., Прохоров В.А., Сердешнов Е.А., Шепеть И.П. Повышение качества оценки параметров траектории движения летательного аппарата // Радиотехника. 2001. № 5. С. 96-100.
Software & Systems Received 24.12.19
DOI: 10.15827/0236-235X.132.661 -670 2020, vol. 33, no. 4, pp. 661-670
Method of synthesis of adaptive radio technical monitoring system S.V. Butsev 1, Dr.Sc. (Engineering), Leading Researcher, butsev57@bk.ru
1 Research Center of Central Research Institute of Air and Space Forces, Tver, 170026, Russian Federation
Abstract. Radio technical monitoring systems for various mission objectives operate in the presence of both uncertainties in the parameters of the monitored process and the uncertainty of a generalized disturbance.
The author of the paper has developed a method for synthesizing algorithms for the functioning of an adaptive radio technical monitoring system, consisting of an adaptive filter and an adaptive control system.
The synthesis of an adaptive filter for a monitoring input signal includes the development of a sliding adaptation algorithm based on a direct estimate of the filter parameters, in particular, the gain weights generated by the adaptation unit and used in the main filter of the monitoring system.
The paper proposes to use a two-level structure, which includes two stages of synthesis: the main control loop (optimal controller) and the adaptation loop (adaptive controller). The optimal controller is synthesized on the basis of the principles of the theory of optimal control of stochastic processes, provided that the parameters of the control object of the radio engineering tracking system are constant and known, and external disturbances do not change (or are absent). The synthesis of the adaptive controller is carried out for the case of the simultaneous presence of the uncertainty of the parameters of the control object of the radio technical monitoring system and external disturbances acting on it (model of disturbances of the control object of the radio technical monitoring system), based on the recurrent modification of the least-squares identification method.
The proposed approach makes it possible to formalize the problem of the functioning of an adaptive radio technical monitoring system under conditions of uncertainty of generalized disturbance parameters. New relationships are obtained for the evaluation of gains of an adaptive filter and an adaptive regulator transmission matrix.
The adaptive radio technical monitoring system developed on the basis of the proposed method ensures efficient functioning within the formalized quality description on the basis of the indicator determined by the quadratic function which characterizes the accuracy of the system operation and its control costs.
The paper provides some efficiency evaluation results of the synthesized adaptive radio tracking system functioning.
Keywords: radio technical monitoring system, control object, generalized disturbance, control device, adaptive filter, adaptive control system, adaptive regulator, operation algorithms.
References
1. Repin V.G., Tartakovsky G.P. Statistical Synthesis with a Priori Uncertainty and Adaptation of Information Systems. Moscow, 1977, 432 p. (in Russ.).
2. Sage A.P., Melsa J.L. System Identification. Moscow, 1974, 246 p. (in Russ.).
3. Pervachev S.V., Perov A.I. Adaptive Message Filtration. Moscow, 1991, 160 p. (in Russ.).
4. Fradkov A.L. Adaptive Control in Complex Systems. Moscow, 1990, 296 p. (in Russ.).
5. Krasovskiy A.A. Automatic Control Theory Reference Guide. Moscow, 1987, 712 p. (in Russ.).
6. Tsypkin L.Z. Adaptively invariant discrete control system. Automation and Remote Control, 1991, no. 5, pp. 96-124 (in Russ.).
7. Kurzhansky A.B. Control and Observation in the Condition of Indeterminacy. Moscow, 1977, 394 p. (in Russ.).
8. Tsi E., Bar-Sholom Ya. Actively adaptive control of nonlinear stochastic systems. Proc. IEEE, 1976, no. 8, pp. 52-61 (in Russ.).
9. Butsev S.V. On an estimate for the sensitivity of tracking systems to uncertainties in the characteristics of the tracked process. Automation and Remote Control, 1992, no. 8, pp. 25-31 (in Russ.).
10. Butsev S.V., Zaytsev A.G., Tikmenov V.N. Evaluation of the adaptive servo systems sensitivity to uncertainties in input signal characteristics. Radiotekhnika, 2004, no. 5, pp. 29-34 (in Russ.).
11. Egorov A.I. Fundamentals of Management Theory. Moscow, 2004, 504 p.
12. Sinitsyn I.N. Kalman andPugachev Filters. Moscow, 2007, 772 p. (in Russ.).
13. Butsev S.V. Algorithms for the operation of an adaptive optical phase conjugation system. Quantum Electronics, 1996, vol. 23, no. 8, pp. 753-758 (in Russ.).
14. Perov A.I. Statistic Theory of Radio Engineering Systems. Moscow, 2003, 400 p. (in Russ.).
15. Мазурова М.А. Kalman filter with range-Doppler coupling compensation: choice of diffusion coefficients. Radiotekhnika, 2019, no. 10, pp. 32-42 (in Russ.). DOI: 10.18127/j00338486-201910(15)-06.
16. Perov A.I. Synthesis and analysis an optimal code delay lock detector algorithm in a satellite navigation receiver. Radiotekhnika, 2019, no. 9, pp. 5-13 (in Russ.). DOI: 10.18127/j00338486-201909(14)-01.
17. Butsev S.V., Goryunov G.A., Prokhorov V.A., Serdeshnov E.A., Shepet I.P. Improving the quality of estimating the parameters of the motion trajectory of the aircraft. Radiotekhnika, 2001, no. 5, pp. 96-100 (in Russ.).
Для цитирования
Буцев С.В. Метод синтеза адаптивной радиотехнической следящей системы // Программные продукты и системы. 2020. Т. 33. № 4. С. 661-670. DOI: 10.15827/0236-235X.132.661-670.
For citation
Butsev S.V. Method of synthesis of adaptive radio technical monitoring system. Software & Systems, 2020, vol. 33, no. 4, pp. 661-670 (in Russ.). DOI: 10.15827/0236-235X.132.661-670.