ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ, 2022, том 58, № 12, с. 1654-1665
УРАВНЕНИЯ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ
УДК 517.95
ПСЕВДОСДВИГ И ФУНДАМЕНТАЛЬНОЕ РЕШЕНИЕ
ΔB-ОПЕРАТОРА КИПРИЯНОВА
© 2022 г. Л. Н. Ляхов, Ю. Н. Булатов,
С. А. Рощупкин, Е. Л. Санина
Изучены решения сингулярных дифференциальных уравнений с оператором Бесселя B-γ
отрицательного порядка -1 < -γ ≤ 0. В связи с этим большой интерес представляют
решения сингулярного дифференциального уравнения Бесселя B-γ J + λ2J = 0, которые
в работе представлены линейно независимыми функциями Jμ и J-μ, μ = (γ + 1)/2.
Рассмотрены некоторые свойства функций Jμ, выраженные через свойства j-функции
Бесселя-Левитана. Введены прямое и обратное Jμ-преобразования Бесселя и определён
оператор T-псевдосдвига, коммутирующий с оператором Бесселя B-γ . Найдено фун-
даментальное решение обыкновенного сингулярного дифференциального оператора B-γ.
Приведено представление фундаментального решения ΔB -оператора Киприянова с осо-
бенностью в точке x = 0 и на конусе |x| = |y| в евклидовом n-мерном полупростран-
стве R+n.
DOI: 10.31857/S037406412212007X, EDN: NCGOAR
1. Оператор Бесселя с отрицательным параметром и ΔB-оператор Киприянова.
Основные результаты. Пусть -γ = -(γ1, . . . , γn) - мультииндекс с отрицательными дроб-
ными координатами -1 < -γi < 0. ΔB -оператором Киприянова называется сингулярный
дифференциальный оператор
∑
∂2
-γi ∂
ΔB =
B-γi , B-γi =
+
,
x ∈ R+n = {x : xi > 0}.
∂x2i
xi ∂xi
i=1
При -1 < -γi ≤ 0 этот оператор естественно назвать оператором Лапласа-Киприянова.
В сферических координатах x = rΘ (Θ ∈ Rn,
|Θ| = 1) он имеет следующее представление
Бельтрами (см. работу [1]):
1
ΔB = Bn-|γ|-1 +
ΔB,Θ, n - |γ| > 0.
(1)
r2
Далее изучается возможность работы с оператором Киприянова с помощью классических
методов математической физики. Отметим, что этот оператор весьма специфичен. Например, с
помощью операторов B-γ формула (1) позволяет произвольно менять размерность евклидова
пространства путём введения скрытых переменных со сферической симметрией. Например, в
области xi > 0 для произвольного m и y ∈ Rm из (1) получим равенство
∑
∂2
f (xi, x′) =
Bγi f(|y|,x′),
∂x2
i
i=1
если |γ| = 0. Тем самым оператор Лапласа в Rn окажется оператором Лапласа-Киприянова
в евклидовом пространстве Rn+m со скрытой сферической симметрией от m переменных.
Более того, оператор Бесселя Bγ с дробным параметром γ = n - 1 + {γ} окажется промежу-
точным между операторами Лапласа в Rn и в Rn+1 в том смысле, что
lim
Bγ = Δn и lim Bγ = Δn+1.
{γ}→0
{γ}→1
1654
ПСЕВДОСДВИГ И ФУНДАМЕНТАЛЬНОЕ РЕШЕНИЕ ΔB -ОПЕРАТОРА
1655
Это позволяет считать сингулярный дифференциальный оператор Бесселя Bγ (с дробным па-
раметром) оператором Лапласа во фрактальной среде со скрытой сферической симметрией∗).
Также принципиально новой по сравнению с фундаментальными решениями, полученными в
работе [3] (см. также монографию [4]), является конструкция фундаментального решения ΔB -
оператора Киприянова: при 0 < n - |γ| < 1 это решение имеет вырождение O(|x|2-n-|γ|) при
x → 0 вместо особенности для ΔB-операторов Лапласа-Бесселя (при n ≥ 2 и положительных
параметрах операторов Бесселя).
Оператор B-γ в “чужой” билинейной форме. В качестве основного пространства
функций рассматриваем подпространство Шварца Sev = Sev[0, ∞), состоящее из быстро убы-
вающих вместе со всеми производными функций, чётных∗∗) по каждой координате аргумента.
Введём весовую билинейную форму
∫∞
(u, v)α = u(x)v(x)xα dx.
0
Соответствующее пространство обобщённых функций будем обозначать S′ev,α. Функцио-
нал Дирака-Киприянова определим равенством (δ-γ , ϕ)α = ϕ(0) для любой функции ϕ ∈ Sev.
Вначале отметим следующее утверждение.
Теорема 1. Пусть γ ∈ [0, 1), α - произвольное действительное число, U и V принад-
лежат основному пространству функций Sev[0, ∞), тогда справедлива формула
(
[
]
)
(γ + α)(α - 1)
(B-γU, V )α = U, Bγ+2α +
V
x2
α
Доказательство легко достигается интегрированием по частям.
Как видим, оператор Bγ+2α + (γ + α)(α - 1)/x2 оказывается сопряжённым к оператору
B-γ в гильбертовом пространстве Lα2(R+1). Конечно, наиболее востребованы случаи α = -γ
и α = 1. Первый представляется наиболее естественным, так как тогда B-γ является само-
сопряжённым оператором в соответствующем скалярном произведении.
В этой работе приводим результаты исследований только для случая α = -γ, 0 < γ < 1.
Это связано не только с традициями, заложенными в научной школе И.А. Киприянова (см. [4]
и имеющуюся в ней библиографию), но и с определением фундаментального решения, которое
жёстко связано с весовой билинейной формой, где определяется это решение.
Перечислим основные результаты.
1. Получены линейно независимые решения сингулярного уравнения Бесселя B-γu(λt) +
+λ2u(λt) = 0 и приведены свойства этих функций: ортогональность, теорема сложения. На ос-
нове одной из этих функций введены взаимно обратные преобразования Бесселя.
2. Введён коммутирующий с оператором B-γ T-псевдосдвиг
)
(√
∫
π
f
x2 + y2 - 2xy cosα
γ+1
Γ((γ + 3)/2)(xy)
Tyxf(x) =
(2)
)γ+1sinγ+1 αdα.
Γ(1/2)Γ((γ + 2)/2)
(√
0
x2 + y2 - 2xy cos α
3. Для γ ∈ [0, 1) фундаментальным решением оператора B-γ в пространстве S′ev,-γ яв-
ляется регулярный S′ev,γ -функционал E(x) = xγ+1/(γ + 1).
4. Фундаментальное решение ΔB -оператора Киприянова при n - |γ| > 0 определено в
теореме 6, а при n - |γ| > 1 - в следствии 2 на основе обобщённого сдвига Пуассона, принад-
лежащего классу обобщённых сдвигов Б.М. Левитана [5, 6].
∗) Ранее этот факт уже отмечен в некоторых прикладных исследованиях (см., например, [2]).
∗∗) Следуя работе [4], функции, определённые на промежутке [0, ∞), называются чётными по xi, если все
её производные нечётного порядка в точке xi равны нулю.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
1656
ЛЯХОВ и др.
Отметим, что конструкция (2) не принадлежит к классу обобщённых сдвигов Левитана,
так как T0xf(x) = f(x) и для TyC = C. Поэтому мы его назвали “псевдосдвигом”. Отметим
также, что конструкции сверток, порождённые двупараметрическим преобразованием Ганке-
ля, изучались В.А. Какичевым и его учениками (см., например, [7-10]). В указанных работах
изучались именно свёртки, полученные на основе одномерной конструкции (2). Будем следо-
вать методу этих работ.
2. Решения сингулярного уравнения Бесселя с отрицательным параметром
-γ ∈ (-1, 0]. Хорошо известное в теории дифференциальных уравнений и в прикладных
задачах математической физики уравнение Бесселя имеет вид вырождающегося уравнения
второго порядка
d2V
dV
t2
+t
+ (t2 - ν2)V = 0, V = V (t),
(3)
dt2
dt
где ν - произвольное число. Линейно независимые решения уравнения (3)
∑
(-1)m
(x)2+ν
Jν(t) =
,
(4)
m!Γ(ν + m + 1)
2
m=0
∑
(-1)m
(t)2m-ν
J-ν(t) =
(5)
m!Γ(m - ν + 1)
2
m=0
называются функциями Бесселя первого рода.
Сингулярное уравнение Бесселя имеет вид
2
d
γ d
BγU + U = 0, Bγ =
+
(6)
dx2
x dx
При γ > 0 решение V классического уравнения Бесселя (3) и решение U сингулярного
уравнения Бесселя (6) в области t > 0 связаны равенством
V (t)
U (t) =
,
γ = 2ν + 1
(ν > -1/2),
tν
которое приводит к следующим линейно независимым решениям уравнения (6):
Γ(1 ± ν)
γ-1
j±ν(t) =
J±ν(t), ν =
,
(t/2)±ν
2
называемым j-функциями Бесселя∗. Для положительного параметра γ оператора Bγ свой-
ства этих функций и их приложения к задачам для сингулярных дифференциальных урав-
нений изучены в работах [5, 4]; на их основе построены взаимно обратные преобразования
Фурье-Бесселя, доказана теорема сложения j-бесселевых функций и другие свойства.
Для решений уравнения (6) с отрицательным параметром -γ ∈ (-1, 0] справедливо сле-
дующее утверждение.
Лемма. Пусть отрицательный параметр оператора Бесселя задан в виде -γ = -2μ + 1
(1/2 < μ < 1). Тогда решение V классического уравнения Бесселя (3) и решение U сингуляр-
ного уравнения Бесселя (6) в области t > 0 связаны равенством U(t) = tμV (t).
Доказательство леммы проводится путём подстановки функции U(t) в уравнение (6) с
оператором B-γ .
Из леммы с очевидностью вытекает
Теорема 2. Пусть -1 < -γ < 0 и μ = (γ + 1)/2. Линейно независимые решения уравне-
ния
2
d
γ d
B-γu(λt) + λ2u(λt) = 0, t > 0, B-γ =
-
,
(7)
dt2
t dt
∗Это название использовал в своих докладах на конференциях в Алма-Ате и в Москве (1975-1980 гг.)
Б.М. Левитан.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
ПСЕВДОСДВИГ И ФУНДАМЕНТАЛЬНОЕ РЕШЕНИЕ ΔB -ОПЕРАТОРА
1657
имеют вид
∑
(-1)m · 2μ
J-μ(t) = tμJ-μ(t) =
(t)2m,
m!Γ(m + 1 - μ)
2
m=0
Jμ(t) = t2μjμ(t) = 2μΓ(1 + μ)xμJμ(t), Jμ(t) = xμJμ(t),
(8)
где Jμ и J-μ - функции Бесселя первого рода (4) и (5) соответственно, а j-μ и jμ - j-
функции Бесселя.
Отметим, что теорема сложения Левитана на основе обобщённого сдвига Пуассона может
применяться к функциям Бесселя, порядок которых строго больше -1/2 (см. [5, § 7, c. 124;
§ 11, с. 137]). Ниже рассмотрим возможности, которые открываются при использовании функ-
ций (8), имеющих особый статус, поскольку, с одной стороны, удовлетворяют сингулярному
уравнению Бесселя с отрицательным параметром -γ (γ > 0), а с другой - определены непо-
средственно через функции Бесселя первого рода с положительным параметром μ = (γ+1)/2.
2.1. Ортогональность Jμ-функций Бесселя. Введём весовую билинейную форму
∫∞
(u, v)-γ = u(x)v(x)x-γ dx,
(9)
0
в рамках которой сингулярный дифференциальный оператор Бесселя B-γ самосопряжён (см.
теорему 1).
Отметим также, что если функция u(x) - решение уравнения (7), то функция u(λx) -
решение уравнения B±γ u(λx) + λ2u(λx) = 0.
Теорема 3. Пусть γ = 2μ - 1 и 0 < μ ≤ 1/2. Имеет место равенство
1
{
∫
0,
k = m,
J∗μ(λkt)J∗μ(λmt)t-γ dt =
(10)
(λkλm)(γ+1)/2M = 0,
k = m,
0
где числа λk являются корнями функций Jμ или корнями функции J′μ.
Действительно, с учётом того, что J∗μ(λkt) = (λkt)(γ+1)/2J(γ+1)/2(λkt), получим следующее
представление левой части равенства (10):
∫
1
∫
1
Jμ(λkt)Jμ(λmt)t-γ dt = (λkλm)(γ+1)/2
Jμ(λkt)Jμ(λmt)tdt.
0
0
Теперь равенство (10) следует из ортогональности функций Бесселя первого рода Jμ поло-
жительных порядков μ (на самом деле для μ > -1, см. [11, с. 265, формулы (12.6) и (13.3)]).
Отметим также, что ряды Фурье по J-функциям Бесселя исследовались, например, в ра-
боте [12].
3. J-преобразование Бесселя. Будем следовать методике В.А. Какичева, опираясь на
теорему сложения Б.М. Левитана для соответствующих j-функций Бесселя [5].
Теорема 4. Пусть γ = 2μ - 1 и 0 < μ ≤ 1/2. Если x-μf(x) ∈ L2(0, ∞), то
∞
∫∞
∫
f (x) = Jμ(xξ)ξ-2μ+1 dξ Jμ(yξ)y-2μ+1f(y) dy,
(11)
0
0
т.е. следующие J-преобразования Бесселя взаимно обратимы:
∞
∫∞
∫
f (ξ) = J(γ+1)/2(yξ)f(y)y-γ dy, F-1
f ](x) = f(x) =
f (ξ)ξ-γ dξ.
FB-γ [f](ξ)
B-γ
J(γ+1)/2(xξ
0
0
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
1658
ЛЯХОВ и др.
Доказательство. Известен следующий интеграл Ханкеля для функций f ∈ L2(0, ∞):
∫∞
∫
∞
f (x) = Jμ(xξ)ξ dξ Jμ(yξ)f(y)y dy.
0
0
Пусть x-μf(x) ∈ L2(0, ∞), μ = (γ + 1)/2 > 0. Тогда выполняются равенства
∫∞
∫
∞
∫
∞
∫
∞
Jμ(xξ)
Jμ(yξ)
f (x)x-μ = Jμ(xξ)ξ dξ Jμ(yξf(y)y-μy dy =
ξ dξ
f (y)y-μy dy =
(xξ)μ
(yξ)μ
0
0
0
0
∫∞
∫
∞
= x-μ Jμ(xξ)ξ-2μ+1 dξ Jμ(yξ)y-2μ+1f(y)dy,
0
0
откуда следует формула (11). Теорема доказана.
4. Обобщённый T-псевдосдвиг и теорема сложения для J-функций Бесселя. Для
j-функций Бесселя известен следующий вариант теоремы сложения [5]: для ν > -1/2
jν(tξ)jν(tτ) = Tτξ jν(ξτ),
где сдвиг Пуассона Tτξ определён равенством
∫
π
)
(√
Γ((γ + 1)/2)
Tτ ϕ(t) =
ϕ
t2 + τ2 - 2tτ cos α sinγ-1 α dα, γ > 0.
(12)
Γ(1/2)Γ(γ/2)
0
Поскольку μ = (γ + 1)/2 > 0, то из (8) имеем
Jμ(tξ)Jμ(tτ) = (tξ)2μjμ(tξ)(tτ)2μjμ(tτ) = (tξ)2μ(tτ)2μ
ξ
jμ(tξ) =
(
)
√
∫π
ξ2 + τ2 - 2ξτ cos α
Jμ t
(tξ)2μ(tτ)2μΓ(μ + 1)
=
(
)2μ sin2μ α dα =
Γ(1/2)Γ(μ + 1/2)
√
0
t
ξ2 + τ2 - 2ξτ cos α
(
)
√
∫
π
Jμ t
ξ2 + τ2 - 2ξτ cos α
(ξτ)2μΓ(μ + 1)
=
)2μ sin2μ αdα = TxJμ(xξ),
2μ = γ + 1.
Γ(1/2)Γ(μ + 1/2)
(√
0
ξ2 + τ2 - 2ξτ cos α
Таким образом, имеет место
Теорема 5 (сложения Jμ-функций Бесселя). Справедлива формула
(13)
J(γ+1)/2(xξ)J(γ+1)/2(yξ) = TxJ(γ+1)/2(xξ), x, y, ξ ∈ R1,
где 0 < γ < 1, μ = (γ + 1)/2 > 0, а T-псевдосдвиг в евклидовом пространстве R+1 имеет
вид (2), где 1/2 ≤ μ < 1.
Отметим, что T-псевдосдвиг не меняет гладкости функции, так как особенность знамена-
теля (возникающая при x = y, α = 0) регулярная.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
ПСЕВДОСДВИГ И ФУНДАМЕНТАЛЬНОЕ РЕШЕНИЕ ΔB -ОПЕРАТОРА
1659
5. Свойства T-псевдосдвига. Весовую билинейную форму в пространстве R+1 опреде-
лим равенством (10).
Свойство 1. Если f и g - функции, суммируемые с весом x-γ, 0 < γ < 1, то
(Tyf, g)-γ = (f, Tyg)-γ .
(14)
Доказательство. Учитывая, что -γ = 1 - 2μ, имеем
)
(√
∞
π
∫
∫
g
x2 + y2 - 2xy cos α
(f, Tyg)γ
=
(xy)2μf(y)(
)2μ sin2μ αdαy1-2μ dy.
C(μ)
√
0
0
x2 + y2 - 2xy cosα
Введём антиполярные координаты
y cos α = z1, y sin α = z2, y dy dα = dz1 dz2 = dz.
Тогда
)
(√
∫
(√
) g (x - z1)2 + z2
(√
2
)-2μ
(f, Tyg)γ
=
x2μf
z21 + z2
z21 + z2
dz.
2
)2μ z2μ
2
2
C(μ)
(√
R+2={z:z2>0}
(z1 - x)2 + z2
2
Теперь сдвиг x - z1 = ξ приведёт к равенству
)
)
(√
(√
∫
f
(x - ξ)2 + z2
g
ξ2 + z2
2
2
(f, Tyg)γ
=x2μ
)2μ (
)2μ z2μ2 dξ dz2.
C(μ)
(√
√
R+2
(x - ξ)2 + z22
ξ2 + z2
2
C помощью формул перехода к полярным координатам
ξ = y cosα, z2 = y sinα, dξ dz2 = dξ dz2 = dαy dy
вернёмся к прежним обозначениям:
)
(√
∫
∞
∫
π
f
x22 + y22 - 2xy cos α g(y)
(f, Tyg)γ
= x2μ
y2μ
sin2μ α dα y-2μ+1 dy =
C(μ)
(x22 + y22 - 2xy cos α)μ
0
0
)
(√
π
∫∞
∫
f
x22 + y22 - 2xy cos α g(y)
(Tyf, g)-γ
= (xy)2μ
sin2μ αy-γ dy =
(x22 + y22 - 2xy cos α)μ
C(μ)
0
0
Отсюда вытекает равенство (15). Свойство доказано.
Свойство 2. Пусть f
- чётная, дважды непрерывно дифференцируемая функция,
x-γf(x) ∈ L2(0,∞), 0 < γ < 1. Тогда
B-γ,xTyxf(x) = TyxB-γ,xf(x).
(15)
Доказательство. Из (11), (13) и (7) имеем
∞
∫
B-γ,xTyxf(x) = B-γ,x Jμ(xξ)Jμ(yξ
f (ξ)ξ-γ dξ =
0
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
1660
ЛЯХОВ и др.
∫∞
∫
∞
=
(-ξ2)Jμ(xξ)Jμ(yξ
f (ξ)ξ-γ dξ =
(-ξ2)TyxJμ(xξ
f (ξ)ξ-γ dξ.
0
0
Здесь множитель (-ξ2) можно внести под знак T-псевдосдвига. В результате получим
∫∞
B-γTyxf(x) = Tyx((-ξ2)Jμ(xξ)
f (ξ)ξ-γ dξ =
0
∫∞
∫
∞
= Tyx(B-γJμ(xξ)
f (ξ)ξ-γ dξ = TyB-γ Jμ(xξ
f (ξ)ξ-γ dξ = TyxB-γf(x).
x
0
0
Свойство доказано.
Свойство 3 (переместительность T-псевдосдвига). Если функция f представлена равно-
мерно сходящимся рядом Фурье по J-функциям Бесселя, то
TyxTzxf(x) = TzxTyxf(x).
Доказательство непосредственно следует из определения T-сдвига и теоремы 5:
TyxTzxJ(λx) = J(λy)J(λz)J(λx) = TzxTyxJ(λx).
Свойство 4 (ассоциативность T-сдвига). Если функция f представлена равномерно схо-
дящимся рядом Фурье по J-функциям Бесселя, то TzyTxf(x) = TzxTxf(x).
Доказательство аналогично доказательству свойства 3.
Свойство 5 (представление многомерного T-псевдосдвига в сферических координатах).
∏n
Если f и g - радиальные в R+n функции, суммируемые с весом x-γ =
i=1
x-γii, μi = (γi +
+ 1)/2, то
∫
∫
∞
μ
s
C(n, μ)
f (|y|),Tx g(|y|)y-γ dy = f(r)
Tρ g(r)rs dr, s = n - |γ| - 1,
Rn
0
μ
где C(n, μ) - константа, не зависящая от подынтегральных функций, через
Tx обозначен
многомерный псевдосдвиг, отвечающий мультииндексу μ = (μ1,
,μn),
1/2 < μi < 1.
Доказательство. В выражении
(√
)
∫
∫
π
g
∏
f (|y|)
y1-2μii sin2μi αi(|x|yi)2μi dαi dy
(. . . + x2
+ y2i - 2xiyi cosαi + ...)μi
i
i=1
Rn
0
по каждой паре переменных (yi, αi) введём антиполярные координаты
yi cos αi = zi,1, yi sin αi = zi,2, yi dyi dαi = dzi,1 dzi,2.
Пусть R+2n = {z : zi,2 > 0}. Тогда справедливо равенство
)
(√
μ
∫
g
... + (xi - zi,1)2 + z2i,2 + ...
(√
)-2μi
∏
(f,Ty
g)-γ
=
f (|z|)
|x|2μi z2μi
z2i,1 + z2
dz.
i,2
i,2
C(μ)
(. . . + (xi - zi,1)2 + z2i,2 + . . .)μi
i=1
R+
2n
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
ПСЕВДОСДВИГ И ФУНДАМЕНТАЛЬНОЕ РЕШЕНИЕ ΔB -ОПЕРАТОРА
1661
Здесь существование интеграла обеспечено выбором функции g и тем, что
)-μi
(z2
i,1
(z2i,1 + z2i,2)-μi =
+1
,2
z2
i,2
Теперь повернём оси координат вокруг гиперплоскости, образованной чётными осями ко-
ординат Ozi,2, так, чтобы орты векторов Oz1,1 и Ox совпали. С учётом того, что радиальные
функции инвариантны относительно вращений системы координат, а точка x в новых коор-
динатах это x = (|x|, 0, . . . , 0), получим
)
(√
μ
∫
g
(|x| - z1,1)2 + z21,2 + . . . z22n,1 + z2
2n,2
∏
(√
)-2μi
(f,Ty
g)-γ
=
f (|z|)
|x|2μ
z2μ
z2i,1+z2
dz =
i,2
i,2
C(μ)
((|x| - z1,1)2 + z21,2 + . . . z22n,2 + z22n,2)μ
i=1
R+
2n
)
(√
∫
g
|x|2 - 2|x|z1,1 + |z|2
(√
)-2μi
∏
=
f (|z|)
|x|2μ
z2μi
z2i,1 + z2
dz.
i,2
i,2
(|x|2 - 2|x|z1,1 + |z|2)μ
i=1
R+
2n
Сферическое преобразование координат
z1,1 = ρcos φ1, z1,2 = ρsin φ1 cos φ2n,
...,
z2n,1 = ρsin φ1 ··· sin φ2n-2 cos φ2n-1, z2n,2 = ρsin φ1 ··· sin φ2n-2 sin φ2n-1
приведёт к равенству
)
(√
μ
∞
∫
∫
(rρ)2|μ|g
r2 + ρ2 - 2rρcos ϕ1
(f,Ty
g)-γ
=
f (ρ)
×
C(μ)
(r2 + ρ2 - 2rρ cos ϕ1)μ
0 S+1(2n)
(√
)-2μi
∏
× Θ2μ
Θ2i,1 + Θ2
dS ρ2(n-|μ|)-1 dρ,
i,2
i,2
i=1
где S+1(2n) = {Θ : |Θ| = 1, Θ2i > 0}. Поскольку на 2n-мерной сфере
dS = sin2n-2 φ dφ1 sin2n-3 φ2 dφ2 . . . sin φ2n-2 dφ2n-2 dφ2n-1,
то
)
(√
μ
∞
π
∫
∫
(rρ)2|μ|g
r2 + ρ2 - 2rρcos α
(f,Ty
g)-γ
=
f (ρ)
sin2|μ| α ×
C(μ)
(r2 + ρ2 - 2rρ cos α)μ
0
0
(∫π
∫
)-1
∏
2μ
i
× sin2n-2 φ dφ
Θ
(Θ2i,1 + Θ2i,2)-μi dS ρ2(n-|μ|-1 dρ.
i,2
i=1
0
S1(2n)
Обозначив
(∫π
∫
)-1
∏
2μ
i
C-1(n,μ) = C(μ)
sin2n-2 ϕ dϕ
Θ
(Θ2i,1 + Θ2i,2)-μi dS,
i,2
i=1
0
S1(2n)
получим утверждение свойства 5.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
1662
ЛЯХОВ и др.
6. Фундаментальное решение оператора B-γ. Пусть u и v - чётные суммируемые
с весом x-γ функции. Весовую билинейную форму со слабой особенностью определим равен-
ством (8).
Предполагая чётные функции u и v суммируемыми вместе со всеми производными и
дважды интегрируя по частям, получаем равенство
(B-gu, v)-γ = (u, B-γ v)-γ .
Здесь внеинтегральные члены исчезли в силу суммируемости и нечётности функций u′ и v′
(см. теорему 1).
Фундаментальное решение оператора B-γ определяется как обобщённая функция E ∈ S′ev,
удовлетворяющая равенству
B-γE = δ-γ
или, что то же самое,
(E, B-γ ϕ)-γ = ϕ(0) для любой ϕ ∈ Sev([0, ∞)).
(16)
Следующая теорема представляет собой одно из возможных обобщений результатов рабо-
ты [6] (см. в ней теорему 1) на операторы Бесселя отрицательного параметра.
Теорема 6. Пусть Z(t) - регулярный функционал из пространства S′ev и пусть
d
lim
t-γ
Z(t) = 1,
0 ≤ γ < 1.
(17)
t→0
dt
Если в области t > 0 функция Z = Z(x) удовлетворяет однородному сингулярному диффе-
ренциальному уравнению второго порядка
B-γZ(t) = 0,
0 ≤ γ < 1,
(18)
то в смысле обобщённых функций S′ev функция Z есть фундаментальное решение оператора
B-γ в пространстве L-γ2.
Доказательство. Пусть ϕ ∈ Sev. Необходимо проверить выполнение равенства (16). Ин-
тегрируя по частям, получаем
(
)
∫
∞
(
)
d
∞
dZ(t)
d
(Z, B-γ ϕ)-γ = Z(t) t-γ
ϕ(t)
-
t-γ
ϕ(t) dt.
dt
dt
dt
0
0
Поскольку функция ϕ быстро убывает вместе со всеми производными, то
(
)
d
lim
Z(t) t-γ
ϕ(t)
= 0.
t→∞
dt
По условию функция ϕ чётная и бесконечно дифференцируемая в окрестности начала коор-
динат, поэтому ϕ′(x) = O(x), x → 0. Учитывая, что 0 ≤ γ < 1, имеем
Z(t)t-γϕ′(t) = O(t1-γ ), t → 0.
Следовательно,
(
)
d
∞
Z(t) t-γ
ϕ(t)
= 0.
dt
0
Вновь проинтегрировав по частям, получим
∫
∞
(
)
∞
dZ(t)
d
dZ(t)
(Z, Bγ ϕ)γ = -t-γ
ϕ(t)
- tγ
t-γ
ϕ(t)t-γ dt.
dt
dt
dt
0
0
Остаётся воспользоваться условиями (17), (18) и получить равенство (16). Теорема доказана.
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
ПСЕВДОСДВИГ И ФУНДАМЕНТАЛЬНОЕ РЕШЕНИЕ ΔB -ОПЕРАТОРА
1663
Следствие 1. Фундаментальным решением оператора B-γ ,
0 ≤ γ < 1, является функ-
ция
1
E (t) =
tγ+1.
γ+1
Действительно, непосредственно проверяется справедливость (18): B-γAtγ+1 = 0, а из
условия (17) следует, что нормирующая константа A = 1/(γ + 1). Отсюда следует уравне-
ние (6).
При γ = 0 сингулярный дифференциальный оператор Бесселя есть вторая производная.
Из (17) вытекает, что фундаментальным решением второй производной в пространстве обоб-
щённых функций S′ev является функция E(t) = t, что легко проверить: если ϕ ∈ Sev, то
∫
∞
∫
∞
tϕ′′(t)dt = tϕ′(t)|∞0 - ϕ′(t)dt = ϕ(0).
0
0
Замечание 1. В работе [6] получено фундаментальное решение Eγ,m натуральных сте-
пеней (Bγ )m сингулярного дифференциального оператора Бесселя при γ > 0. В частности
Eγ,1 = t1-γ/(1-γ), т.е. фундаментальное решение имеет вырождение (вместо привычной син-
гулярности) порядка O(t1-γ ), t → 0. Тем самым результат следствия 1 справедлив в более
общей форме: при γ > -1 фундаментальным решением оператора Bγ является функция
Eγ = t1-γ/(1 - γ).
Замечание 2. Принципиальное отличие фундаментальных решений операторов B-γ и Bγ
заключается в том, что второе (т.е. при γ > 0) может быть представлено с вырождением в
произвольной точке τ обобщённым сдвигом Пуассона (12) (см. в работе [5] формулу Пуассона),
который не определён при γ = 0 и не существует при γ < 0. Последнее приводит к тому, что
обобщённый сдвиг Пуассона не может применяться к фундаментальному решению оператора
Бесселя с отрицательным параметром -γ.
При γ = 0 сдвиг (12), применённый к чётной непрерывной функции f, есть среднее обыч-
ных сдвигов: учитывая, что j-1/2(x) = cos x, и используя формулу обращения преобразования
Фурье-Ганкеля (см. [4, с. 18]), имеем
∞
∫
2
Tyxf(x) = TyxF-1B[FBf](x)|γ=0 =
(cos ξ(x - y) + cos ξ(x + y))[Fcosf](ξ)dξ =
π
0
1
=
(f(x - y) + f(x + y)).
2
7. T-Псевдосдвиг фундаментального решения оператора B-γ . Пусть -1<-γ < 0,
μ = (γ)/2. Согласно формуле (6) фундаментальное решение оператора B-γ - это функция
→ y) = x2 + y2 - 2xy cos α и
определение (14):
(√
)γ+1
∫
π
(xy)γ+1
→ y) sinγ+1 αdα
∫
π
C(g)
C(g)(xy)γ+1
(xy)γ+1
TyE(x) =
(√
)2(γ+1)/2
=
sinγ+1 α dα =
γ+1
γ+1
γ+1
0
→ y)
0
Таким образом, функция TyE(x) - фундаментальное решение в пространстве R2 по каж-
дой из переменных. При этом
(
)
(
)
(yx)γ+1
xγ+1
(B-γ TyE (x), ϕ(x))-γ = B-γ
, ϕ(x)
=yγ+1
B-γ
, ϕ(x)
= yγ+1ϕ(0).
γ+1
γ+1
-γ
-γ
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
1664
ЛЯХОВ и др.
8. Фундаментальное решение ΔB-оператора Киприянова.
Теорема 7. Пусть -γ = -(γ1, . . . , γn), 0 > -γi > -1 и n-|γ| > 0. Тогда фундаменталь-
ным решением ΔB-оператора Киприянова является регулярная S′ev,-γ-обобщённая функция
⎧
⎪
|x|2-n-|γ|
⎨
,
n + |γ| = 2,
(2 - n - |γ|)|S+(n)|γ
1
En,γ(x) =
(19)
⎪
1
1
⎩
ln
,
n + |γ| = 2,
|S+1(n)|γ
|x|
где площадь нагруженной сферы
∫
∏
∏
1
|S+1(n)|γ =
Γ((αk + 1)/2).
xγki dS =
n-1Γ((n + |α|)/2)
2
k=1
i=1
S+1(n)={x:|x|=1,xi>0}
Доказательство. Используем весовую линейную форму (9). Необходимо найти функцию
En,γ, для которой
(ΔBEn,-γ, ϕ)-γ = ϕ(0).
Предполагая функции En,-γ и ϕ = ϕ(|x|) радиальными и учитывая (1), видим, что предыду-
щее соотношение эквивалентно равенству
|S+1|γ (Bn-|γ|-1En,γ(r), ϕ(r))-γ = ϕ(0).
Пусть n + |γ| = 2. Но тогда n + |γ| - 1 = 1, а это даёт возможность воспользоваться след-
ствием 1 при n + |γ| < 1 или замечанием 1 при n + |γ| > 1. Искомая функция имеет вид
2-n-|γ|
|x|
En,γ(x) =
(2 - n - |γ|)|S1(n)|γ
При n+|γ| = 2 применяется теорема о фундаментальном решении оператора Bγ из работы [6],
поскольку фундаментальным решением оператора Бесселя Bn+|γ|-1 при n + |γ| = 2 будет
функция En,γ(x) = ln |x|. Теорема доказана.
Отметим, что для оператора Киприянова не найден коммутирующий с ним обобщённый
сдвиг, а определённый выше T-псевдосдвиг, как это следует из п. 7, не определяет фунда-
ментального решения оператора B-γ с центром в точке y ∈ R+n, не совпадающей с началом
координат.
Следствие 2. Пусть -1 < -γi < 0 и n - |γ| > 1. Фундаментальное решение ΔB -опера-
тора Киприянова с центром на конусе |x| = |y| имеет вид T|y||x|En,-γ, где
r
- обобщённый
сдвиг Пуассона (12), а En,-γ - функция (19).
Доказательство. При n - |γ| > 1 оператор Бесселя Bn-|γ|-1 имеет положительный
параметр, поэтому ему отвечает обобщённый сдвиг Пуассона (12) порядка α = n - |γ| - 1 > 0,
коммутирующий с Bn-|γ|-1.
9. Обобщение для ΔB-оператора Лапласа-Бесселя-Киприянова. Рассмотрим опе-
ратор ΔB с параметрами γ > -1 в евклидовом пространстве R+n. Будем предполагать, что
хотя бы один из операторов Бесселя, входящий в ΔB, имеет отрицательный параметр. Тогда
нельзя определить фундаментальное решение этого оператора с особенностью в точке. Соглас-
но (1) имеем ΔBu(|x|) = Bn+|γ|-1u(r), r = |x|. Тем самым фундаментальное решение ΔB -
оператора найдётся в виде соответствующего фундаментального решения оператора Бесселя
с параметром n + |γ| - 1.
В пространстве R+n = {n ≥ 2, x : xi > 0} введём весовую линейную форму
∫
∏
(u, v)γ =
u(x)v(x)xγ dx, xγ =
xγii, γi > -1.
i=1
n
Rn
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
ПСЕВДОСДВИГ И ФУНДАМЕНТАЛЬНОЕ РЕШЕНИЕ ΔB -ОПЕРАТОРА
1665
Фундаментальное решение ΔB-оператора Лапласа-Бесселя-Киприянова с центром в про-
извольной точке n-мерного полупространства R+n можно получить из фундаментального ре-
шения (19) при дополнительном условии
n + |γ| - 1 > 0.
(20)
Теорема 8. При выполнении условия (20) фундаментальное решение ΔB -оператора Ла-
пласа-Бесселя-Киприянова с особенностью в точке на конусе |x| = |y|, x, y ∈ R+n, имеет вид
⎧
⎪
|x|2-n-|γ|
,
n + |γ| = 2,
⎪
(2 - n - |γ|)|S1(n)|γ
En,γ(|x|,|y|) = T|y||x|
⎪
1
1
⎩
ln
,
n + |γ| = 2,
|S+1(n)|γ
|x|
где T|y| - обобщённый сдвиг Пуассона (12), отвечающий параметру n+|γ|-1 > 0, |S1(n)|α -
площадь нагруженной сферы.
Доказательство. Введём сферические координаты x = rΘ,
|Θ| = 1. Согласно форму-
√
ле (1) имеем ΔBu(|x|) = Bn-|γ|-1u(r), r =
x21 + ... + x2n. Отсюда, предположив функцию
ϕ = ϕ(|x|) радиальной и учтя (20), получим равенства
(ΔBE, ϕ)γ = (ΔBT|y||x|E, ϕ)γ = |S+1|γ (Bn+|α|-1TρrE(r), ϕ(r))-γ =
= |S+1 |γ (Tρr Bn+|α|-1E(r), ϕ(r))-γ = |S+1 |γ (Bn+|α|-1E(r), Tρr ϕ(r))-γ = ϕ(ρ).
Теорема доказана.
СПИСОК ЛИТЕРАТУРЫ
1. Ляхов Л.Н., Санина Е.Л. Оператор Киприянова-Бельтрами с отрицательной размерностью опе-
ратора Бесселя и сингулярная задача Дирихле для B-гармонического уравнения // Дифференц.
уравнения. 2020. Т. 56. № 12. C. 1610-1520.
2. Metzler R., Glöckle W.G., Nonnenmacher T.F. Fractional model equation for anomalous diffusion
// Phys. A. Stat. Mech. and its Appl. 1994. V. 211. № 1. P. 13-24.
3. Киприянов И.А., Кононенко В.И. О фундаментальных решениях некоторых сингулярных уравне-
ний в частных производных // Дифференц. уравнения. 1969. Т. 5. № 8. С. 1470-1483.
4. Киприянов И.А. Сингулярные эллиптические краевые задачи. М., 1997.
5. Левитан Б.М. Разложение по функциям Бесселя в ряды и интегралы Фурье // Успехи мат. наук.
1951. Т. 6. Вып. 2 (42). С. 102-143.
6. Ляхов Л.Н. Фундаментальные решения сингулярных дифференциальных уравнений с DB опера-
тором Бесселя // Тр. Мат. ин-та им. В.А. Стеклова. 2012. Т. 278. С. 148-160.
7. Какичев В.А. О свертках для интегральных преобразований // Изв. АН БССР. Сер. физ-мат. наук.
1967. № 2. С. 48-57.
8. Бритвина Л.Е. Полисвертки преобразования Ханкеля и дифференциальные операторы // Докл.
РАН. 2002. Т. 382. № 3. С. 298-300.
9. Бритвина Л.Е. О некоторых полисвертках, порожденных преобразованием Ханкеля // Мат. замет-
ки. 2004. Т. 76. Вып. 1. С. 20-26.
10. Britvina L.Y. Generalized shift operators generated by convolutions of integral transforms // Current
Trends in Analysis and its Applications. Trends in Mathematics / Eds. V. Mityushev, M. Ruzhansky.
Cham, 2015.
11. Толстов Г.П. Ряды Фурье. М., 1980.
12. Сабитов К.Б., Зайцева Н.В. Вторая начально-граничная задача B-гиперболического уравнения
// Изв. вузов. Математика. 2019. № 10. С. 75-86.
Воронежский государственный университет,
Поступила в редакцию 26.05.2022 г.
Елецкий государственный университет
После доработки 26.05.2022 г.
имени И.А. Бунина,
Принята к публикации 21.10.2022 г.
Липецкий государственный педагогический университет
имени П.П. Семенова-Тян-Шанского
6
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022