Shift Operators and Recurrence Relations
mathsequencerecurrenceshift-operatorlinear-operatorlecture
1Introduction
This lecture unifies differences, indefinite sums, and recurrences through the forward shift operator
(Ea)_n=a_{n+1}.
A constant-coefficient linear recurrence then becomes the linear operator equation P(E)a=b.
2Terminology and Definitions
Work with bilateral sequences indexed by \mathbb Z or unilateral sequences indexed by n\ge0, as specified. Define
(Ea)_n=a_{n+1},\qquad (Ia)_n=a_n,
and the forward difference
(\Delta a)_n=a_{n+1}-a_n,\qquad \Delta=E-I.
This convention agrees with the Z-transform lecture: the bilateral transform maps E to multiplication by z.
3Program
The development proceeds by proving linearity of E, deriving difference identities from \Delta=E-I, treating \Delta F=f as an inverse problem, expressing recurrences as P(E)a=b, and deriving the characteristic equation from the eigen-sequence r^n.
data/lecture/math/sequence/recurrence-types-and-strategies.lecture.n.md
data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md
data/lecture/math/linear-operator/polynomial-operators-and-eigenvalue-problems.lecture.n.md
4Proposition 1: E Is a Linear Operator
4.1Statement
The shift E is linear on the sequence space.
4.2Proof
For sequences a,b and scalars \alpha,\beta,
\begin{aligned}
(E(\alpha a+\beta b))_n
&=(\alpha a+\beta b)_{n+1}\\
&=\alpha a_{n+1}+\beta b_{n+1}\\
&=(\alpha Ea+\beta Eb)_n.
\end{aligned}\begin{aligned}
(E(\alpha a+\beta b))_n
&=(\alpha a+\beta b)_{n+1}\\
&=\alpha a_{n+1}+\beta b_{n+1}\\
&=(\alpha Ea+\beta Eb)_n.
\end{aligned}
The identity holds for every admissible n, hence E(\alpha a+\beta b)=\alpha Ea+\beta Eb. \square
5The Difference Is E-I
Since
((E-I)a)_n=a_{n+1}-a_n=(\Delta a)_n,
\Delta=E-I. Therefore \Delta is linear as the difference of linear operators:
\Delta(\alpha a+\beta b)=\alpha\Delta a+\beta\Delta b.
6Higher Differences from the Binomial Theorem
Because E and I commute,
\Delta^m=(E-I)^m
=\sum_{k=0}^m(-1)^k\binom{m}{k}E^{m-k}.
Consequently,
(\Delta^m a)_n
=\sum_{k=0}^m(-1)^k\binom{m}{k}a_{n+m-k}.
Thus the signs and indices follow directly from the operator binomial expansion.
7Indefinite Summation as the Inverse Difference Problem
An indefinite sum seeks F satisfying \Delta F=f, or F(n+1)-F(n)=f(n). Summing from n=a to b-1 telescopes:
\sum_{n=a}^{b-1}f(n)
=\sum_{n=a}^{b-1}(F(n+1)-F(n))
=F(b)-F(a).
This is the discrete counterpart of \int_a^bF'(x)\,dx=F(b)-F(a).
8The Shift Term in a Product Difference
Direct calculation gives
\begin{aligned}
\Delta(uv)(n)
&=u(n+1)v(n+1)-u(n)v(n)\\
&=v(n+1)(u(n+1)-u(n))+u(n)(v(n+1)-v(n)),
\end{aligned}\begin{aligned}
\Delta(uv)(n)
&=u(n+1)v(n+1)-u(n)v(n)\\
&=v(n+1)(u(n+1)-u(n))+u(n)(v(n+1)-v(n)),
\end{aligned}
so
\Delta(uv)=(Ev)\Delta u+u\Delta v.
The shift Ev cannot be omitted by copying the continuous product rule.
9Summation by Parts Is Discrete Integration by Parts
Rearranging the product rule gives u\Delta v=\Delta(uv)-(Ev)\Delta u. Indefinite summation yields
\sum u\Delta v=uv-\sum(Ev)\Delta u,
the discrete counterpart of \int u\,dv=uv-\int v\,du. The one-step shift explains why the shift operator must remain explicit.
On a finite range, the precise boundary formula is
\sum_{n=a}^{b-1}u(n)\Delta v(n)
=u(b)v(b)-u(a)v(a)
-\sum_{n=a}^{b-1}v(n+1)\Delta u(n).
10A Recurrence Is P(E)a=b
For
a_{n+k}=c_1a_{n+k-1}+\cdots+c_ka_n+b_n,
define P(t)=t^k-c_1t^{k-1}-\cdots-c_k. Then
P(E)a=b.
The operator equation is solvable exactly when b\in\operatorname{Im}P(E). If a_p is one solution, every solution is
a=a_p+h,\qquad h\in\ker P(E).
11Eigen-Sequences and the Characteristic Equation
For bilateral sequences assume r\ne0; for unilateral sequences define r^n on n\ge0. The geometric sequence satisfies
E(r^n)=r\,r^n,
so it is an eigen-sequence of E with eigenvalue r. Hence
P(E)(r^n)=P(r)r^n.
It solves P(E)a=0 precisely when P(r)=0. The characteristic-equation method is therefore the discrete eigenvalue method for constant-coefficient linear recurrences.
12Example: A Fibonacci-Type Recurrence
The recurrence f_{n+2}=f_{n+1}+f_n is
(E^2-E-I)f=0.
Its characteristic polynomial is P(t)=t^2-t-1, with roots
r_1=\frac{1+\sqrt5}{2},\qquad r_2=\frac{1-\sqrt5}{2}.
The sequences r_1^n and r_2^n are linearly independent. Moreover, a solution of this second-order recurrence is uniquely determined by the two initial values f_0,f_1, so its solution space is two-dimensional. Hence these sequences form a basis and
f_n=C_1r_1^n+C_2r_2^n,
and initial conditions determine C_1,C_2. The geometric trial sequences are justified because they are eigen-sequences of E.
13Constant and Nonconstant Coefficients
Constant coefficients commute with E, allowing the polynomial P(E). For the variable-coefficient recurrence a_{n+1}-na_n=0, define (Ma)_n=na_n. The equation is (E-M)a=0, but
(EMa)_n=(n+1)a_{n+1},\qquad (MEa)_n=na_{n+1},
so EM\ne ME. Thus the ordinary polynomial calculus and simple characteristic equation generally fail for nonconstant coefficients. Constancy is an essential commutativity hypothesis.
14Newton's Forward Formula (n\ge0)
Since E=I+\Delta and the operators commute,
E^n=(I+\Delta)^n
=\sum_{k=0}^n\binom{n}{k}\Delta^k.
Applying both sides to the zeroth entry of f gives
f(n)=(E^nf)_0
=\sum_{k=0}^n\binom{n}{k}(\Delta^kf)_0,
which reconstructs a sequence from its initial forward differences.
15Scope
The representation P(E)a=b applies primarily to constant-coefficient linear recurrences. Nonlinear recurrences lack this linear operator structure; variable coefficients generally do not commute with E, so the usual characteristic method is unavailable.
On unilateral sequence spaces, E^{-1} need not have a natural definition. Any use of backward differences or bilateral shifts must specify the index set of the sequence space.
16Canonical Formulas
\boxed{(Ea)_n=a_{n+1}},\qquad
\boxed{\Delta=E-I},
\boxed{\Delta F=f\quad(\text{indefinite sum problem})},
\boxed{P(E)a=b\quad(\text{constant-coefficient linear recurrence})},
\boxed{E(r^n)=r\,r^n,\qquad P(E)(r^n)=P(r)r^n},
\boxed{P(r)=0\quad(\text{characteristic equation})}.
17Summary
A recurrence relation is a linear operator equation in the shift E when it is linear with constant coefficients. Its characteristic equation arises from the eigen-sequences r^n of E.
18Related Lectures
data/lecture/math/analysis/introduction-to-z-transform.lecture.n.md
data/lecture/math/sequence/sequences-and-recurrences.lecture.n.md
data/lecture/math/sequence/recurrence-types-and-strategies.lecture.n.md
data/lecture/math/sequence/first-order-recurrences.lecture.n.md
data/lecture/math/sequence/difference-equation-basics.lecture.n.md
data/lecture/math/linear-operator/linear-operator-equation-basics.lecture.n.md
data/lecture/math/linear-operator/polynomial-operators-and-eigenvalue-problems.lecture.n.md