Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

ADAPTIVE ALGORITHM FOR IDENTIFICATION OF NONLINEAR SYSTEMS BY VOLTERRA SERIES

Tsibizova T.Yu. 1
1 Bauman Moscow State Technical University (National Research University)
This article is devoted to identification of nonlinear systems of automatic control. A constructive approach to solving this problem is the use of filter structure as a sequence of Volterra. Describes the implementation of the Volterra filter of 2nd order, its mathematical description, are defined by linear and quadratic weight filter that minimizes the mean square error. The proposed adaptive algorithm for implementation of the Volterra filter of 2nd order in which the operators are linear and the quadratic filter can be adjusted as an extension of the algorithm of least mean square. On the basis of the mean square error defined by its asymptotic remainder in the adaptive implementation. It is concluded that the application of the adaptive algorithm can give an effective, flexible approaches to development and implementation of the Volterra filter of 2nd order. On the other hand, the distribution of the results on the Volterra filter of higher order is an interesting area for further research.
identification of nonlinear systems
adaptive algorithm
mean square error
linear filter
quadratic filter
optimal filter
1. Voronov E.M., Karpunin A.A., Vanin A.V. Optimizatsiya upravleniya strukturno slozhnymi sistemami, Inzhenernyy zhurnal: nauka i innovatsii, 2013, no. 10 (22), URL http://engjournal.ru/articles/1080/1080.pdf. (accessed 10.09.2016).
2. Zadorozhnaya N.M. Adaptivnye sistemy avtomaticheskogo upravleniya s dvumya etalonnymi modelyami. Sovremennye problemy nauki i obrazovaniya, 2015, no. 1–2; URL: http://science-education.ru/ru/article/view id=20134.
3. Kobzev G.K. K voprosu ob identifikatsii yader Volterra pri modelirovanii nelineynykh dinamicheskikh system. Vestnik Irkutskogo gos. tekhn. Universiteta, 2006, no. 4–1 (28), pp. 80–83.
4. Lukyanova N.V., Kuznetsov I.A. Identifikatsiya nelineynykh dinamicheskikh sistem na osnove razlozheniya funktsionalov metodom Vinera. Materialy konferentsii «Upravlenie v morskikh i aerokosmicheskikh sistemakh» (Proc. of the conference «Management in the marine and aerospace systems»). Saint-Petersburg, 2014, pp. 633–636.
5. Lukyanova N.V. Modulnyy metod modelirovaniya s ispolzovaniem razlozheniya Vinera. Avtomatizatsiya. Sovremennye tekhnologii, 2015, no 9, pp. 17–22.
6. Pavlenko S.V., Polozhaenko S.A. Optimizatsiya vychislitelnykh algoritmov approksimatsionnogo metoda identifikatsii nelineynykh system v vide modeley Volterra. Informatika i matematicheskie metody v modelirovanii, 2013, Vol. 3, no. 2, pp. 103–112.
7. Proletarskiy A.V., Neusypin K.A., Kuznetsov I.A. Razrabotka kriteriya stepeni identifitsiruemosti parametrov dinamicheskikh system, Trudy FGUP NPTSAP, Sistemy i pribory upravleniya, 2014, no. 4 (30), pp. 87–93.
8. Pupkov K.A., Tsibizova T.Yu. Realizatsiya filtra Volterra vtorogo poryadka dlya identifikatsii nelineynykh system upravleniya, Nauka i obrazovanie, 2006, no. 6, URL http://technomag.edu.ru/doc/58741.html. (accessed 10.09.2016).
9. Tsibizova T.Yu. Identifikatsiya nelineynykh sistem avtomaticheskogo upravleniya pri pomoshchi filtrov Volterra. Fundamentalnye issledovaniya, 2015, no. 2–14, pp. 3070–3074.
10. Tsibizova T.Yu. Metody identifikatsii nelineynykh sistem upravleniya. Sovremennye problemy nauki i obrazovaniya, 2015, no. 1–1; URL: http://science-education.ru/ru/article/view id=17910.
11. Tsibizova T.Yu., Chan N.Kh. Sposoby realizatsii protsedury identifikatsii na osnove filtra Volterra. Avtomatizatsiya. Sovremennye tekhnologii, 2015, no. 8, pp. 31–34.

