Section 1.2

Supplementary Problems

[35/0%]
LaTeX
Problem 1.2.1
?
(i)

Are the four 2×22 \times 2 nonnormal matrices

12(1100),12(001−1),12(1−100),12(0011) \frac{1}{\sqrt{2}} \begin{pmatrix} 1 & 1 \\ 0 & 0 \end{pmatrix}, \quad \frac{1}{\sqrt{2}} \begin{pmatrix} 0 & 0 \\ 1 & -1 \end{pmatrix}, \quad \frac{1}{\sqrt{2}} \begin{pmatrix} 1 & -1 \\ 0 & 0 \end{pmatrix}, \quad \frac{1}{\sqrt{2}} \begin{pmatrix} 0 & 0 \\ 1 & 1 \end{pmatrix}

linearly independent?

(ii)

Show that the nine 3×33 \times 3 matrices

(100000000),(110000000),(111000000),(111001000),(111001001), \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}, \quad \begin{pmatrix} 1 & 1 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}, \quad \begin{pmatrix} 1 & 1 & 1 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}, \quad \begin{pmatrix} 1 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix}, \quad \begin{pmatrix} 1 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 0 & 1 \end{pmatrix}, (111001011),(111001111),(111101111),(111111111) \begin{pmatrix} 1 & 1 & 1 \\ 0 & 0 & 1 \\ 0 & 1 & 1 \end{pmatrix}, \quad \begin{pmatrix} 1 & 1 & 1 \\ 0 & 0 & 1 \\ 1 & 1 & 1 \end{pmatrix}, \quad \begin{pmatrix} 1 & 1 & 1 \\ 1 & 0 & 1 \\ 1 & 1 & 1 \end{pmatrix}, \quad \begin{pmatrix} 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 1 \end{pmatrix}

form a basis in the vector space of the 3×33 \times 3 matrices. This basis is called the spiral basis (and can be extended to any dimension). Which of these matrices are nonnormal?

(iii)

Do the eight 3×33 \times 3 skew-hermitian matrices

Γ1=12(0i0i00000),Γ2=12(010−100000), \Gamma _{1} = \frac{1}{\sqrt{2}} \begin{pmatrix} 0 & i & 0 \\ i & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}, \quad \Gamma _{2} = \frac{1}{\sqrt{2}} \begin{pmatrix} 0 & 1 & 0 \\ -1 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}, Γ3=12(i000−i0000),Γ4=12(00i000i00), \Gamma _{3} = \frac{1}{\sqrt{2}} \begin{pmatrix} i & 0 & 0 \\ 0 & -i & 0 \\ 0 & 0 & 0 \end{pmatrix}, \quad \Gamma _{4} = \frac{1}{\sqrt{2}} \begin{pmatrix} 0 & 0 & i \\ 0 & 0 & 0 \\ i & 0 & 0 \end{pmatrix}, Γ5=12(001000−100),Γ6=12(00000i0i0), \Gamma _{5} = \frac{1}{\sqrt{2}} \begin{pmatrix} 0 & 0 & 1 \\ 0 & 0 & 0 \\ -1 & 0 & 0 \end{pmatrix}, \quad \Gamma _{6} = \frac{1}{\sqrt{2}} \begin{pmatrix} 0 & 0 & 0 \\ 0 & 0 & i \\ 0 & i & 0 \end{pmatrix}, Γ7=12(0000010−10),Γ8=16(i000i000−2i) \Gamma _{7} = \frac{1}{\sqrt{2}} \begin{pmatrix} 0 & 0 & 0 \\ 0 & 0 & 1 \\ 0 & -1 & 0 \end{pmatrix}, \quad \Gamma _{8} = \frac{1}{\sqrt{6}} \begin{pmatrix} i & 0 & 0 \\ 0 & i & 0 \\ 0 & 0 & -2i \end{pmatrix}

together with the 3×33 \times 3 unit matrix form an orthonormal basis in the vector space of 3×33 \times 3 matrices over the complex number. Find all pairwise scalar products ⟨A,B⟩:=tr⁡(AB∗)\left\langle A, B \right\rangle := \operatorname {tr}(AB^{*}). Discuss.

Problem 1.2.2
?
(i)

One can describe a tetrahedron in the vector space R3\mathbb {R}^{3} by specifying vectors v→1,v→2,v→3,v→4\overrightarrow {v}_{1}, \overrightarrow {v}_{2}, \overrightarrow {v}_{3}, \overrightarrow {v}_{4} normal to its faces with lengths equal to the faces' area. Give an example.

(ii)

Consider a tetrahedron defined by the triple of linearly independent vectors v→j∈R3\overrightarrow {v}_{j} \in \mathbb {R}^{3}, j=1,2,3j = 1,2,3. Show that the normal vectors to the faces defined by two of these vectors, normalized to the area of the face, is given by

n→1=12v→2×v→3,n→2=12v→3×v→1,n→3=12v→1×v→2. \overrightarrow {n}_{1} = \frac{1}{2}\overrightarrow {v}_{2} \times \overrightarrow {v}_{3}, \quad \overrightarrow {n}_{2} = \frac{1}{2}\overrightarrow {v}_{3} \times \overrightarrow {v}_{1}, \quad \overrightarrow {n}_{3} = \frac{1}{2}\overrightarrow {v}_{1} \times \overrightarrow {v}_{2}.

Show that

v→1=23Vn→2×n→3,v→2=23Vn→3×n→1,v→3=23Vn→1×n→2 \overrightarrow {v}_{1} = \frac{2}{3V}\overrightarrow {n}_{2} \times \overrightarrow {n}_{3}, \quad \overrightarrow {v}_{2} = \frac{2}{3V}\overrightarrow {n}_{3} \times \overrightarrow {n}_{1}, \quad \overrightarrow {v}_{3} = \frac{2}{3V}\overrightarrow {n}_{1} \times \overrightarrow {n}_{2}

