差分模型精选PPT.ppt

《差分模型精选PPT.ppt》由会员分享，可在线阅读，更多相关《差分模型精选PPT.ppt（45页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、关于差分模型第1页，讲稿共45张，创作于星期日 以以t 表示时间，规表示时间，规 定定t只取非负整数。只取非负整数。t=0表示第一周期初，表示第一周期初，t=1表示第二周期初。表示第二周期初。记记 为变量为变量y在时刻在时刻t 时的取值，则称时的取值，则称 为为 的的一阶差分一阶差分，称，称 为的为的二阶差分二阶差分。第2页，讲稿共45张，创作于星期日一般地，函数 的n-1阶差分的差分称为n阶差分阶差分，记为 ，即 例例1设 求第3页，讲稿共45张，创作于星期日 差分的性质：（1）（2）（3）（4）例例2求 的差分第4页，讲稿共45张，创作于星期日 二、差分方程的概念二、差分方程的概念含有未知

2、函数含有未知函数 的差分的方程为差分方程。的差分的方程为差分方程。差分方程的一般形式差分方程的一般形式:其中含的最高阶差分的阶数称为该差分方程其中含的最高阶差分的阶数称为该差分方程的的阶阶。例如，二阶差分方程。例如，二阶差分方程第5页，讲稿共45张，创作于星期日也可改写成也可改写成一般地一般地第6页，讲稿共45张，创作于星期日满足差分方程的序列满足差分方程的序列yt称为此差分方程的称为此差分方程的解解。类似于微分方。类似于微分方程情况，若解中含有的独立常数的个数等于差分方程的阶数程情况，若解中含有的独立常数的个数等于差分方程的阶数时，称此解为该差分方程的时，称此解为该差分方程的通解通解。若解中

3、不含任意常数，则。若解中不含任意常数，则称此解为满足某些初值条件的称此解为满足某些初值条件的 特解特解，例如，考察两阶差分方，例如，考察两阶差分方程程 易易见见与与均是它的特解，而均是它的特解，而则为则为它的通解，其它的通解，其中中c1，c2为为两个两个任意常数。任意常数。类类似于微分方程，称差分方程似于微分方程，称差分方程 为为n阶线性差分方程，阶线性差分方程，当当 0时称其为时称其为n阶非齐次线性差分阶非齐次线性差分方程，而方程，而 第7页，讲稿共45张，创作于星期日则则被称被称为为方程方程对应对应的的齐齐次次线线性差分方程性差分方程。若所有的若所有的 ai(t)均为与均为与t 无关的常数

4、，则称其为无关的常数，则称其为 常系数差分方程常系数差分方程，即即n 阶常系数线性差分方程可分成阶常系数线性差分方程可分成（1）的形式，其对应的齐次方程为的形式，其对应的齐次方程为（2）容易容易证证明，若序列明，若序列与与均均为为方程（方程（2）的解，）的解，则则也是方程（也是方程（2）的解，其）的解，其中中c1、c2为为任意常数，任意常数，这说这说明，明，齐齐次方程的解构成一个次方程的解构成一个线线性空性空间间（解空（解空间间）。）。此规律对于（此规律对于（2）也成立。）也成立。第8页，讲稿共45张，创作于星期日例例3 试改变差分方程 的形式.例例4试确定下列差分方程的阶.例例5指出下列等式

5、哪一个是差分方程,若是,进一步指出是否为线性方程.第9页，讲稿共45张，创作于星期日第10页，讲稿共45张，创作于星期日方程（方程（1）可用如下的代数方法求其通解：）可用如下的代数方法求其通解：（步一步一）先求解）先求解对应对应的特征方程的特征方程 （3）（步二步二）根据特征根的不同情况，求）根据特征根的不同情况，求齐齐次方次方程程(2)的通解的通解情况情况1若特征方程（若特征方程（2）有）有n个互不相同的个互不相同的实实根根,则齐则齐次方程（次方程（2）的通解）的通解为为(C1,Cn为为任意常数任意常数)，情况情况2若若是特征方程（是特征方程（3）的）的k重根，通解中重根，通解中对应对应于于

