在《數學傳播》2020 年第1 期刊登的的探源與推廣 一文[19] 中, 張進安老師推導了與Fibonacci 數列(按定義有關的幾個無窮級數的運算式, 例如他得到了標題中的等式 (改用無窮級數的常用記號)
張進安老師還得到, 對一切正整數 k, 有
他認為(2)式中的k不必限於正整數, 例如取 k=log2, 就有
更進一步, 張進安老師認為可以令 k=logr(注意, 這相當於令), 這裡r是一個正數, 從而(2)式成為
[在原文中,(4)式右邊分母中的r被誤寫為。]
張進安老師在文章末尾進一步提出一個問題。他注意到, (4) 中的 r 不必限於正數, 對 r=−2,−3也成立。於是問:有誰能將這個問題說得更清楚?臺灣大學數學系張鎮華教授在發表於2020 年第2 期《數學傳播》的文章中給出了一個清楚的回答, 我們在此給出另一個解決方案並推而廣之。
整個文章所討論的問題, 似乎以生成函數(generating function) 的觀點最為直接。令Fibonacci 數列的生成函數為
則一個簡單的推理告訴我們, 如果對某個複數z收斂, 則其無窮和G(z)為
推理如下 :
從而
是方程的兩根。現在, 通過做變數替換, 張進安老師的問題轉化為:
問題: 冪級數(5)究竟對那些複數z收斂呢? (也許張進安老師只考慮實數。)
回答這類問題的一個基本結果, 是著名的Cauchy--Hadamard 公式。我們引述如下(參見[2][p.73] 定理1.2, 作者稱之為Abel 定理, 並將Cauchy--Hadamard 公式稱作Hadamard 公式):
定理1 (Cauchy--Hadamard公式):冪級數在以R為半徑的圓盤 |z|<R內收斂,其中R叫做斂半徑,按照下述公式確定:
且冪級數在這個圓的外部任何點處都發散。
回到的具體情況, 看來我們是要計算這樣一個麻煩的上極限
如果我們知道的通項公式, 其實是可以求出上述極限(準確的說, 可以直接算出極限)。當然如果知道的通項公式, 張老師的問題完全可以直接解決(我們留給有興趣的讀者。) 如果不利用的通項公式, 可否求出這個收斂半徑R呢?回答是肯定的。
在複變函數中, 我們有下述基本結果(參見 [3][pp. 73--74], 作者概括為函數f(z)在其冪級數收斂圓周上必有奇點, 我們重新表述如下) :
定理2 (Riemann定理):函數f(z)在原點處的冪級數的收斂圓盤, 恰好是使得函數f(z)有定義並且解析的最大圓盤 |z|<R。
據 [6] 講, Riemann 在1856 年關於複分析的講義中提出定理2, 並(早於Hadamard) 重新發現了定理1中的公式(9)。
回到冪級數(5), 我們知道, 它是在原點的冪級數,於是根據定理2,(5)的收斂圓盤是使得有定義且解析的最大圓盤|z|<R。因為僅僅在在處沒有定義且不解析, 從而容易看出
於是根據定理1," 級數(5)對對一切收斂,而對一切發散。進一步, 利用的通項公式(參見 [20][p.61]) 可以確定, 冪級數(5)在收斂圓周|z|=上均發散。綜合起來, 我們得到
收斂⟺|z|<,且此時有=. (11)於是在(11)中令z=, (r≠0), 就得到收斂⟺|r|>,且此時有=.(12)特別地, 在(12)中分別令 r=10,(k為正整數),2 就得到 (1)(2)(3)各式。而(12)所框定的r之取值範圍 |r|>(也許他只是考慮實數值的r) 就是張老師最後的問題的答案。
注1:這裡介紹的方法適用於由高階遞推關係(給定初值,)
(其中是常數) 給出的數列所對應的冪級數
的收斂區域問題之討論。此時容易確定, 若(14)對某複數z收斂, 則必收斂於其生成函數(推導如前, 也可參見 [1][pp.337--338])
而且可以確定出使得(16)成立的z的一個範圍為|z|<R, (17)
其中R可以如下確定。設d(x)是P(x),Q(x)的最大公因式,
則R=多項式各個根模長之最小者。 (18)
級數(14)對滿足的z發散, 對|z|=R上的點則需具體分析[我們在附記中以更簡明的方式證明了, 級數(14)的收斂域恰好是|z|<R
我們相信, 這一結果不是新的, 但絕不是人人都了解。值得注意的是, 以上結果包含了 是一個 d−1(d≥1)次多項式的特殊情形, 此時結果更簡單。考慮到這一問題張進安老師曾在另一篇文章[18]提及, 我們特別介紹一下。有下述結果:
定理3:設f(x)是x的d次多項式,則收斂⟺|z|<1,且此時有
其中是f在0處的k階差分。特別地, 當z=時,有
考慮到這個情形特別簡單, 我們給出一個直接的證明。
證明:根據級數收斂的必要條件,我們推出,當對某z收斂時有
由此容易推出
從而我們有
對冪級數等式
在收斂圓|z|<1|內求k階導, 我們有
變形即得
代入(22)式即得(19)。在(19)式中令, 即得(20),證畢。
注2:我們特別要指明,定理2可以使我們冪級數(23) 的收斂範圍是 |z|<1獲得更深刻的理解。事實上, 從複變函數的觀點來看, 根本的原因就在於複變函數對z=1恰好無定義且不解析, 從而其冪級數收斂半徑R=1。講到這裡, 我要與讀者分享自學成才的數學家IM Gelfand (1913∼2009) 在16歲期間自學數學分析的一段經歷 [16]:
我開始參加大學的研討班, 在那裡我發現自己壓力很大。我做數學的方式不合適宜。當時, 數學界掀起了一股風氣:對嚴謹的證明的要求、 對實變函數論的濃厚興趣。(今天看來, 這種嚴格性和這個特殊理論已陳腐過時, 但在那時⋯⋯)
直到那時, 我才認識到, 很重要的是:函數未必是連續的, 連續的函數未必是可微的, 一階可微的函數未必是二階可微的, 如此等等;甚至一個無窮次可微的函數, 其Taylor 級數也未必是收斂的;即便收斂, 也未必收斂到函數本身。如果函數的Taylor 級數剛巧收斂到它本身, 這個函數就稱為解析的。(實變函數論愛好者認為) 這類函數是如此狹窄, 以至於它被排除在主流數學之外。而在此之前, 我就只見過這類函數。
在這種觀點的影響下,我讀了de la Vallée-Poussin的"現代化的、嚴格化的"分析教材( Cours d'analyse )。它類似於目前莫斯科大學數學力學系用的教材,但更好一些。因此我很同情那些大一學生,他們只有在歷經長達一年的強調"嚴格基礎"的痛苦考驗之後,才能體會到數學分析的美妙。
即便如此, 我也是幸運的, 我讀了II Privalov 關於單複變函數的卓越教材(指1927 年出版的《複變函數引論》), 有閔嗣鶴等譯的中譯本)。讀這本書時, 我理解了, 為什麼函數的Taylor 級數在x=1發散, 雖然它的圖像是連續的。(這根源於下述事實:它所對應的複函數在z=i有奇性。)
讀完前100 頁, 我感到一陣清風拂過。我發現, 如果一個複變函數有一階導數, 那麼就自動有任意階的導數, 並且其Taylor 級數在某個區域內收斂到函數本身。每樣東西都找到了自己的位置, 又恢復了和諧。
張老師 [18],[19]提出的這類問題, 恰好可以利用複變函數的結果(定理2)獲得近乎圓滿的解決。已故享譽全球之大數學家陳省身先生(1911∼∼2004) 多次強調複數的美妙與重要。比如, 他曾在美國數學會的Notices 訪談中[5] 提到:"My main idea is that you should do topology or global geometry in the complex case. The complex case has more structure and is in many ways simpler than the real case. So I introduced the complex Chern classes." 相信讀者已從本文之討論對複數之美妙有所體會。
注3:最後,我要指出,在上文寫成後,為給讀者指引一個關於討論Fibonacci數列通項公式的文獻,我從《數學傳播》通過關鍵字搜索,得到意外收穫,即發現本刊刊登了幾篇與張老師所論問題相關的文獻,見[8],[9],[14], 特別是 [14]。這些作者的討論正是基於Fibonacci 數列的通項公式, 而且只考慮實數。其結果均可用此處的概念性方法而非計算得到。本文的重點不在於重新推導, 而在於從複變函數的眼光來理解其本質。正如德國大數學家Kronecker (1823∼∼1891)所說:"Analysis does not owe its really significant successes of the last century to any mysterious use of, but to the quite natural circumstances that one has infinitely more freedom of mathematical movement if he lets quantities vary in a plane instead of only on a line."
附記
本文投稿之後, 作者進一步思考了級數(14)在|z|=R上的收斂性問題。蒙昔日恩師指點, 解決了這一問題。結論是:在|z|=R上, 級數 (14)總是發散的。我們表述成以下定理(相信它不是新的)。
定理4:設多項式P(z),Q(z)互素, P(z)的次數小於Q(z)的次數, 且Q(0)≠0。設有理函數在 z=0處有冪級數展開
令人驚訝的是, 以下的證明完全避開了冪級數的收斂半徑(Abel定理), 也不需要複變函數的Riemann 定理, 只需要一個巧妙的引理(下面的引理1) 和關於有理生成函數的一般展開定理(我們有推導)。
其中是次數為的多項式(這一結果在 [1][p.341] 中被稱為有理生成函數的一般展開定理, 作者在那裡並未給出詳細證明)。因此, 我們只要討論冪級數,其中是次數為的多項式(40)收斂範圍。(38)式表明, 當 |z|<r時, 級數(40)收斂到。下面將證明,當|z|≥r時級數 (40)發散。為此,只要證明級數(40)的通項(39)不趨於零。這是下述引理的直接推論。
引理1,設是不同的複數, 模長均大於等於1, 是復係數多項式, 令,n=0,1,…,若, 則是零數列。
再附記
作者在2020 年8 月8 日了解到, 有理函數之冪級數在收斂圓周上不收斂(定理4) 這一結果確實不是新的, 它蘊含於Pólya 和Szegö 的名著《分析中的問題與定理》[15][p.152] 的問題246:
命題1: 設冪級數在收斂圓周上有一個極點,則冪級數在收斂圓周的每一點都不收斂。
其證明用到了E. Cesàro 的一個結果 ([15][p.20] 問題85), 如下:
命題2:設數列{},{}滿足以下條件:
最後, 作者還想補充一個從Fibonacci 數列可能衍出的問題(最初將我引向這個問題的, 是河南大學的陳敏茹博士與學友楊凡)。眾所周知, Fibonacci 數列FnFn 前後兩項之比的極限恰好為黃金分割數(參見[13]), 即我們有:
所以, 我們也許會問:
命題2:對由一般的常係數線性遞推關係(13)所定義的數列{} , 其前後兩項之比的極限是否存在;如果存在, 如何確定?
林鳳美老師在發表於《數學傳播》2019年第4期的文章 [10] 中考慮了這一問題, 但遺憾的是, 其結果有瑕疵(原因在於, [10][p.103] 倒數第二行公式中的係數可能等於0)。例如, 根據該頁定理3, 只要數列{}滿足(n≥2), 就有
事實上, 問題2已蘊含於Pólya--Szegö[15][p.152] 的問題242, 而且那裡的結果更一般:
命題3:設冪級數在收斂圓周上有唯一的奇點, 而且它還是極點, 則有
在常係數線性遞推關系所定義的數列的特殊情形, 命題3可以表述為以下形式:
定理5:設多項式P(z), Q(z)互素, P(z)的次數小於Q(z)的次數, 且 Q(0)≠0。設有理函數在z=0處有冪級數展開
若Q(z)的模長最小的根是唯一的, 記為, 則(49)成立。
實際上, 無需藉助命題3, 從 (39)(39) 式不難推出定理5, 我們留給有興趣的讀者。
注意, 若 Q(z)最小模長的根不止一個, 那麼 (49)不一定成立。例如, 任取非零複數a, 對由
1,a,1,a,… (50)
構成的數列, 容易算出,其生成函數為.
注意, 當 a≠±1時, 這個有理函數的分子與分母互素, 而分母的兩個根模長相等。由前後兩項之比構成的數列為
它存在極限若且唯若, 即 a=±1。因此, 若 a≠0,±1, 則數列 (50)的前後兩項之比就無極限。
從另一種角度看, 等式 (49)意味著, 這給出了求多項式 Q(z)的最小模長零點(假定它唯一) 的一個近似方法, 這本質上就是線性代數中求矩陣的佔優特徵值與特徵向量(dominant eigenvalue and eigenvector) 的冪法(power method), 參見 [17] 6.2 節或 [4][p.81] 1.4 節問題7。
致謝
作者感謝高雄市中正高中張進安老師和《數學傳播》諸位編輯老師的鼓勵推動!感謝復旦大學邵美悅博士指正初稿的一處誤拼並指引文獻[17], 並建議我們改用生成函數而非母函數的稱謂。感謝審稿人對初稿提出寶貴建議。
參考文獻
[1] Ronald Graham, Donald Knuth, Oren Patashnik, Concrete Mathematics , (2nd ed.). Reading, MA: Addison-Wesley. 1994.中譯本《具體數學》,張明堯、張凡譯。人民郵電出版社, 2013。
[2] 龔升、張德健。 複分析五講--第一講。數學傳播季刊, 34(2), 52--75, 2010。
[3] 龔升、張德健。 複分析五講--第三講。數學傳播季刊, 34(4), 50--76, 2010。
[4] Roger A. Horn and Charles R. Johnson, Matrix Analysis(second edition), Cambridge University Press, 2013. 有中譯本:《矩陣分析》, 楊奇譯。機械工業出版社, 2005年;《矩陣分析》, 張明堯、 張凡譯。機械工業出版社, 2014年。
[5] Allyn Jackson, Interview with Shiing Shen Chern , Notices of the AMS , 45(7), 860--865, 1998.
[6] D. Laugwitz, E. Neuenschwander, Riemann and the Cauchy--Hadamard formula for the convergence of power series. Historia Mathematica , 21(1), 64--70, 1994.
[7] 李克大、 李尹裕。《有趣的差分方程》(第2版)。合肥:中國科學技術大學出版社, 2019 年。
[8] 林炳炎。奇妙的費氏數列之一。數學傳播季刊, 6(2), 34--38, 1982。
[9]林炳炎。有關費氏數之無窮級數的分數和。數學傳播季刊, 7(1), 27--33, 1983。
[10] 林鳳美。高階線性遞迴數列的一般化費氏螺線。數學傳播季刊, 43(4), 99--109, 2019。
[11]林開亮。微積分之前奏(或變奏) :高階等差數列的求和。數學傳播季刊, 41(1), 61--79, 2017。
[12] 林開亮。 解常係數線性微分方程和遞推關係的新方法---秦九韶和亥維賽的遺產。數學傳播季刊, 43 (2), 63--79, 2019。
[13] 林開亮。從鹹道老師夢得的數學題談起。好玩的數學, 2020年3月24日。
[14] 林智勇、易正明、許天維。以費氏數列表示的無窮級數和與收斂半徑。數學傳播季刊, 30(3), 53--65, 2006。
[15] George Pólya and Gabor Szegö, Problems and Theorems in Analysis I, Springer, 1978.有中譯本, 《數學分析中的問題和定理》, 張奠宙、宋國棟等譯。上海:上海科學技術出版社, 1981
[16] VS Retakh and AB Sosinsky, A talk with IM Gelfand , Quantum , Jan-Feb, 1989.有中譯文,數學譯林,第4期,李錕譯, 340--347, 1990.
[17] 徐樹方、高立、張平文。《數值線性代數》(第二版)。北京: 北京大學出版社, 2013。
[18] 張進安。2n2n在分母的級數收斂性質。數學傳播季刊, 43(3), 56--59, 2019。
[19] 張進安。limn→∞∑ni=0Fi10i=1089limn→∞∑i=0nFi10i=1089的探源與推廣 。數學傳播季刊, 44(1), 89--93, 2020。
[20] 張鎮華。費氏數列與等比數列的交會處。數學傳播季刊, 44(2), 58--61, 2020。