Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология
Образование Политология Производство Психология Стандартизация Технологии


Тогда справедливы все утверждения теоремы 1.



Для доказательства заметим, что по лемме 1 п. 32.2 опера­тор Ф удовлетворяет в Q условию Липшица с постоянной q е е(0, 1), т. е. является в Q сжатием.

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

Пусть Н — т-мерное унитарное пространство, у — заданный вектор в Я, а Для нахождения решения уравнения

х — Вх = у (П

часто применяется так называемый метод простой итерации (см. [1]). При этом решение уравнения (1) разыскивается как предел последовательности

хп+1 = Вхп + у, п = 0, 1,..., (2)

а начальное приближение хо задано. Если хп-»-х при л-> оо (х„ — решение (2), а х — решение (1)), то говорят, что метод простой итерации сходится.


С точки зрения принципа сжимающих отображений уравне­ние (1) следует записать в виде дг = Ф(дг), где Ф (х) = Вх + у. При этом

|| Ф ([6]') - Ф (*" ) || = || Вх' - Вх" || < || В || || - х" ||.

Если ||В||< 1, то оператор Ф является сжимающим и метод простых итераций сходится. Более общее условие его сходимо­сти дает теорема 2 п. 33.2, когда сжимающим является опера­тор Ф* при некотором натуральном k. Предлагаем читателю проверить, что в данном случае

k

ф* (х) = в*х+Е в1 у

/ = 0

и, следовательно,

|| Фк (*') - Фк (х" ) || < || Вк || || - х" ||.

Таким образом, если удастся показать, что при каком-либо k ||Bfe||< 1, то, согласно теореме 2 п. 33.2, метод простой ите­рации сходится. На этом пути удается получить следующее не­обходимое и достаточное условие (см. [1]).

Теорема. Для сходимости метода простой итерации при любом начальном приближении необходимо и достаточно, что­бы все собственные значения оператора В по модулю были меньше 1.

Доказательство достаточности. Пусть Яь,....., Хт — собственные значения оператора В (среди них могут быть равные). Согласно условию теоремы существует число <? е(0, 1) такое, что

г К<?, i= 1, .... т.

Фиксируем число б > 0 настолько малое, что < 7 + 6 < 1, и рассмотрим оператор 6Его собственные значения равны б-1Я„ г = 1, т. По теореме Жордана (см. п. 23.2) суще­ствует в Н базис {et}™, состоящий из собственных и присоеди­ненных векторов оператора б~1В. При этом выполняются сле­дующие соотношения:

& ~lBei = Я16-1е1, 6-1Ве, = Я; б-1ег + г = 2, ..., т,

где либо аг-i = 0, либо — 1. После умножения на б при­ходим к равенствам

Вех = Ххеь Bei = %iel-{-bai_lei_u г = 2, т. (3)

Возьмем х е Н, и пусть его разложение по базису имеет

вид

т

Пользуясь формулами (3), можно показать, согласно мате­матической индукции, что

т т— 1

Bx—Y, Kli^i + 6 £ а; 1г + 1<? «,

г=i i=i

m m—1 m—2

b2x = Z Aft*, + 26 2 A + б2 У i+qet, (4)

l — l I I » /

m m—I m—&

Фа + С», I + Z

где Cm — биномиальные коэффициенты. При этом, если m^s> то полагаем

Вспомним теперь, что | ^ < 7, и, пользуясь неравенством Коши — Буняковского, получим

I

I m—s II m—s

11 s < 3)

где

C

m Sl/2 / m ч 1/2

Slid2] У^И2; • (6)

Оценка (4) с использованием (5) приводит к следующему не­равенству:

|]B**ll< fo + 6)*a(*)- (7)

Заметим теперь, что функционал а(х), определяемый фор­мулой (6), удовлетворяет всем аксиомам нормы. Но Н конеч­номерно, и потому все нормы в Я эквивалентны. Поэтому най­дется постоянная О 0 такая, что для всех л: е// имеем c||jc||. Обращаясь к (7), получаем оценку

т. е. c(q + b)k. Так как < 7 + б< 1, то при достаточно

большом k будем иметь ||£ ft]|< 1, и метод простой итерации сходится.

Доказательство необходимости. Пусть X соб­ственное значение В и 1. Пусть < р — собственный вектор В, отвечающий %. Выберем в (2) начальное приближение ха = = х + ф, где х — решение уравнения (1). Тогда

х\ — Вх о + У — Вх + У + М> = х + Яф,

х2 = Вхх + у = х + №у, *rt = * + A" q>.


