Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
metods / Управление в системах.pdf
Скачиваний:
112
Добавлен:
26.03.2015
Размер:
676.4 Кб
Скачать

Для рассмотренного выше примера найдем передаточную функцию замкнутой системы по ошибке от задающего воздействия:

 

 

 

Фx (z)

=

(z 1)(z d)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(z 1)(z d) + kγT0

(1 d)z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для g(t)=g0 1(t) с учетом табл. 1 и (6.1) получим:

 

 

 

 

 

 

xус т

= lim

(z 1)

 

 

(z 1)(z d)

 

 

 

 

z

 

= 0 .

 

z

[(z 1)(z d) + kγT0 (1 d)z] (z 1)

 

 

z1

 

 

 

Данная система обладает астатизмом по задающему воздействию.

 

Для g(t)=vt1(t) с учетом табл. 1 и (11.7) получим:

 

 

 

 

 

 

xус т =

lim

(z 1)

 

 

 

(z 1)(z d)

 

 

 

vT0 z

 

=

v

.

 

z [(z 1)(z d) + kγT0 (1

d)z] (z 1)

2

kγ

 

z1

 

 

 

Лекция 12. Особенности построения математических моделей компьютерных систем управления. Методы линеаризации

В системах компьютерного управления, или цифровых автоматических системах (ЦАС), кроме квантования по времени выполняется квантование сиг- нала по уровню. Таким образом, ЦАС являются нелинейными импульсными системами. Достаточно точный анализ процессов в таких системах может быть выполнен только путем имитационного моделирования. Однако для проекти- руемой системы такое моделирование невозможно без предварительного ана- литического исследования. На основе достаточно простых аналитических ме- тодов должны быть определены структура системы, ее основные параметры и показатели качества, выбраны элементы системы.

Рассмотренные выше методы учета квантования по времени в рамках теории импульсных систем могут применяться для ЦАС, если характерный для

них эффект квантования по уровню удается учесть в рамках линеаризованной модели. На основе такой модели выполняется основной объем работ по анали- зу и синтезу системы. Уточнение результатов с учетом влияния квантования по уровню выполняется на заключительном этапе на основе имитационного моде- лирования или соответствующих аналитических оценок.

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

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

45

Предположим, что для нелинейного звена (рис. 31) получена модель в виде нелинейного дифференциального уравнения:

æ

.

.

..

ö

æ

. ö

 

Fç x1, x 2

, x 2

, x 3, x 3

, x

3 ÷

= jç f , f ÷ ,

(12.1)

è

 

 

 

ø

è

ø

 

где x1 и x2 - входные сигналы, x3 - выходной сигнал, f - возмущение.

Допустим, установившийся процесс в сис-

теме имеет место при некоторых постоянных значениях x1 = x10 , x 2 = x20 , x 3 = x30 и

f = f 0 . Тогда для установившегося процесса уравнение (12.1) примет вид:

F(x10 , x20 ,0, x 30 ,0,0) = ϕ(f 0 ,0). (12.2)

В основе линеаризации лежит допущение о том, что в исследуемом динамическом процессе переменные, характеризующие состояние системы (x1 , x2 , x3 ), являются не-

прерывными функциями времени и отклоняются от установившихся значений (x10 , x 20 , x30 ) достаточно мало. Предполагается также, что F - непрерывная и

дифференцируемая по всем аргументам функция. Обозначим эти отклонения x1 , x2 , x3 . Тогда в динамическом процессе будем иметь:

x1 = x10 + Dx1, x 2 = x20 + Dx2 , x 3 = x30 + Dx3 ,

.

. .

. ..

..

x 2

= D x 2 , x 3

= D x 3 , x 3

= D x 3 .

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

Разложим нелинейную функцию из правой части (12.1) в ряд Тейлора по всем ее аргументам. Тогда уравнение (12.1) с учетом первых членов разложе- ния примет вид:

F(x 0

 

,0, x 0 ,0,0)+

æ

 

F

ö

0

 

 

æ

F

ö

0

 

 

æ

 

F

ö 0 .

 

, x 0

ç

 

÷

Dx

 

+ ç

 

÷

Dx

 

+

ç

 

÷

D x 2

