Matrix Polynomial Evaluation
You know how to evaluate a polynomial like p(x)=2x2−3x+5 at a number: plug in x, get a number out. Now plug in a square matrix A instead. The variable becomes A, and — crucially — the constant term becomes a multiple of the identity matrix I, because you cannot add a bare number to a matrix.
The Definition
For p(x)=anxn+⋯+a1x+a0 and a square matrix A,
p(A)=anAn+an−1An−1+⋯+a1A+a0I.
Here Ak is k-fold matrix multiplication, akAk is scalar multiplication, and a0I replaces the constant. The result is a square matrix of the same size as A.
There is no ambiguity from non-commutativity here: a polynomial only ever multiplies A by itself, and A always commutes with A.
A Worked Example
Let p(x)=x2−4x+3 and A=(2013).
A2=(4059),−4A=(−80−4−12),3I=(3003).
Adding term by term,
p(A)=(−1010).
A Shortcut for Diagonal Matrices
If A=(λ100λ2), then Ak=(λ1k00λ2k), so
p(A)=(p(λ1)00p(λ2)).
You simply evaluate p at each diagonal entry. …