Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Брикман, М. С. Аналитическая идентификация управляемых систем

.pdf
Скачиваний:
12
Добавлен:
20.10.2023
Размер:
7.91 Mб
Скачать

21

ОСНОВНЫЕ СВЕДЕНИЯ ИЗ ФУНКЦИОНАЛЬНОГО АНАЛИЗА

сходимости эквивалентны. Однако, если пространство бесконеч­ номерно, то мы можем лишь утверждать, что любая сильно схо­ дящаяся последовательность сходится и слабо. Обратное утвер­ ждение имеет место лишь в некоторых бесконечномерных про­ странствах.

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

Определение 1.25. Множество Т в линейном нормированном пространстве Е называется слабо компактным, если из любой бесконечной последовательности элементов х\ ,. . ., хп , . .. этого множества можно выделить слабо сходящуюся подпоследова­ тельность Хщ , . . ., xnk .

Теорема 1.1 [1.13]. Если линейное нормированное простран­ ство сепарабельно, то любой шар в сопряженном пространстве Е* слабо компактен.

Таким образом, из теоремы следует, что построение слабо компактных множеств в сепарабельных пространствах может быть осуществлено элементарными методами.

Пространство Е* является банаховым пространством, следо­ вательно, можно ввести понятие пространства (£*)*, сопряжен­ ного к Е*.

Определение 1.26. Банахово пространство В называется реф­ лексивным, если оно совпадает с пространством (В*)*.

Конечномерное пространство, а также пространства 1Р, Lp(a,b) при р > 1 и все гильбертовы пространства рефлексивны. В гильбертовом пространстве понятие слабой сходимости фор­ мулируется особенно просто.

Последовательность элементов {xn}cz Н слабо сходится к элементу х е Я , если (хп, г/)->(х, у) для любого i/е Я . Отсюда и из рефлексивности гильбертова пространства вытекает, что из любого ограниченного по норме множества элементов простран­ ства Н можно выделить слабо сходящуюся последовательность.

Ниже, по аналогии с гильбертовым пространством, мы будем пользоваться обозначением < х , /> для значения линейного функционала f<=E* на элементе х е £ .

1.4. ЛИНЕЙНЫЕ ОПЕРАТОРЫ1

Определение 1.27. Если каждому элементу х множества D про­ странства X поставлен в соответствие некоторый элемент у = А х из пространства У, то говорят, что на D задан оператор А со значением в У. Множество D ( A ) = D называется областью опре-

1 Излагается по работам [1.1, 1.9, 1.18, 1.28, 1.37].

ГЛАВА I

22

деления оператора Л, а

совокупность R(A) всех элементов у

из У вида у = А х (xeZ))

называется областью значения опера­

тора А.

 

Из этого определения ясно, что функционал является част­ ным случаем оператора, когда пространство У представляет собой действительную ось. В любом линейном пространстве имеется тождественный или единичный оператор /, ставящий элементу х в соответствие сам этот элемент: /х = х.

Определение 1.28. Оператор А называется непрерывным ли­ нейным оператором, действующим из линейного метрического пространства X в линейное метрическое пространство У, если

1. А линеен, т. е.

а) D(A) — линейное подпространство;

б) A (|3Xi+ Xx2) = ($Axi + XAx2) — для любых Xu Х2^D(A) и

любых действительных чисел [3 и X.

2. А непрерывен, т. е. для любой точки х0еО (Л ) и произ­ вольного числа в> 0 найдется число 5>0 такое, что ру(Лх, Лх0)< е для всякой точки x^D(A) и удовлетворяющей неравенству рх(х, х0)< 6 . Здесь ру и р х — метрики в простран­ ствах У и X соответственно.

Определение 1.29. Линейный оператор Л, действующий из нормированного пространства Е в нормированное пространство F, называется ограниченным, если ЦЛхЦр^сЦхЦд, причем вели­ чина с не зависит от выбора

Наименьшая из постоянных с, удовлетворяющих этому нера­ венству, называется нормой оператора Л и обозначается ||Л||.

\\А\\= sup —

IHH-E

— ~

SUP

X(= E

Ы

Е = 1

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

Между понятиями ограниченности и непрерывности линей­ ных операторов существует тесная связь.

Теорема 1.2 [1.13]. Для того чтобы линейный оператор, дей­ ствующий из Е в F, был ограничен, необходимо и достаточно, чтобы он был непрерывен.

Рассмотрим пространство всех линейных операторов. В этом пространстве можно ввести понятие сходимости последователь­ ности операторов.

Определение 1.30. Последовательность ограниченных опера­ торов {Ап}, действующих из нормированного пространства Е в

23 ОСНОВНЫЕ СВЕДЕНИЯ ИЗ ФУНКЦИОНАЛЬНОГО АНАЛИЗА

нормированное пространство F, сходится к линейному опера­ тору А из Е в F:

1)

по норме, если lim ||Л—ЛП|= 0;

2)

