Технические
науки. Автоматизированные
системы
управления на производстве
Тимонин Д.В.
Самарский Государственный Университет Путей Сообщения
Алгоритм параметрической идентификации в условиях
минимальной априорной информации об особо аддитивных помехах наблюдений
Рассмотрим стационарную нелинейную
динамическую систему, которая описывается следующим разностным уравнением:
, (1)
где выходная переменная наблюдается с
аддитивными помехами в виде:
Требуется
по наблюдаемым конечным выборочным реализациям последовательностей и при известных
порядках и (1) определить
оценки истинных значений параметров.
В [1] показано, что для получения
состоятельных оценок параметров (1) применим следующий критерий:
(2)
где - компакт, ,
- скалярное произведение, ,
, , ,
,,
, ,
Для получения численного метода вычисления
оценок параметров из критерия (2) рассмотрим функцию:
,
,
тогда
(3)
Заметим, что на интервале непрерывна, где
наименьший
корень уравнения:
Далее
,
где
(4)
И что, на
Следовательно, корень существует,
находится на интервале и является
наименьшим, среди всех корней уравнения .
Уравнение (3) можно записать в виде:
Рассмотрим особенности численных методов
определения оценок параметров:
Лемма 1. Для функции справедливы
следующие утверждения:
1. Все корни уравнения (если они
существуют) неотрицательны
2. Уравнение на полусегменте
имеет не более
одного корня, если , где - наименьшее
обобщенное собственное число матрицы (4)
(5)
3. Непрерывность матрицы и существование
корня на интервале является
необходимым и достаточным условием существования единственного решения (2), при
этом:
,
где определяется
уравнением (4)
Эта лемма практически доказана, если
учесть условие утверждения 1 [1].
Утверждение
1.
Пусть последовательность определяется
следующим алгоритмом:
Шаг 0. .
Шаг 1. , определяется из
(5).
Шаг 2. Вычислить и из системы
уравнений:
(6)
Шаг 3. Вычислить
Шаг 4. Проверить условие
.
Тогда, если уравнение имеет корень , то последовательность конечна и , в противном случае она бесконечна.
Доказательство утверждения 1 немедленно
следует из леммы 1.
Этот алгоритм позволяет определить
начальное приближение , необходимое для дальнейшего применения метода Ньютона
или определить, что корень не существует.
Утверждение 2.
Пусть существует , тогда
, , ,
где , и определяются
совместно следующим алгоритмом.
Шаг 1. Вычислить и из системы
уравнений (6).
Шаг 2. Вычислить
Шаг 3. Переход к шагу 1.
Вычисления заканчиваются, если выполняется
условие:
,
где - априорно
заданная точность нахождения оценки.
Это утверждение непосредственно следует из
метода Ньютона:
Обоснованность использования метода
Ньютона вытекает из того, что функция непрерывна на
интервале и, и на интервале .
На основе вышеописанного алгоритма создано
программное обеспечение, позволяющее получать оценки матриц параметров с
наперед заданной точностью.
Литература:
1. Кацюба О.А., Тимонин Д.В. “Математические
методы в технике и технологиях ММТТ-21” – Численный метод идентификации
параметров нелинейных динамических систем при наличии помех наблюдений //
Сборник трудов «Математические методы в технике и технологиях - ММТТ-21». –
Саратов: Изд-во Сарат. гос., 2008. – Т. 2.