-
具有交错扩散的种群模型是一类重要的生物模型,此类模型最初由Shigesada等人提出[1],主要用于描述种群栖息地的分化现象.此后,有很多生物数学工作者对此类模型的动力学行为、平衡态模式进行了研究[2-7].最近,文献[8]研究了如下一类带有分数型交错扩散项的捕食-食饵模型:
其中:
$ \varOmega \subset {\mathbb{R}^N}\left( {1 \le N \le 5} \right)$ 是有界区域并具有光滑边界$ \partial \varOmega $ ;系数a,c,d,r2>0,α,β≥0,r1≤1,$b \in \mathbb{R} $ 是实常数.模型(1) 的相关生物背景可参见文献[8].在该文献中作者对此模型正解的存在性条件及其动力学行为进行了研究,所谓正解是指满足模型(1) 并在Ω内取正的解.下面我们介绍一些本文所需的符号和结果.令λ1(q)表示如下特征值问题的主特征值:
其中
$q\left( x \right) \in C\left( {\overline \Omega } \right) $ ,由文献[8]可知λ1(q)是简单的,其对应的特征函数在Ω内不变号,并且如下结论成立:若(q)1(x)≤(q)2(x),$ {q_1}\left( x \right)\not \equiv {q_2}\left( x \right)$ ,则(λ)1((q)1)<(λ)1((q)2).我们记(λ)1=(λ)1(0),并令${\phi _1} $ 为其对应的满足${\left\| {{\phi _1}} \right\|_{{L^2}\left( \varOmega \right)}} = 1 $ 的正特征函数.类似于文献[8],当常数a>(λ)1时,我们记如下Logistic方程
的唯一正解为θa.由文献[9]可知:如果a≤(λ)1,u=0是(3) 式的唯一非负解;θa是关于a的递增函数且对任意的x∈Ω有0<θa<a;
$ \mathop {\lim }\limits_{a \to \lambda _1^ + } {\theta _a}\left( x \right) = 0$ 是一致的;$\mathop {\lim }\limits_{a \to \infty } {\theta _a} = \infty $ 和$\mathop {\lim }\limits_{a \to \infty } \frac{{{\theta _a}}}{a} = 1 $ 在Ω的任意紧子集上一致成立.设E是一个Banach空间,W是E的一个闭凸子集,
$ \mathit{\boldsymbol{A}}:W \to W$ 是一个紧的Frechet可微的算子,y∈W是A的一个孤立不动点,G是W的一个开子集且A在G的边界$ \partial G$ 上没有不动点,类似于文献[5],我们可以定义Leray-Schauder度degW(I-A,G,0).孤立不动点y的不动点指数index(A,y)定义为degW(I-A,N(y),0),其中N(y)表示y在W中的一个小邻域.进一步,如果degW(I-A,G,0)≠0,则有A在G中至少有一个不动点.近年来,有许多作者研究了具有扩散项的种群模型的正解的多重性,如文献[10-14]等.文献[8]没有考虑模型(1) 正解的多重性.在本文中,我们将利用分支理论和拓扑度理论研究模型(1) 正解的多重性条件.
全文HTML
-
由文献[8]可知,模型(1) 等价于如下问题:
其中,
显然模型(4) 具有平凡解
$\left( {\overline U, \overline V } \right) = \left( {0, 0} \right) $ ;如果a>(λ)1,模型(4) 具有半平凡解$\left( {\overline U, \overline V } \right) = \left( {{{\left( {2 - {r_1}} \right)}^2}{\theta _a}, 0} \right) $ ;如果b>(λ)1,模型(4) 具有半平凡解$\left( {\overline U, \overline V } \right) = \left( {0, {{\left( {1 + {r_2}} \right)}^2}{\theta _b}} \right) $ .由此可以看出模型(1) 有两条半平凡解曲线:为了研究以上各平凡解和半平凡解的不动点指数,我们首先定义以下集合:
其中,
由文献[8]中引理4.1可知(4) 式的非负解都在D内,我们记D的内部为
$ \overset{{}^\circ }{\mathop{D}}\, $ .定义E上的算子A为
其中:f(u,v)=u((2-r1)a-u-cv),g(u,v)=v((1+r2)b+du-v),
$ u\left( {\overline U, \overline V } \right)$ 和$v\left( {\overline U, \overline V } \right) $ 由(5) 式给出.常数p充分大使得$\mathit{\boldsymbol{A}}:\bar D \to W $ .我们令根据文献[8]中的引理5.1和引理5.2,如下结论成立:
(¡)
${{\deg }_{w}}\left( \mathit{\boldsymbol{I}}-\mathit{\boldsymbol{A}}, \overset{{}^\circ }{\mathop{D}}\, \right)=1 $ ,(¡¡)如果a>(λ)1且b≠(λ)1,则有
(¡¡¡)如果a>(λ)1,b≠(λ)1且(λ)1(φ(θa,b))>0,则
(¡v)如果a>(λ)1,b>(λ)1且(λ)1(ψ(a,θb))<0,则
为了给出本文的主要结果,我们在(a, b)平面内定义两条曲线S1和S2如下:
由文献[8]中的引理6.1和引理6.2可知:
其中a*是一个大于(λ)1的正常数且满足
$f\left( {{a}^{*}} \right)=-\frac{d}{\beta {{r}_{2}}}, f\left( a \right) $ 是[(λ)1,a*]上关于a递减的连续可微函数,g(a)是[(λ)1,+∞)上关于a递增的连续可微函数,并且f(a)和g(a)满足利用分支理论,我们可以得到如下定理:
定理1 固定a>(λ)1,模型(4) 曲线在Γ1附近有分支解当且仅当b=f(a),并且对充分小的δ>0,模型(4) 在((2-r1)2θa,0,f(a))附近的分支正解
${{\varGamma }_{3}}=\left\{ \left( \overline{U}, \overline{V}, b \right)=\left( \overline{U}\left( s \right), \overline{V}\left( s \right), b\left( s \right) \right):0<s\le \delta \right\} $ 可表示为:并且有
其中ξ和η的定义见(18) 和(19) 式,{Φ1(s),Ψ1(s),b(s)}关于s是光滑的并满足
证 设
$X={{\left[{{W}^{2, p}}\left( \varOmega \right)\bigcap W_{0}^{1, p}\left( \varOmega \right) \right]}^{2}}, Y={{\left[{{L}^{p}}\left( \varOmega \right) \right]}^{2}} $ ,定义映射$\mathit{\boldsymbol{S}}:X\times \mathbb{R}\to Y $ 如下:其中f和g由(8) 式给出.通过计算,我们可以得到:
由于b=f(a)时,(λ)1(φ(θa,b))=0,令ξ是其对应的满足
${{\left\| \xi \right\|}_{{{L}^{2}}\left( \varOmega \right)}}=1 $ 的正特征函数,即ξ满足方程由于算子-Δ+2θa-a在Dirichlet边界条件下是可逆的,我们令η是问题
的唯一解.经计算得
在
$ \left( \overline{U}, \overline{V}, b \right)=\left( {{\left( 2-{{r}_{1}} \right)}^{2}}{{\theta }_{a}}, 0, f \right)$ 时,我们可以得到此外,由于
$\int{_{\Omega }\left( 1+{{r}_{2}} \right)\frac{1+\beta \left( 2-{{r}_{1}} \right){{\theta }_{a}}}{1+{{r}_{2}}+\beta \left( 2-{{r}_{1}} \right){{\theta }_{a}}}{{\xi }^{2}}\text{d}}x>0 $ ,可知从而由Rabinowitz局部分支理论[15-16]可知定理成立.再根据文献[17]可知,
其中:A,B由(15) 式定义,L是Y上的线性泛函满足〈L,(χ,ψ)〉=∫Ωψξdx.定理证明完毕.
显然我们可以知道
$ {{\left\{ \left( \overline{U}\left( s \right), \overline{V}\left( s \right), b\left( s \right) \right) \right\}}_{0<s<\delta }}$ 是(4) 式的正解,下面我们研究这组正解$\left( \overline{U}\left( s \right), \overline{V}\left( s \right) \right) $ 的稳定性.定理2 假设a>(λ)1且b1≠0,则存在常数δ′∈(0,δ],当0<s<δ′时,解
$\left( \overline{U}\left( s \right), \overline{V}\left( s \right), b\left( s \right) \right) $ 是非退化的.此外,若b1<0,则$\left( \overline{U}\left( s \right), \overline{V}\left( s \right), b\left( s \right) \right)$ 是不稳定的;若b1>0,则$ \left( \overline{U}\left( s \right), \overline{V}\left( s \right), b\left( s \right) \right)$ 是稳定的.证 设X,Y和S是定理1的证明中所定义的集合和映射,则对任意的(φ,ψ)∈X,有
考虑特征值问题
其中
当s→0时,
$\left( \overline{U}\left( s \right), \overline{V}\left( s \right), b\left( s \right) \right)\to {{\left( \left( 2-{{r}_{1}} \right) \right)}^{2}}{{\theta }_{a}}, 0, f $ .因此,由于
故算子
$ {{\mathscr{L}}_{0}}$ 有最小特征值0,并且$ {{\mathscr{L}}_{0}}$ 的其他特征值都大于0.由摄动定理知,存在δ1∈(0,δ]使得当s∈(0,δ1)时,$\mathscr{L} $ 存在唯一的特征值ω(s)满足并且
$\mathscr{L} $ 的其他特征值的实部都大于0.下面讨论当s>0适当小时ω(s)的符号,考虑特征值问题
类似
$\mathscr{L} $ 的讨论,可得存在一个适当小的常数δ2∈(0,δ]使得当|b-f|<δ2时,$-{{\mathit{\boldsymbol{S}}}_{\left( \bar{U},\bar{V} \right)}}\left( {{\left( 2-{{r}_{1}} \right)}^{2}}{{\theta }_{a}},0,b \right)$ 存在唯一的特征值m(b)满足并且
$-{{\mathrm{S}}_{\left( \overline{U}, \overline{V} \right)}}\left( {{\left( 2-{{r}_{1}} \right)}^{2}}{{\theta }_{a}}, 0, b \right)$ 的其他特征值都有正实部.另外,m(b)满足线性化特征值问题由于m(f)=0,ξ(f)=ξ,将(21) 式两边关于b求导,
令ζ=ξ′(f)并取b=f得
将(22) 式两边乘以ξ并在Ω上积分得
由(23) 式知,
因为b1≠0,所以存在一个正数δ3≤min{δ1,δ2}使得当0<s<δ3时ω(s)≠0并且有
此时可知存在一个正数
$\tilde{\delta }\le {{\tilde{\delta }}_{3}} $ 使得当0$0<s<\tilde{\delta } $ 时,ω(s)与b1符号相同.定理得证.基于上述讨论,模型(4) 正解的多重性结果如下:
定理3 固定a>(λ)1,则如下列结论成立:如果b1<0,则存在一个常数ε=ε(b)<(λ)1-f使得
1) 当f-ε<b<f时,模型(4) 至少有两个正解;
2) 当b≥f-ε时,模型(4) 至少有一个正解.
其中f和b1的定义可参见(12) 和(14) 式.
证 由定理1知,方程(4) 在((2-r1)2θa,0,b)附近有正解曲线Γ3.由于b1<0,则当s>0适当小时b(s)<f.
下面用反证法证明此定理.假设在f附近且b<f时方程(4) 有唯一解
$\left( \hat{U}, \hat{V} \right) $ ,显然$ \left( \hat{U}, \hat{V} \right)=\left( \overline{U}\left( s \right), \overline{V}\left( s \right) \right)$ ,并且由定理2知它是非退化的正解.因而$ \mathit{\boldsymbol{I}}-{{\mathit{\boldsymbol{A}}}_{\left( \hat{U},\hat{V} \right)}}\left( \hat{U},\hat{V} \right):\overline{{{W}_{\left( \hat{U},\hat{V} \right)}}}\to \overline{{{W}_{\left( \hat{U},\hat{V} \right)}}}$ 是可逆的.因为$ \left( \hat{U}, \hat{V} \right)\in \overset{{}^\circ }{\mathop{D}}\, $ ,根据文献[19]中定理4.9.3知$\rm{inde}{{\rm{x}}_{w}}\left( \mathit{\boldsymbol{A}},\left( \hat{U},\hat{V} \right) \right)=\pm 1 $ .注意到当$0<s\ll 1, a>{{\lambda }_{1}} $ 时,b<f.因此由(10) 式可知上面的等式不可能成立.故假设错误,定理得证.
-
本文研究了文献[8]中提出的一类具有分数型交错扩散项的捕食-食饵模型.首先利用分支理论研究了模型在半平凡解曲线附近的分支解的存在条件和稳定性条件,其次利用以上结论并结合拓扑度理论得到了模型正解的多重性条件.此结论推广并完善了文献[8]中的相关结果.