李路遙，申建華

(杭州師范大學(xué)理學(xué)院，浙江 杭州 311121)

0 引言

任意一個二維常系數(shù)線性微分系統(tǒng)u′=Au,其中A為非退化的常數(shù)矩陣，經(jīng)過適當(dāng)?shù)淖鴺?biāo)變換之后，就可以轉(zhuǎn)化成3種標(biāo)準(zhǔn)形式之一.本文討論如下標(biāo)準(zhǔn)形式的脈沖系統(tǒng)：

(1)

近年來，關(guān)于脈沖微分方程的研究已有較多深刻結(jié)果[1-7]，但有關(guān)脈沖微分方程幾何理論的研究還不多，主要困難是不連續(xù)動力系統(tǒng)的相空間理論還不完善.本文旨在研究系統(tǒng)(1)的穩(wěn)定性和相圖特征，得到的結(jié)果顯示系統(tǒng)(1)有比系統(tǒng)u′(t)=A1u(t) 更加復(fù)雜的動力學(xué)性質(zhì).

1 主要結(jié)果

經(jīng)過計算，不難得到系統(tǒng)(1)的解為

(2)

其中i=1,2,u1=x,u2=y,σ1=α,σ2=β,u01=x0,u02=y0.為了得到系統(tǒng)(1)更加復(fù)雜的相圖特征，需要引入一個引理，由于證明過程較為簡單，這里不做證明.

令Sk={(xk(t),yk(t))|xk=(1+α)k-1x0eλ1(t-t0),yk=(1+β)k-1y0eλ2(t-t0),k=1,2,3,…}, 則{Sk}是曲線族.事實上，系統(tǒng)(1)的相圖在曲線族{Sk}上.

若1+α=0或1+β=0，當(dāng)t>t1,x(t)=0，或y(t)=0，系統(tǒng)(1)的相圖在y-軸或x-軸上.因此，僅需考慮下列兩種情形：Ⅰ)α>-1,β>-1；Ⅱ)α>-1,β<-1，或α<-1,β>-1，或α<-1,β<-1.

對于情形Ⅰ，容易得到

(3)

從而可以看出，點(x(tk+),y(tk+))和曲線Sk在同一象限上，且(x(tk+),y(tk+))在曲線Sk+1上.因此，曲線Sk和Sk+1在同一象限內(nèi).由于xk(t),yk(t)都是關(guān)于t的單調(diào)函數(shù)，從而由引理1知：

(4)

此時，對系統(tǒng)(1)相圖的特征及解的穩(wěn)定性，有以下結(jié)論.

定理1對系統(tǒng)(1)有：

3) 當(dāng)Δ1>0時，系統(tǒng)(1)的相圖是間斷連續(xù)的，且λ1<0,λ2<0時，系統(tǒng)(1)的解是漸近穩(wěn)定的；當(dāng)λ1>0,λ2>0時，系統(tǒng)(1)是無窮擴散的.

(5)

當(dāng)λ1<0，可得

(6)

(7)

于是，當(dāng)m→+∞時，

因此，當(dāng)m→+∞時有x(t)→0.同理可得，當(dāng)λ2<0時，y(t)是有界的，而且y(t)→0.從而當(dāng)λ1<0,λ2<0時系統(tǒng)(1)是漸近穩(wěn)定的.

當(dāng)λ1>0，從式(5)可以得到

(8)

于是對?t∈(tm,tm+1]，當(dāng)m→+∞時，

即x(t)→+∞(t→+∞).同理可得，當(dāng)λ2>0時，有y(t)→+∞(t→+∞).從而系統(tǒng)(1)是無窮擴散的.

下證穩(wěn)定性.由于

(9)

類似于1)的證明，可以得到當(dāng)λ1>0,λ2>0時，系統(tǒng)(1)是漸近穩(wěn)定的.當(dāng)λ1<0,λ2<0時，系統(tǒng)(1)是無窮擴散的.

3) Δ1>0，由于tk+Δ1>tk，系統(tǒng)(1)在任意兩個相鄰區(qū)間(tk,tk+1]和(tk-1,tk]上沒有重疊部分，因此相圖間斷連續(xù).下證穩(wěn)定性.由于

所以，當(dāng)λ1<0,λ2<0,ln|1+α|<0,ln|1+β|<0時，有

