葉凱莉,王秋燕 ，宋 強

(1.信陽師范學(xué)院 a.工商管理學(xué)院;b.數(shù)學(xué)與信息科學(xué)學(xué)院, 河南 信陽 464000;2.信陽職業(yè)技術(shù)學(xué)院 數(shù)學(xué)與計算機科學(xué)系， 河南 信陽 464000)

0 引言

伴隨著工業(yè)化進程的推進和社會經(jīng)濟的發(fā)展,城市化現(xiàn)象越來越受到人們的關(guān)注.城市化在帶動區(qū)域經(jīng)濟發(fā)展的同時，也引起了城市和農(nóng)村人口結(jié)構(gòu)的變化.而城市外來人口的增加和農(nóng)村勞動力的減少對社會也造成了一定的消極影響，因此，控制城市化進程中的人口流動是必要的.

(1)

其中：x,y分別代表城市和鄉(xiāng)村人口密度；a,b,c,d為正參數(shù).這里假設(shè)x(0)>0,y(0)>0.顯而易見，該模型有非平凡平衡點.

本文對系統(tǒng)(1)進行了全局性分析.利用類似于文獻[5]的方法,首先給出了解的有界性，然后討論了平衡點的存在性和穩(wěn)定性.特別地，對正平衡點的動力學(xué)性質(zhì)進行了系統(tǒng)地研究，得到了其存在性和局部漸近穩(wěn)定性的條件.然后利用Bendixon-Dulac定理[6]得出了系統(tǒng)不存在非平凡正周期解，進而得到正平衡點是全局漸近穩(wěn)定的.最后對所得結(jié)論進行了討論和總結(jié).

1 解的有界性

根據(jù)解的存在唯一性定理和簡單的討論可知，系統(tǒng)(1)的解總是存在且為正.事實上，由系統(tǒng)(1)，易得

這分別意味著

(2)

首先考慮系統(tǒng)(1)中x(t)的有界性.顯然，根據(jù)解的正性，下面假設(shè)c>d.于是，有以下結(jié)論：

定理1 對于系統(tǒng)(1)，如果c>d，那么

(3)

證明由以上的討論可知，存在一個T>0，使得當(dāng)t>T時，y(t)<1.由系統(tǒng)(1)的第二個方程可知，對于t>T，有

于是，得到

這說明，對于c>d,有

通過以上討論，可以得到下面的結(jié)論.

定理2 如果c>d,那么系統(tǒng)(1)是持久的.

2 平衡點分析

2.1 平衡點的存在性

系統(tǒng)(1)有平衡點E0(0,0)和邊界平衡點E1(1,0)與E2(0,1).

假設(shè)存在正平衡點E*(x*,y*),則有

(4)

結(jié)合方程組(4),可得

-(abc+a2d)(x*)2+2ad(a+b)x*+

b2c+abc-a2d-b2d-2abd=0，

(5)

(6)

定義

F(x)=-(abc+a2d)x2+2ad(a+b)x+

b2c+abc-a2d-b2d-2abd,

(7)

(8)

那么(x*,y*)是系統(tǒng)(2)的正平衡點當(dāng)且僅當(dāng)x*是F(x)=0的一個正解和y*是正的，其滿足方程(8).

注意到方程(5)是一個二次方程，通常來講，其最多有兩個正根.但是，下面將說明多項式(7)僅有一個正根滿足式(6).

由式(7)可得，如果Δ=ab2c(a+b)(c-d)≥0,那么方程F(x)=0有正根.

當(dāng)Δ=0時，F(xiàn)(x)=0只有一個正根x=1，此時y=0，顯然這是一個邊界平衡點.

當(dāng)Δ>0即c>d時，F(xiàn)(x)=0有兩個正根x1和x2，其中x1>x2>0.由式(8)可得y1>0，y2<0.因此，根據(jù)解的正性，僅有一根x*=x1滿足式(6)，也就是說，E*(x*,y*)是系統(tǒng)(1)唯一的正平衡點.

綜上所述，對于系統(tǒng)(1)正平衡點的存在性，有下面的定理.

定理3 如果c≤d,則系統(tǒng)(1)存在一個平衡點E0(0,0)和兩個邊界平衡點E1(1,0)，E2(0,1)；不存在正平衡點.如果c>d，則系統(tǒng)(1)存在唯一的正平衡點E*(x*,y*).

2.2 平衡點的穩(wěn)定性

平衡點E0,E1,E2和E*的穩(wěn)定性分別是由系統(tǒng)在這些點處的雅克比矩陣J(E0),J(E1),J(E2)和J(E*)的特征值所決定的.易得，系統(tǒng)(2)的雅克比矩陣為