При | Я хп-/*хпри п-*- оо, и, следовательно, метод ите­рации расходится. Теорема доказана.

Замечание I. Если В — самосопряженный оператор в Я,

то ||S]| =max 11 ^ q < 1 (см. следствие к теореме п. 23.4) и 1

оператор Ф (х) = Вх + у — сжимающий.

Замечание 2. Теорему можно перефразировать так: для того чтобы спектральный радиус г0(В) оператора В был мень­ше 1 (см. п. 24.2), необходимо и достаточно, чтобы все соб­ственнее значения В были по модулю меньше 1.

Рассмотрим теперь метод верхней релаксации (см. [14]). Для решения линейного уравнения

Ах = у, (8)

где у е Н, А ^ 2(H), Н — m-мерное унитарное пространство, обычно используются различные варианты его записи в форме (1). Опишем один из них. Пусть в Н фиксирован ортонормиро- ванный базис, в котором (4) является матричной записью си­стемы линейных уравнений:

*=(ЫГ> г/ = ЫГ> A = (an)7.i=v

Предположим, далее, что оператор А вещественный самосо­пряженный, т. е. Oil вещественны и ач = а; 1, и положительно определенный, т. е. (Ах, х) > 0 при х ф 0. Разложим матрицу А на сумму трех матриц:

Л = С + Л + 0, (9)

где матрица С — нижняя левая треугольная, матрица D — верхняя правая треугольная, а Л — диагональная матрица. Точнее, с£ / = ац при i > j, сц — 0 при i ^ /; Я„ = аи, %ц = О при i ф у, dij = Chi при i < }, а^ = 0 при i ^ /'. Запишем уравнение (8) в виде

(Л + «вС)ж = [(1 - < в) Л - < nD] л: + coy, (10)

где со > 0. Поскольку матрица Л + со С обратима, то уравнению (10) можно придать вид

х = Вах -f со (Л + мС)-1 у, (11)

где Вш = (Л + соС) [ (1 — со) Л — coD].

Покажем, что при 0 < со < 2 все собственные значения мат­рицы В(о по модулю меньше 1. Действительно, если Я— ее соб­ственное значение, а ф — отвечающий Я ее собственный вектор, то

[(1 — со) Л — ©D] ф = Я (Л + юС) ф. (12)

Учитывая представление (9), уравнения (12) можно записать в виде

[(2 — ©) Л — юЛ + ю (С -Чр)] ф = Я ((2 - ©) Л + coy* + со (С - £ > )] ф. J)96

Умножим это равенство скалярно на ф и, пользуясь тем, что ((С — D)ф, ф) = 0, получим

, _ (2 — ю) (Лф, ф) — to (Лф, ф) (2 — (в) (Лф, ф) + (о (Лф, ф) 1

Так как (Лф, ф)> 0 и (Лф, ф) > 0, то при со е (0, 2) имеем ]Л| < 1. Согласно методу простых итераций, метод верхней ре­лаксации, т. е. метод последовательных приближений для урав­нения (11), сходится при любом начальном приближении, если О <; со < 2. Число со называется релаксационным множителем. Его выбор может существенно ускорить сходимость последова­тельных приближений (см. [14]).

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

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

^ = F(t, x), (1)

где F(t, x) —нелинейный оператор от двух переменных: веще­ственного переменного /^0 и переменного х из веществен­ного банахова пространства X; значения F также лежат в X. Для уравнения (1) ставится задача Коши, т. е. задается на­чальное условие (а е X)

х 1г_о = а. (2)

Теорема 1. Пусть F(t, x) непрерывен по t на [0, 0] при каждом фиксированном х с \\х —. хоИ^'" " при t е [0, 0] их таких, что IU — удовлетворяет следующим условиям:

JI F (/, х) || ^ с, (3)

||F(t, Xl)-F(t, х2)||< /я||*а-дс, ||. (4)

Тогда на [0, 8j], где

01 = min± 0), (5)

существует единственное решение x(t) задачи Коши (1) — (2), При этом на [0, 0i] ||jc (/)— *о1К г.

Доказательство. Пусть x'(f)—решение задачи (1)—> (2), т. е. непрерывно дифференцируемая на некотором [0, 0i], Gi ^ 0, абстрактная функция, обращающая (1) в тождество на [0, 01 ] и удовлетворяющая начальному условию (2). Под- ставим x{t) в (1) и полученное тождество проинтегрируем на [О, /] с учетом (2). Получим, что x(t) удовлетворяет следую­щему интегральному уравнению:

t

х (0 = а + J F (s, х (s)) ds. (6)

о

Обратно, пусть x(t)—непрерывное на [0, 9i] решение ин­тегрального уравнения (6). Заметим сначала, что при se[0, 0i] абстрактная функция F(s, x(s)) непрерывна на [0, 0i]. Это вы­текает из непрерывности x(s), непрерывности F(t, x) по t и следующей оценки, использующей условие Липшица (4):

|| F (s, х (s)) - F (s0, х (so)) IK IIF (s, x (s)) -F(s, x (s0)) || +

+ || F (s, x (so)) - F (so, x (s0)) IK m || * (s) - * (s0) II +

+ || F(s, x(so))-F(so, x(s0))f

Если s, s0e[0, 0i] и s-> s0, то правая часть в последнем неравенстве стремится к нулю и, значит, F(s, x(s)) непрерывна в каждой точке [О, Эi]. Отсюда вытекает, что согласно (6) лс(0)=а, а также дифференцируемость х(/) на [0, 6i] и равен ство x'(t) = F(t, x(t)) (см. свойство 8 п. 25.2). Итак, x(i) — ре­шение задачи (1) — (2). Таким образом, задача нахождения решения задачи Коши (1) — (2) эквивалентна задаче нахожде ния непрерывного решения интегрального уравнения (6). Поль­зуясь этим, введем банахово пространство Cx[0, 9i] абстракт ных, непрерывных на [0, 0i] функций x{t) со значениями в X и с нормой

(01 = max ||х (0II.

[о, 0, ]

Рассмотрим в Cx[0, 9i] замкнутый шар

Sr (а) = {х (0 е Сх [0, 6J; ||х (/) - a IK г). Нелинейный оператор

Ф(х) = а + \F(s, x(s))ds (7)

о

переводит Sr(a) в себя, ибо

t I t

[F(s, х (s)) ds < шах [ || F (s, x (s)) || ds< 0^ < r J 10.6, 1 J

| Ф (x) —a | = max

о « №.ej

10, 8, ]

вследствие неравенства (3) и определения 0[ равенством (5). Кроме того, Ф(х) является на Sr{a) оператором сжатия, так


как согласно (4)

| Ф (jci) — Ф(лг2)| = тах [о, ej

t

\[F{s, xl (s))-F(s, x2(s))]ds <

0 II


 

 


< max \ II F (s, xx (s)) — F (s, x2 (s)) || ds < [o, ад J

< 0i«|*i (s) — x2(s)\ = q\xx —x2\,

где q = 0i, m < 1 (см. (5)).

Согласно принципу сжимающих отображений в шаре Sr(a)1 существует единственное решение С*[0, 0i] уравнения

(6). Таким образом, теорема 1 доказана.

Определенным недостатком полученного результата является то обстоятельство, что решение x(t) задачи (1) — (2) опреде­лено не на всем [0, 0], а лишь локально на [0, 0i], где 01 удов­летворяет равенству (5). Пример задачи Коши для обыкновен­ного дифференциального уравнения dxjdt — x2, х(0)=1, точ­ное решение которой имеет вид х (I) — _ -, показывает, что

в общем случае нельзя гарантировать существования решения задачи (1) — (2) «в целом»—на всем [0, 0].

При более жестких предположениях относительно опера­тора F{t, x) удается доказать теорему существования и един­ственности решения задачи Коши (1) — (2) на [0, 0].

Теорема 2. Пусть оператор F(t, x) непрерывен по t на [0, 0] при каждом фиксированном х е X и удовлетворяет усло­вию Липшица (4) при этих же значениях переменных. Тогда на [0, 0] существует единственное решение задачи Коши (1) —■ (2).

Мы приведем ниже два доказательства теоремы 2, каждое из которых достаточно просто и поучительно.

Доказательство 1. Как и в доказательстве теоремы 1, сведем задачу (1) — (2) к эквивалентному интегральному урав­нению (6). Рассмотрим в банаховом пространстве С*[0, 0] ин­тегральный оператор Ф(х), определяемый равенством (7), Имеем следующую оценку:

t

[Ф (*, ) — Ф (х2}| < m ^ || (s) — х2 ( s ) || ds < mt | xi — х21, о

используя которую находим

t

|| Ф2 (*, ) - Ф2 2) II < m $ IIФ (*, (s)) - Ф (х2 (s)) II ds <

о

t

< m2 ^ s II Xl (s) — х2 (s) II ds < ~2j I Xl — x2 [.

Продолжая такие оценки, методом полной математической индукции устанавливаем оценку

II Ф" (*, ) - Ф" (х2) [К-^ |*i - х2|. Переходя к максимуму на [0, 0], получим

|Ф" (х, ) - Ф" (х2)К|х, - х21.

Поскольку (mQ)n/nl -*• 0 при я-»- оо, то, фиксируя настолько большое п, чтобы ( mQ)a/n\ = q < 1, придем к выводу, что Фя — сжатие в С* [0, 0]. Утверждение теоремы 2 следует те­перь из теоремы 2 п. 33.2.

Доказательство 2. В банаховом пространстве С*[0, 0] введем эквивалентную норму

I * if) lm = max II e~mix (0 [|.

[0. ei

Покажем теперь, что интегральный оператор Ф является сжатием в смысле этой новой нормы. Действительно, t

|| Ф (х, ) — Ф (х2) IK т $ II Xi (s) — х2 (s) I, e-msems ds < о

t

< m J ds || x, - x2 ||m = (emi - 1) |x, - x2|m.

о

Умножив полученное неравенство на e~mt, получим

e~mt IIФ (х.) - Ф (х2) IK (1 - e~mt) | х, - х2|т.

Наконец, перейдем к шах на [0, 0], что дает неравенство

| Ф (х, ) - Ф (х2) < (1 - е-9) IX! - х2 \т.

Итак, Ф — сжимающий оператор. Теорема 2 доказана. Упражнение. Перефразируйте теоремы 1 и 2:

а) на случай обыкновенного дифференциального уравнения,

б) на случай системы I обыкновенных дифференциальных уравнений с I неизвестными функциями.

Рассмотрим задачу Коши для линейного дифференциаль- t ного уравнения в банаховом пространстве:

=A(t)x + y(t), t> 0, х (0) = а. (8)

Пусть y{t)] и A(t) —соответственно абстрактная функция со значениями в X и оператор-функция из & (Х), непрерывные на m+.ooh,

Поскольку на каждом [0, 8] шах|]Д (0 И< то из тео­ремы 2 вытекает существование и единственность решения за­дачи (8) на полуоси [0, + оо). Доказательство этого утвержде­ния мы предлагаем провести читателю.

§ 34. Итерационный процесс Ньютона

34.1. Описание итерационного процесса Ньютона. Ньютон предложил эффективный метод вычисления решений уравнения

F(x) = 0 ' (1)

для случая функции F(x) с вещественными значениями, зави­сящей от вещественной переменной х. Впоследствии метод Нью­тона был перенесен на системы уравнений (когда F(x)^Em, х^Ет), а затем обобщен в работах акад. Л. В. Канторовича на уравнения в банаховых пространствах.

Пусть F(x) — нелинейный оператор, определенный в окре­стности S решения х* уравнения (1), непрерывно дифференци­руемый в S в смысле Фреше. Пусть, далее, в S оператор F'(x) непрерывно обратим. Итерационный процесс Ньютона состоит в следующем. Выбирается начальное приближение Хо е 5 и ле­жащее достаточно близко к решению х*. Дальнейшие прибли­жения хп, п= 1, 2, ..., предлагается вычислять по формуле

Хп = Хп-1 — [F' (*n-i)r' F (*„_i), п= 1, 2,... (2)

В книге [11] (гл. XVIII) читатель может найти подробней­шим образом разработанную теорию метода Ньютона, различ­ные варианты теорем о сходимости итерационного процесса (2) и близких к нему процессов, а также многочисленные прило­жения. Метод Ньютона является в настоящее время одним из наиболее употребительных вычислительных методов. Главное его достоинство — это (в определенных предположениях) очень быстрая сходимость последовательных приближений (2) к ре­шению х*. Метод применим также и в случае, когда уравнение (1) имеет несколько решений.

Дадим описание метода Ньютона в простейших случаях. Пусть сначала х — вещественная переменная и значения F{x)' также вещественны. Предположим, что в окрестности решения х? (рис. 21) функция F(x) строго возрастает и выпукла вверх (аналогично рассматриваются и другие случаи строго монотон­ной и выпуклой функции). Выберем начальное приближение Хо достаточно близко к х*. Запишем уравнение касательной к кри­вой y = F(x) в точке (xq, F(xq)): у = F(х0) + F'(х0) {х — х0). Точка пересечения касательной с осью абсцисс получится при у = 0. Ее абсцисса равна значению, вычисляемому по фор­муле (2) при п = 1. Далее проводим касательную к кривой у = F(x) в точке (xi, F(xi)) и получаем х2 (см. (2) при п = 2) и т. д. В результате получается последовательность {х„}, очень быстро сходящаяся к х*, если — достаточно мало. В рас­сматриваемом случае метод Ньютона обычно называют мето­дом касательных Ньютона.

Пусть теперь х, F(x)^Em. Если ft(xь..., хт), £ = 1, ......, п, — координатные функции F(x), то уравнение (1) яв­ляется краткой записью системы уравнений

* „.. fi(x ь.... Хт) = 0, i= 1.......................... т.

Производная F' (я) = ~ это матрица Якоби, а


 

|[Р(л: )]-1 — обратная матрица. Формулы (2) представляют со­бою, таким образом, матричную запись итерационного процес­са Ньютона в Ет.

34.2. Теорема о сходимости итерационного процесса Ньюто­на. Приведем один из наиболее удобных в приложениях вари­антов теорем о методе Ньютона. Существование решения х* здесь не предполагается, а доказывается. Вопрос о единствен­ности решения в рассматриваемом шаре здесь не обсуждается.

Теорема. Пусть в шаре Sr(x0) оператор F(x) дифферен­цируем и его производная удовлетворяет в этом шаре условию Липшица с постоянной I. Пусть, далее, в S < -( jco ) оператор F'(х) непрерывно обратим и существует постоянная тп ~> 0 такая, что в Sr(xо)

|[F'Wr1< m. (1)

Пусть, далее, || F (*o) |К ч. Тогда, если q = -^m2lr\< 1 и

r' = m4+fq2k-\< r, (2)

А-О

то уравнение F(x) = 0 имеет решение х* е Sr' (*о)> к которому сходится итерационный процесс Ньютона, начатый с х0. Ско­рость сходимости хп к х* дается неравенством

2П-1

II хя — х* ||< mil —-у-.

1 — Q

Доказательство. Введем для краткости следующие обозначения: '

r(x) = [F'(x)r\ Г„ = Г(4 Fn = F(xn), Fi = F'(xn).

Итерационный процесс Ньютона в этих обозначениях записы­вается так:

xn+i = x„ — YnFn, п = 0, 1, 2, ... (3)

Покажем сначала, что {х„} с: Sr> {ха), xi — xq = TqFq и, сле- довательно, || хх х0 IK тц. Из (2) следует, что xi< =Sr'(x0). Далее, F0 F" 0{ х0) — 0, поэтому

F^Fi — Fo — F'0 (jc, -x0) = F (ж, ) - F (лг0) - F' (xQ) (xt - x0). Пользуясь неравенством (3) п. 32.2, получаем

11^1 IK у/II xi — xo IP.

Дальнейшие рассуждения проведем методом полной математи- ческой индукции. Пусть уже доказано, что хп s Sr' (хо) и что справедливы оценки

\\Xn-Xn-iW< mi\q^-1-1, (4)

II fJK^ 1\\Хп-Хя-Х II2. (5)

Покажем, что тогда

lUn+i — *n IK ттде2" -1, (о) откуда хп е Sr' (*о), и что

||/-„+11К|Л1*„+1-лгп||2. (7)

Действительно, из (3), (1) и (4) имеем

II *„+! - *!. II = II № IK Ш || F „ IK у ml || *„ - ||2 <

< {ml (т'п)2 з2" -2 == mr\ тЧг\^ q2" -2 = тщ~\

Формула (6) доказана. Далее, из (3) имеем

Fn + F'n {Хп+1~Хп) = 0.

Это позволяет оценить Fn+i: Fn+i: Fn+i — Fn — F'n (xn+i — xn) =

= F (*«+i) - F (x„) - F' (xn) (xn+i - xn). Следовательно, по неравенству (3) п. 32.2

II II II *„ +! -*„ IP,

и неравенство (7) тоже доказано.

Теперь установим фундаментальность { хп }. Из неравенства треугольника и оценок (6) имеем

II Хп + р Хп II ^ II хп Я/г + р-1 II 4" II — Хп+р_2 II +...

п+р—1

... +||*„+1 — *„||< тт| £ q2 (8)

k=n

Отсюда ||х„ — лг„ || —> 0 при п-* оо равномерно по р, так как

оо

ряд 2 Я2*'1 сходится. Итак, {*„} — фундаментальная, а в силу

полноты X — сходящаяся. Обозначим через х* ее предел. Вслед­ствие замкнутости Srr 0), так как {х„} cz Sr' ( х0 ), тои/е Sr' (*о).

Докажем, что х*— решение уравнения F(x) = 0. Для этого достаточно перейти к пределу при гс-> оо в равенстве (3). По­лучаем х* — х* — Г(х*) F(x*), откуда F(x*)= 0, так как Г(х*) обратим. Переходя в оценке (8) к пределу при р -*■ оо, полу­чим следующую оценку скорости сходимости {хл} к х*\

||х%-х„||< тл Е q2k~l. (9)

Воспользуемся теперь элементарным неравенством 2I+S > > 1 + s, справедливым для s = 0, 1, ... Полагая в (9) k = п + + s, учитывая, что <? е(0, 1), и используя формулу суммы гео­метрической прогрессии, имеем

II х- - х„ || < тц £ Я2" *'-1 < mm2*-1 ][У= •

s=0 s-0 1 1

Теорема полностью доказана.

34.3. Модифицированный метод Ньютона. Рассмотрим видо­изменение, или, как говорят, модификацию, итерационного ме­тода Ньютона:

= xn-[F'(x0)rlF(xn), « = 0, 1, 2,... (I)

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

Теорема. Пусть в шаре Sr(xо) оператор F(x) дифферен­цируем и его производная удовлетворяет в Sr(xo) условию Лип­шица с постоянной I. Пусть F'(xо) непрерывно обратим и II [f (*о)]~М1 ^ т. Пусть, кроме того, ll/^xo) IK л- Тогда, если 2 т2 < 1 и

(2)

то уравнение F(x) = 0 имеет единственное решение х* е Sr> (xq), к которому сходится модифицированный итерационный процесс Ньютона (1), начатый с Хо. Скорость сходимости хп к х* дается неравенством

и. ____ < 1 — Vl — 2 тг1г\)п

Доказательство. Покажем сначала, что оператор Ф(х) = х— [/*" ' (лг0) ]'Z7 (jc) отображает шар Srr (хо) в себя. Дей­ствительно, используя оценку для [F' 0) ]неравенство (3) п. 32.2, неравенство треугольника и оценку для F(xo), получаем

II Ф (х) - хо II = | [F'0)Г' [F (х) - F' (хо) (х - х0)] II <

< т || F (х) - F (хо) - F' 0) (* - х0) + F (х0) II <

< т || F (х) - F (х0) - F' (х0) (х - х0) II + т || F0) ||<

< -j ml || * — х01|2 + тц < -i- mlr'2 + mr) = г',

ибо г' — наименьший корень квадратного уравнения (см. фор­мулу (2))

-J- mlr'2 - г' ■ + тц = 0.

Покажем теперь, что в Sr'{x0) оператор является сжатием. Имеем следующую оценку:

IIФ' (*) II ■ = II / - \F' (хо)Г1 F' (х)|| = J [F' (хо)]" 1 [F' (xq) -F' (x)]fl<

< /и/1| х — Хо IK ftilr' = 1 - V1 - 2m2lx] < 1.

Итак, в Sr' (хо) выполнены условия принципа сжимающих ото­бражений, откуда и вытекают.утверждения теоремы о существо­вании и единственности решения х* в Sr 0) и о сходимости к нему {хл}. Далее, < 7 = 1 — У1 — 2тЧг\, 1 — q — У1 — 2тЧх[, а IIФ (хо) - хо 11 = I [F' (хо)]" 1 F 0)1 < тц.

Согласно оценке (2) п. 33.2 скорости сходимости метода последовательных приближений имеем

IU, Хо и < Г тп,

1 — < 7 V1

и, таким образом, теорема доказана.

34.4. Пример к методу Ньютона. Рассмотрим следующую краевую задачу

— х" + f (t, х) = 0, 0 < / < 1, (1)

*(0) = ао, д: (1) = а1. (2)

Предположим для простоты, что функция f{t, x) непрерывна вместе со своей частной производной fx(t, x) в полосе 1],

л: е(—со, +оо). Задачу (1) — (2) представим в операторном виде. Для этого введем сначала банаховы пространства X = = С2[а, Ь] и У = С[а, Ь] 4- Е2. Норму в X и в У зададим фор­мулами

|] х ILr = max | л: (о I + max | х" (/) |, 10. 11 10, 1]

|| у Ну = max \h (t) | + [0, 11 v

где у = {h(t); ао, a, }< = У.

Пусть далее F(x) —нелинейный оператор, действующий из X в У по формуле

F (*) = {- (t) + f(t, x (0); * (0) - Оо, X (1) - a , }. (3)

Теперь задача (1) — (2) записывается в виде операторного урав­нения F{x) = 0. Далее нетрудно убедиться, что уравнение F'(x) -z = у, где F'(x) — производная Фреше оператора F в точ­ке х, представляет собой следующую линейную краевую за­дачу:

~z" + fx{t, x{t))z = h(t), z (0) = cto, z(\) = a,. (4)

Следовательно, выражение F'(x)z задается формулой

F' (х) z = {- z" (t) + fx {t, x (/)) z (/); z (0), z (1)} (5)

Но тогда

. F' (u) z-F'(v)z = {[fx (i >, и (/)) - fx (t, v (/))] z (0; 0, 0} и, следовательно,

|| F' (u) - F' (v) ||= тах|*/Л (t, и (/)) - fx (t, v < /)) |. (6)

