下载文档原格式
差分方程对连续型变量而言，我们常常回到微分方程的问题。
对离散型变量将导致一类的问题。
一、差分的定义定义：设()t y y t =是一个函数，自变量从t 变化到1t +，这时函数的增量记为(1)()t y y t y t ∇=+-，我们称这个量为()y t 在点t 步长为1的一阶差分，简称为()y t 的一阶差分。
为了方便我们也记1(1),()t t y y t y y t +=+=，即1t t t y y y +∇=-。
称21121()()()2t t t t t t t t y y y y y y y y +++++∇∇=---=-+为()y t 的二阶差分，简记为2t y ∇。
同样记2()t y ∇∇为3t y ∇，并称为三阶差分。
一般记1()n n t t y y -∇=∇∇，称为n 阶差分，且有0(1)nni it n t n i i y Cy +-=∇=-∑。
性质：当,,a b C 是常数，t y 和t z 是函数时， （1）()0C ∇=； （2）()()t t C y C y ∇=∇；（3）()()()t t t t ay bz a y b z ∇±=∇±∇；（4）11()()()()()t t t t t t t t t t y z y z z y y z z y ++∇⋅=∇+∇=∇+∇；（5）1111()()()()t t t t t t t t t t t t t t y z y y z z y y z z z z z z ++++⎛⎫∇-∇∇-∇∇== ⎪⋅⋅⎝⎭，（其中，0t z ≠）。
例1：已知,(0)nt y t t =≠，求()t y ∇。
解：()(1)n n t y t t ∇=+-。
特别，当n 为正整数时，1()ni n i t ni y Ct-=∇=∑，阶数降了一阶。
推论：若,m n 为正整数且m n >时，()P t 为n 次多项式，则()0m P t ∇=。
差分方程百科内容来自于:差分方程是含有未知函数及其导数的方程，满足该方程的函数称为差分方程的解。
基本概念一、差分的概念设函数yt=f(t)在t=…，-2，-1，0，1，2，…处有定义,对应的函数值为…，y-2，y-1，y0，y1，y2，…，则函数yt=f(t)在时间t的一阶差分定义为Dyt=yt+1-yt=f(t+1)-f(t)。
依此定义类推，有Dyt+1=yt+2-yt+1=f(t+2)-f(t+1)，Dyt+2=yt+3-yt+2=f(t+3)-f(t+2)，………………一阶差分的性质(1) 若yt=C(C为常数)，则Dyt=0；(2) 对于任意常数k，D(kyt)=kDyt；(3) D(yt+zt)=Dyt+Dzt。
函数yt=f(t)在时刻t的二阶差分定义为一阶差分的差分，即D2yt= D (D yt)= D yt+1- D yt=(yt+2-yt+1)-(yt+1-yt)=yt+2-2yt+1+yt．依此定义类推，有D2yt+1= Dyt+2- Dyt+1=yt+3-2yt+2+yt+1，D2yt+2= Dyt+3-Dyt+2=yt+4-2yt+3+yt+2，………………类推，计算两个相继的二阶差分之差，便得到三阶差分D3yt= D2yt+1- D2yt=yt+3-3yt+2+3yt+1-yt，D3yt+1= D2yt+2- D2yt+1=yt+4-3yt+3+3yt+2-yt+1，………………一般地,k阶差分(k为正整数)定义为这里二、差分方程含有未知函数yt=f(t)以及yt的差分Dyt，D2yt，…的函数方程，称为常差分方程(简称差分方程)；出现在差分方程中的差分的最高阶数，称为差分方程的阶。
n阶差分方程的一般形式为F(t，yt，Dyt，…，Dnyt)=0，其中F是t，yt, Dyt，…，Dnyt的已知函数，且Dnyt一定要在方程中出现。
含有两个或两个以上函数值yt，yt+1，…的函数方程，称为(常)差分方程，出现在差分方程中未知函数下标的最大差，称为差分方程的阶。
差分方程是含有未知函数及其导数的方程，满足该方程的函数称为差分方程的解。
数学意义及性质意义差分方程是微分方程的离散化。
一个微分方程不一定可以解出精确的解，把它变成差分方程，就可以求出近似的解来。
比如dy+y*dx=0,y(0)=1 是一个微分方程，x取值[0,1] (注：解为y(x)=e^(-x)); 要实现微分方程的离散化，可以把x的区间分割为许多小区间[0,1/n],[1/n,2/n],...[(n-1)/n,1] 这样上述微分方程可以离散化为：差分方程y((k+1)/n)-y(k/n)+y(k/n)*(1/n)=0，k=0,1,2,...,n-1 （n 个离散方程组）利用y(0)=1的条件，以及上面的差分方程，就可以计算出y(k/n) 的近似值了。
§1 基本理论差分方程1. 差分2. 任意数列{xn },定义差分算子Δ如下：Δxn=xn+1-xn 对新数列再应用差分算子，有Δ2xn=Δ（Δkxn).性质性质1 Δk(xn+yn)=Δkxn+Δkyn性质2 Δk(cxn)=cΔkxn性质3 Δkxn=∑（-1)jCjkXn+k-j性质4 数列的通项为n的无限次可导函数，对任意k>=1，存在η，有Δkxn=f(k）（η）差分方程定义8。
1 方程关于数列的k阶差分方程：xn-a1xn-1-a2xn-2-……aBxn-k=b(n=k,k+1，……）其中a1,a2,------ak 为常数，ak≠0. 若b=0，则该方程是齐次方程关于λ的代数方程λk-a1λk-1-------ak-1λ-ak=0 为对应的特征方程，根为特征值。
例题1．实验内容与练习2．1 插分例1 Xn={n3},求各阶差分数列：xn△xn △2xn △3xn △4xn 1 7 12 6 0 8 19 18 6 0 27 37 24 6 0 64 61 30 6 125 91 36 216 127 343 可见，{n3},三阶差分数列为常数数列，四阶为0。
什么叫差分方程？给我举几个例子呗§1 基本理论1. 差分2. 任意数列{xn },定义差分算子Δ如下：Δxn=xn+1-xn对新数列再应用差分算子，有Δ2xn=Δ(Δkxn).性质性质1 Δk(xn+yn)=Δkxn+Δkyn性质2 Δk(cxn)=cΔkxn性质3 Δkxn=∑(-1)jCjkXn+k-j性质4 数列的通项为n的无限次可导函数，对任意k>=1,存在η，有Δkxn=f(k)(η) 差分方程定义8。
1 方程关于数列的k阶差分方程：xn-a1xn-1-a2xn-2-……aBxn-k=b(n=k,k+1,……)其中a1,a2,------ak 为常数，ak≠0. 若b=0,则该方程是齐次方程关于λ 的代数方程λk-a1λk-1-------ak-1λ-ak=0为对应的特征方程，根为特征值。
1．实验内容与练习2．1 插分例1 Xn={n3},求各阶差分数列：xn △xn △2xn △3xn △4xn1 7 12 6 08 19 18 6 027 37 24 6 064 61 30 6125 91 36216 127343可见，{n3},三阶差分数列为常数数列，四阶为0。
练习1 对{1}，{n},{n2},{n4},{n5}, 分别求各阶差分数列。
练习2 {C0n-1}{C1n-1}{C2n-1},{C4n-1},分别求各阶差分数列.{Xn}的通项为n的三次函数，Xn=a3n3+a2n2+a1n+a0证明它为常数数列。
证明由Xn=a3n3+a2n2+a1n+a0可直接计算。
定理8。
1 若数列的通项是关于n 的k次多项式，则k 阶差分数列为非零数列，k+1阶差分数列为0。
练习3 证明定理8。
1 。
定理8。
2 若{Xn}的k 阶插分为非零常数列，则{Xn}是n的k次多项式，练习4 根据插分的性质证明定理8。
2例2。
求∑i3例3例4解设Sn=∑i3 表Sn △Sn △2Sn △3Sn △4Sn △5Sn1 8 19 18 6 09 27 37 24 6 036 64 61 30 6 0100 125 91 36 6 0225 216 127 42441 343 169784 5121296设Sn=a4n4+a3n3+a2n2+a1n+a0, s1=1,s2=9,s3=36,s4=100,s5=225,得a0=0, a1=0, a2=1/4, a3=1/2, a4=1/4.所以，Sn=(1/4)n4+(1/2)n3+(1/4)n2.练习{Xn}的通项Xn为n的k次多项式，证明∑xi为n的k+1次多项式；求∑i4.由练习2 {Crn-1}可得。
2.2差分方程对于一个差分方程，如果能找出这样的数列通项，将它带入差分方程后，该方程成为恒等式，这个通项叫做差分方程的解。
例3 对差分方程xn-5xn-1+6xn-2=0,可直接验证xn=c13n+c22n是该方程的解。
例3中的解中含有任意常数，且任意常数的个数与差分方程的阶数相同。
这样的解叫做差分方程的通解。
若k阶差分方程给定了数列前k项的取值，则可以确定通解的任意常数，得到差分的特解。
例4对差分方程xn-5xn-1+6xn-2=0,若已知x1=1,x2=5,则可以得到该差分方程的特解为xn=3n-2n.我们首先研究齐次线性差分方程的求解。
xn=rxn-1对一阶差分方程x1=a显然有xn=arn-1。
因此，若数列满足一阶差分方程，则该数列为一个等比数列。
例5 求Fibonacci数列{Fn}的通项，其中F1=1,F2=1,Fn=Fn-1+Fn-2.Fibonacci数列的前几项为：1，1，2，3，5，8，13，21，34，55，89，…。
该数列有着非常广泛的应用。
Fibonacci数列所满足的差分方程为Fn-Fn-1-Fn-2=0,其特征方程为λ2-λ-1=0其根为λ1= ,λ2= .利用λ1λ2可将差分方程写为Fn-(λ1+λ2)Fn-1+λ1λ2Fn-2=0,即Fn-λ1Fn-1=λ2(Fn-1-λ1Fn-2)数列｛Fn-λ1Fn-1｝满足一个一阶差分方程．显然( )同理可得( )由以上两式可解出的通项。
练习9 证明若数列{ }满足二阶差分方程，其特征方程由两个不相等的根，则为该差分方程的两个特解。
从而其通解为。
由练习9，若二阶差分方程的特征方程有两个不相等的根，可写出其通解的一般性式。
再由的值可解出其中的系数，从而写出差分方程的特解。
练习10 具体求出Fibonacci数列的通项，并证明。
那么,若二阶线性齐次差分方程有两个相等的根，其解有如何来求呢？设二阶线性齐次差分方程的特征方程有两个相等的根，则差分方程可写为。
差分方程的两边同时除以，有。
设,则(n>=3)。
由于该式在n>=3式均成立,我们将它改写为(n>=1)。
(8.2)方程(8.2)的左边是的二阶差分,从而有,于是是n的一次函数,设为则有。
上是即为差分方程的通解。
练习11 证明：若数列{ } 所满足的三阶差分方程的特征方程由三个相等的根，则差分方程的通解为。
一般的，设•••, 为差分方程的特征方程所有不同的解,其重数分别为•••，，则差分方程对应于其中的根（i=1，2，•••，l）的特解••• 。
对于一般的k阶齐次线性差分方程，我们可以通过其特征方程得到上述形式的k个特解，进而得到差分方程的通解。
练习12 若数列{ } 满足差分方程且求{ }的通项。
例6 若实系数差分方程的根为虚数，则其解也是用虚数表示的，这给讨论问题带来不便。
差分方程xn-2xn-1+4xn-2=0的特征值为i.若x1=1,x2=3,由下面的程序易求出其特解为：xn=( )(1+ i)n+(－)(1－i)nClear[x1,x2,c1,c2,l1,l2,solution];x1=1;x2=3;solution=Solve[1^2-2l+4==0,1];l1=l/.solution[[1,1]];l2=l/.solution[[2,1]];c=Solve[{c1*l1+c2*l2==x1,c1*l1^2+c2*l2^2==x2},{c1,c2}];c1=Simplify[Re[c1]]+Simplify[Im[c1]]*I;c2=Simplify[Re[c2]]+Simplify[Im[c2]]*I;Print[“xn=(“,c1,”)(“,l1,”)^n+(“,c2,”)(“,l2,”)^n”]解的形式相当复杂，是否可以将它们用实数表示呢？设=rei ,则＝re ,我们可将（8．4）中的表达式改写为xn=re (2e )n+re (2e )\n=r=2r Cos( )=(2rCos )=可以看出，通项可以写成的形式．那么，与是不是差分方程的特解呢？练习13 验证与是差分方程(8.3)的特解.对于差分方程(8.3),我们找出了它的两个实型的特解,从而可以将通解表示成实数的形式.这一方法对于一般的方程也是成立的.练习14 设的两个特征值为 .证明该差分方程的通解可表示为 .练习15 用实数表示差分方程的特解.上次我们讨论了其次线性差分方程的求解方法.那么,非齐次线性差分方程是否可以化为齐次线性差分方程呢?练习16 若已知非齐次线性差分方程••• （8.5）的一个特解为求证:若令则满足齐次差分方程•••由练习16,若已知非齐次线性差分方程(8.5)的一个特解,就可以将它化为齐次线性差分方程. 显然方程(8.5)的最简单的形式为(其中p为常数),代入(8.5)得•••若••• 则有称p = 为非齐次线性差分方程（8.5）的平衡值。
在（8.5）中，令则有由，得.从而可将原来的非齐次线性差分方程化为齐次线性差分方程.如果方程（8.5）的平衡值不存在，可以将方程（8.5）中所有的n换为n+1,得到（8.6）方程（8.6）和（8.5）相减得.于是可将原来的非齐次线性差分方程化为高一阶的齐次线性差分方程.练习17 分别求差分方程及的通解.2.3 代数方程求根由Fibonacci数列的性质，我们可以用来逼近，用这一性质可以来计算的近似值。
一般地，对a>0，可以用构造差分方程的方法来求的近似值.对给定的正数a，设λ1= ，λ2= ，则λ1 ，λ2是方程λ2-2λ+(1+a)=0的根.该方程是差分方程的特征方程。
于是，选定,利用差分方程可以构造一个数列{ }.练习18 证明：若a>1,对任意的>0, >0,若≠ ，则按上述法构造的数列{ }满足.这样，我们得到了计算的一个方法：1．给定（作为误差控制），任取初始值，令n=1;2．若，则终止计算，输出结果；否则，令n ：=n+1,转第3步；3．令,转第2步.练习19 对a=1.5,10,12345,用上述方法求 .上述方法的收敛速度不够快，我们可以加以改进设整数u满足,令,则，是方程的两个根.练习20 根据上面的差分方程的构件数列{ x },使得.练习21 对练习19中的a，用上面的方法来计算，并比较两种方法的收敛速度.代数方程（8.7）是差分方程（8.1）的特征方程，是否可以用此差分方程来求解方程（8.7）呢？设方程（8.7）有k个互不相同的根满足，（8.8）则对应的差分方程的通解形式为.练习22 设方程（8.7）的根满足条件（8.8），任取初始值用差分方程（8.1）（取b=0）构造数列{ }.若通解中的系数≠0，证明：.利用练习22得到的结论，我们可以求多项式方程的绝对值最大的根.练习23 求方程的绝对值最大的根.事实上，若方程（8.7）的互不相同的根满足≥ ≥…≥(其重数分别为)，则练习22中的结论仍然成立.2.4 国民收入4 国民收入的稳定问题一个国家的国民收入可用于消费，再生产的投资等。
一般地说，消费与再生产投资都不应该没有限制。
合理的控制各部分投资，能够使国民经济处于一种良性循环之中。
如何配各部分投资的比例，才能使国民经济处于稳定状态呢？这就是本节要讨论的问题。
我们首先给出一些假设条件：1．国民收入用于消费、再生产投资和公共设施建设三部分。
2．记分别为第k个周期的国民收入水平和消费水平。
的值与前一个周期的国民收入成正比例。
即=A , (8.9)其中A为常数（0<A<1）.3．用表示第k个周期内用于再生产的投资水平，它取决于消费水平的变化，即 . (8.10) 4．G表示政府用于公共设施的开支，设G为常数.由假设1有 . (8.11)上式是一个差分方程，当给定的值后，可直接计算出国民收入水平（k=2,3,…）来观察其是否稳定。
例7 若，计算可得表8.3中数据。
表8.3 Y 的值的变化k 2 3 4 5 6 7 8 9 10 1111.0 24.5 35.8 39.1 32.9 20.3 7.48 0.95 3.93 15.0k 12 13 14 15 16 17 18 19 20 2128.5 37.8 38.2 29.5 16.0 4.58 0.82 6.65 19.2 32.1我们可以画出的散点图来观察其变化。