由系統(tǒng)(1)相應(yīng)的線性系統(tǒng)易得，原點(0,0)處的雅克比矩陣為

顯然，無論正平衡點E*是否存在，E0總是一個不穩(wěn)定的結(jié)點.

下面分別將雅克比矩陣J(x,y)的跡和行列式記為A和B.

平衡點(1,0)處的雅克比矩陣J1=J(E1)為

容易看出，它的特征根為λ1=-a和λ2=c-d.當(dāng)c>d時，即正平衡點存在時，E1是不穩(wěn)定的，同時，J(E1)的行列式為負，因此，E1是一個鞍點.當(dāng)c≤d時，正平衡點不存在，此時，A2-4B=(a+c-d)2≥0,且J(E1)的行列式B是非負的，因此E1是一個結(jié)點.如果c0，此時E1是一個局部漸近穩(wěn)定的結(jié)點.如果c=d，那么B=0,A<0，此時，E1是一個不穩(wěn)定的結(jié)點.

平衡點(0,1)處的雅克比矩陣J2=J(E2)為

很明顯，它的特征根為λ1=a+b和λ2=-c.因此，E2是一個鞍點.

平衡點(x*,y*)處的雅克比矩陣J*=J(E*)為

J(E*)=

經(jīng)過簡單的計算，可得A<0，B>0，這表明E*是局部漸近穩(wěn)定的.

綜上所述，我們得到下面的定理：

定理4 對于系統(tǒng)(1)，無論正平衡點是否存在，總有一個不穩(wěn)定的結(jié)點E0和一個鞍點E2.如果正平衡點存在，即c>d，那么E1是一個鞍點，E*是局部漸近穩(wěn)定的.如果正平衡點不存在，即c

3 非平凡正周期解的不存在性

下面將證明系統(tǒng)不存在非平凡正周期解.顯然，要在正平衡點存在的情況下討論.

若正平衡點外圍存在極限環(huán)Γ，則Γ必全部位于xOy坐標(biāo)平面的第一象限內(nèi)部.下面用Dulac函數(shù)法證明這種極限環(huán)不可能存在.事實上，在第一象限內(nèi)取

B(x,y)=xα-1yβ-1,

其中α，β為待定常數(shù).容易算得

xα-1yβ-1((-aα-a)x+(-βc-c)y+

由于x和y定號，欲使D不變號,只需令

從而可確定α=-1,β=-1.于是，

(9)

其中

f(s)=(b-d)s2+2ds-a-b-c,

(10)

以上的討論表明，D不變號等價于f(s)在區(qū)間[0,1]上不變號.下面將證明，只要正平衡點存在，就有f(s)≤0，s∈[0,1].

由式(10)可得，f(0)=-a-b-c<0,f(1)=-a-c+d<0.當(dāng)b-d≠0時,f(s)曲線頂點的橫坐標(biāo)為

如果b-d>0，那么x0<0，這表明在區(qū)間[0,1]上f(s)≤0；如果b-d<0，那么x0>1，這也表明在區(qū)間[0,1]上f(s)≤0.

當(dāng)b-d=0時，以上結(jié)論仍然成立.

因此，在正平衡點存在的情況下，式(9)滿足f(s)≤0，即總存在一個Dulac函數(shù)

因此，根據(jù)Bendixson-Dulac定理，可得到以下非平凡正周期解不存在的結(jié)論.

定理5 如果系統(tǒng)存在正平衡點，那么它沒有非平凡的正周期解.

通過前面的平衡點分析可知，當(dāng)E*局部漸近穩(wěn)定時，E0、E1和E2都是不穩(wěn)定的，同時，系統(tǒng)不存在非平凡正周期解.因此，由定理5可知，所有的正解都趨近于E*，于是，E*是全局漸近穩(wěn)定的.

定理6 如果c>d，那么E*是全局漸近穩(wěn)定的.

4 討論

從上面討論可以看到，相對于原模型，改進的模型(1)能夠表現(xiàn)出更豐富、更合理的動力學(xué)行為.原因在于該模型將環(huán)境資源因素考慮進來，因而能更準(zhǔn)確地體現(xiàn)人口變化規(guī)律.

當(dāng)然，本文并沒有進行數(shù)值模擬.動力學(xué)分析結(jié)果表明，該模型在參數(shù)的合理范圍內(nèi)并沒有表現(xiàn)出復(fù)雜的動力學(xué)行為.要表現(xiàn)出復(fù)雜的行為，它的參數(shù)必須突破合理的界限.因此，我們將會在今后的工作中對該模型進行模擬，確定出合理的參數(shù).