看一个例子:
∑i=0ai2
现在我们尝试使用多种方法来求这个式子的封闭形式——即通项公式。
方法0:可以查找公式。最好的资料来源是Sloane所著的《Handbook of Integer Sequences》,里面列出了数以千计的序列的封闭形式。可以去参考。还有些软件如:Axiom、MACSYMA、Maple或者Mathematica可以处理一些庞大的公式。其中,Mathematica很有名,大家可以去试用一下。
方法1:打表,猜测答案,然后用归纳法法证明。
方法2:对和式用扰动法
按照一般的做法,我们在其后加上一项:
Sn+1=Sn+(n+1)2=0≤k≤n∑(k+1)2=0≤k≤n∑(k2+2k+1)=0≤k≤n∑k2+20≤k≤n∑k+n+1=Sn+2∑0≤k≤n+(n+1)
发现Sn消掉了,扰动法失效。
那我们尝试从整数的立方和出发,来用一下扰动法。
设G(n)表示整数的立方和。
G(n)+(n+1)3=0≤k≤n∑(k+1)3=0≤k≤n∑(k3+3k2+3k+1)=0≤k≤n∑k3+30≤k≤n∑k2+23n(n+1)+n+1
发现G(n)消掉了,里面刚好有S(n).
3S(n)=(n+1)3−23n(n+1)+n+1
化简可得:
3S(n)=(n+1)(n+21)n
方法三:建立成套方法
设递归式
R0Rn=α;=Rn−1+β+γn+δn2,n>0(1)
的解的一般形式是
Rn=A(n)α+B(n)β+C(n)γ+D(n)δ(2)
则
Sn对应于递归式取
α=0,β=0,γ=0,δ=1时
Rn的值。.
设
Rn=1,则可得
α=1,β=0,γ=0,δ=0,于是
A(n)=1
设
Rn=n,则可得
α=0,β=1,γ=0,δ=0,于是
B(n)=n
设
Rn=n2,则可得
α=0,β=−1,γ=2,δ=0,于是
C(n)=2n2+n
设
Rn=n3,则可得
α=0,β=1,γ=−3,δ=3,于是
n3=B(n)−3C(n)+3D(n),有
D(n)=62n3−2n+3n(n+1)=6n(n+1)(2n+1)
现在
A(n),B(n),C(n),D(n)都已经求出来了。将
α=0,β=0,γ=0,δ=1代入式子
(2)中,则可求得
Sn=D(n)=6n(n+1)(2n+1)
成套方法是先从特殊到一般,再从一般到特殊的一个经典应用。
首先对一个具体的递归式,将其表示为一个系数待确定的方程,递归式变量为
α,β,γ,δ,各变量的系数分别为
A(n),B(n),C(n),D(n),该方程是一个通用的方程,可以表示符合这种形式的任意递归式(f(n)与f(n-1)的倍数关系不能变)。然后再从一般到特殊,通过一些特殊的数列,求出系数的值。于是该方程的编程了一个系数确定的方程。只要代入
α,β,γ,δ的值,就可以确定方程在那一组变量下的值,即为原递归式的封闭形式。
方法四:用积分替换和式。
∑0≤k≤nk2可以看做宽度为1,高度为
k2的小矩形的面积之和,它可以近似的表示为:
∫0nx2dx=n3/3
我们检查这个误差,设误差
En=Sn−n3/3.
因为Sn满足递归式,所以有:
En=Sn−n3/3=Sn−1+n2−n3/3=En−1+(n−1)3/3+n2−n3/3=En−1+n−1/3
可以很容易计算出
En=(n+1)n/2−n/3
则
Sn=En+n3/3=(n+1)n/2−n/3+n3/3=6n(n+1)(2n+1)
方法五:展开和收缩
Sn=1≤k≤n∑k2=1≤j≤k≤n∑k=1≤j≤n∑j≤k≤n∑k=1≤j≤n∑(2j+n)(n−j+1)=211≤j≤n∑(n(n+1)+j−j2)=21n2(n+1)+41n(n+1)−21Sn=21n(n+21)(n+1)−21Sn
从单重和式变换到多重和式,初看是一种倒推,但实际上是进步。
我们不能指望不断的简化简化再简化,就像你不可能只向上爬坡而登上最高的山峰!
方法六:有限微积分
方法七:无限微积分