Linear Algebra and Optimization

Linear Algebra and Optimization

Chapter One

Problem 2.1.

(1) To prove that (H,)({\cal H},\cdot) is a group we must show that the following properties hold:

  1. Closure: if A,BHA,B \in {\cal H} then ABHA \cdot B \in {\cal H}. To show this note that any upper triangular matrix AA has the property
    [1ab01c001]=[1abd1cef1][111011001] \begin{bmatrix} 1 & a & b \\ 0 & 1 & c \\ 0 & 0 & 1 \end{bmatrix} = \begin{bmatrix} 1 & a & b \\ d & 1 & c \\ e & f & 1 \end{bmatrix} \odot \begin{bmatrix} 1 & 1 & 1 \\ 0 & 1 & 1 \\ 0 & 0 & 1 \end{bmatrix}
    where \odot denotes the Hadamard product of two matrices and the first matrix is non-singular. In index notation we write this as
    Aij=aij1ji A_{ij} = a_{ij} \mathbb{1}_{j \geq i}
    where 1statement=1\mathbb{1}_{\text{statement}}=1 if the statement is true and zero otherwise. Then
    (AB)ij=k=1nAikBkj=k=1naikbkj1ki1jk=k=1naikbkj1ji=(ab)ij1ji \begin{align*} (A \cdot B)_{ij}&=\sum_{k=1}^n A_{ik}B_{kj}=\sum_{k=1}^n a_{ik}b_{kj} \mathbb{1}_{k \geq i}\mathbb{1}_{j \geq k} \\ & =\sum_{k=1}^n a_{ik}b_{kj}\mathbb{1}_{j \geq i} =(a \cdot b)_{ij}\mathbb{1}_{j \geq i} \end{align*}
    and so, ABHA \cdot B \in {\cal H}.
  2. Since H\cal H is a subspace of the vector space of 3×33 \times 3 matrices it inherits associativity and the identity element.
  3. The inverse A1A^{-1} of AHA \in \cal H must have the property
    AA1=A1A=I,A1H. A \cdot A^{-1} = A^{-1} \cdot A = I, \quad A^{-1} \in \cal H.
    We note that
    (AA1)ij=k=1naikakj11ki1jk=k=1nδij1ki1jk=δij (A \cdot A^{-1})_{ij}=\sum_{k=1}^n a_{ik}a^{-1}_{kj} \mathbb{1}_{k \geq i} \mathbb{1}_{j \geq k} = \sum_{k=1}^n \delta_{ij}\mathbb{1}_{k \geq i} \mathbb{1}_{j \geq k} =\delta_{ij}
    where δij=1\delta_{ij}=1 if i=ji=j and zero otherwise. Therefore the inverse exists and is an element of H\cal H. The following explicit form of the inverse
    [1ab01c001] \begin{bmatrix} 1 & a' & b' \\ 0 & 1 & c' \\ 0 & 0 & 1 \end{bmatrix}
    satisfies the diagonal and lower triangular conditions. For the upper triangular elements we must have
    a+a=0b+ac+b=0c+c=0 \begin{align*} a'+a & = 0 \\ b'+ac'+b & = 0 \\ c'+c & = 0 \end{align*}
    and so
    A1=[1aacb01c001]. A^{-1}=\left[\begin{matrix}1 & - a & a c - b\\0 & 1 & - c\\0 & 0 & 1\end{matrix}\right].
    If (H,)({\cal H},\cdot) is abelian then AB=BAA \cdot B = B \cdot A or
    (AB)ij=k=1naikbkj1ki1jk=1jik=1naikbkj=1ji(ab)ij,(BA)ij=k=1nbikakj1ki1jk=1jik=1nbikakj=1ji(ba). \begin{align*} (A\cdot B)_{ij} & = \sum_{k=1}^n a_{ik} b_{kj} \mathbb{1}_{k \geq i} \mathbb{1}_{j \geq k} = \mathbb{1}_{j \geq i}\sum_{k=1}^n a_{ik} b_{kj} = \mathbb{1}_{j \geq i} (a \cdot b)_{ij},\\ (B \cdot A)_{ij} & = \sum_{k=1}^n b_{ik} a_{kj} \mathbb{1}_{k \geq i} \mathbb{1}_{j \geq k} = \mathbb{1}_{j \geq i}\sum_{k=1}^n b_{ik} a_{kj} =\mathbb{1}_{j \geq i} (b \cdot a). \end{align*}
    Since abbaa \cdot b \neq b \cdot a we conclude that H\cal H is not abelian.

