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

книги из ГПНТБ / Цейтлин Я.М. Проектирование оптимальных линейных систем

.pdf
Скачиваний:
2
Добавлен:
24.10.2023
Размер:
9.35 Mб
Скачать

Поэтому

It

f W = 4 r § F

или

*+ajf!!? (

П р и м е р 2. Пусть

F(s)

(s ) e S t d s = S

Res., [ F (s)].

ft=i

 

(1.42)

 

 

ft;-l

 

1

 

d kі— [ ( s - a / ^ s ) ] .

1

i }

ds'

 

=

 

1

 

(a+*)(b-s)

Здесь F (s) —функция, совпадающая с F (s) примера 1,—уже не преобразование Фурье, а преобразование Лапласа.

По формуле (1.42)

/ (О = Res . a

(o +

s)

( 6 - s )

+

Res6 L (a + s)(ft-s) J'

где согласно (1.38)

 

 

 

 

 

 

 

Res _а [•••]

=

I

\b

— s

 

a+

b

 

 

0!

 

R e s 6 [ . . - ]

 

1

 

 

 

 

 

 

0!

\ a

+ s / s = i

 

 

 

 

 

 

 

 

Поэтому

 

 

 

 

 

 

 

 

/(0

=

1

( e w + e - e ' ) .

^ 0 .

 

a + 6

 

 

 

График этой функции изображен на рис. 6. Видно, что эта функция значительно отличается от функции / (і) в примере 1.

Полученные результаты с учетом свойства единственности преобразований Фурье и Лапласа свидетельствуют о том, что обрат­ ные преобразования Фурье и Лапласа одной и той же функции комплексного переменного F (s) различны, если F (s) имеет полюсы как в левой, так и в правой полуплоскостях, и одинаковы (иден­ тичны), если F (s) имеет полюсы только в левой полуплоскости.

7. ОСНОВНЫЕ ПОНЯТИЯ ТЕОРИИ ЛИНЕЙНЫХ СИСТЕМ

Каждая система автоматического управления является уст­ ройством, осуществляющим преобразование информации на входе системы в информацию на ее выходе.

Информация на входе системы представляется сигналом той или иной физической природы, пропорциональным функции вре-

Мени х (t), информация на выходе системы — также некоторым физическим сигналом, пропорциональным функции времени у (t). Функции х (t) и у (t) могут быть как скалярами, так и векторами. Эти функции называют иногда сокращенно входом и выходом системы.

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

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

а)

хШ

F

Рис. 6 Рис. 7

На рис. 7, б схематично изображена многомерная система, имею­ щая N входов и М выходов.

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

Зависимость между функцией х (t) на входе системы и функ­ цией у (t) на ее выходе осуществляется оператором системы F, который является наиболее общей математической характеристи­ кой системы,

Известно, что оператором называется совокупность математи­ ческих (а в общем случае и логических) операций, в результате выполнения которых каждой функции заданного множества при­

водится в соответствие некоторая функция другого

множества.

Тот факт, что функция у {t) является результатом

преобразо­

вания функции х (t) системой, характеризуемой оператором F, записывается в виде операторного равенства

У (0 = Fx {t)

или в более общем случае

у (s) = Fx (t).

Примером операторов этого типа является оператор преобра­ зования Лапласа, приводящий в соответствие множеству ориги­ налов, т. е. функций времени х (t), преобразуемых по Лапласу, множество функций комплексного переменного X (s), представ­ ляющих собой изображения функций х (t). Это преобразование записывается в виде

 

 

 

 

 

со

 

 

 

 

 

X(s)

= L[x(t)]

= j

 

x(t)e-stdt.

 

Приведем примеры

операторов.

 

 

 

 

1.

Оператор

возведения

в квадрат

 

 

 

 

 

 

у = Ах =

х2.

 

 

 

2.

Оператор умножения функции х (t)

на заданную

функцию

 

 

у (t) = Fx (t) =

ф (t) х (О-

 

3.

Оператор

дифференцирования

 

 

 

 

 

 

 

y(t) = Fx(t) =

^

L .

 

В теории автоматического управления оператор дифференци­

рования принято обозначать буквой р.

 

 

Оператору рк соответствует оператор

определения

производ­

ной

k-то порядка. В

частности,

 

 

 

 

 

 

 

y(t)=*px(t)

=

-$rx(t)

=

*(t);

 

 

у (t) = рН (t) = р [рх (01 = -~x(t)

= х" (*);

 

 

y(t) = pkx (t) = р [рк - i x

(*)] = -JL-

х (t) = х" (t).

 

4.

Оператор

интегрирования

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

y(t) =

Fx(t)

=

 

\x(l)dt.

 

 

 

 

 

 

о

 

 

 

Этот оператор, обратный оператору дифференцирования р, принято обозначать ~ . Естественно, что оператор —^- опреде-

ляет

^-кратное повторное

интегрирование. В частности,

 

 

 

 

 

 

 

н

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

/

t

 

 

 

 

 

 

y(t)

=

-jrx(t)

=

\\x(mdr\;

 

 

 

 

 

 

 

 

 

 

 

о

о

 

 

 

 

y{t)

=

\ x

(

t )

= \ . . .

Ix&Jd^di,-

 

-dlk.

 

 

 

 

Р

 

 

0

 

0

 

 

 

 

5.

Оператор

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(p)

=

bJL±hE

 

 

 

где p — оператор

дифференцирования.

 

 

 

Этот оператор ставит в соответствие функции х

(t)

функцию у (t),

определяемую решением дифференциального

уравнения

 

(а„ +

агр

+

а2 р2 ) у

(t)

=

(b0 +

btf)

х

(t),

которое с учетом свойств оператора

дифференцирования может

быть

записано

в

виде

 

 

 

 

 

 

 

 

 

а2у"

(t)

+

аіУ'

(t) + а0у

(t)

= blX'(t)

+

b0x {t).

6.

Нелинейный

интегральный

оператор

 

 

 

 

 

 

y(s)

=

Fx(t)==

J <p[*(/), t,s]dt,

 

 

T

где ф [. . .] — заданная функция, нелинейная относительно пере­ менной х.

7. Оператор решения нелинейного дифференциального урав­

нения у (t) .==

Fx

(t)

 

 

 

 

 

у" (t)

+

a (t)

sin

[b

(і) y(i)]=x

со­

существуют и значительно более сложные операторы.

Оператор F называется линейным, если при произвольных

элементах х и

х2,

. . ., хп,

принадлежащих

множеству М, и

любых постоянных

kx,

k2,

. . .,

kn

справедливо соотношение

 

 

F

( І

ktxt^

=

Д k,Fxt.

(1.43)

Примерами линейных операторов являются операторы, ука­ занные в пунктах 2, 3, 4, 5. Операторы, приведенные в пунктах 1, 6, 7, являются нелинейными, так как для них не удовлетворяется соотношение (1.43), в чем нетрудно убедиться.

З Я. М. Цейтлин

33

8. ОПЕРАТОРЫ ЛИНЕЙНЫХ ДИНАМИЧЕСКИХ

СИСТЕМ. ФУНКЦИЯ ВЕСА

Системы автоматического управления, часто называемые ди­ намическими, бывают линейными и нелинейными. Наиболее ис­ следованным классом динамических систем является класс линей­ ных систем.

Динамическая система называется линейной, если она обладает свойством суперпозиции сигналов (удовлетворяет принципу су­ перпозиции). Этот принцип (свойство) предусматривает соблю­ дение соотношения (1.43), из которого непосредственно следует, что при изменении масштаба входного сигнала в k раз без измене­ ния его формы масштаб выходного сигнала также изменяется в k раз (форма его остается неизменной) и что при подаче на вход суммы нескольких воздействий сигнал на выходе равен аналогичной сумме выходных сигналов, обусловленных каждым из воздействий на входе. Таким образом, оператор линейной динамической системы является линейным оператором.

Рассмотрим динамическую систему, работа которой описывается линейным дифференциальным уравнением с переменными (зави­

сящими от времени)

коэффициентами

 

а„ (0 yw (0 + a«-i (0 У(п~Х) ( * ) + • • • + % (О У (0 *f-

 

+ fl0 (Г) у {t) =

ьт

(о *<"*> (о +

ът_х х-1) (ОН

V

 

+

(t) +

b0(t)x(t).

 

Если воспользоваться оператором дифференцирования р и ввести в рассмотрение операторы А (р, t) и В (р, t), определяемые равенствами

п

А(р,

*) =

M ' ) +

fli(')p+--

- + а « ( ' ) р п

= £

«ЛОР';

 

 

 

 

 

 

 

i*=0

 

В(р,

t) =

b0(f) + b1(t)p-\

 

\- bm (t)pm

т

bt (t)p\

 

= S

 

 

 

 

 

 

 

=o

дифференциальное уравнение

примет вид

 

 

 

 

 

Л (р,

t)y(t)

=

В (р,

t)x(t),

 

откуда следует

формальное равенство

 

 

 

 

 

^ > = 4 о Н Н ' ) = / : > '

о*(о,

 

где

 

 

 

 

 

 

 

 

является

оператором системы.

 

 

 

 

 

34

В том случае, если коэффициенты дифференциального уравне­ ния постоянны:

ak (t) = ak, bk

(t)

=

bk,

выражение для оператора системы F

(р,

t)

принимает вид

Такой оператор, не зависящий от времени (переменной, по которой производится дифференцирование), называется стационарным.

Оператор F (р, і), который явно зависит от времени (пере­ менной, по которой производится дифференцирование), назы­ вается нестационарным. .

Операторы F (р, t) или F (р) характеризуют линейную систему. Однако обратное соображение, вообще говоря, не имеет места, так как существуют линейные системы, т. е. такие системы, для которых справедлив принцип суперпозиции, которые не описы­ ваются линейными дифференциальными уравнениями. Напри­ мер, система, осуществляющая операцию чистого запаздывания,

описываемая уравнением

-

 

 

 

 

у

(t)

=

х (t у),

 

линейна.

 

 

 

 

 

Оператором этой системы, как

известно [7,

11, 13], является

 

F

(р)

=

e-w.

 

Этот оператор не относится к классу операторов,

соответствующих

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

Поскольку рассматриваемая система, характеризуемая опера­ тором F (р, t), является линейной, то к ней применим принцип суперпозиции (1.43).

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

Импульсной функцией называется четная функция б (t — т), обладающая следующими свойствами:

( 0, t + x

(1.44)

J 6(* —т)гіт = -1,

где ЄІ и є2 — любые сколь угодно малые положительные числа. Иными словами, это такая четная функция, которая равна нулю всюду, кроме одной точки т = t} где она обращается в бес-

3*

35

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

Пользуясь

б-функцией,

можно

для

любой непрерывной

в точке t функции /

(т)

написать

 

 

 

t +

E2

 

 

 

 

 

 

J f(x)6(x

— t)dx=f(t).

(1.45)

Действительно,

б ( т — t)

равно нулю

всюду,

кроме точки X = t.

В бесконечно малой окрестности этой точки функция / (т) при­

близительно постоянна и равна f (t). Вынося ее за знак

интеграла

и принимая во внимание (1.44), получим формулу (1.45).

Важной для практических приложений является также фор­

мула

 

 

 

 

\f{x)bw{x-t)dx

= {-\ffw

(1.46)

а

 

 

 

 

или, так как

 

 

 

 

6<*>(т-0 =

( - 1 ) * 6 ( < - т ) , .

 

ь

 

 

 

 

lf(x)bw(t-x)dx^f^(t).

 

 

 

а

 

 

 

 

Формула (1.45) может быть записана в виде

 

ь

 

 

 

 

\f'(x)b(x-t)dx

=

f'(t),

(1.47)

поскольку любая производная

функции

/ (f)

является

некоторой

функцией, а формула (1.45) справедлива для любой функции,

непрерывной в точке t-

 

 

 

 

 

Проинтегрировав

по частям левую часть (1.47),

получим

 

 

 

б

 

 

 

 

 

/ (т) б(т -

0| -

J f (т)б' (x-t)

dx = Г (О-

 

 

 

 

 

а

 

 

 

 

Первое слагаемое левой части

равно нулю, поскольку

б-функ-

ция в силу

свойства

(1.44)

равна

нулю

в точках а

и 6.

Поэтому

 

ь

 

 

 

 

 

 

 

lf{x)&'(x-f)dx

=

-f'(t);

 

 

 

и

 

 

 

 

 

(1.48)

 

 

 

 

 

 

 

 

lf(x)8'(t-x)dx

=

f'(t).

 

 

Из (1.48)

следует

 

 

 

 

 

 

\f'(x)6'{x-t)dx=;-r(t)-

Проинтегрировав

по частям

левую часть,

получим

 

'

ь

ь

 

 

 

 

 

f (т) б' (т -

0 | -

J f (т) 6" (т - 1 ) dx =

-

f" (t),

 

откуда

 

 

 

 

 

 

 

j f (т) 6" (x — t)dx

=

j

f (T) 6" (* — T ) Л

=

/" (/),

(1.49)

поскольку

 

 

 

 

 

 

 

 

б(*) (т —

t)

= 0, т Ф t.

 

 

 

Выражение (1.49) соответствует формуле (1.46) при k = 2.

At 2At

зat

(n-vat nut

 

Рис.

8

Аналогичным образом можно показать справедливость фор­

мулы (1.46) при любом

k.

 

Вернемся к принципу суперпозиции. Этот принцип при задан­

ном операторе линейной

динамической системы позволяет свести

исследование реакции системы у

(t) на произвольное воздействие

х {t) к исследованию реакции системы на типовое воздействие, если воздействие на входе представляется линейной комбинацией (суперпозицией) типовых воздействий. В качестве типового воз­ действия целесообразно принять б-функцию б (t — т).

Представим функцию х (t), воздействующую на систему, как последовательность элементарных прямоугольных импульсов, при­ кладываемых к системе через промежутки времени At (рис. 8). Очевидно, что

я

Это выражение следует непосредственно из формулы (1.45), ко­ торую можно представить в виде

t п

* ( *) = J * ( T ) 6 ( * —T)dT = lim J x(kAt)8(t

— kAt)

At.

о

2 £ ° 0 * = о

 

 

 

Элементарные импульсы х (k

At) At 8 (t —

k At)

при приложении

их к системе вызывают элементарные реакции на выходе,

которые

в силу линейности системы

складываются

(накладываются друг

на друга), образуя суммарную результирующую реакцию на выходе.

Суммарная реакция системы к моменту времени t приближенно равна сумме элементарных реакций, возникших ранее рассмо­ тренного момента времени.

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

мент времени tit

дает реакцию на выходе, отличную от нуля для

t ^5 t{ и равную

нулю для t <^ th поскольку реакция (следствие)

не может появиться раньше, чем обусловивший ее импульс (при­ чина).

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

импульс, соответствующий моменту времени 4 =

k At, представ­

ляется в виде

 

 

х (4)

Ш (t — tk) = х (k At) At6 (t — k

At),

где x (k At) At

площадь данного элементарного прямоугольного

импульса.

 

 

Обозначим через Р (t, | ) реакцию на выходе системы в момент времени t, обусловленную воздействием на входе, представляющим

собой

б (t — ! ) ,

т. е. б-функцию,

приложенную в

момент вре­

мени

Иными словами, функция Р

(t, £) является решением диф­

ференциального

уравнения, если х

(t) — б (t— £).

Эта функция

называется весовой функцией или функцией веса системы и, как следует из приведенного определения, показывает, с каким весом в момент времени t сказывается в реакции системы на выходе воз­ действие вида б-функцйи на входе, приложенное в момент вре­ мени \ . Таким образом, результирующая суммарная реакция на

выходе

системы

будет

 

 

 

 

 

 

 

 

у it)

^ х

(О АО AtP

(t,

0 At)

+

х (1 АО A* Р

(t, 1 АО

+

 

 

 

+

х

(2А0 А* Р

(t,

2At) +. . .

 

 

 

 

 

 

 

1- х (п At)

At

Р (t, n

At).

 

(1.50)

Обозначим

k

At

=

т

и

будем считать,

что по

мере уменьше­

ния At увеличивается

k таким образом, что их произведение

равно

величине т, которая принимает значения от 0 до t. Если At

стре­

мится к

нулю,

то

реакция

на выходе будет представлять

собой.

Предел суммы (1.50), которая с учетом принятых обозначений примет вид

п t

y(t) = \\m J ] x(kM)AtP{t,

kAt) = J x'{x)P(t,

x)dx. (1.51)

M-+0k=0

о

 

rt->-00

Выражение (1.51), играющее исключительно важную роль в тео­ рии автоматического управления и особенно в теории оптимальных систем, является интегральной формой принципа суперпозиции,

если в качестве элементарных

воздействий принята

(выбрана)

б-функция.

 

 

Таким образом, функция веса Р (t, т) является

важнейшей

динамической характеристикой

системы.

 

В соответствии со сказанным ранее для любой реальной системы

Р (t, т) = 0,

t < т .

(1.52)

Выражение (1.52), отражающее причинно-следственную связь, часто называют условием физической осуществимости системы.

Таким образом, функция веса реальной системы, равная, нулю при t < т , может быть записана в виде

Р (t, т) =g(t,

т) 1 (* — т ) ,

(1.53)

где 1 (t — т) — единичная

функция, определенная

известным

равенством

 

 

 

 

( 1,

* S S S T

 

' « - ' Н о ,

« < , .

 

Поэтому для физически реализуемой системы функция у (t), опре­ деляющая реакцию системы на выходе, должна быть записана в виде

t

 

y{t) =

\x(x)g(t,x)l(t-x)dx.

о

 

Заметим, что в подавляющем

большинстве работ (в том числе

во всех учебниках и руководствах) используется формула без единичной функции 1 (t — т), что допустимо, если функция X (t) является обычной (не обобщенной) функцией. В этом случае до­

бавление множителя

1 (t — т) не изменяет результата, определяе­

мого пределами

интегрирования.

 

Если же функция х (t) содержит производные от б-функции,

пренебрежение

множителем

1 — т) приводит к неправильным

результатам.

 

 

 

 

 

х (t) = б' (t).

П р и м е р .

Пусть

функция

на входе

Без учета единичной

функции,

т. е. если

Р (t, х) = g (t, т), получим

 

 

ч

 

 

t

 

y { t )

=

\x(x)g{t,

x)dx=^g(t,

х) 6' (т)dx = g' (t, 0).

Соседние файлы в папке книги из ГПНТБ