差分方程方法

1、第三章差分方程方法3.1差分方程的平衡點及其穩(wěn)定性設有未知序列 xn ,稱F(n；Xn,Xn 1, Xn k)0(3.1)為k階差分方程。若有xnx(n),滿足F(n; x(n), x(n 1), , x(n k) 0則稱Xnx(n)是差分方程(3.1)的解，包含k個任意常數(shù)的解稱為(3.1)的通解，X0,X1,xk1為已知時，稱其為(3.1)的初始條件，通解中的任意常數(shù)都由初始條件確定后的解稱為(3.1)的特解。形如Xn k a(n)Xnk1ak(n)Xnf (n)(3.2)的差分方程，稱為 k階線性差分方程。ai(n)為已知系數(shù)，且ak(n) 0。若差分方程(3.2)中的f(n) 0,則稱

2、差分方程(3.2)為k階齊次線性差分方程，否則 稱為k階非齊次線性差分方程。若有常數(shù) 是差分方程(3.1)的解，即F(n; , , , ) 0,則稱是差分方程(3.1)的平衡點，又對差分方程(3.1)的任意由初始條件確定的解xn x(n),都有Xn(n ),則稱這個平衡點是穩(wěn)定的。若X0,X1,Xk1已知，則形如 Xn k g(n；Xn,Xn1, , Xn k 1)的差分方程的解可以在計算機上實現(xiàn)。下面給出理論上需要的一些特殊差分方程的解。一階常系數(shù)線性差分方程Xn 1Xnb ,(3.3)(其中 2為常數(shù)，且 1,0)的通解為XnC( )n b (a 1)(3.4)易知b/(1)是方程(3.3

3、)的平衡點，由(3.4)式知，當且僅當1時，b/(1)是穩(wěn)定的平衡點。二階常系數(shù)線性差分方程Xn 2Xn 1 bXnr ,(3.5)其中a,b,r為常數(shù)，當r 0時，它有一特解x 0；當r 0,且a b 1 0時，它有一特解X r/(a b 1)。不管是哪種情形，X是方程(3.5)的平衡點。設方程(3.5)的特征方程2 a b 0的兩個根分別為1,2，則當1,2是兩個不同實根時，方程(3.5)的通解為XnX* Ci( i)n C2( 2尸；當1,2是兩個相同實根時，萬程 (3.5)的通解為Xnx(C1 C2n) n ；當1,2 (cos isin )是一對共軻復本|時，方程 (3.5)的通解為

4、* nXn x(C1 cosnC2sin n )一* 一 >，.、 一一 一一* . *易知，當且僅當特征萬程的任一特征根i 1時，平衡點X是穩(wěn)定的。二階方程的上述結果可以推廣到k階線性方程，即k階線性方程平衡點穩(wěn)定的條件是特征方程白根i(i 1,2, ,k)均滿足：i 1 ,即均在復平面上的單位圓內。下面討論一階非線性差分方程Xn 1f(Xn)(3.6)，一，一、一，> * *的平衡點的穩(wěn)定性。 其平衡點X由代數(shù)方程X f(x)解出，為分析平衡點X的穩(wěn)定性，將*萬程(3.6)的右端f(Xn)在X點作泰勒展開，只取一次項，(3.6)近似為_ * * * _ * _Xn 1f (X

5、)(Xn X ) f (X )(3.7)*(3.7)是(3.6)的近似線性萬程，X也是(3.7)的平衡點。關于線性方程平衡點穩(wěn)定的條件上面已給出，而當 f (X )1時，方程(3.6)與(3.7)平衡點的穩(wěn)定性相同。于是得到，當_ '* . I . *. ._ '* .f (X )1時，對于非線性方程(3.6) , X是穩(wěn)定的；當f (X )1時，對于方程6),* 一 .一 、一X是不穩(wěn)止的。3.2市場經(jīng)濟中的蛛網(wǎng)模型在自由貿易的集市上你注意過這樣的現(xiàn)象嗎？ 一個時期由于豬肉的上市量遠大于需求，銷售不暢導致價格下降，農(nóng)民覺得養(yǎng)豬賠錢， 于是轉而經(jīng)營其它農(nóng)副業(yè)。過一段時間后豬肉上

6、市量大減，供不應求導致價格上漲，原來的飼養(yǎng)戶看到有利可圖，又重操舊業(yè)。這樣下一個時期會重視供大于求、價格下降的局面，在沒有外界干擾的情況下，這種現(xiàn)象將如此循環(huán) 下去。在完全自由競爭的市場經(jīng)濟中上述現(xiàn)象通常是不可避免的，因為商品的價格是由消費者的需求關系決定的，商品數(shù)量越多價格越低。而下一時段商品的數(shù)量由生產(chǎn)者的供應關系決 定，商品價格越低生產(chǎn)的數(shù)量就越小。這樣的需求和供應關系決定了市場經(jīng)濟中商品的價格 和數(shù)量必然是振蕩的。這種振蕩越小越好，如果振蕩太大就會影響人們群眾的正常生活?，F(xiàn)在我們用差分方程建模，討論市場經(jīng)濟趨于穩(wěn)定的條件，再用圖形方法建立“蛛網(wǎng)模型”對上述現(xiàn)象進行分析，對結果進行解釋，

7、然后作適當推廣。記第n時段商品數(shù)量為Xn,價格為yn , n 1,2,。這里我們把時間離散化為時段，1個時段相當于商品的1個生產(chǎn)周期，如蔬菜、水果可以是一年，肉類則是1個飼養(yǎng)周期。在n時段商品的價格yn取決于數(shù)量Xn ,設yf (Xn) o它反映消費者對這種商品的需求關系，稱需求函數(shù)。因為商品的數(shù)量越多價格越低，所以在圖3-1中用一條下降曲線 f表示它，f為需求曲線。在n 1時段商品的數(shù)量 xn 1由上一時段價格yn決定，用xn1g(yn)表示，它反映生產(chǎn)者的供應關系，稱為供應函數(shù)。因為價格越高生產(chǎn)量才越大，所以在圖圖4-1需求函數(shù)和供應函數(shù)3-1中供應曲線g是一條上升曲線。設圖3-1中兩條曲

8、線f和g相交于Po(x0, y0)點，在P0點附近取函數(shù)f和g的線性近似，即需求函數(shù)(f): yn yo(Xn Xo) ,0(3.8)供應函數(shù)(g)： Xn 1 Xo (yn yo),0(3.9)由式(3.8) , (3.9)消去yn得一階線性差分方程Xn 1Xn (1)X0 , n 1,2, .(3.10)因此Xo是其平衡點，即 Po是平衡點，對n遞推不難得到Xn 1 ()nX1 1 ()nXo，n 1,2,.由此可得，平衡點 P0穩(wěn)定的條件是：1;不穩(wěn)定的條件是：1。下面用圖形解釋此模型。若對某一個k有xk x0,遞推可得，當n k時xn x0 ,從而yn y0 ,即商品的數(shù)量和價格將永遠

9、保持在 P0(x0, y0)點。但是實際生活中的種種干擾使得xn, yn不可能停止在P。點上。不妨設xi偏離x0 ,我們分析隨著n的增加xn , yn的變化。數(shù)量xi給定后，價格yi由曲線f上的Pi點決定，下一時段的數(shù)量x2由曲線g上的P2點決定，這樣得到圖案稀刖的甲(, y1) , P2(x21y1),圖田口 2海第白屯0>3, y2),，在圖3-2上這些點將按箭頭所示方向趨向P0(x0,y0),表明P。是穩(wěn)定平衡點，意味著市場經(jīng)濟(商品的數(shù)量和價格) 將趨向穩(wěn)定。但是如果需求函數(shù)和供應函數(shù)由圖3-3的曲線所示，則類似的分析發(fā)現(xiàn)，市場經(jīng)濟將按照P1 ,P2,P3,P4,，的規(guī)律變化而遠

