標籤: 暫無標籤

差分方程又稱遞推關係式,是含有未知函數及其差分,但不含有導數的方程。滿足該方程的函數稱為差分方程的解。差分方程是微分方程的離散化。

1簡介

在數學上,遞推關係(recurrence relation),也就是差分方程(difference equation),是一種遞推地定義一個序列的方程式:序列的每一項目是定義為前一項的函數。某些簡單定義的遞推關係式可能會表現出非常複雜的(混沌的)性質,他們屬於數學中的非線性分析領域。
所謂解一個遞推關係式,也就是求其解析解,即關於n的非遞歸函數。

2意義性質

性質
性質1 Δk(xn+yn)=Δkxn+Δkyn
性質2 Δk(cxn)=cΔkxn
性質3 Δkxn=∑(-1)jCjkXn+k-j
性質4 數列的通項為n的無限次可導函數,對任意k>=1,存在η,有 Δkxn=f(k)(η)
差分方程
定義8。1 方程關於數列的k階差分方程:
xn-a1xn-1-a2xn-2-……aBxn-k=b(n=k,k+1,……)
其中a1,a2,------ak 為常數, ak≠0. 若b=0,則該 方程是齊次方程
關於λ 的代數方程
λk-a1λk-1-------ak-1λ-ak=0
為對應的特徵方程,根為特徵值。

3例題

