Methods for identification and optimization in intelligent systems

Abstract

The paper presents several algorithms for identification, describing complex objects or control system, including: an algorithm for identification of linear time-varying control systems, multi-row algorithm for decision-making by the method of group record of arguments and neuropackage to solve the problem of functions approximation

Full Text

Процесс разработки системы управления начинается с создания конструктором концептуальной модели этой системы, при этом ключевая информация при решении задач анализа, синтеза и оптимизации САУ сосредоточена в математической модели (ММ) объекта управления[1]. Существует достаточно много сложных объектов, для которых протекающие физические процессы мало изучены, а значит, применение фундаментальных законов физики для построения ММ неэффективно. В этом случае проблема построения ММ относится к проблеме идентификации математических моделей объектов управления. Рассмотрим следующие примеры алгоритмов и методов идентификации: алгоритм идентификации линейных переменных во времени систем управления, метод Группового Учета Аргументов (МГУА) на примере Байесовского алгоритма принятия решения и нейросетевая аппроксимация функции. Алгоритм идентификации линейных переменных во времени систем управления (алгоритм разработан К.А. Пупковым и описан в [1]) Если динамику системы управления описать с помощью обыкновенного линейного дифференциального уравнения и если в этом уравнении хотя бы один коэффициент (параметр) является функцией времени, то такую систему будем называть переменной во времени. Физически это означает, что в процессе функционирования некоторые элементы системы изменяют свойства под воздействием внешних факторов или каких-либо внутренних эффектов. Рассмотрим блок схему устройства, приведенную на рисунке 1. Рисунок 1 Рисунок 2 На этом рисунки показано, что на вход системы подается сигнал , представляющий собой сумму синусоидальных гармоник с частотой и произвольной амплитудой . На выходе получаем некоторый сигнал: , в структуре которого сказалась переменность параметров системы, так как амплитуда и фаза являются функцией времени. На входы устройств умножения (в общем случае их может быть ) подается выходной сигнал и соответственно сигналы: и . Такой прием называется синхронным детектирование. На выходе фильтров низких частот будем получать соответственно мгновенные значения вещественной и мнимой составляющих вектора частотной характеристики системы управления, переменной во времени. Покажем это для случая частоты и амплитуды (1) (2) Чтобы определить амплитуду и фазу на всех частотах входного сигнала, необходимо провести одновременно 2n операций умножения. Если входной сигнал состоит из n синусоидальных компонентов , то в каждом из приведенных выше выражений содержится n членов с частотой ωi (плюс частота каждого компонента входного сигнала ωi и n членов с частотами Если параметры системы не меняются во времени , то член в (1) и соответствующий член в (2) будут постоянными. Если происходит медленное по отношению к изменениям входного и выходного сигналов изменение параметров системы, то эти члены будут характеризоваться существенно более низкой частотой, чем остальные (2n-1) членов каждого уравнения. Физически выделение компонентов отражающих переменность во времени системы, достигается применением фильтров низких частот (рисунок 3.1). Обозначим эти члены и соответственно, причем в рассмотренном случае . Интересующая нас информация содержится именно в указанных членах, так как является мгновенным значением действительной части вектора частотной характеристики на частоте - его мнимой части. Непрерывное определение действительной и мнимой частей вектора эквивалентно непрерывному определению амплитуды и фазы на частоте ω1. На рисунке 2 показано, каким образом на основе измерений R и I можно построить комплексную вектор-функцию для переменных во времени систем в текущие моменты времени t1, t2,… . Метод Группового Учета Аргументов (МГУА) теории эвристической самоорганизации [2] Автор этого метода А.Г. Ивахненко пишет: «Осуществляется целенаправленный перебор многих моделей-претендентов различной сложности по ряду критериев. В результате находится модель оптимальной структуры в виде одного уравнения или системы уравнений. Минимум критерия селекции определяет модель оптимальной структуры». Рассмотрим технологию МГУА [3] на примере работы многорядного алгоритма принятия решения для формулы Байеса из теории статистических решений: (1) где: dj - решение в пользу j-го класса; Ri - класс распознавания ( где i=1,2);Р(Ri) - априорная вероятность класса; Х1, Х2,…,Хn - дискретные аргументы (где Хi=1,2,…,q, где q - число уровней дискретизации); Р(Х1, Х2,…,Хn/Ri) - экспериментальная матрица, представляющая собой вероятностные распределения по классам; L[Ri,dj(Х1, Х2, …, Хn)] - матрица потерь. Возможность преобразования к многорядному представлению функции, заданной выражением (1), заключается прежде всего в наличии сомножителя Р(Х1, Х2, …, Хn/Ri). Заменяя полное описание (1) системой «частных» описаний, составленных для всевозможных пар аргументов, в первом ряду селекции получим следующие выражения: (2) Лучшие решения, отобранные из всех решений системы (2) на первом ряду селекции, принимаются в качестве аргументов второго ряда; лучшие решения, отобранные на втором ряду - в качестве аргументов для третьего ряда селекции и т.д. Таким образом реализуется многорядная вероятность модель принятия решения с иерархической структурой - для принятия решения предыдущего ряда. В рассмотренном алгоритме решение о классе Ri, к которому следует относить пару дискретных переменных a1am К-го ряда селекции, принимается по минимуму риска: , определяемому на отдельной проверочной последовательности данных. Как видно из выражения (3), адаптация (изменение) матрицы потерь L[•] позволяет ускорить процесс достижения максимальной точности принятия решений по малому числу рядов селекции. С помощью математической модели осуществляется прогноз состояния системы на некотором интервале времени, который выбирается из практических соображений в соответствии с режимом функционирования объекта. Нейросетевое моделирование В области нейрокибернетики развиваются три основных направления, первое из которых связано с построением и применением искусственных нейронных сетей. Далее для их программной реализации разрабатываются нейропакеты, а для их физической реализации создаются нейрокомпьютеры [4]. Задача аппроксимации функции заключается в определении по экспериментальным данным функции, наилучшим образом приближающейся к неизвестной и удовлетворяющей определенным критериям. Эта задача актуальна при моделировании сложных динамических систем, при решении многочисленных инженерных и научных задач моделирования. Одной из таких задач является аппроксимация (подгонка) выборки входного сигнала полиномом. В нейропакете Neural Network Toolbox for Matlab для решения аппроксимации функции можно использовать обобщенно-регрессионную нейронную сеть (НС), сеть с радиальными базисными элементами и линейную НС. [5]. Результаты эксперимента показали, что при использовании сетей с радиальными базисными элементами качество аппроксимации улучшено, а при использовании линейных НС результаты моделирования почти такие же. Однако при использовании альтернативной полиномиальной сети по архитектуре МГУА (Group Method of Data Handling, GMDH) из пакета NeuroShell2 качество аппроксимации может быть улучшено. Выводы В статье рассмотрены несколько алгоритмов методов идентификации, описывающих сложные объекты или системы (в том числе интеллектуальные). К ним относятся разработанный [1] частотный метод анализа нелинейной динамической системы управления и на примере метода МГУА [3] показано эффективное испольэование селекционного процесса самоорганизации для определения модели оптимальной структуры [4].