Loading [MathJax]/jax/output/CommonHTML/autoload/mtable.js

「具体数学」二:和式

记号#

和式和递归式#

Sn=nk=0akS0=a0,Sn=Sn1+anSn=nk=0akS0=a0,Sn=Sn1+an

如此,就能用第一章中用封闭形式求解递归式的方法

例如:

R0=α,Rn=Rn1+β+γnR0=α,Rn=Rn1+β+γn

R1=α+β+γR1=α+β+γ

R2=α+2β+3γR2=α+2β+3γ

一般来说

Rn=A(n)α+B(n)β+C(n)γRn=A(n)α+B(n)β+C(n)γ

用成套方法

R1=1R1=1,则

1=α1=1+β+nγ1=α1=1+β+nγ

得出α=1β=γ=0α=1β=γ=0

从而A(n)=1A(n)=1

Rn=nRn=n得到α=0β=1γ=0α=0β=1γ=0

从而B(n)=nB(n)=n

Rn=n2Rn=n2得到α=0β=1γ=2α=0β=1γ=2

从而2C(n)B(n)=n22C(n)B(n)=n2

C(n)=n+n22C(n)=n+n22

Rn=α+nβ+n+n22γRn=α+nβ+n+n22γ

若想计算ni=0(a+bi)ni=0(a+bi)

则转成递归式其中α=β=aγ=bα=β=aγ=b,解得Sn=a(n+1)+bn(n+1)2Sn=a(n+1)+bn(n+1)2

相反,有些递归式可以转成和式来求和

例如:

T0=0T0=0

Tn=2Tn1+1Tn=2Tn1+1

两边同除2n2n

T020=0T020=0

Tn2n=Tn12n1+12nTn2n=Tn12n1+12n

Sn=Tn2nSn=Tn2n,则

S0=0S0=0

Sn=Sn1+nnSn=Sn1+nn

得到和式

Sn=ni=12i=12nSn=ni=12i=12n

Tn=2nSn=2n1Tn=2nSn=2n1

一般地,形如

anTn=bnTn1+cnanTn=bnTn1+cn

转化为和式的思想是用一个求和因子snsn来乘两边

snanTn=snbnTn1+sncnsnanTn=snbnTn1+sncn

求和因子要恰当地选取使得

snbn=sn1an1snbn=sn1an1

这样一来,记Sn=snanTnSn=snanTn,则

Sn=Sn1+sncnSn=Sn1+sncn

Sn=s0a0T0+ni=1sici=s1b1T0+ni=1siciSn=s0a0T0+ni=1sici=s1b1T0+ni=1sici

Tn=1snan(s1b1T0+ni=1sici)Tn=1snan(s1b1T0+ni=1sici)

关于求snsn

展开sn=sn1an1/bnsn=sn1an1/bn

sn=kn1i=1aini=2bi,ksn=kn1i=1aini=2bi,k

在上述例子中,an=1,bn=2an=1,bn=2,所以sn=2nsn=2n

只要aabb不为0。。。

考虑解决快速排序比较步骤平均次数递归式

C0=C1=0C0=C1=0

Cn=n+1+2nn1i=0Ci,n>1Cn=n+1+2nn1i=0Ci,n>1

先规避掉除法

nCn=n2+n+2n1i=0CinCn=n2+n+2n1i=0Ci

n1n1代替nn

(n1)Cn1=(n1)2+(n1)+2n2i=0Ci(n1)Cn1=(n1)2+(n1)+2n2i=0Ci

两式相减

nCn(n1)Cn1=2n+2Cn1,n>2nCn(n1)Cn1=2n+2Cn1,n>2

关于CnCn的和式就变成了一个简单的递归式

C0=C1=0,C2=3C0=C1=0,C2=3

nCn=(n+1)Cn1+2n,n>2nCn=(n+1)Cn1+2n,n>2

根据求和因子的方法,递归式中an=n,bn=n+1,cn=2n2[n=1]+2[n=2]an=n,bn=n+1,cn=2n2[n=1]+2[n=2]

sn=kn1i=1aini=2bi=k(n1)!(n+1)!2!=k2n(n+1)sn=kn1i=1aini=2bi=k(n1)!(n+1)!2!=k2n(n+1)

k14k14

