Теория функций действительного переменного/Применение принципа сжимающихся отображений

Материал из testwiki
Перейти к навигации Перейти к поиску

Шаблон:Содержание «Теория функций действительного переменного»

В данном разделе приведены примеры применения принципа сжимающих отображений к уравнениям различных видов.

Системы линейных алгебраических уравнений

Рассмотрим отображение A n-мерного арифметического эвклидова пространства в себя

yi=j=1naijxj+bi,i=1,...,n.

Если данное отображение является сжатием, то для решения уравнения Ax=x можно применить принцип сжимающих отображений. Определим, при каких условиях отображение A будет сжимающим. В рассматриваемом пространстве метрику можно ввести несколькими способами. Рассмотрим три варианта.

Равномерная метрика

ρ(x,y)=max1in|xiyi|. Рассмотрим два вектора: x и x' — преобразование A ставит им в соответствие вектора y и y':

ρ(y,y)=max1in|yiy'i|=max1in|j=1naij(xjx'j)|

Применив аксиому треугольника для модуля, получим:

ρ(y,y)max1inj=1n|aij||xjx'j|.

Так как

|xjx'j|max1in|xjx'j|=ρ(x,x),

то

ρ(y,y)(j=1n|aij|)ρ(x,x).

Отсюда следует условие сжимаемости:

j=1n|aij|α<1,i=1,...,n

«Городская метрика»

Пусть теперь метрика задана следующим образом:

ρ(x,y)=i=1n|xiyi|,

тогда

ρ(y,y)=i=1n|yiy'i|=i=1n|j=1naij(xjx'j)|.

Применим необходимое число раз аксиому треугольника:

ρ(y,y)i=1nj=1n|aij||xjx'j|=j=1n(|xjx'j|i=1n|aij|).

Так как

i=1n|aij|max1jni=1n|aij|,

то ρ(y,y)max1jni=1n|aij|j=1n|xjx'j|=(max1jni=1n|aij|)ρ(x,x).

Таким образом, в данном случае, условие сжимаемости имеет вид:

max1jni=1n|aij|α<1.

Евклидова метрика

В данном случае метрика задаётся формулой

ρ(x,y)=i=1n(xiyi)2.

В силу неравенства Коши-Буняковского:

ρ2(y,y)=i=1n(j=1naij(xjx'j))2(i=1nj=1naij2)ρ2(x,x).

Следовательно, в данной метрике, условие сжимаемости принимает вид:

i=1nj=1naij2α<1.

Обыкновенные дифференциальные уравнения и их системы

Задача Коши для одного уравнения

Пусть дано обыкновенное дифференциальное уравнение

dydx=f(x,y)

и задано начальное условие

y(x0)=y0,

причём функция f определена и непрерывна в некоторой плоской области G, содержащей точку (x0,y0), и удовлетворяет в этой области условию Липшица по y:

|f(x,y1)f(x,y2)|M|y1y2|.

Задача решения дифференциального уравнения с начальным условием (задача Коши) эквивалентна задаче решения интегрального уравнения

φ(x)=y0+x0xf(t,φ(t))dt.

Так как функция f непрерывна, то в некоторой плоской области GG,(x0,y0)G будет иметь место

|f(x,y)|K.

Подберём d>0 таким образом, чтобы выполнялись следующие условия:

  1. Если |xx0|d и |yy0|Kd, то (x,y)G;
  2. Md<1.

Обозначим пространство всех непрерывных функций φ, определённых на отрезке |xx0|d и удовлетворяющих условию |φ(x)y0|Kd, как C* и введём в нём метрику следующим образом

ρ(φ1,φ2)=max|xx0|d|φ1(x)φ2(x)|.

Данное метрическое пространство полно, так как является замкнутым подпространством полного метрического пространства (пространства всех непрерывных функций, заданных на том же отрезке). Докажем, что отображение ψ=Aφ, заданное формулой

ψ(x)=y0+x0xf(t,φ(t))dt,

переводит пространство C* в себя и является сжатием этого пространства. Пусть φC* и |xx0|d, тогда

|ψ(x)y0|=|x0xf(t,φ(t))dt|Kd,

а значит A(C*)C*. Теперь докажем, что отображение является сжатием, действительно:

|ψ1ψ2|x0x|f(t,φ1(t))f(t,φ2(t))|dtMdmax|xx0|d|φ1(x)φ2(x)|.

Так как Md<1, то отображение A — является сжатием. Таким образом, задача Коши имеет единственное решение в пространстве C*.

Задача Коши для систем уравнений

Интегральные уравнения второго рода

Уравнения Фредгольма

Уравнение Фредгольма второго рода — это интегральное уравнение вида

x(t)=λabK(t,s)x(s)ds+f(t).

Функции K(t,s) (ядро) и f(t) (правая часть) заданы, λ — произвольный параметр, необходимо найти функцию x.

Будем считать, что функции K(t,s) (ядро) и f(t) является непрерывными при

atb,asb.

Кроме того, существует такое число M, что

|K(t,s)|M.

Рассмотрим отображение y=Ax полного пространства непрерывных на отрезке [a;b] функций C[a;b] в себя, заданное формулой

y(t)=λabK(t,s)x(s)ds+f(t).

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

ρ(x1,x2)=maxt[a;b]|x1(t)x2(t)|.

Рассмотрим две функции x1 и x2 и оценим

ρ(y1,y2)=ρ(Ax1,Ax2)=maxt[a;b]|λabK(t,s)[x1(s)x2(s)]ds|.

По свойствам интеграла:

ρ(y1,y2)|λ|(ba)maxt[a;b]maxs[a;b]|K(t,s)[x1(s)x2(s)]|.

Так как ядро ограничено числом M, то:

ρ(y1,y2)|λ|(ba)Mmaxs[a;b]|x1(s)x2(s)|=|λ|(ba)Mρ(x1,x2).

Таким образом, отображение A является сжимающим при условии

|λ|1(ba)M,

а значит, что при выполнении данного условия последовательность

xn+1(t)=λabK(t,s)xn(s)ds+f(t)

будет сходится к решению уравнения. В качестве начального приближения можно взять произвольную непрерывную функцию x0(t).

Нелинейные интегральные уравнения

Рассмотрим уравнение вида

x(t)=λabK(t,s,x(s))ds+f(t),

где функции K и f непрерывны, а ядро ещё и удовлетворяет условию Липшица по третьему аргументу, то есть выполняется неравенство

|K(t,s,y1)K(t,s,y2)|M|y1y2|.

Рассмотрим отображение A, заданное формулой

Ax=λabK(t,s,x(s))ds+f(t)

Если x1 и x2 — две непрерывные функции и

y1=Ax1,y2=Ax2,

то справедлива оценка

|y1(t)y2(t)||λ|(ba)maxs[a;b]|K(t,s,x1(s))K(t,s,x2(s))||λ|(ba)Mmaxs[a;b]|x1(s)y2(s)|,

таким образом

ρ(Ax1,Ax2)|λ|(ba)Mρ(Ax1,Ax2),

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

|λ|<1(ba)M.

Уравнения Вольтерры

Рассмотрим теперь неоднородное интегральное уравнение типа Вольтерра

x(t)=λtbK(t,s)f(s)ds+f(t).

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

Как и при рассмотрении линейного уравнению Фредгольма, будем считать, что функции K(t,s) (ядро) и f(t) является непрерывными при

atb,asb,

а K(t,s) ещё и ограниченной, то есть существует такое число M, что

|K(t,s)|M.

Рассмотрим отображение A, определённое следующим образом:

y=Ax=λtbK(t,s)f(s)ds+f(t).

Введём обозначение

g0=A(x1x2)=Ax1Ax2

и рассмотрим последовательность

gn+1=Agn.

Рассуждая так же как в случае уравнения Фредгольма, получим оценку:

|g1(t)||λ|M(ta)maxs[a;b]|x1(s)x2(s)|.

Из этой оценки следует

|g2(t)||λ|2M2(ta)22maxs[a;b]|x1(s)x2(s)|,

и вообще

|gn(t)||λ|nMn(ta)nn!maxs[a;b]|x1(s)x2(s)|.

Для любого значения параметра λ можно указать, такое целое число n0, что будет выполняться условие

|λ|n0Mn(ba)n0n0!

и, так как atb, отображение B=An0 будет сжимающим.

Возьмём некоторую непрерывную функцию x0(t) и построим последовательность

{Anx0}.

Рассмотрим подпоследовательность этой последовательности

{Akn0x0},

так как Akn0x0=Bkx0, то данная последовательность будет, в силу принципа сжимающих отображений, сходится к решению уравнения, а следовательно к решению уравнению будет сходится и вся последовательность

{Anx0}.