Решение задач управления динамическими объектами предполагает использование математической модели исследуемого процесса. Математические модели, полученные на основе физических или каких-либо других законов, в практических приложениях, как правило, не всегда точно отражают исследуемые процессы. Поэтому для уточнения структуры и параметров математической модели применяются различные алгоритмы идентификации и алгоритмы построения моделей [1, 2, 4, 5]. Алгоритмы построения моделей позволяют получить высокоточные математические модели исследуемых объектов, однако их использование для управления затруднительно, так как структура модели заранее неизвестна и вопрос проверки устойчивости требует дополнительных исследований. Алгоритмы идентификации позволяют определить отдельные параметры матрицы модели, структура которой задана априори. Для решения задачи идентификации нелинейных объектов разработано довольно много подходов и методов [3, 4, 6, 7]. На современном этапе возросли требования к точностным характеристикам применяемых алгоритмов идентификации. В связи с этим модифицируются классические подходы к решению задачи идентификации нелинейных систем с целью повышения их точности и уменьшения ограничений применения [1, 6, 10], а также универсальные поисковые методы, которые требуют минимальной априорной информации об идентифицируемой системе, но сложны в реализации.

В настоящей работе исследована проблема идентификации нелинейных систем. В качестве основной трудности данной проблемы можно назвать необходимость обработки большой базы данных, характеризующей работу идентифицируемой системы. Конструктивным подходом в решении данной задачи является использование фильтрующей структуры в виде последовательности Вольтерра [3, 6, 8]. Однако одной из главных причин достаточно редкого применения методики фильтрации Вольтерра на практике является значительная сложность, связанная с реализацией фильтров Вольтерра [9, 11]. Таким образом, главной задачей является нахождение упрощений в разработке и реализации фильтра Вольтерра.

Фильтр Вольтерра 2-го порядка

Возьмем фильтр Вольтерра 2-го порядка (ФВ2), который состоит из параллельной комбинации линейного и квадратичного фильтров [8, 9]:

Tsibizova01.wmf (1)

где {a(j)} и {b(j, k)} называются линейным и квадратичным весом соответственно, а N указывает длину фильтра (предполагается симметричность квадратичных весов фильтра, т.е. b(j,k) = b(k,j)).

Далее возьмем, что x(n) и s(n) – это случайные процессы с нулевым математическим ожиданием с дискретным параметром n. Нужно найти веса фильтра, которые минимизируют среднюю квадратичную ошибку (СКОШ) между s(n) и выходом фильтра y(n), т.е.

Tsibizova02.wmf (2)

где предполагается строгая стационарность s(n) и x(n) с нулевым математическим ожиданием.

Первым шагом в определении минимума среднеквадратичной ошибки является требование бездрейфового выхода фильтра. Другими словами, должно быть E[y(n)] = 0, так как основной сигнал имеет нулевое математическое ожидание.

Тогда имеем следующее соотношение между h0 и b(j, k):

Tsibizova03.wmf (3)

где Tsibizova04.wmf обозначает автокорреляционную функцию x(n). Важно включение члена нулевого порядка h0. Некоторые из предыдущих исследователей не имели выхода нулевого порядка, но без этого выход минимальной средней квадратичной ошибки фильтра Вольтерра не является обязательно бездрейфовым и ошибка имеет, следовательно, тенденцию к увеличению, в отличие от фильтра Вольтерра 2-го порядка с h0, выраженном в (3).

Подставляя в (1) выражение (3), получим формулу для определения фильтра Вольтера 2-го порядка:

Tsibizova05.wmf (4)

Следующий шаг – определение линейных и квадратичных весов фильтра a(j) и b(j, k), которые определяют минимум среднеквадратичной ошибки.

Выведем простое решение для оптимального фильтра Вольтерра 2-го порядка в предположении, что на входе фильтра гауссион. Во-первых, заметим, что (4) может быть переписано в виде [8, 9]:

Tsibizova06.wmf (5)

где tr – след матрицы В;

Tsibizova07.wmf

Tsibizova08.wmf

Tsibizova09.wmf

Rx указывает на N×N матрицу от x(n), где rx(j, k) = rx(j – k) – автокорреляционная функция входного сигнала x(n). A и B – операторы линейного и квадратичного фильтра соответственно.

Перед выводом решения определим кросс-корреляционную rsx(j) и кросс-бикорреляционную tsx(j, k) функции между x(n) и s(n) следующим образом:

Tsibizova10.wmf

Tsibizova11.wmf

Поскольку предполагается, что s(n) и x(n) – строго стационарны, то как rsx(j), так и tsx(j, k) являются независимыми от переменной n. Кросс-бикорреляционная функция tsx(j, k) определяет статистическую зависимость между s(n) и x(n), которая является критичной при нахождении оптимального квадратичного оператора фильтра. Кроме того, надо сказать, что кросс-бикорреляционная функция является симметричной, т.е. tsx(j, k) = tsx(k, j).

В матричной форме кросс-корреляционную и кросс-бикорреляционную функцию можно записать следующим образом:

Rsx = [rsx(0), …, rsx(N – 1)]T;

Tsibizova12.wmf

Отсюда линейный и квадратичный операторы фильтра Вольтерра 2-го порядка с минимальной среднеквадратичной ошибкой определяются следующим образом:

Tsibizova13.wmf (6)

