Файл: Уткин, В. И. Скользящие режимы и их применения в системах с переменной структурой.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 127
Скачиваний: 3
§ 2] ПРОЦЕДУРА НАСТРОЙКИ ПАРАМЕТРОВ МОДЕЛИ |
259 |
а вероятность попасть на решение, целиком принадлежа щее этому подмножеству меры нуль, при произвольном выборе начальных данных равна нулю. Поэтому у пас 'тк_возникает задачи сходимости процесса восстановлепия 'оператора. Заметим также, что предлагаемый метод прин ципиально пригоден и в том случае, когда система (13.1) неустойчива (например, система (13.2) с негурвицевой матрицей А). Требуется лишь, чтобы скорость процесса позволила реально наблюдать и измерять его на протяже нии некоторого времени, достаточного для выхода в сколь зящий режим системы «объект — модель» и, если нужно, для организации процесса поиска.
§ 2. Процедура настройки параметров модели
Опишем одну из процедур поиска таких параметров
/модели, при которых все величины щ экв окажутся постоян ными, и следовательно, искомые параметры объекта могут быть определены из соотношений (13.9).
Дополним уравнение (13.8) аналогичными уравнения ми, в которых в левых частях вместо фу (х, it) стоят функ
ции фу (х, t) = Lk {фу (х, i)}, т. е. функции, полученные в результате /с-кратного применения какого-либо линей ного оператора L к фу (х, t). В этом случае правые части, вообще говоря, являются функциями времени.
Введем обозначения
N |
|
|
|
|
|
2 |
(*, 0 = |
е? (О |
(к = |
0,. .. , N — 1). |
(13.11) |
i=i |
|
|
|
|
|
(Обратим |
внимание |
на то, |
что |
функции ef не |
равны |
Lk{е?}, так как в процессе поиска параметров модели коэф фициенты l(j будут изменяться во времени и тогда левая 'часть в уравнении (13.11) не является результатом /с-кратиого применения оператора L к линейной комбина-
N
ции 2 k№h которая согласно (13.7) равна е°.)
3=1
Убедимся теперь, что все функции sf могут быть изме рены непосредственно, несмотря на то, что параметры объекта а гу, от которых зависят коэффициенты li}- в (13.11),
260 ИДЕНТИФИКАЦИЯ ДИНАМИЧЕСКИХ ОБЪЕКТОВ [ГЛ. XIII
неизвестны. Из условий st = 0 при щ = г4ВКв следует
N N N
2 |
аиЪ (*> 0 = 2 |
0 + ( 2 da I фi (*> 0 1) “ ibkb- (13.12)^ |
|
3=1 |
3=1 |
'3=1 |
' |
Применяя линейный оператор Ькк правойи левой части уравнения (13.12), получаем
N N N
2 а«Фз- (я,о = £к'{2 ьрФз(ж. о +(2 dpiФз(*>oiWbkd}.
3=1 |
|
3=1 |
|
'з = 1 |
' |
|
|
|
|
|
(13.13) |
Соотношения |
(13.13) |
позволяют переписать |
уравнения |
||
(13.11) |
следующим образом: |
|
|
||
N |
|
N |
|
|
|
^ |2 h |
m {.Х : |
0 + (2 |
I Ф; { х , 01) М1экв} |
|
|
4 = 1 |
|
3=1 |
' |
1 |
|
|
|
|
N |
|
|
|
|
|
— 2 |
ь«зФ- ix >0 = |
Si • (13.14) |
|
|
|
3=1 |
|
|
Как следует из (13.14), все функции в левой части известны
и это позволяет вычислить величины ef.
Линейный оператор L подберем таким образом, чтобы
для функций фу |
(х, it) выполнялось условие |
|
det ПФ/(зг, 0 1|¥= |
0 (/ = 1 ,..., iV; к = 0 , . . . , N - |
1). (13.15) |
Рассмотрим поведение вспомогательных |
функций |
N
1V)
ц= — 2 1Ф если коэффициенты Ьц модели изменяются
3=1
по следующему закону:
N -1 |
|
Ъхj —- с 2 е£фj(x,t) (с — const, с)> 0). |
(13.16) |
к=0 |
|
Запишем производные по времени для функций гзр.
N
V i = 2
3=1
§ 3] ПРИМЕРЫ ВЫБОРА ЛИНЕЙНЫХ ОПЕРАТОРОВ 261
или с улетом (13.16)
N N -1
v i — |
~ с |
Zj; |
2 eicP; |
(^1 Z). |
(13.17) |
' — |
j = l |
k = 0 |
|
|
|
Уравнение (13.17) |
согласно |
(13.11) |
приводится |
к виду |
JV-1
Щ= — С 2 (8i)2I
к=о
т. е. функции vt являются ограниченными снизу (^ > 0) и монотонно убывающими при всех е£, одновременно не
равных нулю. Процесс прекратится, когда все ef обра тятся в нуль. Из уравнений (13.11) и условия (13.15) *) следует, что при этом
lim lij — 0,
( —►СХЭ
ИЛИ
lim Ъц = &'хj, /—►00
Таким образом, процедура поиска (13.16) для линейных операторов типа (13.15) позволяет определить все пара метры идентифицируемого объекта.
§ 3. Примеры выбора линейных операторов
Описанная процедура определения параметров объекта правомерна для линейных операторов типа (13.15). Приве дем примеры таких операторов.
П р и м е р 1. Пусть объект является линейной системой (13.2), т. е. ер,- = Xj, И является квадратной матрицей раз мерности п X п с действительными и различными собствен ными числами. Покажем, что для линейного оператора L типа «элемент запаздывания» условие (13.15) выполняется.
*) Предполагается, что процесс настройки параметров моде ли происходит намного быстрее процессов в объекте и если
lim del |tpf |= 0, то Н т (е^ / clet |tpf ||) = 0.
/ —►оо |
/ —►СО |
262 ИДЕНТИФИКАЦИЯ ДИНАМИЧЕСКИХ ОБЪЕКТОВ ЕГЛ. XIII
В рассматриваемом случае вектор х |
имеет вид |
|
||||||||
|
|
|
|
х = ПФ (t) с, |
|
|
|
(3.18) |
||
где D — матрица, столбцами которой являются собствен |
||||||||||
ные векторы матрицы А, |
Ф (t) — диагональная матрица с |
|||||||||
элементами |
е“*\ оц (£ = 1, |
. . ., |
п) |
— собственные |
числа |
|||||
матрицы |
А, |
с — вектор-столбец с компонентами с1; .. . , |
||||||||
сп, |
зависящий |
от начальных условий. |
Для выбранного |
|||||||
оператора |
L с |
запаздыванием т задача |
сводится |
к до |
||||||
казательству следующего условия: |
|
|
|
|||||||
|
det |xKj |
(t) I = det |Xj (t — (к — 1) т) |=j= О |
|
|||||||
|
|
|
|
(/, |
к = |
1, . . ., п). |
|
(3.19) |
||
Как |
следует |
из (3.18), |
этот |
определитель имеет вид |
||||||
|
|
|
|
|
|
t |
П |
|
|
|
|
|
|
|
|
|
Z a |
|
|
||
|
det ||хкI — схсч. . . спdet |D |е |
i= l |
det II e |
|
По сделанному ранее предположению, компоненты век тора х в (3.18) должны быть линейно независимыми, поэтому все с, отличны от нуля; матрица D состоит из собственных векторов матрицы А с различными собствен ными числами, поэтому det||Х>| 0; наконец, det |е_аА;~1)т| является определителем Вандермонда и поэтому также отличен от нуля. Таким образом, неравенство (3.19) до
казано. |
|
П р и м е р 2. Рассмотрим теперь применение к функ |
|
циям |
X j ( t ) линейного оператора L типа «инерционное |
звено» |
с постоянной времени Т и коэффициентом усиле |
ния, равным единице, для линейной системы, которая
изучалась |
в примере 1. |
Предположим, |
что |
собственные |
||
движения |
инерционных звеньев затухают |
столь быстро, |
||||
что ими можно пренебречь. |
|
|
|
|||
|
В результате ^-кратного применения оператора L к |
|||||
вектору х (3.18) получим: |
|
|
|
|||
|
|
хк = |
D Ф?: (if) с, |
|
|
|
где |
Фк (г) — диагональная |
матрица |
с |
элементами |
||
еаз1/{Ta,j + |
1),:. Аналогичный |
(3.19) определитель для рас |
||||
сматриваемого оператора, составленный |
из |
компонентов |
§ 3] |
ПРИМЕРЫ ВЫБОРА ЛИНЕЙНЫХ ОПЕРАТОРОВ |
263 |
||||
векторов |
хк |
(к = |
0, . . |
п — 1), запишется следующим |
||
образом: |
|
|
|
П |
|
|
|
|
|
|
,, |
|
|
|
|
|
|
1.s а |
|
|
det II хЧII = с,С9. .. сп det IID II е г~1 |
d e t ------------г |
• |
||||
II |
; l l |
1 - |
п |
1 1 |
I {Tctj + l ) fc |
При изучении оператора типа «элемент запаздывания» в примере 1 было показано, что все с,- и det | | не равны нулю. Последний сомножитель, как и в примере 1, являет ся определителем Вандермонда и поэтому ои также отли чен от нуля. Это означает, что и для оператора L типа «инерционное звено» выполняется условие (13.15), при котором процедура поиска (13.16) позволяет найти пара метры объекта.
П р и м е р 3. Покажем, что для систем произвольного вида (13.1) в качестве оператора L можно выбрать опера тор дифференцирования
кdk<Pj (*. О
фу “ |
dtk |
' |
Докажем, что найдутся такие интервалы времени, на которых условие (13.15) выполняется. Доказательство проведем методом математической индукции.
Рассмотрим определитель второго порядка
|
А 2(t) |
Ф 1 |
Ф г |
|
|
|
|
ф1 |
фг |
|
|
||
|
|
|
|
|
||
Допустим, |
что А 2 |
(t) |
= 0. |
Тогда |
в каждый момент |
|
его столбцы |
линейно |
зависимы, |
т. |
е. |
||
|
< h ( t ) ф1 + сь(£) ср2 = 0, |
(13.20) |
||||
|
ai (*0 ф1 + |
а2.(<0 Фг = |
0. |
|||
|
|
Функция а2(£) не равна тождественно нулю, поэтому
найдется интервал времени, |
на котором а2(t) не обра |
|
щается в нуль и система (13.20) представима в виде |
||
ai (t) Ф1+ |
Фг = |
0, |
«1 (£) ф ! + |
фг = |
(13.21) |
0. |
||
После дифференцирования |
первого уравнения в (13.21) |
|
и вычитания из него второго |
получим |
Я1 Ф1 = 0. |
( 13. 22) |