伯努利數是如何定義到 x/(e^x-1)的冪級數展開式的係數上的?

伯努利數是伯努利在做等冪和的時候定義的,之後又是如何被定義到上述的冪級數展開式的呢


通常,大多數這類問題都是有不顯然的原因的,而且關於Bernoulli數在各種不同場合的出現是有非常大量的文獻以及討論的。

下面是一個很知名的heuristic,能夠幾乎完全給題主的問題一個解答。如果你在讀Concrete Mathematics的話跳到9.5章,如果沒在讀這是一個簡單的介紹:

考慮隨便一個什麼函數f(x),我們可以定義一個離散的差分運算元Delta

(Delta f)(n)=f(n+1)-f(n)

在離散的狀況下這幾乎就是你能定義的導數,但是不全是.

假如差分運算元有逆運算元,就只能是求和(Sigma f)(n) = sum_{k=m}^{n-1}f(k),這樣一來我們就有:

(SigmaDelta f)(n) = sum_{k=m}^{n-1}Delta f(k)=f(m+1)-f(m)+ldots+f(n)-f(n-1)=f(n)-f(m)

這個f(m)有點煩人,就當作m不變好了,就像變上限積分一樣把它當作常數.

所以,用比較瘋狂的寫法,我們可以認為Sigma Delta=1,也就是說Sigma = frac{1}{Delta}。用中文說,意思就是求和是插分的逆運算。下面我們就來換個角度看這個「逆」frac{1}{Delta}的問題:

泰勒公式告訴我們,f(n+1)=sum_{k=0}^infty frac{f^{(k)}(n)}{k!},這當然通常是錯的,但是我們只是在進行heuristic,所以並不重要。這個公式是說,如果我們把導數運算元寫成(Df)(x)=f,泰勒公式就讓我們相信,f(n+1)=sum_{k=0}^infty frac{D^k f}{k!}(n)=e^D f(n).

可是這有什麼用呢?把這個奇怪的f(n+1)=e^Df (n)帶入到插分的定義中,我們可以看到Delta f (n)= (e^D-1) f(n), 也就是說Delta = e^D-1,這也就告訴我們

frac{1}{Delta}=frac{1}{D}cdotfrac{D}{e^D-1}

看到那個熟悉的Bernoulli數的生成函數了吧?在這個背景下,題主好奇的frac{x}{e^x-1}幾乎給出了差分運算元的逆。

用Bernoulli數作為係數展開:

frac{1}{Delta}=frac{1}{D}sum_{k=0}^infty B_k frac{D^k}{k!}=frac{1}{D}+sum_{k=1}^infty B_k frac{D^{k-1}}{k!}

這個結論是什麼意思呢?我們把frac{1}{Delta} f(n)=[frac{1}{D}+sum_{k=1}^infty B_k frac{D^{k-1}}{k!}]f(n)的具體形式寫下來看:

sum_{l=m}^{n-1}f(l)=int_m^n f(x)dx+sum_{k=1}^infty B_k frac{f^{(k-1)}(n)-f^{(k-1)}(m)}{k!}

嘿!這看上去是個非常好的公式,雖然不知道對不對。

那麼我們用個例子看看這個公式對不對:

sum_{l=1}^{n-1}l=int_1^n xdx+ B_1(n-1)=frac{1}{2}(n^2-1)-frac{1}{2}(n-1)=frac{1}{2}n(n-1)

好像還不錯。

但是很遺憾,很多情況下(除了多項式,也許?),第二個無窮和不收斂,所以這個公式不太對。參見Concrete Mathematics p457。 不過,我們計算出來的形式恰好是Euler最初得到的形式. 正確的形式稱為Euler-Maclaurin求和公式:

sum_{l=m}^{n-1}f(l)=int_m^n f(x)dx+sum_{k=1}^pB_k frac{f^{(k-1)}(n)-f^{(k-1)}(m)}{k!}+R_p

那個余項R_p中有Bernoulli多項式和我們的函數的p階導數乘積的積分。這個余項是非常重要的。

關於這個公式,如下的鏈接是一個比較不錯的故事書:

Dances between continuous and discrete:
Euler』s summation formula

https://www.math.nmsu.edu/~davidp/euler2k2.pdf

Bernoulli數本身是個非常複雜的話題,它出現在(但不限於)如下的場合:

1. 解析數論里zeta函數在正偶數的值

2. von Staudt-Clausen定理以及代數數論里的Kummer"s criterion

3. 矩陣里Campbell-Baker-Hausdorff公式:log (e^A e^B)的展開,參見http://www.math.harvard.edu/~shlomo/docs/lie_algebras.pdf

4. 正切、餘切函數的泰勒展開

5. 拓撲 里的Todd class

6. 某個我不懂的東西:http://ocw.mit.edu/courses/mathematics/18-238-geometry-and-quantum-field-theory-fall-2002/lecture-notes/sec5.pdf

7. Mathoverflow上的http://mathoverflow.net/questions/61252/why-do-bernoulli-numbers-arise-everywhere/187599#187599

和這個討論https://golem.ph.utexas.edu/category/2008/02/metric_spaces.html#c014936

