近年来, 随机分数阶微分方程引起了国内外学者的关注, 并被广泛的应用在物理、生物、工程、金融等领域[1-5].然而, 分数阶非线性系统还可能存在与整数阶非线性系统相似的非线性现象, 如自振、多平衡状态、混沌以及更复杂的过渡过程等[6, 7].
另一方面, 分数阶微分方程与分形密切相关, 而分形富含于生物系统中.研究表明, 在传统的整数阶微分方程不能建模的现象中, 分数阶微分方程为此提供了可能性.这里, 特别强调的是, 分数阶与整数阶模型最大区别在于分数阶模型拥有记忆, 而死亡率有长程相似性, 它的主要特征恰恰包含了记忆.因此, 本文首次尝试建立与年龄相关的随机分数阶动力学模型.另外, 已经证实由分数阶微分方程建立的某些生物学模型比整数阶更有优势[8, 9].
现在将分数阶引入与年龄相关的随机种群模型[10].新系统可描述为如下分数阶非线性系统:
其中 $0 < \alpha < 1$, $a\in [0, A]$表示年龄, $t\in [0, T]$表示时间, $0 < T < \infty$, $P_t:=P(t, a)$、 $\beta(t, a)$和 $\mu(t, a)$分别表示时刻 $t$年龄为 $a$的种群密度、生育率和死亡率. $f(t, P)$为外部环境对种群系统的影响, 如迁移、地震等突发性灾害对种群引起的影响, $g(t, P)\frac{d \omega_t}{dt}$为随机外界环境对系统的扰动.
在 $\alpha=1$的情况下, 方程(1.1) 成为经典的与年龄相关的随机种群模型, 该模型通过随机分析理论和数值方法已被广泛研究.例如张启敏等[10]研究了与年龄相关的随机种群方程解的存在性, 唯一性和指数稳定性.李荣华等[11]考虑了带Markovian跳的与年龄相关的随机种群方程数值解的收敛性.马维军和张启敏等[12, 13]讨论了与年龄相关的随机种群方程带分数布朗运动的数值解和Markovian跳的渐近稳定性.杨洪福等在文献[14]中介绍了与年龄相关的随机两种群方程在POD基下的数值解.
本文首次将分数阶引入与年龄相关的随机种群动力学模型, 讨论在满足一定的假设条件下证明了温和解的存在性、唯一性, 得到的结论是文献[10]的推广.
令
$V$是一个Sobolev空间, $\{\omega_t\}_{t\in J}$是定义在完备的概率空间 $(\Omega, \mathcal{A}, \mathbb{P})$取值于可分的Hilbert空间 $K$上的Brown运动,
下面, 首先给出两个相关的基本定义.
定义2.1[4] 函数 $f$的 $\alpha$分数阶积分定义为
其中 $t\geq t_0$且 $\alpha>0$, $\Gamma(\cdot)$为Gamma函数.
定义2.2[15, 16] 函数 $f$的 $\alpha$阶Caputo分数阶微分定义为
其中 $t\geq t_0$, $n$是一个正整数满足 $n-1 < \alpha < n$.特殊地, 当 $0 < \alpha < 1$时,
把与年龄相关的随机分数阶种群方程(1.1) 写成如下形式的积分方程[16]
其中
$S(t)$是由一个 $V$中的有界线性算子 $B=\frac{\partial}{\partial a}$生成的 $C_0$ -半群, $\xi_{\alpha}$是定义在 $(0, \infty)$的概率密度函数, 即 $\xi_{\alpha}(\theta)\geq 0$, $\theta\in(0, \infty)$并且
定义2.3[17] 如果一个 $\cal F$ $_t$ -循序可测的随机变量 $\{P(t, a)\}_{t\in J}$是方程(1.1) 的一个温和解, 则 $\{P(t, a)\}_{t\in J}$满足相应的积分方程(2.1).
引理2.4[18] 如果算子 $\{\mathscr{T}(t)\}_{t\in J}$和 $\{\mathscr{S}(t)\}_{t\in J}$是强连续的, 则对每个 $P_t \in H$, 当 $t_1\rightarrow t_2$, 有
为了证明本文的主要结论, 给出以下假设条件:
(ⅰ) $\mu(t, a)$非负可测, 并且
(ⅱ) $\forall t\geq 0$, $\{\mathscr{T}(t)\}_{t\in J}$和 $\{\mathscr{S}(t)\}_{t\in J}$是有界线性算子, 即
并且满足 $E\displaystyle\int_{0}^{T}\parallel P_t\parallel^2dt < \infty$, 则存在一个常数 $\gamma>0$, 使得
(ⅲ)设 $f(t, \cdot ) : L^2_H\rightarrow H$是对几乎所有 $t$有定义的非线性算子, 满足
(a.1) $f(t, 0)=0$;
(a.2) $\exists N>0$, 使得
设 $g(t, \cdot ) : L^2_H\rightarrow $ $\cal{L}$ $ (K, H)$是对几乎所有 $t$有定义的非线性算子, $g(t, x)\in \cal{L}$ $ (K, H)$且满足
(b.1) $g(t, 0)=0$;
(b.2) $\exists L>0$, 使得
在本节中, 利用第二部分给出的定义和引理, 证明与年龄相关的随机分数阶种群方程(1.1) 的温和解的存在性、唯一性.定义算子 $\Psi: H\rightarrow H$, 满足如下方程
为了证明文章的主要结论, 先证明如下引理.
引理3.1 对任意的 $P_t\in H$, $\Psi P_t$在 $[0, T]$上是 $L^2$ -连续的.
证 令 $0 < t_1 < t_2 < T$.对于方程(3.1)的任意固定点 $P_t\in H$, 有
由 $\mathscr{T}(t)$是强连续的, 则当 $t_2-t_1\rightarrow 0$时, 方程(3.2) 右侧第一项也趋于零.接下来, 应用Holder不等式和假设条件(ⅰ)-(ⅲ), 有
进一步, 有
类似的, 应用Holder不等式、Itô积分和假设条件(ⅰ)-(ⅲ), 可得
因此由 $\mathscr{S}(t)$是强连续的和Lebesgue控制收敛定理, 得出这样的结论:当 $t_2-t_1 \rightarrow 0$时, 方程(3.2)右侧不等式趋于零.因此得出在 $[0, T)$上 $\Psi P_t$是右连续.
类似的可以证明在 $[0, T)$上 $\Psi P_t$是左连续.因此引理得证.
定理3.2 如果条件(ⅰ)-(ⅲ)是成立的, 并且满足下面的不等式, 则方程(1.1) 存在唯一的温和解 $ P_t\in H$.
证 应用压缩映射原理证明算子 $\Psi$存在一个不动点.首先, 要证 $\Psi(H)\subset H$.令 $P_t\in H$, 由(3.1)式, 可得
现在估计方程(3.3) 右边的项.第一项, 由条件(ⅱ)可得
第二项应用Holder不等式和假设条件(ⅰ)-(ⅲ)有
进一步又可以得到
类似的, 应用Itô积分有
所以从(3.4) 式到(3.7)式可以得到 $E\parallel \Psi P_{t}\parallel^2 < \infty$.由引理3.1知 $\Psi P_{t}\in H$.因此算子 $\Psi$是从 $H$到 $H$的自映射.接下来, 将证明 $\Psi $是 $H$上的连续映射.事实上, 对任意的 $P_{1t}, P_{2t} \in H$, 由方程(3.1) 和Itô等距公式, 可得
现在估计方程(3.8)右边的项.第一项由Holder不等式和假设条件(ⅰ)-(ⅲ), 有
类似的, 第二项为
接下来, 对方程(3.7) 的第三项应用Itô积分可得
因此 $\forall t \in J$由方程(3.9)-(3.11)可以得到
即
从方程(3.12) 可以得出 $\Psi$是连续映射.由Banach压缩映射原理, 在 $H$上存在唯一不动点 $P_{t}$, 使得 $\Psi P_{t}=P_{t}$.因此可以得到
即 $P_{t}\in H$是方程(1. 1)的温和解.定理得证.
本文用分数阶导数代替古典的随机种群系统中关于时间的整数阶导数, 得到了一类随机分数阶种群模型, 该模型较传统的随机种群系统能更好地刻画现实中的种群问题, 尤其对具有遗传性的种群现象.在引入分数阶概念的基础上, 利用不动点定理、随机分析和算子半群理论, 讨论了随机分数阶种群系统温和解的存在性、唯一性.所得到的结论为种群未来的研究提供了一定的理论依据.