Cn=2(n+1)ni=11i+123(n+1)Cn=2(n+1)ni=11i+123(n+1)

其中

Hn=i1iHn=i1i

HnHn被称为调和级数,因为小提琴弦产生的第kk个泛音是弦长1k1k处所产生的基音

可以通过将CnCn表示为封闭形式研究快速排序递归式,CnCn中的和式是

nk=11k+1=1kn1k+1nk=11k+1=1kn1k+1

kk改为k1k1

1kn1k+1=1k1n1k1kn1k+1=1k1n1k

=2kn+11k=2kn+11k

=(1kn1k)1+1n+1=(1kn1k)1+1n+1

=Hnnn+1=Hnnn+1

Cn=2(n+1)Hn83n23Cn=2(n+1)Hn83n23

和式的处理#

分配律:

icai=ciaiicai=ciai

结合律:

i(ai+bi)=iai+ibii(ai+bi)=iai+ibi

交换律:

kKak=p(k)Kap(k),p(k)kKak=p(k)Kap(k),p(k)

假设要计算一个等差级数

S=nk=0(a+bk)S=nk=0(a+bk)

根据交换律,可以用nknk代替kk

S=nnk=0(a+b(nk))=nk=0(a+bnbk)S=nnk=0(a+b(nk))=nk=0(a+bnbk)

利用结合律将二者相加

2S=nk=0((a+bk)+(a+bnbk))=nk=0(2a+bn)2S=nk=0((a+bk)+(a+bnbk))=nk=0(2a+bn)

再利用分配律计算平凡和式

2S=(2a+bn)nk=01=(2a+bn)(n+1)2S=(2a+bn)nk=01=(2a+bn)(n+1)

nk=0(a+bk)=(a+12bk)(n+1)nk=0(a+bk)=(a+12bk)(n+1)

艾弗森约定#

kKak+kKak=kKKak+kKKakkKak+kKak=kKKak+kKKak

这是由一般公式

kKak=kak[kK]kKak=kak[kK]

[kK]+[kK]=[kKK]+[kKK][kK]+[kK]=[kKK]+[kKK]

推出的

扰动法#

将一项分出去的运算是扰动法的基础,利用扰动法,常常可以用封闭形式计算和式,思想是从一个未知的和式开始,记为SnSn

Sn=0knSn=0kn

然后,将第一项和最后一项分离

Sn+an+1=0kn+1ak=a0+1kn+1akSn+an+1=0kn+1ak=a0+1kn+1ak

=a0+1k+1n+1ak+1=a0+1k+1n+1ak+1

=a0+0knak+1=a0+0knak+1

如果能对最后的和式用SnSn表示,就会得到一个方程

例如用扰动法求一般几何级数的和

Sn=0knaxkSn=0knaxk

Sn+axn+1=ax0+nk=0axk+1=a0+xSnSn+axn+1=ax0+nk=0axk+1=a0+xSn

Sn=aaxn+11xSn=aaxn+11x

若是更难一点的和式

Sn=0knk2kSn=0knk2k

Sn+(n+1)2n+1=0kn(k+1)2k+1Sn+(n+1)2n+1=0kn(k+1)2k+1

借助结合律将右边和式分离

0knk2k+1+0kn2k+10knk2k+1+0kn2k+1

第一个和式等于2Sn2Sn,另一个等于2n+222n+22

Sn+(n+1)2n+1=2Sn+2n+22Sn+(n+1)2n+1=2Sn+2n+22

Sn=(n1)2n+1+2Sn=(n1)2n+1+2

xx代替22可以得到

Sn+(n+1)xn+1=xSn+xxn+21xSn+(n+1)xn+1=xSn+xxn+21x

Sn=x(n+1)xn+1+nxn+2(1x)2,x1Sn=x(n+1)xn+1+nxn+2(1x)2,x1

然而如果用微积分的初等技巧来推导这个封闭形式

nk=0xk=1xn+11xnk=0xk=1xn+11x

两边同时求导

nk=0kxk1=(1x)((n+1)xn)+1xn+1(1x)2=1(n+1)xn+nxn+1(1x)2nk=0kxk1=(1x)((n+1)xn)+1xn+1(1x)2=1(n+1)xn+nxn+1(1x)2

会得到同样的结论

