Section 7.11

Minimum Polynomials and Krylov Methods

[12/67%]
Show
LaTeX
Problem 7.11.1

Determine the minimum polynomial for A=(512−40−2−4−1−1)\mathbf{A} = \begin{pmatrix} 5 & 1 & 2 \\ -4 & 0 & -2 \\ -4 & -1 & -1 \end{pmatrix}.

?
Problem 7.11.2

Find the minimum polynomial of b→=(−1,1,1)T\overrightarrow {b} = (-1,1,1)^T with respect to the matrix A\mathbf{A} given in Exercise 7.11.1.

?
Problem 7.11.3

Use Krylov's method to determine the characteristic polynomial for the matrix A\mathbf{A} given in Exercise 7.11.1.

?
Problem 7.11.4

What is the Jordan form for a matrix whose minimum polynomial is m(x)=(x−λ)(x−μ)2m(x) = (x-\lambda )(x-\mu )^2 and whose characteristic polynomial is c(x)=(x−λ)2(x−μ)4c(x) = (x-\lambda )^2(x-\mu )^4?

?
Problem 7.11.5

Use the technique for computing a matrix's minimum polynomial via an orthogonalized Krylov sequence (successively Gram--Schmidt-orthogonalizing I,A,A2,…\mathbf{I}, \mathbf{A}, \mathbf{A}^2, \ldots in the Frobenius inner product ⟨X∣Y⟩=trace⁡(XTY)\langle \mathbf{X} \mid \mathbf{Y} \rangle = \operatorname {trace}(\mathbf{X}^T\mathbf{Y}) until a power Ak\mathbf{A}^k is found to be a linear combination of the earlier orthogonalized terms U0,…,Uk−1\mathbf{U}_0, \ldots , \mathbf{U}_{k-1}, with coefficients recovered by solving Rc→=R−1\mathbf{R}\overrightarrow {c} = \mathbf{R}^{-1} triangular system for the αi\alpha_i's) to determine the minimum polynomial for

A=(−7−48−8−4−14−4−16−817−16−6−36−5). \mathbf{A} = \begin{pmatrix} -7 & -4 & 8 & -8 \\ -4 & -1 & 4 & -4 \\ -16 & -8 & 17 & -16 \\ -6 & -3 & 6 & -5 \end{pmatrix}.
?
Problem 7.11.6

Explain why similar matrices have the same minimum and characteristic polynomials.

?
Problem 7.11.7

Show that two matrices can have the same minimum and characteristic polynomials without being similar by considering A=(N00N)\mathbf{A} = \begin{pmatrix} \mathbf{N} & \mathbf{0} \\ \mathbf{0} & \mathbf{N} \end{pmatrix} and B=(N000)\mathbf{B} = \begin{pmatrix} \mathbf{N} & \mathbf{0} \\ \mathbf{0} & \mathbf{0} \end{pmatrix}, where N=(0100)\mathbf{N} = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}.

?
Problem 7.11.8

Prove that if A\mathbf{A} and B\mathbf{B} are nonderogatory matrices that have the same characteristic polynomial, then A\mathbf{A} is similar to B\mathbf{B}.

?
Problem 7.11.9

Use the Lanczos algorithm to find an orthogonal matrix P\mathbf{P} such that PTAP=T\mathbf{P}^T\mathbf{A}\mathbf{P} = \mathbf{T} is tridiagonal, where A=(211121112)\mathbf{A} = \begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 2 \end{pmatrix}.

?
Problem 7.11.10

Starting with x→0=0→\overrightarrow {x}_0 = \overrightarrow {0}, apply the conjugate gradient algorithm to solve Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b}, where A=(211121112)\mathbf{A} = \begin{pmatrix} 2 & 1 & 1 \\ 1 & 2 & 1 \\ 1 & 1 & 2 \end{pmatrix} and b→=(400)\overrightarrow {b} = \begin{pmatrix} 4 \\ 0 \\ 0 \end{pmatrix}.

?
Problem 7.11.11

Use Arnoldi's algorithm to find an orthogonal matrix Q\mathbf{Q} such that QTAQ=H\mathbf{Q}^T\mathbf{A}\mathbf{Q} = \mathbf{H} is upper Hessenberg, where A=(512−40−2−4−1−1)\mathbf{A} = \begin{pmatrix} 5 & 1 & 2 \\ -4 & 0 & -2 \\ -4 & -1 & -1 \end{pmatrix}.

?
Problem 7.11.12

Use GMRES to solve Ax→=b→\mathbf{A}\overrightarrow {x} = \overrightarrow {b} for A=(512−40−2−4−1−1)\mathbf{A} = \begin{pmatrix} 5 & 1 & 2 \\ -4 & 0 & -2 \\ -4 & -1 & -1 \end{pmatrix} and b→=(121)\overrightarrow {b} = \begin{pmatrix} 1 \\ 2 \\ 1 \end{pmatrix}.

?