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

u_course

.pdf
Скачиваний:
18
Добавлен:
27.02.2016
Размер:
605.11 Кб
Скачать

получим k, m > n

1( uk 12 + um 12)

 

 

 

 

 

 

 

 

 

 

 

uk

− um

2

=e

 

 

uk + um

 

=

 

 

 

 

 

2

1

 

 

2 k k k k −

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u +u

 

) + 2(f,

u +u

 

) =

2

(Φ(uk) − 2(f, uk) + Φ(um) − 2(f, um)) − Φ(

 

2

m

2

m

1

 

 

 

 

 

 

 

 

 

 

 

 

k

 

k

 

 

1

(Φ(uk) + Φ(um)) − Φ(

uk+um

)

1

· 2(d + ε) − d = ε.

 

 

 

= 2

2

 

2

 

 

 

 

 

 

 

 

последовательность

uk

 

 

 

 

 

 

 

 

 

 

 

 

n

фундаментальна и, следо-

Доказали, что

 

 

 

 

2 ok=1

 

 

 

 

 

 

ˆ

вательно, последовательность {uk}k=1 также фундаментальна. Так как H -

полное пространство, то существует элемент u

ˆ

 

 

 

 

 

H такой, что klim uk = u

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

→∞

в H.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Покажем, что u есть элемент, реализующий минимум функционала Φ

ˆ

на H.

Замечание 2. Из анализа известно, что если lim vm = v в B (B -

m→∞

банахово пространство), то lim kvmkB = kvkB. Из сильной сходимости vm

m→∞

кv в H следует слабая сходимость vm к v.

Всилу замечания 2

kumk12 → kuk12 при m → ∞.

(14.5)

Ясно, что um → u в L2(Ω) слабо:

(f, um) → (f, u), m → ∞.

(14.6)

В силу (14.5), ( 14.6)

Φ(um) = kumk21 + 2(f, um) → Φ(u) = kuk21 + 2(f, u), m → ∞.

C другой стороны Φ(um) → d, m → ∞. В силу единственности предела числовой последовательности имеет место равенство Φ(u) = d.

Докажем, что элемент, реализующий минимум функционала Φ(u) единственный.

Пусть v, w - два элемента, реализующие минимум функционала Φ(u):

Φ(v) = Φ(w) = d и v 6= w. Рассмотрим последовательность v, w, v, w, . . . .

81

ˆ

Данная последовательность минимизирует функционал Φ на H. По доказанному выше она должна сходится к элементу, реализующему минимум функционала, а по построению она расходится. Следовательно, v = w = u. Теорема 1 доказана.

Далее мы получим необходимое условие минимума функционала. Пусть

ˆ 1 ˆ

H H (Ω) и w - произвольный элемент из H, u - элемент, реализую-

ˆ

щий минимум функционала Φ на H. Рассмотрим элемент u + tw, где t

-действительная переменная, и рассмотрим функцию

F (t) = Φ(u + tw) = ((u + tw, u + tw)) + 2(f, u + tw) =

kuk21 + t2kwk21 + 2t((u, w)) + 2(f, u) + 2t(f, w) =

Φ(u) + 2t[((u, w)) + (f, w)] + t2kwk21].

Функция F (t) принимает минимальное значение при t = 0, так как

F (0) = Φ(u) ≤ F (t) t E1.

Следовательно, F 0(t)|t=0 = 0. Отсюда 0 = 2[((u, w)) + (f, w)] и

ˆ

(14.7)

((u, w)) + (f, w) = 0 w H.

Соотношение (14.7) - необходимое условие, которому удовлетворяет эле-

ˆ

мент, реализующий минимум функционала Φ на H.

ˆ

Cвязь элемента, реализующего минимум функционала Φ на H, с решением краевой задачи для эллиптического уравнения.

Рассмотрим краевую задачу

n

∂xi

k(x)∂xi

+ a(x)u = −f(x),

(14.8)

i=1

X

 

∂u

 

 

 

 

 

u|∂Ω = 0.

(14.9)

Относительно коэффициентов уравнения (14.8) предположим, что

