斐波那契数

来自testwiki
跳转到导航 跳转到搜索

Template:Refimprove Template:NoteTA

以斐波那契數為邊的正方形拼成的近似的黃金矩形(1:1.618)

-{zh-hant:費波那契數;zh-hans:斐波那契数}-意大利语:Numero di Fibonacci),又譯為菲波拿契數菲波那西數斐氏數黃金分割數、費氏數列。所形成的數列稱為-{zh-hant:費波那契數列;zh-hans:斐波那契数列}-意大利语:Successione di Fibonacci),又譯為菲波拿契數列菲波那西數列斐氏數列黃金分割數列、費氏數列。這個數列是由意大利數學家斐波那契在他的《算盤書》中提出。

數學上,斐波那契數是以遞歸的方法來定義:

  • F0=0
  • F1=1
  • Fn=Fn1+Fn2n2

用白話文來說,就是斐波那契數列由0和1開始,之後的斐波那契數就是由之前的兩數相加而得出。首幾個斐波那契數是:

Template:數列987……Template:OEIS

特別指出0 不是第一項,而是第零項(F0)。 Template:TOClimit

起源

公元1150年印度數學家Gopala金月在研究箱子包裝物件長宽剛好為1和2的可行方法數目時,首先描述這個數列。在西方,最先研究這個數列的人是比薩的李奧納多(意大利人斐波那契Leonardo Fibonacci, 1175-1250),他描述兔子生長的數目時用上了這數列:

兔子对的数量就是斐波那契数列
  • 第一個月初有一對剛誕生的兔子
  • 第二個月之後(第三個月初)牠們可以生育
  • 每月每對可生育的兔子會誕生下一對新兔子
  • 兔子永不死去

假設在n月有兔子總共a對,n+1月總共有b對。在n+2月必定總共有a+b對:因為在n+2月的時候,前一月(n+1月)的b對兔子可以存留至第n+2月(在當月屬於新誕生的兔子尚不能生育)。而新生育出的兔子對數等於所有在n月就已存在的a

斐波纳契数是杨辉三角的每一条红色对角线上数字的和。

斐波纳契数也是杨辉三角形(即帕斯卡三角形)的每一条红色对角线上数字的和。

表達式

為求得斐波那契數列的一般表達式,可以藉助線性代數的方法。高中的初等數學知識也能求出。

初等代數解法

已知:

  • a1=1
  • a2=1
  • an=an1+an2(n≥3)

首先構建等比數列

