差分方程模型ppt课件.ppt

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

1、 差分方程模型差分方程模型第四讲第四讲 养老保险养老保险 第一讲第一讲 差分方程差分方程第二讲第二讲 蛛网模型蛛网模型第三讲第三讲 商品销售量预测商品销售量预测1、差分方程简介、差分方程简介t 规定 只取非负整数，记 为变量在 点的取值，则称 为 的一阶向前差分，称为 的二阶差分。 由 、 及 的差分给出的方程称为 的差分方程。其中含 的最高阶差分的阶数称为该差分方程的阶数。差分方程也可以写成不显含差分的形式，例如二阶差分方程 可以写成ttytttyyy1tytttttttyyyyyyy12122)(tyttttytytyty02tttyyy012tttyyy 满足一阶差分方程的序列 称为差分

2、方程的解，若解中含有独立的常数的个数等于差分方程的阶数时，称此解为该差分方程的通解。 称如下形式的差分方程 为 阶常系数线性差分方程，其中 是常数， 。其对应的齐次方程为ty)(110tbyayayatnntntnnaaa,1000a0110tnntntyayaya 求非齐次常系数线性差分方程的通解的步鄹： 1.先求解对应的特征方程 2.根据特征根的不同情况，求解齐次方程的通解 若特征方程有 个不同的实根 ，则齐次方程的通解为 ； 若 是特征方程的 重实根，则齐次方程的通解为 ； 若特 征方程有单重复根 ，则齐次方程的通解为 ，其中 为 的模， 为 的幅角； 0110nnnaaann,1tnn

3、tcc11ktkktcc)(11itctcttsincos1122arctan 若特征方程有 重复根 ，则齐次方程的通解为 3.求非齐次方程的一个特解 ，若 为齐次方程的通解，则非齐次方程的通解为 。 对特殊形式的特解 可以使用待定系数法求非齐次方程的特解。例如 ， 为 的 次多项式时可以证明：若 不是特征根，则非齐次方程有形如 的特解， 也是 的 次多项式；若 是 重特征根，则非齐次方程有形如 的特解。进而可以用待定系数法求出 ，从而得到非齐次方程的一个特解。kittccttcctkktkksin)(cos)(1111tytyttyy )(tb)()(tpbtbkt)(tpktkb)(tqb

4、kt)(tqktkbr)(1tqtbkrt)(tqk例例1. 1. 求解两阶差分方程求解两阶差分方程 解 对应齐次方程的特征方程为 ，其特征根为 ，故齐次方程的通解为，原方程有形如 的特解，带入原方程求得 ，所以原方程的通解为在应用差分方程研究问题时，需要讨论解的稳定性。对常系数非齐次线性差分方程，若不论其对应齐次方程的通解中的任意常数如何取值，当 时， ，则称方程的解是稳定的。tyytt2012i2, 1tctcyt2sin2cos21bat 2/1, 2/1ba21212sin2cos21ttctct0ty2、常系数线性差分方程的、常系数线性差分方程的 变换解法变换解法 采用上述解析解法求

5、解常系数线性非齐次差分方程比较繁琐，下面介绍 变换，将差分方程转化为代数方程去求解 设有离散序列 ，则 的 变换定义为其中 是复变量，上式右端的级数的收敛域是某个圆的外部 的 反变换记作 ZZ, 1 , 0),(kkx)(kxZ0)()()(kkzkxkxZzXz)(zXZ)()(1zXZkx几个常用离散函数的几个常用离散函数的 变换变换 (3) 单边指数函数 的 变换( 为不等于1的正常数)(1) 单位冲激函数 的 变换(2) 单位阶跃函数 的 变换Z0011 )()(kkkkzzkkZ)(kZ)(kUZ) 1|(|11)()(00zzzzzkUkUZkkkkkakf)(Za)|(|0aza

6、zzzaaZkkkk(1)线性性质 设 ，则 变换的性质变换的性质(2)平移性质：设 ，则Z)()(),()(2211zXkxZzXkxZ)()()()(2121zbXzaXkbxkaxZ)()(zXkxZ)0()()1(xzXzkxZ)()()(10NkkNzkxzXzNkxZ) 1()()1(1zxzXzkxZ)()()(1NkkNzkxzXzNkxZ例例2. 2. 求解齐次差分方程求解齐次差分方程 解 令 ，对差分方程取 变换得 对上式取 反变换，便得差分方程的解为1) 1 (, 0)0(, 0)(2) 1(3)2(xxkxkxkx)()(zXkxZZ0)(2)(3)(2zXzzXzzX

7、z2123)(2zzzzzzzzXZkkkx)2() 1()(1、问题的提出、问题的提出 在自由竞争的社会中，很多领域会出现循环波动的现象。在经济领域中，可以从自由集市上某种商品的价格的变化看到如下现象：在某一时期，商品的上市量大于需求，引起价格下跌，生产者觉得该商品无利可图，转而经营其他商品，一段时间之后，随着产量的下降，供不应求又会导致价格上升，又会有很多生产商进行该商品的生产，随之而来的是商品过剩，价格下降。在没有外界干预的情况下，这种现象会反复出现。 如何从数学的角度来描述上述现象呢？2、模型假设、模型假设 (1)设 时段商品数量为 ，其价格为 ，这里把时间离散化为时段，一个时期相当于

8、商品的一个生产周期。 (2)同一时段的商品价格取决于该时段商品的数量，称为需求函数。出于对自由经济的理解，商品的数量越多，其价格就越低。故可以假设需求函数为一个单调递减函数。 (3) 下一时段的商品数量由上一时段的商品价格决定，称为供应函数，由于价格越高可导致产量越大，所以可以假设供应函数是一个单调递增的函数。kkxky)(kkxfy )(1kkygx3、模型求解、模型求解 在同一坐标系中同时做出供应函数和需求函数的图形，设两条曲线相交于 则 为平衡点。因为此时 若某个 有 ，则可推出即商品的数量保持在 ，价格保持在 。不妨假设 下面考虑 在图上的变化如右图所示。),(000yxP0P)(),

