ВУЗ:
Составители:
Рубрика:
Наконец, оценка (1.45) следует из (1.40) и (1.54) при ψ = Sx
∗
с учетом
того факта, что E
T
n
(Ix)
2
= E
T
n
(x)
2
, x ∈ L
2
.
Теорема 1.3 и ее следствие полностью доказаны.
Замечание 1.1. Для полной непрерывности оператора
R : L
2
−→ W
1
2
достаточно, напр., чтобы h(s, σ) имела производную
по s, удовлетворяющую при σ − s → 0 условию
∂h(s, σ)
∂s
= O
µ
ln
m
|s − σ|
|σ − s|
α
¶
, 0 6 α < 1, m + 1 ∈ N, m + α > 0.
1.4. Метод коллокации. Пользуясь (1.25), сначала приведем
две вычислительные схемы метода коллокации (кратко: м.к.). Если
приближенное решение ищется в виде полинома (1.21), то соотношения
(1.25) и (2.41) приводят к СЛАУ
n
X
k=−n
{c
k
(g)e
iks
j
+ h
jk
}α
k
= y(s
j
), j = −n, n, (1.58)
где {s
j
} – любая система из 2n + 1 попарно неэквивалентных узлов,
c
k
(g) определены в (1.12), а h
jk
= R(e
ikσ
; s
j
). В дальнейшем в качестве
узлов будем использовать
s
j
= 2jπ/(2n + 1), j = −n, n; (1.59)
это связано, с одной стороны, с получением наиболее простых вычисли-
тельных схем, а с другой – с получением наиболее точных результатов
для рассматриваемых методов.
Для вывода второй вычислительной схемы м.к. воспользуемся
формулами (1.21), (1.23), (1.30):
x
n
(s) =
n
X
k=−n
α
k
e
iks
=
2
2n + 1
n
X
k=−n
x
n
(s
k
) D
n
(s − s
k
). (1.60)
Тогда соотношения (1.25) приводят к СЛАУ
2
2n + 1
n
X
k=−n
x
n
(s
k
)
n
S(D
n
(σ − s
k
); s
j
) +
1
2π
2π
Z
0
h(s
j
, σ)D
n
(σ − s
k
) dσ
o
=
= y(s
j
), j = −n, n. (1.61)
Интегралы S(D
n
(σ − s
k
); s
j
) вычисляются точно, т.к. в силу (1.1),
(1.30), (1.11) и (1.12) имеем (см. также с.45 [69])
S(D
n
(σ − s
k
); s
j
) = −
1
2π
2π
Z
0
ln
¯
¯
¯
sin
σ
2
¯
¯
¯
· D
n
(s
j
− s
k
− σ) dσ =
19
Наконец, оценка (1.45) следует из (1.40) и (1.54) при ψ = Sx∗ с учетом того факта, что EnT (Ix)2 = EnT (x)2 , x ∈ L2 . Теорема 1.3 и ее следствие полностью доказаны. Замечание 1.1. Для полной непрерывности оператора R : L2 −→ W21 достаточно, напр., чтобы h(s, σ) имела производную по s, удовлетворяющую при σ − s → 0 условию µ m ¶ ∂h(s, σ) ln |s − σ| =O , 0 6 α < 1, m + 1 ∈ N, m + α > 0. ∂s |σ − s|α 1.4. Метод коллокации. Пользуясь (1.25), сначала приведем две вычислительные схемы метода коллокации (кратко: м.к.). Если приближенное решение ищется в виде полинома (1.21), то соотношения (1.25) и (2.41) приводят к СЛАУ X n {ck (g)eiksj + hjk }αk = y(sj ), j = −n, n, (1.58) k=−n где {sj } – любая система из 2n + 1 попарно неэквивалентных узлов, ck (g) определены в (1.12), а hjk = R(eikσ ; sj ). В дальнейшем в качестве узлов будем использовать sj = 2jπ/(2n + 1), j = −n, n; (1.59) это связано, с одной стороны, с получением наиболее простых вычисли- тельных схем, а с другой – с получением наиболее точных результатов для рассматриваемых методов. Для вывода второй вычислительной схемы м.к. воспользуемся формулами (1.21), (1.23), (1.30): Xn Xn iks 2 xn (s) = αk e = xn (sk ) Dn (s − sk ). (1.60) 2n + 1 k=−n k=−n Тогда соотношения (1.25) приводят к СЛАУ Xn n Z2π o 2 1 xn (sk ) S(Dn (σ − sk ); sj ) + h(sj , σ)Dn (σ − sk ) dσ = 2n + 1 2π k=−n 0 = y(sj ), j = −n, n. (1.61) Интегралы S(Dn (σ − sk ); sj ) вычисляются точно, т.к. в силу (1.1), (1.30), (1.11) и (1.12) имеем (см. также с.45 [69]) Z2π ¯ ¯ 1 ¯ σ¯ S(Dn (σ − sk ); sj ) = − ln ¯sin ¯ · Dn (sj − sk − σ) dσ = 2π 2 0 19
Страницы
- « первая
- ‹ предыдущая
- …
- 18
- 19
- 20
- 21
- 22
- …
- следующая ›
- последняя »