讲差分方程模型.pptVIP

数学建模之动态模型 讲授主要内容 一、差分方程知识简介 二、差分方程模型讲解 三、微分方程知识简介 四、微分方程模型讲解 一、 差分方程知识简介 二、差分方程模型讲解 差分方程建模方法的思想 差分方程建模方法的思想与一般数学建模的思想是一致的，也需要经历背景分析、确定目标、预想结果、引入必要的数值表示（变量、常量、函数、积分、导数、差分、取最等）概念和记号、几何形式（事物形状、过程轨迹、坐标系统等）也就是说要把事物的性态、结构、过程、成份等用数学概念、原理、方法来表现、分析、求解。 差分方程的特殊性： 首先应该把系统或过程进行特别分解，形成表现整个系统的各个部分的离散取值形式或形成变化运动过程的时间或距离的分化而得到离散变量。 然后通过内在的机理分析，找出变量所能满足的平衡关系、增量或减量关系及规律，从而建立差分方程。 特别地，有时可能通过多个离散变量的关系得到我们关心的变量的关系，这实际上建立的是离散向量方程，它有着非常重要的意义。有时需要找出决定变量的初始条件；有时还需要将问题适当分成几个子部分并分别求解。 * * 动态 模型 特征 描述对象特征随时间(空间)的演变过程 分析对象特征的变化规律 预报对象特征的未来性态 研究控制对象特征的手段 在研究实际问题时, 我们常常不能直接得出变量之间的关系,但却能容易得出包含变量导数在内的关系式,这就是微分方程. 在现实社会中,又有许多变量是离散变化的,如人口数、生产周期与商品价格等, 而且离散的运算具有可操作性, 差分正是联系连续与离散变量的一座桥梁. 不管是微分方程还是差分方程模型，有时无法得到其解析解(必要时,可以利用计算机求其数值解),既使得到其解析解,尚有未知参数需要估计(这时可利用参数估计方法). 而在实际问题中,讨论问题的解的变化趋势很重要，因此，以下只对其平衡点的稳定性加以讨论. 对于k阶差分方程 F( n; xn, xn+1, … , xn+k ) = 0 (1) 若有xn = x (n), 满足 F(n; x(n), x(n + 1) , … , x(n + k )) = 0, 则称xn = x (n)是差分方程(1)的解, 包含任意常数的解称为(1)的通解, x0, x1, … , xk-1为已知时称为(1)的初始条件,通解中的任意常数都由初始条件确定后的解称为(1)的特解. 若x0, x1, … , xk-1已知, 则形如 xn+k = g(n; xn, xn+1, … , xn+k-1 ) 的差分方程的解可以在计算机上实现. 若有常数a是差分方程(1)的解, 即 F (n; a, a, … , a ) = 0, 则称 a是差分方程(1)的平衡点. 又对差分方程(1)的任意由初始条件确定的解 xn= x(n)都有 xn→a (n→∞), 则称这个平衡点a是稳定的. 一阶常系数线性差分方程 xn+1 + axn= b, (其中a, b为常数, 且a ≠-1, 0)的通解为 xn=C(- a) n + b/(a + 1) 易知b/(a+1)是其平衡点, 由上式知, 当且仅当|a|＜1时, b/(a +1)是稳定的平衡点. 二阶常系数线性差分方程 xn+2 + axn+1 + bxn = r, 其中a, b, r为常数. 当r = 0时, 它有一特解 x* = 0； 当r ≠ 0, 且a + b + 1≠ 0时, 它有一特解 x*=r/( a + b +1). 不管是哪种情形, x*是其平衡点. 设其特征方程 ?2 + a? + b = 0 的两个根分别为 ?=?1, ?=?2. ① 当?1, ?2是两个不同实根时,二阶常系数线性差分方程的通解为 xn= x*+ C1(?1)n + C2(?2)n ; ② 当?1, 2=?是两个相同实根时,二阶常系数线性差分方程的通解为 xn= x* + (C1 + C2 n)?n; ③ 当?1, 2= ? (cos?( +)(-) i sin? ) 是一对共轭复根时,二阶常系数线性差分方程的通解为 xn = x*+ ? n (C1cosn? + C2sinn? ). 易知,当且仅当特征方程的任一特征根 |?i |＜1时, 平衡点x*是稳定的. 则 对于一阶非线性差分方程 xn+1 = f (xn ) 其平衡点x*由代数方程 x = f (x) 解出. 为分析平衡点x*的稳定性, 将上述差分方程近似为一阶常系数线性差分方程 时,上