Составление исходных дифференциальных уравнений САР

Два подхода к получению исходных дифференциальных уравнений систем.
Истинные и ложные модели

Рабочие файлы: [Истинная и ложная модели] [Принцип измерения ЧХ]

Если поставлена задача составления исходных ДУ САР, то возможны две ситуации. Либо детальная декомпозиция системы на модули и отдельные звенья возможна, либо нет.

Если декомпозиция возможна, то, опираясь на постулаты о сохранении материи и энергии (для соответствующего энергетического домена) и на закон Ома (в соответствующей формулировке), приступают к составлению исходных ДУ САР, т.е. к созданию истинной модели системы. Истинной будем называть такую модель или такое математическое описание, о которых известно, что они детально соответствуют физической природе системы.

Если декомпозиция на модули и звенья для системы невозможна, то, не имея детальной информации о ее физической природе, можно получить лишь ложную модель или ложное математическое описание, которые, однако, позволят исследовать систему и получить адекватные результаты. В этом случае совокупность исходных ДУ САР получают через частотный домен, путем экспериментального снятия частотных характеристик.

Для физической системы порядок системы ДУ ее истинной модели обычно в десять и более раз выше порядка системы ДУ ее ложной модели (например, для моделей ОУ). Тем обусловлена широкая популярность ложных моделей, и типовых звеньев, как структурных элементов для их создания.

Общая форма записи систем ДУ

В целях формализации процесса составления исходных ДУ систем используют такие методы, как "Метод контурных токов", "Метод узловых потенциалов" и их аналоги, имеющиеся во всех энергетических доменах. В результате их применения получается единая система:




(1)

a11(p) x1 + a12(p) x2 + ... + a1k(p) xk = f1(t)
a21(p) x1 + a22(p) x2 + ... + a2k(p) xk = f2(t)
       ...
ak1(p) x1 + ak2(p) x2 + ... + akk(p) xk = fk(t)

где:

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

  1. в форме Коши;
  2. в пространстве состояний;
  3. решенную относительно регулируемой величины - y(t);
  4. решенную относительно ошибки - x(t);
  5. в виде передаточных функций - W(p), F(p), Fx(p).

Форма Коши

Форма Коши
Матричная форма записи системы ДУ решенных исключительно относительно первой производной координат САР.

Структура блок-схемы решающая систему ДУ представленную в форме Коши

Изображенная на рисунке блок-схема позволяет решить систему ДУ представленную в форме Коши:

(2)

где:

О форме Коши:

Пространство состояний

Рабочие файлы: [ABCD.vsm] [ABCD-ПФ и W(jw)]

Пространство состояний (ABCD-форма)
Матричная форма записи системы ДУ САР адаптированная для теории управления путем выделения из формы Коши алгебраических уравнений связывающих внутренние координаты САР с выходной(ыми). Применяется для описания САР большого порядка, как правило, с несколькими входами / выходами и с перекрестными связями.

Структура блок-схемы решающая систему ДУ представленную в форме "Пространства состояний"

Изображенная на рисунке блок-схема позволяет решить систему ДУ представленную в форме "Пространства состояний":

(3)

u' = A u + B x
y = C u + D x

где:

О форме "Пространство состояний":

ДУ решенное относительно регулируемой величины y(t) - уравнение движения

Рабочие файлы: [САР]

Система ДУ (1) может быть преобразована к одному уравнению путем исключения промежуточных координат (обычно выходную координату выражают через координату задания):

.

Результатом подобного преобразования является уравнение движения системы:

(4)

D(p) y(t) = R(p) g(t) - N(p) f(t) ,

где:

ДУ решенное относительно ошибки x(t) - уравнение ошибки

Рабочие файлы: [САР]

Если система ДУ (1) решается относительно ошибки системы, то получается уравнение ошибки замкнутой системы:

(5)

D(p) x(t) = Q(p) g(t) + N(p) f(t)

где:

Передаточные функции САР

Передаточная функция
Функция, связывающая один входной и один выходной сигналы САР. Является формой записи системы ДУ САР решённой относительно требуемой выходной координаты. Обычно ПФ записывается не для временного домена, а для домена Лапласа, связывая в этом варианте не сигналы (т.е. не функции времени), а их изображения.

ПФ-ии получают из ДУ решенного относительно требуемой координаты системы (уравнение (4) или (5)). Для чего правую часть уравнения делят на характеристический полином D(p). Отношения полиномов в правой части при возмущающих воздействиях и есть ПФ-ии.

Типовая структура замкнутой САР

