Моделирование сложных динамических объектов Текст научной статьи по специальности «Строительство и архитектура»

Д.Д. Мищенко

МОДЕЛИРОВАНИЕ СЛОЖНЫХ ДИНАМИЧЕСКИХ ОБЪЕКТОВ

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

Ключевые слова: модель, динамическая система, непараметрическая оценка.

D.D. Mishchenko

MODELING THE COMPLEX DYNAMIC OBJECTS

The development of accurate and adequate model of the object in the conditions of non-parametric uncertainty, the task of techniques description that can be applied to solve a particular type of identification tasks are considered in the article.

Key words: model, dynamic system, non-parametric assessment.

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

Цель и задачи исследований. Построение точной и адекватной модели объекта в определенных условиях, описание методики, которая может быть применима для решения различного типа задач с небольшими издержками ресурсов. Допустим, что в ходе изучения объекта пользователь владеет информацией о том, что он относится к классу детерминированных, стационарных, линейных, динамических объектов и объектов с запаздыванием, то есть «качественной» информацией о свойствах поведения объекта. Допускается, что пользователь способен воздействовать на объект, контролируя уровень своего воздействия. Так, он может компенсировать недостаток информации об объекте с удовлетворительной корректностью. Возьмем во внимание, что в силу ограниченности средств изучения, математического аппарата, используемого пользователем для описания объекта, не все факторы, влияющие на объект, будут учтены в модели. Как следствие, в ней могут появиться случайные помехи, причем информации о характере помехи нет: закон распределения случайных помех, действующих при взаимодействии объекта с окружающей средой, неизвестен [1]. Иначе говоря, любая другая информация об объекте и его моделях недоступна для пользователя, либо не удовлетворяет поставленным задачам. В таких рамках пользователь может потратить большое количество ресурсов (вычислительных, материальных, производственных, временных, научных, финансовых), прежде чем найти решение. Выбор пути в поисках решения предполагает необходимость вложения различного количества ресурсов для решения задачи. Между такими ресурсами, как, например, точность и время (быстрота нахождения решения), обычно существует противоречие: для повышения точности уменьшают шаг интегрирования, но при этом увеличивается время счёта [4].

Методика и результаты исследований. Для решения поставленной цели предлагается взять за основу модели свертки характеристик объекта характеристики, воздействующие на объект, и характеристики, которые описывают реакции объекта на такие воздействия, то есть модели в виде интеграла Дюамеля [3] и интеграла Коши-Лагранжа. Математическая запись интеграла Дюамеля приведена в формуле (1), интеграла Коши-Лагранжа - в формуле (2):

і

х() = І И'( — т)л {т)Лт

(1)

х() = / { — ї0 )+| к'{і — т)и{т)Лт,

(2)

где х(І) - выходная реакция объекта, который наблюдается в течение времени 1, начиная с момента времени Ь; И'^) - импульсно-переходная характеристика объекта, связанная с переходной характеристикой объекта Іі^) оператором дифференцирования ^)/Л = іф); и(І) - входное воздействие на объект, которое приводит к изменению его реакции х(І); ВД - свободная составляющая движения объекта - характеристика, описывающая поведение объекта в нестабильном начальном состоянии. Модель (1) является описанием объекта в стабильном начальном состоянии в отличие от модели (2), которая позволяет учитывать нестабильное начальное состояние объекта.

Характеристики и(І) и х(І) могут наблюдаться в ходе экспериментов с помощью разнообразных средств контроля (органы чувств пользователей, контрольно-измерительные приборы и электронные измерительные датчики). Характеристики Іі(І) и ВД - временные характеристики объекта, и они являются характеристиками, идентифицирующими сам объект, то есть информация о поведении этих характеристик позволяет получить достаточное описание поведения объекта. Иными словами, чтобы построить модель объекта в определенных пользователем условиях, достаточно описать временные характеристики объекта и использовать эти описания при поиске модели.

Метод, основанный на наблюдении за поведением объекта в условиях, когда на него воздействуют при помощи специальных воздействий в тестовом режиме, является одним из наиболее используемых методов описания временных характеристик объекта. Нами рассматривается метод временных характеристик, предполагающий воздействие на объект с помощью характеристики и(І) таким образом, чтобы на выходе объекта наблюдались реакции х(І), соответствующие временным характеристикам объекта (или И'(0) и ОД.

Для нахождения описания переходной характеристики объекта при помощи метода временных характеристик необходимо на вход объекта подать тестовый сигнал, описывающийся функцией Хэвисайда 1(1). На выходе объекта будет наблюдаться реакция, приближенно повторяющая поведение переходной характеристики объекта (рис. 1), где и(І) - характеристика, воздействующая на объект; х(І) - реакция объекта на воздействие; із - время запаздывания объекта.

А

/ V « ’ В*1 Тїї?«ят> ■'"'ї'Г'і* І" ^#1'

"Вход" и©

* * * "Выход" х(ї)

50

100

Рис. 1. Идентификация Ь(1)

На входе объекта подается воздействие и наблюдается запаздывающая реакция объекта на это воздействие. Можно заметить наличие случайных помех, оказывающих влияние на измерения реакции объекта. Для нахождения описания поведения свободной составляющей движения объекта необходимо привести его в интересующее нестабильное начальное состояние с последующим прекращением воздействий на объект [5]. Тогда на его выходе будет наблюдаться реакция объекта, приближенно повторяющая свободную составляющую движения объекта (рис. 2).

Рис. 2. Идентификация Щ

На рисунке 2 показано, как при помощи входного воздействия и(1) в интервале времени меньше момента времени 1=0 происходит направление объекта к интересующему нестабильному начальному состоянию. Далее в интервале времени с момента 1=0 до 1=100 воздействие на него прекращается и на выходе объекта наблюдается запаздывающая по времени на величину 1з реакция объекта х(1), соответствующая свободной составляющей его движения, описывающая объект в этом нестабильном начальном состоянии.

Далее необходимо в условиях неопределенности восстановить поведение временной характеристики объекта по наблюдениям характеристик, полученных в результате проведения экспериментов с объектом и построить модель самого объекта. В условиях непараметрической неопределенности и ограниченности ресурсов логично использовать непараметрические оценку Надарая-Ватсона [6], оценку Пристли [7] и оценку регрессии Пристли [2].

Модифицированная оценка Пристли:

К(і )=■

1

1=1

і — і.

\

(3)

где Н[(і) - непараметрическая оценка импульсной переходной характеристики объекта, построенная по выборке значений переходной характеристики объекта іі^), в моменты времени ^ 5 - объем выборки значений іі^); с5Й- параметр размытости непараметрической оценки (3); Н' - колоколообразная функция непараметрической оценки.