6、的的项为项为为为任意常数，任意常数，i=1,k。情况情况3 若特征方程（若特征方程（3）有单重复根）有单重复根 通解中对应它们的项为通解中对应它们的项为 为为的模，的模，为为的幅角。的幅角。第11页，讲稿共45张，创作于星期日情况情况4若若为为特征方程（特征方程（3）的）的k重复根，重复根，则则通通解解对应对应于它于它们们的的项为项为为为任意常数，任意常数，i=1,2k。.若若yt为为方程方程(2)的通的通解解,则则非非齐齐次方程次方程(1)的通解的通解为为（步三步三）求非求非齐齐次方程次方程(1)的一个特解的一个特解第12页，讲稿共45张，创作于星期日例例1求解两求解两阶阶差分方程差分方程解

7、解对应齐对应齐次方程特征方程次方程特征方程:,其特征根其特征根为为，对应齐对应齐次方程的通解次方程的通解为为 原方程有形如原方程有形如的特解。代入原方程求得的特解。代入原方程求得，故原方程的通解，故原方程的通解为为在应用差分方程研究问题时，一般不需要求出方程的通解，在应用差分方程研究问题时，一般不需要求出方程的通解，在给定初值后，通常可用在给定初值后，通常可用 计算机迭代计算机迭代求解，但我们常常需要求解，但我们常常需要讨论解的稳定性。对讨论解的稳定性。对 差分方程差分方程(1),若不论其对应齐次方程的通若不论其对应齐次方程的通解中任意常解中任意常 数数C1,Cn如何取值如何取值,在在 时总有

