实验八__差分方程

1、实验八 差分方程实验目的1. 掌握差分的性质，多项式求和；2. 差分方程的解法；3. 用差分方程解代数方程； 4. 用差分方程分析国民经济。1 基本理论1. 差分2. 任意数列xn ,定义差分算子如下：xn=xn+1-xn对新数列再应用差分算子，有2xn=(kxn).性质性质1 k(xn+yn)=kxn+kyn性质2 k(cxn)=ckxn性质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,

2、a2,-ak 为常数， ak0. 若b=0,则该 方程是齐次方程关于 的代数方程k-a1k-1-ak-1-ak=0为对应的特征方程，根为特征值。 1 实验内容与练习21 差分例1 Xn=n3,求各阶差分数列：xnxn2xn3xn4xn17126081918602737246064613061259136216127343可见，n3,三阶差分数列为常数数列，四阶为0。 练习1 对1，n,n2,n4,n5, 分别求各阶差分数列。练习2 C0n-1C1n-1C2n-1,C4n-1,分别求各阶差分数列. Xn的通项为n的三次函数，Xn=a3n3+a2n2+a1n+a0证明它为常数数列。证明 由Xn=a

3、3n3+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 表SnSn2Sn3Sn4Sn5Sn181918609273724603664613060100125913660225216127424413431697845121296设Sn=a4n4+a3n3+a2n2+a1n+a0, s1=1,s2=9,s3=36,s4=100,s5=

4、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项

5、的取值，则可以确定通解的任意常数，得到差分 的特解。例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-

6、2=0,其特征方程为2-1=0其根为1=,2=.利用可将差分方程写为Fn-(1+2)Fn-1+12Fn-2=0,即Fn-1Fn-1=2(Fn-1-1Fn-2)数列Fn-1Fn-1满足一个一阶差分方程显然() 同理可得 () 由以上两式可解出的通项。 练习 证明若数列满足二阶差分方程，其特征方程由两个不相等的根，则为该差分方程的两个特解。从而其通解为。由练习9，若二阶差分方程的特征方程有两个不相等的根，可写出其通解的一般性式。再由的值可解出其中的系数，从而写出差分方程的特解。练习10 具体求出Fibonacci数列的通项，并证明。那么,若二阶线性齐次差分方程有两个相等的根，其解有如何来求呢？设二

7、阶线性齐次差分方程的特征方程有两个相等的根，则差分方程可写为。差分方程的两边同时除以，有。设,则(n=3)。由于该式在 n=3式均成立,我们将它改写为(n=1)。 (8.2) 方程(8.2)的左边是的二阶差分,从而有,于是是n的 一次函数,设为则有。上是即为差分方程的通解。 练习11 证明：若数列 所满足的三阶差分方程的特征方程由三个相等的根，则差分方程的通解为。 一般的，设,为差分方程的特征方程所有不同的解,其重数分别为，则差分方程对应于其中的根（i=1，2，l）的特解。 对于一般的k阶齐次线性差分方程，我们可以通过其特征方程得到上述形式的k个特解，进而得到差分方程的通解。 练习12 若数列

8、满足差分方程 且求的通项。例6 若实系数差分方程的根为虚数，则其解也是用虚数表示的，这给讨论问题带来不便。差分方程 xn-2xn-1+4xn-2=0的特征值为i.若x1=1,x2=3,由下面的程序易求出其特解为： xn=()(1+i)n+()(1i)n Clearx1,x2,c1,c2,l1,l2,solution; x1=1;x2=3; solution=Solve12-2l+4=0,1; l1=l/.solution1,1; l2=l/.solution2,1; c=Solvec1*l1+c2*l2=x1,c1*l12+c2*l22=x2,c1,c2; c1=SimplifyRec1+Si

9、mplifyImc1*I; c2=SimplifyRec2+SimplifyImc2*I; Print“xn=(“,c1,”)(“,l1,”)n+(“,c2,”)(“,l2,”)n”解的形式相当复杂，是否可以将它们用实数表示呢？设=rei,则re,我们可将（）中的表达式改写为xn=re (2e)n+re (2e)n =r =2rCos() =(2rCos) =可以看出，通项可以写成的形式那么，与是不是差分方程的特解呢？练习13 验证与是差分方程(8.3)的特解.对于差分方程(8.3),我们找出了它的两个实型的特解,从而可以将通解表示成实数的形式.这一方法对于一般的方程也是成立的.练习14 设的