(2) If aG1a \in G_1 then from the definition of e1e_1
ϕ(ae1)=ϕ(a)G2 \phi(a e_1 )=\phi(a) \in G_2
and since ϕ\phi is a homomorphism
ϕ(ae1)=ϕ(a)ϕ(e1). \phi(a e_1)=\phi(a) \phi(e_1).
Hence
ϕ(a)ϕ(e1)=ϕ(a) \phi(a) \phi(e_1)=\phi(a)
and we conclude that ϕ(e1)=e2\phi(e_1)=e_2.

Since aa1=e1a a^{-1} = e_1
ϕ(aa1)=ϕ(a)ϕ(a1)=ϕ(e1). \phi(a a^{-1})=\phi(a) \phi(a^{-1})=\phi(e_1).
We have already shown that ϕ(e1)=e2\phi(e_1)=e_2. Therefore ϕ(a1)=ϕ(a)1\phi(a^{-1})=\phi(a)^{-1}.

(3) For any AHA \in \cal H, elements a,b,cRa, b , c \in \mathbb{R}. Since if bRb \in \mathbb{R}, eibe^{ib} generates all elements of S1S^1 we conclude that ϕ\phi is surjective.

To prove that (G,)(G,\cdot) first note that closure is a consequence of the closure of R\mathbb{R} w.r.t. addition and multiplication and the fact that if u1,u2S1u_1,u_2 \in S^1 then, since eix1y2S1e^{ix_1y_2} \in S_1, eix1y2u1u2S1e^{ix_1y_2}u_1 u_2 \in S_1 (closure of S1S^1 under multiplication). To prove associativity write
((x1,y1,u1)(x2,y2,u2))(x3,y3,u3)==(x1+x2,y1+y2,eix1y2u1u2)(x3,y3,u3)==(x1+x2+x3,y1+y2+y3,ei(x1+x2)y3eix1y2u1u2u3)==(x1+(x2+x3),y1+(y2+y3),eix1(y2+y3)u1(eix2y3u2u3))==(x1,y1,u1)(x2+x3,y2+y3,eix2y3u2u3)==(x1,y1,u1)((x2,y2,u2)(x3,y3,u3)). \begin{align*} & ((x_1,y_1,u_1)\cdot(x_2,y_2,u_2))\cdot(x_3,y_3,u_3)=\\ & \kern3ex= (x_1+x_2,y_1+y_2,e^{ix_1y_2}u_1u_2)\cdot(x_3,y_3,u_3)=\\ & \kern3ex= (x_1+x_2+x_3,y_1+y_2+y_3,e^{i(x_1+x_2)y_3}e^{ix_1y_2}u_1u_2u_3)=\\ & \kern3ex= (x_1+(x_2+x_3),y_1+(y_2+y_3),e^{ix_1(y_2+y_3)}u_1(e^{ix_2y_3}u_2u_3))=\\ & \kern3ex= (x_1,y_1,u_1)\cdot(x_2+x_3,y_2+y_3,e^{ix_2y_3}u_2u_3)=\\ & \kern3ex= (x_1,y_1,u_1)\cdot((x_2,y_2,u_2)\cdot(x_3,y_3,u_3)). \end{align*}
The identity element of GG, (0,0,1)(0,0,1) borrows the identity elements of(R,+)(\mathbb{R},+) and (S1,)(S^1,\cdot). The inverse of (x,y,u)(x,y,u) must have the property
(x,y,u)(x,y,u)=(x+x,y+y,eixyuu)=(0,0,1) (x,y,u)\cdot(x',y',u')=(x+x',y+y',e^{ixy'}uu')=(0,0,1)
and so x=x, y=y, u=eixyu1x'=-x,\ y'=-y, \ u'=e^{ixy}u^{-1}. ϕ\phi is a homomorphism if for A,BHA,B \in \cal H
ϕ(AB)=ϕ(A)ϕ(B). \phi(A \cdot B) = \phi(A) \cdot \phi(B).
Since
[1ab01c001][1de01f001]=[1a+daf+b+e01c+f001], \begin{bmatrix} 1 & a & b \\ 0 & 1 & c \\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} 1 & d & e \\ 0 & 1 & f \\ 0 & 0 & 1 \end{bmatrix}= \begin{bmatrix} 1& a+ d & af +b + e \\ 0 & 1 & c+ f \\ 0 & 0 & 1 \end{bmatrix},
we have
ϕ(AB)=(a+d,c+f,ei(af+b+e))=(a,c,eib)(d,f,eie)=ϕ(A)ϕ(B). \begin{align*} \phi(A \cdot B) & =(a+d, c+f,e^{i(af+b+e)}) \\ & =(a,c,e^{ib}) \cdot (d,f,e^{ie}) \\ & = \phi(A) \cdot \phi(B). \end{align*}
(4) We are looking for a set of matrices AHA \in \cal H such that
ϕ(A)=(0,0,1) \phi(A)=(0,0,1)
the identity element in GG. From the definition of ϕ\phi we must have
a=0,c=0,b=2πn a=0, c=0, b =2\pi n
where nZn \in \mathbb{Z}. For any AkerϕA \in \ker\phi we can write A=I+2πne1e3TA=I + 2\pi n e_1 e_3^{\sf T} where e1T=(1,0,0)e_1^{\sf T}=(1,0,0) and e3T=(0,0,1)e_3^{\sf T}=(0,0,1). Closure follows since if A,BkerϕA,B\in\ker\phi
AB=(I+2πne1e3T)(I+2πme1e3T)=I+2π(n+m)e1e3T+(2π)2nme1e3Te1e3T,e3Te1=0=I+2π(n+m)e1e3Tkerϕ. \begin{align*} A \cdot B & = (I + 2 \pi n e_1 e_3^{\sf T})\cdot(I + 2\pi m e_1 e_3^{\sf T})\\ & = I +2 \pi (n+m) e_1 e_3^{\sf T} + (2\pi)^2nm e_1 e_3^{\sf T} e_1 e_3^{\sf T} , \quad e_3^{\sf T} e_1 =0\\ & =I +2 \pi (n+m) e_1 e_3^{\sf T} \in \ker\phi. \end{align*}
Associativity is inherited; the identity element is II (n=0n=0) while the inverse of A=I+2πne1e3TA=I + 2\pi n e_1 e_3^{\sf T} is A1=I2πne1e3TA^{-1}=I - 2\pi n e_1 e_3^{\sf T}. Hence kerϕ\ker\phi is a subgroup of H\cal H.

