Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология Образование Политология Производство Психология Стандартизация Технологии |
Тогда справедливы все утверждения теоремы 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 У a£ i+qet, (4) l — l I I » / m m—I m—& Фа + С», I + Z где Cm — биномиальные коэффициенты. При этом, если m^s> то полагаем Вспомним теперь, что | ^ < 7, и, пользуясь неравенством Коши — Буняковского, получим
I m—s II m—s 11 s < 3) где
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
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 = тщ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 т21ц < 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 + т || F (х0) ||< < -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. Принцип неподвижной точки Шаудера. С помощью теоремы Броуэра можно доказывать различные теоремы о неподвижных точках нелинейных операторов в бесконечномерных банаховых пространствах. В этом пункте будет доказан принцип неподвижной точки Шаудера. Популярное:
|
Последнее изменение этой страницы: 2016-04-11; Просмотров: 875; Нарушение авторского права страницы