ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ, 2022, том 58, № 5, с.703-716
ЧИСЛЕННЫЕ МЕТОДЫ
УДК 519.63+519.651+517.956.32
ОПТИМАЛЬНЫЕ ВЫЧИСЛИТЕЛЬНЫЕ АГРЕГАТЫ
В ЗАДАЧE ДИСКРЕТИЗАЦИИ РЕШЕНИЙ
УРАВНЕНИЯ КЛЕЙНА-ГОРДОНА
И ИХ ПРЕДЕЛЬНЫЕ ПОГРЕШНОСТИ
© 2022 г. А. Б. Утесов, А. А. Базарханова
Указан конкретный вычислительный агрегат, реализующий точный порядок погрешности,
возникающей при дискретизации решения уравнения Клейна-Гордона вычислительными
агрегатами, построенными по точной числовой информации о начальных условиях, при-
надлежащих многомерным 1-периодическим классам Никольского. Найдена предельная
погрешность указанного оптимального вычислительного агрегата.
DOI: 10.31857/S0374064122050107, EDN: CBUUGO
1. Постановка задачи дискретизации решений уравнения Клейна-Гордона.Пусть
символом u(x, t; f1, f2) обозначено классическое решение уравнения Клейна-Гордона
2u
2u
2u
=
+...+
-u
(0 t < +∞, x = (x1, . . . , xs) Rs)
(1)
∂t2
∂x21
∂x2s
с начальными условиями
∂u(x, 0)
u(x, 0) = f1(x) и
= f2(x).
(2)
∂t
Решение u(x, t; f1, f2) описывает, в частности, свободную релятивистскую (псевдо) скалярную
частицу массы, равной единице (см., например, [1, с. 54]).
Задача дискретизации бесконечного объекта (в нашем случае решения u(x, t; f1, f2), пред-
ставимого кратным рядом (см. приведённую ниже лемму 1)) состоит в его приближении прос-
тым (в некотором смысле) конечным объектом и в указании точности предложенного прибли-
жения.
Впервые задача дискретизации решений была рассмотрена Н.М. Коробовым для уравне-
ния Пуассона [2, с. 185-190]. Далее эта задача изучалась в работах [3] и [4, § 10.5, с. 217] для
уравнения Лапласа и уравнения теплопроводности соответственно. Следует отметить статью
С.А. Смоляка [5], в которой был предложен метод построения оптимальных вычислительных
агрегатов в задачах численного интегрирования, восстановления функций и дискретизации ре-
шений уравнения в частных производных. Впоследствии метод Смоляка и его различные при-
менения стали объектами изучения многих математиков (см., например, [6-8]). Отметим также
работу [9], содержащую ряд результатов, полученных при дискретизации решений классиче-
ских уравнений математической физики, в частности, решений уравнения Клейна-Гордона в
рамках постановки “Компьютерный (вычислительный) поперечник” (подробную информацию
об этой постановке можно найти, например, в [10]).
В данной работе дискретизация решений u(x, t; f1, f2) производится вычислительными аг-
регатами, построенными по точной числовой информации о начальных условиях f1 и f2, и
находится предельная погрешность оптимального вычислительного агрегата.
Следуя работам [11-13], сформулируем задачу дискретизации решений уравнения Клейна-
Гордона вычислительными агрегатами, построенными по точной числовой информации о на-
чальных условиях.
Пусть даны целые числа N1 1 и N2 1, классы F1 и F2 функций, заданных на Ω1
и Ω2 соответственно, нормированное пространство Y функций, заданных на ΩY , и пусть
703
9
704
УТЕСОВ, БАЗАРХАНОВА
N = N1 + N2. Числовая информация l(N1,N2) = (l(1)1(f1),... ,l(N1)1(f1),l(1)2(f2),... ,l(N2)2(f2))
объёма N о начальных условиях f1 ∈ F1 и f2 ∈ F2 снимается с функционалов
l(1)1 : F1 C, ... , l(N1)1 : F1 C и l(1)2 : F2 C, ... , l(N2)2 : F2 C.
Пусть также дана функция ϕN (z1, . . . , zN ; y) : CN × ΩY C - алгоритм переработки
числовой информации l(N1,N2) о начальных условиях f1 ∈ F1 и f2 ∈ F2 такая, что при
всяком фиксированном (z1, . . . , zN ) CN как функция от y принадлежит пространству Y.
Далее, черезN }Y будем обозначать множество функций ϕN , удовлетворяющих всем пе-
речисленным выше условиям, через (l(N1,N2), ϕN ) - вычислительный агрегат дискретизации,
действующий по правилу ϕN (l(1)1(f1), . . . , l(N1)1(f1), l(1)2(f2), . . . , l(N2)2(f2); ·), а через D(N
-
1,N2)
множество всех вычислительных агрегатов (l(N1,N2), ϕN ).
Для заданных F1, F2, Y и DN =N1+N2=N
D(N1,N2) положим
N1,N2=1,2,...
δN (DN ,F1,F2)Y = min
inf
δN ((l(N1,N2)N ),F1,F2)Y ,
N1+N2=N (l(N1,N2)N )∈DN
где
δN ((l(N1,N2)N ),F1,F2)Y =
= sup
∥u(· ; f1, f2) - ϕN (l(1)1(f1), . . . , l(N1)1(f1), l(1)2(f2), . . . , l(N2)2(f2); ·)Y .
(f1,f2)∈F1×F2
Всюду ниже для любого числа A и положительного числа B запись A
≪ B будет
α,β,...
означать существование постоянной C(α, β, . . .) > 0 такой, что |A| C(α, β, . . .)B, а для
положительных чисел A и B одновременное выполнение соотношений A
B и B ≪ A
α,β,...
α,β,...
записывается в виде A ≻≺ B.
α,β,...
В рамках приведённых обозначений и определений задача дискретизации решений урав-
нения Клейна-Гордона вычислительными агрегатами, построенными по точной числовой ин-
формации о начальных условиях, заключается в установлении точного порядка величины
δN (DN ,F1,F2)Y (т.е. в нахождении положительной последовательностиN }N1, удовлетво-
ряющей соотношению δN (DN , F1, F2)Y
≻≺ ψN, здесь в качестве α,β,... берутся параметры
α,β,...
классов F1, F2 и пространства Y ) и в указании такого вычислительного агрегата
(l(N1,N2), ϕN ) ≡ ϕN (l(1)1(f1), . . . ,l(N1)1(f1),l(1)2(f2), . . . ,l(N2)2(f2); ·),
для которого δN ((l(N1,N2), ϕN ), F1, F2)Y
≻≺ ψN (в этом случае говорят, что вычислительный
α,β,...
агрегат (l(N1,N2), ϕN ) реализует точный порядок дискретизации).
2. Определения оптимального вычислительного агрегата и его предельной по-
грешности.
Определение 1. Вычислительный агрегат (l(N1,N2)N), реализующий точный порядок
погрешности дискретизации, называется оптимальным вычислительным агрегатом.
В работе [14] было дано определение предельной погрешности оптимального вычислитель-
ного агрегата ϕN (l(1)(f), . . . ,l(N)(f); ·) в задаче дискретизации решений u(t, x; f) уравнения
теплопроводности
∂u
2u
2u
=
+...+
(0 t < +∞, x = (x1, . . . , xs) Rs)
∂t
∂x21
∂x2s
с начальным условием u(0, x; f) = f(x). Здесь сформулируем определение предельной по-
грешности оптимального вычислительного агрегата (l(N1,N2), ϕN ), строящегося по двум точ-
ным числовым информациям (l(1)1(f1), . . . ,l(N1)1(f1)) и (l(1)2(f2), . . . ,l(N2)2(f2)), полученным от
начальных условий f1 и f2.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
ОПТИМАЛЬНЫЕ ВЫЧИСЛИТЕЛЬНЫЕ АГРЕГАТЫ
705
Определение 2. Предельной погрешностью оптимального вычислительного агрегата
(l(N1,N2), ϕN ) называется пара (εN1 , εN2 ) положительных последовательностейN1 } иN2 }
таких, что, во-первых, для всех N1 N, N2 N имеет место соотношение
N((εN1N2 ),(l(N1,N2)N),F1,F2)Y ≻≺ δN(DN ,F1,F2)Y ,
(3)
где
N((εN1N2),(l(N1,N2)N),F1,F2)Y =
= sup
sup
{∥u(· ; f1, f2) - ϕN (z(1)1, . . . , z(N1)1, z(1)2, . . . , z(N2)2; ·)Y
:
(f1,f2)∈F1×F2 z(1)1,...,z(N1)1
z(1)2,...,z(N2)
2
|z(1)1 -l(1)1(f1)| < εN1 , . . . , |z(N1)1 -l(N1)1| < εN1 , |z(1)2 -l(1)2(f1)| < εN2 , . . . , |z(N2)2 -l(N2)2| < εN2 } =
= sup
sup
∥u(· ; f1, f2) - ϕN (l(1)1(f1) + γ(1)1εN1 , . . .
(f1,f2)∈F1×F2(1)1|1, ..., |γ(N1)1|1
(1)2|1, ..., |γ(N2)2|1
... ,l(N1)1(f1) + γ(N1)1εN1,l(1)2(f2) + γ(1)2εN2,... ,l(N2)2(f2) + γ(N2)2εN2;·)Y ,
во-вторых, справедливо равенство
N((ηN1 εN1N2 εN2),(l(N1,N2)N),F1,F2)Y
lim
= +
(4)
N1→∞
δN (DN ,F1,F2)Y
N2→∞
для всех сколь угодно медленно возрастающих к + положительных последовательностей
N1 } иN2 }.
Соотношение (3) означает, что при вычислении значений оптимального вычислительного
l(τ)
агрегата ϕN (l(1)1(f1), . . . ,l(N1)1(f1),l(1)2(f2), . . . ,l(N2)2(f2); ·) числовые информации
(f1), τ =
1
= 1, N1, иl(ν)2(f2), ν = 1, N2, можно заменить неточными значениями z(τ)1 и z(ν)2 такими,
что |z(τ)1 -l(τ)1(f)| < εN1 (τ
= 1, N1) и |z(ν)2 -l(ν)2(f)| < εN2 (ν = 1, N2) соответственно,
сохранив при этом точный порядок погрешности оптимальной дискретизации. Выполнение
же равенства (4) означает неулучшаемость порядков погрешностей εN1 и εN2 , так как сколь
угодно медленное возрастание к бесконечности величин εN1 и εN2 нарушает точный порядок
погрешности дискретизации.
3. Необходимые обозначения и определения для формулировки полученных ре-
зультатов. В дальнейшем через |E| будем обозначать количество элементов конечного мно-
жества E. Как обычно, [a] - целая часть числа a. Для упрощения записи будем использовать
обозначение
sin(t
4π2(m, m) + 1)
ρm(t) =
,
4π2(m, m) + 1
где m ∈ Zs, t 0.
Пусть s = 1, 2, . . . и r - положительное число. Класс Никольского Hr2(0, 1)s по определе-
нию состоит из всех суммируемых 1-периодических по каждой переменной функций f(x) =
= f(x1,...,xs), тригонометрические коэффициенты Фурье
f (m), m ∈ Zs, которых удовле-
творяют условию
sup
22jr
f (m)|2 1,
j=0,1,2,...
[2j-1]∥m∥<2j
где m = (m1, . . . , ms) Zs и ∥m∥ = max |mj |.
j=1,s
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
706
УТЕСОВ, БАЗАРХАНОВА
Ниже нам понадобится класс Соболева Wr2(0, 1)s, состоящий из всех суммируемых 1-пе-
риодических по каждой переменной функций f(x) = f(x1, . . . , xs), тригонометрические коэф-
фициенты Фурье
f (m), m ∈ Zs, которых удовлетворяют условию
∥f∥W r
(m2r1 + ... + m2rs ) 1,
2
m∈Zs
где m = max{1;|mj|} для каждого j = 1,s.
Нормированное пространство L2,∞ ≡ L2,∞((0, 1)s × [0; +]) определяется как линейное
пространство всех функций g : Rs × [0; +) C таких, что для каждого t 0 функция
gt(x) = g(x,t) как функция аргумента x ∈ Rs является измеримой, 1-периодической по каж-
дой из своих s переменных и удовлетворяет неравенству
∥g∥L2,∞ sup ∥g( · , t)L2(0,1)s < +∞.
t0
Символом Φ(N1,N2) обозначается множество всех вычислительных агрегатов (l(N1,N2), ϕN )
с функционалами l(1)1(f1) =
f1(m(1)1), ..., l(N1)1(f1)
= fˆ1(m(N1)1), l(1)2(f2) =
f2(m(1)2), ...,
l(N2)2(f2)
= fˆ2(m(N2)2) и функцией ϕN
∈ {ϕN }Y , а символом L(N1,N2) - множество всех
вычислительных агрегатов (l(N1,N2), ϕN ) c линейными функционалами l(1)1 : F1 C, . . . ,
l(N1)1 : F1 C, l(1)2 : F2 C, ..., l(N2)2 : F2 C и функцией ϕN ∈ {ϕN }Y .
Конкретизировав в δN (DN , F1, F2)Y множество DN , классы F1 и F2, пространство Y,
получим различные задачи дискретизации решений уравнения Клейна-Гордона, а также за-
дачи нахождения предельной погрешности оптимального вычислительного агрегата.
В статьях [15, 16] при DN = LN , F1 Hr12 (0, 1)s, F2 = H22 (0, 1)s DN = ΦN , F1 =
= Wr12(0, 1)s, F2 =
2
(0, 1)s, где LN =N1+N2=N,
L(N1,N2) и ΦN = N1+N2=N,Φ(N1,N2)
N1,N2=1,2,...
N1,N2=1,2,...
соответственно, были установлены точные порядки погрешности дискретизации, но не бы-
ли указаны конкретные вычислительные агрегаты, построенные по числовым информациям
объёма N и реализующие точные порядки.
В настоящей работе при конкретизации из [15] доказана оптимальность некоторого вычис-
лительного агрегата (l(N1,N2), ϕN ) и найдена его предельная погрешность.
Справедливы следующие утверждения.
Теорема 1. Пусть даны целое положительное число s, положительные числа r1 > 2 +
+s/2, r2 > 1+s/2 и пусть для каждого целого Ki 1, i ∈ {1,2} : Ni ≡ Ni(Ki) = (2Ki +1)s.
Тогда для любого N ≡ N(K1, K2) = N1 + N2 имеют место соотношения
N1,
1
≻≺
δN ((l
≻≺
,
(5)
δN(LN ,Hr12,H22 )L2,∞
N2)N),Hr12,H22 )
s,r1,r2
s,r1,r2
Nmin{r1/s,(r2+1)/s}
N1,
здесь вычислительный агрегат (l
N2)N) состоит из функционаловl(1)1(f1) =
f1(m(1)1),
N1)
N1)
N2)
N2)
l(
l(1)
l(
...,
(f1)
f1(m
),
(f2)
f2(m(1)2), ... ,
(f2)
f2(m
) и функции
1
1
2
2
2
N1)
N2)
ϕN(z(1)1,... ,z
,z(1)2,...,z
;x,t) =
1
2
N1
N2
=
z(τ)1ρ′m
(τ )
(t) exp{2πi(m(τ)1, x)} +
z(ν)2ρ
(t) exp{2πi(m(ν)2, x)},
1
m(ν)
2
τ=1
ν=1
где целые положительные числа
N1
и
N2
такие, что
(
)
1
1
1
1
min
+
=
+
,
Nr1/s
N(r2+1)/s
N1+N2=N Nr1/s1
N(r2+1)/s2
1
2
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
ОПТИМАЛЬНЫЕ ВЫЧИСЛИТЕЛЬНЫЕ АГРЕГАТЫ
707
Ni)
{m(1)i, . . . , m
}, i ∈ {1, 2}, - некоторое упорядочение множества Ai = {m ∈ Zs : ∥m∥
i
Ki}, i ∈ {1,2}.
N1
N2
Далее, для краткости, символы
и
заменим соответственно на N1 и N2.
Теорема 2. Пусть даны целое положительное число s, положительные числа r1 > 2 +
+ s/2, r2 > 1 + s/2 и пусть для каждого целого Ki 1, i ∈ {1,2}, выполняется Ni
≡ Ni(Ki) = (2Ki + 1)s.
Тогда пара (εN1 , εN2 ) с компонентами εN1 = (Nr1/s+1/21)-1 и
(Nr2+12)-1,
если s = 1,
εN2 =
(N(r2+1)/22(ln N2)1/2)-1, если s = 2,
(Nr2/s+1/22)-1,
если s > 2
является предельной погрешностью оптимального вычислительного агрегата (l(N1,N2), ϕN )
из теоремы 1, т.е. имеет место соотношение
(6)
N((εN1N2),(l(N1,N2)N),Hr12,H22 )L2,∞ ≻≺ δN (LN ,H21,H22 )L2,∞ ,
и для любых сколь угодно медленно возрастающих к +∞ положительных последовательно-
стей {ηN1 } и {ηN2 } справедливо равенство
r1
N((ηN1 εN1N2 εN2 ),(l(N1,N2)N),H2
, Hr22 )L2,∞
lim
= +∞.
(7)
K1→∞
2)L2,∞
δN (LN ,Hr12,H2
K2→∞
4. Вспомогательные утверждения.
Лемма 1 (см. [15, § 2, лемма 1]). Пусть даны целое положительное число s, положи-
тельные числа r1 и r2 такие, что r1 > 2 + s/2 и r2 > 1 + s/2. Тогда для каждой пары
(f1, f2) функций f1 ∈ Hr12 (0, 1)s, f2 ∈ H22 (0, 1)s ряд
f1(m)ρ′m(t)
f2(m)ρm(t))exp{2πi(m,x)}
m∈Zs
сходится абсолютно и является решением уравнения Клейна-Гордона (1) с начальными усло-
виями (2).
Лемма 2 (см. [15, § 2, лемма 3]). Пусть даны положительные числа α1, α2, c1 и c2.
Тогда для каждого N = 2, 3, . . . верно соотношение
min
)
≻≺ N-min12}.
(c1N11 + c2
2
N1+N2=N,
α12,c1,c2
N11,N21
Лемма 3. Пусть p = 2, 3, . . . и Bp = {m ∈ Zs : 1 ∥m∥ p}. Тогда справедливо
соотношение
1,
если s = 1,
1
≻≺
ln p, если s = 2,
m21 + ... + m2s
s
m∈BP
ps-2, если s > 2.
Доказательство. Пусть s = 1. Тогда
1
1
=2
m21 + ... + m2
m2
s
m1=1
1
m∈Bp
Поэтому, согласно неравенству
1
1
dx,
m2
x2
m1=1
1
1
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
708
УТЕСОВ, БАЗАРХАНОВА
имеем
1
1
4
2
dx = 2(1 - 1/(p + 1))
(8)
m21 + ... + m2s
x2
3
m∈Bp
1
Так как
1
1
1
1
dx
+
+...+
,
x2
22
32
p2
1
то
1
1
1+
dx
,
x2
m2
m1=1
1
1
откуда получим цепочку неравенств
)
1
1
2 1+
dx
4.
m21 + ... + m2s
x2
m∈Bp
1
Следовательно, в силу (8) при s = 1 имеет место соотношение
1
≻≺
1.
(9)
m21 + ... + m2s
s
m∈Bp
Перейдём к случаю s 2. Рассмотрим множества
Mp = {x ∈ Rs : 1 x1 p, ... , 1 xs p},
Up = {x ∈ Rs : x1 0, ... , xs 0, 1 x21 + ... + x2s sp2},
Tp = {x ∈ Rs : x1 1, ... , xs 1, (x1 - 1)2 + ... + (xs - 1)2 (p - 1)2}.
Поскольку Tp ⊂ Mp ⊂ Up, то
dx1 · · · dx
s
dx1 · · · dxs
dx1 · · · dxs
(10)
x2
+...+x2s
x2
+...+x2s
x2
+...+x2s
1
1
1
Tp
Mp
Up
Перейдём к сферическим (полярным в случае s = 2) координатам в интегралах, рассмат-
риваемых на множествах Tp и Up, и соответственно получим
{
{
dx1 · · · dx
s
ln p, если s = 2,
dx1 · · · dxs
ln p, если s = 2,
≻≺
и
≻≺
x21 + ... + x2s
s
ps-2, если s > 2
x21 + ... + x2s
s
ps-2, если s > 2,
Tp
Up
откуда, согласно (10), вытекает соотношение
{
dx1 · · · dxs
ln p, если s = 2,
≻≺
(11)
x21 + ... + x2s
s
ps-2, если s > 2.
Mp
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
ОПТИМАЛЬНЫЕ ВЫЧИСЛИТЕЛЬНЫЕ АГРЕГАТЫ
709
Так как
dx1 · · · dx
s
1
≻≺
и
(12)
x21 + ... + x2s
s
m2
+...+m2s
m1=1
ms=1
1
Mp
)
(s-1
1
1
1
τ
≻≺
C
s
+
,
m21 + ... + m2s
s
m2
+...+m2
m2
+...+m2s
τ+1
s
1
m∈Bp
τ=1
1mτ+1
m1=1
ms=1
msp
то в силу (9), (11) и (12) имеем
{
1
ln p, если s = 2,
≻≺
m21 + ... + m2s
s
ps-2, если s > 2.
m∈Bp
Лемма доказана.
5. Доказательство теоремы 1. Сначала оценим величину δN ((l(N1,N2), ϕN ), Hr12 , H22 )L2,∞
сверху. Так как
ϕN(l(1)1(f1),... ,l(N1)1(f1),l(1)2(f2),... ,l(N2)2(f2);x,t) =
=
f1(m(τ)1)ρ′m
(τ )
(t) exp{2πi(m(τ)1, x)} +
f2(m(ν)2)ρ
(t) exp{2πi(m(ν)2, x)} =
1
m(ν)2
τ=1
ν=1
=
f1(m)ρ′m(t)exp{2πi(m,x)} +
f2(m)ρm(t)exp{2πi(m,x)},
m∈A1
m∈A2
то в силу леммы 1 запишем равенство
u(x, t; f1, f2) - ϕN (l(1)1(f1), . . . ,l(N1)1(f1),l(1)2(f2), . . . ,l(N2)2(f2); x, t) =
=
f1(m)ρ′m(t)exp{2πi(m,x)} +
f2(m)ρm(t)exp{2πi(m,x)}.
(13)
m∈Zs\A1
m∈Zs\A2
Обозначим через li, i ∈ {1, 2}, некоторое целое число, удовлетворяющее неравенствам
2li
Ki < 2li+1. Тогда для всех m ∈ Zs таких, что ∥m∥ 2li, при fi ∈ Hri2(0,1)s имеет
место неравенство
)1/2
( ∑
1
fi(m)|2
,
i ∈ {1,2}.
(14)
s,r
i Nri/s
∥m∥2li
i
Действительно, учитывая очевидные соотношения N1/s1 ≻≺
K1,N1/s2 ≻≺
K2, для каждого
s
s
i ∈ {1,2} получим
)1/2
)1/2
( ∑
fi(m)|2
=
fi
(m)|2
=
∥m∥2li
τ=li+1 2τ-1∥m∥<2τ
)1/2
=
2-2τri · 22τri
fi(m)|2
τ=li+1
2τ-1∥m∥<2τ
(
{
})1/2(
)1/2
1
sup
22τri
fi(m)|2
2-2τri
τ=li+1,li+2,...
s,ri Nri/s
2τ-1∥m∥<2τ
τ=li+1
i
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
710
УТЕСОВ, БАЗАРХАНОВА
Так как
(
)1/2
f1(m)ρ′m(t)exp{2πi(m,·)}
=
f1(m)|2′m(t)|2
=
L2
m∈Zs\A1
m∈Zs\A1
)1/2
)1/2
( ∑
( ∑
=
f1(m)|2′m(t)|2
f1(m)|2′m(t)|2
∥m∥>K1
∥m∥2l1
и′m(t)| 1 для всех m ∈ Zs, то
)1/2
( ∑
f1(m)ρ′m(t)exp{2πi(m,·)}
f1(m)|2
,
L2
m∈Zs\A1
∥m∥2l1
откуда в силу (14) следует оценка сверху для L2-нормы первой суммы правой части равен-
ства (13):
1
f1(m)ρ′m(t)exp{2πi(m,·)}
(15)
L2
s,r1 Nr1/s
m∈Zs\A1
1
Теперь оценим L2-норму второй суммы правой части (13) сверху. Из равенств Парсеваля
)1/2
( ∑
f2(m)ρm(t)exp{2πi(m,·)}
=
f2(m)|2m(t)|2
L2
m∈Zs\A2
∥m∥>K2
с учётом неравенства 2l2 K2 имеем
)1/2
( ∑
f2(m)ρm(t)exp{2πi(m,·)}
f2(m)|2m(t)|2
(16)
L2
m∈Zs\A2
∥m∥2l2
Поскольку
in(t
4π2(m21 + . . . + m2s) + 1)
1
1
1
m(t)| =s
≤
4π2(m21 + . . . + m2s) + 1
2π
m21 + ... + m2s
∥m∥
2l2
для всех ∥m∥ 2l2 , то
)1/2
)1/2
( ∑
( ∑
1
f2(m)|2m(t)|2
f2(m)|2
(17)
l2
2
∥m∥2l2
∥m∥2l2
Вследствие неравенств (16), (17), (14) и 2l2
N1/s2 получаем
s
1
f2(m)ρm(t)exp{2πi(m,·)}
(18)
L2
s,r2 N(r2+1)/s
m∈Zs\A2
2
При любом t 0 из (13), (15) и (18) следует неравенство
1
1
∥u( · , t; f1, f2) - ϕN (l(1)1(f1), . . . ,l(N1)1(f1),l(1)2(f2), . . . ,l(N2)2(f2); · , t)L2
+
,
s,r1,r2
Nr1/s1
N(r2+1)/s
2
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
ОПТИМАЛЬНЫЕ ВЫЧИСЛИТЕЛЬНЫЕ АГРЕГАТЫ
711
откуда, поскольку переменная t и упорядоченная пара (f1, f2) произвольны, в силу опреде-
ления пространства L2,∞ и леммы 2 заключаем, что
≪ N-min{r1/s,(r2+1)/s}.
(19)
δN((l(N1,N2), ϕN), Hr12, H22 )L2,∞
s,r1,r2
В статье [15, теорема 4] установлено соотношение
≫ N-min{r1/s,(r2+1)/s}.
(20)
δN (LN ,Hr12,H22 )L2,∞
s,r1,r2
Наконец, сопоставив неравенства δN (LN , Hr12 , H22 )L2,∞ δN ((l(N1,N2), ϕN ), H21 , H22 )L2,∞ ,
(19) и (20), получим соотношения (5). Теорема 1 доказана.
6. Доказательство теоремы 2. Доказательство проведём для случая s > 2 (в случаях
s = 1 и s = 2 доказательство теоремы 2 проводится аналогично).
Так как для произвольно заданных чисел γ(1)1, . . . , γ(N1)1, γ(1)2, . . . , γ(N2)2 таких, что
(1)1| 1, . . . ,
(N1)1| 1,
(1)2| 1, . . . ,
(N2)2| 1,
имеет место равенство
ϕN(l(1)1(f1) + γ(1)1εN1,... ,l(N1)1(f1) + γ(N1)1εN1 ,l(1)2(f2) + γ(1)2εN2 ,... ,l(N2)2(f2) + γ(N2)2εN2 ;x,t) =
=
f1(m)ρ′m(t)exp{2πi(m,x)} +
f2(m)ρm(t)exp{2πi(m,x)} +
m∈A1
m∈A2
+εN1
γ(τ)1ρ′m(τ)(t)exp{2πi(m(τ)
,x)} + εN2
γ(ν)2ρ
(t) exp{2πi(m(ν)2, x)},
1
m(ν)
1
2
τ=1
ν=1
то в силу леммы 1 справедливо выражение
u(x, t; f1, f2) - ϕN (l(1)1(f1) + γ(1)1εN1 , . . . ,l(N1)1(f1) + γ(N1)1εN1 ,l(1)2(f2) +
+ γ(1)2εN2,... ,l(N2)2(f2) + γ(N2)2εN2;x,t) =
=
f1(m)ρ′m(t)exp{2πi(m,x)} +
f2(m)ρm(t)exp{2πi(m,x)} +
Zs\A1
Zs\A2
+εN1
((τ)1)ρ′m(τ) (t) exp{2πi(m(τ)1, x)} +
1
τ=1
+εN2
((ν)2)ρ
(t) exp{2πi(m(ν)2, x)}.
(21)
m(ν)2
ν=1
Оценим L2-норму третьей суммы правой части (21) сверху. Так как
)1/2
(N1
N1
((τ)1)ρ′m(τ)
(t) exp{2πi(m(τ,·)1)}
=εN1
(τ)1|2′m(τ)
(t)|2
,
ε
1
1
τ=1
L2
τ=1
1)1/2 =
√N1 получим
=1
1
N1
((τ)1)ρ′m(τ)
(t) exp{2πi(m(τ)1, · )}
(22)
ε
1
L2
Nr1/s
τ=1
1
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
712
УТЕСОВ, БАЗАРХАНОВА
Оценим L2-норму четвертой суммы правой части (21) сверху:
)1/2
(N2
((ν)2)ρ
(t) exp{2πi(m(τ)2, · )}
=εN2
(ν)2|2
(t)|2
ε
N2
m(ν)2
m(ν)
2
ν=1
L2
ν=1
)1/2
(N2
εN2
(t)|2
(23)
m(ν)
2
ν=1
Теперь используем лемму 3:
1
(t)|2 =
m(t)|2 =0(t)|2 +
m(t)|2 1 +
N(s-2)/s2.
2
m(ν)
m21 + ... + m2s
s
ν=1
m∈A2
m∈A2\{0}
m∈BK2
Следовательно, продолжив неравенство (23), будем иметь оценку сверху для L2-нормы
четвертой суммы правой части (21):
1
1
N2
((ν)2)ρ
(t) exp{2πi(m(ν)2, · )}
N(s-2)/(2s)2
(24)
ε
m(ν)2
L2
s,r2 Nr2/s+1/22
s,r2 N(r2+1)/s
ν=1
2
При фиксированном t 0 из (15), (18), (21), (22) и (24) следует неравенство
∥u( · , t; f1, f2) - ϕN (l(1)1(f1) + γ(1)1εN1 , . . . ,l(N1)1(f1) + γ(N1)1εN1 ,l(1)2(f2) + γ(1)2εN2 ,
1
1
... ,l(N2)2(f2) + γ(N2)2εN2;·,t)L2
+
,
s,r1,r2
Nr1/s1
N(r2+1)/s
2
откуда в силу определения пространства L2,∞ и леммы 2 имеем
≪ N-min{r1/s,(r2+1)/s}.
(25)
N((εN1N2 ),(l(N1,N2)N),Hr12,H22 )L2,∞
s,r1,r2
Поскольку δN (LN , Hr12 , H22 )L2,∞ δN ((l(N1,N2), ϕN ), H21 , H22 )L2,∞
N((εN1N2),(l(N1,N2)N),Hr12,H22)L2,∞,
то вследствие неравенств (20) и (25) получим соотношение (6).
Далее убедимся в справедливости (7). Для каждого Kj N, j ∈ {1, 2}, определим мно-
жество
GKj = {m ∈ Zs : [
β-1/rj ]},K
j
β-1/rjKj ] ∥m∥ 5 · [
j
j
где βKj = minNj , ln Nj }, j ∈ {1, 2}.
Так как lim
βKj = +∞, то существует номер K(0)j такой, что для всех Kj K(0)j
Kj+
имеет место неравенство
βKj 1, j ∈ {1,2}.
(26)
При каждом Kj K(0)j для функций
-rj/s-1/2
βKj Nj
hKj (x) =
exp{2πi(m, x)}
s · 11s · 5rj
m∈GKj
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
ОПТИМАЛЬНЫЕ ВЫЧИСЛИТЕЛЬНЫЕ АГРЕГАТЫ
713
из проверяемого неравенства |GKj | 11sNjβ-s/rj и условий (26) получим цепочку неравенствK
j
∥hKjW rj =
β2K
N-2rj/s-1 m1rj +...+msrjj
j
2
s · 11s · 52rj
m∈GKj
1
β2K
N-2rj/s-1
j
j
s(5N1/sjβ-1/rjKj )2rj |GKj | βKs/rj 1.j
s · 11s · 52rj
С учётом |GKj | Njβ-s/rjKj оценим L2-норму функции hKj , j ∈ {1,2} снизу:
βKj
βKj
1
∥hKjL2
|GKj |1/2
Njβ-s/(2rj )K
β(2rj-s)/(2rj).
(27)
j
Kj
s,rj Nrj/sj
Nj
s,rj Nrj/s
Nj
s,rj Nrj/s
Теперь рассмотрим N -мерные векторы
(α(1)1, . . . , α(N1)1, α(1)2, . . . , α(N2)2),
(α(1)1,... , α(N1)1, α(1)2,..., α(N2)2),
(ω(1)1, . . . , ω(N1)1, ω(1)2, . . . , ω(N2)2),
(ω(1)1,... , ω(N1)1, ω(1)2,... , ω(N2)2)
с компонентами, определёнными равенствами:
(τ)
ĥK
1
(m
1
1) α(τ)1 = -
), τ = 1,N1, и α(ν)2 = 0, ν = 1,N2;
εN1ηN1
(τ)
1
2)
α(τ)1 = -(K1 )(m
), τ = 1,N1, и
α(ν)2 = 0, ν = 1,N2;
εN1 ηN1
(ν)
ĥK
(m
2
2
3) ω(τ)1 = 0, τ = 1, N1, и ω(ν)2 = -
), ν = 1,N2;
εN2ηN2
(ν)
2
4)
ω(τ)1 = 0, τ = 1,N1, и
ω(ν)2 = -(K2 )(m
), ν = 1,N2.
εN2 ηN2
Так как все компоненты указанных выше векторов по модулю не превосходят единицы и
при всех значениях τ = 1, N1, ν = 1, N2 имеют место равенства
ĥK
1
(m(τ)1) + α(τ)1εN1 ηN1 = 0,
(K1 )(m(τ)1) + α(τ)1εN1 ηN1 = 0,
ĥK
(m(ν)2) + ω(ν)2εN2 ηN2 = 0,
(K2 )(m(ν)2) + ω(ν)2εN2 ηN2 = 0,
2
то для всякого вычислительного агрегата (l(N1,N2), ϕN ) Φ(N1,N2), в силу известного вложения
Wr2(0,1)s ⊂ Hr2(0,1)s (см., например, [17, c. 230]), имеем
sup
sup
sup∥u(· ; t; f1, f2) - ϕN
f1(m(1)1) + γ(1)1εN1 ηN1,...
(f1,f2)∈Hr12 ×Hr2
t0
2
(1)1|1, ..., |γ(N1)1|1
(1)2|1, ..., |γ(N2)2|1
...,
f1(m(N1)1) + γ(N1)1εN1ηN1
f2(m(1)2) + γ(1)2εN2ηN2 ,...
f2(m(N2)2) + γ(N2)2εN2ηN2 ;·,t)L2
max{∥u(·,0;hK1 ,0) - ϕN (ĥK1 (m(1)1) + α(1)1εN1 ηN1,
...,ĥK1(m(N1)1) + αN11εN1,0,0,... ,0 ;·,0)L2,∥u(·,0;(-hK1),0) - ϕN((K1)(m11)) +
N2
+ α(1)1εN1ηN1,... ,(K1)(m(N1)1) + α(N1)1εN1ηN1,0,0,... ,0;·,0)L2} =
N2
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
714
УТЕСОВ, БАЗАРХАНОВА
= max{∥u(·,0;hK1 ,0) - ϕN(0,... ,0;·,0)L2 ,∥u(·,0;(-hK1 ),0) - ϕN (0,... ,0;·,0)L2 }
(∑
)1/2
∥u( · , 0; hK1 , 0)L2
K1(m)|2′m(0)|2
= ∥hK1L2 .
m∈HK1
Следовательно, в силу (27) получаем
sup
sup
∥u(· ; f1, f2) -
(f1,f2)∈Hr12 ×Hr2
2
(1)1|1, ..., |γ(N1)1|1
(1)2|1,...,|γ(N2)2|1
N
f1(m(1)1) + γ(1)1εN1ηN1 ,...
f1(m(N1)1) + γ(N1)1εN1 ηN1
f2(m(1)2) + γ(1)2εN2ηN2 ,...
1
...,
f2(m(N2)2) + γ(N2)2εN2 ηN2;·)L2,∞
β(2r1-s)/(2r1).
(28)
K1
s,r1 Nr1/s
1
). Тогда для всех m ∈ GK2 имеют место неравенства
Пусть t0 = β1/r2K2 /(15π√s
2
1
4π2(m21 + . . . + m2s) + 1 1,
15√st0
откуда следует
)
( √
1
sin t0
4π2(m21 + . . . + m2s) + 1
sin
(29)
15√s,m∈GK2.
В силу равенства Парсеваля с учётом (27) и (29) имеем
(∑
)1/2
∥u( · , t0; 0, hK2 )L2 =
K2 (m)|2m(t0)|2
=
m∈GK2
(∑
)1/2
4π2(m21 + . . . + m2s) + 1)
=
K2 (m)|2 sin2(t0
4π2(m21 + . . . + m2s) + 1
s,r2
m∈GK2
1/r2
βK
2
β(2r2-s)/(2r2)+1/r2K2
∥hK2L2
(30)
s,r2
s,r2 N1/s2
N(r2+1)/s
2
Поскольку
sup
sup
sup∥u( · , t; f1, f2) - ϕN
f1(m(1)1) + γ(1)1εN1 ηN1 ,...
(f1,f2)∈Hr12 ×Hr2
t0
2
(1)1|1, ..., |γ(N1)1|1
(1)2|1, ..., |γ(N2)2|1
...,
f1(m(N1)1) + γ(N1)1εN1ηN1
f2(m(1)2) + γ(1)2εN2ηN2 ,...
f2(m(N2)2) + γ(N2)2εN2ηN2 ;·,t)L2
max{∥u(·,t0;0,hK2 ) - ϕN (0,0,... ,0
K2 (m(1)2) + ω(1)2εN2 ηN2 ,...
N1
... ,ĥK2(m(N2)2) + ωN22εN2;·,t0)L2,∥u(·,t0;0,(-hK2)) - ϕN(0,0,... ,0 ,(K2)(m21)) +
N1
+ ω(1)2εN2ηN2,... ,(K2)(m(N2)2) + ω(N2)2εN2ηN2;·,t0)L2} =
= max{∥u(·,t0;0,hK2 ) - ϕN (0,... ,0;·,t0)L2 ,∥u(·,t0;0,(-hK2 )) - ϕN (0,... ,0;·,t0)L2 }
∥u( · , t0; 0, hK2 )L2 ,
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
ОПТИМАЛЬНЫЕ ВЫЧИСЛИТЕЛЬНЫЕ АГРЕГАТЫ
715
то, приняв во внимание (30), получаем
sup
sup
∥u(· ; f1, f2) - ϕN
f1(m(1)1) + γ(1)1εN1 ηN1,...
(f1,f2)∈Hr12 ×Hr2
2
(1)1|1, ..., |γ(N1)1|1
(1)2|1, ..., |γ(N2)2|1
...,
f1(m(N1)1) + γ(N1)1εN1 ηN1
f2(m(1)2) + γ(1)2εN2 ηN2,...
f2(m(N2)2) + γ(N2)2εN2 ηN2 ;·)L2,∞
s,r2
(2r2-s)/(2r2)+1/r2
βK
2
(31)
s,r2
N(r2+1)/s
2
Из неравенств (28) и (31) следует
N((ηN1 εN1N2 εN2 ),(l(N1,N2)N ),Hr12,H22 )L2,∞
s,r1,r2
(32)
(N-r1/s1 + N-(r2+1)/s2) minθ1K1K },
2
s,r1,r2
где θ1 = (2r1 - s)/(2r1), θ2 = (2r2 - s)/(2r2) + 1/r2.
(N-r1/s1 + N-(r2+1)/s2), то с учётом (32) имеем
Поскольку δN (LN , Hr12 , H22 )L2,∞
s,r1,r2
N((ηN1 εN1N2 εN2 ),(l(N1,N2)N ),Hr12,H22 )L2,∞
s,r1,r2
δN (LN,Hr12,H22 )L2,∞ minK11K2 },2
s,r1,r2
отсюда, в силу равенства lim
minθ1K1K2 }=+∞,вытекаетсоотношение(7).Теорема22
K1→∞
K2→∞
доказана.
СПИСОК ЛИТЕРАТУРЫ
1. Владимиров В.С. Уравнения математической физики. М., 1981.
2. Коробов Н.М. Теоретико-числовые методы в приближённом анализе. M., 1963.
3. Жилейкин Я.М. О приближённом решении задачи Дирихле для уравнения Лапласа // Докл. АН
СССР. 1964. Т. 155. № 5. С. 999-1002.
4. Loo Keng Hua, Yang Wang. Application of Number Theory to Numerical Analysis. Berlin; Heidelberg;
New York, 1981.
5. Смоляк С.А. Квадратурные и интерполяционные формулы на тензорных произведениях некоторых
классов функций // Докл. АН СССР. 1963. Т. 148. № 5. С. 1042-1045.
6. Sickel W., Ullrich T. The Smolyak’s algorithm, sampling, on spars grids and function spaces of dominating
mixed smoothness // East J. Approx. 2007. V. 13. № 4. P. 287-425.
7. Naurizbayev N.,Temirgaliyev N. An exact order of discrepancy of the Smolyak grid and some general
conclusions in the theory of numerical integration // Found. Comput. Math. 2012. V. 12. P. 139-172.
8. Кудайбергенов С.С., Сабитова С.Г. О дискретизации решений уравнения Пуассона на классе Ко-
робова // Журн. вычислит. математики. и мат. физики. 2013. Т. 13. № 7. С. 1082-1093.
9. Темиргалиев Н., Таугынбаева Г.Е., Абикенова Ш.К. Дискретизация решений уравнений в частных
производных в контексте Компьютерного (вычислительного) поперечника // Вестн. Евразийск.
нац. ун-та. Сер. Математика. Информатика. Механика. 2019. Т. 126. № 1. С. 8-51.
10. Темиргалиев Н., Жубанышева А.Ж. Компьютерный (вычислительный) поперечник в контексте
общей теории восстановления // Изв. вузов. Математика. 2019. № 1. С. 89-97.
11. Абикенова Ш.К.,Темиргалиев Н. О точном порядке информативной мощности всех возможных ли-
нейных функционалов при дискретизации решений волнового уравнения // Дифференц. уравнения.
2010. Т. 46. № 8. С. 1201-1204.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022
716
УТЕСОВ, БАЗАРХАНОВА
12. Абикенова Ш.К., Утесов А.Б., Темиргалиев Н. О дискретизации решений волнового уравнения с
начальными условиями из обобщённых классов Соболева // Мат. заметки. 2012. Т. 91. № 3. С. 459-
463.
13. Абикенова Ш.К. Дискретизация периодических решений волнового уравнения с начальными усло-
канд. физ.-мат. наук. Астана, 2010.
виями из классов Wr2(0, 1)s, Wωr1 ,...,ωrs2 (0, 1)s и Es : дис
14. Утесов А.Б., Базарханова А.А. Об оптимальной дискретизации решений уравнения теплопроводно-
сти и предельной погрешности оптимального вычислительного агрегата // Дифференц. уравнения.
2021. Т. 57. № 12. С. 1705-1714.
15. Ибатулин И.Ж., Темиргалиев Н. Об информативной мощности всех возможных линейных функ-
ционалов при дискретизации решений уравнения Клейна-Гордона в метрике L2,∞ // Дифференц.
уравнения. 2008. Т. 44. № 4. С. 491-506.
16. Темиргалиев Н., Шерниязов К.Е., Берикханова М.Е. Точные порядки компьютерных (вычисли-
тельных) поперечников в задачах восстановления функций и дискретизации решений уравнения
Клейна-Гордона по коэффициентам Фурье // Совр. проблемы математики. 2013. Вып. 17. Матема-
тика и информатика. 2. К 75-летию со дня рождения А.А. Карацубы. С. 179-207.
17. Никольский С.М. Приближение функций многих переменных и теоремы вложения. М., 1977.
Актюбинский региональный университет
Поступила в редакцию 24.01.2022 г.
им. К. Жубанова, г. Актобе, Казахстан,
После доработки 25.03.2022 г.
Назарбаев университет,
Принята к публикации 21.04.2022 г.
г. Нур-Султан, Казахстан
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№5
2022