若将经典的弦振动方程与热传导方程中关于时间的变化率替换成记忆依赖型导数形式,其解的性态会发生怎样的变化呢?相对于普通导数而言,记忆依赖型导数可反映物理过程对过去状态的依赖性。与分数阶导数相比,其核函数的选取更自由,依赖区间也不会随时间的增长而增大,故而记忆依赖型偏微分方程也应具有更强表现力。本文针对核函数为线性函数的情况进行探讨。数值结果显示:1) 其解的性态介于弦振动方程和热传导方程之间,既有波动性又有衰减性,振幅随时滞和扩散系数的增大而减小。2) 与Caputo型分数阶偏微分方程相比,其波动性更强,振幅衰减更慢。 If the rate of change respect to time in the classical string-vibration equation and heat-conduction equation is replaced by the type of memory-dependent derivative, what is the difference for the behavior of the solution? To compare with the ordinary derivative, the memory-dependent type can reflect clearly the dependence of physical process on their past states. To compare with the fractional derivative, the kernel function can be chosen freely and the interval for dependence doesn’t increase with time; so the memory-dependent partial differential equation should have strong expressive force. In this study, the case of kernel function with linear form is considered. Numerical results show that: 1) The characteristics of the solution lie between the string-vibration equation and the heat-conduction equation. It has both fluctuating and decaying properties. The amplitude of it decreases along with the increasing of time-delay and diffusion-coefficient. 2) To compare with the Caputo type of fractional partial differential equation, the fluctuation of the so-lution is stronger and the decay of amplitude is slower.
孙雯雯,王金良
青岛理工大学理学院,ESMD方法及其应用研究所,山东 青岛
收稿日期:2017年7月7日;录用日期:2017年7月24日;发布日期:2017年7月27日
若将经典的弦振动方程与热传导方程中关于时间的变化率替换成记忆依赖型导数形式,其解的性态会发生怎样的变化呢?相对于普通导数而言,记忆依赖型导数可反映物理过程对过去状态的依赖性。与分数阶导数相比,其核函数的选取更自由,依赖区间也不会随时间的增长而增大,故而记忆依赖型偏微分方程也应具有更强表现力。本文针对核函数为线性函数的情况进行探讨。数值结果显示:1) 其解的性态介于弦振动方程和热传导方程之间,既有波动性又有衰减性,振幅随时滞和扩散系数的增大而减小。2) 与Caputo型分数阶偏微分方程相比,其波动性更强,振幅衰减更慢。
关键词 :记忆依赖型导数,记忆依赖型偏微分方程,Caputo型分数阶导数,波动方程,热传导方程
Copyright © 2017 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
分数阶微积分是数学领域的一个重要分支,最初的分数阶微分算子主要有三种形式:Grumwald- Letnikov定义、Riemann-Liouville定义和Caputo定义。Caputo型分数阶导数在解决实际问题时应用更为方便。于2011年,Wang & Li [
经典的弦振动方程和热传导方程的解都有明显的特性,前者解的振幅具有随时间不减的波动性,后者解的振幅具有衰减性。当将关于时间的变化率换成记忆依赖型导数形式时其解的性态到底会发生怎样的变化?本文采用数值方法进行探索。先分析扩散系数及时滞对其数值解的影响;后将其与经典的弦振动方程、一维热传导方程及Caputo型分数阶偏微分方程进行了比较。
基于Caputo型分数阶微分算子,即:
其中核函数
定义1 [
被称为m阶记忆依赖型导数,
与(1)式中导数相比,(2)式中导数的核函数可根据实际情况选择,在 [
此处取
对(2)应用积分的平均值公式得:
其中
定义2:设m是一个自然数(
被称为m阶记忆依赖型导数,
在 [
探讨如下初边值问题:
其中
通过计算得,
利用有限差分法将问题(4)进行离散化,取∆t为时间步长,∆x为空间步长,令
其中一式中
固定
由图1可见,记忆依赖型偏微方程的数值解具有波动性,且随着时间增加其振幅逐渐减小,时滞为15个步长不变时,其振幅随
进一步,固定
由图2可见,在
将问题(4)与经典的热传导方程:
图1. 不同扩散系数对解的影响情况
图2. 不同时滞对解的影响情况
及弦振动方程:
进行比较。利用有限差分法将上述问题进行离散化,并由MATLAB软件实现得到图3。
由图3可见,波动方程随时间以波的形式传播,同时能量不衰减;热传导方程随时间为能量衰减型,并且衰减速度很快;记忆依赖型偏微分方程介于波动方程及热传导方程之间,既有波动性又有衰减性(见图3(d),这种变化符合“加权周期”函数 [
记忆依赖型导数在Caputo型分数阶微分算子的基础上提出,两者依赖区间及核函数都有所不同,故而两者分别构成的偏微分方程数值解也应有所差异。将问题(4)与以下Caputo型分数阶微分算子构成的分数阶偏微分方程进行比较:
图3. a, k = 1的弦振动方程(a)、热传导方程(b)、记忆依赖型偏微分方程(c)及记忆依赖型偏微分方程于
利用以下引理及有限差分法对问题(7)进行离散 [
引理 [
其中
在此我们选取
对比图3(c),图3(d)与图4,从整体形态上发现,记忆依赖型导数与Caputo型微分算子虽都为分数阶的,但所构成的偏微分方程数值解明显有所不同,两者都具有波动性,以波的形式向前传播,并且在传播过程中能量都在衰减,在扩散系数相同时,相对于分数阶偏微分方程数值解的快速衰减,记忆依赖型偏微分方程数值解的衰减速度明显慢很多。
在记忆依赖型导数基础上提出记忆依赖型偏微分方程,通过数值方法与弦振动方程、一维热传导方程及分数阶偏微分方程相比,其解介于两者之间,既具有波动性又具有衰减性,且振幅衰减速度慢于热传导方程。其振幅及趋于稳定的时间随时滞和扩散系数的增大而减少,这更近似于现实生活中声音的传播等现象。另外,此处实验结果显示,其振幅远大于弦振动方程和热传导方程的振幅,这一现象应与记
图4.
忆依赖型导数的时滞及核函数有关,值得进一步探讨。与分数阶偏微分方程数值解相比,其解的衰减速度更缓慢,且波动性更明显。
国家自然科学基金项目(No.41376030)。
孙雯雯,王金良. 记忆依赖型偏微分方程数值解研究On Numerical Solution of the Memory Dependent Partial Differential Equations[J]. 应用数学进展, 2017, 06(04): 637-643. http://dx.doi.org/10.12677/AAM.2017.64074