Problem 2.2.

To check that mZ={mkkZ}m\mathbb{Z}=\{mk \vert k\in \mathbb{Z}\} is an abelian subgroup of Z\mathbb{Z} first note that mk1+mk2=m(k1+k2)mk_1+mk_2=m(k_1+k_2) (closure). Addition is commutative/associative in Z\mathbb{Z} and mZm\mathbb{Z} inherits these properties. It also inherits the identity element (k=0k=0) and the inverse (k1=kk^{-1}=-k).

(1) ϕ(mk):=k\phi(mk):=k; this is an homomorphism since
ϕ(mk1+mk2)=ϕ(m(k1+k2))=k1+k2=ϕ(mk1)+ϕ(mk2). \phi(mk_1+mk_2)=\phi(m(k_1+k_2))=k_1+k_2=\phi(mk_1)+\phi(mk_2).
Since kerϕ={0}\ker\phi=\{0\} it is also an isomorphism.
(2) If i(mk):=mki(mk):=mk then
i(mk1+mk2)=i(m(k1+k2))=m(k1+k2)=mk1+mk2=i(mk1)+i(mk2) i(mk_1+mk_2)=i(m(k_1+k_2))=m(k_1+k_2)=\\mk_1+mk_2= i(mk_1)+i(mk_2)
hence the inclusion map is a group homomorphism. If p:ZmZp:\mathbb{Z} \to m\mathbb{Z} is a homomorphism then
p(l1+l2)=p(l1)+p(l2)=mk1+mk2 p (l_1+l_2)=p(l_1)+p(l_2)=mk_1+m k_2
where l1,l2,k1,k2Zl_1,l_2,k_1,k_2 \in \mathbb{Z} (for example k=l/mk=\lfloor l/m \rfloor). If pi=idp \circ i={\rm id} then pp is the left inverse of ii. By Proposition 2.16 ii must be an isomorphism. However, unless m=1m=1, ii is not an isomorphism since it is not surjective.