8、时总有 ,则称方程则称方程(1)的解是稳定的解是稳定 的的,否则称其解为不稳定否则称其解为不稳定 的的.根据通根据通解的结构不难看出解的结构不难看出,非齐次方程非齐次方程(1)稳定的充要条件为其所有特稳定的充要条件为其所有特征根的模均小于征根的模均小于1。第13页，讲稿共45张，创作于星期日n rsolve(方程，解函数,选项)，rsolve(方程组，初始条件，解函数，选项)选项为genfunc(x)解以x为自变量；选项为makeproc解为过程函数。rsolve(F(n)=F(n-1)+F(n-2),F(1.2)=1,F,genfunc(x);rsolve(s(n)=2*s(n-1),s(0

9、)=1,s,makeproc)第14页，讲稿共45张，创作于星期日第15页，讲稿共45张，创作于星期日MatlabMatlab求解求解n濒危物种的自然演变和人工孵化n问题 Florida沙丘鹤属于濒危物种，它在较好自然环境下，年均增长率仅为1.94%，而在中等和较差环境下年均增长率分别为-3.24%和 -3.82%，如果在某自然保护区内开始有100只鹤，建立描述其数量变化规律的模型，并作 数值计算。第16页，讲稿共45张，创作于星期日模型建立模型建立n xk+1=(1+r)xk k=0,1,2n记第k年沙丘鹤的数量为xk,年均增长率为r，则第k+1年鹤的数量为n已知x0=100,在较好，中等和

10、较差的自然环境下 r=0.0194,-0.0324,和-0.0382 我们利用Matlab编程，递推20年后观察沙丘鹤的数量变化情况第17页，讲稿共45张，创作于星期日Matlab实现实现n首先建立一个关于变量n，r的m函数nfunction x=sqh(n,r)na=1+r;nx=100;nfor k=1:nn x(k+1)=a*x(k);nend第18页，讲稿共45张，创作于星期日n在command窗口里调用sqh函数 k=(0:20);y1=sqh(20,0.0194);y2=sqh(20,-0.0324);y3=sqh(20,-0.0382);round(k,y1,y2,y3)第19页

11、，讲稿共45张，创作于星期日利用plot 绘图观察数量变化趋势n可以用不同线型和颜色绘图nr g b c m y k w 分别表示 红绿兰兰绿洋红黄黑白色：+o*.X s d 表示不同的线型 第20页，讲稿共45张，创作于星期日n plot(k,y1,k,y2,k,y3)在同一坐标系下画图 plot(k,y2,:k)离散点：黑色 plot(k,y2,-)plot(k,y2,r)plot(k,y2,y)plot(k,y2,y,k,y1,:)plot(k,y2,k,y1,:)plot(k,y2,oy,k,y1,:)用gtext(r=0.0194),gtext(r=-0.0324),gtext(r=

12、-0.0382)在图上做标记。第21页，讲稿共45张，创作于星期日 筹措教育经费模型筹措教育经费模型某家庭从现在着手,从每月工资中拿出一部分资金存入银行,用于投资子女的教育,并计算20年后开始从投资账户中每月支取1000元,直到10年后子女大学毕业并用完全部资金.要实现这个投资目标,20年内要总共筹措多少资金?每月要在银行存入多少钱?假设投资的月利率为0.5%,为此,设第t个月,投资账户资金为每月存资金为b元,于 是2 0年 后,关 于 的 差 分 方 程 模 型 为 第22页，讲稿共45张，创作于星期日例例（市场经济的蛛网模型市场经济的蛛网模型）在自由竞争的市场经济中，商品的价格是由市场上该

13、在自由竞争的市场经济中，商品的价格是由市场上该商品的供应量决定的，供应量越大，价格就越低。另商品的供应量决定的，供应量越大，价格就越低。另一方面，生产者提供的商品数量又是由该商品的价格一方面，生产者提供的商品数量又是由该商品的价格决定的，价格上升将刺激生产者的生产积极性，导致决定的，价格上升将刺激生产者的生产积极性，导致商品生产量的增加。反之，价格降低会影响生产者的商品生产量的增加。反之，价格降低会影响生产者的积极性，导致商品生产量的下降。积极性，导致商品生产量的下降。在市场经济中，对每一商品事实上存在着两个不同的在市场经济中，对每一商品事实上存在着两个不同的函数：函数：（1）供应函数）供应函

14、数x=f(P)，它是价格它是价格P的单增函数，其曲的单增函数，其曲线称为供应曲线。线称为供应曲线。（2）需求函数）需求函数x=g(P)，它是价格它是价格P的单降函数，其曲的单降函数，其曲线称为需求曲线，供应曲线与需求曲线的线称为需求曲线，供应曲线与需求曲线的 形状形状如图所示。如图所示。第23页，讲稿共45张，创作于星期日记记t时时段初市段初市场场上的供上的供应应量量(即上即上一一时时段的生段的生产产量量)为为xt，市，市场场上上该该商品的价格商品的价格为为Pt。商品成交的。商品成交的价格是由需求曲价格是由需求曲线线决定的，决定的，即即随着随着,Mt将将趋趋于平衡点于平衡点M*,即商品量将即商

15、品量将趋趋于平衡于平衡量量x*,价价格将格将趋趋于平衡价于平衡价格格P*。图图中的箭中的箭线线反映了在市反映了在市场经济场经济下下该该商品的商品的供供应应量与价格的量与价格的发发展展趋势趋势。xoPP0P2P*P1xx1x2x0 x*需求曲线需求曲线供应曲线供应曲线M0M2M1M*PoM3M2M1第24页，讲稿共45张，创作于星期日图图和图和图的区别在哪里，的区别在哪里，如何判定平衡点的稳定如何判定平衡点的稳定 性呢？性呢？但是，如果供应曲线和需求曲线呈但是，如果供应曲线和需求曲线呈 图图中的形状，则平衡点中的形状，则平衡点M*是不稳定的，是不稳定的，Mt将越来越远离平衡点。即使初始时刻的供将

16、越来越远离平衡点。即使初始时刻的供应量和价格对应于平衡点，一点微小的波动也会导致市场供应量和价格对应于平衡点，一点微小的波动也会导致市场供求出现越来越大的混乱。上述用图示法分析市场经济稳定性求出现越来越大的混乱。上述用图示法分析市场经济稳定性的讨论在经济学中被称为市场经济的的讨论在经济学中被称为市场经济的 蛛网模型蛛网模型。不难看出，在不难看出，在 图图中平衡点中平衡点M*处处供应曲线的切线斜率大于供应曲线的切线斜率大于需求曲线切线斜率的绝对值，需求曲线切线斜率的绝对值，而在图而在图中情况恰好相反。中情况恰好相反。第25页，讲稿共45张，创作于星期日现在利用差现在利用差 分方程方法来研究蛛网模

17、型，以验证上述猜测是否正确。分方程方法来研究蛛网模型，以验证上述猜测是否正确。我们知道，平衡点我们知道，平衡点M*是否稳定取决于在是否稳定取决于在M*附近供、需曲线的局部附近供、需曲线的局部性态。为此，性态。为此，用用M*处供、需曲线的线性近似来代替它们，并讨论此处供、需曲线的线性近似来代替它们，并讨论此线性近似模型线性近似模型 中中M*的稳定性。的稳定性。设设供供应应曲曲线线与需求曲与需求曲线线的的线线性近似分性近似分别为别为 和和式中，式中，a、b分分别为别为供供应应曲曲线线在在M*处处的切的切线线斜率与需求曲斜率与需求曲线线在在M*处处切切线线斜率的斜率的绝对值绝对值。根据市根据市场经济

18、场经济的的规规律，当供律，当供应应量量为为xt时时，现时现时段的价格段的价格，又，又对对价格价格，由供，由供应应曲曲线线解得下一解得下一时时段的商品量段的商品量 第26页，讲稿共45张，创作于星期日由此由此导导出一出一阶阶差分方程：差分方程：（4）此差分方程的解在此差分方程的解在(b/a)b时，顾客需求对价格的敏感度较小（小于生产时，顾客需求对价格的敏感度较小（小于生产者的敏感程度）商品供应量和价格会自行调节而逐步趋于稳者的敏感程度）商品供应量和价格会自行调节而逐步趋于稳定；反之若定；反之若ab（商品紧缺易引起顾客抢购）该商品供售市（商品紧缺易引起顾客抢购）该商品供售市场易造成混乱场易造成混乱

19、Remark:需求曲线越平,供应曲线越陡，越有利于经济稳定如果生产者对市场经济的蛛网模型有所了解，为了减少因价如果生产者对市场经济的蛛网模型有所了解，为了减少因价格波动而造成的经济损失，他应当提高自己的经营水平，不格波动而造成的经济损失，他应当提高自己的经营水平，不应当仅根据上一周期的价格来决定现阶段的生产量。例如可应当仅根据上一周期的价格来决定现阶段的生产量。例如可以根据本时段与前一时段价格的平均值来确定生产量。此时，以根据本时段与前一时段价格的平均值来确定生产量。此时，若若t 时段的商品量为时段的商品量为 xt 时，仍有时，仍有第27页，讲稿共45张，创作于星期日（7）将（将（5）式、（）

20、式、（7）式代入（）式代入（6）式，整理得）式，整理得（5）但但t+1时时段的商品量段的商品量则则不再不再为为而被修正而被修正为为（6）由（由（5）式得）式得（8）（8）式是一个常系数二）式是一个常系数二阶线阶线性差分方程，特征方程性差分方程，特征方程为为其特征根为其特征根为第28页，讲稿共45张，创作于星期日记记。若。若，则则此此时时差分方程（差分方程（8）是不）是不稳稳定的。定的。，若若，此，此时时特征根特征根为为一一对对共共轭轭复数，复数，。由线性差分方程稳定的条件，由线性差分方程稳定的条件，当当r2即即b2a时（时（9）式是稳定的，）式是稳定的，从从 而而M*是稳定的平衡点。是稳定的平

21、衡点。不难发现，生产者管理方式的不难发现，生产者管理方式的这一更动不仅使自己减少了因这一更动不仅使自己减少了因价格波动而带来的损失，而且价格波动而带来的损失，而且大大消除了市场的不稳定性。大大消除了市场的不稳定性。生产者在采取上述方式来确定生产者在采取上述方式来确定各时段的生产量后，如发现市各时段的生产量后，如发现市场仍不稳定（场仍不稳定（b2a），可按类），可按类似方法试图再改变确定生产量似方法试图再改变确定生产量的方式，此时可得到更高阶的的方式，此时可得到更高阶的差分方程。对这些方程稳定性差分方程。对这些方程稳定性条件的研究很可能会导出进一条件的研究很可能会导出进一步稳定市场经济的新措施。

22、步稳定市场经济的新措施。第29页，讲稿共45张，创作于星期日例例3 国民经济的稳定性国民经济的稳定性 国民收入的主要来源是生产，国民收入的开支主要用于消费资国民收入的主要来源是生产，国民收入的开支主要用于消费资金、投入再生产的积累资金及政府用于公共设施的开支。现在金、投入再生产的积累资金及政府用于公共设施的开支。现在我们用差分方程方法建立一个简略的模型，粗略地分析一下国我们用差分方程方法建立一个简略的模型，粗略地分析一下国民经济的稳定性问题。民经济的稳定性问题。再生再生产产的投的投资资水水平平It取决于消取决于消费费水平的水平的变变化量，化量，设设政府用于公共政府用于公共设设施的开支在一个不太

23、大的施的开支在一个不太大的时时期内期内变动变动不大，不大，设设为为常数常数G。故由。故由可得出可得出。将。将及及代入代入。记记yt为为第第t周期的国民收入，周期的国民收入，Ct为为第第t周期的消周期的消费资费资金。金。Ct的的值值决决定于前一周期的国民收入，定于前一周期的国民收入，设设第30页，讲稿共45张，创作于星期日（10）（9）式是一个二）式是一个二阶阶常系数差分方程，其特征方程常系数差分方程，其特征方程为为，相，相应应特征根特征根为为（11）成立成立时时才是才是稳稳定的。定的。（11）式可用于）式可用于预报经济发预报经济发展展趋势趋势。现现用待定系数法求方程用待定系数法求方程（10）的

24、一个特解）的一个特解。令。令代入（代入（10)式，得式，得故当（故当（11）式成立）式成立时时，差分方程，差分方程（10）的通解）的通解为为其中其中为为的模，的模，为为其幅角。其幅角。第31页，讲稿共45张，创作于星期日例如，若取例如，若取，易易见见，此，此时时关系式关系式（10）成立，）成立，又若又若取取y0=1600，y1=1700，G=550，则则由迭代公式由迭代公式求得求得 y2=1862.5,y3=2007.8,y4=2110.3,y5=2171.2,y6=2201.2,y7=2212.15,y8=2213.22,y9=2210.3,。易见易见第32页，讲稿共45张，创作于星期日例例

