在一条直线上,一个人从原点出发,每一次可以往左或者往右走一步.每次往左走的概率为p,往右走的概率为q.设他一共走了总步数N, 设最终所处的位置为m,设p=q,求m,m2,m3,m4.
任何物理量的平均值都是这样算的:X=∑iP(X=xi)xi.在行走问题里概率该怎么算?分步乘法、排列组合。如果往左走了n1步,往右走了n2步,那么平均值就是m=∑(n1−n2)⋅CNn1pn1qn2.因为m和n1或n2的取值是一一对应的,所以原则上我们已经"做完"了。
但是,这不是一个很informative的形式.我们是否能化简?
往左走和往右走是对称的。因此:m=n1−n2,考虑其中一个,n1=∑n1=0Nn1⋅CNn1pn1qN−n1.你可能会注意到,p=q=21,带入得到
n1=n1=0∑Nn1⋅CNn1(21)N
但很可惜,对这个问题,先带入数值会给进一步化简带来困难。相反,我们引入“量纲修改”的技巧:常数能进行加减乘除等初等运算。但是,函数还可以求导、积分。只要自变量取特定值,函数值就是一个普通的常数。所以,我们可以先把情形推广,用函数来运算,然后取特定自变量值得到常数下的结果。
我们把p看成自变量。如果你足够幸运,就可以联想到:∑n1=0NCNn1pn1qN−n1=(p+q)N. 对比f(p)=∑n1=0Nn1⋅CNn1pn1qN−n1只需要改成
f(p)=p∂p∂n1=0∑NCNn1p1nqN−n1
所以当p+q=1,
f(p)=Npq
特别的,对于p=q,n1=41N,进而m=0.
我们如法炮制,求解m2:f(p)=∑n1=0Nn12⋅CNn1pn1qN−n1.如何消去n1的平方项呢?
注意到二次求导(∂p∂)2pn1=n12pn1−2−n1p1n1−2,我们可以改写f(p)为
f(p)=n1=0∑NCNn1(n1pn1+p2(∂p∂)2pn1)qN−n1
注意到括号里的第一项简单给出n1,第二项给出二项式的二阶导数。 带入走路问题
n12=41N2+41N
换一种思路,我们把求导看成算符。算符可以移动位置,于是
f(p)=(p∂p∂)2n1=0∑NCNn1pn1qN−n1
上面的计算中,我们是这样打括号的:(p∂p∂),而不是这样:p(∂p∂).你可能回想,这有什么区别吗? 在我们的推导中,p∂p∂的意思是"对pn1求导,然后乘上p"。我们把这个整体的操作移到了求和号外面。显然,这里求导和乘法的顺序是不能交换的。如果我们多次执行"求导,然后乘上p",就应该写成
(p∂p∂)(p∂p∂)(p∂p∂)⋯
很显然,p(∂p∂)很容易诱导我们写出错误的表达式:
p⋅p⋅p⋯(∂p∂)⋅(∂p∂)⋅(∂p∂)⋯
对于三次方和四次方,做法是类似的。
在求解位置均值时我们使用了"量纲"修改的技巧
常数→自变量→常数
这种变换思想应用广泛。在热统中,我们把离散量看成是连续的,进而可以使用求导运算,其合法性由热力学极限保证;在电动力学中,我们把向量改写成指标,进而可以使用指标运算规则同时操作多个分量。 进一步推广,我们可以把直角坐标改写成球坐标来适配球形边界,可以把矩阵看成张量进行运算...每当卡壳,"量纲修改"给我们提供了一种不断变通、寻找出路方向。
在求解高次均值时我们写出了
(∂p∂)2pn1=n12pn1−2−n1p1n1−2
这种环环相扣的计算思路出现在各种高次方求和问题中。例如,求解∑i=1Nik,k∈Z+。对k=2,考虑到
(i+1)3=i3+3i2+3i+1
移项求和
(N+1)3−1=3i=1∑Ni2+3i=1∑Ni+N
于是可以求∑i=1Ni2.