Для типовой структурной схемы замкнутой САР различают 3 основные ПФ, применяемые для исследований:

  1. W(p) = y(t)/x(tЧ Wос(p) = Wрег(pWо(pWос(p) - ПФ разомкнутой системы;
  2. F(p) = y(t)/g(t) - ПФ замкнутой системы;
  3. Fx(p) = x(t)/g(t) - ПФ замкнутой системы по ошибке.

Запишем по структурной схеме уравнение движения для разомкнутой системы:

Уравнение движения для разомкнутой системы.

где: Wпрк(p) - ПФ прямого канала системы.

Замкнем систему с помощью уравнения замыкания:

Уравнение замыкания.

Тогда совместное решение даст уравнение движения замкнутой системы:

Уравнение движения замкнутой системы

и уравнение ошибки замкнутой системы:

Уравнение ошибки замкнутой системы.

При отсутствии помехи  f (t) выходная величина связана с задающим воздействием ПФ замкнутой системы:

Передаточная функция замкнутой системы .

А ошибка - с задающим воздействием ПФ замкнутой системы по ошибке:

Передаточная функция замкнутой системы по ошибке .

При этом:

Другие связывающие отношения

Разделим уравнение движения (4) на уравнение ошибки (5), считая, что f(t)=0 и Wос(p)=1:

y(t) / x(t) = R(p) / Q(p) ,     =>     W(p) = R(p) / Q(p) .

В соответствии с теми же уравнениями и уравнением замыкания характеристический полином D(p) = R(p) + Q(p). Добавим 1 к W(p):

1 + W(p) = Q(p) / Q(p) + R(p) / Q(p) = D(p) / Q(p) .

При исследованиях характеристический полином приравнивают к нулю, т.е. вместо него можно использовать W(p):

1 + W(p) = 0 , - характеристическое уравнение.

А так же:

W(p) = [D(p) - Q(p)] / Q(p) = D(p)/Q(p) - 1 = R(p) / [D(p) - R(p)] .

и

W(p) = F(p) / [1 - F(p)] ,                   W(p) = [1 - Fx(p)] / Fx(p) .

Линеаризация ДУ САР

Достаточно часто встречаются звенья, имеющие нелинейную зависимость между входной и выходной координатами. Если для малых отклонений от установившегося режима нелинейность несущественна, то в этом случае до составления исходных ДУ САР (1), ..., (5) выполняют процедуру линеаризации.

Суть линеаризации

Представление сигналов на входах и выходах САР для решения задачи линеаризации

  1. Пусть нелинейное динамическое уравнение звена имеет вид:

    F(x1, x2, x2', y, y', y'', y''') = j( f, f ') .

  2. Тогда уравнение установившегося состояния, всилу равенства нулю всех производных, имеет вид:

    Fo(x1o, x2o, 0, yo, 0, 0, 0) = j( f o, 0) .

  3. Перейдем к уравнению динамики для отклонений, выполнив подстановки:

    x1 = x1o+Dx1(t),    x2 = x2o+Dx2(t),    x2' = Dx2'(t),
    y = yo+Dy(t),    y' = Dy'(t),    y'' = Dy''(t),    y''' = Dy'''(t);

    и разложив функцию F в ряд:

    .

  4. Завершая линеаризацию, вычтем из левой и правой части уравнение установившегося состояния:

    (*)

    Линеаризованное уравнение модуля системы.

Особенности линеаризованного уравнения

  1. Оно является приближенным - отброшены члены высшего порядка малости.
  2. Неизвестными функциями являются не полные величины, а их отклонения D... от установившихся значений.
  3. Уравнение является линейным относительно отклонений D..., при этом масштабирующие коэффициенты (частные производные) могут быть постоянными или переменными во времени.
  4. Внешнее воздействие линеаризации не подлежит.

Геометрическая трактовка линеаризации

Рабочие файлы: [de_lin.vsm]

Запись линеаризованных уравнений в стандартных для ТАУ формах

Представим линеаризованное уравнение (*) в форме уравнения движения и в виде ПФ.

Уравнение движения предполагает: а) выходную величину и ее производные в левой части уравнения, а входную и все остальные члены - в правой; б) так же, принято, чтобы сама выходная величина входила в уравнение с коэффициентом единица. Чтобы привести линеаризованное уравнение (*) к такому виду введем обозначения:

тогда:

T33Dy''' + T22Dy'' + T1Dy' + Dy = k1Dx1 + k2Dx2 + k3Dx'2 + k4 f1 .

Знак D обычно опускают и записывают уравнение в символьном виде:

(**)

(T33 p3 + T22 p2 + T1 p + 1) y = k1x1 + (k2 + k3 p) x2 + k4 f1 ,

где:  
 
 

T3, T2, T1 - постоянные времени;
k4, k3, k2, k1 - коэффициенты усиления;
p=d.../dt - оператор дифференцирования.

Для вывода ПФ решим уравнение движения (**) относительно выходной величины:

y = W1(p) x1 + W2(p) x2 + W(p)  f1 ,

Более строго передаточные функции определяются через изображения Лапласа или Карсона-Хевисайда, как отношение изображений выходной и входной величин:Пример отображения структурной блок-схемы модуля системы

.

Запись ПФ для переменных во времени величин и для их изображений совпадает до оператора. В первом случае ПФ зависит от оператора дифференцирования p=d.../dt. Во втором случае - от оператора Лапласа s=c+jw.