25、4 商品销售量预测商品销售量预测(实例实例)某商品前某商品前5年的销售量见表年的销售量见表。现希望根据。现希望根据 前前5年的统计年的统计数据预测数据预测 第第6年起该商品在各季度中的销售量。年起该商品在各季度中的销售量。从表中可以看出，该商品在从表中可以看出，该商品在 前前5年相同季节里的销售量呈增长趋势，年相同季节里的销售量呈增长趋势，而在同一年中销售量先增后减，第一季度的销售量最小而第三季度的而在同一年中销售量先增后减，第一季度的销售量最小而第三季度的销售量最大。预测该商品以后的销售情况，一种办法是应销售量最大。预测该商品以后的销售情况，一种办法是应 用用最小二最小二乘法乘法建立经验模型

26、。即根据本例中数据的特征，可以按季度建立四个经建立经验模型。即根据本例中数据的特征，可以按季度建立四个经验公式，分别用来预测以后各年同一季度的销售量。例如，如认为第一验公式，分别用来预测以后各年同一季度的销售量。例如，如认为第一季度的销售量大体按线性增长，可设销售量季度的销售量大体按线性增长，可设销售量 由由15253217152430151320271512182614111625121234第五年第五年第四年第四年第三年第三年第二年第二年第一年第一年销销售量售量季度季度 年份年份第33页，讲稿共45张，创作于星期日 求得求得 a=1.3，b=9.5。根据根据 预测第六年起第一季度的销售量预