10、離P0,即P0是不穩(wěn)定平衡點，市場經(jīng)濟趨向不穩(wěn)定。由此可見，需求曲線越平，供應曲線越陡，越有利 于經(jīng)濟穩(wěn)定。圖3-2和圖3-3中折線P1F2F3F4 形似蛛網(wǎng)，于是這種用需求曲線和供應曲線分析市場經(jīng)濟穩(wěn)定性的圖示法在經(jīng)濟學中被稱為蛛網(wǎng)模型。實際上，需求曲線f和供應曲線g的具體形式通常是根據(jù)各個時段商品的數(shù)量和價格的一系列統(tǒng)計資料得到的。一般地說，f取決于消費者對這種商品的需要程度和他們的消費水平，g則與生產(chǎn)者的生產(chǎn)能力，經(jīng)營水平等因素有關。 下面再來解釋此模型的實際意義。首先考察參數(shù)，的含義，需求函數(shù) f的斜率(取絕對值)表示商品供應量減少1個單位時價格的上漲幅度；供應函數(shù)g的斜率 表示價格上

11、漲1個單位時(下一時期)商品供應增加量。所以的數(shù)值反映消費者對商品需求的敏感程度。如果這種商品是生活必需品，消費者處于持幣待購狀態(tài)，商品數(shù)量稍缺，人們立即蜂擁搶購，那么會比較大；較小。 的反之，若這種商品非必需品，消費者購物心理穩(wěn)定，或者消費水平低下，則數(shù)值反映生產(chǎn)經(jīng)營者對商品價格的敏感程度，如果他們目光短淺，熱衷于追逐一時的高利潤,價格稍有上漲立即大量增加生產(chǎn)，那么會比較大；反之，若他們素質較高，有長遠的計劃，則 較小。根據(jù) ，的意義很容易對市場經(jīng)濟穩(wěn)定與否的條件作出解釋。當供應函數(shù)g ,即固定時， 越小，需求曲線越平，表明消費者對商品需求的敏感程度越小，越利于經(jīng)濟穩(wěn)定。當需求函數(shù) f ,即

