- •Министерство Образования Украины
- •Задачи идентификации
- •Классификация методов идентификации
- •Основные типы моделей в теории идентификации
- •Основные типы сигналов
- •Методы идентификации моделей объектов типовых звеньев по временным и частотным характеристикам
- •1. Математическая обработка динамических характеристик объектов управления
- •Идентификация параметров модели апериодического звена 1-го порядка по временным характеристикам
- •3. Идентификация моделей в виде апериодических звеньев II-го порядка
- •4. Идентификация моделей в виде передаточной функции колебательного звена II-го порядка по временным характеристикам.
- •5. Особенности идентификации моделей в виде типовых динамических звеньев по частотным характеристикам
- •Методика идентификации моделей в виде передаточной функции по кривым разгона (метод площадей, метод Симою) Постановка задачи
- •Последовательность расчета коэффициентов моделей:
- •Методика идентификации моделей объектов III-го порядка по их временным характеристикам
- •Типы моделей
- •Модель первого типа
- •Модель второго типа
- •Модель третьего типа
- •6. Идентификация линейных динамических систем
- •1. Условия идентифицируемости линейных динамических систем
- •2. Определение весовой функции из уравнения свертки
- •3. Оценивание весовой функции по методу наименьших квадратов
- •7. Регрессионный метод идентификации линейных систем (Метод наименьших квадратов)
- •1. Градиентные методы идентификации нелинейных систем
- •2. Рекуррентное оценивание параметров по методу наименьших квадратов
Основные типы сигналов
Входные сигналы:
ступенчатое единичное воздействие
U(t)=1(t)=
δ-функция:
δ(t)=
гармонический синусоидальный сигнал:
U(t)=a·sin(ωt+φ0)
линейные сигналы:
U(t)=kt+b0
случайный сигнал “Белый шум”.
Методы идентификации моделей объектов типовых звеньев по временным и частотным характеристикам
1. Математическая обработка динамических характеристик объектов управления
Рассматривается объект с одним входом и одним выходом со свойствами: стационарности, линейности, сосредоточенности параметров. На вход подается ступенчатое воздействие и на выходе снимается кривая разгона. Необходимо решить обратную задачу: по известной кривой разгона определить коэффициенты уравнения.
Для представления уравнений в безразмерной форме выполняется математическая обработка кривой разгона. Пересчитывается ордината кривой разгона (операция тарирования) по формуле , где- экстремальные значения выходной величины.
При описании динамических свойств статических промышленных объектов ограничиваются одним из следующих дифференциальных уравнений
T1, T2, T3 - коэффициенты левой части дифференциального уравнения;
T - коэффициент при первой производной в правой части дифференциального уравнения;
R0 - коэффициент усиления объекта.
В уравнении 3-го порядка могут быть T3, T2, = 0 , тогда получаем частные случаи уравнений 1-го и 2-го порядков, и без запаздывания.
Для описания динамических свойств астатических объектов используются дифференциальные уравнения не содержащие члена y(t) и статического коэффициента усиления k0, т.е. имеющих вид:
.
Величина запаздывания может быть определена графически следующим образом (см. рис).
Идентификация параметров модели апериодического звена 1-го порядка по временным характеристикам
Для апериодического звена 1-го порядка
коэффициент усиления k0 определяется следующим образом
Для апериодического звена первого порядка:
- передаточная функция
- переходная функция;
- импульсная переходная функция.
Реакция на единичное входное воздействие
Для определения коэффициента звена:
Временные характеристики.
Если зависимость экспоненциальная, то постоянная времени определяется одним из двух способов. При первом, чисто графическом способе проводится касательная в любой точке графика и берётся разность абсцисс точек: а) касания с графиком и б) пересечения с линией установившегося уровня. При втором способе необязательно производить графические построения. Рассмотрим момент времени t=T: , тогдаh(T)≈0,63k.
С соответствующими корректировками те же способы применимы и при наличии импульсной переходной характеристики. Для второго метода при t=T: , тогда ω(T)≈0,37k.
Если переходная характеристика имеет вид:
,то нужно в передаточную функцию вводить чистое запаздывание:
3. Идентификация моделей в виде апериодических звеньев II-го порядка
;
a0p2+a1p+1=0;
;
;
;
Тогда ,;
Считаем, что процесс апериодический I-го порядка. Для этого начальный участок аппроксимируем.
Получаем кривую, приближенную к зависимости I-го порядка. Считаем точку А началом координат. Взяв 0,63 относительно т. А определяем Т1.
Для определения Т2 строим зеркальную кривую и определяем Т2 на начальном участке кривой по уже известной методике.