希望有幫助


先簡單地說一下吧~伯努利數

Bn的第一次發現是與下列數列求和的公式有關:

sum_{k=0}^{m-1}{k^n}=0^n+1^n+2^n+..........+(m-1)^n

其中n為固定的任意正整數。

這數列的求和公式中,只看最後一項,可以知道左邊的通項可以寫成:

(m-i)^n  其中,iin (1,m)

下面我們看一個生成函數

e^x-1=(1+x+frac{x^2}{2!}+frac{x^3}{3!}+.....+frac{x^n}{n!}+....... )-1=x+frac{x^2}{2!}+frac{x^3}{3!}+......

所以

我們有

x=(B_{0}+B_{1}x+B_{2}frac{x^2}{2!}+B_{3}frac{x^3}{3!}+.... +B_{n}frac{x^n}{n!}+........)	imes (x+frac{x^2}{2!}+frac{x^3}{3!}+frac{x^4}{4!}+.......+frac{x^n}{n!})

如果記a_{0}=0,a_{n}=1 (ngeq 1)

我們有

x=(sum_{n=0}^{infty}{frac{a_{n}}{n!}} x^n)	imes (sum_{n=0}^{infty}{frac{B_{n}}{n!}} x^n)

我們再定義一個項數C_{n}

x=(sum_{n=0}^{infty}{frac{c_{n}}{n!}} x^n)

再看看多項式展開,x的係數相等,所以Bn和an係數,下標和為n就可以啦~

我們有:

c_{n}=sum_{k=0}^{n}{C_{n}^{k}B_ka_{n-k} }

因為

a_{0}=0,a_{n}=1 (ngeq 1)

所以我們有:

c_{n}=sum_{k=0}^{n-1}{C_{n}^{k}B_k}

x=(sum_{n=0}^{infty}{frac{c_{n}}{n!}} x^n)兩端關於x的冪相等,所以最後只保留n=1的項。於是,我們有:

sum_{k=0}^{n-1}{C_{n}^{k}B_k}=0 n=(2,3,....)

在上面式子兩端加上B_n

我們有:

sum_{k=0}^{n}{C_{n}^{k}B_k}=B_n

這就是B_n的遞推公式。

令n=2,B0=1

B0+2B1+B2=B2

由此B1=1/2

再帶入回去:

x=(sum_{n=0}^{infty}{frac{a_{n}}{n!}} x^n)	imes (sum_{n=0}^{infty}{frac{B_{n}}{n!}} x^n)

就變成了:

frac{x}{e^x-1} +frac{x}{2}=1+(sum_{n=2}^{infty}{frac{B_{n}}{n!}} x^n)

很容易知道,左邊的函數是f(x)=f(-x)

所以,等式右邊只有關於x的偶次項。

B_{2n-1}=0

這便是用生成函數的方法來表示數~

現在生成函數的方法已經進一步拓展到其他。

數學是描述邏輯的語言~比如說,你要證明一個邏輯,把N表示成m個數的平方和,這種表示方法的個數。這種邏輯怎麼證明?用其他的工具?實驗?都無法來直接的證明。所以,數學就告訴你,我獨創了一種語言。比如,我就定義一個數論函數

這個函數就是把N表示成m個數的平方和的表法個數~於是,橢圓函數論中一個著名的等式就誕生了~那就是

(1+2x+2x^4+2x^9+....)^2=1+4(frac{x}{1-x}-frac{x^3}{1-x^3}+frac{x^5}{1-x^5}-..... )

有人說,為什麼會有左邊的公式。你仔細想想

sum_{m=-infty }^{infty} x^{m^2} =1+2x+2x^2+2x^4+2x^9+......

這是恆成立的~因為根據級數的定義嘛~

這就是解決數學問題的第一步,把常見的邏輯問題「翻譯」成為數學語言。而平方和的表法個數,我們用初等數論沒法解決,我們考慮引入一個新的東西,比如把平方弄到指數上去。當然你弄到對數上也可以,不過對數做起來不怎麼方便,這是第一步。

第二步,根據等式尋找邏輯關係。平方表示方法個數。要有m_{1}^2+m_2^2

把上面等式兩邊平方,是不是左邊x的指數就變成了m_{1}^2+m_2^2?

而右邊,令n=m_{1}^2+m_2^2

是不是就是表示把n表示成兩個平方和相加有幾種表示方法了?

數學類似這樣的問題都可以用生成函數的方法解決~


那麼, 題主應該知道,

所以, 有遞推式

[m=0]代表m=0時值為1, m不為零時值為零.

即直接展開再用遞推式化簡就可以了.答主直接複製粘貼了Concrete Mathematics和Analysis by Its History中的內容, 皆可參見, 兩書均是不錯的.


推薦閱讀:

可不可以說雖然數學歸納法本身是嚴謹的,但我們使用時卻不可能做到完全嚴謹?
中國為什麼沒有翻譯法國數學家的著作?
物理學與數學中有哪些特別有意義的名詞?
如何理解banach tarski悖論?
幾何平均數的空間意義是什麼啊?

TAG:數學 | 數論 | 高等數學 |