Отбор методов по целевой направленности
Как уже указывалось выше, число работ и разнообразие методов идентификации делает практически невозможной достаточно полную их характеристику. Одним из рациональных подходов в этих условиях является отбор методов параметрической идентификации по их целевой направленности, т. е. зависимости от свойств объектов, отражением которых являются модели определенных классов. Растригиным Л.А. предложена, например, следующая классификация моделей под этим углом зрения:
- статические или динамические;
- детерминированные или стохастические;
- линейные или нелинейные;
- непрерывные или дискретные.
Признаки классификации методов
При определении вида оператора связи между входом и выходом объекта в зависимости от его свойств выбирается либо один из приведенных выше типов моделей, либо некоторая их комбинация, что, в свою очередь, используется при выборе наиболее приемлемых методов параметрической идентификации, в основу классификации которых положены следующие признаки:
- активность (пассивные и активные методы);
- адаптивность (неадаптивные и адаптивные);
- дискретность (непрерывные и дискретные, т. е. шаговые).
Как видно из приведенных классификации, число возможных сочетаний моделей и методов довольно велико, но и оно не исчерпывает всего многообразия реальных ситуаций, хотя, безусловно, вносит определенную целенаправленность в процесс выбора методов. Например, идентификация объекта, описываемого статической, детерминированной, линейной моделью осуществляется более простыми методами, чем для случая динамической стохастической нелинейной модели. Из-за невозможности рассмотреть все промежуточные случаи ниже мы остановимся лишь на некоторых из них, наиболее характерных применительно к металлургическим объектам.
Параметрическая идентификация для случая статической детерминированной модели
Математическое описание модели
Допустим, что поведение объекта описывается регулярной зависимостью, связывающей вход и выход объекта:.
Тогда модель объекта также должна представлять собой некоторую регулярную функцию.
Рассмотрим сначала случай линейной модели объекта с входами и выходами, которая имеет единственно возможную структуру и описывается системой линейных алгебраических уравнений:
В векторной форме эта система имеет вид:
,
где– скалярное произведение векторов и:
Рассмотрим сначала неадаптивный шаговый метод применительно к решению этой задачи, для чего приравниваются выходы модели и объекта в каждом изопытов (отсчетов):
Критерий идентификации
В качестве критерия идентификации чаще всего используется суммарная невязка (остаток) модели и объекта
Некоторым преимуществом рассмотренного метода по сравнению с предыдущим является возможность использования текущей информации, но при этом возникают проблемы сходимости процесса адаптации, связанные с выбором параметра .
Скорость сходимости может быть значительно повышена при использовании активных методов, например, путем выбора и реализации на объекте такой последовательности состояний входов, при которой значения идентифицируемых параметров могут быть найдены за минимальное число шагов. Наилучшим образом это достигается при условии взаимной ортогональности векторов входных параметров.
Для непрерывного случая процесс адаптивной идентификации представляется дифференциальным уравнением вида:
- ,
- ,
- .
Обозначив вектор-функцию времени
Рис. 5.14 Структурная схема непрерывной адаптивной идентификации
Нелинейные модели
Функциональное представление обьекта
Определение вида функции , связывающей вход и выход объекта , является, как уже отмечалось выше, предметом структурной идентификации.
Существующие представления об объекте отражаются в виде определенной функции с неизвестными параметрами . Процедуры же конкретной реализации функции основываются на кусочно – линейном или функциональном представлении зависимости выхода модели от входа. Остановимся на функциональном представлении, как имеющем большую степень общности [35].
Неизвестная функция объектав виде известной функции с неизвестными параметрами .
Для определения неизвестных параметров приравнивают состояния модели и объекта для каждого из наблюдений
Аналитическое решение
Решение такой системы (в общем случае трансцендентных уравнений) сводится, как и в линейном случае, к задаче минимизации суммарной невязки
При правильном выборе алгоритмаэтот процесс должен сходиться к точным значениям параметров, т. е. к решению задачи (5.76)
Особенности адаптивных методов идентификации применительно к нелинейным моделям
Рассмотрим особенности адаптивных методов идентификации применительно к нелинейным моделям.
Локальная невязка выходов модели и объекта имеет для непрерывного случая вид:
Структурная схема, реализующая этот алгоритм, представлена на рис. 5.15, сравнение которого с рис. 5.14 показывает, что последний алгоритм отличается от случая линейной модели наличием функционального преобразователя , предназначенного для определения вектора
Рис. 5.15 Схема адаптивной идентификации нелинейных объектов
Адаптивно шаговый метод
Рассмотренный случай касался адаптивного алгоритма идентификации для непрерывных объектов. Остановимся теперь на особенностях адаптивного шагового метода, который целесообразно применять для дискретных объектов (при дискретном способе получения информации о состоянии объекта). Локальная невязка в этом случае имеет вид:
Параметр выбирается из соображений оптимизации работы алгоритма. Определим его для предельного случая, когда структура модели и объекта совпадают, т. е.
т. е. при локальная невязка на каждом шаге идентификации уменьшается до нуля.
При активной идентификации векторы на каждом шаге выбираются таким образом, чтобы они были ортогональны друг другу, т. е.
Параметрическая идентификация стохастических объектов
Случай статического стохастического объекта
Сначала рассмотрим случай статического стохастического объекта, который может быть представлен в виде:
Рис. 5.16 Схема объекта с аддитивным наложением помех
Для простоты остановимся на таких объектах (рис. 5.16), у которых регулярная и случайная составляющие могут быть разделены, т. е. представлены в виде:
- – регулярная составляющая объекта;
- – случайная составляющая.
При этом предполагается, что свойства случайной составляющей не зависят от входа , т. е. полностью оцениваются определенной плотностью вероятности , в качестве которой часто принимают нормальный закон.
Тогда для объекта с одним входом , когда , плотность нормального распределения характеризуется двумя параметрами: математическим ожиданием и дисперсией
Процедура декорреляции
Идентификация объектов с несколькими выходами значительно затрудняется при наличии корреляции помех, действующих на разные выходы. Преодолеть этот недостаток можно с помощью процедуры «декорреляции», смысл которой состоит в следующем.
Пусть и – коррелированные случайные величины с нулевыми средними, дисперсиями , и корреляционным моментом , который предполагается известным. Образуем путем линейного преобразования новые случайные величины:
Определим коэффициент корреляции этих величин
При этом необходимо подчеркнуть, что оценки идентифицируемых параметров в среднем совпадает с их точными значениями при выполнении двух следующих условий:
В реальных же условиях такая корреляция часто имеет место. В этом случае необходимо получить значение корреляций , для чего нужно привлечь всю дополнительную информацию о свойствах помехи и ее взаимодействии со входами. При известных значениях (по крайней мере, для случая линейной модели) могут быть получены поправки , к значениям параметров .
Оптимизационная задача
Невязки
Для определения параметров после этапа структурной идентификации, на котором выбирается вид функции , задача сводится, как и в предыдущем подразделе, к оптимизационной. Образуют невязки выхода модели и объекта на каждом – том измерении
Оцениваемыепараметры
Оцениваемые параметры выбираются таким образом, чтобы все эти невязки были минимальны по модулю, т. е. решается задача минимизации функций:
Методы свертывания критерия невязки
Эта задача является многокритериальной и может быть решена только при свертывании критерия, что можно сделать, например, одним из следующих способов:
- модульный критерий;
- квадратичный критерий;
- показательный критерий;
- минимаксный критерий, .(осуществляется минимизация этого критерия);
- взвешенный критерий;
где .
Последний критерий обобщает все предыдущие.
Задача идентификации теперь представляется в виде
Выбор критерия, прежде всего, зависит от характера помех , т. е. от их вероятностных свойств. Например, при нормальном распределении помех наибольшую точность дает минимизация по квадратичному критерию, который удовлетворяет также требованию простоты решения задачи ввиду гладкости функций вида .
Выбор параметра
Рассмотрим теперь особенности выбора оптимального значения параметра в алгоритме адаптивной шаговой идентификации для объекта типа (230), по одному их выходов, к которому приложена помеха , имеющая следующие параметры закона распределения , .
Рассмотрим, как и в предыдущих случаях, локальную невязку на -том шаге
Опуская эти преобразования, имеем [35]
Следует заметить, что прямо воспользоваться этим выражением для адаптивной идентификации нельзя, так как неизвестен вектор невязки параметров , поскольку неизвестны параметры модели, которые мы и пытаемся оценить с помощью искомого алгоритма. Получается как бы замкнутый круг. Однако это препятствие можно обойти путем замены на , которые отличаются на и, следовательно, в среднем совпадают. Тогда приближенно имеем:
Вопрос сходимости для непрерывного объекта
Выше рассмотрен случай пассивной адаптивной идентификации. Если имеется возможность активного воздействия на объект, то выбор вектора ортогонального предыдущим векторам способствует ускорению сходимости процесса, но из-за наличия помехи все же не гарантирует окончания идентификации за конечное число шагов.
Для случая непрерывного объекта и метода вопрос сходимости алгоритма решается проще, чем для шагового метода. Сходимость здесь обеспечивается при условии и достаточной вариабельности вектора . Скорость сходимости прямо пропорциональна , но при очень большом его значении может возникать неустойчивость. Структурная схема алгоритма не отличается от приведенного на рис. 4.14 и 5.15 соответственно для линейной и нелинейной моделей.
Динамические модели
Параметрические и непараметрические модели
Рассмотрим методы идентификации объектов, оператор которых имеет память, т. е. выход в момент отражает не столько состояние входа в этот момент, сколько его значения в предыдущие моменты времени.
Следует различать параметрические и непараметрические модели объектов. В первом случае модель определяется набором параметров (коэффициентов), которые оцениваются в процессе идентификации. Непараметрическая же модель определяется, в общем случае, непрерывной функцией (чаще всего функцией времени). Однако она может быть задана также точками или в виде разложения в ряд по некоторой системе функций. В последнем случае мы опять приходим к параметрической модели.
Линейная параметрическая модель для одномерного случая
Линейная параметрическая модель для одномерного случая представляет обыкновенное дифференциальное уравнение вида
Линейная параметрическая модель с числом входов больше единицы
При числе входов больше единицы для случая линейной модели вектор состояния образуется как сумма векторов:
Тогда уравнение для линейной динамической системы с несколькими входами имеет вид:
Как и в предыдущих случаях, задачу идентификации можно свести к минимизации функции невязки в виде квадрата разности правой и левой частей уравнения (5.89) при подстановке в него функций и наблюдений объекта
Сложности связанные с дискретностью объекта
Решение этой задачи усложняется при дискретном характере объекта или способа получения информации. Здесь возникают дополнительные погрешности, связанные с численным интегрированием и дифференцированием, для чего используются специальные подпрограммы, входящие в состав математического обеспечения современных ЭВМ. Погрешность численного дифференцирования можно снизить применением специальных методов сглаживания исходных сигналов.
Адаптивная идентификация динамической нелинейной модели
Математическое описание модели
Остановимся теперь на случае адаптивной идентификации динамической нелинейной модели. Пусть– известная функция с неизвестными параметрами.
Система уравненийрешается либо в непрерывном виде на АВМ либо численно интегрируется на ЦВМ (например, методом Рунге – Кутта) при заданных начальных условиях и фиксированных значениях идентифицируемых параметров. Полученное решение
При адаптивном методе идентификации невязка в момент представляется в виде и минимизируется в каждый момент времени .
Методы решения, структура системы
Эта задача может быть решена с использованием поисковых методов оптимизации (рис. 5.17).
Рис. 5.17 Структурная схема поисковой идентификации
При этом мы не касаемся конкретной реализации алгоритмов блока минимизации, так как они зависят от постановки задачи, свойств реальных объектов и могут требовать в сложных случаях больших затрат вычислительных средств. В отношении же общей структуры системы (рис. 5.17) следует отметить, что наличие в ней операторов дифференцированияи, способных усиливать помехи (особенно высокочастотные), приводит к необходимости применения соответствующих методов и устройств и для сглаживания исходных сигналов и .
Проблема сглаживания и фильтрации
Проблема в общем случае
Проблема сглаживания и фильтрации в общем случае является достаточно сложной, ей посвящено много работ, обзор которых выходит за рамки данного пособия. Основной задачей в этой проблеме является выбор такого оператора сглаживающего фильтра, который бы наилучшим образом подавлял помеху (ее статистические свойства, в том числе частотные, должны быть известны), в наименьшей степени искажая полезный сигнал.
Случай значительного различия между помехами и полезным сигналом
В самом простом случае, когда частотные характеристики полезного сигнала и помехи существенно различны, эту задачу можно решить, пропуская исходный сигнал, например , через инерционное звено, описываемое дифференциальным уравнением первого порядка
Рис. 5.18 К вопросу фильтрации и сглаживания
Свойства фильтра, а значит, и вид сглаженного сигнала можно менять выбором соответствующей постоянной времени .
Дискретным аналогом такого фильтра является оператор экспоненциального сглаживания, который записывается в следующем виде:
- – новое сглаженное значение (оценка среднего значения );
- – предыдущее сглаженное значение;
- – параметр сглаживания.
Выбор параметра сглаживания
При простоте этого алгоритма выбор параметра сглаживания не является столь же простой задачей. Показано, например, что экспоненциальное сглаживание дает результаты, близкие к методу скользящего среднего, если параметр удовлетворяет следующему соотношению
Идентификация непараметрической модели
Описание поведения объекта
Остановимся теперь на методе идентификации для случая непараметрической модели, которая может представляться в виде импульсной (весовой) переходной функции, амплитудной и фазовой характеристик. Известно, что свойства линейного динамического объекта однозначно определяются его реакцией на единичное импульсное возмущение (см. рис. 5.6,а). В этом случае поведение объекта описывается интегралом свертки
Алгоритм для стохастического случая
Аналогом этого уравнения для стохастического случая, когда ко входу объекта наряду с полезным сигналом приложена помеха, является уравнение статистической динамики (5.57) , в котором роль входного и выходного сигналов играют соответственно автокорреляционная функция входаи взаимная корреляционная функция выхода и входа .
Более подробно эти вопросы, в том числе для случая перехода от непараметрической к параметрической форме модели рассмотрены в разделе 3; здесь мы обратились к ним из соображения цельности представления о методах идентификации.
Представление нелинейных динамических объектов линейными моделями
Примечательные модели подобного рода
Ниже рассмотрим случаи, когда нелинейные динамические объекты могут быть представлены моделями линейными относительно идентифицируемых параметров. В этом случае удается построить относительно простые и в то же время достаточно эффективные алгоритмы идентификации. Среди моделей такого рода следует выделить модели Вольтерра и Гаммерштейна.
Модель Вольтерра
Модель Вольтерра связана с рядом, названным его именем, которым представляется выход нелинейной динамической системы
Для удобства идентификации целесообразно представить эту задачу в параметрической форме в виде разложения по определенной системе функций. Линейную часть представим следующим образом:
- ;
- ;
- .
Задача идентификации сводится, таким образом, к определению параметров разложения (5.97), общее число которых.
Интегральная невязка в этом случае имеет следующий вид:
Если идентификация осуществляется в темпе с процессом или используется последовательно поступающая информация, то более удобен алгоритм адаптивной идентификации. В этом случае образуется локальная невязка
Структурная схема реализации этого алгоритма аналогична изображенной на рис. 5.15.
Модель Гаммерштейна
Модель Гаммерштейна основывается на предположении, что нелинейность и динамику объекта можно разделить и представить в виде последовательного соединения нелинейного безинерционногои линейного динамического звеньев:
Описание объекта динамической стохастической моделью
Математическое описание модели
В заключение этого подраздела рассмотрим один из наиболее сложных случаев, когда объект описывается динамической стохастической моделью. Объекты такого рода могут быть представлены в виде:
- – идентифицируемый оператор;
- и – фактические (истинные) значения параметров входа и выхода объекта.
Можно предположить, что измеряемые значения входа и выхода отличаются от фактических на величину ошибок измерения, т. е.
Структурная схема
Структурная схема такого объекта приведена на рис. 5.19, где и – измерители входа и выхода.
Рис. 5.19 Представление стохастического объекта
Связь стохастичности с точностью измерений
Стохастичность объекта связывается, таким образом, с точностью измерений его входа и выхода, в то время как идентифицируемый оператор предполагается детерминированным. В действительности же случайная помеха может входить и в оператор объекта (внутренняя помеха), т. е.
Фильтрация помех
Распространенным подходом к решению задачи синтеза такого оператора является фильтрация помех, что позволяет в ряде случаев приблизить измеренные значения и к истинным , . Для этой цели можно воспользоваться одним из методов фильтрации, которые выбираются, главным образом, в зависимости от характера помехи и ее отношения к полезному сигналу.
Если характер помехи позволяет отфильтровать ее до определенного допустимого уровня, то дальнейшая идентификация может производиться таким же образом, как и в случае детерминированных динамических объектов. Но такой подход не всегда возможен, так как для эффективной фильтрации помех необходимо знать их свойства. При этом, если спектральные характеристики полезного сигнала и помехи близки, то вместе с помехой можно отфильтровать и полезный сигнал, несущий информацию о свойствах объекта. В подобных случаях необходимо использовать другие методы, которые основываются либо на усреднении помех, либо на выборе определенного вида пробных воздействий. Один из таких методов, относящийся к первому подходу, был рассмотрен в разделе 5.3 (см. уравнение статистической динамики и его приложения), там же мы коснулись вопроса нанесения пробных воздействий для детерминированного случая (или близкого к таковому). При наличии же существенного уровня помех эта проблема значительно усложняется, так как для выбора характера пробных сигналов необходимо учитывать свойства объектов и помех.
Добрый день. Подскажите пожалуйста, из какого источника представлена информации на странице: https://bookaa.ru/matematicheskoe-modelirovanie/metody-parametricheskoy-identifikaci.html
Заранее спасибо, с уважением Андрей Александров