an+αan1=β(an1+αan2)
化簡得
an=(βα)an1+αβan2
比較係數可得:
{βα=1αβ=1
不妨設β>0,α>0
解得:

{α=512β=5+12
又因为有an+αan1=β(an1+αan2), 即{an+αan1}為等比數列。

求出數列{an+αan1}

由以上可得:
an+1+αan=(a2+αa1)βn1=(1+α)βn1=βn

變形得: an+1βn+1+αβanβn=1β。 令bn=anβn

求數列{bn}進而得到{an}

bn+1+αβbn=1β
bn+1+λ=αβ(bn+λ),解得λ=1α+β。 故數列{bn+λ}為等比數列
bn+λ=(αβ)n1(b1+λ)。而b1=a1β=1β, 故有bn+λ=(αβ)n1(1β+λ)
又有{α=512β=5+12bn=anβn
可得an=55[(1+52)n(152)n]

得出an表達式

an=55[(1+52)n(152)n]

用數學歸納法證明表達式

證明Fn=15[φn(1φ)n],其中φ黃金比例1+52n為任意整數
  • n為非負整數
n=0時,15[φ0(1φ)0]=15[11]=0=F0,成立
n=1時,15[φ1(1φ)1]=15[φ1+φ]=15[2φ1]=15×5=1=F1,成立
設當n=kn=k+1時皆成立,即Fk=15[φk(1φ)k]Fk+1=15[φk+1(1φ)k+1]
n=k+2
Fk+2=Fk+1+Fk=15[φk+1(1φ)k+1]+15[φk(1φ)k]=15[φk+1+φk(1φ)k+1(1φ)k]=15{φk(φ+1)(1φ)k[(1φ)+1]}=15{φk(φ2)(1φ)k[(1φ)2]}=15{φk+2(1φ)k+2}
亦成立
  • n為非正整數
n=0時,成立
n=1時,15[φ1(1φ)1]=15[(φ1)(φ)]=15[2φ1]=15×5=1=F1,成立
設當n=kn=k1時皆成立,即Fk=15[φk(1φ)k]Fk1=15[φk1(1φ)k1]
n=k2
Fk2=FkFk1=15[φk(1φ)k]15[φk1(1φ)k1]=15[φkφk1(1φ)k+(1φ)k1]=15{φk1(φ1)(1φ)k1[(1φ)1]}=15{φk1(φ1)(1φ)k1[(1φ)1]}=15{φk2(1φ)k2}
亦成立

因此,根據數學歸納法原理,此表達式對於任意整數n皆成立

線性代數解法

(Fn+2Fn+1)=(1110)(Fn+1Fn)

(Fn+2Fn+1Fn+1Fn)=(1110)n+1

構建一個矩陣方程

Jn為第n個月有生育能力的兔子數量,An為這一月份的兔子數量。

(Jn+1An+1)=(0111)(JnAn),

上式表達了兩個月之間,兔子數目之間的關係。而要求的是,An+1的表達式。

求矩陣的特徵值λ

根据特征值的计算公式,我们需要算出来 |λ111λ|=0 所对应的解。

展开行列式有:λ(1λ)1×1=λ2λ1

故當行列式的值為 0,解得 λ1=12(1+5)λ2=12(15)

將兩個特徵值代入

((0111)λE)x=0


求特徵向量x

x1=(112(1+5))

x2=(112(15))

分解首向量

第一個月的情況是兔子一對,新生0對。

(J1A1)=(01)

將它分解為用特徵向量表示。

(01)=15(112(1+5))15(112(15)) (4)

(Jn+1An+1)=(0111)(JnAn)=λ(JnAn)

可得到

(Jn+1An+1)=(0111)n(J1A1)=λn(J1A1) (5)

化簡矩陣方程

將(4) 代入 (5)

(Jn+1An+1)=λn[15(112(1+5))15(112(15))]

根據3

(Jn+1An+1)=15λ1n(112(1+5))15λ2n(112(15))

求A的表達式

現在在6的基礎上,可以很快求出An+1的表達式,將兩個特徵值代入6中

An+1=15λ1n+115λ2n+1
An+1=15(λ1n+1λ2n+1)
An+1=15{[12(1+5)]n+1[12(15)]n+1}(7)

(7)即為An+1的表達式

數論解法

實際上,如果將斐波那契數列的通項公式寫成anan1an2=0,即可利用解二階線性齊次遞迴關係式的方法,寫出其特徵多項式λ2λ1=0(該式和表達斐波那契數列的矩陣的特徵多項式一致),然後解出λ1=12(1+5)λ2=12(15),即有an=c1λ1n+c2λ2n,其中c1,c2为常数。我们知道a0=0,a1=1,因此{c1+c2=0c1(1+5)2+c2(15)2=1,解得c1=15,c2=15

組合數解法

Fn=i=0(nii)[1]

Fn1+Fn=i=0(n1ii)+i=0(nii)=1+i=1(nii1)+i=1(nii)=1+i=1(n+1ii)=i=0(n+1ii)=Fn+1

黃金比例恆等式解法

φ黃金比例1+52,則有恆等式φn=Fn1+φFn(1φ)n=Fn+1φFn,其中n為任意整數Template:NoteTag,則

φn(1φ)n=(Fn1+φFn)(Fn+1φFn)=(Fn1Fn+1)+2φFn=Fn+2φFn=Fn(2φ1)=Fn×5

因此得到Fn的一般式:

Fn=15[φn(1φ)n]=15[(1+52)n(152)n]

此一般式對任意整數n成立

近似值

n為足夠大的正整數時,则

Fn15φn=15[12(1+5)]n0.44721359551.61803398875n
Fn15(1φ)n=15[12(15)]n0.4472135955(0.61803398875)n

用計算機求解

可通過編程觀察斐波那契數列。分為兩類問題,一種已知數列中的某一項,求序數。第二種是已知序數,求該項的值。

可通過遞歸遞推的算法解決此兩個問題。 事實上當n相當巨大的時候,O(n)的遞推/遞歸非常慢……這時候要用到矩陣快速幂這一技巧,可以使遞迴加速到O(logn)。

和黃金分割的關係

開普勒發現數列前、後兩項之比12,23,35,58,813,1321,2134,,也組成了一個數列,會趨近黃金分割

fn+1fna=12(1+5)=φ1.618...

斐波那契數亦可以用連分數來表示:

11=121=1+1132=1+11+1153=1+11+11+1185=1+11+11+11+11

Fn=15[(1+52)n(152)n]=φn5(1φ)n5

而黃金分割數亦可以用無限連分數表示:

φ=1+11+11+11+11+...

而黃金分割數也可以用無限多重根號表示:

φ=1+1+1+1+...

和自然的關係

Template:Further Template:See also

春黄菊頭狀花序上,小花呈螺旋狀排列,從不同方向可以數出21(深藍)和13(淺藍)條旋臂,為相鄰的斐氏數。類似的螺旋狀排列見於多種植物。

斐氏數列見於不同的生物學現象[2],如樹的分枝、Template:Link-en、菠蘿聚花果上小單果的排列、[3]雅枝竹的花蕾、正在舒展的蕨葉、松毬的鱗的排列[4]、蜜蜂的家族樹[5][6]开普勒曾指出斐氏數列存在於自然界,並以此解釋某些花的五邊形形態(與黄金分割率相關)。Template:Sfn法國菊的「瓣」(舌狀花)數通常為斐氏數。Template:Sfn1830年,K. F. Schimper和A. Braun發現植物的旋生葉序中,連續兩塊葉之間轉過的角度與周角之比,約成整數比時,常出現斐氏數[7],如2/55/13[8]

恆等式

資料來源:[9]

證明以下的恆等式有很多方法。以下會用組合論述來證明。

  • Fn可以表示用多個1和多個2相加令其和等於n的方法的數目。

不失一般性,我們假設n1Fn+1是計算了將1和2加到n的方法的數目。若第一個被加數是1,有Fn種方法來完成對n1的計算;若第一個被加數是2,有Fn1來完成對n2的計算。因此,共有Fn+Fn1種方法來計算n的值。

  • F0+F1+F2+F3+...+Fn=Fn+21

計算用多個1和多個2相加令其和等於n+1的方法的數目,同時至少一個加數是2的情況。

如前所述,當n>0,有Fn+2種這樣的方法。因為當中只有一種方法不用使用2,就即1+1+...+1 (n+1項),於是我們從Fn+2減去1。

  1. 若第1個被加數是2,有Fn種方法來計算加至n1的方法的數目;
  2. 若第2個被加數是2、第1個被加數是1,有Fn1種方法來計算加至n2的方法的數目。
  3. 重複以上動作。
  4. 若第n+1個被加數為2,它之前的被加數均為1,就有F0種方法來計算加至0的數目。

若該數式包含2為被加數,2的首次出現位置必然在第1和n+1的被加數之間。2在不同位置的情況都考慮到後,得出Fn+Fn1+...+F0為要求的數目。

  • F1+2F2+3F3+...+nFn=nFn+2Fn+3+2
  • F1+F3+F5+...+F2n1=F2n
  • F2+F4+F6+...+F2n=F2n+11
  • F12+F22+F32+...+Fn2=FnFn+1
  • FnFmkFmFnk=(1)nkFmnFk,其中m,n,kF的序數皆不限於正整數。Template:NoteTag
    • 特別地,當n=mk時,Fn2Fn+kFnk=(1)nkFk2
      • 更特別地,當k=1k=1時,對於數列連續三項,有Fn2Fn1Fn+1=(1)n1
    • 另一方面,當(m,n,k)=(n+1,n,2)時,對於數列連續四項,有FnFn+3Fn+1Fn+2=(1)n+1Template:NoteTag
  • φn=Fn1+φFn(1φ)n=Fn+1φFn,其中φ黃金比例1+52n為任意整數Template:NoteTag
藉由上述公式,又可推得以下恆等式Template:NoteTag
    • FmFn+Fm1Fn1=Fm+n1[9]
    • FmFn+1+Fm1Fn=Fm+n[9]

      特別地,當m=n時,F2n1=Fn2+Fn12F2n=(Fn1+Fn+1)Fn=(2Fn1+Fn)Fn

數論性質

公因數和整除關係

  • Fn整除Fm,若且唯若n整除m,其中n3
  • gcd(Fm,Fn)=Fgcd(m,n)
  • 任意連續三個菲波那契數兩兩互質,亦即,對於每一個n
gcd(Fn,Fn+1)=gcd(Fn,Fn+2)=gcd(Fn+1,Fn+2)=1

斐波那契质数

Template:Main 在斐波那契數列中,有質數[10] 2, 3, 5, 13, 89, 233, 1597, 28657, 514229, 433494437, 2971215073, 99194853094755497, 1066340417491710595814572169, 19134702400093278081449423917……

截至2015年,已知最大的斐波那契質數是第104911個斐波那契數,一共有21925個十進制位。不过,人们仍不知道是不是有无限个斐波那契质数。[11]

Template:Section link所述,Fkn總能被Fn整除,故除F4=3之外,任何斐氏質數的下標必同為質數。由於存在Template:Link-en的一列連續合数,斐氏數列中亦能找到連續任意多項全為合數。

大於F6=8的斐氏數,必不等於質數加一或減一。[12]

與其他數列的交集

斐波那契数列中,只有3個平方數01144[13][14]2001年,Template:Link-hu證明,斐氏數中的次方數衹有有限多個。[15]2006年,Y. Bugeaud、M. Mignotte、S. Siksek三人證明,斐波那契数中的次方數只有0、1、8、144。[16]

1、3、21、55為僅有的斐氏三角形數Template:Link-en曾猜想此結論,後來由罗明證明。[17]

斐波那契數不能為完全数[18]推而廣之,除1之外,其他斐氏數皆非多重完全數[19],任兩個斐氏數之比亦不能是完全數[20]

n的週期性

Template:Main

斐波那契數列各項模n的餘數構成Template:Le,其最小正週期稱為皮萨诺周期[21],至多為6n[22]。皮薩諾週期對不同n值的通項公式仍是未解問題,其中一步需要求出某個整數(同餘意義下)或二次有限域元素的Template:Link-en。不過,對固定的n,求解模n的皮薩諾週期是Template:Link-en問題的特例。

推廣

斐波那西數列是斐波那西n步數列步數為2的特殊情況,也和盧卡斯數列有關。

和盧卡斯數列的關係

FnLn=F2n

反費波那西數列

反費波那西數列的遞歸公式如下:

Gn+2=GnGn+1

如果它以1,-1開始,之後的數是:1,-1,2,-3,5,-8, ...

即是F2n+1=G2n+1=F(2n+1),F2n=G2n=F2n

亦可寫成Fm=(1)m+1Gm=(1)m+1Fm,其中m是非負整數。

反費波那西數列兩項之間的比會趨近1φ0.618

證明關係式

證明Fm=(1)m+1Fm,其中m是非負整數

φ表示黃金分割數1+52,則有φ(1φ)=1
(1)m=[φ(1φ)]m=φm(1φ)m,因此
(1)m+1Fm=(1)m+1×15[φm(1φ)m]=(1)×(1)m×15[φm(1φ)m]=(1)×φm(1φ)m×15[φm(1φ)m]=(1)×15[φm+m(1φ)m(1φ)m+mφm]=(1)×15[(1φ)mφm]=15[φm(1φ)m]=Fm

巴都萬數列

費波那西數列可以用一個接一個的正方形來表現,巴都萬數列則是用一個接一個的等邊三角形來表現,它有Pn=Pn2+Pn3的關係。

佩爾數列

佩爾數列的遞歸公式為Pn=2Pn1+Pn2,前幾項為0,1,2,5,12,29,70,169,408,...

應用

1970年,尤裏·馬季亞謝維奇指出了偶角標的斐波那契函數

y=F2x

正是滿足Julia Robison假設的丟番圖函數,因而證明了希爾伯特第十問題是不可解的。

電腦科學

高為6的斐波那契樹。平衡因子以綠色標記,節點的高度則為紅色。

最左一條路徑上的鍵值全為斐氏數。

延伸閱讀

  • KNUTH, D. E. 1997. The Art of Computer ProgrammingArt of Computer Programming, Volume 1: Fundamental Algorithms, Third Edition. Addison-Wesley. Chapter 1.2.8.
  • Arakelian, Hrant (2014). Mathematics and History of the Golden Section. Logos, 404 p. ISBN 978-5-98704-663-0, (rus.)
  • 克裏福德A皮科夫.數學之戀.湖南科技出版社.

參考文獻

Template:Reflist

註釋

Template:NoteFoot

參見

外部連結

Template:級數 Template:貴金屬比例 Template:Authority control