Файл: Живоглядов, В. П. Адаптация в автоматизированных системах управления технологическими процессами.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