因为和式的导数等于它的各项导数之和,微积分与离散数学有很多联系

多重和式#

一个和式的项可以用两个或多个指标来控制

可以用艾弗森约定表示

P(j,k)aj,k=j,kaj,k[P(j,k)]P(j,k)aj,k=j,kaj,k[P(j,k)]

有时也会用两个求和符号表示

jkaj,k[P(j,k)]jkaj,k[P(j,k)]

对于多重和式,我们有交换求和次序的基本法则,它推广了和式的结合律

jkaj,k[P(j,k)]=kjaj,k[P(j,k)]jkaj,k[P(j,k)]=kjaj,k[P(j,k)]

通常先对某个指标求和更容易些

1j,k3ajbk=1j31k3ajbk=(1j3aj)(1k3bk)1j,k3ajbk=1j31k3ajbk=(1j3aj)(1k3bk)

交换求和次序的基本法则有很多变形,如简易型:

jJkKaj,k=kKjJaj,kjJkKaj,k=kKjJaj,k

复杂型,适用于内和的范围与外和的指标变量有关的情形

jJkK(j)aj,k=kKjJ(k)aj,kjJkK(j)aj,k=kKjJ(k)aj,k

这里须满足

[jJ][kK(j)]=[kK][jJ(k)][jJ][kK(j)]=[kK][jJ(k)]

考虑交换求和次序的应用

有由n2n2个乘积ajakajak组成的阵列

[a1,1a1,nan,1an,n]

目的是求Sn=1jknajak

即这个阵列主对角线及其上方所有元素之和

由于ajak=akaj,所以Sn近似等于所有元素和一半,除了对主对角线计算时加入的修正因子

启发我们想到

Sn=1jknajak=1kjnakaj=Sn

由于

[1jkn]+[1kjn]=[1j,kn]+[1j=kn]

所以

2Sn=1j,knajak+1j=knajak

Sn=12((nk=1ak)2+nk=1a2k)

研究另一相似二重和式

S=1j<kn(akaj)(bkbj)

交换k,j,

=1k<jn(ajak)(bjbk)=1k<jn(akaj)(bkaj)

所以可以和自己求和

[1j<kn]+[1k<jn]=[1j,kn][1j=kn]

2S=1j,kn(ajak)(bjbk)1j=kn(ajak)(bjbk)

第二个和式显然为0,第一个和式可以展开

1j,knajbj1j,knaj,bk1j,knak,bj+1j,knakbk

=21j,knak,bk21j,knajbk

=2n1knakbk2(nk=1ak)(nk=1bk)

抛开原题,我们得到了一个有趣的公式:

(nk=1ak)(nk=1bk)=nnk=1akbk1j<kn(akaj)(bkbj)

这是切比雪夫单调不等式的一个特例

(nk=1ak)(nk=1bk)nnk=1ak,bk,a1anb1bn

(nk=1ak)(nk=1bk)nnk=1akbk,a1anb1bn

根据交换律可知

kKak=p(k)Kak

如果f是一个任意函数

f:JK

将整数jJ变成整数f(j)K,那么用f(j)替换k会发生什么?

指标替换的一般公式是

jJaf(j)=kKak#f(k)

这里#f(k)表示集合

f(k)={j|f(j)=k}

中的元素个数,即使得f(j)=kjJ的值的个数

由于jJ[f(j)=k]=#f(k),容易通过交换求和次序证明

jJaf(j)=jJ,kKak[f(j)=k]=kKakjJ[f(j)=k]

来看一个包含实际数字的多重和式:

Sn=1j<kn1kj

计算二重和式的正规方法是先对j或者先对k求和,分别讨论一下

Sn=1kn1j<k1kjj

=1kn1kj<k1jkjj

=1kn0<jk11jj

=1knHk1

=0k<nHk

然而调和数的和并不会表示为封闭形式!

尝试先对k求和会发现

Sn=0j<nHj

仍然是死胡同

若用kj替换k

Sn=1j<kn1kj

=1j<k+jn1k

=1kn1jnk1k

=1knnkk

=1knnk1kn1

=nHnn

