一类具p(t)-Laplacian算子的Hadamard分数阶微分方程耦合系统边值问题解的唯一性
Uniqueness of Solutions for Boundary Value Problems of Hadamard Fractional Differential Equation Coupled Systems with p(t)-Laplacian Operator
摘要: 本文研究一类含p(t)-Laplacian算子的Hadamard分数阶微分方程耦合系统非局部边值问题。通过Hadamard分数阶积分与微分的性质,将边值问题转化为等价的积分算子方程;利用Banach压缩映射原理,建立系统解存在唯一的充分条件;构造具体耦合系统实例,通过验证条件及数值计算,证实结论的有效性。研究结果为同类分数阶耦合系统的定性分析提供了新的思路与方法。
Abstract: This paper investigates a class of nonlocal boundary value problems for Hadamard fractional differential equation coupled systems involving the p(t)-Laplacian operator. By virtue of the properties of Hadamard fractional integration and differentiation, the boundary value problem is transformed into an equivalent integral operator equation. Using the Banach contraction mapping principle, sufficient conditions for the existence and uniqueness of solutions to the system are established. A specific example of the coupled system is constructed, and the validity of the conclusions is verified through condition checking and numerical calculation. The research results provide new ideas and methods for the qualitative analysis of similar fractional coupled systems.
文章引用:张婷婷, 王来全, 周璞铉, 鲁辉. 一类具p(t)-Laplacian算子的Hadamard分数阶微分方程耦合系统边值问题解的唯一性[J]. 应用数学进展, 2026, 15(1): 13-19. https://doi.org/10.12677/aam.2026.151003

1. 引言

分数阶微分方程因其在物理、工程、生物等领域的广泛应用,成为近年来非线性分析领域的研究热点[1]。不同于整数阶微分方程仅能刻画局部变化率,分数阶导数具备非局部记忆特性,能够精准描述黏弹性材料变形、反常扩散、生物种群演化等具有历史依赖性的复杂过程[2]。Hadamard分数阶导数以对数函数为核,在处理定义于正实数域且具有幂律型尺度不变性的系统时展现出独特优势[3] [4],例如在研究电力系统中绝缘介质的老化特性、地质工程中岩石的蠕变损伤过程等问题时,Hadamard分数阶模型的拟合精度显著高于传统整数阶模型。同时,p(t)-Laplacian算子能够刻画非均匀介质中的扩散、渗透等现象,其变指数特性可有效反映介质属性随空间位置变化的规律,例如在多孔介质流体输运问题中,介质孔隙率的空间异质性会导致扩散系数随位置改变,此时p(t)-Laplacian算子对应的方程模型更贴合实际场景[5]

在实际工程与物理系统中,单一方程往往难以完整描述多场耦合的复杂现象,耦合系统的研究更具现实意义。例如在智能材料与结构领域,压电陶瓷–黏弹性基体复合结构的力–电耦合振动问题,需要同时考虑机械形变的黏弹性记忆效应与电场分布的非均匀扩散特性,这类问题可抽象为含p(t)-Laplacian算子的Hadamard分数阶微分方程耦合系统。目前,关于分数阶微分方程耦合系统的研究已取得一定成果[6] [7],但现有成果多聚焦于Riemann-Liouville或Caputo型分数阶导数,且算子多为常指数P-Laplacian算子,结合p(t)-Laplacian算子与Hadamard分数阶导数的耦合系统边值问题研究尚不多见。此外,现有研究多采用单调迭代法或上下解方法,针对非局部边值条件下解的唯一性分析相对匮乏。本文针对一类具p(t)-Laplacian算子的Hadamard分数阶微分方程耦合系统非局部边值问题,通过构造等价积分算子、应用Banach压缩映射原理,建立解的存在唯一性条件,并通过具体实例验证结论的正确性。本文的研究不仅丰富了分数阶微分方程耦合系统的理论体系,也为相关工程问题的建模与分析提供了新的思路与方法。

2. 预备知识

定义2.1 [4] α>0 ,函数 u( t ) C 1 ( [ 1,T ],R ) ,则 α 阶Hadamard分数阶积分定义为:

H 1 + α u( t )= 1 Γ( α ) 1 t ( ln t s ) α1 u( s ) s ds

其中Gamma(∙)为Gamma函数。

定义2.2 [4] α( 0,1 ) ,函数 u( t ) C 1 ( [ 1,T ],R ) ,则 α 阶Hadamard分数阶导数定义为:

D 1 + α u( t )= 1 Γ( 1α ) ( t d dt ) 1 t ( ln t s ) α u( s ) s ds (2.2)

引理2.1 [5] p(t)-Laplacian算子性质

p( t ):[ 1,T ]( 1,+ ) 为连续函数, q( t )= p( t ) p( t )1 p( t ) 的共轭指数函数, p( t ) -Laplacian算子 Δ p( t ) ( u )= | u | p( t )2 u ,满足以下性质:

(1) 单调性:若 uv0 ,则 ( Δ p( t ) ( u ) Δ p( t ) ( v ) )( uv )0

(2) 逆算子存在性: Δ p( t ) 1 ( z )= | z | q( t )2 z ,且 Δ p( t ) Δ p( t ) 1 ( z )=z

引理2.2 [8] Banach压缩映射原理

( X,d ) 为完备度量空间, T:XX 为压缩映射,即存在常数 θ( 0,1 ) ,使得对任意 x,yX ,有 d( Tx,Ty )θd( x,y ) TX中存在唯一不动点。

3. 边值问题及等价积分算子构造

3.1. 问题描述

考虑如下具p(t)-Laplacian算子的Hadamard分数阶微分方程耦合系统非局部边值问题:

{ D 1 + α ( Δ p 1 ( t ) ( D 1 + α u( t ) ) )+f( t,v( t ), D 1 + α v( t ) )=0,t( 1,T ) D 1 + β ( Δ p 2 ( t ) ( D 1 + β v( t ) ) )+g( t,u( t ), D 1 + β u( t ) )=0,t( 1,T ) u( 1 )=0, Δ p 1 ( T ) ( D 1 + α u( T ) )= 1 T h 1 ( s )u( s )ds v( 1 )=0, Δ p 2 ( T ) ( D 1 + β v( T ) )= 1 T h 2 ( s )v( s )ds (3.1)

其中 0<α,β<1 p 1 ( t ), p 2 ( t ):[ 1,T ]( 1,+ ) 连续, h 1 ( t ), h 2 ( t )C( [ 1,T ],[ 0,+ ) ) fg:[ 1,T ]×R×RR 连续。

3.2. 等价积分算子构造

对任意满足可积性与可微性条件的函数 ϕ( t ) α 阶与 β 阶Hadamard分数阶微分、积分算子满足:

D 1 + α ( H 1 + α ϕ( t ) )=ϕ( t ), H 1 + α ( D 1 + α ϕ( t ) )=ϕ( t )ϕ( 1 )

D 1 + β ( H 1 + β ϕ( t ) )=ϕ( t ), H 1 + β ( D 1 + β ϕ( t ) )=ϕ( t )ϕ( 1 )

方程两边对 u( t ) 的方程两端作用 H 1 + α ,由积分线性得:

H 1 + α ( D 1 + α ( Δ p 1 ( t ) ( D 1 + α u( t ) ) ) )= H 1 + α f( t,v( t ), D 1 + α v( t ) ) (3.2)

根据 α 阶Hadamard微积分互逆性,左端展开为:

H 1 + α ( D 1 + α ( Δ p 1 ( t ) ( D 1 + α u( t ) ) ) )= Δ p 1 ( t ) ( D 1 + α u( t ) ) Δ p 1 ( 1 ) ( D 1 + α u( 1 ) )

其中将 Δ p 1 ( t ) ( D 1 + α u( t ) ) 视为整体函数 ϕ( t ) ,代入互逆性公式得到边界项 Δ p 1 ( 1 ) ( D 1 + α u( 1 ) )

D 1 + α u( 1 )=0 ,得 Δ p 1 ( 1 ) ( D 1 + α u( 1 ) )= | 0 | p 1 ( 1 )2 0=0 ,因此左端化简为 Δ p 1 ( t ) ( D 1 + α u( t ) )

将化简结果代回式(3.1),得:

Δ p 1 ( t ) ( D 1 + α u( t ) )= H 1 + α f( t,v( t ), D 1 + α v( t ) ) (3.3)

再应用p(t)-Laplacian逆算子,得逆算子形式: Δ p( t ) 1 ( z )= | z | q( t )2 z ,其中 q( t )= p( t ) p( t )1 p( t ) 的共轭指数,满足 1 p( t ) + 1 q( t ) =1 ;对任意z,有 Δ p( t ) 1 ( Δ p( t ) ( z ) )=z Δ p( t ) ( Δ p( t ) 1 ( z ) ) 为简化书写,记:

F( t,v )= H 1 + α f( t,v( t ), D 1 + α v( t ) ),

则式(3.3)简写为

Δ p 1 ( t ) ( D 1 + α u( t ) )=F( t,v )

对上式两端同时作用 Δ p 1 ( t ) 1 ,得到:

Δ p 1 ( t ) 1 ( Δ p 1 ( t ) ( D 1 + α u( t ) ) )= Δ p 1 ( t ) 1 ( F( t,v ) )

p(t)-Laplacian算子的互逆性 Δ p( t ) 1 ( Δ p( t ) ( z ) )=z ,左端可直接化简为 D 1 + α u( t ) ,即:

D 1 + α u( t )= Δ p 1 ( t ) 1 ( F( t,v ) )

D 1 + α u( t )= Δ p 1 ( t ) 1 ( H 1 + α f( t,v( t ), D 1 + α v( t ) ) )

再次应用 α 阶Hadamard分数阶积分,结合 u( 1 )=0 ,得:

u( t )= H 1 + α ( Δ p 1 ( t ) 1 ( H 1 + α f( t,v( t ), D 1 + α v( t ) ) ) ) (3.4)

同理,对(3.1)中第二个方程重复上述步骤,结合边值条件 v( 1 )=0 ,得:

v( t )= H 1 + β ( Δ p 2 ( t ) 1 ( H 1 + β g( t,u( t ), D 1 + β u( t ) ) ) ) (3.5)

将非局部边值条件代入(3.4) (3.5),最终得到(3.1)的等价积分算子形式:

{ u( t )= T 1 ( u,v )( t ) v( t )= T 2 ( u,v )( t ) (3.6)

其中:

T 1 ( u,v )( t )= 1 Γ 2 ( α ) 1 t ( ln t s ) α1 1 s | 1 s ( ln s τ ) α1 f( τ,v( τ ), D 1 + α v( τ ) ) τ dτ | q 1 ( s )2 1 s ( ln s τ ) α1 f( τ,v( τ ), D 1 + α v( τ ) ) τ dτds (3.7)

T 2 ( u,v )( t )= 1 Γ 2 ( β ) 1 t ( ln t s ) β1 1 s | 1 s ( ln s τ ) β1 g( τ,u( τ ), D 1 + β u( τ ) ) τ dτ | q 2 ( s )2 1 s ( ln s τ ) β1 g( τ,u( τ ), D 1 + β u( τ ) ) τ dτds (3.8)

其中 q 1 ( t )= p 1 ( t ) p 1 ( t )1 q 2 ( t )= p 2 ( t ) p 2 ( t )1 为共轭指数函数。

4. 解的存在唯一性分析

X=C( [ 1,T ],R )×C( [ 1,T ],R ) ,对任意 ( u,v ),( u * , v * )X ,定义度量:

d( ( u,v ),( u * , v * ) )= max t[ 1,T ] | u( t ) u * ( t ) |+ max t[ 1,T ] | v( t ) v * ( t ) | (4.1)

易证 ( X,d ) 为完备度量空间。

定义算子 T:XX T( u,v )=( T 1 ( u,v ), T 2 ( u,v ) ) ,则(3.1)的解等价于算子T的不动点。

假设以下条件成立:

(H1) 存在常数 L 1 , L 2 >0 ,使得对任意 t[ 1,T ] v 1 , v 2 , w 1 , w 2 R ,有:

| f( t, v 1 , w 1 )f( t, v 2 , w 2 ) | L 1 | v 1 v 2 |+ L 2 | w 1 w 2 |

(H2) 存在常数 M 1 , M 2 >0 ,使得对任意 t[ 1,T ] u 1 , u 2 , z 1 , z 2 R ,有:

| g( t, u 1 , z 1 )g( t, u 2 , z 2 ) | M 1 | u 1 u 2 |+ M 2 | z 1 z 2 |

(H3) 记 C 1 = ( T1 ) α L 1 Γ 2 ( α+1 ) max t[ 1,T ] q 1 ( t ) C 2 = ( T1 ) β M 1 Γ 2 ( β+1 ) max t[ 1,T ] q 2 ( t ) ,且 C 1 + C 2 <1

定理4.1若条件(H1)~(H3)成立,则耦合系统(3.1)在X中存在唯一解。

证明:对任意 ( u 1 , v 1 ),( u 2 , v 2 )X ,计算 | T 1 ( u 1 , v 1 )( t ) T 1 ( u 2 , v 2 )( t ) |

T 1 ( u 1 , v 1 )( t ) T 1 ( u 2 , v 2 )( t ) 1 Γ 2 ( α ) 1 t ( ln t s ) α1 1 s | Δ p 1 ( s ) 1 ( H 1 + α f 1 ) Δ p 1 ( s ) 1 ( H 1 + α f 2 ) |ds max t[ 1,T ] q 1 ( t ) Γ 2 ( α ) 1 t ( ln t s ) α1 1 s | H 1 + α f 1 H 1 + α f 2 |ds ( T1 ) α L 1 Γ 2 ( α+1 ) max t[ 1,T ] q 1 ( t ) max t[ 1,T ] | v 1 ( t ) v 2 ( t ) | C 1 d( ( u 1 , v 1 ),( u 2 , v 2 ) )

同理, | T 2 ( u 1 , v 1 )( t ) T 2 ( u 2 , v 2 )( t ) | C 2 d( ( u 1 , v 1 ),( u 2 , v 2 ) )

因此:

d( T( u 1 , v 1 ),T( u 2 , v 2 ) )( C 1 + C 2 )d( ( u 1 , v 1 ),( u 2 , v 2 ) ) (4.2)

由条件(H3), C 1 + C 2 <1 ,故T为压缩映射。根据Banach压缩映射原理,TX中存在唯一不动点,即耦合系统(2.1)存在唯一解。

注记

条件(H3)作为解存在唯一性的核心约束条件,其本质是对算子压缩性进行量化限制,关键参数的取值直接决定该条件是否成立:其一,分数阶导数阶数 α,β 通过Gamma函数 Γ( α+1 ),Γ( β+1 ) 影响条件中的系数,阶数越大, Γ( α+1 ),Γ( β+1 ) 的值越大, C 1 , C 2 的值越小,越容易满足 C 1 + C 2 <1 的要求;其二,区间长度 T1 C 1 , C 2 呈正相关,当T取值过大时, C 1 + C 2 可能超过1,此时压缩映射条件不再成立,这表明系统解的唯一性仅在一定的时间或空间范围内有效;其三,变指数函数 p 1 ( t ), p 2 ( t ) 的范围决定了共轭指数 q 1 ( t ), q 2 ( t ) 的上界, max q 1 ( t ),max q 2 ( t ) 越小, C 1 , C 2 的值越小,更有利于唯一性条件的满足,这体现了非均匀介质属性对系统解的定性行为的调控作用。

本文的研究结果为分数阶耦合系统的定性分析提供了理论支撑,可进一步推广至更一般形式的分数阶微分方程系统,例如,采用Krasnoselskii不动点定理或Schaefer不动点定理,弱化Lipschitz条件的限制,研究非线性项满足次线性或超线性增长时系统解的存在性。

5. 实例验证

5.1. 实例构造

T=e α=β= 1 2 p 1 ( t )= p 2 ( t )=2 (此时p(t)-Laplacian算子退化为普通Laplacian算子, Δ 2 ( u )=u ,逆算子 Δ 2 1 ( z )=z h 1 ( t )= h 2 ( t )=0 (边值条件退化为齐次边值),构造耦合系统:

{ D 1 + 1 2 ( D 1 + 1 2 u( t ) )+ 1 4 v( t )+ 1 4 D 1 + 1 2 v( t )=0,t( 1,e ) D 1 + 1 2 ( D 1 + 1 2 v( t ) )+ 1 4 u( t )+ 1 4 D 1 + 1 2 u( t )=0,t( 1,e ) u( 1 )=0, D 1 + 1 2 u( e )=0 v( 1 )=0, D 1 + 1 2 v( e )=0 (5.1)

5.2. 条件验证

对比(4.1),此处 f( t,v,w )= 1 4 v+ 1 4 w g( t,u,z )= 1 4 u+ 1 4 z

条件(H1):取 L 1 = L 2 = 1 4 ,满足Lipschitz条件;

条件(H2):取 M 1 = M 2 = 1 4 ,满足Lipschitz条件;

条件(H3): Γ( 1 2 +1 )= π 2 ( e1 ) 1 2 1.284 max q 1 ( t )=max q 2 ( t )=2 ,计算得:

C 1 = C 2 = 1.284× 1 4 ( π 2 ) 2 ×20.32

C 1 + C 2 0.64<1 ,满足条件(H3)。

基金项目

昌吉职业技术学院院级科研课题《具P(t)-Laplacian算子的分数阶微分方程边值问题解的存在性研究》CJZY2025026。

NOTES

*通讯作者。

参考文献

[1] Podlubny, I. (1999) Fractional Differential Equations. Academic Press.
[2] Kilbas, A.A., Srivastava, H.M. and Trujillo, J.J. (2006) Theory and Applications of Fractional Differential Equations. Elsevier.
[3] Hadamard, J. (1982) Essai sur l’étude des fonctions données par leurs dérivées successives. Journal de Mathématiques Pures et Appliquées, 8, 101-186.
[4] 王书超, 胡卫敏, 胡芳芳. 具p-Laplacian算子的Hadamard型分数阶微分方程边值问题的解[J]. 数学的实践与认识, 2025, 55(1): 233-244.
[5] 张迪, 刘文斌, 张伟. 一类带p(t)-Laplacian算子的分数阶微分方程边值问题解的存在性[J]. 河南师范大学学报, 2018, 46(2): 17-21.
[6] 王宁, 周宗福. 具有p-Laplacian算子的分数阶微分耦合系统边值问题的正解(英文) [J]. 应用数学, 2023, 36(2): 530-539.
[7] 潘欣媛, 何小飞, 陈国平, 等. 一类Caputo-Hadamard型分数阶微分方程耦合系统边值问题[J]. 应用数学, 2023, 36(4): 987-996.
[8] 张恭庆, 林源渠. 泛函分析讲义[M]. 北京: 北京大学出版社, 2004.