「具体数学」二:和式
记号#
略
和式和递归式#
如此,就能用第一章中用封闭形式求解递归式的方法
例如:
一般来说
用成套方法
令R1=1R1=1,则
得出α=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γ
若想计算n∑i=0(a+bi)n∑i=0(a+bi)
则转成递归式其中α=β=a,γ=bα=β=a,γ=b,解得Sn=a(n+1)+bn(n+1)2Sn=a(n+1)+bn(n+1)2
相反,有些递归式可以转成和式来求和
例如:
两边同除2n2n
令Sn=Tn2nSn=Tn2n,则
得到和式
则
一般地,形如
转化为和式的思想是用一个求和因子snsn来乘两边
求和因子要恰当地选取使得
这样一来,记Sn=snanTnSn=snanTn,则
关于求snsn
展开sn=sn−1an−1/bnsn=sn−1an−1/bn
在上述例子中,an=1,bn=2an=1,bn=2,所以sn=2−nsn=2−n
只要aa和bb不为0。。。
考虑解决快速排序比较步骤平均次数递归式
先规避掉除法
用n−1n−1代替nn
两式相减
关于CnCn的和式就变成了一个简单的递归式
根据求和因子的方法,递归式中an=n,bn=n+1,cn=2n−2[n=1]+2[n=2]an=n,bn=n+1,cn=2n−2[n=1]+2[n=2]
其中
HnHn被称为调和级数,因为小提琴弦产生的第kk个泛音是弦长1k1k处所产生的基音
可以通过将CnCn表示为封闭形式研究快速排序递归式,CnCn中的和式是
将kk改为k−1k−1
和式的处理#
分配律:
结合律:
交换律:
假设要计算一个等差级数
根据交换律,可以用n−kn−k代替kk
利用结合律将二者相加
再利用分配律计算平凡和式
即
艾弗森约定#
这是由一般公式
和
推出的
扰动法#
将一项分出去的运算是扰动法的基础,利用扰动法,常常可以用封闭形式计算和式,思想是从一个未知的和式开始,记为SnSn
然后,将第一项和最后一项分离
如果能对最后的和式用SnSn表示,就会得到一个方程
例如用扰动法求一般几何级数的和
若是更难一点的和式
借助结合律将右边和式分离
第一个和式等于2Sn2Sn,另一个等于2n+2−22n+2−2
用xx代替22可以得到
然而如果用微积分的初等技巧来推导这个封闭形式
两边同时求导
会得到同样的结论
因为和式的导数等于它的各项导数之和,微积分与离散数学有很多联系
多重和式#
一个和式的项可以用两个或多个指标来控制
可以用艾弗森约定表示
有时也会用两个求和符号表示
对于多重和式,我们有交换求和次序的基本法则,它推广了和式的结合律
通常先对某个指标求和更容易些
交换求和次序的基本法则有很多变形,如简易型:
复杂型,适用于内和的范围与外和的指标变量有关的情形
这里须满足
考虑交换求和次序的应用
有由n2n2个乘积ajakajak组成的阵列
目的是求Sn=∑1≤j≤k≤najak
即这个阵列主对角线及其上方所有元素之和
由于ajak=akaj,所以Sn近似等于所有元素和一半,除了对主对角线计算时加入的修正因子
启发我们想到
由于
所以
研究另一相似二重和式
交换k,j,
所以可以和自己求和
第二个和式显然为0,第一个和式可以展开
抛开原题,我们得到了一个有趣的公式:
这是切比雪夫单调不等式的一个特例
根据交换律可知
如果f是一个任意函数
将整数j∈J变成整数f(j)∈K,那么用f(j)替换k会发生什么?
指标替换的一般公式是
这里#f−(k)表示集合
中的元素个数,即使得f(j)=k的j∈J的值的个数
由于∑j∈J[f(j)=k]=#f−(k),容易通过交换求和次序证明
来看一个包含实际数字的多重和式:
计算二重和式的正规方法是先对j或者先对k求和,分别讨论一下
然而调和数的和并不会表示为封闭形式!
尝试先对k求和会发现
仍然是死胡同
若用kj替换k
这样就算是解决了吧(
在跟前面的死胡同结合起来,我们还得到了一个恒等式:
从代数角度理解这一技巧:
如果我们又一个包含k+f(j)的二重和式,其中f是一个任意的函数,那么这个例子指出,用k−f(j)替换k并对j求和比较好
从几何角度理解:
我们首先尝试,对j求和(按行)和对k求和(按列),给出H1+H2+H3=H3+H2+H1
获得成功的方案本质上是按对角线求和31+22+13
一般性方法#
若要求
方法零:查公式#
方法一:归纳法#
假设我们不知道从哪得知了猜想
n=0时显然成立
那么当n成立时,对n+1有
但还是太依赖神来猜想了
方法二:扰动法#
不幸的是S(2)n抵消了
不过我们得到了
显然∑0≤k≤nk=S(1)n
对S(2)n扰动得到了S(1)n,如果对S(3)n扰动呢?
显然我们已经得到了S(2)n
方法三:成套方法#
列出递归式
解的一般形式是
代入Rn=1
代入Rn=n
代入Rn=n2
代入Rn=n3
a=0,b=1,c=−3,d=3
后略
方法四:积分替换和式#
由积分替换和式得到
接下来检查近似的误差
然后就几乎是白给了
另一种求误差方法是对每块的楔形面积求和
方法五:展开和放缩#
用一个看起来更复杂的二重和式替换原来的和式,二重和式操作得到(超神)是可以化简的
后略
方法六:有限微积分#
方法七:生成函数#
有限微积分和无限微积分#
类比传统的无限微积分,数学家们发展了有限微积分
无限微积分基于由
定义的微分算子D的性质,有限微积分则是基于
所定义的差分算子Δ的性质,差分算子是微分的有限模拟,其中限制了取h的正整数值,于是h→0时,1是能达到的极限,而Δf(x)则是f(x+h)−f(x)h当h=1时的值
符号D和Δ称为算子,因为它们作用在函数上给出新的函数,它们是函数的函数,用于产生函数,如果f是从实数到实数的一个适当光滑的函数,那么Df也是一个从实数到实数的函数,又如果f是任意一个从实数到实数的函数,那么Δf亦然,函数Df和Δf在点x处的值由上面的定义给出
我们知道对于f(x)=xm,Df(x)=mxm−1,非正式地可以写成D(xm)=mxm−1
可惜
但有一类m次幂,它在Δ的作用下可以很好地变换,这就是有限微积分的意义,那就是下降幂和上升幂
对他们有
无限微分算子有一个逆运算,即逆微分算子(或积分算子)∫,微积分基本定理把D和∫联系起来
这里∫g(x)dx是g(x)的不定积分,它是导数等于g(x)的一个函数类,类似地,Δ也有一个逆运算,即逆差分算子(或求和算子)∑,而且有一个类似的基本定理:
这里∑g(x)δx是g(x)的不定和式,是差分等于g(x)的一个函数类,不定积分中的C是一个任意常数,而不定和中的C则是满足p(x+1)=p(x)的任意一个函数p(x),例如C可以是周期函数a+bsin2πx,当取差分时会被消去
无限微积分也有定积分:
于是有限积分,也有确定的和式
如果b=a+1,结果是
更一般地,如果b增加1
通过观察和数学归纳,我们可以得出
也就是说定和式和通常的有界和式是相同的
假如给定一个未知的和式,希望用封闭形式对它的值进行计算,并假设它可以写成形式∑a≤k<bg(k)=∑bag(x)δx,只要能求得一个补丁和式或者逆差分函数f,使得g(x)=f(x+1)−f(x),那么由有限微积分的理论,就可以将答案表示成f(b)−f(a),理解这个原理的一种方法是利用三点记号将∑a≤k<bg(k)完全写开
除了f(b)−f(a)以外,右边的每一项都抵消了,所以f(b)−f(a)就是该和式的值,该种和式常被称为叠缩
但该法则仅a≤b时适用,若a>b
类似方法可以论证
确定的求和法为我们提供了一种计算下降幂的和式的简单方法
定和式的方法也暗示了,对0≤k<n求和通常比对1≤k≤n求和更简单,前者是f(n)−f(0),后者是f(n+1)−f(1)
通常幂也可以用新方法来求和,如果我们先用下降幂表示出通常幂
例如
所以
用n+1替换n又给出另一种将我们的老朋友S(2)n=∑0≤k≤nk2计算成封闭形式的方法
从平方到立方
后面讨论的斯特林数可以做到通常幂和阶乘幂的转换
从而
下降幂还有性质,对应(x+y)2=x2+2xy+y2,也有(x+y)2_=x2_+2xy+y2_
实际上对(x+y)m和(x+y)m_之间也有相似的结果
看上去除以x−1会从x2_得到x1_,除以x会从x1_得到x0_,若接下来除以x+1会从x0_得到x−1_
所以x−m=1(x+1)(x+2)…(x+m),m>0
类比通常幂的指数法则
下降幂的指数法则形式是
对负的n有
我们来确定负数下降幂关于差分性质的成立
成立了
如此只要零不做除数
当m=−1时,回想积分
再回到和式,我们希望找到f(x)使得
当x是整数时
所以f(x)即调和数Hx,于是Hx就是连续的lnx的离散模拟
现在能对下降幂的和式给出完整定义了
这个公式指出了为什么快速排序这样的离散问题解中会冒出调和数,恰如再连续性问题的解中自然会出现自然对数一样
既然我们找到了一个lnx的一个类似物,我们也来研究下是否有一个ex的类似物存在,什么样的函数f(x)有与恒等式D(ex)=ex对应的性质Df(x)=f(x)?
可以取f(x)=2x作为离散指数函数
cx的差分也很简单
故如果c≠1,cx的逆差分是cxc−1,这一事实与前面的基本法则结合起来,就给出了理解几何级数和的一般公式的方法
每当遇到一个可能用作封闭形式的函数f,我们就能计算出它的差分Δf=g,然后就有一个函数g,它的不定和式∑g(x)δx是已知的
但有些连续概念并没有离散的相似概念,例如无限微积分的链式法则是对函数的函数求导的方便法则,但有限微积分没有对应的链式法则,因为Δf(g(x))没有很好的形式
除了像用c±x替换x这样的情形之外,离散的变量变换很难
然而Δ(f(x)+g(x))的确有比较好的形式,而且提供了一个分布求和的法则,它是无限微积分中称为分部积分法则的一个有限相似结果
回顾以下无限微积分中的公式
在积分并重新排列各项的次序之后,它引导出分部积分法
可以在有限微积分中做到类似的事情
这个公式可以由Ef(x)=f(x+1)
所定义的移位算子E表示成方便的形式,用Ev(x)替换v(x+1)就得到乘积的差分的一个紧凑法则:
在这个方程两边取不定和,并重新排列它的项,就得到广为宣扬的分部求和法则:
(无限微积分在此令1→0)来避开了E
如同对无限微积分那样,可以在所有三项上加上界限,从而使得不定和式变成确定的和式
当左边的和式比右边的和式更难以处理时,这个法则是有用的,例如∫xexdx是分部积分的一个典型例子,它的离散模拟是∑x2xδx
为了用分部求和,我们用u(x)=x,Δv(x)=2x,从而Δu(x)=1,v(x)=2x,且Ev(x)=2x+1
加上界限,我们可以用此式来计算以前做过的和式
用这个方法比用扰动法更容易求出这个和,因为不用动脑子
在之前我们碰巧发现了∑0≤k<nHk的一个公式,但是如果知道分部求和,那么就能遵循规则发现公式
例如有和式∑0≤k<nkHk,如果用∫xlnxdx的模拟来引导,解就很好求了
令u(x)=Hx,Δv(x)=x=x1_,从而Δu(x)=x−1_,v(x)=x2_2,Ev(x)=(x+1)2_2
第一行到第二行利用了指数法则将(x+1)2_x−1_结合了起来
现在添加界限并得到结论
无限和式#
对于无限和式,处理∑时的方法并不总是有效的,但存在一大类容易理解的无限和式,我们对它们做过的所有运算都是完全合法的,在我们更密切审视求和的基本意义之后,这两个消息背后所隐藏的理由将会变得清晰起来
无限和式可以很自然的这样定义:
等于2,因为如果将它加倍就得到
另一方面,同样的推理提示我们应该定义
是−1,因为
笑死,其实应该说成T=∞,因为T中相加的项会大于任何指定的有限数
我们尝试对一般和式∑k∈Kak的值构想一个好的定义,其中K可以是无限的
首先,让我们假设所有的项ak都是非负的,这样就不难找到一个合适的定义:如果有一个常数A为界,使得对所有有限的子集F⊂K都有
那么我们就定义∑k∈Kak是最小的这样的A(它由实数众所周知的性质得出:所有这样的A总包含一个最小的元素)
如果没有常数A为界限,我们就说∑k∈Kak=∞
上面的定义与指标集K中可能存在的次序无关,下面我们打算要讨论适用于带有多个指标k1,k2的多重和式,而不仅仅是指标集为(单一)整数集合的和式
在K是非负整数集合的特殊情形下,我们对于非负项ak的定义就意味着
理由是:实数的任何一个非减序列都有极限(可能是∞),如果极限是A,又如果F是任意一个非负整数的有限集合,其元素全都≤n,我们就又∑k∈Fak≤∑nk=0ak≤A,因为A=∞或者A是一个有界常数
又如果A′是任何一个小于所述极限A的数,那么就存在一个n,使得∑nk=0ak>A′,该有限集合f={0,1,…,n}证明A′不是有界常数
按照刚刚的定义,我们很容易地计算某种无限和式的值,例如如果ak=xk就有
特别地,刚刚考虑过的无限和S和T分别取值2和∞,恰如我们所猜测的,另一个有趣的例子是
现在考虑和式中既可能有非负项又可能有负项的情形,例如
对他进行分组
但如果延迟一项再分组
我们或许也可以尝试再公式∑k≥0xk=11−x中令x=−1,因为我们已经证明了这个公式当0≤x<1时成立,但这样就迫使我们得出结论:这个无限和是12,尽管它是由整数组成的和式
另一个有趣的例子是双向无限的∑kak,其中当k≥0时ak=1k+1,而当k<0时,ak=1k−1,我们可以把它写成
如果我们从位于中心的元素开始往外计算这个和式
如果将所有括号左移一步
因为最内层的n个括号中,所有数之和是
类似的讨论声明,如果将这些括号向左或者向右任意移动固定步,它的值都是1,这使我有勇气相信这个和式的确是1
另一方面,如果我们以下述方式对项分组
那么从内往外第n对括号包含数
在后面将证明limn→∞(H2n−Hn+1)=ln2
于是这一组合方式表面,这个双向无限和式实际上应该等于1+ln2
按照不同方式对其项相加而得出不同值的和式,有某些不同寻常之处。不同教材有不同定义
事实上,我们关于无限和式的定义相当简单,设K是任意一个集合,而ak是对每一个k∈K定义的实值项(这里k实际上可以代表若干个指标k1,k2,…,因而K可以是多维的),任何实数x都可以写成其正的部分减去负的部分
(或者x+=0,或者x−=0,或者两者皆成立),我们已经说明了怎样来定义无限和式∑k∈Ka+k和∑k∈Ka−k的值,因为a+k和a−k都是非负的,这样一来,一般性的定义是
除非右边的两个和式都等于∞,在后面这种情形,我们不定义∑k∈Kak
设A+=∑k∈Ka+k,A−=∑k∈Ka−k,如果A+和A−两者都是有限的,就说和式∑k∈Kak绝对收敛于值A=A+−A−。如果A+=∞,A−是有限的,就说和式∑k∈Kak发散于+∞,类似地,如果A−=∞而A+是有限的,就说∑k∈Kak发散于−∞,如果两者都是∞,结果还很难说
我们从非负项所做的定义开始,然后再将它推广到实值的项,如果ak是复数,我们可以再次用显然的方式将定义推广:和式∑k∈Kak定义成∑k∈Kℜak+i∑k∈Kℑak,其中,ℜak和ℑak是ak的实部和虚部,这两个和式都有定义。反之,没有定义
坏消息是某些无限和式必须不予定义,因为我们所做的操作再所有这样的情形中都可能产生矛盾,好消息是,只要我们处理的是刚才所定义的绝对收敛的和式,这一章里所有操作都完全成立
通过证明每一个变换法则都保持所有绝对收敛的和式的值不变,我们可以验证好消息
需要证明的主要结果是多重和式的基本原理:经过两个或多个指标集的绝对收敛的和式永远可以对这些指标中的任何一个首先求和
将来会正式地证明:如果J是任意的指标集,并且{Kj|j∈J}的元素是任意的指标集,使得
那么对每一个j∈J都存在复数Aj,使得
只要再所有的项都非负时证明时就够了,因为可以和前面一样把每一项分解成实部和虚部,正的部分和负的部分来证明一般情形
于是,我们假设对所有的指标对(j,k)∈M都有aj,k≥0,其中M是主指标集{(j,k)|j∈J,k∈Kj}
给定∑(j,k)∈Maj,k是有限的,即对所有有限子集F⊆M有
而A这样的最小上界,如果j是J的任意一个元素,形如∑k∈Fjaj,k的每一个和都以A为上界,其中Fj是Kj的一个有限子集。从而这些有限和式有一个最小的上界Aj≥0,且根据定义有∑k∈Kjaj,k=Aj
我们仍需要证明:对所有有限子集G∈J,A是∑j∈GAj的最小上界
假设G是J的满足∑j∈GAj=A′>A的有限子集。我们至少可以求出一个有限子集Fj⊆Kj,使得对每个满足Aj>0的j∈G均有∑k∈Fjaj,k>AA′Aj。
但是此时有∑j∈G,k∈Fjaj,k>AA′∑j∈GAj=A,这与如下事实矛盾:对所有有限子集F⊆M有∑(j,k)∈Faj,k≤A。从而对所有有限子集G∈j都有∑Aj≤A
(我觉得上面全是废话。。)
最后,设A′是小于A的任何一个实数。如果我们能找到一个有限集合G⊆J,使得∑j∈GAj>A′,证明就完成了
我们知道存在一个有限集合f⊆M,使得∑(j,k)∈Faj,k>A′,设G是这个F中那些j组成的集合,又设Fj={k|(j,k)∈F},那么就有∑j∈GAj≥∑j∈G∑k Fjaj,k=∑(j,k)∈Faj,k>A′
证明完毕
所以现在关于无限和式的一切法则证明完毕,只要它的项的绝对值组成的所有有限和式都有一个有限的界。
由于双向无限和式再我们用两种不同方法计算时给出两个不同的答案,因为它的正项1+12+13+…必定发散于∞,否则,不论我们如何对它的项进行分组都会得到同一个答案
习题#
或
a.
b.
两个不同的指标变量用了同一指标
右边还可能是
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 如何编写易于单元测试的代码
· 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】