更新时间:2024-01-29 04:07
包含未知函数的差分及自变数的方程。在求微分方程*的数值解时,常把其中的微分用相应的差分来近似,所导出的方程就是差分方程。通过解差分方程来求微分方程的近似解,是连续问题离散化*的一个例子。
在数学上,递推关系(recurrence relation),也就是差分方程(difference equation),是一种递推地定义一个序列的方程式:序列的每一项目是定义为前一项的函数。某些简单定义的递推关系式可能会表现出非常复杂的(混沌的)性质,他们属于数学中的非线性分析领域。
所谓解一个递推关系式,也就是求其解析解,即关于n的非递归函数。
在数值分析中首先遇到的问题是如何把微分方程化成相应的差分方程 ,使得差分方程的解能最好地近似表示原来的微分方程的解 ,其次才是进行计算。
比如 dy+y*dx=0,y(0)=1 是一个微分方程, x取值[0,1]
(注:解为y(x)=e^(-x));
要实现微分方程的离散化,可以把x的区间分割为许多小区间 [0,1/n],[1/n,2/n],...[(n-1)/n,1]
这样上述微分方程可以离散化为:y((k+1)/n)-y(k/n)+y(k/n)*(1/n)=0, k=0,1,2,...,n-1 (n 个离散方程组)
利用y(0)=1的条件,以及上面的差分方程,就可以计算出 y(k/n) 的近似值了。
性质1 Δk(xn+yn)=Δkxn+Δkyn
性质2 Δk(cxn)=cΔkxn
性质3 Δkxn=∑(-1)jCjkXn+k-j
性质4 数列的通项为n的无限次可导函数,对任意k>=1,存在η,有 Δkxn=f(k)(η)
设{ut,t=0,±1…}为实序列,若满足如下关系式ut-ᵠ1ut-1-…-ᵠput-p=h(t),其中ᵠ1,ᵠ2…,ᵠp为实数,h(t)为t的已知实函数,则称上式为{ut}所满足的线性差分方程。
如将上式中的确定性函数ut,h (t)代之以统计特性已知的随机序列,于是便得到线性随机差分方程。在时间序列分析中并不讨论这样广泛的模型,只涉及一种特殊的线性随机差分方程:
xt-ᵠ1xt-1-…-ᵠpxt-p=εt-θ1εt-1-…-θqεt-g
其中ᵠ1, …,ᵠp, 及θ1, …,θg为实数, {xt}是零均值平稳序列,{εt}是平稳白噪声序列,且当s>t时Eεsxt=0上述特定的线性随机差分方程就是时间序列分析中的ARMA (p,g) 模型。
形如yt+n+a1(t)yt+n-1+a2(t)yt+n-2+…+an-1(t)yt+1+an(t)yt=f(t)的差分方程,称为n阶非齐次线性差分方程。其中a1(t),a2(t),…,an-1(t),an(t)和f(
而形如yt+n+a1(t)yt+n-1+…+an-1(t)yt+1+an(t)yt=0的差分方程,称为n阶齐次线性差分方程。其中ai(t)(i=1,2,…,n)为t的已知函数,且an(
如果ai(t)=ai(i=1,2,…,n)均为常数(an≠0),则有
yt+n+a1yt+n-1+a2yt+n-2+…+an-1yt+1+anyt=f(t),
yt+n+a1yt+n-1+a2yt+n-2+…+an-1yt+1+anyt=0。
分别称为n阶常系数非齐次线性差分方程和n阶常系数齐次线性差分方程。
定理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个任意(独立)常数。
齐次差分方程的通解
将方程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为对应的齐次方程的通解。
存款模型
设St为t期存款总额,i为存款利率,则St与i有如下关系式:
St+1=St+iSt=(1+i)Si, t=0,1,2,…
其中S0为初始存款总额。
动态供需均衡模型(蛛网定理)
用来分析商品当价格和产品失去平衡时,经济状态发生不同波动情况的理论模型。设需求方程:Pt=D(Qdt),即Qdt=c+bpt
供给方程:Qst=S(Pt-1),即Qst=g+hPt-1
其中:P——价格;
D——需求函数;
S——供给函数;
Q——产量;
b——边际需求倾向;
h——边际供给倾向。
现在的供应同前期的价格相关。均衡条件是:
a+bPt=g+hPt-1。
化为一阶线性差分方程的通式:
在正常供求条件下,b<0,h>0,故h/b≠1。解得:
由于b<0,h>0,得h/b<0。故时间轨迹为振荡(其图形参看一阶线性差分方程的通式条目)。
若|h|<|b|,即|h/b|<1,价格变动对供应量的影响小于对需求量的影响,波动逐渐减弱,经济状况趋于均衡。PQ图如(a)所示。
若|h|>|b|,即|h/b|>1,价格变动对供应量的影响大于需求量的影响,时间轨迹为发散。PQ图如(b)所示。
若|h|=|b|,即h/b=-1,时间轨线的振荡型态相同,波动将一直循环下去。PQ图如(C)所示。
三种情况的图形均似蛛网,故该模型称为蛛网模型。
哈罗德(Harrod.R.H)经济增长模型
哈罗德把加速原理引入凯恩斯的投资——储蓄分析,企图说明动态的经济增长。哈罗德模型假设为:
St=SYt,It=a(Yt-Yt-1)
式中:St——t期实际储蓄;
Yt——t期实际收入;
S——平均储蓄倾向为常量;
It——t期计划投资;
a——单位产品需要的投资,也称加速系数。其平衡条件是计划投资等于计划储蓄。即:St=It
有
G为均衡增长系数,即为哈罗德所说的保证增长率。如果收入按这个比率增长,全部储蓄将转化为投资,从而使经济稳定增长。把上式整理后得一阶齐次线性差分方程:,其解为,在正常情况下,故。
发散而不振荡,它意味着收入无限扩展。
萨缪尔森(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。