0 < k0 ≤ k(x) ≤ K, ≤ a(x) ≤ A,

(14.10)

f(x) L2(Ω).

(14.11)

82

ˆ

1

(Ω) со скалярным произведением

 

 

 

Положим H =H

 

 

 

 

 

((u, v)) = Z (k(x)rurv + auv) dx,

 

(14.12)

 

 

Ω

 

 

 

 

 

 

 

 

 

 

 

 

порождающим в силу (14.10) норму, эквивалентную исходной норме в H1

(Ω). Рассмотрим функционал

 

 

 

 

 

Φ(u) = kuk12 + 2(f, u) = Z (k(x)|ru|2 + au2) dx + 2 Z

fu dx.

 

 

 

 

Ω

 

Ω

 

 

 

 

 

 

 

 

ˆ

1

(Ω).

Пусть u - элемент, реализующий минимум функционала Φ на H =H

Необходимое условие, которому удовлетворяет u:

 

 

 

Z (k(x)rurw + auw) dx + Z

 

 

 

 

fw dx = 0

w H1 (Ω),

 

(14.13)

Ω

 

Ω

 

 

 

 

 

то есть u - обобщенное решение класса H1 (Ω) задачи (14.8), (14.9) (см. определение обобщенного решения первой краевой задачи (12.6)).

Рассмотрим краевое условие

∂u

= 0 (14.14)

∂n ∂Ω

и задачу 14.8), (14.14). Предположим выполнение условий (14.10) и

0 < a0 ≤ a(x) ≤ A.

(14.15)

Рассмотрим ˆ 1 и скалярное произведение

H = H (Ω)

Z

((u, v)) = (k(x)rurv + auv) dx,

Ω

эквивалентное в силу (14.10), (14.15) исходному скалярному произведению в H1(Ω).

Элемент u, реализующий минимум функционала Φ(u) на H1(Ω), есть обощенное решение задачи (14.8), (14.14) в классе H1(Ω) (в силу (14.7)).

Доказаны:

83

ˆ 1

Теорема 2. Пусть H =H (Ω) и выполняются условия (14.10), (14.11).

ˆ

Тогда элемент u, реализующий минимум функционала Φ(u) на H со скалярным произведением (14.12) есть обобщенное решение задачи (14.8), (14.9).

ˆ 1

Теорема 3. Пусть H = H (Ω) и выполняются условия (14.10), (14.15).

ˆ

Тогда элемент u, реализующий минимум функционала Φ(u) на H есть обобщенное решение задачи (14.8), (14.14) в классе H1(Ω).

15. Метод Ритца

Мы рассмотрим метод Ритца - конструктивный метод нахождения элементов, реализующих минимум функционала

Φ(u) = kuk12 + 2(f, u)

 

(15.1)

ˆ

1

(Ω). В (15.1)функция f L2(Ω),

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

 

(f, u) - скалярное произведение в L2(Ω), kuk1 = ((u, u))

1/2

ˆ

 

- норма в H,

порожденная скалярным произведением ((·, ·)) и эквивалентная исходной

норме в H1(Ω).

k

X

{ }ˆ { | }

Пусть wj j=1 - базис в H и Xk = u u = ciwi , k - мерное про-

i=1

странство, натянутое на первые k элементов базиса {wj}j=1. Ясно, что

{ }ˆ 0

Xk Xk+1. Так как wj j=1 - базис в H, то Xk H и функционал Φ определен на Xk. По теореме 1 п.14. существует элемент vk Xk, реализующий минимум функционала Φ(u) на Xk. Ясно, что

Φ(vk) ≤ Φ(u) u Xk.

Запишем необходимое условие, которому удовлетворяет элемент, реализующий минимум функционала на Xk: ((vk, v)) + (f, v) = 0 v Xk, в частности,

((vk, wj)) + (f, wj) = 0, j = 1, . . . , k.

(15.2)

84

k

Xi

 

Так как vk Xk, то vk = ciwi, тогда (15.2) можно записать в виде

=1

 

k

 

Xi

(15.3)

ci((wi, wj)) + (f, wj) = 0, j = 1, . . . , k.