这样就算是解决了吧(

在跟前面的死胡同结合起来,我们还得到了一个恒等式:

0k<nHk=nHnn

从代数角度理解这一技巧:

如果我们又一个包含k+f(j)的二重和式,其中f是一个任意的函数,那么这个例子指出,用kf(j)替换k并对j求和比较好

从几何角度理解:

nk=1k=2k=3k=4j=10111213j=2001112j=300011j=40000

我们首先尝试,对j求和(按行)和对k求和(按列),给出H1+H2+H3=H3+H2+H1

获得成功的方案本质上是按对角线求和31+22+13

一般性方法#

若要求

S(2)n=0knk2,n0

方法零:查公式#

S(2)n=n(n+1)(2n+1)6

方法一:归纳法#

假设我们不知道从哪得知了猜想

S(2)n=n(n+1)(n+2)6

n=0时显然成立

那么当n成立时,对n+1

S(2)n+1=S(2)n+(n+1)2

=n(n+1)(2n+1)6+6(n+1)26

=(n+1)(n+2)(2n+3)6

但还是太依赖神来猜想了

方法二:扰动法#

S(2)n+(n+1)2=0kn(k+1)2=0kn(k2+2k+1)

=0knk2+20knk+0kn1

=S(2)n+20knk+(n+1)

不幸的是S(2)n抵消了

不过我们得到了

20knk=(n+1)2(n+1)

显然0knk=S(1)n

S(2)n扰动得到了S(1)n,如果对S(3)n扰动呢?

S(3)n+(n+1)3=0kn(k+1)3=0kn(k3+3k2+3k+1)

=S(3)n+3S(2)n+3S(1)n+n+1

显然我们已经得到了S(2)n

方法三:成套方法#

列出递归式

R0=a

Rn=Rn1+b+cn+dn2

解的一般形式是

Rn=A(n)a+B(n)b+C(n)c+D(n)d

代入Rn=1

a=1,b=c=d=0

A(n)=1

代入Rn=n

a=0,b=1,c=d=0

B(n)=n

代入Rn=n2

a=0,b=1,c=2,d=0

2C(n)B(n)=n2

代入Rn=n3

a=0,b=1,c=3,d=3

3D(n)3C(n)+B(n)=n3

后略

方法四:积分替换和式#

由积分替换和式得到

S(2)nn0x2dx=n33

接下来检查近似的误差

En=S(2)nn33

=S(n1)+n2n33

=En1+13(n1)3+n2n33

=En1+n13

然后就几乎是白给了

另一种求误差方法是对每块的楔形面积求和

En=S(2)nn0x2dx=nk=1(k2kk1x2dx)

=nk=1(k2k3(k1)33)=nk=1(k13)

方法五:展开和放缩#

用一个看起来更复杂的二重和式替换原来的和式,二重和式操作得到(超神)是可以化简的

S(2)n=1knk2=1jknk

=1jnjknk

=1jnj+n2(nj+1)

=121jn(n(n+1)+jj2)

=12n2(n+1)+14n(n+1)12S(2)n

后略

方法六:有限微积分#

方法七:生成函数#

有限微积分和无限微积分#

类比传统的无限微积分,数学家们发展了有限微积分

无限微积分基于由

Df(x)=limh0f(x+h)f(x)h

定义的微分算子D的性质,有限微积分则是基于

Δf(x)=f(x+1)f(x)

所定义的差分算子Δ的性质,差分算子是微分的有限模拟,其中限制了取h的正整数值,于是h0时,1是能达到的极限,而Δf(x)则是f(x+h)f(x)hh=1时的值

符号DΔ称为算子,因为它们作用在函数上给出新的函数,它们是函数的函数,用于产生函数,如果f是从实数到实数的一个适当光滑的函数,那么Df也是一个从实数到实数的函数,又如果f是任意一个从实数到实数的函数,那么Δf亦然,函数DfΔf在点x处的值由上面的定义给出

我们知道对于f(x)=xm,Df(x)=mxm1,非正式地可以写成D(xm)=mxm1

可惜

Δ(x3)=(x+1)3x3=3x2+3x+1

但有一类m次幂,它在Δ的作用下可以很好地变换,这就是有限微积分的意义,那就是下降幂和上升幂

xm_=x(x1)(xm+1),m0

x¯m=x(x+1)(x+m1),m0

对他们有

Δ(xm_)=(x+1)m_xm_=mxm1_

无限微分算子有一个逆运算,即逆微分算子(或积分算子),微积分基本定理把D联系起来

g(x)=Df(x)g(x)dx=f(x)+C

这里g(x)dxg(x)的不定积分,它是导数等于g(x)的一个函数类,类似地,Δ也有一个逆运算,即逆差分算子(或求和算子),而且有一个类似的基本定理:

f(x)=Δf(x)g(x)δx=f(x)+C

这里g(x)δxg(x)的不定和式,是差分等于g(x)的一个函数类,不定积分中的C是一个任意常数,而不定和中的C则是满足p(x+1)=p(x)的任意一个函数p(x),例如C可以是周期函数a+bsin2πx,当取差分时会被消去

无限微积分也有定积分:

bag(x)dx=f(x)|ba=f(b)f(a)

于是有限积分,也有确定的和式

ba$g(x)δx=f(x)|ba=f(b)f(a)$$bag(x)δx$$$aag(x)δx=f(a)f(a)=0

如果b=a+1,结果是

a+1ag(x)δx=f(a+1)f(a)=g(a)

更一般地,如果b增加1

b+1ag(x)δxbag(x)δx=g(b)

通过观察和数学归纳,我们可以得出

bag(x)δx=ak<bg(k)

也就是说定和式和通常的有界和式是相同的

假如给定一个未知的和式,希望用封闭形式对它的值进行计算,并假设它可以写成形式ak<bg(k)=bag(x)δx,只要能求得一个补丁和式或者逆差分函数f,使得g(x)=f(x+1)f(x),那么由有限微积分的理论,就可以将答案表示成f(b)f(a),理解这个原理的一种方法是利用三点记号将ak<bg(k)完全写开

ak<b(f(k+1)f(k))=(f(a+1)f(a))+(f(a+2)f(a+1))++(f(b)f(b1))

除了f(b)f(a)以外,右边的每一项都抵消了,所以f(b)f(a)就是该和式的值,该种和式常被称为叠缩

但该法则仅ab时适用,若a>b

bag(x)δx=f(b)f(a)=(f(a)f(b))=abg(x)δx

类似方法可以论证

ba+cb=ca

确定的求和法为我们提供了一种计算下降幂的和式的简单方法

0k\<nkm_=km+1_m+1|n0=nm+1_m+1,n,m0

定和式的方法也暗示了,对0k<n求和通常比对1kn求和更简单,前者是f(n)f(0),后者是f(n+1)f(1)

通常幂也可以用新方法来求和,如果我们先用下降幂表示出通常幂

例如

k2=k2_+k1_

所以

0k<nk2=n3_3+n2_2=13n(n12)(n1)

n+1替换n又给出另一种将我们的老朋友S(2)n=0knk2计算成封闭形式的方法

从平方到立方

k3=k3_+3k2_+k1_

后面讨论的斯特林数可以做到通常幂和阶乘幂的转换

从而

ak<b=k4_4+k3_+k2_2|ba

下降幂还有性质,对应(x+y)2=x2+2xy+y2,也有(x+y)2_=x2_+2xy+y2_

实际上对(x+y)m(x+y)m_之间也有相似的结果

x2_=x(x1)

x1_=x

x0_=1

看上去除以x1会从x2_得到x1_,除以x会从x1_得到x0_,若接下来除以x+1会从x0_得到x1_

所以xm=1(x+1)(x+2)(x+m),m>0

类比通常幂的指数法则

xn+m=xnxm

下降幂的指数法则形式是

xn+m_=xm_(xm)n_

对负的n

x23_=x2_(x2)3_=x(x1)1(x1)x(x+1)=1x+1=x1_

我们来确定负数下降幂关于差分性质的成立

Δx2_=1(x+2)(x+3)1(x+1)(x+2)

=(x+1)(x+3)(x+1)(x+2)(x+3)

=2x3_

成立了

如此只要零不做除数

baxm_δx=xm+1_m+1|ba,m1

m=1时,回想积分

bax1dx=lnx|ba

再回到和式,我们希望找到f(x)使得

x1_=1x+1=Δf(x)=f(x+1)f(x)

x是整数时

f(x)=11+12++1x

所以f(x)即调和数Hx,于是Hx就是连续的lnx的离散模拟

现在能对下降幂的和式给出完整定义了

baxm_δx={xm_m+1|ba,m1Hx|ba,m=1

这个公式指出了为什么快速排序这样的离散问题解中会冒出调和数,恰如再连续性问题的解中自然会出现自然对数一样

既然我们找到了一个lnx的一个类似物,我们也来研究下是否有一个ex的类似物存在,什么样的函数f(x)有与恒等式D(ex)=ex对应的性质Df(x)=f(x)

f(x+1)f(x)=f(x)f(x+1)=2f(x)

可以取f(x)=2x作为离散指数函数

cx的差分也很简单

Δ(cx)=cx+1cx=(c1)cx

故如果c1,cx的逆差分是cxc1,这一事实与前面的基本法则结合起来,就给出了理解几何级数和的一般公式的方法

ak<bck=bacxδx=cxc1|ba=cbcac1,c1

每当遇到一个可能用作封闭形式的函数f,我们就能计算出它的差分Δf=g,然后就有一个函数g,它的不定和式g(x)δx是已知的

f=gΔf=gf=gΔf=g1x0_=102x2x2x1_=x1cx(c1)cx3x2_=x(x1)2xcxc1cx4xm_mxm1_cucΔu5xm+1_m+1xm_u+vΔu+Δv6Hxx1_=1x+1uvuΔv+EvΔu

但有些连续概念并没有离散的相似概念,例如无限微积分的链式法则是对函数的函数求导的方便法则,但有限微积分没有对应的链式法则,因为Δf(g(x))没有很好的形式

除了像用c±x替换x这样的情形之外,离散的变量变换很难

然而Δ(f(x)+g(x))的确有比较好的形式,而且提供了一个分布求和的法则,它是无限微积分中称为分部积分法则的一个有限相似结果

回顾以下无限微积分中的公式

D(uv)=uDv+vDu

在积分并重新排列各项的次序之后,它引导出分部积分法

uDv=uvvDu

可以在有限微积分中做到类似的事情

Δ(u(x)v(x))=u(x+1)v(x+1)u(x)v(x)

=u(x+1)v(x+1)u(x)v(x+1)+u(x)v(x+1)u(x)v(x)

=u(x)Δv(x)+v(x+1)Δu(x)

这个公式可以由Ef(x)=f(x+1)

所定义的移位算子E表示成方便的形式,用Ev(x)替换v(x+1)就得到乘积的差分的一个紧凑法则:

Δ(uv)=uΔv+EvΔu

在这个方程两边取不定和,并重新排列它的项,就得到广为宣扬的分部求和法则:

uΔv=uvEvΔu

(无限微积分在此令10)来避开了E

如同对无限微积分那样,可以在所有三项上加上界限,从而使得不定和式变成确定的和式

当左边的和式比右边的和式更难以处理时,这个法则是有用的,例如xexdx是分部积分的一个典型例子,它的离散模拟是x2xδx

为了用分部求和,我们用u(x)=x,Δv(x)=2x,从而Δu(x)=1,v(x)=2x,且Ev(x)=2x+1

x2xδx=x2x2xδx=x2x2x+1+C

加上界限,我们可以用此式来计算以前做过的和式

nk=0k2k=n+10x2xδx

=x2x2x+1|n+10

=((n+1)2n+12n+2)(02021)=(n1)2n+1+2

用这个方法比用扰动法更容易求出这个和,因为不用动脑子

在之前我们碰巧发现了0k<nHk的一个公式,但是如果知道分部求和,那么就能遵循规则发现公式

例如有和式0k<nkHk,如果用xlnxdx的模拟来引导,解就很好求了

u(x)=Hx,Δv(x)=x=x1_,从而Δu(x)=x1_,v(x)=x2_2,Ev(x)=(x+1)2_2

xHxδx=x2_2Hx(x+1)2_2x1_δx

=x2_xHx12x1_δx

=x2_xHxx2_4+C

第一行到第二行利用了指数法则将(x+1)2_x1_结合了起来

现在添加界限并得到结论

0k<nkHk=n0xHxδx=n2_2(Hn12)

无限和式#

对于无限和式,处理时的方法并不总是有效的,但存在一大类容易理解的无限和式,我们对它们做过的所有运算都是完全合法的,在我们更密切审视求和的基本意义之后,这两个消息背后所隐藏的理由将会变得清晰起来

无限和式可以很自然的这样定义:

S=1+12+14+18+

等于2,因为如果将它加倍就得到

2S=2+1+12+14+=2+S

另一方面,同样的推理提示我们应该定义

T=1+2+4+8+16+

1,因为

2T=2+4+8+16+32+=T1

笑死,其实应该说成T=,因为T中相加的项会大于任何指定的有限数

我们尝试对一般和式kKak的值构想一个好的定义,其中K可以是无限的

首先,让我们假设所有的项ak都是非负的,这样就不难找到一个合适的定义:如果有一个常数A为界,使得对所有有限的子集FK都有

kFakA

那么我们就定义kKak是最小的这样的A(它由实数众所周知的性质得出:所有这样的A总包含一个最小的元素)

如果没有常数A为界限,我们就说kKak=

上面的定义与指标集K中可能存在的次序无关,下面我们打算要讨论适用于带有多个指标k1,k2的多重和式,而不仅仅是指标集为(单一)整数集合的和式

K是非负整数集合的特殊情形下,我们对于非负项ak的定义就意味着

k0ak=limnnk=0ak

理由是:实数的任何一个非减序列都有极限(可能是∞),如果极限是A,又如果F是任意一个非负整数的有限集合,其元素全都n,我们就又kFaknk=0akA,因为A=或者A是一个有界常数

又如果A是任何一个小于所述极限A的数,那么就存在一个n,使得nk=0ak>A,该有限集合f={0,1,,n}证明A不是有界常数

按照刚刚的定义,我们很容易地计算某种无限和式的值,例如如果ak=xk就有

k0xk=limn1xn+11x={11x,0x<1,x1

特别地,刚刚考虑过的无限和ST分别取值2,恰如我们所猜测的,另一个有趣的例子是

k01(k+1)(k+2)=k0k2_

=limnn1k=0k2_=limnk1_1|n0=1

现在考虑和式中既可能有非负项又可能有负项的情形,例如

k0(1)k=11+11+11+

对他进行分组

k0(1)k=(11)+(11)+(11)+=0

但如果延迟一项再分组

k0(1)k=1(11)(11)(11)=1

我们或许也可以尝试再公式k0xk=11x中令x=1,因为我们已经证明了这个公式当0x<1时成立,但这样就迫使我们得出结论:这个无限和是12,尽管它是由整数组成的和式

另一个有趣的例子是双向无限的kak,其中当k0ak=1k+1,而当k<0时,ak=1k1,我们可以把它写成

+(14)+(13)+(12)+1+12+13+14+

如果我们从位于中心的元素开始往外计算这个和式

+(14+(13+(12+(1)+12)+13)+14)+=1

如果将所有括号左移一步

+15+(14+(13+(12)+1)+12)+13)+=1

因为最内层的n个括号中,所有数之和是

1n+11n12+1+12++1n1=11n1n1

类似的讨论声明,如果将这些括号向左或者向右任意移动固定步,它的值都是1,这使我有勇气相信这个和式的确是1

另一方面,如果我们以下述方式对项分组

+(14+(13+(12+1+12)+13+14)+15+16)+

那么从内往外第n对括号包含数

1n+11n12+1+12++12n1+12n=1+H2nHn+1

在后面将证明limn(H2nHn+1)=ln2

于是这一组合方式表面,这个双向无限和式实际上应该等于1+ln2

按照不同方式对其项相加而得出不同值的和式,有某些不同寻常之处。不同教材有不同定义

事实上,我们关于无限和式的定义相当简单,设K是任意一个集合,而ak是对每一个kK定义的实值项(这里k实际上可以代表若干个指标k1,k2,,因而K可以是多维的),任何实数x都可以写成其正的部分减去负的部分

x=x+x,x+=x[x>0],x=x[x<0]

(或者x+=0,或者x=0,或者两者皆成立),我们已经说明了怎样来定义无限和式kKa+kkKak的值,因为a+kak都是非负的,这样一来,一般性的定义是

kKak=kKa+k+kKak

除非右边的两个和式都等于,在后面这种情形,我们不定义kKak

A+=kKa+k,A=kKak,如果A+A两者都是有限的,就说和式kKak绝对收敛于值A=A+A。如果A+=,A是有限的,就说和式kKak发散于+,类似地,如果A=A+是有限的,就说kKak发散于,如果两者都是,结果还很难说

我们从非负项所做的定义开始,然后再将它推广到实值的项,如果ak是复数,我们可以再次用显然的方式将定义推广:和式kKak定义成kKak+ikKak,其中,akakak的实部和虚部,这两个和式都有定义。反之,没有定义

坏消息是某些无限和式必须不予定义,因为我们所做的操作再所有这样的情形中都可能产生矛盾,好消息是,只要我们处理的是刚才所定义的绝对收敛的和式,这一章里所有操作都完全成立

通过证明每一个变换法则都保持所有绝对收敛的和式的值不变,我们可以验证好消息

需要证明的主要结果是多重和式的基本原理:经过两个或多个指标集的绝对收敛的和式永远可以对这些指标中的任何一个首先求和

将来会正式地证明:如果J是任意的指标集,并且{Kj|jJ}的元素是任意的指标集,使得

jJ,kKkaj,kA

那么对每一个jJ都存在复数Aj,使得

kKjaj,kAj,

jJAjA

只要再所有的项都非负时证明时就够了,因为可以和前面一样把每一项分解成实部和虚部,正的部分和负的部分来证明一般情形

于是,我们假设对所有的指标对(j,k)M都有aj,k0,其中M是主指标集{(j,k)|jJ,kKj}

给定(j,k)Maj,k是有限的,即对所有有限子集FM

(j,k)Faj,kA

A这样的最小上界,如果jJ的任意一个元素,形如kFjaj,k的每一个和都以A为上界,其中FjKj的一个有限子集。从而这些有限和式有一个最小的上界Aj0,且根据定义有kKjaj,k=Aj

我们仍需要证明:对所有有限子集GJAjGAj的最小上界

假设GJ的满足jGAj=A>A的有限子集。我们至少可以求出一个有限子集FjKj,使得对每个满足Aj>0jG均有kFjaj,k>AAAj
但是此时有jG,kFjaj,k>AAjGAj=A,这与如下事实矛盾:对所有有限子集FM(j,k)Faj,kA。从而对所有有限子集Gj都有AjA

(我觉得上面全是废话。。)

最后,设A是小于A的任何一个实数。如果我们能找到一个有限集合GJ,使得jGAj>A,证明就完成了

我们知道存在一个有限集合fM,使得(j,k)Faj,k>A,设G是这个F中那些j组成的集合,又设Fj={k|(j,k)F},那么就有jGAjjGk Fjaj,k=(j,k)Faj,k>A

证明完毕

所以现在关于无限和式的一切法则证明完毕,只要它的项的绝对值组成的所有有限和式都有一个有限的界。

由于双向无限和式再我们用两种不同方法计算时给出两个不同的答案,因为它的正项1+12+13+必定发散于,否则,不论我们如何对它的项进行分组都会得到同一个答案

习题#

0

0k=4qk=k0qkk<4qk=q1q2q3

|x|

a0+a1+a2+a3+a4+a5

a4+a1+a0+a1+a4

a.

4i=14j=i+14k=j+1aijk=(((a123+a124)+a134)+a234)

b.

4k=1k1j=1j1i=1aijk=(a123+(a124+(a134+a234)))

两个不同的指标变量用了同一指标

[1jn](nj+1)

(xm_)=xm_(x1)m_

=x(x1)(x2)(xm+1)(x1)(x2)(x3)(xm)

=m(x1)m1_

0m_={0,m11|m|!,m0

x¯n+m=x¯n(x+n)¯m

xn=1(x1)n_

Δ(u(x)v(x))=u(x+1)v(x+1)u(x)v(x)

=u(x+1)v(x+1)u(x)v(x+1)+u(x)v(x+1)u(x)v(x)

=uΔv+EvΔu

右边还可能是

=EuΔv+vΔu

posted @   lovelyred  阅读(1713)  评论(0编辑  收藏  举报
编辑推荐:
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
阅读排行:
· 地球OL攻略 —— 某应届生求职总结
· 周边上新:园子的第一款马克杯温暖上架
· Open-Sora 2.0 重磅开源!
· 提示词工程——AI应用必不可少的技术
· .NET周刊【3月第1期 2025-03-02】
点击右上角即可分享
微信分享提示
CONTENTS