Файл: Методические указания по дисциплине Идентификация и диагностика систем к контрольной работе идентификация динамических систем по методу мнк.doc

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 12.12.2023

Просмотров: 72

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ

ТЮМЕНСКИЙ ГОСУДАРСТВЕННЫЙ НЕФТЕГАЗОВЫЙ УНИВЕРСИТЕТ

ИНСТИТУТ ГЕОЛОГИИ И НЕФТЕГАЗОДОБЫЧИ
Кафедра «Автоматизации и управления»

МЕТОДИЧЕСКИЕ УКАЗАНИЯ

по дисциплине

«Идентификация и диагностика систем»

к контрольной работе

«ИДЕНТИФИКАЦИЯ ДИНАМИЧЕСКИХ СИСТЕМ ПО МЕТОДУ МНК»

Тюмень 2021

Содержание
1. Обратная задача моделирования. Метод наименьших квадратов 3

2. Приведение модели системы к линейно-регрессионному виду 8

3. Построение алгоритма идентификации и проверка результатов 10

4. Задания к лабораторной работе 14

1. Обратная задача моделирования. Метод наименьших квадратов

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

Прямая задача моделирования: пусть сигнал управления задан на некотором интервале времени , известна функция модели F и определенны ее параметры с.

Необходимо определить значение выхода на указанном периоде времени при заданных начальных условиях (н.у.) выходной функции и ее производных (если потребуется).




Рис. 1 Структурная схема задачи прямого моделирования

Обратная задача моделирования (идентификация): известно значение сигнала управления на интервале времени и известна реакция системы на это управление – т.е. измеренный сигнал переходного процесса на выходе .

Необходимо определить вид функции F (идентификация в большом) или параметры модели с (идентификация в малом). Для идентификации в малом необходимо построить алгоритм оценивания параметров модели на основе данных измерений входного и выходного сигналов модели.



Рис. 2 Структурная схема идентификации

Одним из распространенных методов идентификации динамических систем является метод наименьших квадратов (МНК). Рассмотрим динамическую систему, заданную дифференциальным уравнением (иначе - в виде передаточной функции):

(1)

( , - коэффициенты, , ) или в виде:



где - оператор производной, Перепишем это уравнение, оставляя в левой части выходной сигнал :



иначе:


(2)

где введены следующие обозначения: ,

- вектор параметров,

- вектор регрессионных переменных

По-сути мы получили модель в виде , при этом параметры описываются вектором , а функция описывает линейную зависимость между параметрами и производными входного и выходного сигналов. Запись модели в виде (2) называется записью в линейно-регрессионном виде.

При решении задачи идентификации, параметры модели (1) (и соответственно (2)) являются неизвестными и подлежащими определению. Пусть сигнал и регрессионные переменные - доступны к измерению (известны) на периоде . Будем искать оценку параметров . Для этого построим оценку выхода модели (2) с использованием оценки параметров и измеряемых регрессоров : . Запишем квадратичное уравнение невязки и будем искать его минимум:



Данная запись означает, что минимальное отклонение (по квадратичному критерию) измеряемого выхода от его оценки для всех значений времени
достигается при точных оценках параметров .

Для поиска минимума критерия находим его производную по вектору и приравниваем его к нулю:





по правилу:

()

отсюда по правилу дифференцирования сложных функций :





по правилу:



и окончательно:



или в матричном виде (т.е. в виде системы линейных алгебраических уравнений):

(3)

где ( означает запись симметричного относительно диагонали элемента).

.

Тем самым решение задачи идентификации по схеме МНК будет получаться на основе решения системы уравнений (3), которое дает искомые оценки вектора параметров в виде:

, (4)

где - обратная матрица к .

Пример № 1. Пусть динамическая система задана уравнением второго порядка:




или в операторном виде:



Оставляя в левой части переменную , запишем:



или в линейно-регрессионном виде (2):

,

Переходя к дискретному времени, получим:

,

Решающее уравнение схемы МНК запишется в виде ,





2. Приведение модели системы к линейно-регрессионному виду

Для использования схемы МНК необходимо уравнение модели привести к линейно-регрессионному виду (2):

Необходимо отметить, что такое приведение всегда неоднозначно (может быть получено множеством способов).

Отметим некоторые особенности построения схемы МНК:

1) Приведение модели к линейно-регрессионному виду всегда начинается с составления уравнения модели, включающего только измеряемые переменные. В данном уравнении производится группировка относительно переменных, после чего выбирается выходная переменная , и формируется вектор регрессоров в правой части.

Пример: в уравнении

После группировки получаем ;

Выбираем в качестве выходной переменную : ;

Заметим, что в этом уравнении можно сгруппировать два последних слагаемых относительно общего параметра