сильно,

71-»сО

если lim \\Ах — Апх\\Р= 0 для всех х е £ ;

3)

слабо,

п —*оо

если при любом х е £ последовательность {Апх}

слабо сходится к Ах.

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

Если линейные операторы A t и А2 действуют в одном и том же нормированном пространстве Е, то для них вводится опера­ ция умножения по правилу Л =Л 1Л2, если Л х=Л 1Х (Л 2х).

В общем случае Л1Л2=Л=Л2ЛЬ т. е. умножение операторов некоммутативно. Имеет место неравенство

1!Л1Л2|!^ЦЛ1!М|Л2||. (1.7)

При решении операторных уравнений фундаментальное зна­ чение имеет понятие обратного оператора.

Определение 1.31. Если оператор Л действует из простран­ ства Е в пространство F и существует оператор из F в Е та­ кой, что для всех х е £ и у е £ имеют место соотношения:

А~1А х = х , АА-^у— у,

(1.8)

то операторы А и А~1 называются взаимно-обратными. Используя определение единичного оператора из (1.8), по­

лучим Л-1 Л= ЛЛ~1= /.

Допустим, что пространства Е и F нормированы, и рассмот­ рим ограниченный линейный оператор Л, действующий из Е в F,

и линейный функционал g(y) = < y , g > ,

y^F, g^F*, опреде­

ленный на F. Допустим, что у=Ах, хеЯ . Тогда

g ( y ) = < A x , g > = f ( x ) = < x , f y ,

f<=E*.

Это равенство ставит в соответствие каждому линейному функ­

ционалу g^F*

некоторый линейный функционал f^E*.

Определение

 

1.32. Оператор Л*, определяемый равенством

<Ах,

g > = <x,

f> = < x,

A*g>,

называется оператором,

сопря­

женным к оператору Л.

Если оператор Л действует из £ в £*

и для

любых

х, у ^ Е

имеет

место равенство <х,

А у > =

= <у ,

Ах >, то оператор Л называется симметрическим.

||Л||=

Сопряженный

оператор линеен и ограничен, причем

= |(Л*[| и (Л*)—1=

(А~1)* тогда и только тогда, когда Л-1 сущест­

вует.

В гильбертовом пространстве Н сопряженный оператор вво­ дится при помощи скалярного произведения (Ах, у) — (х. А*у) ,

ГЛАВА I

24

x,yi=H. Симметрический оператор в данном случае определя­ ется условием (Ах, у) = (х, Ау) — (Ау, х) и называется самосо­ пряженным. Очевидно, что для самосопряженного оператора А выполняется тождество А=А*.

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

Из полной непрерывности линейного оператора вытекает его непрерывность. Обратное утверждение в общем случае места не имеет. В частности, единичный оператор / непрерывен, но не компактен, если В бесконечномерно.

Теорема 1.3 [1.13]. Справедливы следующие утверждения:

1)линейная комбинация компактных операторов является компактным оператором;

2)произведение компактного и ограниченного операторов есть компактный оператор;

3) оператор, сопряженный к компактному, компактен;

4) предел последовательности компактных операторов, схо­ дящихся по норме, есть компактный оператор;

5) оператор, обратный компактному, не ограничен. Сформулированная теорема определяет основные свойства

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

Среди самосопряженных операторов наиболее простую структуру имеют проекционные операторы.

Теорема 1.4 [1.13]. Если L является подпространством Я, то любой элемент х е Я единственным образом может быть пред­ ставлен в виде

x = y + z ,

(1.9)

где y^L, Z-LL (определение 1.18).

(1.9) называется

Определение 1.34. Элемент у из формулы

проекцией х на L, а оператор PL, ставящий в соответствие каж­ дому элементу х его проекцию у на L : y = PLx, называется опе­ ратором ортогонального проектирования, или ортопроектором.