1. 實驗內容與練習
2.1 差分
例1 Xn={n3},求各階差分數列:
xn △xn △2xn △3xn △4xn
1 7 12 6 0
8 19 18 6 0
27 37 24 6 0
64 61 30 6
125 91 36
216 127
343
可見,{n3},三階差分數列為常數數列,四階為0。
練習1 對{1},{n},{n2},{n4},{n5},分別求各階差分數列。
練習2 {C0n-1}{C1n-1}{C2n-1},{C4n-1},分別求各階差分數列.
{Xn}的通項為n的三次函數,
Xn=a3n3+a2n2+a1n+a0
證明它為常數數列。
證明 由Xn=a3n3+a2n2+a1n+a0可直接計算。
定理8。1 若數列的通項是關於n 的k次多項式,則 k 階差分數列為非零數列,k+1階差分數列為0。
練習3 證明定理8。1。
定理8。2 若{Xn}的 k 階插分為非零常數列,則{Xn}是 n的 k次多項式,
練習4 根據差分的性質證明定理8。2
例2。求∑i3
例4
解 設Sn=∑i3 表
Sn △Sn △2Sn △3Sn △4Sn △5Sn
1 8 19 18 6 0
9 27 37 24 6 0
36 64 61 30 6 0
100 125 91 36 6 0
225 216 127 42
441 343 169
784 512
1296
設Sn=a4n4+a3n3+a2n2+a1n+a0,s1=1,s2=9,s3=36,s4=100,s5=225,得
a0=0,a1=0,a2=1/4,a3=1/2,a4=1/4.
所以, Sn=(1/4)n4+(1/2)n3+(1/4)n2.
練習 {Xn}的通項Xn為n的k次多項式,證明∑xi為n的 k+1次多項式;求 ∑i4.
由練習 2 {Crn-1}可得。
2.2差分方程
對於一個差分方程,如果能找出這樣的數列通項,將它帶入差分方程后,該方程成為恆等式,這個通項叫做差分方程的解。
例3 對差分方程 xn-5xn-1+6xn-2=0,可直接驗證xn=c13n+c22n是該方程的解。
例3中的解中含有任意常數,且任意常數的個數與差分方程的階數相同。這樣的解叫做差分方程的通解。
若k階差分方程給定了數列前k項的取值,則可以確定通解的任意常數,得到差分
的特解。
例4對差分方程xn-5xn-1+6xn-2=0,若已知x1=1,x2=5,則可以得到該差分方程的特解為xn=3n-2n.
我們首先研究齊次線性差分方程的求解。
xn=rxn-1
對一階差分方程
x1=a
顯然有xn=arn-1。因此,若數列滿足一階差分方程,則該數列為一個等比數列。
例5 求Fibonacci數列{Fn}的通項,其中F1=1,F2=1,Fn=Fn-1+Fn-2.
Fibonacci數列的前幾項為:1,1,2,3,5,8,13,21,34,55,89,…。該數列有著非常廣泛的應用。
Fibonacci數列所滿足的差分方程為 Fn-Fn-1-Fn-2=0,
其特徵方程為 λ2-λ-1=0
其根為λ1=,λ2= .利用λ1λ2可將差分方程寫為
Fn-(λ1+λ2)Fn-1+λ1λ2Fn-2=0,
即Fn-λ1Fn-1=λ2(Fn-1-λ1Fn-2)
數列{Fn-λ1Fn-1}滿足一個一階差分方程.顯然 ( )
同理可得 ( )
由以上兩式可解出 的通項。
練習9 證明若數列{ }滿足二階差分方程 ,其特徵方程由兩個不相等的根 ,則 為該差分方程的兩個特解。從而其通解為。
由練習9,若二階差分方程的特徵方程有兩個不相等的根,可寫出其通解的一般性式。再由 的值可解出其中的係數,從而寫出差分方程的特解。
練習10 具體求出 Fibonacci數列的通項,並證明。那麼,若二階線性齊次差分方程有兩個相等的根,其解有如何來求呢?
設二階線性齊次差分方程的特徵方程有兩個相等的根 ,則差分方程可寫為。差分方程的兩邊同時除以 ,有。設,則 (n>=3)。由於該式在 n>=3式均成立,我們將它改寫為 (n>=1)。(8.2)
方程(8.2)的左邊是 的二階差分,從而有,於是 是n的一次函數,設為 則有。上是即為差分方程的通解。
練習11 證明:若數列{ } 所滿足的三階差分方程的特徵方程由三個相等的根 ,則差分方程的通解為。
一般的,設 ···,為差分方程的特徵方程所有不同的解,其重數分別為 ···, ,則差分方程對應於其中的根 (i=1,2,···,l)的特解 ···。
對於一般的k階齊次線性差分方程,我們可以通過其特徵方程得到上述形式的k個特解,進而得到差分方程的通解。
練習12 若數列{ } 滿足差分方程
且 求{ }的通項。
例6 若實係數差分方程的根為虛數,則其解也是用虛數表示的,這給討論問題帶來不便。差分方程
xn-2xn-1+4xn-2=0
的特徵值為 i.若x1=1,x2=3,由下面的程序易求出其特解為:
xn=( )(1+ i)n+(- )(1- i)n
Clear[x1,x2,c1,c2,l1,l2,solution];
x1=1;x2=3;
solution=Solve[1^2-2l+4==0,1];
l1=l/.solution[[1,1]];
l2=l/.solution[[2,1]];
c=Solve[ {c1*l1+c2*l2==x1,c1*l1^2+c2*l2^2==x2},{c1,c2}];
c1=Simplify[ Re[c1]]+Simplify]*I;
c2=Simplify[Re[c2]]+Simplify]*I;
Print[「xn=(「,c1,」)(「,l1,」)^n+(「,c2,」)(「,l2,」)^n」]
解的形式相當複雜,是否可以將它們用實數表示呢?
設 =rei,則 =re,我們可將(8.4)中的表達式改寫為
xn=re (2e )n+re (2e )\n
=r
=2r Cos( )
=(2rCos )
=
可以看出,通項可以寫成 的形式.那麼, 與 是不是差分方程的特解呢?
練習13 驗證 與 是差分方程(8.3)的特解.
對於差分方程(8.3),我們找出了它的兩個實型的特解,從而可以將通解表示成實數的形式.這一方法對於一般的方程也是成立的.
練習14 設 的兩個特徵值為 .證明該差分方程的通解可表示為 .
練習 15 用實數表示差分方程 的特解.
上次我們討論了其次線性差分方程的求解方法.那麼,非齊次線性差分方程是否可以化為齊次線性差分方程呢?
練習16 若已知非齊次線性差分方程
··· (8.5)
的一個特解為 求證:若令 則 滿足齊次差分方程
···
由練習16,若已知非齊次線性差分方程(8.5)的一個特解,就可以將它化為齊次線性差分方程.
顯然方程(8.5)的最簡單的形式為 (其中p為常數),代入(8.5)得
···
若 ··· 則有
稱p = 為非齊次線性差分方程(8.5)的平衡值。在(8.5)中, 令 則有
由 ,得
.
從而可將原來的非齊次線性差分方程化為齊次線性差分方程.
如果方程(8.5)的平衡值不存在,可以將方程(8.5)中所有的n換為n+1,得到
(8.6)
方程(8.6)和(8.5)相減得
.
於是可將原來的非齊次線性差分方程化為高一階的齊次線性差分方程.
練習17 分別求差分方程 及 的通解.
2.3代數方程求根
由 Fibonacci數列的性質,我們可以用 來逼近 ,用這一性質可以來計算 的近似值。一般地,對a>0,可以用構造差分方程的方法來求 的近似值.
對給定的正數a,設λ1= ,λ2= ,則λ1 ,λ2是方程λ2-2λ+(1+a)=0的根.該方程是差分方程 的特徵方程。於是,選定,利用差分方程 可以構造一個數列{ }.
練習 18 證明:若a>1,對任意的 >0,>0,若 ≠ ,則按上述法構造的數列{ }滿足
.這樣,我們得到了計算 的一個方法:
1. 給定 (作為誤差控制),任取初始值 ,令n=1;
2. 若
則終止計算,輸出結果;否則 ,令n :=n+1,轉第3步;
3. 令,轉第2步.
練習 19 對a=1.5,10,12345,用上述方法求 .
上述方法的收斂速度不夠快,我們可以加以改進
設整數u滿足,令,則 , 是方程 的兩個根.
練習 20 根據上面的差分方程的構件數列{ x },使得
.
練習 21 對練習19中的a,用上面的方法來計算 ,並比較兩種方法的收斂速度.
代數方程
(8.7)
是差分方程(8.1)的特徵方程,是否可以用此差分方程來求解方程(8.7)呢?
設方程(8.7)有k個互不相同的根滿足
, (8.8)
則對應的差分方程的通解形式為
.
練習 22 設方程(8.7)的根滿足條件(8.8),任取初始值 用差分方程(8.1)(取b=0)構造數列{ }.若通解中 的係數 ≠0,證明:
.
利用練習22得到的結論,我們可以求多項式方程的絕對值最大的根.
練習 23 求方程 的絕對值最大的根.
事實上,若方程(8.7)的互不相同的根滿足
≥ ≥…≥
(其重數分別為 ),則練習22中的結論仍然成立.
2.4 國民收入4 國民收入的穩定問題
一個國家的國民收入可用於消費,再生產的投資等。一般地說,消費與再生產投資都不應該沒有限制。合理的控制各部分投資,能夠使國民經濟處於一種良性循環之中。如何配各部分投資的比例,才能使國民經濟處於穩定狀態呢?這就是本節要討論的問題。
我們首先給出一些假設條件:
1. 國民收入用於消費、再生產投資和公共設施建設三部分。
2. 記 分別為第k個周期的國民收入水平和消費水平。的值與前一個周期的國民收入成正比例。即 =A,(8.9)其中A為常數(0 3. 用 表示第k個周期內用於再生產的投資水平,它取決於消費水平的變化,即 . (8.10)
4. G表示政府用於公共設施的開支,設G為常數.由假設1有 . (8.11)上式是一個差分方程,當給定 的值后,可直接計算出國民收入水平 (k=2,3,…)來觀察其是否穩定。
例7 若 ,計算可得表8.3中數據。
表8.3 Y 的值的變化
k 2 3 4 5 6 7 8 9 10 11
11.0 24.5 35.8 39.1 32.9 20.3 7.48 0.95 3.93 15.0
k 12 13 14 15 16 17 18 19 20 21
28.5 37.8 38.2 29.5 16.0 4.58 0.82 6.65 19.2 32.1
我們可以畫出 的散點圖來觀察其變化。其計算及畫圖的程序如下:
y0=2;y1=2;a=0.5;b=2;g=10;
y={y0,y1};
For[k=1,k<=20,k++,
Y2=a(1+b)*y1-b*a*y0+g;
Y=Append[y,y2];
Y0=y1,y1=y2]
YListPlot[y,PlotJoined True,
PlotStyle Thickness[0.012]]
圖8.1 國民收入 的變化
由圖8.1利用發現,又例7的數據得出的 的呈現出周期變化的跡象。
練習 24設 ,對於表8.4中的參數A,B,分別計算 (k=2,3,…)並畫圖觀察 的變化。
表8.4 參數A,B的取值
A 1/2 1/2 1/2 8/9 9/10 3/4 4/5
B 1 2 3 1/2 1/2 3 3
可以看出,隨著參數的值不同,國民收入水平 (k=2,3,…)的穩定性呈現出不同的狀態。
那麼,參數滿足什麼條件時,國民收入水平才處於穩定發展之中呢?
差分方程(8.11)是一個常係數非齊次線性差分方程。由A<1容易求出其平衡值為
令 可得
.
其特徵值為
若 則
其中 為 的幅角。
從而可的差分方程的解為
其中 為常數。
若 易見{ }為一周期函數在 ---的取值,從而{ }呈周期變化的狀態。正如在例7中所見到的。
練習25 若 在 及 的情形下,討論{ }的變化趨勢。國民收入會穩定發展嗎?
練習26 若 ,國民收入在什麼條件下會穩定發展?
本實驗涉及的Mathematica軟體語句說明
1. solution=Solve[1^2-2l+4==0,1];
l1=1/.solution[[1,1]];
l2=l/.solution[[2,1]];
將方程l^2-2l+4==0的兩根分別賦值給l1及l2.
2. c=Solve[{c1*l1+c2*l2==x1,c1*l1^2+c2*l2^2==x2},{c1,c2}];
{c1,c2}={c1,c2}/.c[[1]];
將方程組{c1*l1+c2*l2==x1,c1*l1^2+c2*l2^2==x2}的解賦值給c1及c2.
3. c1=Simplify[Re[c1]]+Simplify]*I
將複數c1化簡.