where VV is the volume of the tetrahedron given by

V=13!(v→1×v→2)⋅v→3=29(n→1×n→2)⋅n→3. V = \frac{1}{3!}(\overrightarrow {v}_{1} \times \overrightarrow {v}_{2}) \cdot \overrightarrow {v}_{3} = \sqrt{\frac{2}{9}(\overrightarrow {n}_{1} \times \overrightarrow {n}_{2}) \cdot \overrightarrow {n}_{3}}.

Apply it to normalized vectors which form an orthonormal basis in R3\mathbb {R}^{3}

n→1=12(101),n→2=(010),n→3=12(10−1). \overrightarrow {n}_{1} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix}, \quad \overrightarrow {n}_{2} = \begin{pmatrix} 0 \\ 1 \\ 0 \end{pmatrix}, \quad \overrightarrow {n}_{3} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 0 \\ -1 \end{pmatrix}.
Problem 1.2.3
?
(i)

Show that any rank one positive semidefinite n×nn \times n matrix MM can be written as M=v→v→TM = \overrightarrow {v}\overrightarrow {v}^{T}, where v→\overrightarrow {v} is some (column) vector in Cn\mathbb {C}^{n}.

(ii)

Let u→\overrightarrow {u}, v→\overrightarrow {v} be normalized (column) vectors in Cn\mathbb {C}^{n}. Let AA be an n×nn \times n positive semidefinite matrix over C\mathbb {C}. Show that (u→∗v→)(u→∗Av→)≥0(\overrightarrow {u}^{*}\overrightarrow {v})(\overrightarrow {u}^{*}A\overrightarrow {v}) \geq 0.

(iii)

Show that if hermitian matrices SS and TT are positive semidefinite and commute (ST=TSST = TS), then their product STST is also positive semidefinite. We have to show that

(STv→)∗v→≥0 (ST\overrightarrow {v})^{*}\overrightarrow {v} \geq 0

for all v→∈Cn\overrightarrow {v} \in \mathbb {C}^{n}. Note that if S=0nS = 0_{n}, then obviously (STu)∗u→=0(STu)^{*}\overrightarrow {u} = 0. So assume that S≠0nS \neq 0_{n}, i.e. ∥S∥>0\left\| S\right\| > 0.

Problem 1.2.4

Show that the symmetric matrices

(1111),(101000101),(1001000000001001) \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix}, \quad \begin{pmatrix} 1 & 0 & 1 \\ 0 & 0 & 0 \\ 1 & 0 & 1 \end{pmatrix}, \quad \begin{pmatrix} 1 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 1 \end{pmatrix}

are positive semidefinite. Extend to the n×nn \times n case.

?
Problem 1.2.5

Let AA be a hermitian n×nn \times n matrix and A≠0nA \neq 0_{n}. Show that Am≠0nA^{m} \neq 0_{n} for all m∈Nm \in \mathbb {N}.

?
Problem 1.2.6

Let A∈Rm×nA \in \mathbb {R}^{m \times n} be a nonzero matrix. Let x→∈Rn\overrightarrow {x} \in \mathbb {R}^{n}, y→∈Rm\overrightarrow {y} \in \mathbb {R}^{m} be vectors such that c:=y→TAx→≠0c := \overrightarrow {y}^{T}A\overrightarrow {x} \neq 0. Show that the matrix B:=A−c−1Ax→y→TAB := A-c^{-1}A\overrightarrow {x}\overrightarrow {y}^{T}A has rank exactly one less than the rank of AA.

?
Problem 1.2.7
?
(i)

Consider the two-dimensional Euclidean space and let e→1\overrightarrow {e}_{1}, e→2\overrightarrow {e}_{2} be the standard basis. Consider the seven vectors

v→0=0→,v→1=12e→1+32e→2,v→2=−12e→1−32e→2, \overrightarrow {v}_{0} = \overrightarrow {0}, \quad \overrightarrow {v}_{1} = \frac{1}{2}\overrightarrow {e}_{1}+\frac{\sqrt{3}}{2}\overrightarrow {e}_{2}, \quad \overrightarrow {v}_{2} = -\frac{1}{2}\overrightarrow {e}_{1}-\frac{\sqrt{3}}{2}\overrightarrow {e}_{2}, v→3=−12e→1+32e→2,v→4=12e→1−32e→2,v→5=−e→1,v→6=e→1. \overrightarrow {v}_{3} = -\frac{1}{2}\overrightarrow {e}_{1}+\frac{\sqrt{3}}{2}\overrightarrow {e}_{2}, \quad \overrightarrow {v}_{4} = \frac{1}{2}\overrightarrow {e}_{1}-\frac{\sqrt{3}}{2}\overrightarrow {e}_{2}, \quad \overrightarrow {v}_{5} = -\overrightarrow {e}_{1}, \quad \overrightarrow {v}_{6} = \overrightarrow {e}_{1}.

Find the distance between the vectors and select the vectors pairs with the shortest distance.

(ii)

Consider the vectors in R2\mathbb {R}^{2}

v→k=(cos⁡((k−1)π/4)sin⁡((k−1)π/4)),k=1,3,5,7 \overrightarrow {v}_{k} = \begin{pmatrix} \cos ((k-1)\pi /4) \\ \sin ((k-1)\pi /4) \end{pmatrix}, \quad k = 1,3,5,7

and