Problem 2.3.

As defined EE is an abelian group w.r.t. ++.
From the definition of \cdot we have
λ((x1,y1)+(x2,y2))=λ(x1+x2,y1+y2)=(λ(x1+x2),y1+y2)=(λx1+λx2,y1+y2)=(λx1,y1)+(λx2,y2)=λ(x1,y1)+λ(x2,y2) \lambda\cdot((x_1,y_1)+(x_2,y_2))=\lambda\cdot(x_1+x_2,y_1+y_2) \\ =(\lambda(x_1+x_2),y_1+y_2)=(\lambda x_1+\lambda x_2,y_1+y_2)\\ =(\lambda x_1,y_1)+(\lambda x_2,y_2)=\lambda\cdot(x_1,y_1)+\lambda\cdot(x_2,y_2)
so axiom (V1)\text{(V1)} holds.
(λ+μ)(x,y)=((λ+μ)x,y)=(λx+μx,y)=(λx,y)+(μx,0)=λ(x,y)+μ(x,0)λ(x,y)+μ(x,y) (\lambda + \mu) \cdot (x,y)=((\lambda+\mu)x,y)=(\lambda x + \mu x,y) \\=(\lambda x,y)+(\mu x,0)=\lambda\cdot(x,y)+\mu\cdot(x,0)\neq\lambda\cdot(x,y)+\mu\cdot(x,y)
so axiom (V2)\text{(V2)} fails.
(λμ)(x,y)=((λμ)x,y)=(λ(μx),y)=λ(μx,y)=λ(μ(x,y)) (\lambda \mu)\cdot(x,y)=((\lambda\mu)x,y)=(\lambda(\mu x),y)\\=\lambda\cdot(\mu x,y)=\lambda\cdot(\mu\cdot(x,y))
so axiom (V3)\text{(V3)} holds.
1(x,y)=(1x,y)=(x,y) 1\cdot(x,y)=(1*x,y)=(x,y)
so axiom (V4)\text{(V4)} holds.

Problem 2.4.

(1) Below we (mostly) omit steps that use the properties of a field.
α0=α(uu),any vector u has an additive inverse=αuαu,(V1)=0,the additive inverse of αu is αu \begin{align*} \alpha\cdot 0&=\alpha\cdot(u-u), \quad \text{\footnotesize any vector $u$ has an additive inverse} \\ &=\alpha \cdot u - \alpha \cdot u, \quad \text{\footnotesize (V1)} \\ & =0, \quad \text{\footnotesize the additive inverse of $\alpha \cdot u$ is $-\alpha \cdot u$} \end{align*}
0v=(αα)v,K is a field=αvαv,(V2)=0,the additive inverse of αv is αv \begin{align*} 0 \cdot v &= (\alpha - \alpha) \cdot v, \quad \text{\footnotesize$K$ is a field} \\ & = \alpha \cdot v - \alpha \cdot v, \quad \text{\footnotesize (V2)} \\ & =0, \quad \text{\footnotesize the additive inverse of $\alpha \cdot v$ is $-\alpha \cdot v$} \end{align*}
α(v)=α(1v)=(α1)v,(V3)=(αv) \begin{align*} \alpha \cdot (-v) & = \alpha \cdot (-1 \cdot v) \\ & = (\alpha * -1) \cdot v, \quad \text{\footnotesize (V3)} \\ & = - (\alpha \cdot v) \end{align*}
(α)v=(1α)v=1(αv),(V3)=(αv),the additive inverse of αv is 1(αv) \begin{align*} (-\alpha) \cdot v & = (-1 *\alpha) \cdot v \\ & = -1\cdot (\alpha\cdot v), \quad \text{\footnotesize (V3)} \\ & = - (\alpha \cdot v), \quad \text{\footnotesize the additive inverse of $\alpha \cdot v$ is $-1\cdot(\alpha \cdot v)$} \end{align*}

