Hint $ $ This method works for any sequence $f_i\,$ satisfying a $\rm\color{#0a0}{constant}$ coefficient linear $\,k$'th order recurrence, since then the shift map $\, S\, (f_k,\,\ldots,f_2,f_1)\, =\, (f_{k+1},\ldots,f_3,f_2)\,$ has matrix being a $\rm\color{#0a0}{constant}$ coefficient (companion) matrix. So $\,S^n,\,$ a shift by $\,n,\,$ has matrix an $\,n$'th $\rm\color{#0a0}{power}$ of said shift matrix. If instead the linear recurrence had $\rm\color{#c00}{variable}$ (nonconstant) coef's (i.e. coef's that depend on the index $\,k),\,$ then the shift $\,S^n$ would not be an $\:\!n$'th $\rm\color{#0a0}{power}$ of a $\rm\color{#0a0}{constant}$ coef matrix but, rather, a (more complicated) $\rm\color{#c00}{product}$ of $\,n\,$ matrices with $\rm\color{#c00}{variable}$ coef's.
Note that this matrix representation yields addition formulas and fast polynomial-time algorithms for computing the (Fibonacci) sequence by computing matrix powers by repeated squaring.