УДК 62.501
Вестник СибГАУ Том 17, № 4. С. 891-898
ОПРЕДЕЛЕНИЕ СТРУКТУРЫ ЛИНЕЙНОГО ДИНАМИЧЕСКОГО ОБЪЕКТА В ЗАДАЧАХ НЕПАРАМЕТРИЧЕСКОЙ ИДЕНТИФИКАЦИИ
А. В. Раскина
Сибирский федеральный университет Российская Федерация, 660041, г. Красноярск, просп. Свободный, 79 Е-mail: raskina.1012@gmail.com
Рассматривается проблема идентификации линейных динамических объектов, в частности, исследуется задача определения параметрической структуры динамического объекта с точностью до параметров при использовании непараметрических моделей. Линейные динамические процессы достаточно часто имеют место в различных контурах управления аэрокосмической техники, в частности, в процессе производства космических аппаратов. В данном случае локальный канал «вибростол - космический аппарат», определяемый вибросигналом и соответствующим сигналом датчика, установленным на аппарате, может описываться динамическими разностными уравнениями. Так как разностное уравнение динамического объекта состоит из запаздывающих на соответствующее число тактов выходных переменных, то задача сводится к задаче определения существенных переменных. Таким образом, метод определения структуры динамического разностного уравнения с точностью до параметров основан на применении правила выделения существенных переменных при непараметрической идентификации. Суть данного правила состоит в следующем: если для каждой переменной, входящей в непараметрическую оценку, поставить в соответствие оптимальный коэффициент размытости ядерной функции, при котором разница между полученной моделью и фактическим выходом объекта будет минимальна, то тогда та переменная, коэффициент при которой будет наименьшим, вносит больший вклад в формирование итоговой оценки и тем самым имеет большее влияние на выходную переменную. Рассматриваются непараметрические модели динамических объектов. Анализируется связь коэффициента размытости ядерной функции с влиянием конкретной переменной, учтенной в непараметрической модели на выход объекта. Приводится алгоритм восстановления структуры разностного уравнения динамического объекта, который включает в себя этапы нахождения оптимальных коэффициентов размытости ядерной функции для каждой переменной, отбор коэффициентов, исключение несущественных переменных, построение модели и расчет относительной ошибки моделирования. Алгоритм будет последовательно производить набор перечисленных выше действий, пока относительная ошибка моделирования не станет больше значения, полученного на предыдущей итерации. Подробно приводятся результаты численного исследования, проводимого методами статистического моделирования, которые показали эффективность предложенного метода как для численных аналогов дифференциальных уравнений, так и для разностных уравнений объектов с памятью.
Ключевые слова: разностное уравнение динамического объекта, выделение существенных переменных, объект с памятью, непараметрическая идентификация.
Sibirskii Gosudarstvennyi Aerokosmicheskii Universitet imeni Akademika M. F. Reshetneva. Vestnik Vol. 17, No. 4, P. 891-898
DETERMINATION OF THE STRUCTURE OF LINEAR DYNAMIC OBJECTS IN NONPARAMETRIC IDENTIFICATION PROBLEMS
A. V. Raskina
Siberian Federal University 79, Svobodny Av., Krasnoyarsk, 660041, Russian Federation E-mail: raskina.1012@gmail.com
The problem of identification of linear dynamic objects is considered. The problem of constructing the parametric structure of dynamic object up to the parameters by using non-parametric models is analyzed. Linear dynamic processes found in various aerospace control loops, for example, in vibration test spacecraft during their production. In this case, the local channel "vibrator - spacecraft" defined by the vibrator and the corresponding sensor signal installed on the spacecraft, can be described by the dynamic difference equations. Since the difference equation of the dynamic object consist of lagging on the appropriate number of cycles of the output variables, the problem reduces to the problem of determining the essential variables. Thus, the method of determining the structure of the dynamic differ-
entiaJ equation with up to parameters is based on the application of the rules of allocation of significant variables for nonparametric identification. The essence of this rule is as follows: if, for each variable, a member of the non-parametric estimation assign a optimum ratio blur kernel function in which the difference between the obtained model and the actual output of the object to be minimal, then, is the variable in which the rate will be lower, will contribute more to the formation of the final assessment, and therefore has a greater influence on the output variable. The article deals wih the non-parametric model of dynamic objects. The relationship of coefficient blurs kernel function and the influence of a particular variable, measured in non-parametric model output object is investigate. The algorithm of identification of the structure of the difference equation of the dynamic object includes the steps of finding the optimal coefficients blur kernel function for each variable sampling rates, elimination of the unimportant variables, modeling and simulation calculation of the relative error. The algorithm will produce consistently set the above steps until the relative errormodeling will not be greaterthan the value obtained in the previous iteration. Detailed results of numerical study conducted by the methods of statistical modeling of the effectiveness of the proposed method for numerical analogs of differential equations and difference equations for objects with memory are given.
Keywords: differential equation of dynamic object, the selection of essential variables, object with memory, non-parametric identification.
Введение. Одной из основных проблем современной теории идентификации и управления является определение структуры модели исследуемого объекта. На сегодняшний день большинство работ по данной проблеме посвящены методам параметрической идентификации [1-3], где структура процесса постулируется с точность до вектора параметров, исходя из имеющейся априорной информации. Причем выбор структуры модели процесса основывается на накопленном опыте работы с подобными объектами, а также зачастую осуществляется методом перебора различных моделей из заданного класса. Далее следует этап настройки параметров модели с использованием различных математических подходов, таких как, например, метод наименьших квадратов, различные рекуррентные оценки и т. д. В условиях недостатка априорной информации об объекте исследования такие подходы малоэффективны из-за наличия существующего многообразия процессов, их сложности и малоизученности.
В научной литературе по данной проблематике можно выделить несколько подходов. Рассмотрим наиболее известные из них. В работе [4] предлагается метод восстановления модели, в котором сначала следует этап выбора предполагаемой структуры из априорной информации, настройка параметров модели с использованием рекуррентных методов, а далее на основе имеющейся выборки измерений выполняется проверка предложенной модели. В случае получения неудовлетворительных результатов моделирования структура модели подлежит доработке, а параметры модели настраиваются заново. Как отмечалось ранее, подавляющее большинство подходов в задачах выбора параметрической структуры модели основаны на переборных методах, к которым относится, в частности, метод группового учета аргументов [5; 6]. Сущность метода состоит в выборе регрессионных моделей оптимальной сложности, где под сложностью понимается число параметров модели. Таким образом, в рамках данного метода осуществляется перебор многих моделей-претендентов с различным числом параметров по ряду критериев. В результате находится модель оптимальной структуры в виде
одного уравнения или системы уравнений [5]. В [7] излагаются методы структурной идентификации динамических систем на основе анализа информационных портретов. Предложена процедура нахождения порядка модели, основанная на минимизации ширины интервала изменения коэффициента структурности [7]. В настоящей статье предполагается новый метод определения структуры линейного динамического объекта, основанный на непараметрической теории идентификации [8-10]. В основу предлагаемого метода ложится правило выделения существенных переменных на основе измерения коэффициента размытости ядерной функции в непараметрической модели. Впервые данная идея была предложена профессором А. В. Медведевым в [11].
Постановка задачи идентификации. Пусть объект представляет собой одномерную динамическую систему и описывается следующим разностным уравнением:
X, = /.....Х(_к ,и,X (1)
где Д.) - неизвестный функционал; х( - выходная переменная процесса; и( - управляющее воздействие; к - глубина памяти динамического объекта (в терминологии А. А. Фельдбаума) [12]. Если проводить аналогию с описанием исследуемого процесса в непрерывном времени в виде дифференциальных уравнений, то к - порядок старшей производной в соответствующем уравнении. Здесь существенным является то, что вид функционала не определен с точностью до параметров.
Классическая для теории идентификации блок-схема моделирования данного процесса представлена на рис. 1.
На рис. 1 приняты следующие обозначения: х,, (г) -выход модели объекта; (г) - непрерывное время; индекс t - дискретное время; , Ь^ - случайные помехи измерений соответствующих переменных процесса; индекс Ь у переменных объекта из соображения простоты опущен: и( = и1, х( = х^; £ (г) - векторная случайная помеха.
Введем обозначения:
Ъ = (.....2к+\) = ((.....Х-к , Ч). (2)
тогда
Х = Г (ъ )•
(3)
При идентификации динамической системы (1) ее параметрическую модель естественно принять в форме
= 4 (
ч -1>
, Х
ч-к
, и(,а),
(4)
где а - вектор параметров, подлежащий оцениванию на основании обучающей выборки. Таким образом в случае линейной динамической системы определение структуры динамического объекта (1) сводится к определению переменных, входящих в состав модели (4). С учетом переобозначений (2) модель (4) может быть показана на схеме (рис. 2), которая иллюстрирует модель динамической системы в дискретном
времени, сведенную к модели статической системы, когда на вход последней поступают не только переменные и(, но и , х(-2.....х(-к и т. д.
На рис. 2 ЭЗ - элемент запаздывания [1]. Схема, изображенная на рис. 2, также является близкой к классической схеме, рассматриваемой в теории (рис. 1), так как дискретное уравнение, описывающее объект в данном случае, имеет соответствующий аналог среди дифференциальных уравнений. Как было замечено ранее [13; 14], рассмотрение динамического процесса с помощью данной схемы не является единственным, так как одна из переменных процесса, например х-2, может отсутствовать, а конечное выходное воздействие может зависеть только от переменных х-1, хч-3 и входного воздействия и: х( = /(хч-1, хч-3, и) . Таким образом, схема примет вид, изображенный на рис. 3.
Рис. 1. Классическая схема моделирования одномерного динамического объекта
Рис. 2. Блок-схема моделирования объекта с памятью
Рис. 3. Блок-схема моделирования объекта с памятью при фиксированных запаздывающих на соответствующие число тактов выходных переменных
В рассматриваемой задаче контроль входных-выходных переменных осуществляется через одинаковые интервалы времени Д^. Производя соответствующие измерения, получаем исходную выборку входных-выходных переменных {х;., щ, 1 = 1,5}, где 5 -объем выборки. Задача состоит в определении запаздывающих на определенные такты выходных переменных, которые необходимо учесть в модели (4), что позволит сделать предположение о параметрической структуре модели динамического объекта (4).
Выделение существенных переменных. В классическом случае, изображённом на рис. 2, задача идентификации состоит в оценивании класса операторов на основе выборки [х1, щ, 1 = 1,5}. Таким образом, в качестве оценки модели объекта может быть принято условное математическое ожидание вида
X (0 = Ы{х(0/ и (0}. (5)
В случае, когда динамический объект описывается дифференциальным уравнением при последовательной дискретизации, итоговое разностное уравнение будет содержать последовательно все переменные:
лгм, х(-2.....х(-к . Этому соответствует блок-схема
моделирования, изображенная на рис. 2. Тогда в качестве непараметрической модели объекта можно использовать модель, в которой все коэффициенты разностного уравнения будут учтены [11]:
1 5
х5 =-Х х •Н 51~1
^ и5 - и ^ к
1
п-7Н
м СХ
( - \ х5-] х1- 7
(6)
=1X XI • Н
^ и5 - и, ^
1=1
7 Н
(
сХ
хД- Н
сХ
(
\
\
(7)
В моделях (6), (7) Н (•) - ядерная колоколообраз-
ная функция, сЩ, сХ..... сК8 - коэффициенты размытости ядерной функции, которые удовлетворяют условиям сходимости [15]. В качестве колоколооб-разной функции Н(•) могут быть использованы различные ядра [15]. Следует учитывать, что модели вида (6), (7) можно использовать только при равных интервалах измерения Д^.
Оптимальные параметры размытости с 5 , с 5 .....
с 5 при наличии обучающей выборки находятся из задачи минимизации показателя соответствия выхода объекта х( и выхода модели Х5, основанного на методе скользящего экзамена, когда в моделях (6), (7) по индексу 1 исключается д-е наблюдение переменной, предъявляемой для экзамена:
к(с"5, с*:2,
* А
- хд)
д=1
(8)
Ш1П , а ф 1,
х 1 х х2 с хк н
Для случая, изображенного на рис. 3, модель может быть принята в виде (7), где учтены только выходные переменные, непосредственно от которых зависит выход процесса:
где индекс 1 фигурирует в формулах (6), (7).
Существенным в оценках (6), (7) является то, что в соответствие каждой выходной переменной х5-1, ..., хв-к, запаздывающей на некоторые величины,
ставится свой коэффициент размытости ядра
^ У
сз ..... с5 •
Из представленных выше моделей видно, что степень вклада той или иной выходной переменной из правой части уравнений в итоговое значение оценки зависит от
7
н
( _ \ Х5_] Х_]
сГ
(9)
\
Выражение (9) состоит из двух частей: —г
и н
( _ \ Х5_] ху_]
—г, то наблюдается следующая зависимость: чем
меньше с5 , тем больший вклад значение —] вносит
в итоговую оценку. Таким образом, можно построить следующую цепочку неравенств:
Х ¡2 ¡к 11 1
если сс < сс < ••• < сс , то —- >—- > ••• •
(10)
- 1 Н
сШ J
О(и)
сСш = 1,
(11)
где ш - некоторая переменная^ Исходя из данного
условия, при учете, что с1 = с 5 , с^ = с 52..... с1 =
* ¡к
= с 5 , можно построить следующую последовательность неравенств:
*Х *Хк
если с „ < с „ < ••• < с „ , то
Н
Л (
¡5-1 Х_1
< Н
Х5_ 2 Х_к
• •• < Н
( \
¡¡5 к _ Х_к
* ¡к
¡1
то Н ХБ_ 1 Х_1 <Н Х5_ 2 Х_ к
•х1 5,2
с V и 5 V с V и 5 V
выходная переменная х5_1 имеет большее влияние на выходную величину, чем х_2 • Графически это выражается так, как показано на рис-На рис 4
Н1 = Н
Х5_1 Х_1
Н 2 = Н
Х5 2 Х_к
•х2
Алгоритм вычисления значимых переменных х_] строится по следующей схеме^ Сначала задаем
начальное значение к. Затем строим модель по формуле (6) и считаем относительную ошибку моделирования Щ:
- =1 §< ^о2/§
Х)
• Что касается первого множителя
где тХ - математическое ожидание^
Далее на каждой У-й итерации выполняем следующий набор действий:
1 к
1 Для каждого коэффициента с5 ..... с5 находится оптимальное значение: с1 = с*^, = с*^.....
к *хк
сХ = с Х
2^ Находим из всех полученных значений максимальное - ] 5 •
Строим модель по формуле (6), исключая мно-
житель Н
( _ \ Х5_] Х(_]
, при учете, что ] - номер
Рассмотрим вторую составляющую уравнения (9) Коэффициенты с5 и ядерная функция Н (•) должны удовлетворять следующему свойству:
Приведем пример^ Пусть с 5 - коэффициент раз* Х2 -Г^
мытости ядра при х5_ 1, а с 5 - при х5_ 2 • Если
( \ ( \
„ „„ ^ ^ . - , и тогда
пРи стах_ 5 •
4^ Считаем относительную ошибку Щ •
Данные действия будут повторяться, пока Щ > Щ_ 1 •
Вычислительный эксперимент Численное исследование алгоритма непараметрической идентификации осуществлялось методами статистического моде-лирования^ Для вычислительного эксперимента был выбран объект, описываемый уравнением вида
0,5 С3Х + 0,5 + 0,8— + 1,5х(/) = (12)
Л2 Л2 Л
В дискретном виде данное уравнение равно Х = 2,69 • х_1 _2,46 х_ 2 + 0,73 х(_ 3 + 0,015^ • (13)
Зададим начальное значение к = 6^ Результаты моделирования с использованием формулы (6) представлены на рис 5^
Относительная ошибка моделирования Щ = 0,017^ Далее находим оптимальные коэффициенты с использованием алгоритма покоординатного спуска^ Найденные коэффициенты равны: с* = 0,675 при и(, с* = 0,131 при хм, с* = 0,186 при х(_2, с* = 0,412 при х_3, с* = 1,023 при Х(_4, с6 = 1,988 при х_5, с7 = 1,876 при х(_6 • Зависимость между полученными коэффициентами и значениями колоколообразной функции представлены на рис^ 6^
На риа 6
Н1 = Н
и -и.
\
(
, Н 2 = Н
(
5Х1
Н 7 = Н
После выполнения всех этапов алгоритма были исключены переменные х(-4, х(_5, х(-6. Таким образом, полученная итоговая структура динамического объекта, описываемого дифференциальным уравнением третьего порядка (12), равна:
ЛТ5 = а х(_1 +а2х(_2 + ах(_3 +р щ,
что соответствует структуре разностного уравнения объекта (13), а непараметрическая модель имеет вид
1 5 =1X X • Н
„X2
1=1
Н
и5 - и.
-Ф
Х5-1 Х1-1
„X2
„X3
Н
„X3
(14)
Результаты моделирования с использованием формулы (14) представлены на рис. 7.
Рис. 4. Зависимость величины значения ядерной колоколообразной функции от величины коэффициента размытости ядра
Рис. 5. Результаты моделирования процесса (13) при значении к = 6
Рис. 6. Зависимость величины значений ядерной колоколообразной функции от величин полученных оптимальных коэффициентов размытости ядра с^с*^, ..., с*^
О 13 30 3D « 50 60 70 80 &Э IQD
*
Рис. 7. Результаты моделирования процесса (13) при использовании модели (14)
Относительная ошибка моделирования равна
Ш3 = 0,005.
Рассмотрим процесс восстановления структуры модели для объекта с памятью. Пусть объект описывается следующим дискретным уравнением:
X = 0,2 • х(-1 + 0,4х,-3 + 0,5«,. (15)
Выберем в качестве начального значения к = 4. Относительная ошибка моделирования, полученная на первой итерации, равна Ш0 = 0,105. Далее следует этап нахождения оптимальных коэффициентов: с* = 0,124 при щ, с* = 0,231 при хм, с* = 0,393 при х,_2, с3 = 0,129 при х,_3, с5 = 0,425 при х,_4. После исключения переменных х^_2 и х(_4 относительная ошибка моделирования сократилась до Ш2 = = 0,012, а параметрическая структура была найдена в виде
X =а1 X-1 +а2 X-3 +$«!,
что соответствует структуре дискретного уравнения (15).
Заключение. В статье предложен алгоритм определения параметрической структуры линейного динамического объекта, основанный на применении правила выделения существенных переменных. Данный подход тесно связан с определением оптимальных коэффициентов размытости ядерной функции при использовании непараметрической модели функции регрессии по наблюдениям. В статье предлагается алгоритм, действие которого сводится к определению глубины памяти динамического объекта, которая представляет собой порядок старшей производной в соответствующем дифференциальном уравнении при проведении аналогии с описанием исследуемого процесса в непрерывном времени. Данный алгоритм состоит из нескольких этапов, включающих в себя определение оптимальных коэффициентов размытости, их отбор и построение итоговой модели. Одним из основных преимуществ предложенного алгоритма по сравнению с доминирующими на сегодняшней день методами восстановления параметрической структуры является тот факт, что разработанный
непараметрический алгоритм более применим к задачам практики, так как способен работать в условиях малой априорной информации об объекте.
Благодарности. Выражаю благодарность профессору А. В. Медведеву за его руководство и помощь в проведении исследований.
Acknowledgments. I would like to express my gratitude to Prof. A. V. Medvedev for continuous guidance and research assistance.
Библиографические ссылки
1. Цыпкин Я. З. Адаптация и обучение в автоматических системах. М. : Наука, 1968. 400 с.
2. Цыпкин Я. З. Информационная теория идентификации. М. : Наука. Физматлит, 1995. 336 с.
3. Эйкхофф П. Основы идентификации систем управления. М. : Мир, 1975. 683 с.
4. Бокс Д., Дженкинс Г. Анализ временных рядов, прогноз и управление. М. : Мир, 1974. Вып. 1. 406 с.
5. Ивахненко А. Г., Мюллер И. А. Самоорганизация прогнозирующих моделей. Киев : Техника, 1984. 350 с.
6. Стрижов В. В., Крымова Е. А. Методы выбора регрессионных моделей. М. : ВЦ РАН, 2010. 60 с.
7. Карабутов Н. Н. Наблюдаемые информационные портреты и задача структурной идентификации // Идентификация систем и задачи управления, SICPRO '07 : тр. VI Междунар. конф. (29 янв.-1 фев. 2007, г. Москва) М. : Институт проблем управления им. В. А. Трапезникова РАН, 2007. С. 89-115.
8. Медведев А. В. Основы теории адаптивных систем / СибГАУ. Красноярск, 2015. 525 с.
9. Медведев А. В. Теория непараметрических систем. Моделирование // Вестник СибГАУ. 2010. № 4 (30). С. 4-9.
10. Медведев А. В. Теория непараметрических систем. Процессы // Вестник СибГАУ. 2010. № 3 (29). С. 4-9.
11. Медведев А. В. Адаптация и обучение в условиях непараметрической неопределенности // Фундаментальные исследования (физико-математические
и технические науки). Новосибирск : Наука. Сиб. отд-ние, 1977. С. 92-97.
12. Фельдбаум А. А. Основы теории оптимальных автоматических систем. М. : Физматгиз, 1963. 552 с.
13. Банникова А. В., Медведев А. В. Об управлении объектами с памятью в условиях непараметрической неопределенности // Вестник СибГАУ. 2014. № 5(57). С. 26-37.
14. Банникова А. В., Сергеева Н. А. О непараметрическом моделировании стохастических объектов с памятью // Вестник СибГАУ. 2014. № 2 (54).
15. Надарая Э. А. Непараметрические оценки плотности вероятности и кривой регрессии. Тбилиси : Изд-во Тбил. ун-та, 1983. 194 с.
References
1. Cypkin Ja. Z. Adaptatsiya i obuchenie v avtoma-ticheskikh sistemakh [Adaptation and learning in automatic systems]. Moscow, Nauka Publ., 1968, 400 p.
2. Cypkin Ja. Z. Informatsionnaya teoriya identifi-katsii [Information theory of identification]. Moscow, Nauka, Fizmatlit Publ., 1995, 336 p.
3. Jejkhoff P. Osnovy identifikatsii sistem upravleniya [Identity-based control systems]. Moscow, Mir Publ., 1975, 683 p.
4. Boks D., Dzhenkins G. Analiz vremennyh rjadov, prognoz i upravlenie [Time series analysis, forecasting and management]. Iss. 1. Moscow, Mir Publ., 1974, 406 p.
5. Ivahnenko A. G., Mjuller I. A. Samoorganizatsiya prognoziruyushchikh modeley [Self-organization of predictive models]. Kiev, Tehnika Publ., 1984, 350 p.
6. Strizhov V. V., Krymova E. A. Metody vybora re-gressionnyh modeley [Methods of selection of regression models]. Moscow, VC RAN Publ., 2010, 60 p.
7. Karabutov N. N. [Observed information portraits and structural identification problem]. Trudy VI
Mezhdunarodnoy Konferentsii "IdentiSkatsiya sistem i zadachi upravleniya" SICPRO '07 [Proceedings of the VI International Conference "System Identification and Control Problems" SICPRO '07.] Moscow 29 January to January 2007. 2007, Moscow, Institut problem upravleniya im. V. A. Trapeznikova RAN Publ., P. 89-115 (In Russ.).
8. Medvedev A. V. Osnovy teorii adaptivnykh sistem [Basic theory of adaptive systems]. Krasnoyarsk, SibSAU Publ., 2015, 525 p.
9. Medvedev A. V. [The theory of nonparametric systems. Simulation]. Vestnik SibGAU. 2010, No. 4 (30), P. 4-9 (In Russ.).
10. Medvedev A. V. [The theory of non-parametric systems. Processes]. Vestnik SibGAU. 2010, No. 3 (29), P. 4-9 (In Russ.).
11. Medvedev A. V. [Adaptation and learning in a non-parametric uncertainty]. Fundamental'nye issledova-niya (ñziko-matematicheskie i tekhnicheskie nauki). Novosibirsk, Nauka Publ., 1977, P. 92-97 (In Russ.).
12. Fel'dbaum A. A. Osnovy teorii optimal'nykh av-tomaticheskikh system [Fundamentals of the theory of optimal automatic systems]. Moscow, Fizmatgiz Publ., 1963, 552 p.
13. Bannikova A. V. Medvedev A. V. [On objects with memory management in a non-parametric uncertainty]. Vestnik SibGAU. 2014, No. 5(57), P. 26-37 (In Russ.).
14. Bannikova A. V., Sergeeva N. A. [Nonparametric Stochastic modeling objects with memory]. Vestnik Sib-GAU 2014, No. 2 (54), P. 6-11 (In Russ.).
15. Nadaraja Je. A. Neparametricheskie otsenki plotnosti veroyatnosti i krivoy regressii [Nonparametric estimation of probability density and the regression curve]. Tbilisi, Tbil. un-t Publ., 1983, 194 p.
© Раскина А. В., 2016