v→k=2(cos⁡((k−1)π/4)sin⁡((k−1)π/4)),k=2,4,6,8 \overrightarrow {v}_{k} = \sqrt{2} \begin{pmatrix} \cos ((k-1)\pi /4) \\ \sin ((k-1)\pi /4) \end{pmatrix}, \quad k = 2,4,6,8

which play a role for the lattice Boltzmann model. Find the angles between the vectors.

Problem 1.2.8

Show that the 4×44 \times 4 matrices

Π=12(11/20−1/21/211/2001/211/2−1/201/21),Π~=I4−Π \Pi = \frac{1}{2} \begin{pmatrix} 1 & 1/\sqrt{2} & 0 & -1/\sqrt{2} \\ 1/\sqrt{2} & 1 & 1/\sqrt{2} & 0 \\ 0 & 1/\sqrt{2} & 1 & 1/\sqrt{2} \\ -1/\sqrt{2} & 0 & 1/\sqrt{2} & 1 \end{pmatrix}, \quad \tilde{\Pi } = I_{4}-\Pi

are projection matrices. Describe the subspaces of R4\mathbb {R}^{4} they project into.

?
Problem 1.2.9
?
(i)

The vectors u→\overrightarrow {u}, v→\overrightarrow {v}, w→\overrightarrow {w} point to the vertices of a equilateral triangle

u→=(1/30),v→=(−1/(23)1/2),w→=(−1/(23)−1/2). \overrightarrow {u} = \begin{pmatrix} 1/\sqrt{3} \\ 0 \end{pmatrix}, \quad \overrightarrow {v} = \begin{pmatrix} -1/(2\sqrt{3}) \\ 1/2 \end{pmatrix}, \quad \overrightarrow {w} = \begin{pmatrix} -1/(2\sqrt{3}) \\ -1/2 \end{pmatrix}.

Find the area of this triangle.

(ii)

Let u→\overrightarrow {u}, v→\overrightarrow {v} be (column) vectors in Rn\mathbb {R}^{n}. Show that

A=∣(u→Tu→)(v→Tv→)−(u→Tv→)2∣ A = \sqrt{\left|(\overrightarrow {u}^{T}\overrightarrow {u})(\overrightarrow {v}^{T}\overrightarrow {v})-(\overrightarrow {u}^{T}\overrightarrow {v})^{2}\right|}

calculates the area spanned by the vectors u→\overrightarrow {u} and v→\overrightarrow {v}.

Problem 1.2.10

Assume that two planes in R3\mathbb {R}^{3} given by

kx1+ℓx2+mx3+n=0,k′x1+ℓ′x2+m′x3+n′=0 kx_{1}+\ell x_{2}+mx_{3}+n = 0, \quad k'x_{1}+\ell 'x_{2}+m'x_{3}+n' = 0

be the mirror images with respect to a third plane in R3\mathbb {R}^{3} given by

ax1+bx2+cx3+d=0. ax_{1}+bx_{2}+cx_{3}+d = 0.

Show that (r2=a2+b2+c2r^{2} = a^{2}+b^{2}+c^{2})

(k′ℓ′m′)=1r2(a2−b2−c22ab2ac2ab−a2+b2−c22bc2ac2bc−a2−b2+c2)(kℓm). \begin{pmatrix} k' \\ \ell ' \\ m' \end{pmatrix} = \frac{1}{r^{2}} \begin{pmatrix} a^{2}-b^{2}-c^{2} & 2ab & 2ac \\ 2ab & -a^{2}+b^{2}-c^{2} & 2bc \\ 2ac & 2bc & -a^{2}-b^{2}+c^{2} \end{pmatrix} \begin{pmatrix} k \\ \ell \\ m \end{pmatrix}.
?
Problem 1.2.11
?
(i)

Let AA be a hermitian n×nn \times n matrix over C\mathbb {C} with A2=InA^{2} = I_{n}. Find the matrix (A−1+iIn)−1(A^{-1}+iI_{n})^{-1}.

(ii)

Let BB be an n×nn \times n matrix over C\mathbb {C}. Find the condition on BB such that

(In+iB)(In−iB)=In. (I_{n}+iB)(I_{n}-iB) = I_{n}.

This means that (In−iB)(I_{n}-iB) is the inverse of (In+iB)(I_{n}+iB).

Problem 1.2.12

Consider the vector space Rd\mathbb {R}^{d}. Suppose that {v→j}j=1d\left\{ \overrightarrow {v}_{j}\right\}_{j=1}^{d} and {w→k}k=1d\left\{ \overrightarrow {w}_{k}\right\}_{k=1}^{d} are two bases in Rd\mathbb {R}^{d}. Then there is an invertible d×dd \times d matrix

T=(tjk)j,k=1d T = (t_{jk})_{j,k=1}^{d}

so that

v→j=∑k=1dtjkw→k,j=1,…,d. \overrightarrow {v}_{j} = \sum _{k=1}^{d} t_{jk}\overrightarrow {w}_{k}, \quad j = 1,\ldots ,d.

The bases {v→j}j=1d\left\{ \overrightarrow {v}_{j}\right\}_{j=1}^{d} and {w→k}k=1d\left\{ \overrightarrow {w}_{k}\right\}_{k=1}^{d} are said to have the same orientation if det⁡(T)>0\operatorname {det}\left(T\right) > 0. If det⁡(T)<0\operatorname {det}\left(T\right) < 0, then they have the opposite orientation. Consider the two bases in R2\mathbb {R}^{2}

{v→1=(10),  v→2=(01)},{w→1=12(11),  w→2=12(1−1)}. \left\{ \overrightarrow {v}_{1} = \begin{pmatrix} 1 \\ 0 \end{pmatrix}, \; \overrightarrow {v}_{2} = \begin{pmatrix} 0 \\ 1 \end{pmatrix}\right\} , \quad \left\{ \overrightarrow {w}_{1} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ 1 \end{pmatrix}, \; \overrightarrow {w}_{2} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ -1 \end{pmatrix}\right\} .