Проекционный оператор обладает следующими свойствами:

1)P*l = Pl ^ P l2-,

2)||PJI = 1;

3) P lP m = 0 в том и только в том случае, когда простран­ ство М является ортогональным дополнением к пространству L; подобные операторы называются ортогональными:

25

ОСНОВНЫЕ СВЕДЕНИЯ ИЗ ФУНКЦИОНАЛЬНОГО АНАЛИЗА

4) ортопроектор PL конечномерен и вполне непрерывен, если подпространство L конечномерно. Если пространство L беско­ нечномерно, то соответствующий ему ортопроектор не вполне непрерывен.

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

1.5. ЛИНЕЙНЫЕ ОПЕРАТОРНЫЕ УРАВНЕНИЯ1

Воставшихся параграфах данной главы (если не указано про­ тивное) символом А обозначается линейный оператор, дейст­ вующий из банахова пространства В в банахово пространство С

и имеющий область определения D(A)czB и область

значений

£(А)с=С.

 

Определение 1.35. Уравнение

 

А х — у

(1.10)

с известными оператором А и свободным членом р е С

называ­

ется линейным операторным уравнением относительно неизвест­ ного элемента х.

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

Определение 1.36. 1) Совокупность всех решений однородного

уравнения

 

(1.11)

 

А х = 0

 

называется нуль-пространством, или ядром N (А)

оператора А.

2)

Уравнение (1.10) называется:

если N ( A ) = 0 ;

а)

однозначно-разрешимым на 7?(А),

б)

корректно-разрешимым на /?(А),

если

для любого

x^D(A) справедливо неравенство ||х||в^/г||Ах||с, причем число k не зависит от х;

в) везде разрешимым, если R (A) = С;

г) плотно разрешимым, если J?(A) = С ;

д) нормально разрешимым, если /?(А) замкнуто: R(A) =

=ЩА).

Очевидно, что для однозначно-разрешимого уравнения на i?(A) существует обратный оператор А~1. Если уравнение кор­ ректно разрешимо, то его решение непрерывно зависит от сво­

1 Излагается по работам [1.9, 1.16, 1.34, 1.37].

ГЛАВА I

26

бодного члена у, так как корректная разрешимость эквивален­ тна ограниченности оператора Л-1. Смысл остальных определе­ ний разрешимости ясен и зависит только от характера области значений -К (Л ).

Определение 1.37. Линейный оператор А называется замкну­ тым, если из того, что хп-^х и и/, следует, что x ^ D ( A ) и

Ах = у.

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

Теорема 1.5 [1.16]. Допустим, что оператор А замкнут, а его область определения всюду плотна в В. Тогда существует сле­

дующая связь между свойствами уравнения (1.10)

и сопряжен­

ного к нему

(g^C*;

f^B*).

 

(1.12)

 

A * g = f

 

1)

Уравнение (1.10)

 

 

Уравнение

(1.12)

однозначная разрешимость1

плотная разрешимость; .

2)

корректная разрешимость

< = >

везде разрешимость;

3)

везде разрешимость

 

<=:>

корректная

разрешимость;

4)

плотная разрешимость

 

< = >

однозначная

разрешимость;

5)

нормальная разрешимость

<д=)

замкнутая разрешимость.

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

В теории систем линейных алгебраических уравнений с квад­ ратной матрицей имеется простая связь между решениями основного и сопряженного уравнений [1.35]. Эта связь имеет место и для операторных уравнений при специальном виде опе­ ратора А.

Определение 1.38. Оператор А называется фредгольмовым, если он представим в виде

А = и~'+Т,

(1.13)

где U — ограниченный оператор, определенный на всем С и имеющий обратный U~l; Т — компактный оператор, действую­ щий из В в С. Если в (1.13) U = I при В = С, то оператор А называется каноническим фредгольмовым.

Очевидно, что уравнение (1.10) с фредгольмовым оператором А корректно разрешимо. Остальные условия разрешимости опи­ саны в следующей теореме.

1 Запись а

