Файл: Брикман, М. С. Аналитическая идентификация управляемых систем.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 22.10.2024
Просмотров: 121
Скачиваний: 0
ГЛАВА III |
76 |
|
Непрерывность h(a) вытекает из полученной оценки для ||/?а~1|
иследующих соотношений:
\h (a l) —h (a 2) | = | ||ша1—ш0||— \\wai—w0\\| <
Wwa—w^W — II (Я аГ 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) однозначно разрешимо