27、测第六年起第一季度的销售量 为为 =17.3，=18.6，如认为销售量并非逐年等量增长而是如认为销售量并非逐年等量增长而是按前一年或前几年同期销按前一年或前几年同期销售量的一定比例增长的售量的一定比例增长的，则可建立相应的差分方程模型。仍以，则可建立相应的差分方程模型。仍以第一季度为例，为简便起见不再引入上标，以表示第一季度为例，为简便起见不再引入上标，以表示 第第t年第一年第一节季度的销售量，建立形式如下的差分方程：节季度的销售量，建立形式如下的差分方程：或或等等。等等。上述差分方程中的系数不一定能使所有上述差分方程中的系数不一定能使所有统计统计数据吻合，数据吻合，较为较为合合理的理的办办法

28、是用最小二乘法求一法是用最小二乘法求一组总组总体吻合体吻合较较好的数据。以建立好的数据。以建立二二阶阶差分方程差分方程为为例，例，为选为选取取a0,a1,a2使使最小，解最小，解线线性方程性方程组组：第34页，讲稿共45张，创作于星期日即求解即求解得得a0=-8，a1=-1，a2=3。即所求二。即所求二阶阶差分方程差分方程为为 第35页，讲稿共45张，创作于星期日虽然这一差分方程恰好使所有统计数据吻合，但这只是一个巧合。根虽然这一差分方程恰好使所有统计数据吻合，但这只是一个巧合。根据这一方程，可迭代求出以后各年第一季度销售量的预测值据这一方程，可迭代求出以后各年第一季度销售量的预测值 y6=2

