Технические науки. Автоматизированные

системы управления на производстве

Тимонин Д.В.

Самарский Государственный Университет Путей Сообщения

Алгоритм параметрической идентификации в условиях минимальной априорной информации об особо аддитивных помехах наблюдений

Рассмотрим стационарную нелинейную динамическую систему, которая описывается следующим разностным уравнением:

,                                 (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.