+

 

 

 

 

 

 

 

 

.

 

1

2

 

3

 

 

è

x1 ø

 

 

1

è

x 2 ø

 

2

 

ç

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ö 0

 

 

 

 

 

 

 

0

è

 

x 2

ø

 

 

 

 

æ

F ö

0

 

 

 

 

æ

 

 

F

.

 

æ

F

ö

 

..

 

 

 

 

 

 

 

 

 

 

 

ç

 

 

÷

 

ç

÷

 

 

 

 

 

 

 

 

 

÷

Dx

 

 

+

 

 

 

 

 

 

D x 3

+

 

 

 

 

D x 3

+

 

 

 

 

 

 

3

ç

 

 

.

÷

ç

..

÷

 

 

 

 

 

è

x 3 ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è ¶ x 3

ø

 

 

 

è

x 3

ø

 

 

 

 

 

 

 

 

 

 

æ

 

 

Dx 2,

 

.

 

 

 

 

 

 

.

 

 

..

ö

 

æ

 

.

ö

 

 

 

 

+0çDx1,

D x 2, Dx 3, D x 3,

D x

3 ÷ = jç f , f ÷ ,

 

(12.3)

 

 

è

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ø

 

è

 

 

 

ø

 

 

 

46

æ F ö 0

где ç ÷ - значение частной производной в точке, соответствующей устано-

è ¶x ø

вившемуся процессу, 0(...) - члены разложения более высокого порядка мало- сти.

После вычитания (12.2) из (12.3) получим уравнение:

æ

F

ö

0

 

æ

F

ö

0

 

 

 

æ

F

 

ö

0

.

æ

F

ö

0

 

 

 

 

 

 

 

 

 

 

 

 

ç

 

÷

 

 

 

 

 

 

 

ç

 

 

÷

Dx

 

+ ç

 

÷

 

Dx

 

 

+

 

 

 

 

D x 2 +

ç

 

 

÷

 

Dx

 

+

 

 

1

x 2

2

ç .

 

÷

x 3

 

 

3

è

x1 ø

 

è

ø

 

 

 

 

 

 

è

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

x 2

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

F

ö

0

.

æ

 

F

 

ö

0

 

..

 

 

æ

. ö

 

 

 

 

0

 

)

 

 

 

 

 

ç

 

÷

 

ç

 

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

D x 3 +

 

 

 

 

 

 

D x 3

= jç f , f ÷

- j

f

 

,0

,

 

 

ç .

÷

ç

..

 

÷

 

 

 

 

 

 

 

 

 

 

 

 

 

è

ø

 