4線性差分

定理
定理1(齊次線性差分方程解的疊加原理)
若y1(t),y2(t),…,ym(t)是齊次線性差分方程yt+n+a1yt+n-1 +a2yt+n-2+…+an-1yt+1+anyt=0的m個特解(m≥2),則其線性組合y(t)=A1y1(t)+A2y2(t)+…+Amym(t)也是方程 的解,其中A1,A2,…,Am為任意常數。
定理2n階齊次線性差分方程yt+n+a1yt+n-1 +a2yt+n-2 +…+an-1yt+1+anyt=0一定存在n個線性無關的特解。
定理3(齊次線性差分方程通解結構定理)
如果y1(t),y2(t),…,yn(t)是齊次線性差分方程yt+n+a1yt+n-1 +a2yt+n-2 +…+an-1yt+1+anyt=0的n個線性無關的特解,則方程 的通解為:
yA(t)=A1y1(t)+A2y2(t)+…+Anyn(t),
其中A1,A2,…,An為n個任意(獨立)常數。
定理4(非齊次線性差分方程通解結構定理)
如果 (t)是非齊次線性方程yt+n+a1(t)yt+n-1+a2(t)yt+n-2 +…+an-1(t)yt+1+an(t)yt=f(t)的一個特解,yA(t)是其對應的齊次線性方程yt+n+a1yt+n-1 +a2yt+n-2 +…+an-1yt+1+anyt=0的通解,那麼,非齊次線性差分方程的通解為:
y(t)=yA(t)+ (t)
y(t)=A1y1(t)+A2y2(t)+…+Anyn(t)+ (t),
這裡A1,A2,…,An為n個任意(獨立)常數。

