![](/user_photo/2706_HbeT2.jpg)
- •1 Лекция 1. Понятие моделирования. Цели моделирования
- •2 Лекция 2. Основные термины в математическом моделировании. Виды математических моделей
- •2.1 Основные термины в математическом моделировании
- •2.2. Основные виды математических моделей
- •3 Лекция 3. Этапы процесса моделирования. Общие принципы построения моделей
- •3.1 Этапы процесса моделирования
- •3.2 Общие принципы построения моделей
- •4 Лекция 4. Аналитические методы определения характеристик объектов
- •4.1 Основные уравнения динамики
- •4.2 Упрощение уравнений динамики
- •4.3 Линеаризация уравнений
- •5 Лекция 5. Аналитические методы моделирования объектов с сосредоточенными параметрами
- •6 Лекция 6. Примеры моделирования объектов с сосредоточенными параметрами
- •7 Лекция 7. Объекты с сосредоточенными параметрами. Динамика теплообменных процессов
- •8 Лекция 8. Аналитические методы определения характеристик тепловых объектов
- •9 Лекция 9. Моделирование объектов с распределенными параметрами
- •10 Лекция 10. Общие сведения о проблеме идентификации
- •10.1 Основные понятия
- •10.2 Классификация методов идентификации
- •11 Лекция 11. Постановка задачи идентификации
- •11.1 Объект идентификации
- •11.2 Постановка задачи идентификации
- •12 Лекция 12. Идентификация линейных динамических объектов. Прямые методы
- •12.1 Прямые методы определения динамических характеристик
- •12.2 Идентификация с помощью переходной функции
- •13 Лекция 13. Идентификация линейных динамических объектов. Прямые методы
- •13.1 Графическая идентификация с помощью переходной функции процессов второго порядка
- •13.2 Графическая идентификация с помощью импульсной переходной функции
- •13.3 Идентификация с помощью частотной характеристики
- •14 Лекция 14. Параметрическая идентификация линейных объектов
- •14.1 Статические детерминированные линейные модели
- •14.2 Динамические детерминированные модели
- •15 Лекция 15. Непараметрическая идентификация линейных динамических объектов. Корреляционные функции
- •15.1 Общий подход к определению непараметрической модели
- •15.2 Определение корреляционных функций сигналов
- •16 Лекция 16. Непараметрическая идентификация линейных динамических объектов. Уравнение Винера-Хопфа
- •16.1 Уравнение Винера-Хопфа для определения импульсной переходной функции
- •16.2 Алгебраический метод решения уравнения Винера-Хопфа
- •17 Лекция 17. Методы идентификации, основанные на
- •17.1 Краткие сведения об аппроксимации функций
- •17.2 Сглаживание дискретных значений импульсной переходной функции
- •17.3 Метод идентификации, основанный на предварительной аппроксимации импульсной переходной функции
- •18 Лекция 18. Методы идентификации, основанные на аппроксимации характеристик объектов и сигналов
- •18.1 Метод идентификации, основанный на совместной аппроксимации импульсной переходной и корреляционных функций
- •18.2 Метод идентификации, основанный на аппроксимации сигналов
- •19 Лекция 19. Идентификация нелинейных объектов
- •19.1 Особенности идентификации нелинейных динамических объектов
- •19.2 Методы идентификации, основанные на линеаризации характеристик объектов
- •19.3 Идентификация нелинейных функций априорно известного вида
- •19.4 Идентификация объекта с нелинейностями общего вида
- •20 Лекция 20. Алгоритмы предварительной обработки и оценки идентичности
- •20.1 Алгоритмы оценки стационарности и линейности объекта
- •20.2 Количественная оценка степени идентичности модели реальному объекту
- •Список литературы
14.2 Динамические детерминированные модели
В одномерном случае связь между входной величиной x =x(t) и выходом y = y(t) представляет собой обыкновенное дифференциальное уравнение:
, (14.5)
вместе
с начальными условиями для
Модель определяется (p+l+1) параметрами c = (a0,…,ap-1,b0,…,bl).
Для описания динамики объектов, которые характеризуются дискретными значениями входных и выходных сигналов, вместо дифференциальных уравнений можно воспользоваться разностными уравнениями. Обозначив дискретные значения входного и выходного сигналов соответственно xk-j =x[(k-j)] и yk-j = y[(k-j)], разностное уравнение (аналог дифференциального) запишем в виде
(14.6)
Необходимо также задать начальные условия. Структурными параметрами модели являются p и l, которые должны быть выбраны в процессе структурной идентификации.
Модель (14.5) часто удобно записывать в виде системы дифференциальных уравнений. Обозначимy1 = y, y2 = y(1), y3 = y(2),… , yk = y(k-l). Тогда (14.5) в виде системы имеет вид:
(14.7)
(14.5) сводится к (14.7), но (14.7) в общем случае не сводится к (14.5).
Можно записать в обобщенной форме
Тогда
векторная форма системы: ,
где Y = (y1,…, yp) – вектор состояний, X’=(x, x(1), …, x(l)) – вектор возмущений, A – квадратная матрица коэффициентов, B – прямоугольная матрица коэффициентов.
Итак, исходной информацией для построения процедуры идентификации является вид идентифицируемой модели (14.5) и наблюдения {Xt, Yt} в промежутке [0,T]. Надо определить ai, bj. В общем случае при подстановке наблюдений в уравнение модели (14.5) равенство в этом уравнении не выполняется. Надо подобрать ai, bj таким образом, чтобы в уравнении (14.5) правая и левая части отличались друг от друга наименьшим образом.
Строим функцию невязки в виде среднего квадрата разности правой и левой части уравнения (14.5) при подстановке туда функций Xt, Yt – наблюдений объекта
.
(14.8)
Минимизируем
(14.8)
по ai и bj: Результат
минимизации – с* дает значения
идентифицируемых
параметров. c*=(a0*,…,ap-1*, b0*,…,bl*)
– решение.
Приравнивая нулю производные функции (14.8) по всем неизвестным параметрам (так как функция гладкая), получим систему линейных уравнений, решение которойявляется решением задачи минимизации.
Полученную систему можно решить стандартными вычислительными методами. Однако для вычисления коэффициентов системы необходимо знать производные входного xt и выходного yt сигналов объекта. Если эти сигналы заданы в аналитической форме – нет никаких проблем. В противном случае надо выбрать метод получения производных; это зависит от конкретных условий и должен определяться каждый раз специально:
1) можно использовать численное дифференцирование, то есть использовать точечную оценку производных:
(14.9)
-
интервал, характеризующий базу оценок.
Этот метод обладает двумя недостатками.
Во-первых, точность оценки производных
быстро падает с ростом порядка производных.
Практически хорошо оценить производные
выше второго порядка не удается.
Во-вторых, для оценки i-той
производной при t=0
необходимо иметь значения функции
приt<0,
то есть z(-∆t),z(-2∆t),
…, z(-i*∆t).
Так как этих значений нет в измерениях,
то приходится в (14.8) брать пределы
интегрирования не [0,T],
а [i*∆t, T],
где i –
номер максимальной производной функции;
2) можно использовать аппарат сглаживания функции с весом ρ(t); либо
3) разложение в ряд по заданной системе функций (из-за наличия операции сглаживания, ошибки в определении производных будут значительно меньше).
Для того чтобы задача идентификации имела решение, необходимо, чтобы определитель системы, полученный в результате решения задачи минимизации, был отличен от 0. Если это условие не выполняется, то следует взять другую реализацию входов и выходов объекта, либо уменьшить количество идентифицируемых параметров (то есть p иl), то есть уменьшить порядок модели. Первая мера не всегда приводит к цели, а вторая – при последовательном применении – всегда.