(

 

 

 

 

 

 

 

 

 

è ¶ x 3

ø

 

 

è

 

x 3

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которое дальше может быть записано в виде, позволяющем перейти к переда- точным функциям:

..

 

.

 

 

 

.

 

 

 

 

 

 

a0 x 3

+ a1D x 3

+ a2 Dx 3

= b1Dx1 + b1' Dx 2

+ b2' D x 2

+ f1.

(12.5)

 

 

 

 

 

 

 

 

 

æ

F

ö 0

 

Получено

линеаризованное уравнение,

где

a

=

ç

÷

,

..

 

 

 

 

 

 

 

0

 

ç

÷

 

 

 

 

 

 

 

 

 

 

è

x 3

ø

 

 

æ

F

ö

0

æ

F

ö

0

a = ç

÷

 

 

, a = ç

÷

, b

 

 

1

ç .

÷

2

è

x 3 ø

1

 

è ¶ x 3

ø

 

 

 

 

 

- постоянные коэффициенты,

æ

ö

0

 

 

æ

F

ö

0

 

æ

ö 0

= -ç

F

÷

 

, b ' =

÷

 

, b ' = -ç

F

÷

 

 

x 2

 

 

è

x1 ø

 

 

1

è

ø

 

2

ç

. ÷

 

 

 

 

 

 

 

 

 

è ¶ x 2 ø

 

æ .

 

- j(f 0

,0)

 

 

 

 

ö

-

эквивалентное воз-

f1 = jç f , f

÷

 

è

 

 

ø

 

 

 

 

 

 

 

 

 

мущающее воздействие. Символы приращения при дальнейшем использова- нии уравнений вида (12.5) отбрасывают.

Из полученных соотношений видно, что рассмотренный метод линеари- зации может применяться только при условии, что функция F не имеет разры- вов первого или второго рода. Для входящих в состав компьютерной системы управления ЦАП и АЦП это условие не соблюдается. Поэтому при математиче- ском описании систем компьютерного управления рассмотренный метод ли- неаризации может применяться только для их непрерывной части. Распростра- нить его на часть системы, где выполняется квантование сигнала по уровню, невозможно.

Для линеаризации ступенчатых нелинейных характеристик ЦАП и АЦП применяется метод, основанный на допущении, противоположном использо-

47

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

Аналого-цифровое преобразование состоит в формировании двоичного кода в зависимости от величины непрерывного входного сигнала. Если на вы- ходе АЦП формируется l-разрядный двоичный код, то количество уровней его релейной статической характеристики будет равно 2l . ЦАП выполняет обрат- ное преобразование. Количество уровней его характеристики определяется раз- рядностью входного двоичного кода аналогичным образом.

На рис. 32а показана начальная часть статической характеристики АЦП, стоящего на входе компьютерной системы управления (например, рис. 20). По горизонтальной оси отложено непрерывное значение задающего воздействия g, по вертикальной оси - его цифровое представление g0 (в десятичной системе счисления). Величина единицы младшего разряда обозначена d1 . Обычно АЦП

строят таким образом, что величина d1 постоянна для всего диапазона возмож- ных значений преобразуемого сигнала. Следовательно, характеристику АЦП можно описать зависимостью:

 

é

g

 

1

ù

 

g0

= Eê

 

+

 

sign gú

,

d

2

 

ë

1

 

 

û

 

где E[...] - целая часть числа.

48

Линеаризация проводится путем замены ступенчатой характеристики ли- нейной, как показано на рис. 32а. Это означает, что АЦП рассматривается как безынерционное звено с передаточной функцией W(p)=k1, где коэффициент передачи определяется через величину младшего разряда преобразователя:

k1 =

1

.

(12.6)

 

 

δ1

 

Наибольшая ошибка от замены нелинейной характеристики линейной не

будет превышать 0,5δ1 .

Аналогично выполняется линеаризация и для аналого-цифрового преоб- разователя управляемой величины y, статическая характеристика которого по- казана на рис. 32б. Коэффициент передачи линеаризованного звена будет ра- вен:

k2 =

1

,

(12.7)

 

 

δ2

 

где δ2 - величина единицы младшего разряда преобразователя. Ошибка линеа- ризации определяется аналогичным образом.

Чаще всего задачей системы является обеспечение равенства y=g. В та-

ких случаях выбирают δ1=δ2, и рассмотренные АЦП при составлении модели могут быть условно объединены в один преобразователь - преобразователь сигнала ошибки x=g-y. Его характеристика (рис. 32в) оказывается неодно- значной. Но линеаризация для нее, выполненная рассмотренным способом, да- ет однозначный результат - безынерционное звено с коэффициентом передачи

(12.6).

На рис. 32г показана статическая характеристика ЦАП, стоящего на вы- ходе цифровой части системы. По горизонтальной оси отложена выходная ве- личина цифровой части (цифровое представление сигнала управления), по вер- тикальной оси - сигнал управления, преобразованный в непрерывную форму (чаще всего - электрическое напряжение). Единица младшего разряда преобра- зователя обозначена δ. При линеаризации ЦАП также заменяется безынерци- онным звеном, коэффициент передачи которого будет равен

k=δ.

(12.8)

Рассмотрим подробнее вопрос о погрешностях описанных моделей ЦАП и АЦП.

Для преобразователей со нелинейными статическими характеристиками, показанными на рис. 32а, б, г, может рассматриваться эквивалентная схема (рис. 33). Здесь звено 1 является линейным, с коэффициентом передачи (12.6), (12.7) или (12.8). Звено 2 обеспечивает учет погрешности линеаризованной ха- рактеристики преобразователя относительно реальной ступенчатой характери- стики. Звено 3 имеет коэффициент передачи, равный единице, и отражает ог- раничение на величину преобразуемого сигнала, связанное с конечным количе-

49