5通解特解

齊次差分方程的通解
將方程yt+1+ayt=0改寫為:yt+1=-ayt,t=0,1,2,…。假定在初始時刻(即t=0)時,函數yt取任意值A,那麼由上式逐次迭代,算得
y1=-ay0=-aA,y2=-ay1=(-a)2A,………………
差分方程

  差分方程

方程的通解為yt =A(-a)t ,t=0,1,2,…
如果給定初始條件t=0時yt=y0,則A=y0,此時特解為:yt =y0(-a)t
差分方程

  差分方程

非齊次方程的通解與特解
迭代法求通
將方程改寫為 yt+1=(-a)yt+f(t), t=0,1,2,…。
逐步迭代,則有
y1=(-a)y0+f(0),y2=(-a)2y0+(-a)f(0)+f(1),y3=(-a)3y0+(-a)2f(0)+(-a)f(1)+f(2),………………
由數學歸納法,可得
差分方程
其中
差分方程
為方程的特解。yA(t)=(-a)ty0為對應的齊次方程的通解。

6應用

存款模型
設St為t期存款總額,i為存款利率,則St與i有如下關係式:
St+1=St+iSt=(1+i)Si, t=0,1,2,…
其中S0為初始存款總額。
動態供需均衡模型(蛛網定理)
差分方程

  差分方程

設Dt表示t期的需求量,St表示t期的供給量,Pt表示商品t期價格,則傳統的動態供需均衡模型為:
差分方程

  差分方程