=1

 

Однозначная разрешимость системы (15.3) доказана нами в п.13, см.

k

 

 

Xi

ciwi определяется единственным

уравнение (13.4). Следовательно, vk =

=1

 

 

образом.

 

 

Определение. Последовательность

 

 

v1, v2, . . . , vk, . . . .

(15.4)

называют последовательностью Ритца.

Покажем, что последовательность Ритца - минимизирующая последо-

ˆ

вательность функционала Φ на H. Так как vk - элемент, реализующий минимум функционала Φ на Xk, и имеет место следующая схема

X1 X2 . . .

Xk

. . .

ˆ

 

H

 

↓ . . .

. . .

 

v1

v2 . . .

vk

. . .

u,

 

то имеют место соотношения

 

 

 

 

Φ(v1) ≥ Φ(v2) ≥ · · · ≥ Φ(vk) ≥ · · · ≥ d = Φ(u).

(15.5)

Соотношения (15.5) имеют место в силу свойства точной нижней грани числового множества: если E и G - два ограниченных снизу числовых множества и E G, то inf E ≥ inf G.

Зафиксируем ε > 0 и найдем число N = N(ε) и элемент vε XN

N(ε)

(vε = X c(iε)wi) такой, что

i=1

ku − vεk1 < ε.

(15.6)

85

Имеют место следующие соотношения

Φ(vε) = Φ(u + vε − u) = ((u + (vε − u), u + (vε − u))) + 2(f, u + (vε − u)) =

kuk21 + kvε − uk21 + 2((u, vε − u)) + 2(f, u) + 2(f, vε − u) =

= Φ(u) + Φ(vε − u) + 2((u, vε − u)) ≤

d + Φ(vε − u) + 2((u, vε − u)) ≤

d + kvε − uk21 + 2kfk · kvε − uk + 2kuk1 · kvε − uk1

≤ d + kvε − uk21 + 2kfk · kvε − uk + 2kuk1 · kvε − uk1

ˆ

1

(Ω) ) ≤

 