當(dāng)λ1>0,λ2>0,ln|1+α|>0,ln|1+β|>0時，有

所以，當(dāng)λ1<0,λ2<0時，系統(tǒng)(1)是漸近穩(wěn)定的.當(dāng)λ1>0,λ2>0時，系統(tǒng)(1)是無窮擴散的.

(10)

可以看出點(x(tk+),y(tk+))和曲線Sk不在同一象限，但與曲線Sk-1在同一象限，這也說明曲線Sk+1與曲線Sk不在同一象限，但與曲線Sk-1在同一象限.由于xk(t),yk(t)都是關(guān)于t的單調(diào)函數(shù)，從而由引理1知：

(11)

此時，對系統(tǒng)(1)相圖的特征及解的穩(wěn)定性，有以下結(jié)論.

定理2對系統(tǒng)(1)有：

系統(tǒng)的熱容Ct是熱水器的固有屬性，文獻[7]中通過實驗數(shù)據(jù)的樣本點進行最小二乘法的擬合并通過線性回歸進行校核分析，得到的最終待定系數(shù)Ct的取值為21.7～25.3 kJ/K。在此，將熱水溫升與加熱時間的數(shù)學(xué)模型導(dǎo)入MATLAB軟件進行模擬計算，系統(tǒng)熱容Ct以2為步長進行劃分，其余參數(shù)取均值且保持不變，具體參數(shù)見表1，熱水溫升與加熱時間的關(guān)系曲線如圖3所示。

于是有

(12)

當(dāng)λ1<0時容易得到

(13)

因此

從而

(14)

因此x(t)→0(t→+∞).同理可得，當(dāng)λ2<0時，y(t)→0(t→+∞).從而系統(tǒng)(1)是漸近穩(wěn)定的.當(dāng)λ1>0，從式(13)可得

(15)

(16)

即x(t)→+∞(t→+∞).同理可得，當(dāng)λ2>0時，有y(t)→+∞(t→+∞).從而系統(tǒng)(1)是無窮擴散的.

于是有

(17)

因此，類似1)的證明，可以得到以下結(jié)論：當(dāng)λ1>0,λ2>0時，系統(tǒng)(1)是漸近穩(wěn)定的；當(dāng)λ1<0,λ2<0時，系統(tǒng)(1)是無窮擴散的.

(18)

于是有

因此，可以得到以下結(jié)論：當(dāng)λ1<0,λ2<0時，

當(dāng)λ1>0,λ2>0時，

所以當(dāng)λ1<0,λ2<0時，系統(tǒng)(1)是漸近穩(wěn)定的；當(dāng)λ1>0,λ2>0時，系統(tǒng)(1)是無窮擴散的.

2 數(shù)值模擬

下面僅給出當(dāng)λ1<0,λ2<0時的數(shù)值模擬.

對于情形Ⅰ，當(dāng)α>-1,β>-1時，取定參數(shù)λ1=-0.6,λ2=-0.3,x0=100,y0=120,tk=5k.

(a)函數(shù)圖

(a)函數(shù)圖

由圖1、2可以看出系統(tǒng)(1)是漸近穩(wěn)定的，且此時相圖在整個區(qū)間上看起來是連續(xù)的，但實際是間斷連續(xù)的.

3) 當(dāng)Δ1>0時，取α=-0.36,β=-0.2，得到的函數(shù)圖和相圖見圖3.

由圖3可以看出系統(tǒng)(1)是漸近穩(wěn)定的，且此時相圖在整個區(qū)間上明顯是間斷連續(xù)的.

考慮情形Ⅱ，即α>-1,β<-1，或α<-1,β>-1，或α<-1,β<-1.

從圖4可以看出系統(tǒng)(1)是漸近穩(wěn)定的，從圖5可以看出相圖在整個區(qū)間上看起來是連續(xù)的，但實際上是間斷連續(xù)的.

從圖6可以看出系統(tǒng)(1)是無窮擴散的，從圖7可以看出相圖在兩個象限上看起來是連續(xù)的，但其實是間斷連續(xù)的.

但在不同象限的相圖見圖9.

由圖8、圖9可以發(fā)現(xiàn)，系統(tǒng)(1)是漸近穩(wěn)定的，而且相圖明顯是間斷連續(xù)的.