9、(0000 xfyygxk0 xxk), 1( ,00kklxxyyll0 x0y01xx kkyx ,当 给定后，价格 由 上的 点决定，下一时段的数量 由 上的 点决定， 又可由 上的点 决定。依此类推，可得一系列的点图上的箭头表示求出 的次序，由图知即市场经济趋于稳定。1x1yf1P1P2xg2P2yff3P),(),(),(),(234223122111yxPyxPyxPyxPkPkP),(),(lim000yxPyxPkk并不是所有的需求函数和供应函数都趋于稳定，若给定的 和 的图形如右图所示，得到的就不趋于 ，此时市场经济趋于不稳定。gfg,21PP0P图1和图2中的折线 形如蛛网

10、，故把这种模型称为蛛网模型。在进行市场经济分析中， 取决于消费者对某种商品的需求程度及其消费水平， 取决于生产者的生产、管理等能力。 当已知需求函数和供应函数之后，可以根据和 的性质判断平衡点 的稳定性。当较小时， 的稳定性取决于 和 在点 的斜率，即当时， 点稳定。当时， 点不稳定。,3221PPPPfgfg0P|01xx 0Pfg0P| )(| )(|00ygxf0P| )(| )(|00ygxf0P 这一结论的直观解释是，需求曲线越平，供应曲线越陡，越有利于经济稳定。 设 在 点附近取 和 的线性近似，于是得从上两式中消去 得 .| )(|1|,)(|00ygxf0Pfg)(00 xxy

11、ykk)(001yyxxkkky02231201201)1 ()()()()1)()()()1 (xxxxxxxxxkkkkkk以上 个式子相加，有若 是稳定点，则应有 ，所以 点稳定的条件是同理 点不稳定的条件是 01121022132)1 ()()()()1 ()()()(xxxxxxkkkkkkk011011)(1 )()()(1 )1 ()(xxxxxkkkkk0P0limxxkk0P10P14、模型修正、模型修正 在上述模型的基础上，对供应函数进行改进。下面在决定商品的生产数量 时，不仅考虑前一时期的价格 ，而且考虑了价格 ，取 ，在 附近取线性近似，则有于是得将上述两式整理得到二阶

12、线性差分方程 1kxky1ky)2(11kkkyygx0P)2(20101yyyxxkkk)(00 xxyykk)(0101xxyykk), 3 , 2( ,)1 (22011kxxxxkkk其特征方程为经计算得其特征根结论：若方程的特征根均在单位圆内，则 为稳定点。当 时，该特征方程有两个实根，因则有 ，故此时 不是稳定点。当 时，特征方程有两个共轭复根，共轭复根的模的绝对值为02248)(22, 10P8448)(222|20P824)(84)(|2/1222, 1要使 点为稳定点只需 与前面的模型的结果相比， 的范围扩大了。这是由于经营管理者的水平提高带来的结果。0P2,商品销售量预测商

13、品销售量预测 在利用差分方程建模究实际问题时，常常需要根据统计数据并用最小二乘法来拟合出差分方程的系数。其系统稳定性的讨论要用到代数方程的求根。 例3 某商品前五年的销售量见右表1。现希望根据前五年的统计数据预测第六年起该商品在各季度中的销售量。第一第一 年年第二第二 年年第三第三 年年第四第四 年年第五第五 年年1111213151621618224253252627303241214151517 由于该问题的数据少，用回归分析效果不一定好。 如果认为销售量并非逐年等量增长而是按前一年或者前几年同期销售量的一定比例增长的，则可建立相应的差分方程模型。以第一季度为例，以 表示第 年第一季度的销

14、售量，建立形式如下的差分方程上述差分方程不一定能使所有统计数据吻合，较为合理的办法是用最小二乘法求一组总体吻合较好的数据。即选取 使得最小。根据这一方程可以迭代求解以后各年第一 季度销售tyt32211ayayayttt321,aaa53232211ttttayayay量的预测值 。第7年销售量预测值居然小于第6年的，稍作分析，不难看出，如分别对第一季度建立差分方程，则根据统计数据拟合出的系数可能会相差甚大，但对同一种商品，这种差异应当是微小的，故应根据统计数据建立一个共用于各个季度的差分方程，为此，将季度编号为 ，令 ，利用全体数据来拟合求拟合得到最好的系数。即求 使得最小。于是得二阶差分方

15、程为由此式可得 ，这个结果还是较为可信的。,19,2176yy20, 2 , 1t38241ayayayttt321,aaa209238241321),(ttttayayayaaaQ)21( ,6957. 01941. 08737. 084tyyyttt1676.19,5869.172521yy1、问题的提出、问题的提出 某保险公司的一份材料指出，在每月交费200元至59岁年底，60岁开始领取养老保险金的约定下，男子若25岁开始投保，届时月养老金2282元；假定人的寿命为75岁，是求出保险公司为了兑现保险责任，每月应至少有多少投资收益率 （也就是投保人的实际收益率）？ 2、模型的建立与求解、模

16、型的建立与求解 设投保人在投保后第 个月所交保险费及利息的累计总额为 ，那么得到数学模型为分段表示的差分方程其中 分别为60岁前所交月保险费和60岁起所领取的月养老金的数目（月）， 是所交保险金获得的利率， 分别是自投保起至停交保险费和至停领取养老金的时间（月），这里 ，可推出差分方程的解（这里 ），得 kkF1, 1,)1 (1, 1 , 0,)1 (11MNNkqrFFNkprFFkkkkqp,rMN,600,420,2282,200MNqp00MFF于是得到由于 ，可得如下方程MNkrqrFNkrprFMkkkk, 1,)1 (1 , 1 , 0, 1)1(rprFN1)1(rqrFMNN)1 (1 11qrFFNN)1 (1qrprrqrNMN 1)1()1 (1 10)1)(1 ()1 (pqrpqrMNMrx1041.1141.12180600 xx%49. 0,0049. 1rx