Непараметрическая идентификация систем класса Винера

Тип работы:
Реферат
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 62−506. 1
НЕПАРАМЕТРИЧЕСКАЯ ИДЕНТИФИКАЦИЯ СИСТЕМ КЛАССА ВИНЕРА Коплярова Н. В.
Сибирский федеральный университет, Красноярск, Россия, (660 041, Красноярск, пр. Свободный, 79), e-mail:
koplyarovanv@mail. ru_
Рассматривается задача идентификации нелинейных динамических процессов, представимых в виде моделей Винера. Исследуются системы, находящиеся в условиях частичной параметризации. В данном случае структура и параметры уравнения, описывающего линейную динамическую часть объекта, неизвестны. Вид нелинейности предполагается известным с точностью до набора параметров. Приводится методика построения моделей для получении прогноза выхода нелинейных систем посредством сочетания моделей линейного динамического и нелинейного статического процессов в общей модели системы. Предлагается алгоритм, который не требует наличия полной априорной информации о структуре объекта. Полученная непараметрическая модель адекватно описывает исследуемые системы при различных видах нелинейной части объектов, в условиях помех в каналах измерения, при различном объеме выборки и виде входных воздействий.
Ключевые слова: непараметрическая идентификация, априорная информация, модель, нелинейная динамическая система, модель Винера.
NONPARAMETRIC IDENTIFICATION OF WIENER SYSTEMS Koplyarova N.V.
Siberian Federal University, Krasnoyarsk, Russia (660 041, Krasnoyarsk, 79 Svobodny Prospect), e-mail:
koplyarovanv@mail. ru_
The problem of nonlinear dynamical systems of Wiener type identification is considered. Systems in the conditions of partial parameterization are considered. In this case structure and parameters of equation, which describes the linear dynamical part of the object, are unknown. The common type of nonlinearity is assumed to be known with set of parameters. The proposed method of dynamic objects modeling is based on the nonparametric estimation of linear and nonlinear parts of the system. Presented algorithm don'-t requires complete a priory information about the object structure. This nonparametric model can adequate describe the investigation systems with different types of nonlinearity, in conditions of noise in the measure channels, and with different sample sizes and input actions.
Keywords: nonparametric identification, priory information, model, nonlinear dynamical system, Wiener model.
Введение. Современные производственные системы являются сложными техническими
объектами, часто имеющими нелинейную динамическую структуру. Для их исследования
необходимо создание математических моделей и алгоритмов, ориентированных на
применение современных средств вычислительной техники. В связи с этим актуальными
являются вопросы, связанные с разработкой методов идентификации нелинейных
динамических объектов. При этом необходимо рассмотрение всей имеющейся априорной
информации при построении моделей процесса, что достигается за счет применения
современных методов параметрической и непараметрической идентификации.
Наиболее важным с точки зрения приложений классом динамических процессов
являются системы типа «вход — выход», представимые в виде черного ящика, допускающие
активный эксперимент при отсутствии полной априорной информации о моделируемом
объекте. Настоящая работа посвящена рассмотрению задачи идентификации динамических
систем «в широком смысле», то есть в условиях, когда параметризация невозможна или
1
удается частично параметризовать модель исследуемого процесса на основании имеющейся априорной информации.
Рассматривается задача идентификации нелинейных динамических систем, представленных в виде последовательного соединения линейного динамического и нелинейного статического блоков (модель Винера), линейная динамическая часть которых находится в условиях непараметрической неопределенности. Вид нелинейности предполагается известным с точностью до параметров. Целью работы является повышение эффективности моделирования нелинейных динамических объектов в условиях малой априорной информации. Таким образом, исследуется задача моделирования нелинейных динамических процессов, находящихся в условиях частичной параметризации.
Постановка задачи. Рассматривается задача идентификации динамических систем класса Винера в условиях частичной параметризации. Пусть исследуемый объект может быть описан как последовательное сочетание линейных динамических и нелинейных статических блоков (модель Винера). Линейный динамический блок модели воспроизводит динамические свойства исследуемого объекта, нелинейный блок — имитирует его нелинейные свойства [5]. Сложность в построении моделей таких объектов заключается в том, что значения промежуточных сигналов (значения выхода блока, находящегося первым в последовательности) недоступны для измерения.
Задача идентификации нелинейной динамической системы поясняется на Рисунке 1:
Рисунок 1. Постановка задачи идентификации где Объект — система, состоящая из ЛЭ — линейной динамической и НЭ — нелинейной статической частей, ИУ — измерительное устройство, ы^) — входная переменная объекта,
X X
х (1) — выходная переменная, и, Х или (щхи Х=1,й} - наблюдения переменных процесса в
дискретный момент времени,) — ненаблюдаемое случайное воздействие, еи (/), ех (/)-случайные факторы, действующие в каналах измерения переменных в дискретные моменты времени 1-, Х (/) — выход модели объекта.
Исходные данные о состоянии исследуемого объекта представляют собой выборку
измерений реакции объекта на входное воздействие u (t): {м, xi, i = 1,5}. Параметры и порядок дифференциального уравнения, которым может быть описана линейная динамическая часть системы, неизвестны. Пусть нелинейность в объекте описывается функцией, вид которой предполагается известным с точностью до набора параметров. Рассмотрим системы, в которых нелинейный элемент представляет собой квадратор: /(w) = а ¦ 2,
w, w? Ь1
или звено насыщения (с порогом насыщения Ь_/): /(w) = & lt- Ь, w & gt- Ь1.
— Ь, w & lt- -Ь1
Требуется на основании имеющейся информации построить модель данной системы, адекватно описывающую ее поведение при подаче на ее вход воздействий различного вида.
Непараметрическая модель системы типа Винера. Будем исследовать системы, которые могут быть представлены в виде модели Винера [2] (Рисунок 1). Пусть порядок и параметры линейной динамической части модели Винера неизвестны, а структура нелинейного элемента задана с точностью до набора параметров а. Основная идея алгоритма, предлагаемого для построения моделей таких систем, заключается в использовании непараметрических оценок для описания связей системы, информация о которых по каким-то причинам неизвестна (в данном случае измерению недоступна линейная динамическая часть системы), а также в применении непараметрического подхода к оцениванию функций.
Как видно из рисунка, значение выхода исследуемого объекта вычисляется как некоторая функция от значения w (t):
X (г) = / ^(г), а}, (1)
где х (г) — выходной сигнал системы- w (t) — выход линейной части системы (не измеряемый) — м (1) — входной сигнал системы-), а} - нелинейный оператор.
Предположим, что выход w (t) может быть измерен, тогда модель линейного динамического звена системы при входном воздействии м (1) может быть оценена как интеграл Дюамеля, то есть:
г
w (г) = | к'-(г -т)м (т)с1т, (2)
о
где к (г) — переходная характеристика системы, т — переменная интегрирования. Тогда в соответствии с (1) модель исследуемого объекта:
Л л
х (|) = / ] И'-(I-т)и (т)Ст, ас (3)
V 0)
Математическая модель нелинейного объекта может быть представлена в виде системы (3), где вместо переходной функции И (1) и параметров, а используются их статистические оценки.
Значение переходной функции И (1) представляет собой реакцию линейного динамического элемента системы на подачу на ее вход воздействия и (1)=1, то есть И (1)=м (1/и (1)=1). При этом значение н'-(1) недоступно для измерения. После подачи на вход единичного воздействия возможно измерить только выход нелинейного объекта х (1), который будет иметь значение, равное х (1)=/(И ($.
Тогда значения переходной функции линейного элемента системы в дискретные моменты времени могут быть оценены следующим образом:
И = /-1(Х1., а), (4)
где х1г- - экспериментально полученные значения выхода исследуемого объекта при подаче на вход единичного ступенчатого воздействия 1(1), И1г- - рассчитанное значение переходной характеристики линейного элемента системы. Алгоритмы оценивания неизвестной функции /-1(х) и ее параметров, а зависят от типа нелинейности в объекте и будут приведены далее.
Далее на основании выборки дискретных значений можно оценить переходную функцию системы в виде стохастической аппроксимации регрессии непараметрического типа следующим образом [4]:
1 s
h (t) =--1 h, H
sc" i=i
ft — t4
V cs У
(5)
При этом колоколообразная функция H (•) и параметр размытости cs должны удовлетворять условиям сходимости [4]:
0- lim sc «= ?
cs & gt- 0- lim cs
S ®? s ® ?
J H u) du = 0
W (u) W (u)
cs J H'- (u)udu = -1
(6) (7)
lim
c ~lH
s ®? s 11
f л t-1
V cs
= d (t-1), u =
t-1
В этом случае непараметрическая оценка весовой функции, которая является производной по времени от переходной, примет следующий вид:
t
c
s
1 *
к (г)=К (0=--1 Кн'-
ЭС, ,=1
'-г -г & gt-
V С у
Подставив оценку весовой функции в интеграл Дюамеля в соответствии с (2), получим непараметрическую модель, оценивающую выход линейного элемента [3]:

1 э г/Ат
1 '-I I Кн'-
'-г -т- О
ЭС* ,=1,=1
и (т,)Ат
(9)
где т — переменная интегрирования, Ат — шаг дискретизации.
Тогда непараметрическая модель нелинейного объекта примет вид [7]:
х
(г) = / | К'-(г — т) и (т)ёт, а = /
V о у
1
Ч А
• IIКН'-
(г — т, — г, & gt-
v эсэ ,=1 ,=1
и (т,)Ат, а
(10)
где х (г) — выходной сигнал системы- и (г) — входной сигнал системы- /{м (г), а} - оценка нелинейной функции, а — оценки параметров нелинейного элемента системы, К (г) — оценка переходной функции линейного элемента системы.
Идентификация нелинейной системы с квадратором. Пусть имеем систему, представленную в виде модели Винера (рисунок 2). Причем нелинейная часть системы
представляет собой квадратор, описываемый функцией вида: /(м) = ам& gt-2, где а=еопБХ. Выход объекта вычисляется следующим образом: х (г) = /(м, а) = ам2. При единичном входном воздействии и (г) = 1(г) выход нелинейной системы равен х1(г) = х (г/и (г) = 1) = ам (г)2. То есть оценку значений переходной характеристики линейного элемента К (г) = м (г) можно рассчитать через выход исследуемого процесса следующим образом: Кг, = ^хх~[а, где х1 — выход исследуемой системы при единичном
входном воздействии в момент времени X [2].
Далее оцениваем переходную характеристику системы с применением непараметрического метода следующим образом:
1 Э «
К (г) = - IКН
ЭС& lt- ,=1
г г — г ^
V Сэ у
О-
ЭСЭ ,=1 а
гг — г, ^
V сэ у
(12)
При произвольном входном воздействии и нулевых начальных условиях выход линейной части системы описывается интегралом Дюамеля (8). С учетом рассчитанного значения переходной функции (16) выход линейного элемента равен:
1 [ХГ л
м. (г)=-• Ц^н'-
ЭСЭ ,=1 ,=л а
г-т- г,

V сэ
и (т,)Ат,
где, а — неизвестная константа (параметр нелинейного элемента системы).
5
С
Э
т
С
Э
Таким образом, модель нелинейного динамического объекта типа Винера имеет вид:
х. (У) = / (М'-(У), а) =) = а
ъЪ^'-
г=1 у=м а
'-У-ъ — У & gt-
V
и (тг)Аъ
(14)
Расчетная формула для оценки выхода нелинейной системы при произвольном входном воздействии примет следующий вид:
X. (У):
. УАъ _
ЪЪл/Х^я'-
1=1 ]=1
V
и (т,)Аъ
(15)
где х1, — реакция нелинейной системы на единичное входное воздействие, и (У) — входное воздействие. Таким образом, получается методика, для реализации которой не требуется точное численное значение параметра квадратора.
Пример. Рассмотрим нелинейную динамическую систему, состоящую из квадратора (параметр а=2) и разностного аналога дифференциального уравнения (имитирующего объект): 3 • х& quot- (у) + 1,3 • X (у) + 1 • х (у) = и (у). Предполагаем, что вид и параметры реального уравнения, описывающего процесс, неизвестны. Результат моделирования такого процесса представлен на рисунке 3.
Рисунок 3. Результат оценки выхода х (1). хшоёе1(1:) — модель нелинейной системы, х -выход системы, объем выборки б=250, Ь=0. 16, помеха 5%, входное воздействие: и (1:)=3со8(0. 51:)+8т (0. 21:), относительная ошибка моделирования 3.7%
Анализируя результаты моделирования, можно сказать, что предлагаемые алгоритмы позволяют достаточно точно оценивать значения выхода исследуемой системы при различных значениях параметров нелинейной части объекта, в условиях зашумленности каналов связи, при различных входных воздействиях.
Идентификация нелинейной системы с насыщением. Пусть нелинейное звено системы описывается функцией вида (2). В данном случае при н'-{У)& lt-Ь выход объекта совпадает с выходом его линейной динамической части. В остальных случаях выход объекта представляет собой константу, которую возможно определить опытным путем в результате
2
2
нескольких статических экспериментов, описанных далее [1].
1. Проводится серия экспериментов, в ходе которых будем подавать на вход системы воздействия различной амплитуды ui j = Cj, Cj = const, в результате получим выборку
{ui j, xi j}, i = 1, s, j = 1, m, на основании которой делается вывод о нелинейности.
2. В результате получается оценка нелинейного элемента системы, параметры которого определяются следующим образом:
— находится расстояние между двумя соседними измерениями: Ah- ¦ =
i, j
h — шаг дискретизации-
— оценка порога нелинейности b = x- f, если hh, & lt-?,?>- 0-
i, j j
— b = M{bj }, где Xi, j = bb и xj = y_!, j.
3. На вход объекта подается ступенчатая функция, амплитуда которой не превышает значение порога b — в результате получаем переходную характеристику и строим модель линейной части объекта в виде интеграла Дюамеля.
4. Строится модель объекта, выход которой равен значению функции, описывающей нелинейное звено, аргументом которой является выход модели линейной части объекта.
Рассмотрим нелинейную динамическую систему, состоящую из звена насыщения вида (2) (параметры b=0. 75, b 1=1) и разностного аналога дифференциального уравнения (имитирующего объект): 3 • x& quot- (t) + 0,57 • X (t) +1 • x (t) = u (t).
— x (t) -xmodel (t)

7| j i I
i 3
0 10 20
Рисунок 4. Результаты оценки выхода х (1), хшоёе1(1-) — модель нелинейной системы, х (1-) — выход системы, объем выборки 8=200, шаг дискретизации Ь=0. 15, помеха 5%, входное воздействие и (1-)=0. 358 т (1-), средняя ошибка моделирования 2%
Делая анализ работы модели нелинейного динамического объекта с видом нелинейности типа звено насыщения и квадратор, можно сделать выводы о том, что непараметрическая модель достаточно точно описывает систему при различных значениях параметров нелинейной части объекта, в условиях зашумленности каналов связи, при различном объеме выборки и различных входных воздействиях.
xi, j xi1, j
h, где
Заключение. В статье рассмотрена задача идентификации нелинейных динамических систем, представимых в виде моделей Винера. Исследуются системы, находящиеся в условиях частичной параметризации. В данном случае структура линейного динамического блока неизвестна, а вид нелинейности предполагается известным с точностью до параметров.
Задача идентификации нелинейной системы рассмотренного типа может быть разделена на две части. Сначала рассматривается непараметрическая идентификация линейного элемента, алгоритм которой связан с тем, что реакция линейной системы на входное воздействие описывается интегралом Дюамеля. Суть метода заключается в построении моделей и получении прогноза выхода нелинейных динамических систем посредством сочетания моделей линейного динамического и нелинейного статического процессов в общей модели системы. Данный метод не предусматривает наличия полной априорной информации о структуре объекта.
Показана работоспособность предлагаемого алгоритма и применимость к задаче построения нелинейных динамических систем типа Винера в различных условиях, что расширяет область применения непараметрических методов к новым объектам исследования, а также закладывают теоретическую основу для дальнейшего проектирования. Описанные методы построения и исследования моделей позволяют производить моделирование технологических процессов и систем, имеющих нелинейную структуру.
Список литературы
1. Коплярова Н. В., Сергеева Н. А. Непараметрические алгоритмы идентификации систем класса Винера и Гаммерштейна. Системы управления и информационные технологии. -Воронеж, 2013. — Вып. 2.1 (52).
2. Медведев А. В. Непараметрические алгоритмы идентификации нелинейных динамических систем: сб. Стохастические системы управления. — Новосибирск: Изд-во «Наука», 1979. — С. 15−22.
3. Медведев, А. В. Теория непараметрических систем. Моделирование / А. В. Медведев // Вестн. Сиб. гос. аэрокосмич. ун-та им. акад. М. Ф. Решетнева. — Красноярск, 2010. — Вып. 4. -С. 4−9.
4. Надарая Э. А. Непараметрическое оценивание плотности вероятностей и кривой регрессии / Э. А. Надарая. — Тбилиси: Изд-во Тбилисского ун-та, 1983.
5. Попков Ю. С. Идентификация и оптимизация нелинейных стохастических систем. — М.: Энергия, 1976. — 440 с.
6. Цыпкин Я З. Основы информационной теории идентификации. — М.: Наука, 1984. — 320 с.
7. Чайка С. Н. К идентификации динамических систем при частично параметризованной структуре модели / Динамика систем: Управление и оптимизация. — Горький: Изд-во Горьковского гос. ун-та, 1989.
Рецензенты:
Охорзин В. А., д.т.н., профессор, профессор кафедры Прикладной математики Сибирского государственного аэрокосмического университета «Сибирский государственный аэрокосмический университет имени академика М.Ф. Решетнева», г. Красноярск. Демиденко Н. Д., д.т.н., профессор, ведущий научный сотрудник СКТБ «Наука» КНЦ СО РАН, Федеральное государственное бюджетное учреждение науки Специальное конструкторско-технологическое бюро «Наука» Красноярского научного центра Сибирского отделения Российской академии наук, г. Красноярск.

ПоказатьСвернуть
Заполнить форму текущей работой