差分方程模型(讲义).pdf

差分方程模型差分方程模型一一.引言引言数学模型按照离散的方法和连续的方法，可以分为离散模型和连续模型。1.1.确定性连续模型确定性连续模型1)微分法建模(静态优化模型)，如森林救火模型、血管分支模型、最优价格模型。2)微分方程建模(动态模型)，如传染病模型、人口控制与预测模型、经济增长模型。3)稳定性方法建模(平衡与稳定状态模型)，如军备竞赛模型、种群的互相竞争模型、种群的互相依存模型、种群弱肉强食模型。4)变分法建模（动态优化模型），如生产计划的制定模型、国民收入的增长模型、渔业资源的开发模型。2.2.确定性离散模型确定性离散模型1)逻辑方法建模，如效益的合理分配模型、价格的指数模型。2)层次分析法建模，如旅游景点的选择模型、科研成果的综合评价模型。3）图的方法建模，如循环比赛的名次模型、红绿灯的调节模型、化学制品的存放模型。4）差分方程建模，如市场经济中的蛛网模型、交通网络控制模型、借贷模型、养老基金设置模型、人口的预测与控制模型、生物种群的数量模型。随着科学技术的发展，人们将愈来愈多的遇到离散动态系统的问题，差分方程就是建立离散动态系统数学模型的有效方法。在一般情况下，动态连续模型用微分方程方法建立，与此相适应，当时间变量离散化以后，可以用差分方程建立动态离散模型。有些实际问题既可以建立连续模型，又可建立离散模型，究竟采用那种模型应视建模的目的而定。例如，人口模型既可建立连续模型(其中有马尔萨斯模型 Malthus、洛杰斯蒂克 Logistic模型)，又可建立人口差分方程模型。这里讲讲差分方程在建立离散动态系统数学模型的的具体应用。二二.差分方程简介差分方程简介在实际中，许多问题所研究的变量都是离散的形式，所建立的数学模型也是离散的，譬如，像政治、经济和社会等领域中的实际问题。有些时候，即使所建立的数学模型是连续形式，例如像常见的微分方程模型、积分方程模型等。但是，往往都需要用计算机求数值解。这就需要将连续变量在一定的条件下进行离散化，从而将连续型模型转化为离散型模型。因此，最后都归结为求解离散形式的差分方程解的问题。关于差分方程理论和求解方法在数学建模和解决实际问题的过程中起着重要作用。1.1.差分方程的定义差分方程的定义给定一个数列xn,把数列中的前n1项xi(i 0,1,2,n)关联起来得到的方程，则称这个方程为差分方程。差分方程。2.2.常系数线性齐次差分方程常系数线性齐次差分方程常系数线性齐次差分方程的一般形式为xn a1xn1 a2xn2 akxnk 0，(1)或者表示为F(n,xn,xn1,xnk)0 (1)其中k为差分方程的阶数为差分方程的阶数，其中a1,a2,ak为差分方程的系数，且ak 0(k n)。对应的代数方程k a1k1 a2k2 ak 0 (2)称为差分方程(1)的对应的特征方程特征方程。(2)式中的k个根1,2,k称为(1)式的特征根特征根。2.12.1 差分方程的解差分方程的解常系数线性齐次差分方程的解主要是由相应的特征根的不同情况有不同的形式。下面分别就特征根为单根、重根和复根的情况给出方程解的形式。2.1.12.1.1特征根为单根特征根为单根(互不相同的根互不相同的根)设差分方程(1)有k个单特征根(互不相同的根)1,2,k，则nnnxn c11c22 ckk为该差分方程(1)的通解。其中c1,c2,ck为任意常数，且当给定初始条件xi xi(0)，(i 1,2,k)(3)时，可以确定一个特解。例例 1 1 在信道上传输三个字母a,b,c且长度为n的词，规定有两个a连续出现的词不能传输，试确定这个信道允许传输的词的个数。解解:令xn表示允许传输且长度为为n的词的个数，n 1,2,3,，通过简单计算可得x1 3，(a,b,c),x2 8(即 ab,ac,bc,bb,cc,ba,ca,cb)。当n 3时，若词的第一个字母是b或c，则词可按xn1种方式完成;若词的第一个字母是a，则第二个字母是b或c，该词剩下的部分可按xn2种方式完成。于是得差分方程xn 2xn1 2xn2(n 3,4,)其特征方程为2 2 2 0，特征根为113，213则通解为xn c1(13)nc2(13)n，(n 3,4,)利用条件x1 3，x2 8求参数c1，c2，即由c1(13)c2(13)3，22c1(13)c2(13)8解得c1故得到原差分方程的通解为232 3，c2232 3xn232 3(13)n232 3(13)n，(n 1,2,3,4,)2.1.22.1.2特征根为重根特征根为重根设1,2,l是k阶 差 分 方 程xn a1xn1 a2xn2 akxnk 0的l(1 l k)个根，重数分别为m1,m2,ml，且mi k，则该差分方程的通解i1l为xnc1inc2in clini1lni1n1i1n2i1i1i1m1m2ml同样的，有给定的初始条件(3)可以唯一确定一个特解。例例 2 2设初始值为x01,x1 0,x21,x3 2，解差分方程xn xn13xn25xn3 2xn4 0，(n 4,5,)解解:该差分方程的特征方程为43325 2 0，解得其根为1,1,1,2，故通解为xn c1(1)nc2n(1)nc3n2(1)nc42n代入初始条件x01,x1 0,x21,x3 2，得c14229710，c1，c3，c452525252故该差分方程的满足初始条件的解为xn2.1.32.1.3特征根为复根特征根为复根4229710(1)nn(1)nn2(1)n2n52525252设k阶 差 分 方 程xn a1xn1 a2xn2 akxnk 0的 一 对 共 轭 复 根1,2i和相异的k 2个单根3,4,k，则该差分方程的通解为nnxn c1ncosn c2nsin n c33 c4n4 ckk其中22，arctan。同样由给定的初始条件(3)可以唯一确定一个特解。另外，对于有多个共轭复根和相异实根，或共轭复根和重根的情况，都可类似的给出差分方程解的形式。3.3.常系数线性非齐次差分方程常系数线性非齐次差分方程常系数线性非齐次差分方程的一般形式为xn a1xn1 a2xn2 akxnk f(n)(4)其中k为差分方程的阶数为差分方程的阶数，其中a1,a2,ak为差分方程的系数，且ak 0(k n)，f(n)为已知函数。在差分方程(4)中，令f(n)0，所得方程xn a1xn1 a2xn2 akxnk 0 (5)称为非齐次差分方程(4)(4)对应的齐次差分方程对应的齐次差分方程，即与差分方程(1)的形式相同。求解非齐次差分方程通解的一般方法：*首先求对应的齐次差分方程首先求对应的齐次差分方程(5)(5)的通解的通解xn，然后求非齐次差分方程，然后求非齐次差分方程(4)(4)的一的一(0)个特解个特解xn，则*(0)xn xn xn为非齐次差分方程(4)的通解。*关于求xn的方法同求差分方程(1)的方法相同。对于求非齐次方程(4)的特(0)解xn的方法，可以用观察法确定，也可以根据f(n)的特性用待定系数法确定，具体方法可参照常系数线性非齐次微分方程求特解的方法。4.4.差分方程的平衡点及其稳定性差分方程的平衡点及其稳定性在应用差分方程研究问题时，一般不需要求出方程的通解，在给定初值后，通常可用计算机迭代求解，但常常需要讨论解的稳定性。对于差分方程F(n,xn,xn1,xnk)0，若有常数a是其解，即有F(n,a,a,a)0则称a是差分方程F(n,xn,xn1,xnk)0的平衡点平衡点，又对该差分方程的任意由初始条件确定的解xn x(n)，均有lim x ann则称这个平衡点a是稳定的稳定的；否则是不稳定的不稳定的。下面给出一些特殊差分方程的平衡点和稳定性。4.14.1一阶常系数线性差分方程一阶常系数线性差分方程一阶常系数线性差分方程的一般形式为xn1 axn b，(6)其中a,b为常数，且a 1,0。它的通解为xn C(a)nba 1 (7)易知ba1是方程(6)的平衡点，由(7)式知，当且仅当a 1时，ba1是方程(6)的稳定的平衡点。4.24.2二阶常系数线性差分方程二阶常系数线性差分方程二阶常系数线性差分方程的一般形式为xn2 axn1bxn r，(8)其中a,b,r为常数，当r 0时，它有一特解x*0，当r 0，且a b1 0时，它有一特解x*ra b 1，不管是哪种情形，x*是方程(8)的平衡点。设方程(8)的特征方程为2 ab 0的两个根分别为1，2，则当1,2是两个不同的实根时，方程(8)的通解为xn x*C1(1)n C2(n2)；当12是两个相同实根时，方程(8)的通解为xn x*(C1 C2n)n当1,2(cos isin)是一对共轭复根时，方程(8)的通解为xn x*n(C1cosn C2sin n)易知，当且仅当特征方程的任一特征根i1时，平衡点x*是稳定的。4.34.3一阶非线性差分方程一阶非线性差分方程一阶非线性差分方程的一般形式为xn1 f(xn)(9)其平衡点x*由代数方程x f(x)解出。为了分析平衡点x*的稳定性，将方程(9)的右端f(xn)在x*点作泰勒展开，只取一次项，得到xn1 f(x*)(xn x*)f(x*)(10)(10)是(9)的近似线性方程，x*是(10)的平衡点，根据一阶常系数线性差分方程(6)xn1 axn b的稳定性判定的相关结论，得：当f(x*)1时，方程(9)的平衡点是稳定的；当f(x*)1时，方程(9)的平衡点是不稳定的。三三 差分方程建模实例差分方程建模实例1 1 贷款买房问题贷款买房问题某居民买房向银行贷款 6 万元，利息为月利率 1%，贷款期为 25 年，要求建立数学模型解决如下问题：1)问该居民每月应定额偿还多少钱？2)假设此居民每月可节余 700 元，是否可以去买房?1.11.1确定参变量：确定参变量：用n表示月份，An表示第n个月欠银行的钱，r表示月利率，x表示每月还钱数，A0表示贷款额。1.21.2模型的建立与求解模型的建立与求解1)1)模型的建立模型的建立时间初始一个月后二个月后三个月后欠银行款A0A1 A0(1 r)xA2 A1(1 r)xA3 A2(1 r)xAn An1(1 r)xn个月后由上表可得相邻两个月的递推关系式An An1(1 r)x1.31.3模型的求解：模型的求解：(1)差分方程求解方法先求其特解。令An An1 y，则y y(1 r)x，得特解为y 再求对应齐次方程An An1(1 r)的通解。对应的特征方程为x。r(1r)0，得(1r)。齐次方程的通解为：c(1r)n因此原方程的通解为：An c(1r)nxr又因为n 0时An A0，得c A0故xrn1r1An A01r xnr(2)递推法：An A0(1 r)x11 rnn1 r11 r A01 r xn1nr令A0=60000，A300 0，n=300，r=0.01得x A01 rrn1 rn1600001 0.010.013001 0.013001 632元因此，该居民每月应偿还 632 元。又 632700，所以该居民可以去买房。2 2借贷问题借贷问题中国建设银行北京市分行个人住房贷款一至二十年“月均还款金额表”（自1998 年 3 月 25 日起执行）的一部分如下：(借款额为一万元)单位：元贷款期限(年)1520年利率(%)10.20610.206还款总额(元)19569.6023488.80利息负担总和(元)9569.6013488.80月均还款额(元)108.7297.87试问他们是怎样算出来的？借贷问题的数学模型借贷问题的数学模型一一.符号说明符号说明以贷款期限 20 年为例:借贷额-A010,000;贷款期限-为 N 年;月利率-r 10.206/12 0.008505;“月均还款额”-表示每月还款额是相同的，记为x;还款总额-记为S.二二.建立模型建立模型一开始借款A010,000，一个月后欠银行本利为A1 A0(1 r)，但为了减少欠款，还了x元，因而A1 A0(1 r)x，第k个月情况也是这样的，即Ak Ak1(1r)x,k 1,2,N注意到了第 N 个月已经不欠银行的钱了，即AN 0，因此，我们得到以下的数学模型:Ak Ak1(1 r)xk 1,2,NA0,x,N KnownFind out such that AN 0三三.数学模型的求解数学模型的求解首先求出用已知量表出的表达式。由A2 A1(1 r)x A0(1 r)x(1 r)x A0(1 r)2 x1(1 r)可以猜想，并用数学归纳法证明:Ak A0(1 r)k x1(1 r)(1 r)2 (1 r)k1由等比数列前k 1项的求和公式知:xAk A0(1 r)k(1 r)k1,k 1,2,Nr再由AN 0，得到:A0r(1 r)Nx(1 r)N1把已知量带入，就得到表中的x。3 3生物种群数量问题生物种群数量问题一问题的提出一问题的提出种群的数量问题是当前世界上引起普遍关注的一个问题。要预测未来种群的数量，最重要的影响因素是当前的种群数量，今后一段时间内种群的增长状况和环境因素。由于随着种群数量增加到一定的程度后，种群在有限的生存空间进行竞争，种群的增长状况会随着种群数量的增加而减少，而且在有限的生存空间，种群数量也不可能无限增长，假设只能达到某一固定的数量值记为xm，称为最大种群容量。又假设单位时间内种群数量的增长量与当时种群数量x的比记为：r(x)r sx,r、s 0，其中r相当于x0时的增长率，称为固有增长率，记当前(即t0时)种群数量为x0，时刻t种群数量为x(t)。若利用统计数据可知xm，r，x0，则1 1）设x(t)为连续、可微函数，请给出未来时间里种群数量满足的数学模型。2 2）由于某些种群是在固定的一段时间内进行繁殖，所以可用种群繁殖周期作为时间段来研究其增长状况。请给出未来时间里这类种群数量应满足的离散数学模型。二二.问题分析与模型建立问题分析与模型建立 1.1.由于r(x)为单位时间内种群数量的增长量与当时种群数量的比，所以t到t t时间内种群数量的增量为x(t t)x(t)r(x)x(t)t (1)又由于r(x)r sx,而当x xm时增长率应为零，即r(xm)0，所以s 则r(x)r r，xmrx，xm把它代入方程(1)得:x(t t)x(t)(r r)x(t)t (2)xm此方程两边同除t，并令t 0，加上初始条件x(0)x0可得未来任意时刻t种群数量所满足的数学模型为：dxx 1 rx (3)dtxmx(0)x0 2.2.由于是利用种群繁殖周期作为时段来研究种群增长状况，则令t 1，t视为整数及r(x)r rx代入方程(1)得：xmr)x(t)(4)xmx(t 1)x(t)(r 加上初始条件x(0)x0得任意时刻t种群数量所满足的离散型数学模型为rx(t 1)(r 1)x(t)xmx(0)x0通过这个差分方程就可以很容易得到任意时刻t种群的数量。三模型求解三模型求解 1 1利用Mathematica求解方程(1)，可得任意时刻t种群数量为x(t)xm xmrt11xe0Mathematica源程序为：DSolvex(t)r*(1 xt/xm)*xt 0,xt,t 2 2 根据方程(2)，只要给出初值x0就可以很容易进行递推而得到任意时刻t种群的数量。四结果分析四结果分析 1 1上面方程(3)有时称为阻滞增长模型或Logistic模型，它有着广泛的应用。例如传染病在封闭地区的传播，耐用消费品在有限的市场上的销售等现象，都可以合理的、简化的用这个模型来进行描述。但它存在不足，因为随着环境的变迁，最大种群容量可能会发生变化，而且最大种群容量也不容易准确得到。2 2一方面，用离散化的时间来研究问题有时是很方便的，尤其出现了计算机以后，人们可以很方便的对问题进行求解；另一方面，对这个种群数量问题，由于许多种群实际上是由单一世代构成的，在相继的世代之间几乎没有重叠，所以种群的增长是分步进行的。这种情况下，为了准确的描述种群的数量动态就不能用微分方程，而应利用离散的模型来描述。4.4.人口的控制与预测模型人口的控制与预测模型一问题的提出一问题的提出常见的两个常微分方程模型(马马尔尔萨萨斯斯(Malthus)(Malthus)模模型型和和洛洛杰杰斯斯蒂蒂克克(Logistic)(Logistic)模型模型)没有考虑到社会成员之间的个体差异，即不同年龄、不同体质的人在死亡、生育方面存在的差异。完全忽略了这些差异显然是不合理的。但我们不可能对每一个人的情况逐个加以考虑，故仅考虑年龄的差异对人口的变动的影响，即假设同一年龄的人具有相同的死亡率和生育能力，这样建立的模型不但使我们能够更细致的预测人口总数，而且能够预测老年人口、劳动力人口、学龄人口等不同年龄组的人口信息.下面来建立离散的差分数学模型来表现人口数量的变化规律。二模型的建立与求解设xk(t)为第t年年龄为k的人口数量，k 0,1,2,100，即忽略百岁以上的人口。如果知道了第t年各年龄组的人口数，各年龄组人口的生育及死亡状态，就可以根据人口发展变化规律推得第t 1年各年龄组的人口数。首先引入k岁人口的死亡率和k岁育龄妇女的年生育率这两个概念，他们的含义和记号如下：k岁人口的年死亡率：dk一年内k岁的死亡人数这年内k岁的人口数k岁妇女的年生育率：bk一年内k岁妇女生育的婴儿数这年内k岁妇女人数第t 1年k 1岁的人口数就是第t年k岁人口数扣除它在该年的死亡人数，即xk1(t 1)(1 dk)xk(t)，令pk1 dk称为k岁人口的存活率，故各年龄组人口随时间的变化规律可用递推公式xk1(t 1)pkxk(t),(k 0,1,99)来表示。再考虑到零岁的人数x0(t 1)bkuk(t)xk(t)，k0100其中uk(t)xk(t)为第t年k岁的妇女人数，uk(t)为第t年k岁人口的女性比(占全部k岁人口数)，bkuk(t)xk(t)就是第t年k岁妇女所生育的婴儿数.由此得到的人口模型是：100 x0(t 1)bkuk(t)xk(t)(1)k0 x(t 1)p x(t),k 0,1,99kkk1根据人的生理特征和人口学中的习惯，妇女的育龄区间一般取为 15 岁至 49 岁之间，即当k 15和k 49时，bk 0，令x(t)(x0(t),x1(t),xk(t),x100(t)Tu0(t)b0p00L 0u1(t)b10p10u2(t)b2000u99(t)b9900p99u100(t)b100000则人口模型（1）的矩阵形式为x(t 1)Lx(t)(2)其中L称为莱斯利（Lwslie）矩阵.当第t0年的人口状况已知时，从式(2)就可以推得第t年的人口为x(t 1)Ltt0 x(t0).5.5.市场经济中的蛛网模型市场经济中的蛛网模型在自由竞争的市场经济中，商品的价格是由市场上该商品的供应量决定的，供应量越大，价格就越低。另一方面，生产者提供的商品数量又是由该商品的价格决定的，价格上升将刺激生产者的生产积极性，导致商品生产量的增加。反之，价格降低会影响生产者的积极性，导致商品生产量的下降。在没有外界干扰的情况下，这种现象将如此反复下去。这样的需求和供应关系决定了市场经济中商品的价格和数量必然是振荡的。这种振荡越小越好，如果振荡太大就会影响人民群众的正常生活。供大于求供大于求价格下降价格下降数量和价格在振荡数量和价格在振荡产量减少产量减少产量增加产量增加价格上涨价格上涨供不应求供不应求(1)商品数量与价格的振荡在什么条件下趋向稳定?(2)当不稳定时政府能采取什么干预手段使之稳定?下面用差分方程理论建模，讨论市场经济趋于稳定的条件，再用图形方法建立“蛛网模型”对上述现象进行分析，对结果进行解释，然后作适当推广。3.13.1 模型的假设和符号说明模型的假设和符号说明 记第n时段商品数量为xn，价格为yn，n 1,2,。这里我们把时间离散化为时段，1 个时段相当于商品的 1 个生产周期，如蔬菜、水果可以是 1 年，肉类可以是一个饲养周期。在n时段商品的价格yn取决于数量xn。设yn f(xn)。它反映消费者对这种商品的需求关系，称为需求函数需求函数。因为商品的数量越多，价格越低。需求函数在图 1 中用一条下降的曲线f表示，f称为需求曲线。在n1时段商品的数量xn1由上一时段的价格yn决定，用xn1 g(yn)表示。它反映生产者的供应关系，称为供应函数供应函数。因为价格越高，生产量越大。供应函数在图 1 中用一条上升的曲线g表示，g称为供应曲线。y yf fg gy y0 0P P0 0O Ox x0 0 x x图图 1 1商品供求关系曲线商品供求关系曲线3.23.2模型的建立与求解模型的建立与求解设需求曲线f和供应曲线g相交于点P0(x0,y0)，在P0附近取函数f和g的线性近似，即需求曲线需求曲线f：yn y0(xn x0)，0 (11)供应曲线供应曲线g：xn1 x0(yn y0)，0 (12)由式(11)(12)消去yn，得到一阶线性差分方程xn1 xn(1)x0，n 1,2,(13)因此x0是其平衡点，即P0是平衡点。对式(13)进行递推，得xn1()nx11()nx0，n 1,2,由此可得，平衡点稳定的条件是：1；不稳定的条件是：1。下面用图形解释此模型。若对某一个k有xk x0，则由(11)式得，当n k时xn x0，从而yn y0，即商品的数量和价格将永远保持在P0(x0,y0)点。但是实际生活中的种种干扰使得xn,yn不可能停止在P0(x0,y0)上。不妨设x1偏离x0(见图 2，图 3)，我们来分析随着n的增加，xn,yn的变化情况。y yf fy y2 2P P3 3y y0 0y y3 3y y1 10 0g gP P4 4f f需求曲线需求曲线g g供应曲线供应曲线P P0 0P P2 2x x2 2x x0 0 x x3 3P P1 1x x1 1x x图图 2 2P0点是稳定的点是稳定的数量x1给定后，价格y1由曲线f上的P1点决定，下一时段的数量x2由曲线g上的P2点决定，这样得到一序列的点P1(x1,y1)，P2(x2,y2)，P3(x3,y3)，P4(x4,y4)，,在图 2 上，这些点将按照箭头所示方向趋向P0(x0,y0)，表明P0(x0,y0)是稳定的平衡点，意味着市场经济(商品的数量和价格)将趋向稳定。但是如果需求函数和供应函数由图 3 的曲线所示，则类似的分析发现，市场将按照P1(x1,y1)，P2(x2,y2)，P3(x3,y3)，P4(x4,y4)，的规律变化为远离P0(x0,y0)，即P0(x0,y0)是不稳定的平衡点，市场经济趋向不稳定。y yP P3 3P P4 4f f需求曲线需求曲线y y0 0g g0 0P P2 2P P0 0P P1 1f fx x0 0图图 3 3P0点是不稳定的点是不稳定的x x图 2 和图 3 中折线P1P2P3P4形似蛛网，于是这种用需求曲线和供应曲线分析市场经济稳定性的图示法在经济学中被称为蛛网模型蛛网模型。实际上，需求曲线f和供应曲线g的具体形式通常是根据各个时段商品的数量和价格的一系列统计资料得到的。一般地说，f取决于消费者对这种商品地需要程度和他们地消费水平，g则与生产者的生产能力，经营水平等因素有关。下面来解释此模型的实际意义。首先来考虑参数首先来考虑参数,的含义。的含义。需求函数需求函数f的斜率的斜率(取绝对值取绝对值)：表示商品供应量减少 1 个单位时价格的上涨幅度；供应函数供应函数g的斜率的斜率：表示价格上涨 1 个单位时(下一时期)商品供应增加量。的值反映消费者对商品需求的敏感程度。如果这种商品是生活必需品，消费者处于持币待购状态，商品数量稍缺，人们立即蜂拥购买，那么会比较大；反之，若这种商品非必需品，消费者购物心理稳定，或者消费水平低下，则会比较小。的数值反映生产经营者对商品价格的敏感程度。如果他们目光短浅，热衷于追逐一时的高利润，价格稍有上涨立即大量增加生产，那么会比较大；反之，若他们目光长远，则会比较小。根据根据,的意义很容易对市场经济稳定与否的条件作出解释。的意义很容易对市场经济稳定与否的条件作出解释。当供应函数g的斜率固定时，越小，需求曲线越平，表明消费者对商品需求的敏感程度越小，越有利于经济稳定。当需求函数f的斜率固定时，越小，供应曲线越陡，表明生产者对价格的敏感程度越小，越有利于经济稳定。反之，当,较大，表明消费者对商品的需求和生产者对商品的价格都很敏感，则会导致经济不稳定。经济不稳定的解决方案经济不稳定的解决方案当市场经济趋向不稳定时，政府有两种干预办法：一种办法是控制价格，无论商品数量多少，命令价格不得改变，于是 0；不管曲线g如何，总是稳定的；另一种办法是控制市场上的商品数量，当上市量小于需求时，政府从外地收购或调拨，投入市场，当上市量多于需求时，政府收购过剩部分，于是 0，不管曲线f如何，也总是稳定的。3.33.3模型的改进和推广模型的改进和推广如果生产者的管理水平更高一些，他们再决定商品生产数量时，不是仅根据前一时期的价格，而是根据前两个时期的价格，为简单起见不妨设根据二者的平均值yn yn12于是供应函数为xn1 g(yn yn1)2在P0点附近取线性近似时，式(12)表示为供应函数(g):xn1 x0(yn yn1 y0)，0 (14)2又设需求函数仍由式(11)表示，则由(11),(14)得到2xn2xn1xn(1)x0，n 1,2,(15)(15)式是二阶线性差分方程。P0点稳定的条件可由特征方程22 0的根1,2()28确定。4结论：若方程的特征根均在单位园内，即11,21，则P0为稳定点。当 8时，显然有()282 2，44从而2 2，故此时P0是不稳定的。当 8时，特征方程有两个共轭复数根1,2此时1i8()2441,21()(8()2)2244221要使P0为稳定点，只需1,21，即有 2这与原有模型中P0点稳定的条件1相比，保持经济稳定的参数,的范围放大了(,的含义未变)。可以想到，这是生产经营者的生产管理水平提高，对市场经济稳定起着有利影响的必然结果。专题训练题：专题训练题：养老金计划养老金计划养老金是指人们在年老失去工作能力后可以按期领取的补偿金，这里假定养老金计划从 20 岁开始至 80 岁结束，年利率为 10%。参加者的责任是，未退休时（60 岁以前）每月初存入一定的金额，其中具体的存款方式为：20 岁29 岁每月存入X1元，30 岁39 岁每月存入X2元，40 岁49 岁每月存入X3元，50 岁59 岁每月存入X4元。参加者的权利是，从退休（60 岁）开始，每月初领取退休金P，一直领取 20 年。试建立养老金计划的数学模型，并计算下列不同年龄的计划参加者的月退休金。(1)从 20 岁开始参加养老金计划，假设X1 X2 X3 X4 200元；X41000(2)从35岁开始参加养老金计划，假设X2 200元，X3 500元，元；(3)从 48 岁开始参加养老金计划，假设X31000元，X4 2000元。