Find the orientation.

?
Problem 1.2.13
?
(i)

Let ω>0\omega > 0 (fixed) and t≥0t \geq 0. Show that the 4×44 \times 4 matrix

S(t)=14(1+3e−ωt1−e−ωt1−e−ωt1−e−ωt1−e−ωt1+3e−ωt1−e−ωt1−e−ωt1−e−ωt1−e−ωt1+3e−ωt1−e−ωt1−e−ωt1−e−ωt1−e−ωt1+3e−ωt) S(t) = \frac{1}{4} \begin{pmatrix} 1+3e^{-\omega t} & 1-e^{-\omega t} & 1-e^{-\omega t} & 1-e^{-\omega t} \\ 1-e^{-\omega t} & 1+3e^{-\omega t} & 1-e^{-\omega t} & 1-e^{-\omega t} \\ 1-e^{-\omega t} & 1-e^{-\omega t} & 1+3e^{-\omega t} & 1-e^{-\omega t} \\ 1-e^{-\omega t} & 1-e^{-\omega t} & 1-e^{-\omega t} & 1+3e^{-\omega t} \end{pmatrix}

is a stochastic matrix (probabilistic matrix). Find the eigenvalues of S(t)S(t). Let

π=12(1001). \pi = \frac{1}{2} \begin{pmatrix} 1 & 0 & 0 & 1 \end{pmatrix}.

Find the vectors πS(t)\pi S(t), (πS(t))S(t)(\pi S(t))S(t) etc.. Is S(t1+t2)=S(t1)S(t2)S(t_{1}+t_{2}) = S(t_{1})S(t_{2})?

(ii)

Let ϵ∈[0,1]\epsilon \in [0,1]. Consider the stochastic matrix

S(ϵ)=(ϵ1−ϵ1−ϵϵ). S(\epsilon ) = \begin{pmatrix} \epsilon & 1-\epsilon \\ 1-\epsilon & \epsilon \end{pmatrix}.

Let n=1,2,…n = 1,2,\ldots. Show that

Sn(ϵ)=12(1+(2ϵ−1)n1−(2ϵ−1)n1−(2ϵ−1)n1+(2ϵ−1)n). S^{n}(\epsilon ) = \frac{1}{2} \begin{pmatrix} 1+(2\epsilon -1)^{n} & 1-(2\epsilon -1)^{n} \\ 1-(2\epsilon -1)^{n} & 1+(2\epsilon -1)^{n} \end{pmatrix}.
Problem 1.2.14

Can one find 4×44 \times 4 matrices AA and BB such that

A(0001001001001000)B=(1001000000001001)  ? A \begin{pmatrix} 0 & 0 & 0 & 1 \\ 0 & 0 & 1 & 0 \\ 0 & 1 & 0 & 0 \\ 1 & 0 & 0 & 0 \end{pmatrix} B = \begin{pmatrix} 1 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 1 \end{pmatrix} \; ?

Of course at least one of the matrices AA and BB must be singular. The underlying field F\mathbb {F} could be R\mathbb {R} or char⁡(F)=2\operatorname {char}(\mathbb {F}) = 2.

?
Problem 1.2.15
?
(i)

Find all 2×22 \times 2 matrices such that A2≠02A^{2} \neq 0_{2} but A3=02A^{3} = 0_{2}.

(ii)

Find all non-hermitian 2×22 \times 2 matrices BB such that BB∗=I2BB^{*} = I_{2}.

(iii)

Find all hermitian 2×22 \times 2 matrices HH such that H2=02H^{2} = 0_{2}.

(iv)

Find all 2×22 \times 2 matrices CC over R\mathbb {R} such that CTC+CCT=I2C^{T}C+CC^{T} = I_{2}, C2=02C^{2} = 0_{2}.

(v)

Find all 2×22 \times 2 matrices A1A_{1}, A2A_{2}, A3A_{3} such that A1A2=A2A3A_{1}A_{2} = A_{2}A_{3}, A3A1=A2A3A_{3}A_{1} = A_{2}A_{3}.

(vi)

Find all 2×22 \times 2 matrices AA and BB such that A2=B2=I2A^{2} = B^{2} = I_{2}, AB=−I2AB = -I_{2}.

Problem 1.2.16

Let AA be an m×nm \times n matrix and BB be a p×qp \times q matrix. Then the direct sum of AA and BB, denoted by A⊕BA \oplus B, is the (m+p)×(n+q)(m+p) \times (n+q) matrix defined by

A⊕B:=(A00B). A \oplus B := \begin{pmatrix} A & 0 \\ 0 & B \end{pmatrix}.

Let A1A_{1}, A2A_{2} be m×mm \times m matrices and B1B_{1}, B2B_{2} be n×nn \times n matrices. Show that

(A1⊕B1)(A2⊕B2)=(A1A2)⊕(B1B2). (A_{1} \oplus B_{1})(A_{2} \oplus B_{2}) = (A_{1}A_{2}) \oplus (B_{1}B_{2}).
?
Problem 1.2.17

Consider the symmetric matrix AA and the normalized vector v→\overrightarrow {v}

A=(010101010),v→=13(111). A = \begin{pmatrix} 0 & 1 & 0 \\ 1 & 0 & 1 \\ 0 & 1 & 0 \end{pmatrix}, \quad \overrightarrow {v} = \frac{1}{\sqrt{3}} \begin{pmatrix} 1 \\ 1 \\ 1 \end{pmatrix}.