≤ ( в силу (15.6) и эквивалентности норм в H и H

 

≤ d + ε2 + ε(2Ckfk + 2kuk1) ≤ d + C1ε.

(15.7)

Здесь C1 = 1 + 2Ckfk + 2kuk1.

Через vN(ε) обозначим элемент, реализующий минимум функционала на

XNε . В силу (15.7)

 

Φ(vN(ε)) ≤ Φ(vε) ≤ d + C1ε.

(15.8)

Из (15.5), (15.8) следует оценка

 

d ≤ Φ(vk) ≤ Φ(vN(ε)) ≤ d + C1ε k ≥ N(ε).

 

Отсюда

 

|d − Φ(vk)| < C1ε k ≥ N(ε).

 

Доказано, что

 

lim Φ(vk) = d

 

k→∞

 

и, следовательно, последовательность Ритца {vk}k=1 минимизирующая.

Доказанное выше сформулируем в виде теоремы.

 

 

ˆ

1

(Ω). Пусть

Теорема 1. Пусть H - подпространство пространства H

 

ˆ

=

в H введено скалярное произведение ((·, ·)), порождающее норму kuk1

86

((u, u))1/2 эквивалентную норме kuk =

Z

(u2 + |ru|2) dx 1/2

и {wi}i=1

ˆ

Ω

 

 

- базис в H. Тогда существует (единственная) последовательность Ритца

ˆ { }

функционала Φ(u) на H по базису wi i=1. Эта последовательность явля-

ˆ

ется минимизирующей для функционала Φ(u) на H.

Из теоремы 1 п.14 и теоремы 1 следует

Теорема 2. При выполнении условий теоремы 1 последовательность Ритца для функционала Φ(u), построенная по базису {wi}i=1, сходится к

ˆ ˆ

элементу u H, реализующему минимум на H функционала Φ(u), сильно в H1(Ω).

16. Параболическое уравнение

Пусть Ω - ограниченная область в En, ∂Ω C1, ST

= (0, T ] × ∂Ω,

T = ST Ω ∂Ω, QT = (0, T ) × Ω.

 

 

 

 

Первая начально-краевая задача

 

Рассмотрим задачу

 

 

 

+ f(x),

 

 

∂u

n

∂u

 

 

j=1

(16.1)

 

 

∂t =

∂xj k(x)

∂xj

 

 

 

 

X

 

 

(16.2)

 

 

 

 

 

u|ST = 0,

 

 

 

 

 

 

u|t=0 = 0.

 

(16.3)

Предположим, что выполняются следующие условия:

 

 

функция k(x) измерима по Лебегу на Ω и

(16.4)

 

 

0 < k0 ≤ k(x) ≤ K,

x Ω,

 

 

 

 

 

 

 

 

f L2(Ω).

 

(16.5)

 

 

 

 

 

 

 

 

Определим пространства H1(QT ), H1ST (QT ).

 

 

def

 

 

 

 

 

H1(QT ) ≡ {u|u, ut, uxi L2(QT ), i = 1, . . . , n}.

 

def

 

 

 

 

 

H1ST

(QT ) ≡ { замыкание в норме H1(QT ) функций класса C(QT )

равных нулю вблизи гиперповерхности ST }.

87

Определение. Функция u H1ST (QT ) называется обобщенным решением задачи (16.1) - (16.3) класса H1(QT ), если выполняются условие (16.3)

и

Z

n

Z

utv dxdt + i=1

QT

XQT

k∂xi ∂xi dxdt = Z

fv dxdt

v H1ST (QT ). (16.6)

 

∂u ∂v

 

 

 

 

 

QT

 

 

Для доказательcтва однозначной разрешимости задачи (16.1) - (16.3) воспользуемся методом Галеркина. Так же как и в случае эллиптических уравнений мы должны построить галеркинские приближения, получить априорные оценки, перейти к пределу, доказать, что предел есть решение

задачи, доказать единственность этого решения.

Пусть {wj(x)}j=1 - базис в H1 (Ω) ортонормированный в L2(Ω).

Галеркинские приближения

m

Определение. Функция um(t, x) = X c(km)(t)wk(x) называется m-ым

k=1

галеркинским приближением решения задачи (16.1) - (16.3), если c(km)(t)

- непрерывно дифференцируемые функции на [0, T ] (c(km)(t) C1[0, T ]) и являются решением следующей задачи:

Z

∂um(t, x)

wj(x) dx + Z

krum(t, x)rwj(x) dx = Z

f(x)wj(x) dx,

(16.7)

∂t

 

Ω

 

 

Ω

 

Ω

 

 

 

 

 

cj(m)(0) = 0,

j = 1, . . . , n.

 

(16.8)

Задача (16.7), (16.8) есть задача Коши для системы обыкновенных дифференциальных уравнений с постоянными коэффициентами. Действитель-

но,

Z

m

(m)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ci

(t)

wi(x)wj(x) dx =

 

 

 

 

 

 

 

 

i=1

∂t

 

 

 

 

 

 

 

 

 

Ω

X

 

 

m

(m)

 

 

Z wi(x)wj(x) dx =

m

(m)

 

 

 

(m)

 

 

 

 

 

= i=1

∂ci

(t)

i=1

∂ci

(t)

δij =

∂cj

(t)

,

 

 

 

∂t

 

 

∂t

 

 

∂t

 

 

 

 

 

X

 

 

 

Ω

X

 

 

 

 

 

 

 

Z

f(x)wj(x) dx = fj,

Ω

88

Z

krum(t, x)rwj(x) dx =

Ω

 

 

 

 

(t) Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

(m)

 

 

n ∂wi(x)

·

∂wj(x)

dx =

m

 

(m)

(t).

 

= i=1 ci

k i=1

∂xi

 

∂xi

i=1

αijci

 

Z

X

 

 

 

 

Ω

 

X

 

 

 

 

X

 

 

 

n

∂wi(x) ∂wj(x)

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

Ω

 

 

 

 

 

 

 

 

 

 

dx), и мы имеем

 

 

 

 

 

 

 

 

 

 

 