29、1，y7=19，等。等。上述为预测各年第一季度销售量而建立的二阶差分方程，虽然其系上述为预测各年第一季度销售量而建立的二阶差分方程，虽然其系数与前数与前 5年第一季度的统计数据完全吻合，但用于预测时预测值与事年第一季度的统计数据完全吻合，但用于预测时预测值与事实不符。凭直觉，第六年估计值明显偏高，第七年销售量预测值甚至实不符。凭直觉，第六年估计值明显偏高，第七年销售量预测值甚至小于第六年。稍作分析，不难看出，如分别对每一季度建立一差分方小于第六年。稍作分析，不难看出，如分别对每一季度建立一差分方程，则根据统计数据拟合出的系数可能会相差甚大，但对同一种商品，程，则根据统计数据拟合出的系数可能会相

30、差甚大，但对同一种商品，这种这种 差异差异 应当是微小的，故应根据统计数据建立一个共用于各个季应当是微小的，故应根据统计数据建立一个共用于各个季度的差分方程。度的差分方程。为为此，将季度此，将季度编编号号为为t=1,2,20，令，令或或等，利用全体数等，利用全体数据来据来拟拟合，求合，求拟拟合得最好的系数。以二合得最好的系数。以二阶阶差分方程差分方程为为例，例，为为求求a0、a1、a2使得使得 最小最小第36页，讲稿共45张，创作于星期日求解求解线线性方程性方程组组即求解三元一次方程即求解三元一次方程组组解得解得a0=0.6937,a1=0.8737,a2=0.1941,故求得二阶差分方程故求

31、得二阶差分方程（t21）根据此式迭代，可求得第六年和第七年第一季度销售量的预测值为根据此式迭代，可求得第六年和第七年第一季度销售量的预测值为y21=17.58，y25=19.16还是较为可信的。还是较为可信的。第37页，讲稿共45张，创作于星期日例例5 人口问题的差分方程模型人口问题的差分方程模型 在在3.2中，我们已经讨论了人口问题的两个常微分方程模型中，我们已经讨论了人口问题的两个常微分方程模型Malthus模型模型和和Logistic模型）。前者可用于人口增长的短模型）。前者可用于人口增长的短期预测，后者在作中、长期预测时期预测，后者在作中、长期预测时 效果较好。效果较好。1、离散时间离

