一類帶beddingon-angelis反應(yīng)項(xiàng)的捕食-食餌模型的共存態(tài)問題

一類帶beddingon-angelis反應(yīng)項(xiàng)的捕食-食餌模型的共存態(tài)問題

1模型1:非生物背景-本研究包括以下反應(yīng)擴(kuò)散系統(tǒng)。其中,Ω為Rn中的有界開區(qū)域,且邊界?Ω充分光滑;u,v分別表示食餌和捕食者的密度;參數(shù)a,b,d均為正常數(shù),m,k為非負(fù)常數(shù),c可正可負(fù),c>0表示捕食者有其他的食物來源;uv/(1+mu+kv)通常稱之為Beddington-DeAngelis(簡稱B-D)反應(yīng)函數(shù),它與Holling-Tanner反應(yīng)函數(shù)非常相似,只是分母中多了一項(xiàng)kv,這一項(xiàng)表示捕食者之間的相互干涉.同時(shí),B-D反應(yīng)函數(shù)保持了比率依賴反應(yīng)函數(shù)的所有優(yōu)點(diǎn),并避免了低密度問題引起的爭議,更好地反映了捕食者對食餌的捕食效應(yīng).模型(1)的生物背景和各參數(shù)的生物意義可參見.本文將考慮(1)對應(yīng)的平衡態(tài)問題正解的存在性和穩(wěn)定性.目前,對于帶B-D反應(yīng)項(xiàng)的捕食-食餌系統(tǒng)的研究更多關(guān)注的是常微分模型,對于偏微分模型的研究并不是很多.如果k,m=0,那么(2)即為經(jīng)典的Lotka-Volterra捕食-食餌模型,此類模型在過去的幾十年里已被廣泛研究.如果k=0且m>0,那么(2)即為帶Holling-Tanner反應(yīng)項(xiàng)的捕食-食餌模型.首次討論了該類模型,利用分歧理論給出了正解的存在性.進(jìn)一步研究了該類模型,詳細(xì)分析了當(dāng)m充分大時(shí)正解的存在性情況.本文主要利用分歧理論和線性特征值擾動理論研究了(2)正解的存在性和穩(wěn)定性.內(nèi)容安排如下:第2節(jié)利用Crandall-Rabinowitz分歧定理給出了局部分支解的存在性;第3節(jié)將局部分支延拓為整體分支,從而得到了正解存在的充分條件;第4節(jié)給出了局部分支解的穩(wěn)定性.2分布式/[a算子t10,x設(shè)q(x)∈C1(),考察如下特征值問題則(3)的所有特征值滿足:λ1(q)<λ2(q)≤λ3(q)≤…→∞,相應(yīng)的特征函數(shù)為φ1,φ2,φ3…,其主特征值是單重的.而且,有如下比較定理:若q1≥q2,則λj(q1)≥λj(q2),且若又有q1q2,則λj(q1)>λj(q2).記λ1=λ1(0),相應(yīng)的特征函數(shù)記為Φ1>0(x∈Ω).考慮非線性邊值問題由知,如果a≤λ1(q),那么u=0是(4)的惟一非負(fù)解;如果a>λ1(q),(4)有惟一正解.特別地,當(dāng)q(x))≡0時(shí),如果a>A1,那么非線性邊值問題有惟一正解,記為Θ[a],且Θ[a]關(guān)于a是單調(diào)遞增的,并有Θ[a]<a.顯然,當(dāng)a>λ1時(shí),(2)存在半平凡解(Θ[a],0);當(dāng)c>λ1時(shí),(2)存在半平凡解(0,Θ[c]).利用極大值原理很容易得到(2)非負(fù)解的一個(gè)先驗(yàn)估計(jì).引理1若(u,v)是(2)的一個(gè)非負(fù)解,則a>λ1,且令,固定c>λ1,取a為分歧參數(shù),從半平凡的非負(fù)解分支{(a;0,Θ[c])}出發(fā),構(gòu)造(2)的正解.設(shè)?為邊值問題的主特征值,相應(yīng)的主特征函數(shù)不妨設(shè)為Φ>0,且滿足||Φ||∞=1.令ω=u,X=則0≤ω≤Θ[a],X≥0,且ω,X滿足其中顯然F=(F1,F2)連續(xù),F(0,0)=0,且F關(guān)于(ω,X)的Frechet導(dǎo)數(shù)D(ω,x)F|(0,0)=0.令K為(-Δ)-1,則(6)等價(jià)于定義算子T:R+×X→X為則T(a;ω,X)為X上的緊可微算子.令G(a;ω,x)=(ω,X)′-T(a;ω,X),則G是C1函數(shù),且G(a;0,0)=0.易知,G(a;ω,X)滿足0≤ω≤Θ[a],X≥0的零點(diǎn)恰好是(2)的非負(fù)解.定理1設(shè)c>λ1,則(a;u,v)=(?;0,Θ[c])為(2)的分歧點(diǎn),且在(?;0,Θ[c])的鄰域內(nèi),(2)存在正解,其中?由(5)給出.證記G(a;ω,X)關(guān)于(ω,X)在(0,0)點(diǎn)的Frechet導(dǎo)數(shù)為L(a;0,0)=D(ω,X)G(a;0,0),則因此L(a;0,0)·(ω若ω≡0,則X≡0.因此ω≠0,那么由(5)知,a=?,ω=Φ.由于算子(-Δ-c+2Θ[c])的特征值均大于0,那么且由廣義極大值原理知Ψ>0.因此N(L(?;0,0))=span{(Φ,Ψ)}.另一方面,L(?;0,0)的伴隨算子因此L*(?;0,0)·(ω,X)=0等價(jià)于容易看出,X≡0.因此ω≡Φ,從而有N(L*(?;0,0))=span{(Φ,0)}.由Fredholm選擇公理,R(L(?;0,0))={(ω,X)∈X:∫ΩωΦdx=0}.因此codimR(L(?;0,0))=1.令,則L1(?;0,0))·(Φ,Ψ)=(-KΦ,0)′.假設(shè)存在(ω0,X0)∈X滿足即(ω0,X0)滿足在上式第一個(gè)方程兩端同乘以Φ,在Ω上積分,并利用分部積分得從而有∫ΩΦ2dx=0,與‖Φ‖∞=1矛盾.因此L1(?;0,0)·(Φ,Ψ)?R(L(?;0,0)).至此,中簡單特征值分歧定理的條件得到滿足,故存在充分小δ>0及C1連續(xù)曲線(a(s);φ(s),φ(s)):(-δ,δ)→R×Z,使得a(0)=?,φ(0)=0,φ(0)=0,且(a(s);ω(s),x(s))=(a(s);s(Φ+φ(s)),s(Ψ+φ(s)))滿足G(a(s);ω(s),X(s))=0,其中Z⊕span{(Φ,Ψ)}=X.因此(a(s);U(s),V(s))(|s|<δ)是(2)的分支解,其中U(s)=s(Φ+φ(s)),V(s)=Θ[c]+s(Ψ+φ(s)).若取0<s<δ,則它恰好為(2)的正解.而且,在分歧點(diǎn)(?;0,Θ[c])附近的非平凡非負(fù)解要么在分支{(a;0,Θ[c]):a∈R+}上,要么在分支{(a(s);U(s),V(s)):0<s<δ}上.設(shè)a>λ1固定,以c為分歧參數(shù),考察(2)在(Θ[a],0)附近的分支解,類似可得定理2設(shè)a>λ1,則(c;u,v)=(;Θ[a],0)為(2)的分歧點(diǎn),且在(;Θ[a],0)的鄰域內(nèi),(2）存在正解,其中.引理2設(shè)λ1<a<λ1+.若(u,v)是(2)的一個(gè)正解,則存在M>0,使得<c<M.證由于0<u<Θ[a],結(jié)合(2)中關(guān)于v的方程得令.由于λ1(qc)是關(guān)于c的單調(diào)遞增函數(shù),且有.因此選取充分大的M>0,如果c>M,那么有λ1(qc)>a.又由于當(dāng)c>λ1時(shí),Θ[c]≤v,結(jié)合(2)中關(guān)于u的方程得矛盾說明c不可能充分大,因此c≤M.注1引理2說明在半平凡解分支{(c;Θ[a],0)}上產(chǎn)生的局部分支解必在分歧點(diǎn)(;Θ[a],0)的右側(cè).同時(shí)也說明,如果該分支解能延拓為整體分支,那么該整體分支不可能隨參數(shù)c到達(dá)無窮遠(yuǎn)處.3有條件區(qū)分本節(jié)主要借助于中的方法(也可參見),將定理1和定理2給出的局部分支延拓為全局分支.令:u(x)>0,x∈Ω,<0,x∈?Ω}.定理3設(shè)c>λ1,則由定理1給出的正分支解在正錐P內(nèi)沿a增大到∞,這里P={(a;U,V)∈R+×X:U,V∈P1}.證令設(shè)μ≥1是T′(a)的一個(gè)特征值,則有假設(shè)ω≡0,由于算子(-μΔ-c+2Θ[c])的特征值均大于0,那么X≡0,因此ω0,于是存在某個(gè)i,使得a=ai(μ),其中ai(μ)是問題的特征值.顯然,ai(μ)關(guān)于μ在[1,+∞)上遞增,且ai(μ)滿足.另一方面,如果μ≥1,則算子(-μΔ-c+2Θ[c])的所有特征值都大于0,從而.因此μ≥1為T′(a)的特征值當(dāng)且僅當(dāng)存在某個(gè)i,使得a=ai(μ).假設(shè)a<?,則對任意μ≥1,i≥1,有a<a1(1)≤ai(μ).因此T′(a)沒有大于或等于1的特征值,即當(dāng)a<?時(shí),index(T(a;·),0)=1.假設(shè)?<a<a2(1),則對任意μ≥1,i>2,有a<ai(μ).又a1(1)=?,,且a1(μ)關(guān)于μ嚴(yán)格遞增,因此存在惟一的μ1>1,使得a=a1(μ1).于是,dimN(μ1I-T′(a))=1,其中>0為如下特征值問題的主特征函數(shù).下面證明R(μ1I-T′(a)∩N(μ1I-T′(a))=0.若不然,不失一般性,設(shè),那么存在(ω,X)∈X,使得,即在上述方程兩端同乘以,在Ω上積分,并利用格林公式得另一方面,結(jié)合(7)式可得,矛盾.因此μ1的代數(shù)重?cái)?shù)是1.故當(dāng)?<a<a2(1)時(shí),index(T(a;·),0)=-1.由全局分歧定理知,在R+×X內(nèi),存在從(?;0,Θ[c])出發(fā)的連通分支C0滿足G(a;ω,X)=0,且在點(diǎn)(?;0,Θ[c])附近,G(a;ω,X)的所有零點(diǎn)都在定理1給出的那條曲線上.令C1=C0-{(a(s);s(Φ+φ(s)),s(Ψ+φ(s))):-δ<s<0},C={(a;U,V):U=ω,V=Θ[c]+X,(a;ω,x)∈C1},則C為(2)由(?;0,Θ[c])出發(fā)的解曲線,在(?;0,Θ[c])的小鄰域內(nèi)有C?P,并且分支C-{(?;0,Θ[c])}滿足下列條件之一:(i)C連接了分歧點(diǎn)(?;0,Θ[c])和(a;0,Θ[c]),其中I-T′(a)不可逆,且?≠a.(ii)C在R+×X內(nèi)由(?;0,Θ[c])延伸到∞.(iii)C包含形如(a;u,Θ[c]+v),(a;-u,Θ[c]-v)的點(diǎn),其中(u,w)(0,0).下面證明C-{(?;0,Θ[c])}?P假設(shè)C-{(?;0,Θ[c])}?P,則存在點(diǎn)(a;?,)∈(C-{(?;0,Θ[c])})∩?P和序列{(an;un,vn)}?C∩P,un>0,vn>0,使得當(dāng)n→∞時(shí),(an;un,vn)→(a;?,).易知?∈?P1或∈?P1.假設(shè)?∈?P1,那么?≥0,x∈,則要么存在x0∈Ω使得?(x0)=0,要么存在x0∈?Ω,使得.由于?滿足方程由極大值原理得?三0.同理,假設(shè)∈?P1,則≡0.因此(?,v)有以下三種可能:(1)假設(shè),則當(dāng)n→∞時(shí),(an;un,vn)→(a;Θ,0).令,那么Vn滿足由LP估計(jì)和Sobolev嵌入定知,存在Vn的收斂子列(不失一般性,仍記為Vn)使得理當(dāng)n→∞時(shí),Vn→V在上成立,且V≥0,0,x∈Ω因此,在方程(8)中,令n→∞得由極大值原理得V>0,x∈Ω因此,與c>λ1矛盾.(2)假設(shè)(?,)≡(0,Θ[c]),則當(dāng)n→∞時(shí),(an;un,vn)→(a;0,Θ[c]).令那么Un滿足由LP估計(jì)和Sobolev嵌入定理知,存在Un的收斂子列(不失一般性,仍記為Un)使得當(dāng)n→∞時(shí),Un→U在上成立,且U≥0,0,x∈Ω因此,在方程(9)中,令n→∞得到由極大值原理得U>0,x∈Ω.因此,與a≠?矛盾.因此C-{(?;0,Θ[c])}?P.由引理1知,0≤U≤a,Θ[c]≤V≤.由LP估計(jì)和Sobolev嵌入定理,存在常數(shù)M>0,使得‖U‖,‖V‖≤M.因此C在正錐內(nèi)只能沿參數(shù)a延伸到∞.類似地,定理2給出的局部分支也可以延拓為全局分支,我們只須對分支的走向作進(jìn)一步說明.定理4設(shè)λ1<a<,則由定理2給出的正分支解與S1={(c;0,Θ[c])}相連接.證令S為(2)由(;Θ[a],0)出發(fā)的解曲線,在(;Θ[a],0)的小鄰域內(nèi)有S?P.那么由整體分歧定理知,分支S-(;Θ[a],0)滿足下列條件之一:(i)S連接了分歧點(diǎn)(;Θ[a],0)和(;Θ[a],0),且.(ii)S在R×X內(nèi)由(;Θ[a],0)延伸到∞.(iii)S包含形如(c;Θ[a]-u,v),(c;Θ[a]+u-v)的點(diǎn),其中(u,v)(0,0).假設(shè)S-(;Θ[a],0)?P,則上述條件(i)、(iii)均不成立.由引理1知‖u‖∞,‖v||∞≤M.由引理2知,c有界,那么(ii)也不成立.因此S-(;Θ[a],0)?P,從而存在點(diǎn)和序列{(cn;un,vn)}?S∩P使得當(dāng)n→∞時(shí),(cn;un,vn)→(;,).由極大值原理易知,(,)有以下三種可能:(1);(2);(3).假設(shè).令,則Vn滿足由LP估計(jì)和Sobolev嵌入定理知,存在Vn的收斂子列(不失一般性,仍記為Vn)使得當(dāng)n→∞時(shí),Vn→V在上成立,且V≥0,0,x∈Ω.因此,在方程(10)中,令n→∞得,V|?Ω=0.由極大值原理得V>0,x∈Ω.因此,與矛盾.假設(shè),類似于上面的方法同樣得出矛盾.因此,即S與S1相連接于點(diǎn)(;0,Θ).注2事實(shí)上(;0,Θ)是S1={(c;0,Θ[c])}上的一個(gè)分歧點(diǎn).由于(?;0,Θ[c])是半平凡解分支{(a;0,Θ[c])}上的分歧點(diǎn),那么恰由來決定.由以上分析我們可以得到(2)正解存在的一個(gè)充分條件:定理5設(shè)λ1<a<.若<c<,則(2)至少存在一個(gè)正解.4特征值aa下面討論該分支解的穩(wěn)定性.令,.其中0<α<1.i:X1→Y為由X1到Y(jié)的包含映射.設(shè)L0是(2)在(?;0,Θ[c])處的線性化算子,由定理1的證明可知,N(L0)=span{(Φ,Ψ)},coLdimR(L0)=1,且R(L0)={(u,v)∈X:∫ΩuΦdx=0}.又因?yàn)閕(Φ,Ψ)?R(L0),所以0是L0的i-單重特征值.引理30是L0實(shí)部最大的特征值,其它的特征值均在左半復(fù)平面上.證假設(shè)λ0是L0實(shí)部大于0的特征值,(φ,ψ)′為相對應(yīng)的特征函數(shù),則L0·(φ,ψ)=λ0(φ,ψ)′,即假設(shè)φ≡0,則λ0是算子(Δ+(c-2Θ[c])I)的一個(gè)特征值,那么λ0∈R,且λ0<0,與假設(shè)矛盾.因此φ0,那么λ0是算子的一個(gè)特征值.因?yàn)?所以0是算子的主特征值,從而有λ0≤0.同樣與假設(shè)條件矛盾,故假設(shè)不成立.因此0是L0實(shí)部最大的特征值,其它特征值均在左半復(fù)平面上.L(a;0,Θ[c])S(a)=γ(a)S(a),其中|a-?|<<1;L(a(s);U(s),V(s))T(s)=η(s)T(s),其中0<S<<1,這里S(a)=(φ1(a),φ2(a))′,T(s)=(ψ1(s),ψ2(s))′.而且γ′(?)≠0.又若η(S)≠0(|s|<<1),則這里a′(s)為a(s)關(guān)于s的導(dǎo)數(shù),γ′(?)為γ(a)關(guān)于a在a=?處的導(dǎo)數(shù).分支解(U(s),V(s))的穩(wěn)定性由η(s)的符號決定.當(dāng)η(s)<0時(shí),分支解是穩(wěn)定的;當(dāng)η(s)>0時(shí),分支解是不穩(wěn)定的.而η(s)的符號與sa′(s)γ′(?)的符號相反,因此可以通過判斷sa′(s)γ′(?)的符號來判斷分支解(U(s),V(s))的穩(wěn)定性.引理4γ′(?)>0.證由L(a;0,Θ[c])S(a)=γ(a)S(a)可知因?yàn)閨a-?|<<1,所以|γ(a)