Файл: Живоглядов, В. П. Адаптация в автоматизированных системах управления технологическими процессами.pdf

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

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

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

Добавлен: 23.10.2024

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

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

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

со

+ slnu>* J y(x,t)gz(x)dx\ - T (t)mf(t),

о

где f(t)

[1 + y (s in 2 co/+cos2co/)].

(2. 78)

Сравним алгоритмы (2 . 78) и (2 . 42). Различные по фор­ ме, при выборе у (/) =/С (t) они совпадают.

Применение настраиваемых моделей

Пусть параллельно объекту подключена модель с фикси­ рованной структурой и настраиваемыми параметрами Ш(, i = 1 ,

2 , п.

Примем в качестве критерия близости объекта и модели математическое ожидание некоторой функции потерь

W = W (y,

qM), где y = q-\-h, h — помеха, q M— выход модели.

Применяя одновременно аппарат

теории чувствительности и

стохастической аппроксимации, получаем следующий

алго­

ритм настройки коэффициентов т ,

модели:

 

dm.(

dW

i= 0 ,1 ,...,

(2. 79)

~ d f

T (*)V „ W -----т (0

дЯм

 

тде 5 /= ддк — функции чувствительности. Помеху h считаем dm i

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

точной функцией ФГ1 (р), учитывающей специфику распреде­ ленной системы

Sl{p) = <Pir(p)qlJ. (р).

(2 . 80)

Рассмотрим для конкретности, какой вид имеет модель 'чувствительности для объекта с передаточной функциейV

V Ш \ '

<2'8"

тде В(р) и D(p)— многочлены.

440


Дифференцируя Ф(о) до параметрам mt, находим

=Ф(Р)Ф1Г(Р)-

(2. 82')

В табл. 1 приведены выражения Ф/ (р) для некоторых типичных случаев. В последней графе содержится ссылка на лите* ратуру, где исследуются системы с функцией Ф(р) соответ­ ствующего вида. Перечень передаточных функций, конечно* может быть значительно расширен. Однако из приведенной таблицы видно, что лишь в некоторых случаях для получе­ ния функций чувствительности можно использовать звенья с рациональными передаточными функциями.

 

 

 

 

 

Таблица I

m i

Ф(р)

 

& г ( Р )

Литера­

 

тура

 

 

 

 

 

X

— Р~

 

- р

 

[В.47,В.54]'

е

 

 

X

— V "F ~

 

 

 

1В. 54]

е

 

 

 

 

 

а

 

 

а

[B.56J

«1

Ч + Р

 

(ai+P)a

е

 

 

 

 

 

 

а1

I& V 1+Р)(а2+Д)

£ .

/

аз+р

[В.54]

2

у

c c i+ p

 

 

Итак, применение методов теории чувствительности и стохаотической аппроксимации может привести к решению по­ ставленной задачи. Но изложенный формальный подход стра­ дает некоторыми недостатками. С одной стороны, построение распределенных моделей и тем более настраиваемых моделей представляет трудную техническую задачу. С другой сторо­ ны, эго не является необходимым. Рассмотрим способ определения неизвестных параметров без использования мо­ дели с распределенными параметрами [2. 32].

141


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

/?(iV/z)=yW{W(y,(7M) I I1."*),

(2. 83)

где р, т — параметры объекта и модели,

соответственно,

имел единственный экстремум в точке m —

= р)

т1пЯ(|А,от)= /ф,/я=/(р)].

(2. 84)

/neQ

 

Здесь Q — область изменения параметров р, и т, f(p) — вза­ имнооднозначная функция. В частном случае f(p )= p . Отсю­ да следует, что необязательно, чтобы структуры объекта и модели совпадали, но они должны быть известными для нахождения f(p). При идентификации распределенного объ­ екта можно использовать модель с сосредоточенными пара­ метрами.

Следующий этап упрощения алгоритма (2. 79) состоит в замене функций чувствительности 5 некоторыми функциями vlt которые вычисляются проще.

Можно показать, что для обеспечения сходимости т к т* достаточно, чтобы выполнялись условия (одномерный случай):

Игл

1 Г

dw

 

(2. 85)

0—>оо

0 J

дЯ\,

 

 

 

О

 

 

 

= 0

при

=

( 2. 86)

 

> 0

при т > т *,

 

 

 

< 0

при т<т*,

 

rip.mi^rip.nii)

при m j> m 2\

 

| r(p,m) | </Vi=const при

| m | </Va, I {a I <iV2-

Аналогичные условия можно выписать и для многомерного случая. Условия (2. 85), (2. 86) позволяют упростить и в от­

142


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

2. 6. СХОДИМОСТЬ АЛГОРИТМОВ ОЦЕНКИ ПАРАМЕТРОВ

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

прерывных

распределенных

систем надо

быть

уверенным,

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

Вначале

предположим,

что решение

задачи

идентифика­

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

ется сходимость алгоритма стохастической аппроксимации.

Для простоты рассмотрим случай, когда ц — скаляр.

 

 

Теорема.

 

 

 

Алгоритм (2. 57) сходится с вероятностью, равной едини­

це, и в среднеквадратическом и дает

несмещенную

оценку,

если:

 

(2.

87)

1) М{(у J —e f | Xs}<&(1+X2S), fe = const>0;

2) 3=Хв/И{ут / - / Д 8» 0 ;

 

 

 

E= 0 лишь при 7s=(m [s]—fi.)= 0 .

 

 

Доказательство.

(2. 57)ц, возведем

обе

Вычтем из левой и правой частей

части в квадрат и возьмем условное математическое ожида­

ние при фиксированном векторе Xs= || Xi,...,Xs ||

Т

\ ty=X s* -2 T[s]XsAf{Vm/ - e 1~М +

 

+ 7г[ф И {[у / ш- е]2 I U -

(2 - 88)

.1.43


Можно показать, что в условиях

теоремы процесс Xs2'

является

нижним

иолумартингалом

и

в

соответствии с

[2. 38] jM{X2s} сходится

к некоторому

пределу Доо с ве­

роятностью

единица,

т.

е- НшЖ(Хг5)= д < с о ,

и с вероятно-

 

 

 

S —*oj

00

 

 

стью единица существует предел

 

 

 

lim(/?z[s]—|а)=(/я до —|j.)2< co .

5—»оо

Используя условие 1) теоремы, получим из (2. 88)

M{X2S+1| Х5}<Х2, - 2 Т[5]Х,М{ у J - e | Xs}+ T2[# [l+ XЛ 1

Ж{Х%+1}<ДИХ12} - 2 2 - г [ / ] ( Л ^ ( У т ^ - е | Х ;))+

/= 1

 

+ b £J i f[i](\+M {\*i]).

 

 

(2.

89)

 

/ = 1

 

 

 

 

 

Отсюда,

учитывая ограничения

(1. 101)

на

переменный

коэффициент, условие 2) теоремы

 

f fs + 1 ]

,

и что 11т ■- м

= 1

 

S

S—>оэ

Tlsl

 

 

 

 

 

 

 

 

выводим, что р я д ^ 7 [/]-/И(Х/-М{ у то/ —g|X/))

сходится,

имеет

предел и

 

 

 

 

 

 

 

11mM{\sM { y mI—e | XsH =0.

 

(2.

90)

 

S'—>00

 

 

 

 

 

Используя

условия 2) теоремы,

выводим,

что

lim

m[s] = ^,

 

 

 

 

5 —>оо

 

 

т. е- алгоритм (2. 57) дает оценку /n[s], сходящуюся к истинному значению ц с вероятностью единица и в средне­ квадратическом.

Если объект линейный, то из (2. 87) следует, что метод стохастической аппроксимации неприменим в том случае, когда функция потерь растет быстрее квадратичной парабо-

144