递归方程组解的渐进阶的求法——差分方程法

原创 2007年09月17日 12:56:00

递归方程组解的渐进阶的求法——差分方程法

这里只考虑形如:

T(n)=c1T(n-1)+c2T(n-2)+…+ ckT(n-k)+f(n),nk (6.18)

的递归方程。其中ci (i=l,2,…,k)为实常数,且ck≠0。它可改写为一个线性常系数k阶非齐次的差分方程:

T(n)-c1T(n-1)- c2T(n-2)-…-ckT(n-k)=f(n),nk (6.19)

(6.19)与线性常系数k阶非齐次常微分方程的结构十分相似,因而解法类同。限于篇幅,这里直接给出(6.19)的解法,略去其正确性的证明。

第一步,求(6.19)所对应的齐次方程:

T(n)-c1T(n-1)- c2T(n-2)-…-ckT(n-k)=0 (6.20)

的基本解系:写出(6.20)的特征方程:

C(t)=tk-c1tk-1-c2tk-2 -…-ck=0 (6.21)

t=r是(6.21)的m重实根,则得(6.20)的m个基础解rnnrnn2rn,…,nm-1rn;若ρeρe-是(6.21)的一对l重的共扼复根,则得(6.20)的2l个基础解ρncosρnsinncosnsin,…,nl-1ρncosnl-1ρncos。如此,求出(6.21)的所有的根,就可以得到(6.20)的k个的基础解。而且,这k个基础解构成了(6.20)的基础解系。即(6.20)的任意一个解都可以表示成这k个基础解的线性组合。

第二步,求(6.19)的一个特解。理论上,(6.19)的特解可以用Lagrange常数变易法得到。但其中要用到(6.20)的通解的显式表达,即(6.20)的基础解系的线性组合,十分麻烦。因此在实际中,常常采用试探法,也就是根据f(n)的特点推测特解的形式,留下若干可调的常数,将推测解代人(6.19)后确定。由于(6.19)的特殊性,可以利用迭加原理,将f(n)线性分解为若干个单项之和并求出各单项相应的特解,然后迭加便得到f(n)相应的特解。这使得试探法更为有效。为了方便,这里对三种特殊形式的f(n),给出(6.19)的相应特解并列在表6-1中,可供直接套用。其中pi,i=1,2,…,s是待定常数。

表6-1 方程(6.19)的常用特解形式

f(n)的形式 条    件 方程(6.19)的特解的形式
an C(a)≠0
aC(t)的m重根
ns C(1)≠0
1是C(t)的m重根
nsan C(a)≠0
aC(t)的m重根

第三步,写出(6.19)即(6.18)的通解

(6.22)

其中{Ti(n),i=0,1,2,…,n}是(6.20)的基础解系,g(n)是(6.19)的一个特解。然后由(6.18)的初始条件

T(i)=Ti ,i=1,2,…,k-1

来确定(6.22)中的待定的组合常数{ai},即依靠线性方程组

解出{ai},并代回(6.22)。其中βj=Tj-g(j),j=0,1,2,…,k-1。

第四步,估计(6.22)的渐近阶,即为所要求。

下面用两个例子加以说明。

例l 考虑递归方程

它的相应特征方程为:

C(t)=t2-t-1=0

解之得两个单根。相应的(6.20)的基础解系为{r0nr1n}。相应的(6.19)的一个特解为F*(n)=-8,因而相应的(6.19)的通解为:

F(n)=a0r0n +a1r1n- 8

令其满足初始条件,得二阶线性方程组:

解之得,从而

于是

例2 考虑递归方程

T(n)=4T(n-1)-4T(n-2)+2nn (6.23)

和初始条件T(0)=0,T(1)=4/3。

它对应的特征方程(6.21)为

C(t)=t2-4t+4=0

有一个两重根r =2。故相应的(6.20)的基础解系为{2n,2nn}。由于f(n)=2nn,利用表6-1,相应的(6.19)的一个特解为

T*(n)=n2(p0+p1n)2n

代人(6.23),定出p0=1/2,p1=1/6。因此相应的(6.19)的通解为:

T(n)=a02n+a1n2n+n2(1/2+n/6)2n

令其满足初始条件得a0=a1=0,从而

T(n)=n2(1/2+n/6)2n

于是T(n)=θ(n32n)。

 

相关文章推荐

递归方程的求解

递归方程的求解 生成函数求解递归方程 生成函数及其性质 生成函数求递归方程 用特征方程求解递归方程 用递推法求解递归方程 用递推法求解常系数递归方程 用递推法求解变系数递归方程 函数换名 递归方程的求...

递归算法时间复杂度的计算方程式一个递归方程:

算法导论上的解递推公式的主函数法 递归算法时间复杂度的计算方程式一个递归方程:     在引入递归树之前可以考虑一个例子:   T(n) = 2T(n/2) + n2   迭代2次可...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

解递归式的方法总结

算法设计中经常会用到递归,利用递归式的方法可以清晰地显示算法的整个过程,而对于分析算法的复杂度,解递归式就有了用处,这里的方法来自于《算法导论》。 (一)代换法: 实质上就是数学归纳法,先对一个小...

Delphi7高级应用开发随书源码

  • 2003年04月30日 00:00
  • 676KB
  • 下载

算法设计关于递归方程T(n)=aT(n/b)+f(n)之通用解法

在算法设计中经常需要通过递归方程估计算法的时间复杂度T(n),本文针对形如T(n)=aT(n/b)+f(n)的递归方程进行讨论,以期望找出通用的递归方程的求解方式。 算法设计教材中给出的Master...
  • lsp1991
  • lsp1991
  • 2014年10月08日 10:57
  • 1783

递归算法的时间复杂度求解

在算法分析中,当一个算法中包含递归调用时,其时间复杂度的分析会转化为一个递归方程求解。实际上,这个问题是数学上求解渐近阶的问题,而递归方程的形式多种多样,其求解方法也是不一而足,比较常用的有以下四种方...
  • scutth
  • scutth
  • 2011年09月29日 12:44
  • 3351

差分方程建模

差分方程建模: 1:差分方程:针对的是离散点 2:微分方程:针对的是连续型的 一般是已知差分序列:y1-y0,y2-y1,y3-y2....... 目标是求:未知序列:y0,y1,y2....... ...

C语言结构体声明中冒号的使用(占位符) & C结构体的乱序初始化

有些信息在存储时,并不需要占用一个完整的字节, 而只需占几个或一个二进制位。例如在存放一个开关量时,只有0和1 两种状态, 用一位二进位即可。为了节省存储空间,并使处理简便,C语言又提供了一种数据结构...

MySQL 数据库的启动与关闭

MySQL数据库服务器通常指的的是mysqld,而命令行mysql则是mysql客户端程序,这两个概念通常容易混淆。通常启动mysql服务器即是启动mysqld进程,mysqld启动后,可以通过mys...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:递归方程组解的渐进阶的求法——差分方程法
举报原因:
原因补充:

(最多只允许输入30个字)