Ludium
Sign In
Linear Algebra
Solving Linear Systems
01Row Picture, Column Picture & Matrix Form02Row Picture vs Column Picture in 3D03Singular vs Non-Singular Matrices04Two Ways to Compute AxProblem set0/10Problem set 20/10MIT problem set0/4Practice∞
01Gaussian Elimination & Back Substitution02Elimination Matrices03Permutation Matrices and InversesProblem set0/10Problem set 20/10MIT problem set0/4Practice∞
015 Views of Matrix Multiplication02Matrix Inverses and Singularity03Gauss-Jordan Elimination for InversesProblem set0/10Problem set 20/10MIT problem set0/5Practice∞
01Product Inverse & Transpose Rules02A = LU Factorization03Gaussian Elimination: ⅓n³ Operation Count04Permutation Matrices: P⁻¹ = PᵀProblem set0/10Problem set 20/10MIT problem set0/6Practice∞
01Permutation Matrices and PA = LU02Why RᵀR Is Always Symmetric03Vector Spaces and Subspaces04The Column Space of a MatrixProblem set0/10Problem set 20/10MIT problem set0/4Practice∞
Vector Spaces & Subspaces
01Vector Spaces and Subspaces02Column Space of a Matrix03Null Space of a MatrixProblem set0/10Problem set 20/10Practice∞
01Rank, Pivots, and Free Variables02Special Solutions of the Null Space03Reading the Null Space from RREFProblem set0/10Problem set 20/10MIT problem set0/3Practice∞
01Solvability of Ax = b02Particular Solutions and the Null Space03Rank and the Four Cases of Ax = bProblem set0/10Problem set 20/10MIT problem set0/8Practice∞
01Linear Independence & the Null Space02Basis and Dimension03The Rank-Nullity TheoremProblem set0/10Problem set 20/10MIT problem set0/5Practice∞
01Rank and the Four Fundamental Subspaces02Four Bases from One Row Reduction03Matrices as Vector SpacesProblem set0/10Problem set 20/10MIT problem set0/4Practice∞
01The Dimension Formula for Subspaces02Differential Equations as Linear Algebra03Rank-One Matrices and Outer Products04The Four Fundamental SubspacesProblem set0/10Problem set 20/10MIT problem set0/1Practice∞
01The Incidence Matrix as Difference Operator02The Incidence Matrix of a Graph03The Equilibrium Equation AᵀCAx = fProblem set0/10Problem set 20/10MIT problem set0/1Practice∞
Orthogonality & Least Squares
01Orthogonal Subspaces02Null Space as Orthogonal Complement03The Normal EquationsProblem set0/10Problem set 20/10MIT problem set0/4Practice∞
01The Projection Matrix P = aaᵀ/aᵀa02Projection onto Subspaces03Projection and Least SquaresProblem set0/10Problem set 20/10MIT problem set0/6Practice∞
01Complementary Projections: P and I−P02Least Squares and the Normal Equations03Least Squares as Projection04Invertibility of AᵀAProblem set0/10Problem set 20/10MIT problem set0/5Practice∞
01Orthogonal Matrices02Projections onto Orthonormal Bases03The Gram-Schmidt Process04QR Factorization from Gram-SchmidtProblem set0/10Problem set 20/10MIT problem set0/3Practice∞
Determinants
01The Three Axioms of the Determinant02Determinant Properties from Three Axioms03Determinants via Elimination04Multiplicative and Transpose PropertiesProblem set0/10Problem set 20/10MIT problem set0/2Practice∞
01The Permutation Formula for Determinants02Cofactor Expansion03Periodic Determinants of Tridiagonal MatricesProblem set0/10Problem set 20/10MIT problem set0/3Practice∞
01The Cofactor Formula for Matrix Inverses02Cramer's Rule03Determinants as VolumeProblem set0/10Problem set 20/10MIT problem set0/4Practice∞
Eigenvalues & Dynamics
01Eigenvalues and Eigenvectors02Finding Eigenvalues and Eigenvectors03Complex Eigenvalues and Defective MatricesProblem set0/10Problem set 20/10MIT problem set0/2Practice∞
01Diagonalization: A = SΛS⁻¹02Matrix Powers and the Stability Theorem03Defective Matrices and Multiplicity04Eigenvalues and the Fibonacci SequenceProblem set0/10Problem set 20/10MIT problem set0/3Practice∞
01Solving du/dt = Au with Eigenvalues02Eigenvalue Stability in the Complex Plane03The Matrix Exponential e^(At)04The Companion MatrixProblem set0/10Problem set 20/10MIT problem set0/5Practice∞
01Why 1 Is Always an Eigenvalue of Markov Matrices02Markov Chain Steady States via Eigenvalues03Markov Matrices and Steady States04Orthonormal Bases and Fourier CoefficientsProblem set0/10Problem set 20/10MIT problem set0/7Practice∞
Symmetric Matrices & the SVD
01The Spectral Theorem02Why Symmetric Matrices Have Real Eigenvalues03Sylvester's Law of Inertia04Testing Positive DefinitenessProblem set0/10Problem set 20/10MIT problem set0/5Practice∞
01The Conjugate Transpose02Inverse of the Fourier Matrix03The FFT FactorizationProblem set0/10Problem set 20/10Practice∞
01Positive Definite Matrices02Positive Definiteness via Pivots03The Hessian and Positive Definiteness04Principal Axis Theorem and the EllipsoidProblem set0/10Problem set 20/10MIT problem set0/6Practice∞
01Closure Properties of Positive Definite Matrices02Matrix Similarity and Invariant Eigenvalues03Why Eigenvalues Don't Classify Similarity04Jordan Canonical FormProblem set0/10Problem set 20/10MIT problem set0/5Practice∞
01The Singular Value Decomposition02Reducing SVD to A^T A03Computing the SVD: Two Worked Examples04SVD and the Four Fundamental SubspacesProblem set0/10Problem set 20/10MIT problem set0/3Practice∞
Linear Transformations & Applications
01Defining Linear Transformations02Basis, Coordinates, and Linear Maps03How a Linear Map Becomes a MatrixProblem set0/10Problem set 20/10Practice∞
01Image Compression as a Change of Basis02Choosing a Basis: JPEG, Fourier, and Wavelets03Change of Basis and Image Compression04Change of Basis and DiagonalizationProblem set0/10Problem set 20/10Practice∞
01Matrix Inverses and the Four Subspaces02One-Sided Inverses and Projections03The Pseudo-Inverse04The Pseudo-Inverse via the SVDProblem set0/10Problem set 20/10Practice∞

