高階線性遞迴數列的一般化費氏螺線
林鳳美
考慮數列 {an}, ≥ 0, 若滿足二階齊次線性遞迴關係式:
an = an−1+ an−2, (1) 則稱為費氏數列 (Fibonacci sequence)。 費氏數列具有許多與黃金比例 (golden ratio) 及費 氏螺線 (Fibonacci spiral) 密切的性質, 皆涉及美學與生物的生長模式。
事實上, 我們利用費氏數列的各項為邊作正方形, 稱為數列方形 (Sequence square), 再 依序以逆時針螺旋排列, 由 0 點出發, 不斷在正方形內逆時針作出四分之一的圓弧, 連結成一 條螺線, 稱為費氏螺線, 這近似於鸚鵡螺的螺線, 參見圖 1。
圖1: 費氏螺線與鸚鵡螺的螺線
圖 1 中的數列方形長與寬的比為費氏數列的後前項的比 an+1/an, 當項數 n 越大時, 其 值正巧為黃金比列 φ, 近似於 1.618, 或特徵正實根。
我們把 (1) 式推廣到更一般 k 階正實數齊次線性遞迴關係式:
an = c1an−1+ c2an−2+· · · + ckan−k, (2) 其中係數 c1, c2, . . . , ck 為正實數且 ck = 0。 顯然, 如果給定首 k 項的初期值, 就可確定 (2) 式唯一的整個數列, 則稱為 k 階正實數齊次線性遞迴數列。
我們給定首 k 項的初期值為 a1 = 1, ai = 0, 其中 i = 0, −1, −2, . . . , −(k − 2), 是為 了結合數列及螺線與幾何間相互間妙趣, 並且盡量以推廣與類推的眼光來開展。 美妙地, 這個 k
99
階問題的結果是豐碩的。
底下是我們要探討的主要三個問題:
(i) 唯一性問題: 數列的特徵正實根是否等於數列的後前項比的極限呢?
(費氏數列的後前項比的極限等於特徵正實根, 也是黃金比列 φ。) (ii) 幾何建構問題 : 依照費氏螺線的樣式如何建構出 k 階的數列螺線呢?
(iii) 收斂性問題 : 建構出 k 階的數列螺線收斂條件為何?
1. 唯一性問題
甲、 唯一特徵正實根
首先要求 (2) 式的一般項 an, 可嘗試用 an = xn 代入, 得 n 次代數方程式 f (x) = xn− c1xn−1− c2xn−2− · · · − ckxn−k = 0.
因為零解不是我們所要的, 所以可消去 xn−1 的因式, 得到 k 次代數方程式
f (x) = xk− c1xk−1− c2xk−2− · · · − ck = 0, (3) 則 (3) 式稱為 (2) 式的特徵方程式, 其中 (3) 式的每一個根都稱為特徵複根。 由代數基本定理 知最多有 k 個相異特徵複根, 可令 k 個特徵複根為 αi (1≤ i ≤ k) 且 αi 為最大特徵正實根, 由於 c1, c2, . . . , ck > 0, 則有 α1 ≥ |αi|, i = 2, 3, . . . , k, 這性質讀者可自行證明。
例如: 數列 {an} : an= an−1+ an−2+ an−3+ 2an−4, 特徵方程式為 f (x) = x4− x3− x2− 2 = 0,
⇒ f(x) = (x − 2)(x + 1)(x2 + 1) = 0,
所以特徵複根為 2, −1, ±i, 即 2 > | − 1| = |i| = | − i| 是滿足 α1 ≥ |α2| ≥ |α3| ≥ |α4|。
定理1: 設 {an} 為 k 階正實係數齊次線性遞迴數列, 若有 k1 個特徵實根為 αi,且 α1 為最大 特徵正實根, 其中 α1 ≥ αi (1≤ i ≤ k1 ≤ k), 則 α1為數列 {an} 的唯一特徵正實根, 且實根 α1, α2, . . . , αk 均相異 (沒有重根的情形)。
證明: 現在利用笛卡兒符號法則 (Descartes’s rule of signs) 來判定特徵方程式的正實根的個 數。 首先觀察 f(x) = 0 的係數 (1, −c1,−c2, . . . ,−ck)之正負符號關係為 (+, −, −, . . . , −), 由於第一項及第二項符號改變, 即係數符號的變化次數 N+ 為 1。 令正實根的個數為 N, 則由 笛卡兒符號法則知 N+(= 1)≥ N 並且 1 − N 為非負偶數, 所以 N = 1, 故 α1 為特徵方程 式 f(x) = 0 的唯一正實根。
另一方面, 若實數 α 為數列的特徵根, 則 f(α) = 0, 故 f(α) = αk−c1αk−1−· · ·−ck = 0,但當 α > 0 時,
f(α) = kαk−1− c1(k− 1)αk−2− · · · − ck−1
> k[αk−1− c1αk−2− · · · − ck−1] = k
f (α) + ck α
= kck α > 0,
又當 α < 0 時, 若 k = 3, α = −1, c1 = 100, c2 = 1, c3 = 2, 則 f(α) = 3α2−2c1α−c2。 考慮 f(α)− 3[α2− c1α− c2] = c1α + 2c2 =−98 < 0, 得 f(α) < 0。
所以 f(α)= 0, 故 f(x) = 0 的實根沒有重根的情形, 即其特徵實根 α1, α2, . . . , αk1 均
相異。
由上可知 f(x) = 0 的特徵實根均相異 (沒有重根), 在探討 (2) 式的一般解時, 分成 k 個 相異特徵實根或有特徵共軛複根來討論。 事實上, 特徵共軛複根的情形只是相異特徵實根的一 個特例。 不失一般性, 底下僅以出現 k − 2 個相異特徵實根 α1, α2, . . . , αk−2 及一組特徵共軛 複根 αk−1, αk 來證明。
定理2: 設 {an} 為 k 階正實係數齊次線性遞迴數列,
(i) 若 α1, α2, . . . , αk 為數列 {an} 的特徵相異實根, 則數列 {an} 的一般解為
an = d1αn1 + d2αn2 +· · · + dkαnk, (4) 其中 di (i = 1, 2, . . . , k) 為不為零的任意常數。
(ii) 若 α1, α2, . . . , αk−2 為數列 {an} 的特徵相異實根且有一組特徵共軛複根 αk−1, αk, 則 數列 {an} 的一般解為
an= d1α1n+ d2αn2 +· · · + dk−2αnk−2+ (dk−1+ dk)γncos nθ + i(dk−1− dk)γnsin nθ, 其中 di (i = 1, 2, . . . , k) 為不為零的任意常數。
證明: (i) 給定首 k 項的初期值為
a1 = 1, a2 = β2, a3 = β3, a4 = β4, . . . , ak= βk (βk 不為零, i = 1, 2, . . . , k) 事實上, 由 (4) 式得知
線性方程組 Γ :
⎧⎪
⎪⎪
⎪⎨
⎪⎪
⎪⎪
⎩
d1α1+ d2α2+· · · + dkαk= 1 d1α21+ d2α22+· · · + dkα2k= β2
... ... d1αk1 + d2αk2 +· · · + dkαkk= βk
為關於 di (i = 1, 2, . . . , k) 的線性方程組, 其係數行列式為
α1 α2 · · · αk α12 α22 · · · α2k ... ... . .. ... αk1 αk2 · · · αkk
= α1· α2· · · αk
1 1 · · · 1 αk−11 αk−12 · · · αk−1k
... ... . .. ... αk1 αk2 · · · αkk
.
由 Vandermonde 行列式且 Cramer 法則得知係數行列式不為零及 di = 0 (因為特徵根 α1, α2, . . . , αk 均相異), 因此, 線性方程組 Γ 必有唯一一組解 di (i = 1, 2, . . . , k),故 (4) 式 為 (3) 式的一般解。
(ii) 事實上, 特徵共軛複根的情形只是相異特徵實根的一個特例。
令 αk−1= δ1+ iδ2, αk = δ1− iδ2, 其中 γ =
δ12+ iδ22, θ = tan−1(δ2/δ1), 則 dk−1αnk−1+ dkαkn= dk−1(δ1+ iδ2)n+ dk(δ1− iδ2)n
= dk−1γn(cos nθ + i sin nθ) + dkγn(cos nθ− i sin nθ)
= (dk−1+ dk)γncos nθ + i(dk−1− dk)γnsin nθ.
因此, 數列 {an} 的一般解
an= d1αn1 + d2αn2 +· · · + dk−2αnk−2+ (dk−1+ dk)γncos nθ + i(dk−1− dk)γnsin nθ, 其中 di (i = 1, 2, . . . , k) 為不為零的任意常數。 註: 由定理 2 知我們要討論的線性遞迴關係式都是 k 階線性遞迴關係式, 也就是在一般解中的 每一個係數 di 均不為 0。
例如: 數列 {an} : an= an−1+ an−2+ an−3+ 2an−4, 特徵複根為 2, −1, ±i。
因為 ±i = cos(π/2) ± i sin(π/2), 所以數列 {an} 的一般解為
an= d1(2)n+ d2(−1) + (d3+ d4) cos(nπ/2) + i(d3− d4) sin(nπ/2). (5) 由於初期值為 a−2 = 0, a−1 = 0, a0 = 0, a1 = 1, 代入 (5) 式得知
⎧⎪
⎪⎪
⎪⎨
⎪⎪
⎪⎪
⎩
a−2= d1(2)−2+ d2(−1)−2− (d3+ d4) = 0 a−1= d1(2)−1+ d2(−1)−1− i(d3− d4) = 0 a0= d1(2)0+ d2(−1)0+ (d3+ d4) = 0 a1= d1(2)1+ d2(−1)1+ i(d3− d4) = 1
⇒
⎧⎪
⎪⎪
⎪⎨
⎪⎪
⎪⎪
⎩
a−2= 14d1+ d2− d3− d4 = 0 a−1=12d1− d2− i(d3− d4) = 0 a0= d1+ d2+ d3+ d4 = 0 a1= 2d1− d2+ i(d3− d4) = 1 .
求解得 d1 = 4
15, d2 = 1
6, d3+ d4 =− 1
10, i(d3− d4) = 3
10, 故數列 {an} 的一般解為 an = 4
15(2)n− 1
6(−1)n− 1
10cos(nπ/2) + 3
10sin(nπ/2).
乙、 特徵正實根等於數列的後前項比的極限
「費氏數列的後前項比的極限等於特徵正實根, 其值也是黃金比列 φ。」 這個性質是德國天 文學家克卜勒 (Johannes Kepler, 1571∼1630) 與費氏數列、 黃金比列 φ 的不期而遇。 如何 觀察出來呢?
他將黃金比列寫成連分數
φ = 1 + 1 1 + 1
1 + 1 1 +· · ·
.
若逐項分斷連分數, 然後一步一步往下算去, 不僅可以估計出黃金比列 φ, 美妙地可發現費氏數 列的後前項比, 如下一系列值:
1 = 1
1 = 1.00000, 1 + 1 1 = 2
1 = 2.00000, 1 + 1
1 + 1 = 3
2 = 1.50000, 1 + 1
1 + 1 1 + 1
= 5
3 = 1.66666, 1 + 1 1 + 1
1 + 1 1 + 1
= 8
5 = 1.60000,
1 + 1
1 + 1
1 + 1 1 + 1
1 + 1
= 13
5 = 1.62500.
從這些逐項近似值中, 發現費氏數列的後前項比的極限等於黃金比列 φ。 推廣至 k 階時, 數列 的後前項比的極限等於唯一特徵正實根, 底下就來證明。
定理3: 設 {an} 為 k 階正實係數齊次線性遞迴數列,
(i) 若 α1, α2, . . . αk為數列 {an} 的特徵相異實根, 其中 α1 為唯一特徵正實根並且 α1 ≥ αi, i = 2, 3 . . . , k, 則 lim
n→∞
an+1
an = α1。
(ii) 若 α1, α2, . . . αk 為數列 {an} 的特徵相異實根且有一組特徵共軛複根 αk−1, αk, 其中 α1為唯一特徵正實根並且 α1 ≥ |αi|, i = 2, 3 . . . , k, 則 lim
n→∞
an+1
an = α1。
證明: (i) 由定理 2 得知數列 {an} 的一般解為 an = d1αn1 + d2α2n+· · · + dkαnk, 其中 di
(i = 2, 3, . . . , k)為不為零的任意常數。
由於 α1 ≥ αi, i = 2, 3, . . . , k,因此,
n→∞lim an+1
an = lim
n→∞
d1αn+11 + d2αn+12 +· · · + dkαn+1k d1αn1 + d2αn2 +· · · + dkαnk
= lim
n→∞
d1α1+ d2α2
α2 α1
n
+ d3α3
α3 α1
n
+· · · + dkαk
αk α1
n
d1+ d2
α2 α1
n + d3
α3 α1
n
+· · · + dk
αk α1
n = α1. (6)
(ii) 令 αk−1 = δ1+ iδ2, αk= δ1− iδ2,其中 γ =
δ12+ iδ22, θ = tan−1(δ2/δ1)。 因為 α1 >|αi|, i = 2, 3, . . . , k, 所以 α1 >|αk−1| = |αk| = γ ⇒ 0 < γ
α1 < 1。 由定理 2 得知數列 {an} 的一般解為
an= d1α1n+ d2αn2 +· · · + dk−2αnk−2+ (dk−1+ dk)γncos nθ + i(dk−1− dk)γnsin nθ 其中 di (i = 1, 2, . . . , k 為不為零的任意常數。 由於
n→∞lim
(dk−1+ dk)γn+1cos(n + 1)θ + i(dk−1− dk)γn+1sin(n + 1)θ (dk−1+ dk)γncos nθ + i(dk−1− dk)γnsin nθ
= lim
n→∞
(dk−1+ dk)γ
γ α1
n
cos(n + 1)θ + (dk−1− dk)γ
γ α1
n
sin(n + 1)θ (dk−1+ dk)
γ α1
n
cos nθ + (dk−1− dk)
γ α1
n sin nθ
= 0, (7)
因此, 由 (6) 式與 (7) 式知
n→∞lim an+1
an
= lim
n→∞
d1αn+11 +· · ·+dk−2αn+1k−2+(dk−1+dk)γn+1cos(n+1)θ+i(dk−1−dk)γn+1sin(n+1)θ d1αn1+· · ·+dk−2αk−2n +(dk−1+dk)γncos nθ+i(dk−1− dk)γnsin nθ
= α1.
推論: 若不僅有一組特徵共軛複根, 且 α1為唯一特徵正實根並且 α1 >|αi|, i = 2, 3, . . . , k, 則極限 lim
n→∞
an+1
an 必存在, 且 lim
n→∞
an+1
an = α1。
2. 幾何建構問題
若依照費氏螺線的樣式 (逆時針旋轉), 推廣 k 階齊次線性遞迴數列時產生的曲線均為螺 線嗎?
本文將 k 階時產生的曲線稱為數列螺線。 顯然地, 數列螺線的建構條件為
數列 {an} 中每一項皆為正實數且為 (嚴格) 遞增數列, 即 0 < a1 ≤ a2 ≤ a3 ≤ · · · 。
那麼當 c1, c2, . . . , ck ≥ 1 時, 顯然是 0 < a1 ≤ a2 ≤ a3 ≤ · · · , 並且數列螺線是保持費 氏螺線的樣式建構方式, 以數列 {an} : an = an−1+ 2an−2 為例, 參見圖 2。
圖2 : 數列螺線 圖3: 第2項後成為數列螺線
事實上, 如果將 c1, c2, . . . , ck ≥ 1 改為 c1 + c2+· · · + ck ≥ 1, 則數列在某項後皆為 (嚴格) 遞增數列, 證明參見定理 5。
例如: 數列 {an} : an = 1
3an−1+1
3an−2+1
2an−3, 是從第 2 項後開始嚴格遞增, 即第 2 項 後成為由內而外逆時鐘旋出的螺線, 參見圖 3, 此情形也列入本文數列螺線的範疇。
此外, 若考慮 0 < c1+ c2+· · · + ck < 1, 則數列在某項後為 (嚴格) 遞減數列, 證明參 見定理 5。 若採用寬廣觀點, 同樣地也可建構出數列螺線, 則數列螺線是在某項後成為由外而內 逆時鐘旋入的螺線
例如: 數列 {an} : an = 1
3an−1+1
3an−2+1
4an−3, 是從第 6 項後開始嚴格遞減, 即第 6 項 後成為由外而內逆時鐘旋入的螺線, 參見圖 4, 此情形也列入本文數列螺線的範疇。
綜合以上, 本文所談數列螺線包含由內而外逆時鐘旋出或由外而內逆時鐘旋入的螺線。
圖4 : 第 6 項後才形成數列螺線
值得一提, 針對 ai (i = 2, 3, . . . , n) 為負實數或 0, 不妨將 「ai 的負號」 代表曲線方向與
「ai 的正號」 相反, 並且凹性是相反的。 而 「ai = 0」 是代表曲線沒有移動, 配合費氏螺線的樣 式建構方式, 則必出現尖點產生, 其對應曲線必不是螺線。
例如: 數列 {an} : an = −an−1− an−2 中 a1 = 1, a2 =−1, a3 = 0, a4 = 1, 參見圖 5, 由於 a1 與 a2 是正負相間, 對應四分之一的圓弧間必出現尖點。 而 a3 = 0 代表曲線沒有移動, 相鄰 a2 與 a4 所對應四分之一的圓弧間必出現尖點, 因此, 其對應曲線必不是螺線。
圖5: 數列非螺線
我們將數列螺線的方程式用參數式來表示, 為了方便, 令數列螺線的起點坐標為 (1, 0) 且 第一個四分之一的圓弧之圓心坐標為 (0, 0)。
定理4: 設 {an} 為 k 階正實數齊次線性遞迴數列, 則數列螺線參數式為
x(t) = 1 + [t]
i=1
δi(ai − ai−1) + a[t]· cost− 1 2 π
,
y(t) = [t]
i=1
δi(ai− ai−1) + a[t]· sint− 1 2 π
,
其中實數 t≥1, δi=1 + (−1)i+1
2 × (−1)[t+12 ], δi=1 + (−1)i
2 × (−1)[t2], [ ] 為高斯符號。
證明: 利用圓的參數式來推導數列螺線參數式。
首先先考慮圓心的移動量 : 數列螺線從 (1, 0) 起, 每個四分之一的圓弧的圓心以 ai− ai−1 單 位移動 (0 ≤ i ≤ n), 但要注意 x, y 坐標是否移動 ai− ai−1 單位, 我們發現除了前四個圓心 外, 之後圓心軌跡就遵循 「向左、 向下、 向右與向上」 等四個循環方向移動, 參見下表, 其中用
−1, 0, 1來表示 x 軸方向-向左、 沒有移動、 向右和 y 軸方向-向下、 沒有移動、 向上。
圓心移動軌跡 向左 向下 向右 向上
x 軸方向 −1 0 +1 0
y 軸方向 0 −1 0 +1
上表可用數學式子來表示, 其中 [ ]為高斯符號:
(i) 由 δi=1 + (−1)i+1
2 ×(−1)[t+12 ] 來決定 x 坐標是否移動 ai−ai−1單位 (i = 1, 2, 3, . . .)。
(ii) 由 δi=1 + (−1)i
2 × (−1)[2t] 來決定 y 坐標是否移動 ai− ai−1 單位 (i = 1, 2, 3, . . .)。
由上可知若參數 t 為大於 1 的實數, 則圓心的坐標為
1 +
[t]
i=1
δi(ai − ai−1), [t]
i=1
δi(ai− ai−1)
, 其中 [ ] 為高斯符號。
又四分之一圓的半徑為 a[t] 且旋轉角度為 t− 1
2 π, 則由圓的參數式知數列螺線參數式為
x(t) = 1 + [t]
i=1
δi(ai− ai−1) + a[t]· cost− 1 2 π
,
y(t) = [t]
i=1
δi(ai− ai−1) + a[t]· sint− 1 2 π
,
其中 t≥1, δi=1 + (−1)i+1
2 × (−1)[t+12 ], δi=1 + (−1)i
2 × (−1)[2t], [ ] 為高斯符號。
3. 收斂性問題
我們發現 : 當係數 c1, c2, . . . , ck ≥ 1 時, 數列 {an} 為嚴格遞增, 則必可建構出數列螺 線, 但係數 c1, c2, . . . , ck < 1 時, 數列 {an} 有遞增也有遞減, 無法區分怎樣的數列螺線, 以 下定理 5 我們提供最佳的充分條件就是由 c1+ c2+· · · + ck 的值來決定。
定理5: 設 {an} 為 k 階正實係數齊次線性遞迴數列, 則數列在某項後, 有以下三種數列曲線:
(i) 當 c1+ c2+ c3+· · · + ck > 1 時, 則 {an} 為遞增數列, 數列曲線是由內而外逆時鐘旋 出的螺線。
(ii) 當 c1+ c2+ c3+· · · + ck = 1 時, 則 {an} 趨近常數數列, 數列曲線是由內而外逆時鐘 旋轉的螺線。 當 n → ∞ 時, 數列螺線會收斂於一圓 (此圓稱為極限圓)。
(iii) 當 0 < c1+ c2+ c3+· · · + ck < 1 時, 則 {an} 為遞減數列, 數列曲線是由外而內逆時 鐘旋入的螺線。 當 n → ∞ 時, 數列螺線會旋入至一點 (此點稱為點圓)。
註: 極限圓及點圓為數列螺線的退化情形。
證明: 由數列的特徵方程式如 (3) 式, 當 x = 1 代入 (3) 式, 得到 f (1) = 1− (c1+ c2+· · · + ck).
(i)當 c1+ c2+ c3+· · · + ck> 1 時, 則
f (1) = 1− (c1+ c2+· · · + ck) < 0;
另一方面, 因為 (3) 式的領導係數為 1, 所以正實根 α1 > 1, 參見圖 6(a), 再由定理 3 可知
n→∞lim an+1
an = α1> 1。 因此, 數列{an}在某項後會是遞增的, 數列曲線是由內而外逆時鐘旋出的 螺線。
圖6: 數列的特徵方程式的增減情形 (ii) 當 c1+ c2 + c3+· · · + ck = 1 時, 則由定理 3 可知:
f (1) = 1− (c1+ c2 +· · · + ck) = 0, 並且正實根 α1 = lim
n→∞
an+1
an = 1, 參見圖 6(b)。
因此, 數列在某項後會是趨近常數, 數列曲線是由內而外逆時鐘旋轉的螺線。 另一方面, 當 n →
∞ 時, 由於 an→ 常數, 數列螺線會收斂於一圓。
(iii) 當 0 < c1 + c2+· · · + ck< 1 時, 則由定理 3 可知:
f (1) = 1− (c1+ c2+· · · + ck) > 0, 並且正實根 α1 = lim
n→∞
an+1
an < 1, 參見圖 6(c)。
因此, 數列在某項後會是遞減的, 數列曲線是由外而內逆時鐘旋入的螺線。 另一方面, 當 n →
∞ 時, 數列螺線會旋入至一點。
4. 結語
費氏螺線的誕生早已與宇宙、 黃金比例共存著, 為大自然偉大的法則, 這樣的生長現象並 非是偶然的, 近代的數學家已研究出其數學模式的成因, 可說:
大自然的奧秘或線索藏在數學裡。
值得一提, 費氏螺線近似於黃金螺線 (是對數螺線的一種), 數學家笛卡兒 (R. Descartes,
1596∼1650) 稱它為等角螺線, 說明了對數螺線的獨一無二的特性, 也因此對應了許多美學與 生物的生長模式。
事實上, 本文數列螺線跟隨圓的半徑的增減擴充而產生, 除了如費氏螺線的模式是由內而 外逆時鐘旋出的螺線外, 還擴充有由外而內逆時鐘旋入的螺線, 反映出多樣對數螺線的新面貌, 皆近似於等角螺線, 或許大自然早已出現此螺線的生長或運動行為模式。
例如: 「飛蛾撲火」 的運動行為模式正是數列螺線會旋入至一點, 參見圖 7, 可說數列螺線早已 出現在大自然裡。 此外, 大自然提供如此多的線索, 等待我們挖掘並且證明之, 我們繼續追隨它, 使之美麗。
圖7: 飛蛾撲火