Непараметрическая оценка Надарая-Ватсона:

^(і)= ~ГІ-Її )• Н

ґі — О

V С У

(4)

где /5 (ї) - непараметрическая оценка свободной составляющей движения объекта; /ї ) - значение свободной составляющей движения объекта в момент времени ^; с{ - параметр размытости непараметрической оценки.

Для построения модели выходной характеристики (5) используются оценки (3) и (4):

с

X()=Л( -10)+1ИМ - Т)• , где (5)

х11 (/) - непараметрическая модель реакции объекта на входное воздействие и(1).

Для оценки качества полученных моделей используется несколько критериев, например, ошибка моделирования выходной реакции объекта (6), средняя абсолютная относительная ошибка моделирования выходной реакции объекта (7):

& = - І (х{і,)— *5 (і, ))2

5і=1

1 5

&2 = - І

5,=1

Х(І, )— *5 (І, )

*(ї, )

•100%, где

(6)

(7)

х(ї,) - значение выходной характеристики объекта в момент времени ^.

Оценка качества будет происходить путем сравнения при помощи критериев (6) и (7) реакций объекта и модели на одинаковое входное воздействие. Работа метода показана на сложном динамическом объекте (8) - параметрической модели объекта, заданной дифференциальным уравнением двенадцатого порядка.

— 12 — 11 79 7 8

х(ї) + 4.598—— х(ї) + 0.926—- х(ї) + 4.007 —- х(ї) + 3.6988х(ї)

—

+

—6 —5 — 4

+ 0.256—— х(ї)+ 6.032—- х(ї) + 2.522—- х(ї) +

—ї6 —ї5 —ї4 W

+ 0.418-7т х(ї)+4.152-7— х(ї) + 0.854—х(ї)+ 0.211х(ї ) = и(ї — ї3).

Жї3 Жї2 сії 4 7 4 3/

(8)

По выражениям (3) и (4) при определенных параметрах было найдено непараметрическое описание временных характеристик объекта и построены графики непараметрических оценок (рис. 3-4).

Рис. 3. Непараметрическая оценка импульсной переходной характеристики объекта

Рис. 4. Непараметрическая оценка свободной составляющей движения объекта

Построенная непараметрическая модель объекта по формуле (5) позволяет описать реакцию объекта на произвольное входное воздействие с учетом нестабильного начального состояния объекта. Сравнение полученной непараметрической модели объекта с реакцией самого объекта показано на рис. 5.

Показатели критериев качества полученной модели объекта, рассчитанные по формулам (6) и (7), составили соответственно 0і~1000, 02^24 %. Данные значения показывают, что полученная модель доста-

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

Рис. 5. Сравнение реакций объекта и модели

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

Литература

1. Бронов С.А. Способы управления двигателем двойного питания // Информатика и системы управления. - Красноярск: НИИ ИПУ, 2000.

2. Рубан А.И. Методы анализа данных: учеб. пособие. - Изд. 2-е, испр. и доп. - Красноярск: ИПЦ КГТУ, 2004. - 319 с.

3. Красовский А.А. Справочник по теории автоматического управления. - М.: Наука, 1987. - 712 с.

4. Бронов С.А. Электропривод для приборов с индукторным двигателем двойного питания в квазиша-

говом режиме // Изв. вузов. - 2004. - № 4. - С. 61-65.

5. Бродовский В.Н., Иванов Е.С. Исполнительные устройства с многофазными электродвигателями переменного тока // Устройства и элементы систем автоматического регулирования и управления. Техническая кибернетика. - М.: Машиностроение, 1976. - Кн. 3.

6. Рубан А.И. Теория вероятностей и математическая статистика: учеб. пособие. - Изд. 2-е, испр. и

доп. - Красноярск: ИПЦ КГТУ, 2002. - 320 с.

7. Надарая Э.А. Непараметрические оценки кривой регрессии // Тр. ВУ АН ГрССР. - Тбилиси, 1965. -Вып. 5. - С. 55-68.