· ∂xm

 

 

 

 

 

 

 

(здесь αji =

k m=1 ∂xm

 

 

 

 

 

 

 

 

 

 

cj(m)(t)

m

 

 

 

(m)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Xi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(16.9)

 

dt

=

αijci

 

(t) = fj,

 

 

j = 1, . . . , m.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из теории обыкновенных дифференциальных [17] уравнений следует, что задача (16.9), (16.8) имеет единственное решение класса C1[0, T ] (c(jm)(t)

C1[0, T ], j = 1, . . . m). Априорные оценки

Умножим (16.7) на c(jm)(t) и просуммируем результат по j от 1 до m.

Получим равенство

 

 

 

 

 

krumrum dx = Z

 

 

 

 

 

 

Z

utmum dx + Z

fum dx.

 

 

 

 

 

Ω

Z

 

 

 

 

Ω

Ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть ((w, v)) =

krumrum dx скалярное произведение в H1

(Ω) ( в

 

 

 

Ω

 

 

 

 

1 = ((u, u))1/2 эквивалентной норме

 

 

 

силу (16.4)) с нормой

k

u

k

k

u

=

 

 

1/2

 

 

 

 

 

 

 

kH1(Ω)

Z

|ru|2 dx

 

. Имеют место соотношения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

(m)

 

 

(m)

 

 

1 (m)

(m)

1 (m) 2

1

Xk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ut

(t), u (t)) = 2 ∂t(u (t), u (t)) = 2 ∂tku (t)k = 2 ∂t

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

(m)

 

2

 

(m)

2

 

(m)

 

 

 

 

 

 

 

2

 

ku

 

(t)k

 

+ ku

 

 

(t)k1 = (f, u

 

(t)).

 

 

 

 

∂t

 

 

 

 

 

 

 

 

 

Проинтегрируем последнее соотношение по интервалу (0, t):

2 c(km)(t) ;

2ku(m)(t)k2

2ku(m)(0)k2

t

ku(m)(τ)k12 dτ = Z0

t

+ Z0

(f, u(m(τ)) dτ.

1

 

 

1

 

 

 

89

Вследствие равенства u(m)(0) нулю

t

 

t

ku(m)(t)k2 + 2 Z0

ku(m)(τ)k12 dτ = 2 Z0

(f, u(m(τ)) dτ ≤

≤ (неравенство Шварца) ≤

t

 

t

2 Z0

kfkku(m)(τ)k dτ ≤ 2C Z0

kfkku(m)(τ)k1 dτ ≤

≤ (неравенство Коши с ε : a, b |ab| ≤ 12 (εa2 + 1ε b2) ε > 0) ≤

ε Z0

t

 

dτ + Cε Z0

t

 

(τ)k1 dτ.

kfk

2

ku

(m)

C

 

 

 

2

Положив ε = C1 , получим неравенство

t

 

 

t

 

ku(m)(t)k2 + Z0

ku(m)(τ)k12 dτ ≤ C2

Z0

kfk2 dτ.

(16.10)

Положив в (16.10) t = T и отбросив первый неотрицательный член в левой части последнего неравенства, приходим к неравенству

T

 

T

 

 

 

 

 

Z

ku(m)(τ)k12 dτ =

Z Z

k(x)|ru(m)(t, x)|2 dxdt =

0

 

0 Ω

 

 

!

2

(16.11)

 

=

k

n

∂xi

 

 

dxdt ≤ C1,

 

 

Z

 

∂u(m)

 

 

 

 

 

i=1

 

 

 

 

 

 

QT

X

 

 

 

 

где C1 = C2kfk2T .

Отбрасывая второй член в левой части неравенства(16.10) и интегрируя полученное неравенство по t от 0 до T , получим, что

T

u(m)(t, x)

 

dxdt = ku(m)kL2

2(QT ) ≤ C1T.

(16.12)

Z Z

 

 

 

 

2

 

 

0 Ω

 

 

 

Неравенства (16.11), (16.12) дают равномерную по m оценку в L2(QT ) функций u(m), u(xmi ), i = 1, . . . , n.

90

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]