Ь обозначает, что из а следует b а запись а ( = > Ь,

что а к Ь эквивалентны.

27

ОСНОВНЫЕ СВЕДЕНИЯ ИЗ ФУНКЦИОНАЛЬНОГО АНАЛИЗА

Теорема 1.6 [1.9]. Допустим, что оператор Л является кано­ ническим фредгольмовым в банаховом пространстве В. Тогда 1) уравнение (1.10) разрешимо в том и только в том случае, когда < £ / , / > = 0 для любого (еЛ * и удовлетворяющего сопря­

женному однородному уравнению

 

A * f = 0;

(1.14)

и

2) число линейно-независимых решений g уравнений (1.10)

(1.14) одинаково.

 

к

Любое фредгольмово уравнение может быть преобразовано

каноническому фредгольмову, поэтому

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

с соответствующими изменениями справедлива для произволь­ ных фредгольмовых уравнений. В связи с этим рассмотрим воз­ можности преобразования уравнения (1.10) к каноническому фредгольмову.

Определение 1.39. 1. Линейный ограниченный оператор S, действующий из С в банахово пространство F, называется ле­ вым эквивалентным регуляризатором для оператора А, если оператор 5Л — канонический Фредгольмов, D (S)z^R(A) и N (S) = 0 .

2. Ограниченный линейный оператор V, действующий из ба­ нахова пространства G в В, называется правым эквивалентным регуляризатором для оператора А, если оператор AV — кано­ нический фредгольмов и D(A)c^R(V).

Из этого определения следует, что в результате эквивалент­ ной регуляризации уравнение (1.10) с произвольным операто­ ром может быть преобразовано к эквивалентному каноничес­ кому фредгольмову уравнению. В частности, для фредгольмова оператора (1.13) оператор U является эквивалентным левым

иправым регуляризатором, если R(U) zdD(A) .

Взаключение данного параграфа рассмотрим понятие спек­

тра оператора, которое связано с разрешимостью уравнения

Ах—к х = (Л—Х1)х=у,

(1.15)

где к — некоторое произвольное число.

Определение 1.40. Число а называется регулярным значе­ нием оператора А, если уравнение (1.15) корректно и плотно разрешимо. Совокупность всех регулярных значений называется резольвентным множеством оператора А, а дополнение к нему на комплексной плоскости — спектром оператора А. Оператор

Rx(A) — {А — kl)~l, определенный для регулярных

а, называ­

ется резольвентой оператора А.

 

Пример 1.1. Интегральное уравнение

 

%x{t)+ JЬ K(t,x)x(x)dx=y{t)

(те)

а

ГЛАВА I

28 '

называется

интегральным уравнением Фредгольма второго рода. Для всех

регулярных X решение уравнения

(1.16), если оно существует,

может быть

определено

при помощи резольвенты r(t, х,Х) ядра K(t,%) по

следующей

формуле [1.34]:

Jь r(t,x,X)y(x)dx].

 

 

=

(1.17)

' а

Спектр ограниченного оператора лежит в круге |Я|^||Л|| и всегда не пуст. Радиус наименьшего круга с центром в начале координат, содержащего спектр оператора А, называется спект­ ральным радиусом г (А) оператора А. Справедлива следующая формула [1.371:

г(Л) = Нт У1И||"^М||.

Если |Я|>г(Л), то резольвента представляется сходящимся по норме рядом

 

« 1(л , = - ф ( ; + ф л +

 

 

 

 

 

г

ь

ь

Пример 1.2.

Допустим, что в уравнении

(1.16) |Я|>

/ /

КЦ1, т)Х

 

 

 

*-п а

 

X dtdzJ

Тогда резольвента f(t, т, Я) представляет собой сильно сходящийся

¥

 

 

 

 

 

в L2[a, b) X (а, Ъ)]

ряд1 [1.34]

 

 

 

 

 

r(t, т Д )= X ( - а)-^КДТ т),

 

(1.18)

 

 

/=1

 

 

 

в котором итерированные ядра Kj{t*т ) определяются рекуррентными фор­ мулами

 

Kl(t,x)=K(t,x)

 

Kn+m(t,x)= Jъ Kn(t,&)Km(@,x)d@

(т, П= 1, 2, . . . ) .

 

а

 

 

Пример 1.3. Интегральным оператором Вольтерра называется оператор

Ax(t) = Jt

K(t,x)x(x)dx .

(1.19)

__________________ __________________________________

a

 

 

1 X x Y называется прямым произведением пространств X и Y и обозна­ чает совокупность всевозможных пар {х, у}, х<=Х, г/еУ.

29