12、 固定時， 越小，供應曲線越陡，表明生產(chǎn)者對價格的敏感程度 越小，越利于經(jīng)濟穩(wěn)定。反之，當 ， 較大，表明消費者對商品的需求和生產(chǎn)者對商品的價格都很敏感，則會導致經(jīng)濟不穩(wěn)定。當市場經(jīng)濟趨向不穩(wěn)定時，政府有兩種干預辦法：一種辦法是控制價格，無論商品數(shù)量多少，命令價格不得改變，于是 0,不管曲線g如何，總是穩(wěn)定的；另一種辦法是控制市場上的商品數(shù)量，當上市量小于需求時，政府從外地收購或調撥，投入市場，當上市量多于需求時，政府收購過乘部分，于是0,不管曲線f如何，也總是穩(wěn)定的。最后我們將模型稍加推廣。如果生產(chǎn)者的管理水平更高一些，他們在決定商品生產(chǎn)數(shù)量時，不是僅根據(jù)前一時期的價格，而是根據(jù)前兩個時期的

13、價格，為簡單起見不妨設根據(jù)二者的平均值(yn yn 1)/2 ,于是供應函數(shù)為xn 1 g (yn yn i/2)。在Po點附近取線性近似時，(3.9)表示為：供應函數(shù)(gXn 1 Xo(yn ynl/2 yo) ,0。(3.11)又設需求函數(shù)仍由(3.8)式表示，則由(3.8) , (3.11)式得到2Xn 2Xn 1Xn (1 )X0, n 1,2, .(3.12)(3.12)式是二階線性差分方程。Po點穩(wěn)定的條件可由其特征方程2 20的根1,2、.()2 84確定。8時，顯然有2,從而22, P0是不穩(wěn)定的。當 8時，可以算出 1,2 J2-,由1,21得到P0點穩(wěn)定的條件為2。這與原有

14、模型中的 P0點穩(wěn)定的條件1相比，保持經(jīng)濟穩(wěn)定的參數(shù)， 的范圍放 大了(， 的含義未變)?？梢韵氲?，這是生產(chǎn)經(jīng)營者的生產(chǎn)管理水平提高，對市場經(jīng)濟 穩(wěn)定起著有利影響的必然結果。3.3差分形式的阻滯增長模型阻滯增長模型dx dtrx(1 )K(3.13)的差分形式是將其化為Yn 1Ynryn(1 Yn . K )yn 1(r 1)yn(1ryn. r 1 K)令b r 1, Xn ryn/ (r 1)K，將上式簡化為Xn 1bXn 1 Xn(3.14)方程(3.14)的解可根據(jù)給定的初值X0,利用計算機算出:根據(jù)r , K的含義，有b 1,0 Xn1 ,再由(3.14)應有b 4。下面主要討論當1

15、 b 4時，序列Xn的收斂性，它比簡單地討論方程(3.13)或(3.14)的穩(wěn)定性復雜，且會看到一個十分有趣的現(xiàn)象。記f(X) bX(1 X),解代數(shù)方程X f(X)得方程(3.14)的平衡點X*0,1 1b.由于f (x) b(1 2x),注意到b 1及平衡點穩(wěn)定的條件f (x ) 1得，0不是穩(wěn)定的平衡點，x 1 1/b的狀況可以通過方程(3.14)的圖解法清楚地表示出來。在xOy平面. 一 . 一'，*、 一 上回出y x和y f(x)的圖形(圖3-4)。由于f (x ) 2 b ,因此當1 b 2時，xn(圖3-4 (1)；當2 b 3時，Xn基本上是衰減震蕩地收基本上是單調遞

16、增地收斂于*八人. * 一斂于 X (圖 3-4 (2)。圖 3-4 (1)圖 3-4 (2)當b 3時，雖然萬程(3.14)仍可形式地求解，但 X不穩(wěn)定，其圖解法如圖3-5所示。事情到此并未完結，由(3.14)式迭代一次可得Xn 2bXn i (1 Xn i)bbXn(1 Xn ) 1 bXn (1Xn)即Xn 2b2Xn(1 Xn)1 b4 (1 4)(3.15)方程(3.15)雖是二階非線性差分方程，但缺少 Xn 1,相當于一階差分方程。解代數(shù)方程x b2X(1x) 1 bXn(1 Xn)得到方程(3.15)的4個平衡點，除了方程(3.14)的2個平衡點0, 1 1/b外，還有兩個*X1

17、,2b 1 b2 2b 32b圖3-6圖3-5可以驗證(根據(jù)一階差分方程的判別方法)，在條件b 3下，平衡點0, 1 1/b不是-4,、,，一* _._._.穩(wěn)定的，而X1,2是穩(wěn)定的條件為b 1 763.449。這就是說，當3b3.449時，雖然序列Xn不收斂，但它的兩個子序列X2n和X2n 1都是收斂的(圖 3-6)。它的生物學意義是，當固有增長率 2 r 2.449時，從一個繁殖周期(即一代)的角度看，其數(shù)量增長 是不穩(wěn)定的，但從兩個繁殖周期(即兩代)的角度看，其數(shù)量增長是穩(wěn)定的，這就是所謂的 2倍周期收斂。讀者不難想到，當b 3.449時，x；2不再是穩(wěn)定白1即方程(3.15)不存在穩(wěn)定的平衡點，從而對于方程(3.14)來說2倍周期也不收斂了，但是將方程(3.15)迭代一次或者將方程 (3.14)迭代2次，得Xn 4f ()(3.16)方