差分方程
其中a,b,a1 ,b1均為已知常數。
(1)式表示t期(現期)需求依賴於同期價格;
(2)式表示t期(現期)供給依賴於(t-1)期(前期)價格。
差分方程

  差分方程

(3)式為供需均衡條件。
若在供需平衡的條件下,而且價格保持不變,即 Pt=Pt-1=Pe,靜態均衡價格
差分方程

  差分方程

差分方程
需求曲線與供給曲線的交點(Pe ,Qe)即為該種商品的靜態均衡點。
差分方程

  差分方程

動態供需均衡模型的等價差分方程
差分方程

  差分方程

差分方
方程的一個特解
差分方
方程的通解為
差分方程
若初始價格P0已知時,將其代入通解,可求得任意常數A=P0-Pe ,此時,通解改寫
差分方程

  差分方程

差分方程
如果初始價格P0=Pe ,那麼Pt=Pe ,這表明沒有外部干擾發生,價格將固定在常數值Pe上,即靜態均衡。如果初始價格P0≠Pe ,那麼價格Pt將隨t的變化而變化。
差分方程
<1時,
差分方程
動態價格Pt隨著t的無限增大逐漸地振蕩趨近於靜態均衡價格Pe。
差分方程
凱恩斯(Keynes.J.M)乘數動力學模
差分方程

  差分方程

設Yt表示t期國民收入,Ct為t期消費,It為t期投資,DI0為自發(固定)投資,I為周期固定投資增量。凱恩斯國民經濟收支動態均衡模型為:
差分方程
(1)式為均衡條件,即國民收入等於同期消費與同期投資之和;(2)式為消費函數,即現期消費水平依賴於前期國民收入(消費滯後於收入一個周期),a(≥0)為基本消費水平,b為邊際消費傾向(0<b<1);(3)式為投資函數,這裡僅考慮為固定投資。
在(1)(2)(3)式中消去Ct和It,得到一階常係數非齊次線性差分方
Yt-bYt-1=a+I0+DI 方程的一個特解
差分方程
方程的通解為
差分方程
其中A為任意常數。稱係數
差分方程
為凱恩斯乘數。
哈羅德(Harrod.R.H)經濟增長模型
差分方程

  差分方程

設St為t期儲蓄,Yt為t期國民收入,It為t期投資,s稱為邊際儲蓄傾向(即平均儲蓄傾向),0<s<1,k為加速係數。哈羅德宏觀經濟增長模型為:
差分方程
其中s,k為已知常數。
(1)式表示t期儲蓄依賴於前期的國民收入;(2)式表示t期投資為前兩期國民收入差的加速,且預期資本加速係數k為常數;(3)式為均衡條件。
差分方程

  差分方程

經整理后得齊次差分方程
差分方程
其通解為
差分方程
其中A為任意常數,
差分方程

  差分方程

差分方程
,哈羅德稱之為「保證增長率」 其經濟意義就是:如果國民收入Yt按保證增長率
差分方程
增長,那麼就能保證t期儲蓄與t期投資達到動態均衡,即It=St , t=0,1,2,…。
薩繆爾森(Samuelson P.A)乘數加速數模型
設Yt為t期國民收入,Ct為t期消費,It為t期投資,G為政府支出(各期均相同)。薩繆爾森將乘數和加速數兩個參數同時引進而得到國民經濟收支均衡模型(也稱為乘數-加速數模型):
差分方程
其中G>0為常數,b稱為邊際消費傾向(常數),k為加速數。
將(2)(3)兩式代入(1)並經整理后得:Yt-b(1+k)Yt-1+bkYt-2=G.其特解
差分方程
其經濟意義為:國民收入的均衡值等於凱恩斯乘數
差分方程
與政府支出自發投資G的乘積。
對應的齊次方程為 Yt-b(1+k)Yt-1+bkYt-2=0,
其特徵方程為 A2-b(1+k)A+bk=0,特徵方程的判別式
差分方程
差分方程
時,特徵方程有兩相異實根
差分方程
齊次方程的通解為:
差分方程
差分方程
時,特徵方程有一對相等實特徵根
差分方程
齊次方程的通解為:
差分方程
差分方程
時,特徵方程有一對共軛復根:
差分方程
齊次方程的通解為:
差分方程

相關評論

同義詞:暫無同義詞