markdown
Shift Operators and Recurrence Relationsmd 6167471
lecture/math/sequence/shift-operators-and-recurrences.lecture.n.md
Download PDF

Shift Operators and Recurrence Relations

mathsequencerecurrenceshift-operatorlinear-operatorlecture

1Introduction

This lecture unifies differences, indefinite sums, and recurrences through the forward shift operator

(Ea)n=an+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 Z or unilateral sequences indexed by n[PARSE ERROR: Undefined("Command(\"ge\")")]0, as specified. Define

(Ea)n=an+1,(Ia)n=an,

and the forward difference

(Δa)n=an+1-an,Δ=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 Δ=E-I, treating ΔF=f as an inverse problem, expressing recurrences as P(E)a=b, and deriving the characteristic equation from the eigen-sequence rn.

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 α,β,

\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(αa+βb)=αEa+βEb.

5The Difference Is E-I

Since

((E-I)a)n=an+1-an=(Δa)n,

Δ=E-I. Therefore Δ is linear as the difference of linear operators:

Δ(αa+βb)=αΔa+βΔb.

6Higher Differences from the Binomial Theorem

Because E and I commute,

Δm=(E-I)m=k=0m(-1)k(mk)Em-k.

Consequently,

(Δma)n=k=0m(-1)k(mk)an+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 ΔF=f, or F(n+1)-F(n)=f(n). Summing from n=a to b-1 telescopes:

n=ab-1f(n)=n=ab-1(F(n+1)-F(n))=F(b)-F(a).

This is the discrete counterpart of abF(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}

so

Δ(uv)=(Ev)Δu+uΔ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Δv=Δ(uv)-(Ev)Δu. Indefinite summation yields

uΔv=uv-(Ev)Δu,

the discrete counterpart of udv=uv-vdu. The one-step shift explains why the shift operator must remain explicit.

On a finite range, the precise boundary formula is

n=ab-1u(n)Δv(n)=u(b)v(b)-u(a)v(a)-n=ab-1v(n+1)Δu(n).

10A Recurrence Is P(E)a=b

For

an+k=c1an+k-1++ckan+bn,

define P(t)=tk-c1tk-1--ck. Then

P(E)a=b.

The operator equation is solvable exactly when bImP(E). If ap is one solution, every solution is

a=ap+h,hkerP(E).

11Eigen-Sequences and the Characteristic Equation

For bilateral sequences assume r0; for unilateral sequences define rn on n[PARSE ERROR: Undefined("Command(\"ge\")")]0. The geometric sequence satisfies

E(rn)=rrn,

so it is an eigen-sequence of E with eigenvalue r. Hence

P(E)(rn)=P(r)rn.

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 fn+2=fn+1+fn is

(E2-E-I)f=0.

Its characteristic polynomial is P(t)=t2-t-1, with roots

r1=1+52,r2=1-52.

The sequences r1n and r2n are linearly independent. Moreover, a solution of this second-order recurrence is uniquely determined by the two initial values f0,f1, so its solution space is two-dimensional. Hence these sequences form a basis and

fn=C1r1n+C2r2n,

and initial conditions determine C1,C2. 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 an+1-nan=0, define (Ma)n=nan. The equation is (E-M)a=0, but

(EMa)n=(n+1)an+1,(MEa)n=nan+1,

so EMME. 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[PARSE ERROR: Undefined("Command(\"ge\")")]0)

Since E=I+Δ and the operators commute,

En=(I+Δ)n=k=0n(nk)Δk.

Applying both sides to the zeroth entry of f gives

f(n)=(Enf)0=k=0n(nk)(Δ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

[PARSE ERROR: Undefined("Command(\"boxed\")")](Ea)n=an+1,[PARSE ERROR: Undefined("Command(\"boxed\")")]Δ=E-I,
[PARSE ERROR: Undefined("Command(\"boxed\")")]ΔF=f(indefinitesumproblem),
[PARSE ERROR: Undefined("Command(\"boxed\")")]P(E)a=b(constantcoefficientlinearrecurrence)[PARSE ERROR: Undefined("RBrace")],
[PARSE ERROR: Undefined("Command(\"boxed\")")]E(rn)=rrn,P(E)(rn)=P(r)rn,
[PARSE ERROR: Undefined("Command(\"boxed\")")]P(r)=0(characteristicequation).

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 rn 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
raw .n.md をコピー
loc をコピー (filepath:line ~ line)
copy share link
copy encoded share link
path をコピー
copy share link
copy encoded share link
copy share link
copy encoded share link
タブを全て閉じる