- •Моделирование систем управления
- •Два аспекта понятия моделирования. Понятие об идентификации.
- •Причины необходимости создания новых моделей
- •Характеристики объектов и процессов, которые надо учитывать при создании моделей
- •Приемы упрощения моделей
- •Этапы построения моделей
- •Определение цели получения модели
- •Определение ограничений и условий, учитываемых при построении моделей
- •Выбор подхода к решению задачи получения модели
- •Определение класса модели. Выбор метода решения задачи и ее решение
- •Общие принципы построения аналитических моделей
- •Модель поплавкового уровнемера
- •Модель процесса теплопередачи
- •Модель смесителя.
- •Модель реактора
- •Модель емкости с изменяющимся уровнем
- •Метод наименьших квадратов в одномерном случае
- •Метод наименьших квадратов в многомерном случае
- •Рекуррентный метод наименьших квадратов
- •Взвешенный мнк и другие разновидности мнк
- •Получение модели по частотным характеристикам
- •Идентификация систем по переходной характеристике
- •Идентификация звена 1-го порядка по переходной функции
- •Идентификация звена 1-го порядка с запаздыванием по переходной функции
- •Идентификация параметров колебательного звена 2-го порядка
- •Определение параметров апериодического звена 2-го порядка
- •4.1. Метод отрезков Ta и Tb
- •4.2. Метод отрезков Tb и Tc
- •Идентификация по апериодической переходной функции с точкой перегиба звена первого порядка с запаздыванием
- •Метод кратных корней
- •Метод площадей
- •Основное уравнение идентификации
- •Решение основного уравнения идентификации
Определение параметров апериодического звена 2-го порядка
При коэффициенте затухания больше единицы звено второго порядка становится апериодическим и его передаточную функцию можно представить в виде:
Соответствующая переходная функция:
Горизонтальная асимптота дает величину К.
Т1и Т2находятся по точке перегиба и положению касательной в точке перегиба.
Существует несколько методов определения T1и Т2.
Рис. Апериодическое звено 2-го порядка.
4.1. Метод отрезков Ta и Tb
Пусть T1<T2и Т1/Т2= d; 0d1.
Находим по кривой h(t) интервалы Таи Tb (см. рис.).
По аналитическому выражению h(t) для звена второго порядка можно вывести соотношения:
Приложение. |
Точка перегиба |
Уравнение касательной в точке перегиба: |
По этим соотношениям можно построить две кривые:
Рис. Метод отрезков Ta, Tb.
Кривые f1(d) и f2(d) используются в следующем порядке:
имея Taи Tb, вычисляют Ta/Tbи по кривой f1находят d. Затем по кривой f2для найденного значения d определяют Т1/Tb. Из Т1/Тbи Tbвычисляют Т1. Затем вычисляется Т2= Т1/d.
Кроме того, по экспериментальной кривой можно определить интервал Tc= Т1+ Т2и сравнить с найденными значениями Т1и Т2. При большом расхождении следует искать другой метод идентификации.
Примечание. Похожий метод (метод отрезков h0и Tb) описан в учебнике по ТАУ под ред. А. В. Нетушила, 1976, т.1, стр.265.
4.2. Метод отрезков Tb и Tc
Этот метод можно рассматривать как более совершенный метод идентификации апериодических объектов 2-го порядка.
На экспериментальной h(t) (см. рис.33) измеряется Тbи Тc= Т1+ Т2и вводятся обозначения
Можно показать, что
.
По этим уравнениям можно построить кривую y = f(x) (для 1 xe/2 и 0y0.25) и использовать ее для нахождения у по экспериментальному значению x.
Рис. Метод отрезков Tbи Tc.
Решая систему уравнений
относительно неизвестных T1и T2, получаем:
Если x > e/2, то система, возможно, не второго порядка.
Идентификация по апериодической переходной функции с точкой перегиба звена первого порядка с запаздыванием
В некоторых случаях (например, при ориентировочных расчетах) можно h(t) апериодического звена второго или большего порядка аппроксимировать передаточной функцией вида:
Проводя касательную в точке перегиба, как показано на рис. 26, в первом (грубом) приближении можно принять = Ta; T = Tb(см. кривую 1 на рис. 28).
Рис. 28. Аппроксимация с чистым запаздыванием
Более подходящие значения и Т могут быть найдены, если потребовать, чтобы аппроксимирующая ha(t) проходила через точку перегиба и чтобы касательная для h(t) в точке перегиба была бы также касательной и для ha(t). [См. ТАУ, часть 1, под ред. А.В.Нетушила, 1976, стр.263-264]
При t > :Положим
тогда будем иметь систему уравнений:
Решение системы:
Метод кратных корней
Определение постоянной времени Т и кратности n для модели системы в виде
можно выполнить с использованием информации о всей кривой переходного процесса h(t).
Метод использует вычисление площади, заключенной между уровнем установившегося значения К и кривой h(t):
Рис. 29. Площадь S.
Поскольку теоретически для искомой модели переходная функция h(t) описывается выражением:
то
В последнем выражении находим табличный интеграл:
После подстановки табличного интеграла получаем:
S = K T n
Алгоритм вычисления T и n по экспериментальной кривой может быть следующим:
1. По экспериментальной кривой или по статической характеристике оценивается величина К.
2. Каким-либо подходящим методом численного интегрирования вычисляется площадь
3. Для последовательных значений n = 1, 2, 3 ... вычисляют соответствующие значения по формуле
4. Для каждой пары значений n и T можно вычислить затем h(t) и оценить тем или иным способом разность h(t) , например:
5. В качестве результата выбирается та пара n и T, которой соответствует минимум .