Tsibizova14.wmf (7)

Адаптивный алгоритм реализации

Алгоритм наименьшего среднего квадратического (НСК) для линейного фильтра хорошо известен и представляется в виде

A(n + 1) = A(n) – 2μAe(n)X(n). (8)

Здесь A(n) – оператор линейного фильтра в момент n и e(n) = AT(n)X(n) – s(n) – остаточная ошибка фильтра. Также μA – положительная константа, от которой зависит устойчивость и сходимость алгоритма.

Рассматривая адаптивную реализацию ФВ2, можно заметить, что оператор линейного фильтра может быть скорректирован, используя алгоритм (8), за исключением остаточной ошибки, представляемой в виде:

Tsibizova15.wmf (9)

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

Однако в (9) следует отметить, что выражение нулевого порядка tr{B(n)Rx} не является константой при адаптивной реализации. Вспомним, что выражение нулевого порядка необходимо для вычитания математического ожидания выходного значения квадратичного фильтра из выходного значения ФВ2 [9]. Следовательно, когда Rx неизвестно, рекурсивный алгоритм оценки (например, НЧ-фильтр) для уровня среднего выходного значения квадратичного фильтра можно заменить выражением нулевого порядка.

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

Tsibizova16.wmf (10)

где μB – положительная константа.

Концептуально данный алгоритм имеет под собой основание, т.к.

Tsibizova17.wmf

и его реализация вполне приемлема.

Первая аппроксимация должна допускать, что A(n) и B(n) – независимы от пары {s(n), X(n)}. Тогда мы имеем

Tsibizova18.wmf

Таким образом, если мы определим δB(n) = B(n) – B0, где B0 – оптимум оператора квадратичного фильтра, определяемого с минимальной среднеквадратичной ошибкой следующим образом: Tsibizova19.wmf то математическое ожидание, взятое от обеих частей (10), даст

Tsibizova20.wmf (11)

Заметим, что математическое ожидание оператора квадратичного фильтра сходится к оптимуму B0 тогда и только тогда, когда Tsibizova21.wmf при n → ∞.

Действительно, если Tsibizova22.wmf для некоторого n, то это остается справедливым и для всех других n. Для Tsibizova23.wmf Tsibizova24.wmf если размер шага μB выбирается следующим образом:

Tsibizova25.wmf (12)

где λmax – наибольшее собственное значение матрицы R.

Отсюда средняя ошибка между B(n) и B0 монотонно снижается до 0 с увеличением времени. Далее при детальном рассмотрении процесса адаптации можно увидеть, что каждый элемент матрицы E[δB(n)] состоит из набора компонент вида e–αt, в которых самое быстрое и самое медленное затухание определяется λmax и λmin.

В адаптивной реализации ФВ2 флуктуации операторов линейного и квадратичного фильтров добавляют некоторую дополнительную СКОШ в выходное значение фильтра даже при устойчивом положении адаптационного процесса. Таким образом, асимптотическая СКОШ адаптивного ФВ2 в общем случае больше, чем СКОШ оптимального ФВ2. Для оценки остатка СКОШ запишем СКОШ адаптивного ФВ2 в виде

ξ(n) = ξopt + ξA(n) + ξB(n);

причем остаток СКОШ линейного и квадратичного представляется:

Tsibizova26.wmf

Tsibizova27.wmf

где δA(n) = A(n) – A0 и δB(n) = B(n) – B0 – суть отклонения A(n) и B(n) от их оптимальных значений.

Затем, получим следующие ограничения для остатка СКОШ:

Tsibizova28.wmf

Tsibizova29.wmf

Так как нас интересуют главным образом асимптотические значения остатка СКОШ после достижения адаптивным процессом своего устойчивого состояния, то мы допускаем, что A(n) и B(n) близки к A0 и B0 соответственно. Отсюда выбираем μA и μB так, чтобы Tsibizova30.wmf и Tsibizova31.wmf

Таким образом, имеем

δA(n) = –2μAe(n)X(n)

и

δB(n) = –μBe(n)X(n)XT(n).

Кроме того, при данных условиях как X(n), так и X(n)XT(n) не коррелирует с остаточной ошибкой e(n).

Следовательно, получаем следующие ограничения для асимптотического остатка СКОШ (АОСКОШ):

ξA ≤ 4μAξoptNrx(0);

ξB ≤ 3μBξopt[Nrx(0)]2.

Заключение

В работе представлены способы реализации фильтра Вольтерра 2-го порядка, связанные с упрощением его реализации. Показано, что простое решение для оптимального ФВ2, его адаптивный алгоритм могут дать эффективные гибкие подходы для разработки и реализации ФВ2. С другой стороны, распространение полученных результатов на фильтр Вольтерра более высокого порядка является интересным объектом дальнейшего исследования.