Calculate Sj=v→∗Ajv→S_{j} = \overrightarrow {v}^{*}A^{j}\overrightarrow {v} for j=1,2,3j = 1,2,3. Can AA be reconstructed from S1,S2,S3S_{1}, S_{2}, S_{3} and the information that AA is real symmetric?

?
Problem 1.2.18

The standard simplex Δn\Delta_{n} is defined by the set in Rn\mathbb {R}^{n}

Δn:={(x1,…,xn)T:xj≥0,  ∑j=1nxj=1}. \Delta _{n} := \left\{ (x_{1},\ldots ,x_{n})^{T} : x_{j} \geq 0, \; \sum _{j=1}^{n} x_{j} = 1\right\} .

Consider nn affinely independent points B1,…,Bn∈ΔnB_{1},\ldots ,B_{n} \in \Delta_{n}. They span an (n−1)(n-1)-simplex denoted by Λ=Con⁡(B1,…,Bn)\Lambda = \operatorname {Con}(B_{1},\ldots ,B_{n}), that is

Λ=Con⁡(B1,…,Bn)={λ1B1+⋯+λnBn:∑j=1nλj=1,  λ1,…,λn≥0}. \Lambda = \operatorname {Con}(B_{1},\ldots ,B_{n}) = \left\{ \lambda _{1}B_{1}+\cdots +\lambda _{n}B_{n} : \sum _{j=1}^{n} \lambda _{j} = 1, \; \lambda _{1},\ldots ,\lambda _{n} \geq 0\right\} .

The set corresponds to an invertible n×nn \times n matrix [Λ][\Lambda ] whose columns are B1,…,BnB_{1},\ldots ,B_{n}. Conversely, consider the matrix C=(bjk)C = (b_{jk}), where Ck=(b1k,…,bnk)TC_{k} = (b_{1k},\ldots ,b_{nk})^{T} (k=1,…,nk = 1,\ldots ,n). If det⁡(C)≠0\operatorname {det}\left(C\right) \neq 0 and the sum of the entries in each column is 1, then the matrix CC corresponds to an (n−1)(n-1)-simplex Con⁡(B1,…,Bn)\operatorname {Con}(B_{1},\ldots ,B_{n}) in Δn\Delta_{n}. Let C1C_{1} and C2C_{2} be n×nn \times n matrices with nonnegative entries and all the columns of each matrix add up to 1.

?
(i)

Show that C1C2C_{1}C_{2} and C2C1C_{2}C_{1} are also such matrices.

(ii)

Are the n2×n2n^{2} \times n^{2} matrices C1⊗C2C_{1} \otimes C_{2}, C2⊗C1C_{2} \otimes C_{1} such matrices?

Problem 1.2.19
?
(i)

Let SS, TT be n×nn \times n matrices over C\mathbb {C} with S2=InS^{2} = I_{n}, (TS)2=In(TS)^{2} = I_{n}. Thus SS and TT are invertible. Show that STS−1=T−1STS^{-1} = T^{-1}, ST−1S=TST^{-1}S = T.

(ii)

Let AA be an invertible n×nn \times n matrix over C\mathbb {C} and BB be an n×nn \times n matrix over C\mathbb {C}. We define the n×nn \times n matrix D:=A−1BAD := A^{-1}BA. Show that Dn=A−1BnAD^{n} = A^{-1}B^{n}A applying AA−1=InAA^{-1} = I_{n}.

Problem 1.2.20

Find all invertible 2×22 \times 2 matrices SS such that

(0100)S(0010)=(0010)S(0100). \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} S \begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix} = \begin{pmatrix} 0 & 0 \\ 1 & 0 \end{pmatrix} S \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}.
?
Problem 1.2.21
?
(i)

Consider the 3×23 \times 2 matrix

A=(111111). A = \begin{pmatrix} 1 & 1 \\ 1 & 1 \\ 1 & 1 \end{pmatrix}.

Can one find 2×32 \times 3 matrices BB such that AB=I2AB = I_{2}?

(ii)

Consider the 3×23 \times 2 matrix

X=(121111). X = \begin{pmatrix} 1 & 2 \\ 1 & 1 \\ 1 & 1 \end{pmatrix}.

Can one find 2×32 \times 3 matrices YY such that XY=I2XY = I_{2}?

Problem 1.2.22

Let x→∈R3\overrightarrow {x} \in \mathbb {R}^{3} and ×\times be the vector product. Find all solutions of

(100)×(x1x2x3)=(001)and(010)×(y1y2y3)=(1/201/2). \begin{pmatrix} 1 \\ 0 \\ 0 \end{pmatrix} \times \begin{pmatrix} x_{1} \\ x_{2} \\ x_{3} \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix} \quad \text{and} \quad \begin{pmatrix} 0 \\ 1 \\ 0 \end{pmatrix} \times \begin{pmatrix} y_{1} \\ y_{2} \\ y_{3} \end{pmatrix} = \begin{pmatrix} 1/\sqrt{2} \\ 0 \\ 1/\sqrt{2} \end{pmatrix}.
?
Problem 1.2.23
?
(i)

Let AA be an 3×33 \times 3 matrix over R\mathbb {R} and u→,v→∈R3\overrightarrow {u}, \overrightarrow {v} \in \mathbb {R}^{3}. Find the conditions on AA, u→\overrightarrow {u}, v→\overrightarrow {v} such that

A(u→×v→)=(Au→)×(Av→). A(\overrightarrow {u} \times \overrightarrow {v}) = (A\overrightarrow {u}) \times (A\overrightarrow {v}).
(ii)