ОСНОВНЫЕ СВЕДЕНИЯ ИЗ ФУНКЦИОНАЛЬНОГО АНАЛИЗА

Если K(t, Т)

непрерывна в

треугольнике

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

С(а,Ь) спектр оператора (1.19)

состоит из одной точки X=0.

Определение 1.41. Число X называется собственным значе­ нием оператора А, если однородное уравнение Ах = лх имеет при данном X нетривиальное (отличное от нуля) решение. Это ре­ шение называется собственным элементом оператора А, соот­ ветствующим собственному числу X.

Смысл собственных чисел и собственных элементов стано­ вится особенно ясным, если допустим, что оператор А явля­ ется самосопряженным в гильбертовом пространстве Н.

Определение 1.42. Нижней и верхней границами самосопря­ женного оператора А называются соответственно числа т —

= inf (Ах, х)

и M = sup(Ax, х).

11x11=1

||х!1 = 1

Самосопряженный оператор А называется положительным,

если m^sO,

и положительно определенным, если т > 0.

Для нормы самосопряженного оператора справедливо соот­ ношение \\А|= т а х { |от , |} =sup |(Ах, х) |.

||Х || = 1

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

1.6. МЕТОДЫ РЕШЕНИЯ ОПЕРАТОРНЫХ УРАВНЕНИЙ1

Допустим, что существование и единственность решения урав­ нения (1.10) установлены. После этого основной задачей явля­ ется приближенное в общем случае отыскание элемента х, удовлетворяющего рассматриваемому уравнению. Для этой цели конструируется алгоритм, позволяющий определять по­ следовательность {хп} а В приближенных решений таких, что

limxra = x. Фактически, мы можем найти лишь конечное число

П-^ОО

членов последовательности {х"}. Поэтому кроме исследования сходимости рассматриваемого алгоритма должна быть по­ строена также оценка погрешности ||хп —х||в, позволяющая определить степень близости приближенного решения к точному.

Затем необходимо исследовать устойчивость используемого алгоритма, т. е. влияние на него ошибок вычислений и погреш­ ностей задания оператора А и свободного члена у.

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

1 Излагается по работам [1.2, 1.20, 1.25, 1.27, 1.37].

ГЛАВА I

30

методов построения последовательности {хп}. Эти методы можно разбить на две группы:

1)методы последовательных приближений (итерационные);

2)методы сведения к более простым уравнениям (проек­ ционные) .

Изложение итерационных методов начнем с того случая,

vкогда уравнение (1.10) представляет собой операторное урав­ нение второго рода, т. е. А = 1+ К.

Последовательные приближения описываются формулой

хп = у —Кх*-\ (п = 1,2, ... ).

(1.20)

Здесь хп — значение приближенного решения на п-м шаге итерации.

Если спектральный радиус оператора К р( К)<1, то при любом начальном приближении х° итерационный процесс (1.20) сходится и имеет место следующая оценка [1.27]:

||*»-*||<с(е) [р(/С) +е]Ч х°+ К х°-у\\,

(п = 0, 1,2, ...),

где е — произвольное число, 0 < 8 < 1—р(/();

с(е) — постоянная, зависящая от е.

Врассматриваемом случае последовательные приближения (1.20) сходятся со скоростью геометрической прогрессии.

Пример 1.4. В условиях примера 1.2 решение уравнения (1.16) может быть определено при помощи следующего сходящегося итерационного про­

цесса:

I

 

*П(0=Нг

’ («=1.2,...); X0(t) = j - y { t ) .

1

о

-*

Пример 1.5. Интегральное уравнение Вольтерра второго рода с непре­

рывным ядром K(t, т)

 

 

 

t

 

 

 

x (t)+ J K(t, x)x(i)dx=y(t),

(a^t,

(1.21)

a

 

 

 

имеет единственное непрерывное решение x(t)

при любой непрерывной функ­

ции y(t) и конечном отрезке

[а, Ь], так как

р(К)=0. Это решение

явля­

ется пределом последовательных приближений

 

 

t

 

 

 

x n(t)=y(t) — j K ( t , x ) x n~i(x)dx,

(n = 0, 1,2, . . . ) ,

 

a

 

 

 

причем

 

 

 

\\xn- x \\c(a

Kn( b - a )"

,(n= 1,2,...).

 

 

n\

 

 

Соседние файлы в папке книги из ГПНТБ