10、两个特征值为.证明该差分方程的通解可表示为.练习 15 用实数表示差分方程的特解.上次我们讨论了其次线性差分方程的求解方法.那么,非齐次线性差分方程是否可以化为齐次线性差分方程呢?练习16 若已知非齐次线性差分方程 （.5）的一个特解为求证:若令则满足齐次差分方程 由练习16,若已知非齐次线性差分方程(8.5)的一个特解,就可以将它化为齐次线性差分方程.显然方程(8.5)的最简单的形式为(其中p为常数),代入(8.5)得 若则有称p=为非齐次线性差分方程（8.5）的平衡值。在（8.5）中， 令则有由，得 .从而可将原来的非齐次线性差分方程化为齐次线性差分方程.如果方程（8.5）的平衡值不存在，

11、可以将方程（8.5）中所有的n换为n+1,得到 （8.6）方程（8.6）和（8.5）相减得.于是可将原来的非齐次线性差分方程化为高一阶的齐次线性差分方程.练习17 分别求差分方程及的通解.2.3 代数方程求根由 Fibonacci数列的性质，我们可以用来逼近，用这一性质可以来计算的近似值。一般地，对a0，可以用构造差分方程的方法来求的近似值. 对给定的正数a，设1=，2= ，则1 ，2是方程2-2+(1+a)=0的根.该方程是差分方程 的特征方程。于是，选定,利用差分方程可以构造一个数列.练习 18 证明：若a1,对任意的0, 0,若，则按上述法构造的数列满足.这样，我们得到了计算的一个方法：

12、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 设

13、方程（8.7）的根满足条件（8.8），任取初始值用差分方程（8.1）（取b=0）构造数列.若通解中的系数0，证明： .利用练习22得到的结论，我们可以求多项式方程的绝对值最大的根.练习 23 求方程的绝对值最大的根.事实上，若方程（8.7）的互不相同的根满足 (其重数分别为)，则练习22中的结论仍然成立.2.4 国民收入4 国民收入的稳定问题 一个国家的国民收入可用于消费，再生产的投资等。一般地说，消费与再生产投资都不应该没有限制。合理的控制各部分投资，能够使国民经济处于一种良性循环之中。如何配各部分投资的比例，才能使国民经济处于稳定状态呢？这就是本节要讨论的问题。 我们首先给出一些假设条件：

14、1 国民收入用于消费、再生产投资和公共设施建设三部分。2 记分别为第k个周期的国民收入水平和消费水平。的值与前一个周期的国民收入成正比例。即 =A, (8.9)其中A为常数（0A1）.3 用表示第k个周期内用于再生产的投资水平，它取决于消费水平的变化，即. (8.10)4 G表示政府用于公共设施的开支，设G为常数.由假设1有. (8.11)上式是一个差分方程，当给定的值后，可直接计算出国民收入水平（k=2,3,）来观察其是否稳定。例7 若，计算可得表8.3中数据。 表8.3 Y的值的变化k23456789101111.024.535.839.132.920.37.480.953.9315.0k

15、1213141516171819202128.537.838.229.516.04.580.826.6519.232.1我们可以画出的散点图来观察其变化。其计算及画图的程序如下：y0=2;y1=2;a=0.5;b=2;g=10;y=y0,y1;Fork=1,k=20,k+,Y2=a(1+b)*y1-b*a*y0+g;Y=Appendy,y2;Y0=y1,y1=y2YListPloty,PlotJoinedTrue, PlotStyleThickness0.012 图.国民收入的变化由图.利用发现，又例的数据得出的的呈现出周期变化的迹象。练习设，对于表.中的参数A,B,分别计算（k=2,3,）并

16、画图观察的变化。表.参数A, B的取值A1/21/21/28/99/103/44/5B1231/21/233可以看出，随着参数的值不同，国民收入水平（k=2,3,）的稳定性呈现出不同的状态。那么，参数满足什么条件时，国民收入水平才处于稳定发展之中呢？差分方程（8.11）是一个常系数非齐次线性差分方程。由A1容易求出其平衡值为令可得.其特征值为 若则其中为的幅角。从而可的差分方程的解为其中为常数。若易见为一周期函数在-的取值，从而呈周期变化的状态。正如在例7中所见到的。 练习25 若在及的情形下，讨论的变化趋势。国民收入会稳定发展吗？练习26 若，国民收入在什么条件下会稳定发展？ 本实验涉及的Mathema