Файл: Основы автоматического управления..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 304

Скачиваний: 15

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

682

ГЛ . 17. САМ ОНАСТРАИВАЮ Щ ИЕСЯ СИСТЕМ Ы

 

Аки . . .,

Акп. В результате получим

 

 

I = I (fcj -f- Ski, . . . , kn-p Akn) Ä# I (ky, . .. ,

kn) -f

 

 

П

 

 

 

+ s

4 k ai sin

(17-5-5)

 

г=і

 

Умножая это равенство на at sin о)г< и осредняя результат по вре­

мени,

 

получим

 

 

 

 

 

 

 

 

-jf

т

Іщ sin coit dt та I

(klt

 

kn) at

г

sin сo;< dt -f

 

I

. . . ,

j

 

 

о

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

n

 

 

T

 

 

 

 

 

 

 

+ 2

~§k aiQ-i--j- j

sin соit sin (Ojt dt. (17.5.6)

Ho

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin coddt =

1 — cos wiT

 

 

 

 

 

T - J si

 

 

(£>;T

 

I

T

 

. .

. ,,

1

Г sin (cof— tOi) T

sin (coi+ coj) T 1

f .

T

J

sin cousin CO| f ü = -r

[

 

 

------------- COt+COj

J

 

0

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

. j.

 

sin 2 согГ

 

 

 

 

 

1

f .

1

 

 

 

 

 

- у -

J S in 2 сог* dt

 

4<Bj7’

 

Если взять время осреднения Т кратным всем периодам 2я/соі, . . .

. . ., 2я/соп, то первые два интеграла будут равны нулю, а тре­ тий — половине. При любом достаточно большом времени осред­ нения Т первые два интеграла будут как угодно малы, а третий как угодно близок к 1/2. Следовательно, формула (17.5.6) может

быть переписана

в виде

 

 

л

2

 

 

-у- j Ійі sinCöi* dt та- у - -щ-

(i= 1, 2, . . . , п).

(17.5.7)

о

 

 

 

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


§ 17.5. ЭК СТРЕМ А Л ЬН А Я НА СТРО Й КА ПА РА М ЕТРО В

683

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

Выѵиеяитеяб кри­ териятѵестеа

/

а, sin со, і

 

a„sinco„t

7

1

1

7+Ts

71-Ts

71-Ts

an svnco„t

 

 

 

 

*

 

У

 

 

/7осяеВоеатеяекая

ОВъекяг

 

 

- корректирующая

уяравяекия

 

 

 

 

цепь

 

 

 

 

 

 

 

 

 

 

корректирующая

 

 

 

 

 

цепе

 

 

 

 

 

 

оВраткой связи

 

 

 

 

 

 

Рис. 17.5.2.

 

 

Для

непрерывного

осреднения по

времени произведений

Ia t sin

соit можно использовать апериодические звенья с большими

постоянными времени Т, например цепочки RC. В этом случае

интеграл в (17.5.7) заменяется интегралом в пределах

от t0 до t

от произведения Ia t sin

соtt на весовую функцию апериодического

звена:

 

і

_ 1-1

 

 

 

 

 

 

 

 

=

dki

a-LT Jf е

т lat sin cojt dx

(i = l, . . . , » ) .

(17.5.8)

 

 

to

 

 

 

 

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

представлена на рис. 17.5.2.

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


684 ГЛ . 17. САМ ОНАСТРАИВАЮ Щ ИЕСЯ СИСТЕМ Ы

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

Рассмотренные методы определения градиента критерия каче­ ства управления используют очень небольшой объем априорной информации об объекте управления и вследствие этого требуют значительного времени на поиск экстремума. Время, необходимое для определения градиента критерия качества, а следовательно и на поиск его экстремума, можно существенно сократить, если использовать значительную априорную информацию об объекте управления, во многих случаях имеющуюся в распоряжении конструктора. Так, например, если известны дифференциальные уравнения объекта, можно легко составить другие, вспомогатель­ ные уравнения, определяющие частные производные выходных переменных системы по параметрам системы управления &і, . . ., кп, решая которые, можно вычислить и производные критерия качества управления по параметрам кі, . . ., кп. Вводя в систему быстро­ действующую вычислительную машину, например используя моделирование в сокращенном масштабе времени, можно суще­ ственно ускорить вычисление градиента критерия качества. Для применения этого метода необходимо, чтобы были известны не только уравнения объекта управления, но и все входящие в эти уравнения параметры, которые заранее не могут быть известными. Поэтому, используя этот метод ускоренного определения градиента критерия качества управления, необходимо ввести в систему управления дополнительные устройства, определяющие неизвест­ ные параметры объекта управления І і , . . 1т. Для этого можно использовать в системе управления модель объекта управления с автоматической настройкой ее параметров с таким расчетом, чтобы эти параметры были возможно более близкими к параметрам объекта управления [44, 22] (так называемая подслежиеающая модель). При этом для наилучшего приближения параметров модели к параметрам объекта управления можно использовать экстремальную настройку, обеспечивающую экстремум некоторого критерия Іі близости поведения модели к поведению объекта. Общая структурная схема такой системы приведена на рис. 17.5.3.

Для поиска экстремума критерия качества I можно также использовать различные другие методы. В частности, можно использовать статистические методы нахождения экстремума или градиента I по случайному набору значений параметров к і , . . ., к п .

В качестве критерия качества управления в системах с экстре­ мальной настройкой часто принимается среднее значение квадрата

параметра управления (сигнала ошибки), характеризующего ошиб­

ку управления:

 

і

 

/ = 4 - J А2('c)dT.

(17.5.9)

t - T


§ 17.5. ЭК С ТРЕМ А Л ЬН А Я Н А СТРО Й К А ПАРАМ ЕТРОВ

685

Ввиду того что организовать непрерывное вычисление по форму­ ле (17.5.9) затруднительно, так же как и при высчилении средних значений произведений Ia t sin юг£ при определении градиента методом синхронного детектирования, для осреднения величины

Рис. 17.5.3.

Д2 (т) часто используют апериодическое звено с большой постоян­ ной времени Т. Соответственно формула (17.5.9) заменяется формулой

t

t -

X

 

7 = - | r j

е~ ~

Д£(г)йт.

(17.5.10)

to

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

t

 

Іі = 4" j [г/м (т) — y{x)\*dx.

(17.5.11)

t-т

 

При использовании для осреднения апериодического звена с боль­ шой постоянной времени Т формула (17.5.11) заменяется

686

Г Л . 17. САМ ОНАСТРАИВАЮ Щ ИЕСЯ СИСТЕМ Ы

 

формулой

1 _ 1-т

 

 

 

 

h = -jr je ~ т [ун (т)-у(т)]Чт.

(17.5.12)

 

ІО

 

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

N

 

 

 

Уі = 2 alh{t)yh-\-Fi(t)

(7=1,

N \

(17.5.13)

h=l

 

 

 

где уі = у — выходная переменная системы, у2 — выходная пере­ менная корректирующей цепи в обратной связи, у 3, . . ., y N — вспомогательные переменные, Ft (t) — внешние возмущения, дей­ ствующие на систему, alk (t) — коэффициенты системы, зависящие

от времени и параметров системы klt

. .

кп. Параметр управления

А (t) выражается

формулой

 

 

 

 

 

 

 

 

А (t)

=

X

(t) — у2 (і),

(17.5.14)

где я (<) =

F, (г) — входной сигнал

управления. Подставляя это

выражение

в

формулу

(17.5.9),

получим

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

I = -jr

j

[Уг (t) — x{x)fdx.

(17.5.15)

 

 

 

 

1 -Г

 

 

 

 

 

Продифференцируем это выражение по параметру kt:

-|£ т = Т -

( f

e

W

-

z

M

(

i =

l .........n). (17.5.16)

1

 

t - т

 

 

 

 

 

 

 

 

Для определения производных dy2ldkt продифференцируем по пара­ метру кі уравнения (17.5.13). В результате получим систему уравнений, определяющую частные производные переменных уі

по

параметрам /сг:

П

 

 

 

d

дуі

 

dyk

 

 

 

2 aihw

daih (0

Ук

1 =

1 ,

dt

dkl

dkl

dki

=

,

 

 

h—i

 

h=l

 

 

 

E

1 1

Эта система уравнений совместно с формулой (17.5.16) определяет

частные производные критерия качества по настраиваемым пара­ метрам. При вычислении частных производных дуі/дкі исполь­ зуется априорная информация о системе в виде ее дифференциаль­ ных уравнений (17.5.13). Система уравнений (17.5.17) должна решаться при нулевых начальных условиях.


§ 17.5. ЭК СТРЕМ А Л ЬН А Я Н А СТРО Й К А ПАРАМ ЕТРОВ

687

П р и м е р 17.5.1. Рассмотрим самонастраивающуюся систему стабили­ зации продольной оси самолета или ракеты в вертикальной плоскости с экстремальной настройкой коэффициентов усиления сигналов от свободного и демпфирующего гироскопов. Движение самолета в прямолинейном гори­ зонтальном полете при малых углах атаки а и малых отклонениях вектора скорости от горизонта Ѳ без учета случайных возмущений описывается урав­ нениями (3.16.10):

а + с.аД -саа = сп—сбб,

Ѳ= Лаа —— , д = Ѳ+ а. (17.5.18)

а

V

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

та д. Предполагая для простоты, что эти измерения осуществляются абсо­ лютно точно и что рулевая машина безынерционна, примем уравнение работы рулевой машины в виде

б = fe4d +• fe2d,

(17.5.19)

где fei, fe2 — коэффициенты усиления. При полете с различными скоростями в большом диапазоне высот коэффициенты уравнений (17.5.18), как было отмечено в примере 17.3.1, изменяются в сотни раз. Чтобы система управле­ ния (автопилот) хорошо стабилизировала ось самолета при любых условиях полета, целесообразно применить экстремальную настройку параметров ки fe2. При этом качество управления естественно считать тем лучшим, чем меньше вектор скорости самолета отклоняется от горизонта. В этом случае за крите­ рий качества управления естественно принять среднее значение квадрата

угла отклонения вектора скорости

самолета от

горизонта:

 

г

 

/ = 4

| ѳ2 (0<**.

(17.5.20)

 

о

 

Для определения градиента величины I применим метод вычисления ее про­ изводных по kt и fe2 с помощью уравнений движения самолета (17.5.18). Дифференцируя формулу (17.5.20) по fet и fe2, получим

т

 

’ • ~ щ г = т J e - S " “

(17'5'211

о

 

Для вычисления производных d&dkt продифференцируем по fe4 и к2 урав­ нения (17.5.18) и (17.5.19). В результате получим

d2

да

,

d

да

да

д&

N

Ш ~ д І й ^ Са‘

'1 Г 'дІц +Са''Ш н ~

С&~дкГ ’

 

d

дѲ

Ла

да

ад

 

дѲ л.

да

 

dt

дкі

дкі

дкі

 

дкі

дкі

(17.5.22)

 

д&

»

,

зд

 

d

ад

 

 

+ ^2

 

 

 

dkt

“ ^ +Аі

дкі

dt

дкі

У

 

 

д&

—,$А-кі ^

J_ Ь

d ад

 

 

 

дк2

'

1 dk2

Т"«2

dt

дк2

 

 

Таким образом, в системе управления должны быть вычислительные устрой­ ства, интегрирующие уравнения (17.5.22) и вычисляющие производные кри­ терия качества I по формуле (17.5.21).