close

Вход

Забыли?

вход по аккаунту

?

Наилучшие квадратурные формулы приближенного вычисления криволинейного интеграла первого рода для некоторых классов функций и кривых.

код для вставкиСкачать
Известия Тульского государственного университета
Естественные науки. 2013. Вып. 2. Ч.1. С. 50–57
Математика
УДК 517.5
Наилучшие квадратурные формулы
приближенного вычисления
криволинейного интеграла первого рода
для некоторых классов функций и кривых
К. Тухлиев
Аннотация. Рассматривается
задача
о
приближенном
вычислении криволинейного интеграла первого рода для некоторых
классов функций и классов пространственных кривых.
Ключевые слова: квадратурная формула, криволинейный
интеграл, градиент, погрешность, узлы.
Пусть функция f (M ) = f (x1 , x2 , ..., xm ) определена и интегрируема вдоль
кривой Γ ⊂ Rm и
Z
Z
J (f ; Γ) :=
f (M )ds =
f (x1 , x2 , ..., xm )ds.
(1)
Γ
Γ
Предположим, что на кривой Γ установлено положительное направление
так, что положение точки M = M (x1 , x2 , ..., xm ) на кривой может быть
⌣
определено длиной дуги s = AM , отсчитываемой от начальной точки А.
Тогда кривая Γ параметрически выразится уравнениями
x1 = ϕ1 (s), x2 = ϕ2 (s), ..., xm = ϕm (s), 0 6 s 6 L,
(2)
а функция f (x1 , x2 , ..., xm ), заданная в точках кривой, сведется к
сложной функции f (ϕ1 (s), ϕ2 (s), ..., ϕm (s)) от переменной s. В этом случае
криволинейный интеграл (1) запишется в виде определенного интеграла
Z L
J (f ; Γ) =
f (ϕ1 (s), ϕ2 (s), ..., ϕm (s))ds.
(3)
0
Всякая квадратурная формула вида
J (f ; Γ) ≈ LN (f ; Γ; P, S) :=
N
X
k=1
pk f (ϕ1 (sk ), ϕ2 (sk ), ..., ϕm (sk ))
(4)
Наилучшие квадратурные формулы для вычисления криволинейного интеграла
51
для приближенного вычисления интеграла (3) задается векторами
N
коэффициентов P = {pk }N
k=1 и узлов S = {sk }k=1 (0 6 s1 < s2 < ... < sN 6 L),
где p1 , p2 , ..., pN — произвольные действительные числа. При фиксированном
N > 1 через A будем обозначать множество векторов коэффициентов
и узлов (Р , S), либо некоторое его подмножество, определяемое теми
или иными ограничениями на коэффициенты и узлы формулы (4)
(например, требованием точности формулы на многочленах заданной
степени, положительностью коэффициентов pk , k = 1, 2, ..., N и др.).
Обозначим через NQ (L) — класс пространственных спрямляемых кривых
Γ, у которых длина равна L и кривизна кусочно-непрерывна. В дальнейшем
предположим, что все кривые Γ ∈ NQ (L) расположены в области
Q = {(x1 , x2 , ..., xm ) :
(1)
m
X
i=1
x2i 6 L2 }.
(1)
Обозначим через Wp (K; Q) := Wp Lp (K; Q), 1 6 p 6 ∞ класс функций
f (M ) := f (x1 , x2 , ..., xm ), у которых почти всюду в области Q существуют
частные производные ∂f /∂xi (i = 1, 2, ..., m) и удовлетворяют условию
k▽f (ϕ1 (.), ϕ2 (.), ..., ϕm (.))kp :=
=
µZ
0
L
¶1/q
|grad f (ϕ1 (s), ϕ2 (s), ..., ϕm (s))| ds
6 K, 1 6 p < ∞,
p
essup{|grad f (ϕ1 (s), ϕ2 (s), ..., ϕm (s))| s ∈ [0, L]} 6 K, p = ∞,
где
▽f (ϕ1 (s), ..., ϕm (s)) =
m
X
∂f dϕi
·
.
∂ϕi ds
i=1
Через
(1)
W0,p (K; Q)
обозначим
класс
функций,
состоящий
из
всех
(1)
Wp (K; Q),
функций f ∈
удовлетворяющих дополнительному условию
(1)
f (ϕ1 (0), ϕ2 (0), ..., ϕm (0)) = 0. Для каждой функции f ∈ W0,p (K; Q) и каждой
кривой Γ ∈ NQ (L) остаток квадратурной формулы (4) имеет вполне
определенное значение
|RN (f ; Γ; P, S)| = |J (f ; Γ) − LN (f ; Γ; P, S)|.
За величину, характеризующую точность квадратурной формулы для
(1)
всех функций f ∈ W0,p (K; Q), определенных вдоль кривой Γ ∈ NQ (L),
примем число
³
´
n
o
(1)
(1)
RN W0,p (K; Q); Γ; P, S = sup |RN (f ; Γ; P, S)| : f ∈ W0,p (K; Q) .
К. Тухлиев
52
Наибольшую погрешность квадратурной формулы (4) всего класса функций
(1)
W0,p (K; Q) на классе кривых NQ (L) обозначим
³
´
n
³
´
o
(1)
(1)
RN W0,p (K; Q); NQ (L); P, S = sup RN W0,p (K; Q); Γ; P, S : Γ ∈ NQ (L) .
Для получения формулы, которую можно было бы считать оптимальной
(1)
для всех функций f ∈ W0,p (K; Q) и кривых Γ ∈ NQ (L), полагаем, что
соотношение (4) является точным для констант:
Z
N
X
ds =
pk = L.
Γ
k=1
Нижнюю грань
³
´
(1)
EN W0,p (K; Q); NQ (L) =
n
³
´
o
(1)
= inf RN W0,p (K; Q); NQ (L); P, S : (P, S) ∈ A
по аналогии с монографией С.М.Никольского [1] будем называть
оптимальной оценкой погрешности квадратурной формулы (4) на
(1)
рассматриваемых классах функций W0,p (K; Q) и кривых NQ (L). Если
существует квадратурная формула, для которой
³
´
³
´
(1)
(1)
EN W0,p (K; Q); NQ (L) = RN W0,p (K; Q); NQ (L); P 0 , S 0 ,
(1)
то будем ее называть оптимальной или наилучшей на классах W0,p (K; Q)
0
0 N
и NQ (L), а векторы P 0 = {p0k }N
k=1 и S = {sk }k=1 назовем соответственно
оптимальные коэффициенты и узлы квадратурной формулы.
(1)
Записав для произвольной функции f ∈ W0,p (K; Q) как функции одного
переменного F (s) := f (ϕ1 (s), ϕ2 (s), ..., ϕm (s)) формулу Тейлора с остаточным
членом в интегральной форме Коши, остаток RN (f ; Γ; P, S) квадратурной
формулы представим в виде
Z L
▽f (ϕ1 (s), ϕ2 (s), ..., ϕm (s)) · Φ(s)ds,
(5)
RN (f ; Γ; P, S) =
0
где
Φ(s) = L − s −
N
X
k=1
pk (sk − s)0+ , t0+ = {1, если t > 0; 0, если t 6 0}.
Используя неравенство Гельдера с учетом тождества
¶
m µ
X
dϕi 2
i=1
(1)
ds
≡ 1 для
произвольной функции f ∈ W0,p (K; Q), 1 6 p 6 ∞ и произвольной кривой
Наилучшие квадратурные формулы для вычисления криволинейного интеграла
53
Γ ∈ NQ (L) из (5), применяя неравенство Гельдера, получаем оценку сверху
Z L
|▽f (ϕ1 (s), ϕ2 (s), ..., ϕm (s))| · |Φ(s)| ds =
|RN (f ; Γ; P, S)| 6
0
=
L
Z
| grad f (ϕ1 (s), ϕ2 (s), ..., ϕm (s))| · |Φ(s)|ds 6
0
6
µZ
L
0
¶1/p µZ
| grad f (ϕ1 (s), ϕ2 (s), ..., ϕm (s))| ds
·
p
L
0
¶1/q
|Φ(s)| ds
6
q
¶1/q
L
1
1
q
,
|Φ(s)| ds
+ = 1, 1 6 q 6 ∞.
p
q
0
µZ
6K·
(6)
Рассмотрим кривую Γ∗ ∈ NQ (L), которая задана параметрическими
уравнениями
√
xi = ϕi (s) := s/ m, 0 6 s 6 L, i = 1, 2, ..., m,
и определим функцию f (ϕ1 (s), ..., ϕm (s)) на кривой Γ∗ ∈ NQ (L) равенством
m Z ϕi (s)
X
∗
ψ(t)dt, f ∗ (ϕ1 (0), ϕ2 (0), ..., ϕm (0)) = 0, (7)
f (ϕ1 (s), ..., ϕm (s)) =
o
i=1
где
µZ
K
ψ(t) = √
m
Φ̃(s) = L −
√
L
0
ms =
¶− p1
|Φ(s)| ds
· |Φ̃(t)|q−1 sgnΦ̃(t),
q
m
X
k=1
pk (sk −
√
ms)0+ ,
Φ̃
µ
s
√
m
¶
≡ Φ(s).
(1)
Покажем, что функция f ∗ ∈ W0,p (K; Q). Из (7) имеем:
∗
▽f (ϕ1 (s), ..., ϕm (s)) =
m
X
i=1
=
m
X
ψ
i=1
=K
µZ
0
L
µ
s
√
m
¶
√
1
· √ = m·ψ
m
¶− p1
|Φ(s)| ds
=K
q
µZ
0
L
ψ(ϕi (s)) · ϕ′i (s) =
µ
s
√
m
¶
=
¯ µ
¶¯q−1
¶
µ
¯
¯
s
s
¯
=
sgnΦ̃ √
· ¯¯Φ̃ √
m ¯
m
¶− p1
· |Φ(s)|q−1 sgnΦ(s),
|Φ(s)| ds
q
(1)
а потому, используя определение класса W0,p (K; Q), будем иметь
k▽f ∗ (ϕ1 (s), ..., ϕm (s))kpLp (Q) =
(8)
К. Тухлиев
54
L
= Kp
µZ
= Kp
µZ
L
0
0
L
0
¶−1 Z
·
|Φ(s)|q ds
¶−1 Z
q
·
|Φ(s)| ds
0
L
|Φ(s)|(q−1)p ds =
|Φ(s)|q ds = Kp ,
(1)
Wp (K; Q)
и этим включение f ∗ ∈
доказано.
Используя соотношение (7) и (8) из равенства (5), при любом векторе
коэффициентов и узлов (P, S) получаем
Z L
∗
∗
▽f ∗ (ϕ1 (s), ..., ϕm (s)) · Φ(s)ds =
RN (f ; Γ ; P, S) =
0
=K
µZ
0
L
¶1/q
|Φ(s)| ds
, 1 6 q 6 ∞.
q
Таким образом, правая часть (6) является точной верхней границей
(1)
квадратурной формулы (4) на множествах функций W0,p (K; Q) и кривых
NQ (L) :
¶1/q
µZ L
´
³
(1)
q
RN W0,p (K; Q); NQ (L); P, S = K
|Φ(s)| ds
, 1 6 q 6 ∞.
(9)
0
Нетрудно заметить, что экстремальная пара f ∗ и Γ∗ не единственная.
Полагая σk = sk /L, αk = pk /L, перепишем функцию Φ(s) в следующем
удобном для дальнейшего виде
#
"
N
s 0
s X
αk (σk − )+ := LΦ1 (s/L),
(10)
Φ(s) = L 1 − −
L
L
i=1
причем
N
X
αk = 1, 0 6 σk 6 1, k = 1, 2, ..., N.
(11)
k=1
Подставляя равенство (10) в правую часть (9) и сделав замену переменой
t = s/L, приходим к равенству
µZ 1
¶1/q
³
´
(1)
1+ 1
.
|Φ1 (t)|q dt
RN W0,p (K; Q); NQ (L); P, S = M L q
0
Из последнего равенства и (11) сразу следует , что
³
´
³
´
(1)
(1)
1+ 1
EN W0,p (K; Q); NQ (L) = inf RN W0,p (K; Q); NQ (L); P, S = KL q ·
(P,S)
· inf
(µZ
0
1
|Φ1 (t)|q dt
¶1/q
N
: {σk }N
k=1 , 0 6 σk 6 1; {αk }k=1 ,
N
X
k=1
)
αk = 1 .
Наилучшие квадратурные формулы для вычисления криволинейного интеграла
Пользуясь схемой рассуждения из [3, c.78], имеем:
¯q
¯q
Z 1 ¯¯
Z 1 ¯¯
N
N
¯
¯
X
X
¯
¯
¯
¯
αk (σk − s)0+ ¯ ds =
λk (t − tk )0+ ¯ dt,
¯1 − s −
¯t −
¯
¯
0 ¯
0 ¯
k=1
55
k=1
где положено
λk = αN −k+1 , tk = 1 − σN −k+1 , tk < tk+1 , k = 1, 2, ..., N, tN +1 = 1.
(12)
Таким образом, вопрос сводится к нахождению минимума интеграла
¯q
Z 1 ¯¯
Z 1¯
N
¯q
¯
X
¯
¯
¯
(1)
0¯
λk (t − tk )+ ¯ dt :=
¯t − ρN (t)¯ dt
¯t −
¯
0 ¯
0
(13)
k=1
среди всевозможных систем чисел λk и tk , где 0 6 t1 < t2 < ... < tN 6 1, и
фиксированном N .
Итак, требуется найти наилучшее приближение функции f (t) = t
(1)
ступенчатыми функциями ρN (t) в метрике Lq [0, 1], 1 6 q 6 ∞, т.е. такими
функциями, которые на каждом из интервалов (tk , tk+1 ) принимают
постоянное значение, причем
ª
©
(1)
ρN (t) = 0, если 0 < t < t1 ; ck , если tk < t < tk+1 , k = 1, N ; tN +1 = 1 . (14)
Из равенств (13) и (14) имеем:
Z
Z 1
(1)
q
|t − ρN (t)| dt =
t1
q
t dt +
0
0
k=1
N
=
N Z
X
X
tq+1
1
+
q+1
k=1
Z
tk+1
tk
tk+1
tk
|t − ck |q dt =
|t − ck |q dt.
Обычными средствами анализа из (15) при условии t1 +
находим минимальные значения t = t∗k :
t∗k =
(15)
N
X
k=1
(tk+1 − tk ) = 1
2(k − 1) + 1
, k = 1, 2, ..., N.
2N + 1
Экстремальная ступенчатая функция примет вид:
½
¾
1
2k
(1)
∗
∗
ρN ∗ = 0, если 0 < t <
;
, если tk 6 t 6 tk+1 .
2N + 1 2N + 1
Из (12) следует, что λk = 2/(2N + 1) (k = 1, 2, ..., N ),
αk∗ = 2/(2N + 1), σk∗ = 2k/(2N + 1) k = 1, 2, ..., N.
Таким образом, справедливо следующее утверждение.
К. Тухлиев
56
Теорема. Среди квадратурных формул вида (4) для приближенного
вычисления криволинейного интеграла первого рода на классе функций
(1)
W0,p (K; Q) и классе кривых NQ (L) наилучшей является формула
Z
f (M )ds =
Γ
где
Mk∗
=
µ
ϕ1
N
X
2L
f (Mk∗ ) + RN (f ; Γ),
·
2N + 1
k=1
µ
2kL
2N + 1
¶
, ϕ2
µ
2kL
2N + 1
¶
, ..., ϕm
µ
2kL
2N + 1
¶¶
;
xi = ϕi (s) (i = 1, 2, ..., m) — параметрические уравнения кривой Γ, L —
ее длина. При этом точная оценка погрешности на указанных классах
функций и кривых равна
³
´
(1)
EN W0,p (K; Q); NQ (L) = √
q
1+ 1
KL q
1
1
+ = 1.
, 1 6 q 6 ∞,
p
q
q + 1(2N + 1)
Список литературы
1. Никольский С.М. Квадратурные формулы. М.: Наука, 1988. 256 с.
2. Вакарчук С.Б.
Оптимальная формула численного интегрирования
криволинейного интеграла первого рода для некоторых классов функций и
кривых // Укр. матем. журнал. 1986. Т. 38. № 5. C. 643–645.
3. Шабозов М.Ш., Мирпоччоев Ф.М. Оптимизация приближенного интегрирования
криволинейного интеграла первого рода для некоторых классов функций и
кривых // ДАН РТ. 2010. T. 53. № 6. C. 415–419.
4. Крылов В.И. Приближенное вычисление интегралов. М.: Наука, 1967.
Тухлиев Камаридин (Kamaridin.t54@mail.ru), к.ф.-м.н., доцент,
кафедра информатики и вычислительной математики, Худжандский
государственный университет им. академика Б. Гафурова, Худжанд,
Республика Таджикистан.
The best quadrature formula of approximate calculation of
curvilinear integral of first kind for some classes of functions and
curves
K. Tukhliev
Abstract. In this paper was considered an approximate calculation of curvilinear integral of first kind for some classes of functions and for classes of curvilinear
spaces.
Keywords: quadrature formula, curvilinear integral, gradient, error, node.
Наилучшие квадратурные формулы для вычисления криволинейного интеграла
57
Tukhliev Kamaridin (Kamaridin.t54@mail.ru), candidate of physical and
mathematical sciences, associate professor, department of informatics and computational mathematics, Gafurov Khujand State University, Khujand, Republic
of Tajikistan,.
Поступила 23.03.2013
1/--страниц
Пожаловаться на содержимое документа