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

GOS for Iphone / mobile / ТАУ / Книги / ТАУ учебник

.pdf
Скачиваний:
228
Добавлен:
10.05.2015
Размер:
4.33 Mб
Скачать

 

 

ШИМ

θ (t)

f(x)

x(t)

x*(t)

 

 

 

W(ck,p)

 

 

 

x

 

 

 

t

0

x* τ 1 = k1T τ 2 = k2T A

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

–A

Рис. 4.1

[7, 9], включив последовательно с импульсным элементом линейный формирующий элемент с необходимой передаточной функцией, несущим импульсам всегда можно придать любую, в том числе и прямоугольную форму.

В общем виде уравнение движения САУ с ШИМ записывается следующим образом:

Q (ck , D ) x (t ) + Q

 

 

 

 

(t ),

(4.2)

 

(ck , D ) x

(t ) = S (ck , D ) f (t ) + S

(ck , D ) f

где x(t), x*(t) – сигналы на входе и выходе широтно-импульсного модулятора, соответственно, относительно которых записано уравнение динамики синтезируемой системы управления; f(t), f*(t) – внешнее входное воздействие на входе и выходе ШИМ, соответственно; Q(ck,D), Q*(ck,D), S(ck,D), S*(ck,D) – полиномы оператора обобщенного дифференцирования D соответствующих степеней с вещественными постоянными коэффициентами [см. формулу (2.17)].

121

В соответствии с общей схемой решения задачи синтеза систем управления, рассмотренной во втором разделе, исходя из уравнения (4.2) получаем целевую функцию

m n

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

J =

 

 

ai (ck ) Aqi

+

 

 

 

 

 

ei (ck)Cqi +

 

 

 

 

 

 

 

 

ai

(ck )Aqi

 

ei

(ck )Cqi , (4.3)

 

 

 

 

 

 

 

 

 

 

i=0

 

 

 

 

 

i=0

i=0

 

 

 

 

 

 

 

q=1 i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

minc

J

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

−ρqt

 

 

 

 

 

D

i

A

 

 

 

 

 

 

1(t

nT ) 1(t (n + kn )T )

 

e

dt,

 

Aqi =

 

 

sign (xn )

 

 

 

 

 

 

 

 

n=

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i = 0,1,…, n

,

q = 1,2,…, m,

 

 

 

 

 

 

 

(4.4)

 

 

 

 

 

 

 

 

 

0

)

 

 

 

 

 

 

 

 

−ρqt

 

 

 

=

 

D

i

A

 

 

 

 

1(t nT ) 1(t (n + kn )T )

 

e

dt,

 

Cqi

 

 

 

sign ( fn

 

 

 

 

 

 

 

 

n=

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i = 0,1,…, v

,

q = 1,2,…, m,

 

 

 

 

 

 

 

(4.5)

а интегралы Aqi, Cqi были определены ранее.

Таким образом, для решения задачи параметрического синтеза САУ с ШИМ необходимо определить интегралы вида (4.4), (4.5) и далее минимизировать функционал J в соответствии с общей схемой решения поставленной задачи.

Прежде чем рассматривать вычисление интегралов (4.4), (4.5) необходимо отметить следующее. Интегралы Aqi, Cqi вычислялись при определенном (заданном) виде программного движения и внешнего входного воздействия, поскольку вид аналитических рекуррентных соотношений непосредственно связан с видом желаемого процесса и входного воздействия. В случае САУ с ШИМ модулятор преобразует входной сигнал любого вида в последовательность прямоугольных импульсов одинаковой амплитуды, следовательно, вид аналитических выражений, определяющих Aqi*, Cqi*, будет одинаковым. Поэтому рассмотрим лишь вычисление интеграла Aqi*.

122

При i = 0 получаем

 

 

 

0

 

 

 

 

−ρqt

 

=

 

 

A

 

1(t nT ) 1(t (n + kn )T )

 

e

dt =

Aq0

 

sign (xn )

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

= A

 

 

 

 

sign (xn0 )

e−ρqtdt e−ρqtdt

=

 

 

 

 

(n+kn )T

 

 

n=0

nT

 

 

 

= A

sign (xn0 )1 e−ρqknT e−ρqnT .

 

(4.6)

 

n=0

 

 

ρq

 

 

Далее, применяя к полученному соотношению (4.6), формулу (3.58), определяющую сумму членов геометрической прогрессии, получим

 

 

 

 

 

 

Aq

0 =

A

 

0

 

n=0

0

 

 

 

 

 

 

−ρqt

 

sign (xn )

1(t nT ) 1(t (n + kn )T )

 

e

 

dt =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= A

sign (xn0 )

1 e−ρqknT

.

 

 

 

(4.7)

 

 

 

 

n=0

 

 

ρq (1 e−ρqT )

 

 

 

 

Применяя аналогичные выкладки, определим интеграл Aqi* при i = 1

 

 

 

 

Aq1 =

A

0

 

n=0

0

 

 

 

 

−ρqt

 

 

sign(xn )

δ(t nT ) δ(t (n + kn )T )

 

e

 

dt.

(4.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Преобразуя соотношение (4.8), аналогично вычисленному ранее (4.6), и используя фильтрующее свойство δ -функции

f (t )δ(k ) (t τ)dt = (1)k f (k ) (τ),

0

получаем следующее:

Aq1 = A

sign(xn0 )1e−ρqknT .

(4.9)

n=0

1e−ρqT

 

123

В случае i = n* получаем

 

 

 

n

 

 

 

 

 

 

0

 

 

 

 

 

 

 

−ρqT

 

 

A

 

=

D

 

 

 

 

Asign(xn ) 1(t nT )1(t (n + kn )T )

 

e

 

 

dt =

qn

0

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−1)

 

 

 

−1)

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