Consider the three vectors v→1,v→2,v→3\overrightarrow {v}_{1}, \overrightarrow {v}_{2}, \overrightarrow {v}_{3} in R3\mathbb {R}^{3}. Show that v→1⋅(v→2×v→3)=0\overrightarrow {v}_{1} \cdot (\overrightarrow {v}_{2} \times \overrightarrow {v}_{3}) = 0 if v→1,v→2,v→3\overrightarrow {v}_{1}, \overrightarrow {v}_{2}, \overrightarrow {v}_{3} are linearly dependent.

(iii)

Let u→,v→∈R3\overrightarrow {u}, \overrightarrow {v} \in \mathbb {R}^{3}. Show that (u→×v→)⋅(u→×v→)≡(u→⋅u→)(v→⋅v→)−(u→⋅v→)2(\overrightarrow {u} \times \overrightarrow {v}) \cdot (\overrightarrow {u} \times \overrightarrow {v}) \equiv (\overrightarrow {u} \cdot \overrightarrow {u})(\overrightarrow {v} \cdot \overrightarrow {v})-(\overrightarrow {u} \cdot \overrightarrow {v})^{2}.

(iv)

Let a→\overrightarrow {a}, b→\overrightarrow {b}, c→\overrightarrow {c}, d→\overrightarrow {d} be vectors in R3\mathbb {R}^{3}. Show that (Lagrange identity)

(a→×b→)⋅(c→×d→)≡det⁡(a→⋅c→b→⋅c→a→⋅d→b→⋅d→). (\overrightarrow {a} \times \overrightarrow {b}) \cdot (\overrightarrow {c} \times \overrightarrow {d}) \equiv \det \begin{pmatrix} \overrightarrow {a} \cdot \overrightarrow {c} & \overrightarrow {b} \cdot \overrightarrow {c} \\ \overrightarrow {a} \cdot \overrightarrow {d} & \overrightarrow {b} \cdot \overrightarrow {d} \end{pmatrix}.
Problem 1.2.24
?
(i)

Consider four nonzero vectors v→1,v→2,v→3,v→4\overrightarrow {v}_{1}, \overrightarrow {v}_{2}, \overrightarrow {v}_{3}, \overrightarrow {v}_{4} in R3\mathbb {R}^{3}. Let

w→:=(v→1×v→2)×(v→3×v→4)≠0→. \overrightarrow {w} := (\overrightarrow {v}_{1} \times \overrightarrow {v}_{2}) \times (\overrightarrow {v}_{3} \times \overrightarrow {v}_{4}) \neq \overrightarrow {0}.

Find w→×(v→1×v→2)\overrightarrow {w} \times (\overrightarrow {v}_{1} \times \overrightarrow {v}_{2}) and w→×(v→3×v→4)\overrightarrow {w} \times (\overrightarrow {v}_{3} \times \overrightarrow {v}_{4}). Discuss. Note that v→1\overrightarrow {v}_{1} and v→2\overrightarrow {v}_{2} span a plane in R3\mathbb {R}^{3} and v→3\overrightarrow {v}_{3} and v→4\overrightarrow {v}_{4} span a plane in R3\mathbb {R}^{3}.

(ii)

Let u→,v→,w→\overrightarrow {u}, \overrightarrow {v}, \overrightarrow {w} be vectors in R3\mathbb {R}^{3}. Show that (Jacobi identity)

u→×(v→×w→)+w→×(u→×v→)+v→×(w→×u→)=0→. \overrightarrow {u} \times (\overrightarrow {v} \times \overrightarrow {w})+\overrightarrow {w} \times (\overrightarrow {u} \times \overrightarrow {v})+\overrightarrow {v} \times (\overrightarrow {w} \times \overrightarrow {u}) = \overrightarrow {0}.
(iii)

Consider the normalized column vectors

v→1=(010)T,v→2=(101)T. \overrightarrow {v}_{1} = \begin{pmatrix} 0 & 1 & 0 \end{pmatrix}^{T}, \quad \overrightarrow {v}_{2} = \begin{pmatrix} 1 & 0 & 1 \end{pmatrix}^{T}.

Do the three vectors v→1,v→2,v→1×v→2\overrightarrow {v}_{1}, \overrightarrow {v}_{2}, \overrightarrow {v}_{1} \times \overrightarrow {v}_{2} form an orthonormal basis in R3\mathbb {R}^{3}?

(iv)

Let v→1,v→2,v→3,v→4\overrightarrow {v}_{1}, \overrightarrow {v}_{2}, \overrightarrow {v}_{3}, \overrightarrow {v}_{4} be column vectors in R3\mathbb {R}^{3}. Show that

(v→1×v→2)T(v→3×v→4)≡(v→1Tv→3)(v→2Tv→4)−(v→2Tv→3)(v→1Tv→4). (\overrightarrow {v}_{1} \times \overrightarrow {v}_{2})^{T}(\overrightarrow {v}_{3} \times \overrightarrow {v}_{4}) \equiv (\overrightarrow {v}_{1}^{T}\overrightarrow {v}_{3})(\overrightarrow {v}_{2}^{T}\overrightarrow {v}_{4})-(\overrightarrow {v}_{2}^{T}\overrightarrow {v}_{3})(\overrightarrow {v}_{1}^{T}\overrightarrow {v}_{4}).

A special case with v→3=v→1\overrightarrow {v}_{3} = \overrightarrow {v}_{1}, v→4=v→2\overrightarrow {v}_{4} = \overrightarrow {v}_{2} is

