《经济数学微分方程.ppt》由会员分享,可在线阅读,更多相关《经济数学微分方程.ppt(30页珍藏版)》请在三一文库上搜索。
1、2019/3/16,1,CH6 微分方程,微分方程是研究动态经济学的基本工具。通过计算微分方程来分析变量的具体时间路径,以及能否收敛于均衡。,2019/3/16,2,一、导论 变量为导数的方程称为微分方程。 例如:,如果只有一个自变量,称为常微分方程(ODE)。 常微分方程的阶是方程中最高导数的阶数。 宏观经济学使用的ODE都是对时间的导数。 例:,若x(t)是常数,方程被称为自控的(一个方程仅通过变量y而依赖于时间)。若x(t)0,方程被称为齐次的。,2019/3/16,3,微分方程的解法,求解微分方程的目的在于找到变量的变化特征。 第一种解法:图解法。只能用于自控方程。 第二种解法:解析法
2、。可以找到精确的解,但只能用于有限的函数形式,如线性函数。 第三种解法:数值分析。使用现存软件,如Matlab的子程序ODE23和ODE45。,2019/3/16,4,二、一阶常微分方程的解法,1、图解法 例1:一阶线性自控常微分方程: 其中a和x是常数且大于0。 以y为横轴,以 为纵轴。由于 是y对时间的导数,因此, 为正时,意味着y随着时间的变化而增加; 为负则减少。,2019/3/16,5,图形为直线。 在纵轴的截距为-x,在横轴的截距为-x/a。 在y*点,0,即y不会随时间而变化,y*称为y的稳态值。 当y(0)y*, 0, y随时间而减少。反之则增加。 练习:当a0的动态。,y,y
3、*,当直线的斜率为负时方程是稳定的:无论初始值y(0)在何处,y(t)都将回到y*。,稳态值,2019/3/16,6,例2:非线性函数的动态。 微分方程: 其中s、和都是正常数,且1。,y,y*,0,当 时,y=0;,因此,方程有两个稳态。稳态0是不稳定的,稳态y*是稳定的。,2019/3/16,7,微分方程稳定性总结,对于微分方程: 当 时,可以找出稳态值y*。 若 ,即函数在稳态处的斜率为正,则y是局部不稳定的。 若 ,即函数在稳态处的斜率为负,则y是局部稳定的。,2019/3/16,8,2、解析解,一些简单的微分方程,解为:,解为:,2019/3/16,9,求解线性常微分方程的一般方法,
4、常系数一阶线性微分方程: 求解步骤: 第一:把所有涉及y及其导数的项放在方程的一边,把其余项放在方程的另一边; 第二:两边同乘以积分因子eat并积分,左边就变成了eaty(t)对时间的导数;如果微分方程的系数是可变系数a(t),那么积分因子为 ,左边变成了 的导数。 第三:计算出y(t)。,2019/3/16,10,练习:求解微分方程,例1: 解答:移项、在两边同乘以e-t并积分: 可以解出通解:y(t)=-1+bet 。要想得到特解,需要知道边界条件。 A:如果已知初始条件:t=0时y(t)=0。 则y(0)=-1+be0=0,得到b=1。 特解为: y(t)=-1+et 即曲线A B:如果
5、已知终端条件:t=1000时y为0。 则y(1000)=-1+be1000=0,得到b=e1000。 特解为: y(t)=-1+e1000et 即曲线B,t,y(t),-1,1000,A,B,微分方程的解,2019/3/16,11,练习:人口增长,例2:已知人口增长率为n,计算人口数量的动态变化。 解答: 增长率的定义: 离散形式: 连续形式: 人口增长: 求解微分方程得到:L(t)=L(0)ent,2019/3/16,12,三、线性常微分方程系统,n个微分方程组成的系统: 其中, 、y(t)和x(t)是n维列向量,A是常系数的nn方阵。,用矩阵表示,2019/3/16,13,微分方程系统解法
6、,第一种:相位图。 简单地提供了定性解,但只适用于22系统以及有稳态的自控方程; 第二种:解析解。 第三种:数值法。用时间消去法来计算非线性系统。,2019/3/16,14,1、相位图,(1)对角系统: 以y1为横轴,以y2为纵轴,平面中的每一点都代表了系统(y1,y2)在任一给定时刻的位置。 相位图的目标:把由两个微分方程所隐含的动态转换为一个描述了经济随时间的定性行为的箭头系统。,2019/3/16,15,情形1:a110且a220:系统不稳定。 情形2:a110:鞍点路径稳定。,y1,y2,鞍点路径稳定的相位图,稳定臂,不稳定臂,原点是稳态。鞍点路径既不是稳定又不是不稳定的。系统只有从稳
7、定臂(横轴)开始才会回到稳态。,结论:对角系统的稳定性依赖于系数的符号。若两者都为正,系统不稳定;若两者都为负,系统稳定;若两者异号,系统是鞍点路径稳定。,2019/3/16,16,(2)非对角系统 边界条件为y1(0)=1和 的轨迹是直线y2=0.06y1+1.4 在直线的下方, y20.06y1+1.4, ,即在该区域y1递增;同理,在直线的上方区域y1递减。 的轨迹是直线y1=10 在直线的左边, , y2递增;右边递减。 将 和 联立求解,可以得到稳态值:y1*=10,y2*=2。,2019/3/16,17,具有鞍点路径稳定性的非对角系统的相位图,y1,y2,y1*=10,稳定臂,不稳
8、定臂,y2*=2,稳定臂和不稳定臂对应于两个特征向量,2019/3/16,18,非对角系统稳定性:结论,1、两个特征值是正实数,系统不稳定。 2、两个特征值是负实数,系统稳定。 3、两个特征值是实数但异号,系统是鞍点路径稳定。稳定臂对应于与负特征值相关的特征向量,不稳定臂对应于与正特征值相关的特征向量。 4、两个特征值是有负实部的复数,系统振荡收敛。 5、两个特征值是有正实部的复数,系统振荡且不收敛。 6、两个特征值是有零实部的复数,系统轨迹是环绕稳态运动的椭圆。 7、两个特征值相等,解为y(t)=(b1+b2t)eat。a是特征值,若a0,解是不稳定的。,2019/3/16,19,两变量时的
9、特征根及稳定性,两变量动态系统,2019/3/16,20,(3)非线性系统 解答: 的轨迹为:c=k0.3; 轨迹为k=10。 将k和c的动态结合到一起,系统的稳态是两条轨迹的交点。稳态值:k*=10,c*=2。 系统是鞍点路径稳定的。,2019/3/16,21,k,c,稳定臂,不稳定臂,具有鞍点路径稳定的非线性系统的相位图,2019/3/16,22,2、解析解,(1)线性齐次系统,y(t)是一个n1列向量,A是nn常系数矩阵。,解法:,假设z(t)=V-1y(t) ,则,其中,V是特征向量矩阵,D是特征值对角矩阵。 先解出z(t),然后根据y(t)=Vz(t)可解出y(t)。,2019/3/
10、16,23,(2)线性非齐次系统,解法与齐次系统相同。,练习1:求解以下线性系统的解。,边界条件:y1(0)=1和limte-0.06ty1(t)=0,2019/3/16,24,2019/3/16,25,根据边界条件可以确定b1和b2的值。,根据初始条件y1(0)=1,得到b1+b2=-9,根据终端条件得到:,当t趋于无穷大时,上式中第一和第三项都趋于0,除非b1等于0,否则第二项将趋于无穷大。因此,b1=0。这意味着b2=-9。,这个ODE系统的精确解为:,y1(t)=10-9e-0.04t y2(t)=2-0.9e-0.04t,当t=0时,y1=1,y2=1.1; 当t时,y1*=10,y
11、2*=2。,2019/3/16,26,练习2:求解非线性系统的解,边界条件:y1(0)=1和lim te-0.06ty1(t)=0,非线性系统,稳态值:k*=10,c*=2。围绕稳态值将上述系统进行线性化:,原系统线性化为:,系统与练习1求解方法完全相同。,2019/3/16,27,高维系统,三维非线性系统,根据,得到稳态值,在稳态附近进行泰勒展开,求解系数矩阵A的特征根以判别系统的稳定性,2019/3/16,28,判别方法: 如果特征根都是正数,系统是不稳定的。 如果特征根都是负数,系统是稳定的。 如果负(具有负实部)特征根的个数等于初始条件的个数,则存在惟一收敛路径(稳定臂);如果大于,则有无数个收敛路径,如果小于,则不稳定。 这意味着经济的稳态具有不定性:有可能存在多个收敛路径。 这可以用来解释各国经济增长的差异性。,2019/3/16,29,解释,根据泰勒展开式可以求解:,2019/3/16,30,
链接地址:https://www.31doc.com/p-2281420.html