(2) Axiom (V2)\text{(V2)} states that
α(u+v)=αu+αv \alpha\cdot(u+v)=\alpha\cdot u + \alpha \cdot v
where u,vu,v are vectors and α\alpha is a scalar. Then
αx=α(x1e1+x2e2++xnen)=α(x1e1+(x2e2++xnen))=α(x1e1)+α(x2e2++xnen)=α(x1e1)++α(xn2en2)+α(xn1en1+xnen)=α(x1e1)++α(xnen). \begin{align*} \alpha \cdot x &= \alpha \cdot (x_1 e_1 + x_2 e_2 + \cdots + x_n e_n) \\ & = \alpha \cdot (x_1 e_1 + ( x_2 e_2 + \cdots + x_n e_n)) \\ & = \alpha \cdot (x_1 e_1)+ \alpha \cdot (x_2 e_2 + \cdots + x_n e_n) \\ & \kern15ex \vdots \\ & = \alpha \cdot (x_1 e_1)+\cdots + \alpha \cdot (x_{n-2} e_{n-2})+ \alpha \cdot(x_{n-1}e_{n-1}+ x_n e_n) \\ & = \alpha\cdot (x_1 e_1)+ \cdots + \alpha\cdot(x_n e_n). \end{align*}
Since for any i=1,2,,ni=1,2,\ldots,n,
α(xiei)=α(xiei)=(αxi)ei=(xiα)ei=xi(αei) \alpha\cdot(x_i e_i)=\alpha\cdot(x_i\cdot e_i)=(\alpha * x_i) \cdot e_i \\ =(x_i * \alpha) \cdot e_i=x_i\cdot(\alpha\cdot e_i)
given {αei}i=1n\{\alpha \cdot e_i\}_{i=1}^n we can determine the result of scalar multiplication on any vector by its action on the basis vectors.
(3) Lets work backwards. We want to define scalar multiplication so that
1uu;1 \cdot u \neq u;
if we denote 1u=ν(u)1\cdot u=\nu(u) then if axiom (V3)\text{(V3)} holds
(11)u=1(1u)=1ν(u)=ν2(u). (1*1)\cdot u=1\cdot(1 \cdot u)=1\cdot \nu(u)=\nu^2(u).
Since 11=11*1=1 we conclude that ν(u)=ν2(u)\nu(u)=\nu^2(u) and so, in this setting, scalar multiplication must be a projection. One example is
αx=(αx1,,αxn1,0) \alpha \cdot x=(\alpha x_1,\ldots,\alpha x_{n-1},0)
which for α=1\alpha=1 gives the desired result. All other axioms follow; for example
α(x+y)=(α(x1+y1),,α(xn1+yn1),0)=(αx1,αxn1,0)+(αy1,,αyn1,0)=αx+αy. \begin{align*} \alpha\cdot(x+y)&=(\alpha(x_1+y_1),\ldots,\alpha(x_{n-1}+y_{n-1}),0) \\ & =(\alpha x_1 ,\alpha x_{n-1},0)+(\alpha y_1,\ldots,\alpha y_{n-1},0) \\ & = \alpha \cdot x + \alpha \cdot y. \end{align*}
(4) If nNn \in \mathbb{N} then according to (V2)\text{(V2)}
nu=1u++1u=n(1u). n\cdot u= 1 \cdot u + \ldots + 1 \cdot u = n(1 \cdot u).
Next use (V3)\text{(V3)}:
1u=(1nn)u=1n(nu)=1n(n(1u)). 1\cdot u =\left( \frac{1}{n}* n\right)\cdot u=\frac{1}{n}\cdot(n \cdot u)=\frac{1}{n}\cdot (n(1 \cdot u)).
This is possible only if
1nu=1n(1u). \frac{1}{n}\cdot u=\frac{1}{n}(1 \cdot u).
A rational number r=m/nr=m/n where m,nNm,n \in \mathbb{N}; then
mnu=(m1n)u=m(1nu)=m(1n(1u))=1n(m(1u))=1n(1(mu))=1n(1m(1u))=mn(1(1u))=r((11)u)=r(1u). \frac{m}{n}\cdot u = \left( m * \frac{1}{n} \right) \cdot u =m \cdot \left( \frac{1}{n} \cdot u\right)=m \cdot \left(\frac{1}{n} (1\cdot u) \right) \\=\frac{1}{n}\left(m \cdot (1 \cdot u)\right)=\frac{1}{n}\left(1 \cdot (m \cdot u)\right)=\frac{1}{n}\left( 1 \cdot m (1\cdot u) \right)\\=\frac{m}{n}(1\cdot(1\cdot u))=r((1*1)\cdot u)=r(1 \cdot u).

Written with StackEdit.

Comments

Popular posts from this blog

Poor man's Daletskii-Krein

PAT 2022

BIO 2025