Экстраполяция в дискретных системах с мультипликативными возмущениями при неполной информации | Вестн. Том. гос. ун-та. Управление, вычислительная техника и информатика. 2019. № 47. DOI: 10.17223/19988605/47/6

Экстраполяция в дискретных системах с мультипликативными возмущениями при неполной информации

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

Extrapolation in discrete systems with multiplicative perturbations at incomplete information.pdf Задачи оценивания для дискретных систем с мультипликативными возмущениями рассматривались в работах [1-3]. Аналогичные задачи оценивания для непрерывных систем с мультипликативными возмущениями изучались в статье [4]. Для моделей систем с неопределенными параметрами такие задачи были рассмотрены в [5]. Задачи фильтрации и экстраполяции для дискретных систем с неизвестным входом и аддитивными гауссовскими возмущениями изучались в [6-13]. В [7-8] для вычисления оценок неизвестного входа использовался МНК, в работах [9-13] для улучшения качества оценивания применялись алгоритмы непараметрического сглаживания. В настоящей статье рассматривается задача экстраполяции в дискретных системах с мультипликативными и аддитивными возмущениями с неизвестными входом и неизвестными параметрами. Задача решается на основе принципа разделения с помощью алгоритмов МНК и с использованием сглаживающих процедур (скользящего среднего и непараметрического сглаживания). 1. Постановка задачи Пусть модель объекта с мультипликативными возмущениями и неполной информацией описывается разностным уравнением: m x(k +1) = (A + АЛ)x(k) + (B + AB)u(k) + £ Asx(k)0S (k) + f (k) + q(k), x(0) = x0, (1) s=1 где x(k) e Rn - вектор состояния, u(k) e Rp - известный вход, f(k) - неизвестный вход; x0 - случайный вектор (предполагаются известными дисперсионная матрица N0 = M{(x0 - x0)(x0 - x0 )T} и математическое ожидание x0 = M{x0}); A, B, As (s = 1...m) - заданные матрицы; АЛ, AB - матрицы неизвестных параметров; q(k) и (k) - векторные гауссовские случайные последовательности со следующими характеристиками: M{q(k)} = 0, M{q(k)qT(J)} = Q5j , M{0X (k)} = 0, M{0X (k)06J (J)} = ®s5j. Здесь 5 j - символ Кронекера. Канал наблюдений имеет вид: y(k) - Sx(k) + v(k), (2) где v(k) - гауссовская случайная последовательность с характеристиками: M(v(k)} - 0, M(v(k)vT (j)} - V5^ . Предполагается, что последовательности q(k), 0s (k), v(k) независимы между собой, пара матриц S, A детектируема, матрицы Q, &s (s - 1...m) - неотрицательно определены, V> 0. По информации, поступившей в момент k e [0; T ], требуется найти оценку прогноза xc(k +1) на основе минимизации следующего критерия: J[0;T ] - M(^£eT(k)R(k)e(k)}, (3) k-0 где R(k) > 0 - весовая матрица, e(k) - x(k) - X(k) - вектор ошибок. 2. Синтез оптимального экстраполятора Представим модель (1) в следующем виде: m x(k +1) - Ax(k) + Bu(k) + £ Asx(k^ (k) + r(k) + q(k), x(0) - x0, (4) s-1 где r(k) - AAx(k) + ABu(k) + f (k) - неизвестный вектор. Для решения задачи экстраполяции воспользуемся принципом разделения. Это означает, что сначала мы построим оценки вектора xc(k) в предположении, что вектор r(k) известен, затем строятся оценки вектора r(k) в предположении, что оценка прогноза вектора состояния x(k) известна. Тогда для построения оценки можно использовать рекуррентный алгоритм экстраполятора Калмана x(k +1) - Ax(k) + Bu(k) + r(k) + K(k)(y(k) - Sx(k)), x(0) - x0, (5) где K(k) - матрица коэффициентов передачи экстраполятора. Найдем матрицу K(k), обеспечивающую минимум критерия (3). Для этого запишем уравнение для вектора ошибок e(k), вычитая из уравнения (4) уравнение (5): e(k +1) - x(k +1) - x(k +1) - m (6) - (A - KS)e(k) + £ Ax(k)0,(k) + q(k) - Kv(k). s-1 Учитывая (6), уравнение для матрицы N (k) - M{e(k )e(k )T} получим в виде следующего разностного уравнения: m N(k +1) - (A -K(k)S)N(k)(A -K(k)S)J +£ASN(k)Aj + s-1 (7) m + £ Asx(k)x(kf Aj + Q + K(k)VK(kf, N(0) - N0. s-1 Представим критерий (3) в виде T J [0;T ] - tr N0R(0) + £ tr N(k) R(k), (8) k-1 где tr - след матрицы. Тогда, подставив в (8) формулу (7) со сдвигом на один такт, в результате получим T m J[0;T] - trN0R(0) + £tr[(A - K(k - 1)S)N(k -1)(A - K(k - 1)S)J + £AsN(k -1)Aj + k-1 s1 (9) m + £ Asx(k -1)x(k - 1)J Aj + Q + K(k - 1)VK(k - 1)J ]R(k). s-1 Используя правила дифференцирования функции след (tr) от произведения матриц [14]: д tr AXB T T д tr ATXBT = A1 B , -= BA, (10) дх дх ( ) из уравнения =0, (,1) 8K учитывая, что матрица R(k) не вырождена, получим выражение для матрицы K(k): K (k) = AN (k )S T(SN (k )S T + V)-1 . (12) Теперь рассмотрим стационарный случай, когда дополнительно предполагается, что пара матриц A, B управляема и существует управление в форме обратной связи, зависящей от x(k), которое осуществляет слежение за некоторым постоянным вектором z. Тогда матрица коэффициентов передачи экстраполятора также будет постоянной и может быть вычислена по формуле K = ANST (SNST + V)-1, (13) где матрица N определяется из решения матричного алгебраического уравнения m m N = (A - KS) N (A - KS )T AsNAj AszzT Aj + Q + KVKT. (14) s=1 s=1 Если решение уравнения (14) существует, N > 0, а также в силу того, что в этом случае матрица m m ZT ^ т т т ASNAS + > Aszz As + Q + KVK положительно определена, следует устойчивость матрицы динаs =1 s=1 мики экстраполятора (A - KS), так как будет выполнено условие устойчивости для дискретных систем в виде следующего матричного неравенства: (A - KS)N(A - KS)T - N < 0. (15) 3. Оценка неизвестного входа В качестве алгоритма оценивания неизвестного входа будем использовать алгоритмы МНК, в этом случае оценку можно построить на основе минимизации дополнительного критерия [7, 8]: (16) t=i где W, W - положительно определенные весовые матрицы, x(t) = A5c(t -1) + Bu(t -1) + r(t -1) . Построенные на основе минимизации (16) оценки неизвестного входа примут вид: rLsM(k) = [STWS+W]-1STW{y(k)-S (Ax(k-1) + Bu(k-1))}. (17) Оценка неизвестного входа, использующая дополнительное непараметрическое сглаживание [11-13], определим по формуле: Г №>(k) = [S TWS + W]-1S TWQ, (18) где j-я компонента вектора Q(k) имеет вид: * [y(i) - S(Ax(i -1) + Bu(i -1))] j fk - i + G I j) = ^----V ^ J. (19) Л 1 „ I k - i +1 I G i=1 Ij I Ij J В соотношении (19), G() является ядерной функцией и | j - коэффициентом сглаживания. Для построения оценки неизвестного входа можно также использовать скользящее среднее Q : Г SM> (k) = [S TWS + W]"1 ST WQ, (20) где j-я компонента вектора Q вычисляется по формуле (21) * y(i)-S(Ax(i -1) + Ba(i -1))|; i=k-l-1 l (в (21) l - ширина окна сглаживания). 4. Результаты моделирования Экстраполяции в дискретных системах с мультипликативными возмущениями реализована для трех вариантов вычисления оценок неизвестного вектора (формулы (17), (18), (20)). Оценки вектора состояния определялись из рекуррентного уравнения X(k +1) = Ax(k) + Bu(k) + r(k) + K(k)(y (k) - Sx(k)), X(0) = x0, (22) где матрица коэффициентов передачи K(k) вычислялась по формуле (12). Моделирование выполнено для следующих исходных данных: ( 0,85 0,1 ^ (1 0^ B = A = ( 0,01 -0,0015^ v0 1 у A = v 0,003 0,01 ( 0,1 0 ^ АЛ = АВ = 0,05 0,1 (1 0^ V 0 1 у (1 0^ C = , Q = v-0,05 0,94у ( 0,01 0,004^ , A2 = -0,002 0,01 (0,1 -0,03^ v 0 0,1 (0,03 0 Л ч 0 0,04,, (0,06 0 ^ S = V = (0,1 V0 1 у (0,1 0 ^ W= R = v 0 0,02у 0 ^ - (10 ^ W = 01 v 0 0,1у u(k) v 0 0,15 у (0,1 -0,1)T if k e [1;30], (0,2 0,1)T if k e [31; 75], (0 0)T if k e [76;126], (0,1 -0,1)T if k e [127;184], (0,2 0,1)T if k e [185; 200]. В (19) использовалась ядерная функция гауссовского вида: exp G (u) = -72л Результаты моделирования представлены на рис. 1, 2 и в табл. 1, 2. a b Рис. 1, a, b. Графики компонент вектора состояния x и оценок вектора состояния, полученные с использованием алгоритмов МНК х~LsM , модифицированного МНК с использованием непараметрического сглаживания х~, модифицированного МНК с использованием метода скользящего среднего x''sM Fig. 1, a, b. Graphs for the state vector x components and state vector estimates components using LSM algorithms x'LsM . modified LSM using nonparametric smoothing x'-Np , modified LSM using the moving average method x X (r (к) - Гг (к))2 В табл. 1 приведены результаты сравнения среднеквадратических ошибок отклонений оценок вектора состояния для алгоритмов МНК и модифицированного МНК, использующего процедуру скользящего среднего и алгоритм непараметрического сглаживания для оценки неизвестного входа. Расчет среднеквадратических ошибок оценивания выполнялся по формулам: X (x, (к) - x, (к))2 (i = 1,2). 1 g = N -1 N -1 Таблица 1 Среднеквадратические ошибки оценок вектора состояния ( axj) Номер компоненты i МНК МНК и скользящее среднее МНК и непараметрическое сглаживание 1 0,732 0,527 0,294 2 0,821 0,425 0,268 Рис. 2, a, b. Графики компонент вектора неизвестного входа r и оценок вектора неизвестного входа, полученных с использованием алгоритмов МНК , модифицированного МНК с использованием непараметрического сглаживания г''^, модифицированного МНК с использованием метода скользящего среднего г'SM Fig. 2, a, b. Graphs for the unknown input vector r components and estimates vector of the unknown input components using the LSM algorithms r, modified LSM using nonparametric smoothing r^, modified LSM using the moving average method r''sM> В табл. 2 показаны среднеквадратические ошибки отклонений оценок неизвестного входа. Таблица 2 Среднеквадратические ошибки оценок r(k) (arj) Номер компоненты i МНК МНК и скользящее среднее МНК и непараметрическое сглаживание 1 0,712 0,454 0,293 2 0,763 0,302 0,099 Из табл. 1 и 2 видно, что применение для оценок неизвестного входа МНК с дополнительным сглаживанием с помощью алгоритмов скользящего среднего и непараметрического сглаживания позволяет уменьшить среднеквадратические ошибки оценок экстраполяции вектора состояния и среднеквадратические ошибки оценок вектора неизвестного входа. Заключение Получено решение задачи синтеза экстраполятора для линейной дискретной модели с аддитивными и мультипликативными возмущениями, с неизвестным входом и неизвестными параметрами. Показано, что применение алгоритмов сглаживания при оценивании неизвестного входа позволяет повысить точность прогнозирования.

Ключевые слова

оценки экстраполяции, дискретная система, мультипликативные возмущения, неизвестный вход, extrapolation estimates, discrete system, multiplicative perturbations, unknown input, неизвестные параметры, unknown parameters

Авторы

ФИООрганизацияДополнительноE-mail
Ким Константин СтаниславовичТомский государственный университетаспирант кафедры прикладной математики Института прикладной математики и компьютерных наукkks93@rambler.ru
Смагин Валерий ИвановичТомский государственный университетпрофессор, доктор технических наук, профессор кафедры прикладной математики Института прикладной математики и компьютерных наукvsm@mail.tsu.ru
Всего: 2

Ссылки

Mohanty N.C., Soong T.T. Linear filtering in multiplicative noise // Information and Control. 1977. V. 34, is. 2. P. 141-147.
Stoica A.-M., Yaesh I. Kalman-type filtering for discrete-time stochastic systems with state-dependent noise // Proc. Mathematical Theory of Networkand Systems-MTNS. 28 July - 1 August, Blacksburg, VA. 2008. P. 1-6.
Wu Y., Zhang Q., Shen Z. Kalman filtering with multiplicative and additive noises // Proc. 12th World Congress on Intelligent Control and Automation (WCICA). 12-15 June, 2016. P. 483-487.
Germani A., Manes C., Palumbo P. Linear filtering for bilinear stochastic differential systems with unknown inputs // IEEE Trans. on Automatic Control. 2002. V. 47, No. 10. P. 1726-1730.
Yang F., Wang Z., Hung Y.S. Robust Kalman filtering for discrete time-varying uncertain systems with multiplicative noises // IEEE Trans. on Automatic Control. 2002. V. 47, No. 7. P. 1179-1183.
Smagin V.I. State estimation for nonstationary discrete systems with unknown input using compensations // Russian Physics Jour nal. 2015. V. 58, is. 7. P. 1010-1017.
Janczak D., Grishin Yu. State estimation of linear dynamic system with unknown input and uncertain observation using dynamic programming // Control and Cybernetics. 2006. No. 4. P. 851-862.
Gillijns S., Moor B. Unbiased minimum-variance input and state estimation for linear discrete-time systems // Automatica. 2007. V. 43. P. 111-116.
Hsien C.-S. On the optimality of two-stage Kalman filter for systems with unknown input // Asian J. of Control. 2010. V. 12, No. 4. P. 510-523.
Witczak M. Fault diagnosis and fault-tolerant control strategies for non-linear systems. Ch. 2: Unknown input observers and filters // Lecture Notes in Electrical Engineering. Springer International Publishing, 2014. P. 19-56.
Koshkin G.M., Smagin V.I. Filtering and prediction for discrete systems with unknown input using nonparametric algorithms // Proc. 10th International Conference on Digital Technologies. Zilina, Slovakia, July 9-11, 2014. P. 120-124.
Smagin V.I., Koshkin G.M. Kalman filtering and conrol algorithms for systems with unknown disturbances and parameters using nonparametric technique // Proc. 20th International Conference on Methods and Models in Automation and Robotics (MMAR 2015). 24-27 August, 2015, Miedzyzdroje, Poland. P. 247-251.
Koshkin G.M., Smagin V.I. Kalman filtering and forecasting algorithms with use of nonparametric functional estimators // Springer Proceeding in Mathematical Statistics / eds. R. Cao et al. 2016. Vol. 175. P. 75-84.
Athans M. The matrix minimum principle // Informat. and Contr. 1968. V. 11. P. 592-606.
 Экстраполяция в дискретных системах с мультипликативными возмущениями при неполной информации | Вестн. Том. гос. ун-та. Управление, вычислительная техника и информатика. 2019. № 47. DOI: 10.17223/19988605/47/6

Экстраполяция в дискретных системах с мультипликативными возмущениями при неполной информации | Вестн. Том. гос. ун-та. Управление, вычислительная техника и информатика. 2019. № 47. DOI: 10.17223/19988605/47/6