ВУЗ:
Составители:
Рубрика:
6. Метод Либмана
Если число узлов сетки S
h
велико, то непосредственное решение системы (5.2) становит-
ся затруднительным. Кроме того, для криволинейной области G значения функции и в гра-
ничных узлах сетки S
h
выбраны слишком грубо. Эти обстоятельства заставляют для реше-
ния указанной системы прибегать к итерационным методам с одновременным исправлением
граничных значений.
Согласно процессу усреднения Либмана, выбрав начальные приближения
)0(
ij
u , последова-
тельные приближения
)(k
ij
u
для внутренних узлов (x
i
,y
i
) сетки S
h
определяем по формуле
[]
)1(
1,
)1(
1,
)1(
,1
)1(
,1
)(
4
1
−
+
−
−
−
+
−
−
+++=
k
ji
k
ji
k
ji
k
ji
k
ij
uuuuu
(k=1,2, ...). (6.1)
Рис. 13.
Что касается граничных узлов А
h
сетки S
h
, то значения функции и(А
h
) в этих узлах после-
довательно исправляем по формулам линейной интерполяции:
u
(o)
(A
h
)=u(A)=
ϕ
(A) ,
u
(k)
(A
h
)=u(A)+
δ
δ
+
−
−
u
AuBu
k
)()(
)1(
, (k=1,2, …), (6.2)
где А - ближайшая к A
h
точка границы Г (u(А)==
ϕ
(A)), В - ближайший к A
h
внутренний
узел сетки S
h
(рис. 13) и
δ
- удаление узла A
h
, от точки А, причем
δ
>0, если A
h
- внутрен-
няя точка области G, и
δ
<0, если A
h
- внешняя точка области G. В частном случае, если
узел A
h
лежит на границе Г (A
h
≡
A,
δ
=0), то имеем точно
u
(k)
(A
h
)=u(A)=
ϕ
(A).
На практике после некоторого шага h можно считать u
(k)
(A
h
) неизменными (например,
если эти значения установятся с заданной точностью).
За начальные значения
)(o
ij
u теоретически можно взять любую систему чисел. Однако сле-
дует иметь в виду, что в силу принципа максимума для значений искомой функции и(х, у)
должны быть выполнены неравенства
m
≤
u
ij
≤
M ,
где границе Г: m=min
ϕ
(P) и M=max
ϕ
(P). Поэтому разумно полагать
m
≤
)(o
ij
u
≤
M .
В [1] и [3] доказаны, что для любого шага сетки h метод Либмана независимо от выбо-
ра начальных значений сходится, т. е. существует
ij
k
ij
k
uu =
∞→
)(
lim
, причем погрешность при-
ближенного решения имеет порядок O(h
2
).
6. Метод Либмана Если число узлов сетки Sh велико, то непосредственное решение системы (5.2) становит- ся затруднительным. Кроме того, для криволинейной области G значения функции и в гра- ничных узлах сетки Sh выбраны слишком грубо. Эти обстоятельства заставляют для реше- ния указанной системы прибегать к итерационным методам с одновременным исправлением граничных значений. Согласно процессу усреднения Либмана, выбрав начальные приближения u ij( 0) , последова- тельные приближения u ij(k ) для внутренних узлов (xi,yi) сетки Sh определяем по формуле u ij( k ) = 4 [ 1 ( k −1) u i −1, j + u i(+k1−,1j) + u i(,kj−−11) + u i(,kj−+11) ] (k=1,2, ...). (6.1) Рис. 13. Что касается граничных узлов Аh сетки Sh , то значения функции и(Аh) в этих узлах после- довательно исправляем по формулам линейной интерполяции: u(o)(Ah)=u(A)=ϕ(A) , u ( k −1) ( B ) − u ( A) u(k)(Ah)=u(A)+ δ , (k=1,2, …), (6.2) u +δ где А - ближайшая к Ah точка границы Г (u(А)==ϕ(A)), В - ближайший к Ah внутренний узел сетки Sh (рис. 13) и δ - удаление узла Ah , от точки А, причем δ>0, если Ah - внутрен- няя точка области G, и δ<0, если Ah - внешняя точка области G. В частном случае, если узел Ah лежит на границе Г (Ah ≡A, δ=0), то имеем точно u(k)(Ah)=u(A)=ϕ(A). На практике после некоторого шага h можно считать u(k)(Ah) неизменными (например, если эти значения установятся с заданной точностью). За начальные значения u ij(o ) теоретически можно взять любую систему чисел. Однако сле- дует иметь в виду, что в силу принципа максимума для значений искомой функции и(х, у) должны быть выполнены неравенства m≤ uij ≤M , где границе Г: m=minϕ(P) и M=maxϕ(P). Поэтому разумно полагать m≤ u ij(o ) ≤M . В [1] и [3] доказаны, что для любого шага сетки h метод Либмана независимо от выбо- ра начальных значений сходится, т. е. существует lim u ij( k ) = u ij , причем погрешность при- k →∞ ближенного решения имеет порядок O(h2).
Страницы
- « первая
- ‹ предыдущая
- …
- 12
- 13
- 14
- 15
- 16
- …
- следующая ›
- последняя »