Eigenvalues and the Fibonacci Sequence

Solve uₖ₊₁ = Auₖ by expanding in eigenvectors; the Fibonacci recurrence has eigenvalues (1±√5)/2, giving the golden-ratio growth rate.


Loading…

←Previous Defective Matrices and MultiplicityNext Solving du/dt = Au with Eigenvalues →

Your summary note

    1. 1

      Solving uk+1=Auku_{k+1} = Au_kuk+1​=Auk​ by eigenvector expansion

      Write u0=c1x1+⋯+cnxn=Scu_0 = c_1x_1 + \cdots + c_nx_n = Scu0​=c1​x1​+⋯+cn​xn​=Sc, apply AAA once to see each term pick up its own λi\lambda_iλi​, and reach uk=∑ciλikxi=SΛkcu_k = \sum c_i\lambda_i^k x_i = S\Lambda^k cuk​=∑ci​λik​xi​=SΛkc.

    2. 2

      Fibonacci as a 2×22 \times 22×2 first-order system

      Set uk=(Fk+1,Fk)Tu_k = (F_{k+1}, F_k)^Tuk​=(Fk+1​,Fk​)T, pair the trivial equation Fk+1=Fk+1F_{k+1} = F_{k+1}Fk+1​=Fk+1​ with Fk+2=Fk+1+FkF_{k+2} = F_{k+1} + F_kFk+2​=Fk+1​+Fk​, and record the matrix with rows (1,1)(1,1)(1,1) and (1,0)(1,0)(1,0).

    3. 3

      The eigenvalues 1±52\frac{1 \pm \sqrt{5}}{2}21±5​​

      Take det⁡(A−λI)=λ2−λ−1\det(A - \lambda I) = \lambda^2 - \lambda - 1det(A−λI)=λ2−λ−1, solve by the quadratic formula, check the roots against trace 111 and determinant −1-1−1, and record x=(λ,1)Tx = (\lambda, 1)^Tx=(λ,1)T.

    4. 4

      Golden ratio growth rate of FkF_kFk​

      Write Fk≈c1λ1kF_k \approx c_1\lambda_1^kFk​≈c1​λ1k​ with λ1≈1.618\lambda_1 \approx 1.618λ1​≈1.618 and the λ2≈−0.618\lambda_2 \approx -0.618λ2​≈−0.618 term shrinking away, then solve u0=(1,0)T=c1x1+c2x2u_0 = (1,0)^T = c_1x_1 + c_2x_2u0​=(1,0)T=c1​x1​+c2​x2​ for the constants.

    Attempt 1 of 2