差分方程模型预案.pptVIP

* 差分方程模型 一、差分方程简介 以t 表示时间，规 定t只取非负整数。t=0表示第一周期初，t=1表示第二周期初等。 记yt 为变量y在时刻t 时的取值，则称 为yt 的一阶差分，称 为的二阶差分。类似地，可以定义yt的n阶差分。 由t、yt及yt的差分给出的方程称 为yt差分方程，其中含的最高阶差分的阶数称为该差分方程的阶。差分方程也可以写成不显含差分的形式。例如，二阶差分方程 也可改写成 满足一差分方程的序 列yt称为此差分方程的解。类似于微分方程情况，若解中含有的独立常数的个数等于差分方程的阶数时，称此解为该差分方程 的通解。若解中不含任意常数，则称此解为满足某些初值条件的 特解，例如，考察两阶差分方程 易见 与 均是它的特解，而 则为它的通解，其 中c1，c2为两个任 意常数。 类似于微分方程，称差分方程 为n阶线性差分方程， 当 ≠0时称其为n阶非齐次线性差分方程，而 则被称为方程对应的 齐次线性差分方程 。 若所有的 ai(t)均为与t无关的常数，则称其为 常系数差分方程，即n阶常系数线性差分方程可分成 （1） 的形式，其对应的齐次方程为 （2） 容易证明，若序列 与 均为方程（2）的解，则 也是方程（2）的解，其 中c1、c2为任意常数，这说明， 齐次方程的解构成一个 线性空间（解空间）。 方程（1）可用如下的代数方法求其通解： （步一）先求解对应的特征方程 （3） （步二）根据特征根的不同情况，求齐次方程(2)的通解 情况1 若特征方程（3）有n个互不相同的实根 ,…, ，则齐次方程（2）的通解为 (C1,…,Cn为任意常数) ， 情况2 若λ?是特征方程（3）的k重根，通解中对应 于λ的项为 为任意常数，i=1,…,k。 情况3 若特征方程（3）有单重复根 通解中对应它们的项为 为λ的模， 为λ的幅角。 情况4 若 为特征方程（3）的k重复根，则通 解对应于它们的项为 为任意常数，i=1,…,2k。 . 若yt为方程(2)的通解,则非齐次方程 (1)的通解为 （步三） 求非齐次方程 (1)的一个特解 求非齐次方程（1）的特解一般要用到 常数变易法，计算较繁。对特殊形式 的b(t)也可使用 待定系数法。 例1 求解两阶差分方程 解 对应齐次方程的特征方程为 ，其特征根为 ，对应齐次方程的通解为 原方程有形如 的特解。代入原方程求得 ， ，故原方程的通解为 在应用差分方程研究问题时，一般不需要求出方程的通解，在给定初值以后，通常可用 计算机迭代求解，但我们常常需要讨论解的稳定性。对差分方程 (1) ,若不论其对应齐次方程的通解中任意常数 C1,…,Cn 如何取值 , 在 时总 有 ,则称方程 (3)的解是稳定的,否则称其解为不稳定 的.根据通解的结构不难看出 ,非齐次方程(1)稳定的充要条件为其所有特征根的模均小于1。 三、差分方程的平衡点及稳定性 1．? 一阶线性常系数差分方程的平衡点及稳定性 一阶线性常系数差分方程 xk+1+axk=b,k=0,1,2,…（1） 的平衡点由x+ax=b解得，为 ，当 时，若xk?x*，则x*是稳定的。 方程（1）的平衡点的稳定性问题可以通过变量代换转换为齐次方程 xk+1+axk=0,k=0,1,2… (2) 的平衡点x*=0的稳定性问题。而对于方程（2），其解可以表示为 xk=(-a)kx0, k=1,2,… (3) 所以当且仅当|a|1时，方程（2）（从而方程（1））的平衡点是稳定的。 对于n维向量x(k)和n*n常数矩阵A构成的方程组 x(k+1)+Ax(k)=0 其平衡点稳定的条件是A的特征根λi,I=1,2,…，均有|λi|1。 2．? 二阶线性差分方程的平衡点及稳定性 考察二阶线性差分方程xk+a1xk+1+a2xk+2=0 (4) 在平衡点x*=0的稳定性。为求（4）