2.6 Аппроксимация моделью второго порядка

2.6.1 Аппроксимация двумя инерционными звеньями первого порядка с разными постоянными времени

Если на графике переходной функции объекта h_o(t) не просматривается характерный для s-образных кривых прямолинейный участок, а сама переходная кривая приближается к установившемуся значению h_o(∞)=k_o сравнительно медленно, то для аппроксимации можно выбрать передаточную функцию второго порядка с запаздыванием и различными постоянными времени:

.

Постоянные времени Т_о1 и Т_о2 определяются по некоторым характерным точкам переходного процесса.

В большинстве случаев модель обеспечивает достаточную точность для практических расчетов, если принять Т_о1/Т_о2=0,5. При этом постоянные времени Т_о1 и Т_о2 определяют следующим образом: по ординате h(t₂)=0,63k_o экспериментальной переходной характеристики находят момент времени t₂, отсчитываемый от окончания времени чистого запаздывания, а затем вычисляют Т_о2=0,64t₂ и Т_о1=0,5Т_о2. Такая аппроксимация целесообразна, когда h(0,5 t₂)>0,3k_o.

Полученная передаточная функция в виде:

задается в программе ТАУ-3 Linsys. Полученный переходный процесс оцифровывается по тому же принципу, что и исходная кривая разгона. И производится построение исходной кривой разгона и аппроксимированной с различными постоянными времени на одном графике для произведения сравнительного анализа.

2.6.2 Аппроксимация двумя инерционными звеньями первого порядка с одинаковыми постоянными времени

Модель второго порядка с одинаковыми постоянными времени и запаздыванием имеет вид:

,

где Тoi = То/2,72;

.

Параметры данной модели однозначно выражаются через параметры Т₀ и τ₀ экспериментальной переходной характеристики (рисунок 4).

Рисунок 4. Переходные характеристики реального объекта (кривая 1)

и его приближенных моделей второго (кривая 2) и первого (кривая 3)

порядков с запаздыванием

Этой передаточной функции соответствует кривая 2, которая начинается в точке D΄. Очевидно, что при данном способе аппроксимации почти весь интервал τ_о должен быть смоделирован как чистое запаздывание, т.е. .

Полученная передаточная функция в виде:

задается в программе ТАУ-3 Linsys. Полученный переходный процесс оцифровывается по тому же принципу, что и исходная кривая разгона. И производится построение исходной кривой разгона и аппроксимированной с одинаковыми постоянными времени на одном графике для произведения сравнительного анализа.

2.7 Определение оптимальной передаточной функции

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

Таблица 3

Определение интегральной квадратичной ошибки

для различных моделей аппроксимации

Время	Переходные функции объектов и моделей				Вычисление ошибки
	Исходная	1 пор.	2 п. разл. Т	2 п. один. Т	1 пор.	2 п. разл. Т	2 п. один. Т
Ошибка:					σ1	σ2p	σ2o

Интегральная квадратичная ошибка вычисляется по формуле:

Подсчитывается квадратичная ошибка для каждой модели и объекта: σ1, σ2o, σ2p. Рассчитанные квадратичные ошибки сравниваются между собой, и в дальнейших расчетах за объект принимается модель с наименьшей ошибкой.