Современные проблемы теории управления
..pdfозначает отсутствие соответствующей ветви подграфа. При n нечетном выход модели положительный. Необходимое число усилителей Z = n + 3.
Причинно-следственный подграф сигналов, преобразованный к виду, удобному для моделирования, и принципиальная схема модели, соответствующая данному подграфу сигналов, показанному на рис. 6.3, справедливы при любых знаках в полиноме числителя и знаменателе передаточной функции, а также для любых значений ее коэффициентов, т.е. построены на общих основаниях и являются общим случаем.
Рассмотрим возможные пути минимизации числа усилителей, реализующих подграф сигналов, который построен на общих основаниях, например, n = 5.
Рис. 6.3. Принципиальная схема электронной модели подграфа сигналов
В случаях, когда порядок полинома числителя равен порядку полинома знаменателя передаточной функции, т.е. an ≠ 0, подграф сигналов можно преобразовать. Для этого необходимо путь Pn = Cn перенести из узла n в узел n + 1. При этом путь Pn перестанет соприкасаться с контурами Ln–k, где k – нечетные числа от 1 до n. Для получения подграфа сигналов, эквивалентного исходному, необходимо коэффициенты путей Pn–k = Cn–k уменьшить на величину Cn dn–k.
121
Если Cn–k – Cn dn–k ≥ 0, где n = 5, то очевидно:
C0 – C5 d0 ≥ 0; C2 – C5 d2 ≥ 0; C4 – C5 d4 ≥ 0,
и все приведенные неравенства справедливы, то подграф сигналов эквивалентен предыдущему (без доказательства). Если хотя бы одно неравенство не выполняется, то минимизация не допустима. Число усилителей на этом этапе минимизации Z = n + 2, а подграф сигналов приведен на рис. 6.4, а.
Рис. 6.4. Минимизированный граф сигналов
Если окажется, что Cn–k – Cn dn–k ≤ 0, где n = 5, то очевидно
C0 – C5 d0 ≤ 0; C2 – C5 d2 ≤ 0; C4 – C5 d4 ≤ 0,
122
и все приведенные неравенства справедливы, то для получения эквивалентного подграфа сигналов (без доказательства) необходимо в нем все знаки коэффициентов передачи прямых путей изменить на положительный, при этом y будет иметь отрицательный знак при нечетном n. Подграф сигналов в этом случае показан на рис. 6.1, б. Число усилителей, реализующих данный подграф сигналов, Z = n + 1.
В случае, если порядок числителя меньше порядка знаменателя передаточной функции, т.е. нет пути Pn (C5 = 0), то подграф сигналов строится на общих основаниях, где общее число усилителей равно n + 2, так как не требуется усилитель, реализующий соприкасаемость пути Pn = C5 сконтурамиподграфа сигналов.
Пример 4. Составить электронную модель реального дифференцирующего звена:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
W ( p) = − |
Tp |
|
= − |
Tp + 0 |
|
= − |
|
T1 |
|
, |
||||||||
|
T1 p + 1 |
T1 p + 1 |
1 |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
+ 1 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T1 p |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где C0 |
= 0; C1 |
= |
T |
; |
d0 = |
1 |
|
; n = 1; |
k = 1, |
|
|
|
|||||||
|
T p |
|
|
|
|||||||||||||||
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
Cn–k – Cn dn–k < 0 или C0 – C1 d0 < 0.
На рис. 6.5 приведены графы сигналов и его электронная версия.
В случае, если необходимо ввести масштаб по времени, то необходимо оператор p передаточной функции заменить на Pм/Mt, где Pм – машинный оператор. Если Mt > 1, то масштаб по времени – ускоряющий, если Mt < 1, то масштаб по времени замедляющий. При Mt = 1 процессы в электронной модели будут протекать в реальном времени.
Если значение коэффициентов схемы модели оказывается за пределами возможностей усилителя, то необходимо перераспределить коэффициенты согласно правилам преобразования графа сигналов.
123
Рис. 6.5. Граф сигналов реального дифференцирующего звена и его электронная версия
Таким образом, данный метод моделирования отличается простотой и наглядностью. Здесь не требуется предварительных преобразований уравнений. Метод графа сигналов позволяет вскрыть структуру САУ, что способствует более ясному пониманию существа процессов, протекающих в системе.
Таким образом, можно любую сложную САР представить аналитически в виде передаточной функции.
6.2. АДАПТИВНЫЕ СИСТЕМЫ С ЖЕСТКОЙ ЭТАЛОННОЙ МОДЕЛЬЮ
1. Построение адаптивной системы с применением жесткой эталонной модели путем сравнения передаточных функций модели W ( p)м и передаточной функции САР Ф( p) . На рис. 6.6
приведена структура адаптивной САР.
Методомграфасигналов (см. рис. 6.6) определяемпутиграфа:
P1 = W ( p)p W ( p)об ;
P2 = W ( p)M K W ( p)oc W ( p)p W ( p)об.
Контур графа сигналов:
L2 = −W ( p)p W ( p)об K W ( p)oc .
124
Рис. 6.6. Структура САР с эталонной моделью
Миноры путей графа сигналов:
1= 1 ;
2= 1.
Общий определитель графа сигналов:
= 1− L1 .
Определяем передаточную функцию САР согласно формуле Мезона:
|
δ |
|
W ( p) |
W ( p) |
об |
+ W ( p) |
|
K W ( p) |
|
W ( p) |
W ( p) |
|
|
Φ( p) = |
= |
p |
|
|
M |
|
oc |
p |
|
об |
. |
||
|
|
|
+W ( p)p W ( p)об K W ( p)oc |
|
|
||||||||
|
δзад |
|
1 |
|
|
|
При условии K >> 1 очевидно
(6.1)
Таким образом, при большом коэффициенте в канале обратной связи K при любых изменениях параметров объекта будет выполняться равенство (6.1).
2. Построение адаптивной системы с применением жесткой эталонной модели путем сравнения амплитудно-частотных
характеристик модели W ( p)м и Ф( p) .
На рис. 6.7 приведена структура адаптивной САР.
125
Рис. 6.7. Структурная схема адаптации с помощью сравнения амплитудно-частотных характеристик модели и САР:
Г – генератор синусоидальных колебаний; В1 и В2 – выпрямители; ЭС1 иЭС2 – элементы сравнения; МН – механизм настройки, ЭМ – эталонная модель; СФ – синхронный фильтр
Амплитудно-частотная характеристика (АЧХ) эталонной модели САР (выход выпрямителя В1 ) сравнивается с АЧХ пер-
вой гармоники САР (выход выпрямителя В2 ). Сигнал разности
подается на вход механизма настройки (МН), который корректирует настройку пропорциональной части закона регулирования с целью придания адаптивности САР.
6.3. МЕТОД АКТИВНОЙ АДАПТАЦИИ
Рассмотрим самонастраивающуюся систему с «гибкой» эталонной моделью, где оценка модели объекта может выполняться спомощью пассивного илиактивногоэксперимента.
При пассивном эксперименте используются для обработки случайные возмущения, для снятия которых в необходимом объеме требуется достаточно большое время.
При активном эксперименте обычно на объект подается дополнительное сигнальное идентифицирующее воздействие (единичное воздействие, синусоидальное воздействие и т.д.), либо дается параметрическое воздействие (меняются настройки в регуля-
126
торах), либо алгебраическое нелинейное воздействие (включение нелинейного элемента), либо структурное воздействие (образование новых замкнутых контуров), либо введение в систему квантования сигналов повремении поуровню.
Пусть за основу принимается оценка модели объекта путем подачи на его вход синусоидального воздействия в определенные интервалы времени.
На рис. 6.8 показана структурная схема адаптивной одноконтурной САР, где блок идентификации (БИ) и блок оптимизации (БО) обеспечивают параметрическую адаптацию САР.
Рис. 6.8. Структурная схема адаптивной САР
В процессе активной адаптации САР БИ постоянно формирует модель объекта, максимально приближенную к реальному объекту в его диапазоне существенных частот.
БО позволяет рассчитывать оптимальные настройки регулятора с корректированной модели объекта с дальнейшей установкой их на регуляторе.
Организация процедуры идентификации объекта
Формирование синусоидального воздействия, оценка модуля и фазы замкнутой САР и объекта
Частота вынужденных колебаний ωГ выбирается экспери-
ментально (из опыта эксплуатации) и реализуется отдельным генератором с постоянной амплитудой выходных колебаний.
127
Оценка модуля и фазы вектора замкнутой САР выполнена корреляционнымметодомпоследующимрекуррентнымформулам:
U К =U К−1 + |
|
2 |
tYK |
; |
|
|||||
(T A )sin |
2πtK |
|
||||||||
|
|
|
|
|
|
|
||||
|
|
|
|
|||||||
|
|
|
|
|
Г |
K |
|
TГ |
|
|
|
|
|
|
|
|
|
|
|
|
|
V K =V |
К−1 + |
|
|
|
|
2 |
tYK |
|
, |
|
(TГ AK )cos(2πtK TГ ) |
||||||||||
где tk = k t ; k 1,n ; |
N = |
TГ |
– число интервалов дискретности |
|||||||
|
||||||||||
|
|
|
t |
|
|
|
|
|
на каждом периоде колебаний генератора;
Значение величин Uk и Vk , полученных при k = N, принимается в качестве текущих оценок U (ω)L и V (ω)L . По этим
данным вычисляются текущие оценки модуля и фазы вектора Φ( jω) замкнутой САР по окончании каждого L периода коле-
баний генератора TГ .
A(ω)L = U 2 (ω)L + V 2 (ω)L ,
φ(ω)L |
= arctg |
V (ω)L |
. |
|
|||
|
U (ω)L |
Далее вычисляем средние арифметические значения (математическое ожидание) оценок модуля и фазы за все L закончившихся периодов колебаний по рекуррентным формулам:
A(ω)L = A(ω)L−1 + {A(ω)l − A(ω)L−1} / L, φ(ω)L = φ(ω)L−1 + {φ(ω)L − φ(ω)L−1} / L.
Необходимая длина реализации (число периодов L) определяется путем сравнения математического ожидания оценки A(ω)L сдопустимымзначением εA спомощью следующегоусловия:
128
DL |
≤ εA , |
(6.2) |
A(ω)опт |
где A(ω)опт – заданное оптимальное значение A(ω) , а
|
= |
1 |
L |
|
|
|
|
2 |
|
|
|
|
|
||||
DL |
|
|
|
|
|
|||
|
|
|
||||||
(L − 1) |
A(ω)i |
− A(ω)L |
|
|||||
|
|
i=1 |
|
|
|
|
|
является оценкой дисперсии случайной величины A(ω)L .
Если условие (6.2) не выполняется, то текущее число периодов считается недостаточным и процесс определения оценок A(ω)L и φ(ω)L продолжается в течение следующих дополни-
тельных периодов колебаний.
Как только условие (6.2) будет выполнено, значения A(ω)L и φ(ω)L , полученные за прошедшие L периодов, принимаются в качестве искомых оценок модуля A(ω) и фазы φ(ω)
замкнутой САР.
По установившимся колебаниям регулируемой величины оцениваются модуль R(ω)об и фаза φ(ω)об вектора
АФЧХ объекта для текущей частоты эксперимента с помощью формулы
|
|
W ( jω)об = |
Ф( jω) |
, |
(6.3) |
||
|
|
W ( jω)p [1− Ф( jω)] |
|||||
где Ф( jω) = A(ω)e jφ(ω) , |
|
|
|||||
A(ω) = |
|
(ω)L – АЧХ замкнутой САР; |
|
|
|||
A |
|
|
|||||
φ(ω) = φ(ω)L – ФЧХ замкнутой САР; |
|
|
|||||
2 |
|
C0 |
|
||||
|
|
W ( jω)P = C12 + |
C0 |
e− jarctg |
|
|
|
|
|
ωC1 |
; |
||||
|
|
|
|||||
|
|
|
|
ω2 |
|
|
– АФЧХ ПИ-регулятора.
129
Тогда согласно (6.3), получим следующие выражения:
R(ω)об = |
|
|
|
|
|
|
A(ω) |
|
; |
(6.4) |
||
|
|
C 2 |
|
|
|
|
||||||
2 |
|
|
2 |
|
|
|
||||||
|
0 |
|
|
|
|
|
|
|||||
|
C1 |
+ |
|
|
|
|
1− 2A(ω)cosφ(ω) + A |
(ω) |
|
|||
ω2 |
|
|
||||||||||
φ(ω)об = arctg |
|
C0 |
+ arctg |
sin φ(ω) |
|
; |
|
(6.5) |
||||
C1ω |
cosφ(ω) − A(ω) |
|
||||||||||
|
|
|
|
|
|
|
|
|||||
где R(ω)об – АЧХ объект регулирования; φ(ω)об |
– ФЧХ объект |
|||||||||||
регулирования. |
|
|
|
|
|
|
|
|
|
|
|
Принятая модель объекта имеет рациональную переда-
точную функцию вида |
|
|
|
|
|
|
W ( p)мод = |
K |
a |
e− pτa |
, |
(6.6) |
|
(Ta p + 1)n |
||||||
об |
|
|
где Ka – коэффициент передачи; Ta |
– постоянная времени; τa – |
|||||||
чистое запаздывание; n – порядок полинома знаменателя. |
|
|||||||
Запишем модуль R(ω)обмод и фазу φ(ω)обмод из (6.6), заменив |
||||||||
р = jω: |
|
|
|
|
|
|
|
|
R(ω)обмод = Ka e− ωτa (1− ωTa )2 |
|
+ ω2Ta 2 −0,5n , |
(6.7) |
|||||
|
|
|
|
|
|
|
||
φ(ω)обмод = −ωτa − n arctg |
|
ωTa |
. |
(6.8) |
||||
|
|
|
||||||
|
|
|
|
1− ωTa |
|
|||
– Приравняв выражения (6.4) и (6.7), а также (6.5) и (6.8), |
||||||||
получим систему двух уравнений: |
|
|
|
|
|
|
|
|
R(ω)об = Ka e− ωτa (1− ωTa )2 |
+ ω2Ta 2 −0,5n , |
(6.9) |
||||||
|
|
|
|
|
|
|
||
φ(ω)об = −ωτa − n arctg |
|
|
|
ωTa |
. |
(6.10) |
||
1 |
|
|||||||
|
|
− ωTa |
|
130