Файл: Уткин, В. И. Скользящие режимы и их применения в системах с переменной структурой.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 если коэффициенты Ьц модели изменяются

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)