(v→1×v→2)T(v→1×v→2)≡(v→1Tv→1)(v→2Tv→2)−(v→2Tv→1)(v→1Tv→2). (\overrightarrow {v}_{1} \times \overrightarrow {v}_{2})^{T}(\overrightarrow {v}_{1} \times \overrightarrow {v}_{2}) \equiv (\overrightarrow {v}_{1}^{T}\overrightarrow {v}_{1})(\overrightarrow {v}_{2}^{T}\overrightarrow {v}_{2})-(\overrightarrow {v}_{2}^{T}\overrightarrow {v}_{1})(\overrightarrow {v}_{1}^{T}\overrightarrow {v}_{2}).
Problem 1.2.25

Consider the coordinates

p→1=(x1,y1,z1)T,p→2=(x2,y2,z2)T,p→3=(x3,y3,z3)T \overrightarrow {p}_{1} = (x_{1},y_{1},z_{1})^{T}, \quad \overrightarrow {p}_{2} = (x_{2},y_{2},z_{2})^{T}, \quad \overrightarrow {p}_{3} = (x_{3},y_{3},z_{3})^{T}

with p→1≠p→2\overrightarrow {p}_{1} \neq \overrightarrow {p}_{2}, p→2≠p→3\overrightarrow {p}_{2} \neq \overrightarrow {p}_{3}, p→3≠p→1\overrightarrow {p}_{3} \neq \overrightarrow {p}_{1}. We form the vectors

v→21=(x2−x1y2−y1z2−z1),v→31=(x3−x1y3−y1z3−z1). \overrightarrow {v}_{21} = \begin{pmatrix} x_{2}-x_{1} \\ y_{2}-y_{1} \\ z_{2}-z_{1} \end{pmatrix}, \quad \overrightarrow {v}_{31} = \begin{pmatrix} x_{3}-x_{1} \\ y_{3}-y_{1} \\ z_{3}-z_{1} \end{pmatrix}.

What does 12∣v→21×v→31∣\frac{1}{2}\left|\overrightarrow {v}_{21} \times \overrightarrow {v}_{31}\right| calculate? Apply it to p→1=(0,0,0)T\overrightarrow {p}_{1} = (0,0,0)^{T}, p→2=(1,0,1)T\overrightarrow {p}_{2} = (1,0,1)^{T}, p→3=(1,1,1)T\overrightarrow {p}_{3} = (1,1,1)^{T}.

?
Problem 1.2.26

Consider the 4×44 \times 4 Haar matrix

K=12(111111−1−12−200002−2). K = \frac{1}{2} \begin{pmatrix} 1 & 1 & 1 & 1 \\ 1 & 1 & -1 & -1 \\ \sqrt{2} & -\sqrt{2} & 0 & 0 \\ 0 & 0 & \sqrt{2} & -\sqrt{2} \end{pmatrix}.

Find all 4×44 \times 4 hermitian matrices HH such that KHKT=HKHK^{T} = H.

?
Problem 1.2.27

Consider the real symmetric 3×33 \times 3 matrix

A=(−1/2−3/66/3−3/6−5/6−2/36/3−2/31/3). A = \begin{pmatrix} -1/2 & -\sqrt{3}/6 & \sqrt{6}/3 \\ -\sqrt{3}/6 & -5/6 & -\sqrt{2}/3 \\ \sqrt{6}/3 & -\sqrt{2}/3 & 1/3 \end{pmatrix}.

Show that AT=A−1A^{T} = A^{-1} by showing that the column of the matrix are normalized and pairwise orthonormal.

?
Problem 1.2.28

Let X∈Rn×nX \in \mathbb {R}^{n \times n}. Show that XX can be written as

X=A+S+cIn X = A+S+cI_{n}

where AA is antisymmetric (AT=−AA^{T} = -A), SS is symmetric (ST=SS^{T} = S) with tr⁡(S)=0\operatorname {tr}(S) = 0 and c∈Rc \in \mathbb {R}.

?
Problem 1.2.29

Let cj>0c_{j} > 0 for j=1,…,nj = 1,\ldots ,n. Show that the n×nn \times n matrices

(cjckcj+ck),(1/cj+1/ckcjck) \left(\frac{\sqrt{c_{j}c_{k}}}{c_{j}+c_{k}}\right), \quad \left(\frac{1/c_{j}+1/c_{k}}{\sqrt{c_{j}c_{k}}}\right)

(k=1,…,nk = 1,\ldots ,n) are positive definite.

?
Problem 1.2.30

Let a,b∈Ra,b \in \mathbb {R} and c2:=a2+b2c^{2} := a^{2}+b^{2} with c2>0c^{2} > 0. Consider the 2×22 \times 2 matrix

M(a,b)=1c(abb−a). M(a,b) = \frac{1}{c} \begin{pmatrix} a & b \\ b & -a \end{pmatrix}.

Show that M(a,b)M(a,b)=I2M(a,b)M(a,b) = I_{2}. Thus M(a,b)M(a,b) can be considered as a square root of I2I_{2}.

?
Problem 1.2.31

Show that a nonzero 3×33 \times 3 matrix AA over R\mathbb {R} which satisfies

A2AT+ATA2=2A,AATA=2A,A3=03 A^{2}A^{T}+A^{T}A^{2} = 2A, \quad AA^{T}A = 2A, \quad A^{3} = 0_{3}

is given by

A=2(000100010)⇒AT=2(010001000). A = \sqrt{2} \begin{pmatrix} 0 & 0 & 0 \\ 1 & 0 & 0 \\ 0 & 1 & 0 \end{pmatrix} \Rightarrow A^{T} = \sqrt{2} \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix}.
?
Problem 1.2.32

Let x→\overrightarrow {x} be a normalized column vector in Rn\mathbb {R}^{n}, i.e. x→Tx→=1\overrightarrow {x}^{T}\overrightarrow {x} = 1. A matrix TT is called a Householder matrix if T:=In−2x→x→TT := I_{n}-2\overrightarrow {x}\overrightarrow {x}^{T}. Show that

