![](/user_photo/2706_HbeT2.jpg)
u_course
.pdfполучим 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
![](/html/2706/1221/html_7IqABxhUXY.RbI_/htmlconvd-raZv7H83x1.jpg)
ˆ |
◦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
![](/html/2706/1221/html_7IqABxhUXY.RbI_/htmlconvd-raZv7H88x1.jpg)
◦
Определение. Функция 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
![](/html/2706/1221/html_7IqABxhUXY.RbI_/htmlconvd-raZv7H89x1.jpg)
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
![](/html/2706/1221/html_7IqABxhUXY.RbI_/htmlconvd-raZv7H90x1.jpg)
Вследствие равенства 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