代数方程与差分方程模型
原创tianguiyuyu 发布于2018-05-29 23:13:18 阅读数 1149 收藏
展开
1 代数方程模型。
所谓的代数方程模型就是有一边变量表示未知量,
代数方程,即由多项式组成的方程。有时也泛指由未知数的代数式所组成的方程,包括整式方程、分式方程和根式方程。
例如:5x+2=7,x=1等。 代数,把algebra翻译成代数,就是用字母代替数的意思,继而推广。随着数学的发展,内在涵义又推广为用群结构或各种结构来代替科学现象中的各种关系。也就是说“代数”本质是个“代”字,通过研究各种抽象结构“代替”直接研究科学现象中的各种关系。
2 差分方程模型
就我个人的观点而言,差分方程模型最重要的作用在于,当我们在解微分方程的时候,有时候微分方程很难直接解,那么这个时候,我们就可以将微分方程的连续化变成离散的。通过找到一个递推式和知道初始条件,那么就可以近似的求解出微分方程的最终解。讲个笑话,高中数学中的等差数列的通项就是差分方程的形式。
在数学上,递推关系(recurrence relation),也就是差分方程(difference equation),是一种递推地定义一个序列的方程式:序列的每一项目是定义为前一项的函数。某些简单定义的递推关系式可能会表现出非常复杂的(混沌的)性质,他们属于数学中的非线性分析领域。
所谓解一个递推关系式,也就是求其解析解,即关于n的非递归函数。
意义性质编辑
意义
在数值分析中首先遇到的问题是如何把微分方程化成相应的差分方程 ,使得差分方程的解能最好地近似表示原来的微分方程的解 ,其次才是进行计算。 [2]
比如 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
差分方程(5张)
的已知实函数,则称上式为{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) 模型。 [3]
形如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(t)都是t的已知函数,且an(t)≠0,f(t)≠0。
而形如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(t)≠0。
如果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图b
其中:P——价格;
D——需求函数;
S——供给函数;
Q——产量;
b——边际需求倾向;
h——边际供给倾向。
现在的供应同前期的价格相关。均衡条件是:图c
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)所示。
三种情况的图形均似蛛网,故该模型称为蛛网模型
————————————————
版权声明:本文为CSDN博主「tianguiyuyu」的原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/tianguiyuyu/article/details/80503590