io, H

Перейдем к рассмотрению итерационного процесса Ньютона. Пусть мы имеем достаточно хорошее начальное приближение

*oW решения задачи (1) —(2), т. е. ||F(x0)IK т], где т|-=достав точно мало. Согласно формуле (3) это означает, что

шах | (() + / (I, хо (/)) I + У(*о(0)-яо)2 + (*о(1)-а1)2< л (7) [0. 1]

Рассмотрим в пространстве X шар S/(jt0)', т. е. множество всех тех x(t)^Xr для которых выполняется неравенство

max | лг (0 - х0 (t) | + max | х" (0 — 4 (t) \ < г, [о. и [о, п

и предположим, что для всех и и а из этого шара, выполняется неравенство

|hit, U(t))-fx(t, v(t))\^l\u(t)-v(t)\. (8)

Согласно равенству (5) тем более в Sr(xo) будет выполняться условие Липшица: l|.F'(u)—Р(и)||^/||и — и||. Потребуем, на­конец, чтобы при всех j: eSr(^o} для решений z задачи (4) вы­полнялась априорная оценка

+ (9)

Тем самым, в Sr(xo) гарантирована оценка || [^'(х)]-1!! ^ т и выяснен смысл постоянных /, г| и т, входящих в формулировки теорем п. п. 34.2 и 34.3 о методе Ньютона.

Выпишем теперь формулы последовательных приближений. Здесь удобнее записать формулы итерационного процесса Нью­тона в виде

F' (х„) xn+i = F' (хп) хп — F (хп.), п = 1, 2, ...

Используя формулы (3) и (5), получаем последовательность линейных краевых задач

— -Си + fx (*> xn(t))xn+l = fs(t, xn(i))xn{t)-f(t, xn(t))

Xn+i(0) = Oo. *n+i(l) = ®i. « = 0, 1, 2,...

Заметим, что полугенная линеаризация задачи (вместо не­линейной задачи решается последовательность линейных задач) еще не решает эту задачу практически. Для' решения получив­шихся линейных задач обычно приходится прибегать к разно­стным методам или к методу Галёркина.

В заключение выпишем формулы модифицированного про­цесса Ньютона:

-< +1 + /„('. X0(t))xn+1=fx(t, X0(t))xn (t) — f(t, *„(/))

^/i+i(0) = a0, xn+i (1) = ab « = 0, 1, 2,...

Упражнение. Для краевой задачи примем xo'(t) = 0. Найдите модифицированным итерационным процессом Ньютона xi(t) и x2(t).

Пусть в этой же задаче xq {t) = t. Найдите x{(t)y.

Принцип Шаудера

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

35.1. Теорема Броуэра о неподвижной точке. Начнем со сле­дующего элементарного примера. Пусть f(x)—непрерывная на [0, 1] функция, отображающая [0, 1] в себя, т. е. для всех х е е[0, 1] 1. Покажем, что в этом случае существует

на [0, 1] неподвижная точка х0 отображения т. е. f(xо) = Хо. Действительно, рассмотрим на [0, 1] функцию < р(х) = х — f(x). Имеем < р(0) = — /(0)^0, ф (1) = 1 — f (1) 2> 0. Если f(0) — 0 или /(1)= 1, то неподвижной точкой f будет соответственно точка л- = 0 или jc = 1. Если же /(0)> 0и /(1)< 1, то ф(0)< 0, ф(1)> 0, и мы можем воспользоваться теоремой о промежу­точных значениях (см. [18]), согласно которой найдется точка £ е(0, 1), в которой ф(Е) = 0, т. е. /(£ )=£. Итак, существо­вание неподвижной точки отображения f в рассмотренном слу­чае доказано.

Упражнение 1. Покажите, что функция Дх) = 2jc(I—х) отображает [0, 1] в себя. Найдите ее неподвижные точки на [0, 1].

Оказывается, рассмотренный пример можно обобщить на л-мерный случай.

JI. Броуэр доказал следующую теорему о неподвижной точке.

Теорема 1 (Броуэр). Пусть оператор А отображает еди­ничный шар S={x^En: 1} n-мерного евклидова про­странства Еп в себя. Тогда в S найдется неподвижная точка оператора А.

Доказательство этой теоремы Броуэра обычно использует тонкие топологические соображения (см. [21], [12]), однако в [8] (см. стр. 505—508) приведено ее доказательство, осно­ванное на методах математического анализа в Еп. Все извест­ные нам доказательства теоремы Броуэра довольно сложны, и останавливаться здесь на них мы не будем. Нам понадобится небольшое усиление теоремы Броуэра. Введем сначала несколь­ко понятий выпуклого анализа (см. п. 1.7).

Определение 1. Пусть в банаховом пространстве X за­дано множество М из конечного числа элементов

М = {х(< =Х: г = 1, я}. (1)

п

Множество всевозможных линейных комбинаций £ А< л: г, где

i-l

л

Бее и = называется выпуклой оболочкой мно-

i~l

жества М и обозначается Со (М).

Упражнение 2. Докажите, что множество Со(М), где А1 задается формулой (1), выпукло и замкнуто.

Упражнение 3. Докажите, что если в (1) е D, i = = 1, ..., л, где множество D выпуклое, то Со(М)а D.

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

Упражнение 4. Покажите, что для множества М = = {0, хь х„}, где элементы xi, х„ образуют базис ба­нахова пространства X, множество Со(тИ) является выпуклым телом в X.

п

Указание. Проверьте, что точка х = -^ является

внутренней точкой множества Со(М).

В заключение приведем усиленный вариант теоремы Броуэ- ра (доказательство см. [11], стр. 573).

Теорема 2 (Броуэр). Пусть Q — выпуклое ограниченное тело n-мерного банахова пространства X. Тогда всякий непре­рывный оператор А, отображающий множество Q в себя, имеет в Q неподвижную точку.

35.2. Нелинейные вполне непрерывные операторы и их ап­проксимации. В § 20 мы изучали линейные вполне непрерыв­ные операторы. В этом пункте будут приведены необходимые сведения о нелинейных вполне непрерывных операторах.

Определение. Нелинейный оператор А с областью опре­деления D(A)czX и с областью значений в У (X и Y — бана­ховы пространства) называется вполне непрерывным (на D(A)), если он непрерывен на D(A) и переводит каждое огра­ниченное множество, лежащее в Ь(Л), в компактное в У мно­жество.

Упражнение I. Покажите, что любой непрерывный опе­ратор Л: X —> У, где X и У — конечномерные банаховы про­странства, является вполне непрерывным.

Упражнение 2. Пусть X, У, Z — банаховы пространства. Пусть оператор F: Х-*- У ограниченный и непрерывный, а Л: У -> Z— линейный вполне непрерывный оператор (Л е Gfl(y, Z)). Докажите, что оператор AF — вполне непрерывный.

Введем операторы Шаудера (нелинейные проекторы), ап­проксимирующие нелинейный вполне непрерывный оператор Л. Пусть оператор Л непрерывен на ограниченном множестве


D а X и множество A(D)aY компактно. По критерию Хаус- дорфа (см. п. 19.3) для каждого е> 0 существует множество

Mz = {yt е A (D), i= 1, .... п),

являющееся конечной е-сетью множества A(D) — замыкания множества A(D). Рассмотрим оператор Аг, отображающий D в У и определяемый следующим правилом: для х е D

п

2 f1! (*) У,

Ае(х) = —«-----------. (О

Z ^ < *)

г=1

где pt(x) = 0, если \\Ах — у, \\ > е, и р, (х) = в — \\Ах — у, \\, если ЦАх — уЛ^е. Оператор Аг часто называют е-проектором Шау­дера.

Упражнение 3. Докажите, что оператор Аг непрерывен на D. Воспользуйтесь непрерывностью нормы.

Отметим теперь некоторые свойства оператора Ае. Свойство 1. || А (х) — Ае (х) е для всех х е D, т. е. опе­ратор Ае аппроксимирует оператор Л на Л с точностью е. Действительно,

п п п

2 ц, (х) А (х) £ Цг (х) yt Yj ^г М 1А (*) - У А А[х)-АЛх) = ^-п =

£ (*) £ 1Ч (х) Z ^г М

i =1 i-l

Отсюда вытекает, что

|| Л (а) - А, (*) к ------------,

i=-i

где суммирование в числителе и знаменателе ведется только по тем индексам i, для которых ||Лх — уМ<. е, ибо если НДх — — уЛ ^ е, то fit(x) = 0. Следовательно,

п

Е (*) е

||Л(х)-Ле(*)||< -Ц; --------------- = 8.

ЁМ*>

г = 1

Тем самым свойство 1 доказано.

Замечание. При построении s-проектора Ле можно вме­сто конечной е-сети множества A(D) взять конечную e-сеть лю­бого компактного множества, содержащего A{D).


Упражнение 4. Покажите, что вполне непрерывный опе­ратор A: Dc.X-+Y (X и Y—банаховы пространства) можно аппроксимировать е-проекторами на любом ограниченном мно­жестве D\ cz D.

В заключение пункта сформулируем еще одно свойство е- проектора Шаудера Ае, доказать которое мы предлагаем чита­телю.

Свойство 2. Область значений оператора Аг содержится в Co(Mg), где Ме — конечная е-сеть множества A{D).

35.3. Принцип неподвижной точки Шаудера. С помощью теоремы Броуэра можно доказывать различные теоремы о не­подвижных точках нелинейных операторов в бесконечномерных банаховых пространствах. В этом пункте будет доказан прин­цип неподвижной точки Шаудера.


Поделиться:



Популярное:

  1. I. Когда все это закончится?
  2. I. Руководство к совершению Всенощного бдения.
  3. Iron-CLUB.ru – Всё о бодибилдинге
  4. MS Word. К объединениям ячеек в таблице относится все ниже перечисленное,
  5. NFMC-30 -инновационный коктейль оказывает комплексное интенсивное воздействие на все аспекты старения, запускает ряд биохимических реакций, восстанавливающих кожу.
  6. VI. ВСЕРОССИЙСКИЙ ПАТРИАРХ ТИХОН И ЕГО ВРЕМЯ
  7. Абсолютно неверно; 2 – едва ли это верно; 3 – скорее всего верно; 4 – совершенно верно»
  8. Авраам встал рано утром, оседлал осла своего, взял с собою двоих из отроков своих и Исаака, сына своего; наколол дров для всесожжения, и встав пошел на место, о котором сказал ему Бог. (Быт. 22:3).
  9. Американский пирог: Все в сборе
  10. Базовые утверждения теории транзактного анализа
  11. Без этого справочника обойтись невозможно: все позиции, включаемые в товарно-сопроводительные документы, акты о выполнении работ и т.п. в обязательном порядке должны быть в него внесены.
  12. Бесполезность всех остальных способов приближения к Аллаху


Последнее изменение этой страницы: 2016-04-11; Просмотров: 875; Нарушение авторского права страницы


lektsia.com 2007 - 2024 год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! (0.228 с.)
Главная | Случайная страница | Обратная связь