低扩散分数阶微分方程有限差分方法-计算数学专业论文.docx

承诺书 本人郑重声明：所呈交的学位论文，是本人在导师指导下，独立 进行研究工作所取得的成果。尽我所知，除文中已经注明引用的内容 外，本学位论文的研究成果不包含任何他人享有著作权的内容。对本 论文所涉及的研究工作做出贡献的其他个人和集体，均已在文中以明 确方式标明。 本人授权南京航空航天大学可以有权保留送交论文的复印件，允 许论文被查阅和借阅,可以将学位论文的全部或部分内容编入有关数 据库进行检索，可以采用影印、缩印或其他复制手段保存论文。 (保密的学位论文在解密后适用本承诺书) 作者签名： 日 期： 南京航空航天大学硕士学位论文 南京航空航天大学硕士学位论文 低扩散分数阶微分方程的有限差分方法 低扩散分数阶微分方程的有限差分方法 PAGE PAGE 11 PAGE PAGE 10 第一章 分数阶微分方程简介 1.1 分数阶导数定义 [1]分数阶微分方程在工程科学领域中的应用逐渐扩大，特别是近年来，在流体力学、粘弹性 力学、流变学、分数控制系统与分数控制器和电分析化学、电子回路、生物系统的电传导等领 [1] 定义[2]分数阶导数有三种定义，分别是： Grunwald ? Letnik 定义 [2] (G ? L) 定义 、 Riemann ? Liouville 广。 R ? L 定义： ( R ? L) [1] 和 Caputo 定义 ，其中 R ? L 定义是 G ? L 定义的推 设 γ ∈ (0,1), a, b ∈ R, a t b, f (t) 在[a, b] 上连续，定义 R ? L 分数微分为下式： D1?γ f (t) = 1 d t f (τ ) dτ a t Γ(γ ) dt ∫ a (t ?τ )1?γ , 其中 Γ(?) 是 Gamma 函数。 Caputo 定义如下： Dγ f (t) = 1  tf (s)(t ? s)?γ ds, 0 ≤ t ≤ T , 0 γ 1. t a t Γ(1? γ ) ∫ a 在 γ 为负实数和正整数时三种定义可互相转换，G ? L 定义一般用于离散化计算；R ? L 定 义和 Caputo 定义常用于对分数阶微分方程的讨论。 1.2 应用举例 例 1 牛顿定律推广应用. Westerlund 等建议把牛顿第二定律 f  = ma 用 f  = mx(a ) , (1 a 2) ，取而代之，其中 x(a) = Da x 是 x 的 a 阶导数，具体参见文献[3-6]。 例 2 生物学传导应用. ?α专家在研究生物电传导实例时，给出了传递函数为： Χ(ω) = Χ0 ?α  , (0 α 1) ，其中 ω 是电流频率，Χ0 ,α 都是常数，其值视具体情况而定，与细胞种类相关。如果视上式为 Laplace 0变换式（L变换），即 Sα G(s) = G ，则 L 逆变换是分数微分方程 0 α 0 Dt g(t) = 0, 0 α 1. 例 3 神经控制应用. [7] 提出了人的眼球随着头的转动方式，并给出了分数阶的反射效应，模型如下： r (t ) ν +1 ν τ1D p P + r(t) = τ1τ 2 D u(t) +τ1D u(t), 其中 D =0 Dt 是 R-L 导数，r(t) 是放电率，u(t) 是头转动的速度，αd 是分数导数的阶， αi 是实际情况中的分数，ν = ad ? ai 。 例 4 熔炉加热问题应用. [8] 中讨论熔炉再加热问题时给出了以下一系列方程： 2 1 0a y (t) + a y (t) + a y(t) = u 2 1 0 (α ) ( β ) b2 y (t) + b1 y (t) + b0 y(t) = u(t), 1b y( β ) 1 (t) + b0 y(t) = u(t), 作者经过详尽的测量计算之后比较结果得出分数阶方程比整数阶方程效果要好。 例 5 电回路的应用. 文献 [1] 中推导出带有特定电阻器的系统方程如下： ? Dq x(t) = α  Dq?1 ( y(t) ? x(t)) ? 2α (4x(t) ? x3 (t)), ? 0 t 0 t 7 ? ? dy(t) x(t) ? y(t) + z(t), ?? dt = ? ? dz(t ) 100 y(t). ?? dt =? 7 ? 例 6 PI λ Dμ 控制器应用. PI λ Dμ 控制器，函数为 G(s) = U (s) = k + k s?λ + k sμ , 其中 λ, μ 0, k , k , k 为常 E(s) P I D P 1 D 数，输出方程为： I D Pk D?λ e