ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ, 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 =
Bi , Bi =
+
,
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
Соответствующее пространство обобщённых функций будем обозначать Sev. Функцио-
нал Дирака-Киприянова определим равенством (δ , ϕ)α = ϕ(0) для любой функции ϕ ∈ Sev.
Вначале отметим следующее утверждение.
Теорема 1. Пусть γ ∈ [0, 1), α - произвольное действительное число, U и V принад-
лежат основному пространству функций Sev[0, ∞), тогда справедлива формула
(
[
]
)
(γ + α)(α - 1)
(BU, V )α = U, Bγ+2α +
V
x2
α
Доказательство легко достигается интегрированием по частям.
Как видим, оператор Bγ+2α + (γ + α)(α - 1)/x2 оказывается сопряжённым к оператору
B в гильбертовом пространстве Lα2(R+1). Конечно, наиболее востребованы случаи α =
и α = 1. Первый представляется наиболее естественным, так как тогда B является само-
сопряжённым оператором в соответствующем скалярном произведении.
В этой работе приводим результаты исследований только для случая α = -γ, 0 < γ < 1.
Это связано не только с традициями, заложенными в научной школе И.А. Киприянова (см. [4]
и имеющуюся в ней библиографию), но и с определением фундаментального решения, которое
жёстко связано с весовой билинейной формой, где определяется это решение.
Перечислим основные результаты.
1. Получены линейно независимые решения сингулярного уравнения Бесселя Bu(λ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 в пространстве Sev,-γ яв-
ляется регулярный Sev -функционал 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
Bu(λ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. Если xf(x) ∈ L2(0, ∞), то
f (x) = Jμ()ξ-2μ+1 Jμ()y-2μ+1f(y) dy,
(11)
0
0
т.е. следующие J-преобразования Бесселя взаимно обратимы:
f (ξ) = J(γ+1)/2()f(y)y dy, F-1
f ](x) = f(x) =
f (ξ)ξ dξ.
FB [f](ξ)
B
J(γ+1)/2(
0
0
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
1658
ЛЯХОВ и др.
Доказательство. Известен следующий интеграл Ханкеля для функций f ∈ L2(0, ∞):
f (x) = Jμ()ξ dξ Jμ()f(y)y dy.
0
0
Пусть xf(x) ∈ L2(0, ∞), μ = (γ + 1)/2 > 0. Тогда выполняются равенства
Jμ()
Jμ()
f (x)x = Jμ()ξ dξ Jμ(yξf(y)yy dy =
ξ dξ
f (y)yy dy =
()μ
()μ
0
0
0
0
= x Jμ()ξ-2μ+1 Jμ()y-2μ+1f(y)dy,
0
0
откуда следует формула (11). Теорема доказана.
4. Обобщённый T-псевдосдвиг и теорема сложения для J-функций Бесселя. Для
j-функций Бесселя известен следующий вариант теоремы сложения [5]: для ν > -1/2
jν()jν() = Tτξ jν(ξτ),
где сдвиг Пуассона Tτξ определён равенством
π
)
(√
Γ((γ + 1)/2)
Tτ ϕ(t) =
ϕ
t2 + τ2 - 2 cos α sinγ-1 α dα, γ > 0.
(12)
Γ(1/2)Γ(γ/2)
0
Поскольку μ = (γ + 1)/2 > 0, то из (8) имеем
Jμ()Jμ() = ()2μjμ()()2μjμ() = ()2μ()2μ
ξ
jμ() =
(
)
π
ξ2 + τ2 - 2ξτ cos α
Jμ t
()2μ()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μ(),
2μ = γ + 1.
Γ(1/2)Γ(μ + 1/2)
(√
0
ξ2 + τ2 - 2ξτ cos α
Таким образом, имеет место
Теорема 5 (сложения Jμ-функций Бесселя). Справедлива формула
(13)
J(γ+1)/2()J(γ+1)/2() = TxJ(γ+1)/2(), 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
- чётная, дважды непрерывно дифференцируемая функция,
xf(x) ∈ L2(0,∞), 0 < γ < 1. Тогда
B-γ,xTyxf(x) = TyxB-γ,xf(x).
(15)
Доказательство. Из (11), (13) и (7) имеем
B-γ,xTyxf(x) = B-γ,x Jμ()Jμ(
f (ξ)ξ =
0
ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ том 58
№ 12
2022
1660
ЛЯХОВ и др.
=
(2)Jμ()Jμ(
f (ξ)ξ =
(2)TyxJμ(
f (ξ)ξ dξ.
0
0
Здесь множитель (2) можно внести под знак T-псевдосдвига. В результате получим
BTyxf(x) = Tyx((2)Jμ()
f (ξ)ξ =
0
= Tyx(BJμ()
f (ξ)ξ = TyB Jμ(
f (ξ)ξ = TyxBf(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
xii, μ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μii 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 dyii = 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
приведёт к равенству
)
(√
μ
()2|μ|g
r2 + ρ2 - 2cos ϕ1
(f,Ty
g)
=
f (ρ)
×
C(μ)
(r2 + ρ2 - 2 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 φ22 . . . sin φ2n-22n-22n-1,
то
)
(√
μ
π
()2|μ|g
r2 + ρ2 - 2cos α
(f,Ty
g)
=
f (ρ)
sin2|μ| α ×
C(μ)
(r2 + ρ2 - 2 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 ∈ Sev,
удовлетворяющая равенству
BE = δ
или, что то же самое,
(E, B ϕ) = ϕ(0) для любой ϕ ∈ Sev([0, ∞)).
(16)
Следующая теорема представляет собой одно из возможных обобщений результатов рабо-
ты [6] (см. в ней теорему 1) на операторы Бесселя отрицательного параметра.
Теорема 6. Пусть Z(t) - регулярный функционал из пространства Sev и пусть
d
lim
t
Z(t) = 1,
0 γ < 1.
(17)
t→0
dt
Если в области t > 0 функция Z = Z(x) удовлетворяет однородному сингулярному диффе-
ренциальному уравнению второго порядка
BZ(t) = 0,
0 γ < 1,
(18)
то в смысле обобщённых функций Sev функция Z есть фундаментальное решение оператора
B в пространстве L2.
Доказательство. Пусть ϕ ∈ 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): BAtγ+1 = 0, а из
условия (17) следует, что нормирующая константа A = 1/(γ + 1). Отсюда следует уравне-
ние (6).
При γ = 0 сингулярный дифференциальный оператор Бесселя есть вторая производная.
Из (17) вытекает, что фундаментальным решением второй производной в пространстве обоб-
щённых функций Sev является функция E(t) = t, что легко проверить: если ϕ ∈ Sev, то
′′(t)dt =(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](ξ) =
π
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-оператора Киприянова является регулярная Sev,-γ-обобщённая функция
|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