T2=(In−2x→x→T)(In−2x→x→T)=In−4x→x→T+4x→x→T=In. T^{2} = (I_{n}-2\overrightarrow {x}\overrightarrow {x}^{T})(I_{n}-2\overrightarrow {x}\overrightarrow {x}^{T}) = I_{n}-4\overrightarrow {x}\overrightarrow {x}^{T}+4\overrightarrow {x}\overrightarrow {x}^{T} = I_{n}.
?
Problem 1.2.33

Let AA be an n×nn \times n matrix over R\mathbb {R}. Show that there exists nonnull vectors x→1\overrightarrow {x}_{1}, x→2\overrightarrow {x}_{2} in Rn\mathbb {R}^{n} such that

x→1TAx→1x→1Tx→1≤x→TAx→x→Tx→≤x→2TAx→2x→2Tx→2 \frac{\overrightarrow {x}_{1}^{T}A\overrightarrow {x}_{1}}{\overrightarrow {x}_{1}^{T}\overrightarrow {x}_{1}} \leq \frac{\overrightarrow {x}^{T}A\overrightarrow {x}}{\overrightarrow {x}^{T}\overrightarrow {x}} \leq \frac{\overrightarrow {x}_{2}^{T}A\overrightarrow {x}_{2}}{\overrightarrow {x}_{2}^{T}\overrightarrow {x}_{2}}

for every nonnull vector x→\overrightarrow {x} in Rn\mathbb {R}^{n}.

?
Problem 1.2.34

A generalized Kronecker delta can be defined as follows

δI,J:={1if J=(j1,…,jr) is an even permutation of I=(i1,…,ir)−1if J is an odd permutation of I0if J is not a permutation of I \delta _{I,J} := \begin{cases} 1 & \text{if } J = (j_{1},\ldots ,j_{r}) \text{ is an even permutation of } I = (i_{1},\ldots ,i_{r}) \\ -1 & \text{if } J \text{ is an odd permutation of } I \\ 0 & \text{if } J \text{ is not a permutation of } I \end{cases}

Show that δ126,621=−1\delta_{126,621} = -1, δ126,651=0\delta_{126,651} = 0, δ125,512=1\delta_{125,512} = 1. Give a computer algebra implementation.

?
Problem 1.2.35
?
(i)

Let j0,j1,k0,k1∈{0,1}j_{0},j_{1},k_{0},k_{1} \in \left\{ 0,1\right\}. Consider the tensor

Tk0,k1j0,j1. T^{j_{0},j_{1}}_{k_{0},k_{1}}.

Give a 1−11-1 map that maps Tk0,k1j0,j1T^{j_{0},j_{1}}_{k_{0},k_{1}} to a 22×222^{2} \times 2^{2} matrix S=(sℓ0,ℓ1)S = (s_{\ell_{0},\ell_{1}}) with ℓ0,ℓ1=0,1,2,3\ell_{0},\ell_{1} = 0,1,2,3, i.e.

(t0000t0100t1000t1100t0001t0101t1001t1101t0010t0110t1010t1110t0011t0111t1011t1111)↦(s00s01s02s03s10s11s12s13s20s21s22s23s30s31s32s33). \begin{pmatrix} t^{00}_{00} & t^{00}_{01} & t^{00}_{10} & t^{00}_{11} \\ t^{01}_{00} & t^{01}_{01} & t^{01}_{10} & t^{01}_{11} \\ t^{10}_{00} & t^{10}_{01} & t^{10}_{10} & t^{10}_{11} \\ t^{11}_{00} & t^{11}_{01} & t^{11}_{10} & t^{11}_{11} \end{pmatrix} \mapsto \begin{pmatrix} s_{00} & s_{01} & s_{02} & s_{03} \\ s_{10} & s_{11} & s_{12} & s_{13} \\ s_{20} & s_{21} & s_{22} & s_{23} \\ s_{30} & s_{31} & s_{32} & s_{33} \end{pmatrix}.
(ii)

Let j0,j1,j2,k0,k1,k2∈{0,1}j_{0},j_{1},j_{2},k_{0},k_{1},k_{2} \in \left\{ 0,1\right\}. Consider the tensor

Tk0,k1,k2j0,j1,j2. T^{j_{0},j_{1},j_{2}}_{k_{0},k_{1},k_{2}}.

Give a 1−11-1 map that maps Tk0,k1,k2j0,j1,j2T^{j_{0},j_{1},j_{2}}_{k_{0},k_{1},k_{2}} to a 23×232^{3} \times 2^{3} matrix S=(sℓ0,ℓ1)S = (s_{\ell_{0},\ell_{1}}) with ℓ0,ℓ1=0,1,…,23−1\ell_{0},\ell_{1} = 0,1,\ldots ,2^{3}-1.

(iii)

Let n≥2n \geq 2 and j0,j1,…,jn,k0,k1,…,kn∈{0,1}j_{0},j_{1},\ldots ,j_{n},k_{0},k_{1},\ldots ,k_{n} \in \left\{ 0,1\right\}. Consider the tensor

Tk0,k1,…,kn−1j0,j1,…,jn−1. T^{j_{0},j_{1},\ldots ,j_{n-1}}_{k_{0},k_{1},\ldots ,k_{n-1}}.

Give a 1−11-1 map that maps the tensor to a 2n×2n2^{n} \times 2^{n} matrix S=(sℓ0,ℓ1)S = (s_{\ell_{0},\ell_{1}}) with ℓ0,ℓ1=0,1,…,2n−1\ell_{0},\ell_{1} = 0,1,\ldots ,2^{n}-1. Give a SymbolicC++ implementation. The user provides nn.