ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 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 всеми изложенными способами можно также применить
§ 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).