Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Современные проблемы теории управления

..pdf
Скачиваний:
9
Добавлен:
15.11.2022
Размер:
6.12 Mб
Скачать

означает отсутствие соответствующей ветви подграфа. При 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

Φ( p) W ( p)M .

Рис. 6.6. Структура САР с эталонной моделью

Миноры путей графа сигналов:

1= 1 ;

2= 1.

Общий определитель графа сигналов:

= 1L1 .

Определяем передаточную функцию САР согласно формуле Мезона:

 

δ

 

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

tK

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

K

 

TГ

 

 

 

 

 

 

 

 

 

 

 

 

V K =V

К1 +

 

 

 

 

2

tYK

 

,

(TГ AK )cos(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(ω)L1 + {A(ω)l A(ω)L1} / L, φ(ω)L = φ(ω)L1 + {φ(ω)L φ(ω)L1} / 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

ejarctg

 

 

 

 

ωC1

;

 

 

 

 

 

 

 

ω2

 

 

– АФЧХ ПИ-регулятора.

129

Тогда согласно (6.3), получим следующие выражения:

R(ω)об =

 

 

 

 

 

 

A(ω)

 

;

(6.4)

 

 

C 2

 

 

 

 

2

 

 

2

 

 

 

 

0

 

 

 

 

 

 

 

C1

+

 

 

 

 

12A(ω)cosφ(ω) + A

(ω)

 

ω2

 

 

φ(ω)об = arctg

 

C0

+ arctg

sin φ(ω)

 

;

 

(6.5)

C1ω

cosφ(ω) A(ω)

 

 

 

 

 

 

 

 

 

где R(ω)об – АЧХ объект регулирования; φ(ω)об

– ФЧХ объект

регулирования.

 

 

 

 

 

 

 

 

 

 

 

Принятая модель объекта имеет рациональную переда-

точную функцию вида

 

 

 

 

 

W ( p)мод =

K

a

epτ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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]