(n

(t nT )δ

(n

(t

 

 

 

 

 

 

A

sign(xn0 )

δ

 

 

 

 

 

(n + kn )T )

e−ρqT dt. (4.10)

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее, приводим соотношение (4.10) к виду

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

0

 

 

 

 

 

 

 

 

−ρqT

 

A

 

=

D

 

 

 

A

sign(xn ) 1(t

nT )1(t (n + kn )T )

e

 

 

dt =

qn

0

 

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= A

sign(xn0 )1e−ρqknT

ρnq −1.

 

 

 

 

 

(4.11)

 

 

 

 

 

 

 

 

 

 

n=0

 

 

1e−ρqT

 

 

 

 

 

 

 

Таким образом, обобщая полученные результаты, получаем следую-

щее рекуррентное выражение, определяющее интеграл Галеркина A

*

 

 

 

 

 

 

 

 

qi

при наличии в САУ широтно-импульсного модулятора [188, 189]

 

 

0

 

 

−ρqT

i

 

 

i

A

 

e

 

Aqi = D

 

sign(xn ) 1(t nT )1(t (n + kn )T )

 

 

dt = Aqρq,

 

0

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

i = 0,1,…,n ;

q = 1,2…,m,

 

(4.12)

где

 

 

 

 

 

Aq = A

sign (xn0 )

1 e−ρqknT

 

(4.13)

(1 e−ρqT )ρq

n=0

 

 

 

 

здесь A – амплитуда импульса на выходе модулятора; sign (xn0) – знак сигнала на входе модулятора в n моменты квантования.

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

124

ствия будет использоваться скачкообразное или линейно нарастающее, т. е. функции постоянные по своему знаку, то формулы (4.12), (4.13) могут быть упрощены

Cqi = Di

0

 

 

A

 

 

n=0

 

0

)

 

 

 

−ρqT

i

sign( fn

1(t − nT ) −1(t − (n + kn )T )

 

e

 

dt = Cqρq,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, q =

 

(4.14)

i = 0,1,

,v

 

1,2m

 

где

 

 

 

 

 

 

 

 

Cq = A

 

 

1 e−ρqknT

 

 

 

 

(1 e

−ρqT

)ρq

(4.15)

 

n=0

 

 

 

 

 

 

 

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

Задача определения sign(xn0) может быть решена довольно просто. Поскольку известно желаемое программное движение x0(t), время переходного процесса Tп.п и период квантования, то для моментов времени, соответствующих n-ым периодам квантования в интервале от 0 до Tп.п определяются значения процесса. По найденным числовым значениям процесса x0(t) формируется массив sign(xn0) в моменты коммутации входного сигнала.

Ширина импульса на выходе широтно-импульсного модулятора определяется значением kn, которое нелинейно зависит от соотношения величин входного x0(t) и опорного y(t)=sign (xn0) t сигналов (рис. 4.1). В качестве опорного необходимо использовать тот сигнал, который применяется в реальном ШИМ проектируемой системы управления. Как правило, в схемах широтно-импульсных модуляторов опорным является сигнал пилообразной формы с периодом T. Исходя из этого, процедура определения массива значений kn может быть следующей: в начале задается приращение ∆ t координаты t и определяется значение разности процессов x0(t) и y(t): ε = x0 (t)− y (t), когда указанная разность не будет превышать заданное значение ε 0, соответствующее значение tj,

запоминается и по нему определяется kn

=

t j

. Далее происходит обну-

 

 

 

T

125

Начало

Ввод исходных данных

 

 

Вычисление

 

 

x0(tj ), x0(tj)

 

 

Вычисление

 

 

 

y(t*j ) = k t*j

 

 

Вычисление

 

ε =

x0(tj) − y(tj)

 

 

 

Нет

 

 

 

ε ≤ ε 0

 

 

 

Д а

 

Формирование масси

 

 

ва значений kj

 

Формирование мас-

 

сива знаков амплитуд

t*j = 0

 

tj+ 1 = tj + ∆ t

 

Нет

tj+ 1 ≤ Tп.п

 

 

 

 

 

Да

 

 

Печать массива

 

 

знаков амплитуд

 

 

и

kj

 

 

 

 

Конец

Рис. 4.2.

Нет

Да

tj+ 1 ≤ T п.п

tj+ 1

=

 

tj + ∆

t

 

t*j

= tj + ∆

t

ление опорного сигнала, и вычисления повторяются до тех пор, пока текущее значение координаты t не превысит заданную длительность

126

переходного процесса. Таким образом, формируется массив значений

kn на интервале от 0 до Tп.п.

Блок-схема алгоритма программы, реализующей данную процедуру, представлена на рис. 4.2.

4.2. САУ с частотно-импульсной модуляцией

Область применения систем автоматического управления с частотноимпульсными модуляторами (ЧИМ) весьма широка. Модуляторы данного вида используются как в системах регулирования инерционных технологических процессов, так и в системах ориентации и стабилизации летательных аппаратов [9]. Системы управления данного класса являются непериодическими и существенно нелинейными в силу принципа частотного управления, и во многих режимах работы не допускают линеаризации даже при малых глубинах модуляции. Отмеченные особенности САУ с ЧИМ затрудняют решение задачи синтеза параметров операторов управления традиционными методами. Поэтому подробно рассмотрим применение обобщенного метода Галеркина к системам данного класса, общая схема которого для решения задачи синтеза импульсных САУ с различными видами модуляции сигнала была показана во втором разделе.

Задача синтеза решается в традиционной для обобщенного метода Галеркина постановке: известна структура САУ с ЧИМ, формирующим сигналы, частота следования которых зависит от сигнала на входе модулятора, и параметры объекта управления. Требуется определить параметры регулятора, структура которого задана в общем виде, исходя из условия приближенного обеспечения в синтезируемой системе заданных показателей качества ее работы в переходном режиме, абсолютной устойчивости и грубости САУ по варьируемым параметрам.

Уравнение движения САУ n-го порядка с ЧИМ, структурная схема которой показана на рис. 4.3, а и б имеет вид (4.2).

Следовательно, в соответствии с общей схемой решения задачи синтеза САУ с различными видами модуляции сигнала, получаем целевую функцию

 

 

n

 

 

 

 

n

 

 

 

 

v

 

 

 

v

 

 

2

 

J =

 

a

(c

) A

+

a

(c

e

(c

)C

 

(c

 

, (4.16)

 

 

) A

e

)C

qi

 

i

k

qi

i

k

qi

 

i

k

qi

 

i

k

 

 

 

 

 

 

 

 

 

i=0

 

 

 

 

i=0

 

 

 

i=0

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

minck J → 0,

127

а)

 

 

 

 

x * ( t )

 

 

 

 

f ( t )

 

 

x ( t )

ЧИМ

 

W ( c k , p )

θ ( t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

x *

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T 1

 

 

T 2

T 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

– A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 4.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qi

 

 

 

 

{

 

 

 

 

 

 

}

ρqt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

=

 

 

 

D

i

 

x

0

(t)

e

dt,

i = 0,1,

 

 

 

 

,m;

(4.17)

 

A

 

 

 

 

 

 

 

 

 

 

 

 

,n ; q = 1,2,

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{

 

 

 

(t)

 

ρqt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

=

D

i

f

e

dt, i =

 

 

 

 

 

q = 1,2,…,m,

 

 

 

 

 

qi

 

 

 

 

 

 

0,1,…,v ;

 

(4.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

здесь x0*(t) и f*(t) – процессы на выходе частотно-импульсного модулятора, определяемые в соответствии с (2.10)

x0 (t) = A

sign(x0 (Tn )) 1(t − Tn )−1(t − Tn − τ ) ;

(4.19)

n=0

 

 

128

f (t) = A

sign( f (Tn )) 1(t Tn )1(t Tn τ ) ,

(4.20)

n=0

 

 

где A = const – амплитуда импульса, длительностью τ = const на выходе модулятора; Tn – период следования импульсов, величина которого зависит от сигнала на входе модулятора; sign(x0(Tn)), sign(f(Tn)) – знаки желаемого программного движения x0(t) и внешнего воздействия f(t) на входе модулятора в момент формирования n-го импульса (длительностью Tn) соответственно.

Интегралы Aqi, Cqi уравнения (4.16) были определены ранее, поэтому для решения поставленной задачи требуется лишь определить интегралы вида (4.17), (4.18). При этом следует отметить, что, как и в случае систем управления с ШИМ (см. формулы (4.3)–(4.5)), вид данных интегралов не зависит от вида сигнала на входе частотноимпульсного модулятора, поскольку последний преобразует сигнал любой формы в последовательность прямоугольных импульсов равной амплитуды. Поэтому рассмотрим только лишь вычисление соотношения вида (4.17).

С учетом соотношения (4.19) интеграл (4.17) принимает вид

 

 

 

sign(x0 (T

))

 

 

 

 

 

−ρ t

A

=

 

A

1(t T

)

(t T

τ )

 

Di

 

e q dt,

qi

 

 

n=0

n

 

n

 

n

 

 

 

 

 

0

 

 

 

 

 

 

 

 

(4.21)

 

 

 

 

 

i = 0,1,…,n ;

 

 

 

 

 

 

 

 

 

 

q = 1,2,…,m

 

 

 

Рассуждая аналогично вычислению соответствующего интеграла для САУ с ШИМ, получаем следующее

Aqi = Di

0

 

 

A

 

 

n=0

 

sign (x

 

(Tn )) 1(t Tn )1(t Tn

 

 

 

−ρ t

 

 

0

τ )

 

e

,

 

 

q

dt = Aqρiq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i = 0,1,…,n ;

q =1,2 …, m,

(4.22)

где

 

 

 

Aq =

A(1 e−ρqτ )

sign (x0 (Tn ))e−ρqTn .

(4.23)

ρq

 

n=0

 

 

129

Для вычисления соотношения (4.23), а, следовательно, и интеграла вида (4.22), необходимо рассмотреть вопрос определения периода следования импульсов Tn на выходе модулятора, как функции входного сигнала ЧИМ. Решение данного вопроса непосредственно связано с конкретными схемными реализациями ЧИМ.

4.3. Частотно-импульсные модуляторы 1-го рода

При импульсной модуляции 1-го рода, как показано в [9], значения модулируемого параметра определяются мгновенными значениями модулирующей функции в моменты срабатывания импульсного элемента, причем между обеими величинами существует функциональная зависимость вида an = f (xn ), где an = a(tn ) – модулируемый параметр; xn = x(tn ) – модулирующая функция; tn (n = 0, 1, 2,...) – моменты срабатывания импульсного элемента.

Рассмотрим функциональные и принципиальные электрические схемы некоторых двухполярных частотно-импульсных модуляторов [9].

Двухполярный ЧИМ 1-го рода на базе аналоговых микросхем

Функциональная схема рассматриваемого модулятора и диаграммы сигналов, поясняющих работу ЧИМ, показаны на рис. 4.4, где приняты следующие обозначения: 1 – блок вычисления абсолютной величины входного сигнала; 2 – запоминающее устройство, реализованное в виде интегратора с коэффициентом интегрирования K1; 3 – времязадающий блок, в состав которого входят интегратор с коэффициентом интегрирования K2 и релейный элемент с зоной нечувствительности 1; 4 – блок выдержки времени, в состав которого входят интегратор с коэффициентом интегрирования K3 и релейный элемент с зоной нечувствительности 2; 5 – блок формирования импульсов, реализованный на релейном элементе с зоной нечувствительности .

Модулятор работает следующим образом. В исходном состоянии сигналы на выходе всех трех интеграторов равны нулю. После подачи на вход модулятора воздействия x(t), величина которого превышает значение зоны нечувствительности релейного элемента, замыкается ключ SA5 и интеграторы (блоки 2 и 4) начинают интегрировать до тех пор, пока сигнал y3(t) не достигнет значения 2, которое определяется следующим образом:

2= E0K3τ.

130