Программная реализация цифрового регулятора

СКОРОСТИ И МОДЕЛИРОВАНИЕ КОНТУРА СКОРОСТИ

 

Цель работы

Формировать умения:

– формулировать цели проекта (программы) решения задач
по программированию цифровых регуляторов;

– использовать информационные технологии при проектировании ЦАСУ;

– применять методы анализа синтезируемых цифровых систем управления;

– использовать компьютерные технологии моделирования ССДМ ЦАСУ с применением векторно-матричного представления рабочих программ цифровых регуляторов и обработки результатов моделирования на этапе проектирования.

 

УРАВНЕНИЯ СОСТОЯНИЯ ЛИНЕЙНЫХ НЕПРЕРЫВНЫХ СТАЦИОНАРНЫХ СИСТЕМ

 

В современной теории управления широко применяется метод пространства состояний, который базируется на математических моделях в переменных состояния [4, 5, 6, 7, 8, 10, 12, 13]. Такая математическая модель представляет собой систему уравнений первого порядка, которую обычно представляют в векторно-матричной форме. Преимущество такого представления заключается в том, что облегчается решение задач анализа и синтеза с применением ЦВМ на этапе проектирования.

Динамические свойства линейной непрерывной стационарной системы описываются векторно-матричными уравнениями:

, (3.1)

, (3.2)

где x ^T = [ x₁ x₂... x_n ] – вектор-состояния размерности (n ´1);

u ^T = [ u₁ u₂... u_m ] – вектор-входа размерности (m ´1);

y ^T = [ y₁ y₂... y_q ] – вектор-выхода размерности (q ´1);

A, B, C и D – матрицы коэффициентов размерности (n ´ n), (n ´ m), (q ´ n), (q ´ m), соответственно.

Уравнение (3.1) связывает переменные состояния и входные переменные и называется уравнением состояния. Уравнение (3.2) связывает выходные переменные с переменными состояния и входными переменными, и называется уравнением выхода.

На рис. 3.1 изображена векторная структурная схема линейной непрерывной стационарной системы, построенная по уравнениям (3.1) и (3.2).

 

 

Рис. 3.1. Векторная структурная схема

 

Из схемы следует, что на выходе интегратора (s ^–1) формируется вектор состояния x. Матрица A называется основной матрицей системы;
B – матрица связи вектора входа системы с переменными состояния;
C – матрица связи переменных состояния с вектором выхода системы;
D – матрица связи вектора входа с вектором выхода.

 

СТРУКТУРНЫЕ СХЕМЫ МОДЕЛИРОВАНИЯ

ЦИФРОВЫХ СИСТЕМ

 

Применительно к цифровым системам управления используют структурные схемы моделирования, которые называются схемами программирования [6, 12, 13]. На основании схем программирования получают разностные уравнения цифровых регуляторов в векторно-матричной форме. Уравнения являются основой рабочей программы, которая может быть реализована в виде вычислительной процедуры микропроцессора.

Основное отличие схем моделирования цифровых систем от непрерывных заключается в наличии элемента задержки. Элемент задержки в преобразовании Лапласа представляется звеном запаздывания , что соответствует значению z ^–1 в теории z – преобразовании. Поэтому в структурных схемах моделирования цифровых систем вместо интегрирующего звена s ^–1 будет стоять элемент задержки z ^–1. При этом переход от непрерывной модели к цифровой осуществляется заменой
в соответствии с формулой трапеций:

,

где Т ₀ – период квантования.

Переход с применением формулы трапеций от аналоговой модели регулятора к цифровой будем называть дискретной аппроксимацией регулятора.

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

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

Далее выполняются следующие операции:

– моделирование синтезированного ЭП с цифровым регулятором и уточнение периода квантования с применением прикладной программы MatLab;

– составление структурной схемы программирования и определение разностных уравнений цифрового регулятора;

Модель цифрового регулятора в переменных состояния описывается векторно-матричными уравнениями следующего вида:

X [(k +1) Т ₀] = AX (kТ ₀) + B D U (kТ ₀);

U _цр(kТ ₀) = СX (kТ ₀) + D D U (kТ ₀),

где U _цр(kТ ₀) – вектор выхода цифрового регулятора;

D U (kТ ₀) – вектор входного сигнала рассогласования.

Непосредственное программирование. Передаточная функция цифрового регулятора имеет вид:

.

После преобразования будем иметь:

. (3.3)

Перепишем (3.3) в виде:

, (3.4)

где

. (3.5)

По выражениям (3.47) и (3.48) находим:

, (3.6)

. (3.7)

Уравнения (3.6) и (3.7) реализуются в виде программы, структурная схема которой изображена на рис. 3.2.

 

 

Рис. 3.2. Структурная схема непосредственного программирования

 

В качестве переменных состояния системы принимаем выходные сигналы элементов задержки х ₁(kT ₀) и х ₂(kT ₀). По схеме (рис. 3.2) составляем уравнения состояния:

(3.8)

(3.9)

и уравнение выхода системы:

. (3.10)

Окончательно

 

(3.11)

. (3.12)

 

Последовательное программирование. Передаточная функция системы записывается в виде произведения:

. (3.13)

Применяя к сомножителям преобразование (3.4) получаем структурную схему в переменных состояния, изображенную на рис. 3.3.

 

 

Рис. 3.3. Структурная схема последовательного программирования

 

Из схемы находим уравнения состояния и выхода системы
в векторно-матричной форме:

Параллельное программирование. Представим передаточную функцию (3.8) в виде суммы передаточных функций:

(3.14)

где К ₁, К ₂ зависят от а ₁, а ₂и К и вычисляются в результате разложения (3.13) на элементарные дроби.

Структурная схема, составленная непосредственно по слагаемым выражения (3.14) показана на рис. 3.4.

 

 

Рис. 3.4. Структурная схема параллельного программирования

 

Обозначив выходные сигналы элементов задержки символами переменных состояния x ₁(kT ₀) и x ₂(kT ₀), получаем:

;

.

Диагональные элементы основной матрицы уравнений состояния являются полюсами заданной передаточной функции, что подтверждает правильность составления уравнений по схеме, показанной на рис. 3.4.

3.3. u -ПРЕОБРАЗОВАНИЕ

u - преобразование предполагает применение хорошо зарекомендовавших себя в инженерной практике методов анализа и синтеза непрерывных систем [2, 8] и позволяет с помощью замены

(3.15)

перейти к характеристическому уравнению относительно новой переменной u. Причем, комплексная u-плоскость корней характеристического уравнения

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

Действительно, если ЦСУ устойчива, то в выражении (3.15) модуль числителя должен быть меньше модуля знаменателя.

Правомочность приведенных рассуждений иллюстрируется расположением вектора u на рис. 3.5.

 

Рис. 3.5. υ-плоскость

Наибольший практический интерес представляет применение логарифмических амплитудно-частотных и фазочастотных характеристик с использованием такого понятия, как псевдочастота. Для определения псевдочастоты рассмотрим -преобразование (3.15) и выделим из него комплексную переменную . Тогда

(3.16)

Полагая, что в (3.16), имеем

.

Домножив числитель и знаменатель последнего выражения на степень , получим

.

После преобразований с применением формулы Эйлера можно записать, что

, (3.17)

где – относительная псевдочастота;

Удобнее использовать абсолютную псевдочастоту

,

которая имеет размерность с^{– 1}.

При малых частотах

и абсолютная псевдочастота

.

При изменении частоты w в пределах от до псевдочастота изменяется от –¥ до +¥.

Таким образом, для построения логарифмических частотных характеристик необходимо в дискретной передаточной функции разомкнутой ЦСУ W (z) сделать подстановку (3.15), получить:

(3.18)

и от (3.18) перейти к псевдочастоте w_υ:

.