算法大全差方程模型.doc

算法大全差分方程模型 第十六章差分方程模型 离散状态转移模型涉及的范围很广，可以用到各种不同的数学工具。下面我们对差 分方程作一简单的介绍，下一章我们将介绍马氏链模型。 §1 差分方程 1.1差分方程简介 规定t只取非负整数。记 yt为变量 y在t点的取值，则称 Δyt = yt+1. yt 为 yt的一 阶向前差分，简称差分，称 Δ2 yt =Δ(Δyt ) =Δyt+1 .Δyt = yt+2 . 2 yt+1 + yt 为 yt的二 阶差分。类似地，可以定义 yt的n阶差分Δnyt 。 由t、yt 及 yt的差分给出的方程称为 yt的差分方程，其中含 yt的最高阶差分的阶 数称为该差分方程的阶。差分方程也可以写成不显含差分的形式。例如，二阶差分方程 Δ2 yt +Δyt + yt = 0 也可改写成 yt+2 . yt+1 + yt = 0 。 满足一差分方程的序列 yt称为差分方程的解。类似于微分方程情况，若解中含有 的独立常数的个数等于差分方程的阶数时，称此解为该差分方程的通解。若解中不含任 意常数，则称此解为满足某些初值条件的特解。 称如下形式的差分方程 a0 yn+t + a1 yn+t.1 +L+ an yt = b(t) （1） 为n阶常系数线性差分方程，其中 a0,a1,L,an 是常数，a0 ≠ 0 。其对应的齐次方程为 a0 yn+t + a1 yn+t.1 +L+ an yt = 0 （2） 容易证明，若序列 yt(1) 与 yt(2) 均为（2）的解，则 yt =c1 yt (1) + c2 yt (2) 也是方程（2）的 解，其中 c1,c2 为任意常数。若 yt(1) 是方程（2）的解， yt(2) 是方程（1）的解，则 yt =yt (1) + yt (2) 也是方程（1）的解。 方程（1）可用如下的代数方法求其通解： （I）先求解对应的特征方程 a0λn + a1λn.1 +L+ a0 = 0 （3） （II）根据特征根的不同情况，求齐次方程（ 2）的通解。 （i）若特征方程（ 3）有n个互不相同的实根 λ1,L,λn，则齐次方程（ 2）的通解 为 c1λ1 t +L+ cnλtn （c1,L,cn 为任意常数） （ii）若λ是特征方程（ 3）的k重根，通解中对应于 λ的项为(c1 +L+ cktk .1)λt ， ci (i =1,L,k) 为任意常数。 （iii）若特征方程（ 3）有单重复根 λ =α±βi ，通解中对应它们的项为 c1 ρ t cos.t +c2 ρ t sin.t ，其中ρ=α 2 +β 2 为λ的模，.= arctgβ 为λ的幅角。 α （iv）若λ =α±βi 是特征方程（3）的k重复根，则通解对应于它们的项为 k .1 tk .1 t (c +L+ ct )ρ cos.t + (c +L+ ct )ρ sin.t 1 kk +12k -192-  ci (i =1,L,2k) 为任意常数。 （III）求非齐次方程（ 1）的一个特解 yt。若 yt为方程（2）的通解，则非齐次方 程（1）的通解为 yt +yt 。 求非齐次方程（ 1）的特解一般要用到常数变易法，计算较繁。对特殊形式的 b(t) 也可使用待定系数法。例如，当b(t) =bt pk (t) ， pk(t) 为t的k次多项式时可以证明： 若b不是特征根，则非齐次方程（ 1）有形如 btqk (t) 的特解， qk(t) 也是t的k次多项 式；若 b是r重特征根，则方程（ 1）有形如 bttrqk (t) 的特解。进而可利用待定系数法 求出qk(t) ，从而得到方程（1）的一个特解 yt。 例 1 求解两阶差分方程 yt+2+ yt = t 。 解对应齐次方程的特征方程为 λ2 +1 = 0 ，其特征根为 λ1,2=±i ，对应齐次方程 的通解为 ππ =c cos t + c sin t yt 12 22 原方程有形如 at +b 的特解。代入原方程求得 a =1 ，b =.1 ，故原方程的通解 22 为 ππ 11 c cos t + c sin t + t . 1 2 2222 例 2 在信道上传输仅用三个字母 a,b,c 且长度为n的词，规定有两个a连续出现 的词不能传输，试确定这个信道容许传输的词的个数。 解令h(n)表示容许传输且长度为 n的词的个数， n =1,2,L ，通过简单计算可 求得：h(1) = 3