Файл: Брикман, М. С. Аналитическая идентификация управляемых систем.pdf

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

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

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

Добавлен: 22.10.2024

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

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

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

ГЛАВА III

76

 

Непрерывность h(a) вытекает из полученной оценки для ||/?а~1|

иследующих соотношений:

\h (a l) —h (a 2) | = | ||ша1—ш0||— \\wai—w0\\| <

Wwa—w^WII (Я аГ 1— /?«*-*) foil ^ | а 2— а 4| х

Х11^аГ1^1--^а2_11' llfoll ^

2— cci) -llfoll

l|/-/?oi?a2- 1ll-

(oCiPl+ Po) «2

Доказательство второго утверждения леммы сводится к построе­ нию оценки для нормы элемента R<s(wa— w0):

\\Ro(Wa,~w0) |= ||a/?i^J|<a||/?ia>0||+a||#i(twa—™o) II —

 

=

a||i?i^oll+all#i#a_1ajRiZ0oll^

 

 

^ a | | i?ity o | | (1 +

Ц / — RoRa~M I).

( 3 . 9 )

Правая часть выражения (3.9) равна нулю при а = 0. Отсюда

l i m

Ro(wa— Wo) = R o ( l i m Wo,—Wo) = 0 .

 

a-»0

 

a-*0

 

Следовательно,

limi0a= ®o,

так как в противном

случае

уравнение (3.8)

 

a->-0

 

 

имело бы более одного решения.

 

Объединяя полученные выше результаты, можем сформули­ ровать условия устойчивости и сходимости регуляризованного метода Галеркина—Петрова при решении задач идентифика­ ции.

Теорема 3.3. Пусть выполнены условия теоремы 3.1 и леммы 3.1, а координатные системы |фг} и {ф*} сильно минимальны в гильбертовом пространстве W. Тогда при любом f0^ W най­ дутся натуральное число п0 и положительные числа а и о та­ кие, что приближенное решение г5аи задачи идентификации по методу Галеркина—Петрова устойчиво по отношению к малым

возмущениям системы (3.4) при п ^ п 0 и ||Г„||^астУЯоро, причем

существует такая

зависимость

аи = а(н, ||Г„||, ||уп||), что

IIWann- w 0II—>-0 при а„~*0, п-+ оо,

IIAlIbn), |[уп|—Д).

Здесь Яо = Пт Ап; ро= Пт ци.

 

П—>аэ

п~>оо

 

Доказательство. Воспользуемся неравенством треугольника для оценки погрешности приближенного решения

\\wan—®ol|sSl|z5a«— Wan\+ \\wan—0>all + l|a>a— Я>о||. (3.10)


77

ПРОЕКЦИОННЫЕ МЕТОДЫ РЕШЕНИЯ УРАВНЕНИЙ ИДЕНТИФИКАЦИИ

Из сильной минимальности координатных последовательнос­ тей следует (определение 1.43.5), что Яо>0, р0>0.

Наименьшие собственные числа матриц Грама ограничены снизу [1.20]: ЯпЗ=^о, рп^ЦоПоэтому из оценки (3.7) получаем

при ЦГпЦ^аоУЯоЦо

\ \ w a n —

- - - \\гп\\+

а2УЯ0ро аУро

Из этого неравенства и теоремы 3.2 вытекает, что при доста­ точно малых по норме возмущениях метода Галеркина—Пет­ рова, удовлетворяющих для произвольного числа е>0 условиям

llrii^s

оьУцое

 

 

Зс2

 

 

 

Ц/,п|| = ; а у Я 0р о - m i n | сг,

а

в

c2llvnll

 

3

5 - ) } , (3.11)

 

 

аУро

существует натуральное число п01 такое,

что при п ^ п 01

11йаи-ш аи||г^-|-.

(3.12)

Воспользовавшись теоремой 3.1, находим натуральное число п02, удовлетворяющее требованию

\\wan—OJcclK - у . (3.13)

Из утверждения леммы 3.1 следует, что существует положи­

тельное число <хо такое, что для всех 0 ^ а ^ а 0

 

 

\\Wa-

(3.14)

Таким образом, для всех возмущений, удовлетворяющих

условиям (3.11),

можно выбрать такое натуральное число п0 =

= шах {п0и «ог},

что в силу (3.12) — (3.14) для всех

п ^ п 0 и

0<аг^сс0 неравенство (3.10) даст следующую оценку:

\\wan—^oll^e.

Отсюда ввиду произвольности е вытекает доказываемая теорема.

Мы свели решение задачи идентификации по методу Галер­ кина—Петрова к задаче аппроксимации элемента Ri~l при


ГЛАВА Ш

78

 

помощи линейной комбинации функций rn = R\~i Яафг

Система функций {%г} линейно независима, так как в против­ ном случае нашлись бы (определение 1.2) такие не равные нулю одновременно постоянные лц , . . . , vn, что

п

ViXi — Rr'Ra ^

VjCpi---0.

2=1

 

Отсюда следует, что однородное

уравнение Raw = 0 имеет

п

 

