-
近几十年来,用于描述生物种群动力学行为的偏微分方程模型引起了许多研究者的兴趣,其中两个重要的建模形式为以年龄结构建模以及大小结构建模,而对所建立的模型最值得关注的问题之一就是如何证明方程(组)的(弱)解的存在唯一性,其中称之为半群理论的一种方法是研究自治系统适定性和通过逼近得到收敛性的常用手段[1-2],除此之外还有其他几种,例如不动点定理[3-5]、差分格式[6]以及比较原理[7-9]. 其中,差分格式与不动点定理则可用于增长率依赖于种群总量的非线性非自治系统中,比较原理一般用于增长率与种群总量无关或常微分与偏微分耦合的系统[10]. 不动点定理的运用很广泛,例如:文献[3]利用不动点定理等方法,给出了线性矩阵不等式条件的反应扩散马尔科夫跳跃周期模糊时滞系统的随机稳定性判据,并通过建立在乘积空间上的压缩映射克服了反应扩散模型带来的数学上的困难;文献[5]利用不动点定理证明了具有非线性自然增长率的大小结构模型的偏微分方程的解的存在唯一性;文献[4]用不动点定理证明了n个种群耦合的大小结构模型方程组的解的存在唯一性并分析了种群的共存与排斥;文献[6]利用差分格式证明了弱解的存在唯一性,值得一提的是,文献[4-5]仅考虑了具有单一阶段的物种模型,而在现实世界中,有很多物种是具有两个阶段的,例如蝌蚪与青蛙,蚕与蚕蛾. 因此对于此类物种,模型的建立与文献[4-5]有所不同. 文献[9]将种群分成了幼年与成年两个阶段,并分别用J(a,t),A(x,t)来表示,其中幼年阶段J(a,t)用年龄结构进行建模,成年阶段A(x,t)用大小结构进行建模. 但值得注意的是,文献[9]利用比较原理以及单调有界定理证明解的存在唯一性时,非线性增长率g通常无法依赖φ(t),而本文主要参考文献[9]中的模型,即仍然以年龄结构和大小结构分别对物种的两个阶段(幼年与成年)进行建模,并假设幼年阶段的死亡率v,成年阶段的增长率g,死亡率μ,繁殖率β均依赖于种群的总数量φ(t):
再利用特征线法和不动点定理证明模型的解的存在唯一性. 显然,对每个参数作如此假设是符合现实意义的,原因是在资源有限(例如有限的光照、有限的水资源、有限的食物等等)的环境中,会同时存在种内竞争与种间竞争.
全文HTML
-
给出模型如下:
先对模型的参数作一些假设.
(H1) 对∀φ∈[0,∞),g(φ)>0且连续可微.
(H2) 对∀φ∈[0,∞),v(φ)>0且连续可微.
(H3) 对∀φ∈[0,∞),μ(φ)>0且连续可微.
(H4) 对∀φ∈[0,∞),β(φ)>0连续可微且有界,即β(φ(t))≤βM.
(H5) J(amax,t)(>0)∈L∞[0,T].
(H6) J0(a)(>0)∈L1[0,amax],A0(x)(>0)∈L1[xmin,xmax].
我们可以通过特征线法[4-5]求出模型(1)J(a,t),A(x,t)的解析表达式为
其中J,A的特征线由下列常微分方程组给出
由方程组(6)可解出a=t+c,进而得到(2)式与(3)式,对于方程组(6)第2式由(H1)中关于g>0的假设可知,x关于t是严格单增的,从而方程组(6)第2式对于任意初值(t0,x0)有唯一解X(t;t0,x0),且对应有唯一反函数Γ(x;x0,t0),而(t,X(t;0,xmin))则表示一条过点(0,xmin)且将平面分为两部分的特征线,因此结合初边值条件可得到(4)式与(5)式. 再对(2)-(3)式与(4)-(5)式分别关于a与x进行积分,得到如下关于φ(t)的表达式:
在(7)式的右端作两次变量代换:x→η(X(t;η,xmin)=x,即η=Γ(xmin;x,t)),以及x→ξ:=X(0;t,x). 因此,如果φ(t)是(7)式的非负连续解,则(1)式有连续解.
-
由于已经建立(1)式与(7)式的对应关系,因此要研究原模型解的存在唯一性,只需要关注问题(7)的可解性. 首先令K>‖J0‖L1+‖A0‖L1=
$\int_0^{{a_{\max }}} {{J_0}} \left( a \right){\rm{d}}a + \int_{{x_{\min }}}^{{x_{\max }}} {{A_0}\left( x \right)} {\rm{d}}x$ ,ST,K= {f∈C[0,T]|f(0) =‖J0‖L1+‖A0‖L1,0≤f(x)≤K },再定义算子O:ST,K→ST,K,由于ST,K是Banach空间C[0,T]的闭子空间,故也是完备的,因此仍可取范数为C[0,T]的范数,最后定义(Oφ)(t) 为等式(7)右端函数组.定理 1 假设(H1)-(H6)成立,则存在充分小的T0满足T>T0>0,使得问题(1)在[0,T0]存在唯一解.
证 首先令(Oφ)(t)∈ST,K. 由于对任意的φ∈ST,K,有
当T充分小时,总可以使得(Oφ)(t)≤K. 因此O:ST,K→ST,K.
注意到(4)-(5)式与(7)式的后两项,得到:
为了简便可记
再证明算子O是压缩映射. 对任意的φ1,φ2∈ST,K. 此时对应有B1,B2,则有
注意到
其中βK=supφ∈[0,K]|β′(φ)|,μK=supφ∈[0,K]|μ′(φ)|. 从而
其中
当T充分小时(即在[0,T0]内),可得出O是压缩的,从而存在唯一的不动点φ*使得Oφ*=φ*,此即说明(7)式有唯一解,故而(1)式有唯一解,证毕.
进一步,为了建立问题(1)在[0,T]上的解的存在唯一性,我们先证明φ(t)有上界. 由于
由Gronwall不等式可知
故而由文献[5]中定理3的证明可知问题(1)在[0,T]内有唯一解.
-
本文中,我们利用不动点定理证明了模型(1)的解的存在唯一性,值得注意的是,我们仅考虑了各种参数依赖于种群总量φ(t),后续值得进一步研究的问题如下:①可将参数设为依赖大小x,年龄a,时间t来建立模型从而考虑用差分方式来证明解的存在唯一性;②选取适当的函数应用到死亡率、增长率、繁殖率中再通过有限差分法进行数值模拟,清晰刻画种群动力学变化,即种群的持久与灭绝;③现实中,种群所在的环境大多存在污染物,可建立有关模型考虑不同的污染物浓度是怎样影响种群的动力学行为的.