Skip to main content

Fibonacci difference equation: limiting ratio 1.6180

Short answer

Solve via z-transform: X(z)=z/(z^2–z–1). Partial fractions give x(k)=(1/√5)[((1+√5)/2)^k – ((1–√5)/2)^k]. The dominant root φ=(1+√5)/2, so lim x(k+1)/x(k)=φ=1.6180.

Find X(k) when x(k+2)=x(k+1)+x(k), x(0)=0 and x(1)=1, and show that the limiting value of x(k+1)/x(k)=1.6180 when k tends to infinity.

This is the Fibonacci equation x(k+2)−x(k+1)−x(k)=0x(k+2) - x(k+1) - x(k) = 0. Solve it by z-transform, then take the limit of the ratio.

Step 1: Take the z-transform

Use Z[x(k+2)]=z2X(z)−z2x(0)−zx(1)\mathcal{Z}[x(k+2)] = z^2X(z) - z^2x(0) - zx(1) and Z[x(k+1)]=zX(z)−zx(0)\mathcal{Z}[x(k+1)] = zX(z) - zx(0), with x(0)=0x(0) = 0 and x(1)=1x(1) = 1:

[z2X(z)−z]−zX(z)−X(z)=0\left[z^2X(z) - z\right] - zX(z) - X(z) = 0 X(z)=zz2−z−1X(z) = \frac{z}{z^2 - z - 1}

Step 2: Factorise the denominator

z2−z−1=0  ⇒  z1,2=1±52,z1=1.6180,  z2=−0.6180z^2 - z - 1 = 0 \;\Rightarrow\; z_{1,2} = \frac{1 \pm \sqrt{5}}{2}, \quad z_1 = 1.6180,\; z_2 = -0.6180

Step 3: Partial fractions of X(z)/zX(z)/z

X(z)z=1(z−z1)(z−z2)=1z1−z2[1z−z1−1z−z2],z1−z2=5\frac{X(z)}{z} = \frac{1}{(z-z_1)(z-z_2)} = \frac{1}{z_1 - z_2}\left[\frac{1}{z-z_1} - \frac{1}{z-z_2}\right], \quad z_1 - z_2 = \sqrt{5} X(z)=15[zz−z1−zz−z2]X(z) = \frac{1}{\sqrt5}\left[\frac{z}{z - z_1} - \frac{z}{z - z_2}\right]

Step 4: General (closed-form) solution

x(k)=15[(1+52)k−(1−52)k],k=0,1,2,…x(k) = \frac{1}{\sqrt5}\left[\left(\frac{1+\sqrt5}{2}\right)^k - \left(\frac{1-\sqrt5}{2}\right)^k\right], \quad k = 0, 1, 2, \dots

Here 1/5=0.44721/\sqrt5 = 0.4472. Check: x(0)=0x(0) = 0, and x(1)=15⋅5=1x(1) = \frac{1}{\sqrt5}\cdot\sqrt5 = 1.

First terms: 0,1,1,2,3,5,8,13,21,34,…0, 1, 1, 2, 3, 5, 8, 13, 21, 34, \dots

Step 5: Limiting value of x(k+1)/x(k)x(k+1)/x(k)

Let z1=1.6180z_1 = 1.6180 and z2=−0.6180z_2 = -0.6180:

x(k+1)x(k)=z1k+1−z2k+1z1k−z2k=z1−z2(z2z1)k1−(z2z1)k\frac{x(k+1)}{x(k)} = \frac{z_1^{k+1} - z_2^{k+1}}{z_1^{k} - z_2^{k}} = \frac{z_1 - z_2\left(\dfrac{z_2}{z_1}\right)^{k}}{1 - \left(\dfrac{z_2}{z_1}\right)^{k}}

Since ∣z2z1∣=0.6181.618=0.382<1\left|\dfrac{z_2}{z_1}\right| = \dfrac{0.618}{1.618} = 0.382 < 1, we get (z2z1)k→0\left(\dfrac{z_2}{z_1}\right)^k \to 0 as k→∞k \to \infty:

lim⁡k→∞x(k+1)x(k)=z1=1+52=1.6180\lim_{k\to\infty}\frac{x(k+1)}{x(k)} = z_1 = \frac{1+\sqrt5}{2} = 1.6180

This is the golden ratio.

Numerical check

kk56789
x(k+1)/x(k)x(k+1)/x(k)8/5 = 1.60013/8 = 1.62521/13 = 1.615434/21 = 1.619055/34 = 1.6176

The ratio oscillates about 1.6180 and settles there.

Answer: x(k)=15[(1.6180)k−(−0.6180)k]x(k) = \dfrac{1}{\sqrt5}\left[(1.6180)^k - (-0.6180)^k\right], and lim⁡k→∞x(k+1)/x(k)=1+52=1.6180\lim_{k\to\infty} x(k+1)/x(k) = \dfrac{1+\sqrt5}{2} = 1.6180.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Digital Control System chapter 2 questions →