1. 引言
分数阶微积分是数学领域的一个重要分支,最初的分数阶微分算子主要有三种形式:Grumwald- Letnikov定义、Riemann-Liouville定义和Caputo定义。Caputo型分数阶导数在解决实际问题时应用更为方便。于2011年,Wang & Li [1] 参照Caputo型分数阶导数提出了一种新的导数——记忆依赖型导数。与分数阶导数相比,其核函数可以根据实际选择,且记忆依赖区间不会随时间t增长而增大,始终将依赖区间集中在与过去状态有关的某个有限的时间段
,其中τ为时滞。记忆依赖型导数概念的提出引起了不少学者的关注,Ezzat等人在 [2] [3] [4] 中将记忆依赖型导数用于广义热粘弹性等问题。此外,Wang & Li还提出了记忆依赖型微分方程,并在 [5] 中对其做了深入研究。若将其引入偏微分方程结果会如何呢?
经典的弦振动方程和热传导方程的解都有明显的特性,前者解的振幅具有随时间不减的波动性,后者解的振幅具有衰减性。当将关于时间的变化率换成记忆依赖型导数形式时其解的性态到底会发生怎样的变化?本文采用数值方法进行探索。先分析扩散系数及时滞对其数值解的影响;后将其与经典的弦振动方程、一维热传导方程及Caputo型分数阶偏微分方程进行了比较。
2. 记忆依赖型偏微分方程
2.1. 记忆依赖型导数(MD)
基于Caputo型分数阶微分算子,即:
(1)
其中核函数
,
表示通常的m阶导数,Wang & Li给出了记忆依赖型导数的定义。
定义1 [1] [6] :设m是一个自然数(
),对m-次可微函数
,
(2)
被称为m阶记忆依赖型导数,
是时滞,它代表记忆依赖时段的长短。权重函数K是一个m阶连续函数。
与(1)式中导数相比,(2)式中导数的核函数可根据实际情况选择,在 [1] [6] 中通过计算发现,Caputo型导数在大的时间t下将导致长时间变化过程的失效,这是此类导数的固有缺陷,故而记忆依赖型导数较之更为合理。
此处取
![](//html.hanspub.org/file/26-2620455x17_hanspub.png)
对(2)应用积分的平均值公式得:
![](//html.hanspub.org/file/26-2620455x18_hanspub.png)
其中
。由此可见原记忆依赖型导数大约等于
的一半。将上式中的
用
代替后量值基本保持不变。由此可得改进后的定义。
定义2:设m是一个自然数(
),对m-次可微函数
,
(3)
被称为m阶记忆依赖型导数,
是时滞,它代表记忆依赖时段的长短。权重函数K是一个m阶连续函数。
2.2. 记忆依赖型偏微分方程
在 [1] [6] 中已指出记忆依赖型导数相比于分数阶导数在刻画记忆效应方面更具表现力,相应的记忆依赖型微分方程也具有更强的表现力,由其构成的偏微分方程也应具有较强的表现力,而经典的弦振动方程和热传导方程数值解有明显的特性,由此,启发我们利用(3)式中的记忆依赖型导数代替原来的时间导数,探讨其解的性态变化。
探讨如下初边值问题:
(4)
其中
为二阶的(3)式中的记忆依赖型导数,
取的是一次函数形式,即:
![](//html.hanspub.org/file/26-2620455x30_hanspub.png)
通过计算得,
。易发现,除扩散系数
外,时滞
是影响问题(4)数值解的主要因素。
利用有限差分法将问题(4)进行离散化,取∆t为时间步长,∆x为空间步长,令
,
,
,
,得下式:
![](//html.hanspub.org/file/26-2620455x38_hanspub.png)
其中一式中
。
固定
时,选取时滞为15个步长不变,分别选取扩散系数
时,得图1。
由图1可见,记忆依赖型偏微方程的数值解具有波动性,且随着时间增加其振幅逐渐减小,时滞为15个步长不变时,其振幅随
的增加而减小。
进一步,固定
及
时,分别选取
(其中
为时间步长)时,由MATLAB软件实现得图2。
由图2可见,在
时,随着时滞的增加,解的振度随之降低,趋于稳定所需的时间也随之减少。由图1和图2可见,记忆依赖型偏微分方程的解既有波动性又有衰减性,其振幅随扩散系数和时滞的增大而减少,趋于稳定的时间也随之缩短。
3. 与经典的热传导方程及波动方程的比较
将问题(4)与经典的热传导方程:
(5)
![](//html.hanspub.org/file/26-2620455x49_hanspub.png)
Figure 1. The influence of different diffusion coefficients on the solutions
图1. 不同扩散系数对解的影响情况
![](//html.hanspub.org/file/26-2620455x50_hanspub.png)
Figure 2. Influence of different time delay on solution
图2. 不同时滞对解的影响情况
及弦振动方程:
(6)
进行比较。利用有限差分法将上述问题进行离散化,并由MATLAB软件实现得到图3。
由图3可见,波动方程随时间以波的形式传播,同时能量不衰减;热传导方程随时间为能量衰减型,并且衰减速度很快;记忆依赖型偏微分方程介于波动方程及热传导方程之间,既有波动性又有衰减性(见图3(d),这种变化符合“加权周期”函数 [7] 的特性。这种现象的产生应与记忆依赖型导数依赖于时滞有关。通过比较分析,发现记忆依赖型导数对于描述现实中与时滞相关的能量逐步衰减的波(如击鼓时产生的波)更有表现力。
4. 与分数阶偏微分方程的对比
记忆依赖型导数在Caputo型分数阶微分算子的基础上提出,两者依赖区间及核函数都有所不同,故而两者分别构成的偏微分方程数值解也应有所差异。将问题(4)与以下Caputo型分数阶微分算子构成的分数阶偏微分方程进行比较:
(7)
![](//html.hanspub.org/file/26-2620455x53_hanspub.png)
Figure 3. When a, k = 1, string vibration equation (a), the heat conduction equation (b), memory dependent partial differential equations (c) and the numerical solution of memory-dependent partial differential equations in
(d)
图3. a, k = 1的弦振动方程(a)、热传导方程(b)、记忆依赖型偏微分方程(c)及记忆依赖型偏微分方程于
的数值解图像(d)
利用以下引理及有限差分法对问题(7)进行离散 [8] :
引理 [9] 设
,则
![](//html.hanspub.org/file/26-2620455x57_hanspub.png)
其中
![](//html.hanspub.org/file/26-2620455x58_hanspub.png)
在此我们选取
,将离散化后的问题(7)利用MATLAB实现得图4。
对比图3(c),图3(d)与图4,从整体形态上发现,记忆依赖型导数与Caputo型微分算子虽都为分数阶的,但所构成的偏微分方程数值解明显有所不同,两者都具有波动性,以波的形式向前传播,并且在传播过程中能量都在衰减,在扩散系数相同时,相对于分数阶偏微分方程数值解的快速衰减,记忆依赖型偏微分方程数值解的衰减速度明显慢很多。
5. 结论
在记忆依赖型导数基础上提出记忆依赖型偏微分方程,通过数值方法与弦振动方程、一维热传导方程及分数阶偏微分方程相比,其解介于两者之间,既具有波动性又具有衰减性,且振幅衰减速度慢于热传导方程。其振幅及趋于稳定的时间随时滞和扩散系数的增大而减少,这更近似于现实生活中声音的传播等现象。另外,此处实验结果显示,其振幅远大于弦振动方程和热传导方程的振幅,这一现象应与记
![](//html.hanspub.org/file/26-2620455x60_hanspub.png)
Figure 4. Numerical solutions of FPDE composed of Caputo type differential operators when
and ![](//html.hanspub.org/file/26-2620455x62_hanspub.png)
图4.
时,Caputo型微分算子构成的分数阶偏微方程数值解
忆依赖型导数的时滞及核函数有关,值得进一步探讨。与分数阶偏微分方程数值解相比,其解的衰减速度更缓慢,且波动性更明显。
基金项目
国家自然科学基金项目(No.41376030)。