32、散时间 的的Logistic模型模型在研究人口或种群数量的实际增长情况时，有时采用离散化在研究人口或种群数量的实际增长情况时，有时采用离散化的时间变量更为方便。例如，有些种群具有相对较为固定的的时间变量更为方便。例如，有些种群具有相对较为固定的繁殖期，按时段统计种群数量更接近种群的实际增长方式。繁殖期，按时段统计种群数量更接近种群的实际增长方式。人口增长虽无这种特征，但人口普查不可能连续统计，任何人口增长虽无这种特征，但人口普查不可能连续统计，任何方式的普查都只能得到一些离散时刻的人口总量（指较大范方式的普查都只能得到一些离散时刻的人口总量（指较大范围的普查）。这样，如何建立人口问题的离散模型

33、的问题十围的普查）。这样，如何建立人口问题的离散模型的问题十分自然地提了出来。分自然地提了出来。第38页，讲稿共45张，创作于星期日建立离散模型的一条直接途径是建立离散模型的一条直接途径是用用差分代替微分差分代替微分。从人口。从人口问问题题的的Logistic模型模型可可导导出一出一阶阶差分方程差分方程（12）（12）式中右端的因子）式中右端的因子 常被称为常被称为阻尼因子阻尼因子。当当PtN时，种群增长接时，种群增长接 近近Malthus模型；当模型；当Pt接近接近N时，这一因时，这一因子将越来越明显地发挥阻尼作用，子将越来越明显地发挥阻尼作用，若若PtN，它将使种群增，它将使种群增长速度长

34、速度 在在Pt接近接近N时变得越来越慢，若时变得越来越慢，若 PN，它将使种群呈，它将使种群呈负增长。负增长。（12）式可改写）式可改写为为（13）记记,于是于是(13)式又可改写式又可改写为为（14）第39页，讲稿共45张，创作于星期日虽然，（虽然，（14）式是一个非线性差分方程，但对确定的初值）式是一个非线性差分方程，但对确定的初值x0，其后的其后的 x1可利用方程确定的递推关系迭代求出。可利用方程确定的递推关系迭代求出。差分方程（差分方程（14）有两个平衡点，）有两个平衡点，即即x*=0和和 。类似。类似于微分方程稳定性的讨论，非线性差分方程平衡点的稳定性也于微分方程稳定性的讨论，非线性

35、差分方程平衡点的稳定性也可通过对其线性近似方程平衡点稳定性的讨论部分地得到确定可通过对其线性近似方程平衡点稳定性的讨论部分地得到确定（时不能确定除外）。例如，对时不能确定除外）。例如，对 ，讨论讨论 在在x*处的线性近似方程处的线性近似方程可知，当可知，当（即（即）时时平衡点平衡点是是稳稳定的，此定的，此时时（）若当若当 ，则平稳点，则平稳点 是不稳定的，（这与对是不稳定的，（这与对 一切一切a，p*=N均为均为Logistic方程的稳定平衡点不同）。方程的稳定平衡点不同）。第40页，讲稿共45张，创作于星期日 作业第41页，讲稿共45张，创作于星期日设差分方程求第42页，讲稿共45张，创作于星期日n参参 考考 文文 献献n1.数学模型数学模型 姜启源（第二版）姜启源（第二版）P271-303 高等教育出版社高等教育出版社n2.数学建模方法数学建模方法 刘承平刘承平 P95-105 高等教高等教育出版社育出版社n3.数学模型引论数学模型引论 唐焕文唐焕文、贺明峰（第二版）、贺明峰（第二版）P239-247 高等教育出版社高等教育出版社第43页，讲稿共45张，创作于星期日例设差分方程求解法1：特征方程为有根：故：为方程的解.由条件得：第44页，讲稿共45张，创作于星期日感感谢谢大大家家观观看看第45页，讲稿共45张，创作于星期日