нетривиальное решение w = V f(p,, ЧТО НеВОЗМ ОЖ НО .

1=1

Обозначим через Хп подпространство, натянутое на элементы

Xl >•••>%П'

Теорема 3.4. Пусть

Q(Gn,Hn) < 1. Тогда при а > 0

решение

задачи идентификации

по

методу Галеркина— Петрова

опреде­

ляется формулой

 

 

 

 

п __

 

(RrH, Ф О

Ror'RiXi,

(3.14)

1

Т«

 

 

 

гДе {%г} и {г|ц} — ортонормальные базисы в подпространствах Хп и Нп соответственно такие, что

(Хь Хг) =

(%й to)

-8{j

(Фг>to)

(i,

0 < т г^ 1 ) .

(3.15)

Доказательство. Система уравнений (3.4) однозначно раз­

решима и эквивалентна следующей:

 

 

X i (Хо tyk)wi= (RrH, Фа),

{k=\,n) .

(3.16)

2=1

Применяя к (3.16) теорему 1.13, получаем, что

0(Х„, Нп) < 1.

Таким образом, подпространства Хп и Нп /г-мерные, а их раствор меньше единицы. В этих условиях, согласно лемме 17.3

работы [1.27], существуют ортонормальные базисы {%i}a:Xn и

{tyi}czHn, удовлетворяющие условиям (3.15). Выбирая постро­ енные базисы в качестве координатных систем и используя ме­ тод Галеркина—Петрова, получаем приближенное решение за­ дачи в виде (3.14).


79 ПРОЕКЦИОННЫЕ МЕТОДЫ РЕШЕНИЯ УРАВНЕНИЙ ИДЕНТИФИКАЦИИ

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

пользовании

метода оказывается достаточным

в

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

со_ структурой уравнений

(3.16)

построить такие

базисы

{%*} и

{фг}, чтобы

выполнялось

лишь

последнее из

равенств

(3.15),

т. е.

 

 

 

 

 

 

Ф ; } = т Л >

Формула (3.14) для решения остается в силе и в этом случае.

3.2. ПРИМЕНЕНИЕ МЕТОДА БУБНОВА— ГАЛЕРКИНА ДЛЯ РЕШЕНИЯ УРАВНЕНИЙ ИДЕНТИФИКАЦИИ

В ряде случаев условия устойчивости и сходимости метода Га- леркина—Петрова (теорема 3.3) могут оказаться ограничитель­ ными при практической реализации. Например, требование сильной минимальности координатной системы не выполняется

для системы

степеней

{th} (& = 0, 1 ,2 ,... ) в

пространстве

L2 (0,1) [1.21],

поэтому

использование последней

для решения

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

Подобное положение имеет место и для некоторых других часто употребляемых систем функций. Кроме того, во многих случаях априорная информация об исследуемой системе не со­ держит ограничений на координатную систему {фг}. В подобных случаях для определения коэффициентов разложения (3.1) следует применять метод Бубнова— Галеркина.

Теорема 3.5. Пусть координатная система {фг} полна и ли­ нейно независима в сепарабельном гильбертовом пространстве W, а параметр регуляризации а>0 . Тогда невязка RaWa.n — f и погрешность wan— RaTlf приближенного решения задачи иден­ тификации по методу Бубнова—Галеркина стремятся к нулю при п-*-оо и любом f e F , если коэффициенты {да*} выражения (3.1) определяются из следующей системы линейных алгебраи­ ческих уравнений:

П

 

У*, ( R r lRa(fi,(f>h)Wi=(Rrif,(fh), (k— l,n).

(3.17)

2=1

 

Доказательство. Оператор R r lRa= aI + R r 1Ro является пра­ вильным (определение 1.43.3), так как оператор R r lRo компак­ тен. Кроме того, из полноты координатной системы {фг} выте­


ГЛАВА HI

8 0

кает предельная плотность последовательности конечномерных подпространств {Gn} в W. Следовательно, выполнены все усло­ вия теоремы 1.11 и теорема 3.5 доказана.

Сравнивая условия теорем 3.5 и 3.1, мы видим, что метод Бубнова— Галеркина обеспечивает сходимость приближенных решений уравнений идентификации при менее ограничительных условиях, чем метод Галеркина— Петрова. Однако в общем слу­ чае порядок системы (3.17), необходимый для получения задан­ ной точности решения, оказывается выше, чем порядок системы (3.4) при тех же условиях.

Предположим теперь, что уравнения (3.17) вычисляются приближенно, т. е. нам приходится решать уравнения методом

Бубнова— Галеркина с возмущениями

 

awan+ (P e ^ r ^ o + S n ) w*n+ g n= P o nR r %

(3.18)

Здесь wan— / ,Wj(pi — решение возмущенного уравнения (3.18) ;

i=i

PGn — ортопроектор на подпространство Gn (определение

1.34).

Оператор Sn и элемент gn описывают произвольные возму­ щения галеркинского уравнения

aWan+PGnR r 1RoWan = PGnR r if,

соответствующего системе (3.17).

Теорема 3.6. Пусть выполнены условия теоремы 3.5. Тогда найдутся не зависящие от элемента f^ F натуральное число п0 и положительное число ап такие, что при всех п ^ п 0 и ||Sn||<an уравнение (3.18) однозначно разрешимо и

l\wan—®a"||^C32l|5„|M|PGrai?r1/||+c3||g-n||,

(3.19)

где

 

 

 

a"

WI-RaT'RoW

II(Pg" /)/?1

 

СЗ a„-||S„|| ;

числа Pi и Po определены в лемме 3.1.

Доказательство. Оператор R^Ra по условию непрерывно обратим, а оператор (I — PGn) вследствие полноты координат­ ной системы стремится по норме к нулю при п—>-оо. Следова­ тельно, для всех п ^ п 0 уравнение (3.18) однозначно разрешимо