16

Произведения Кронекера и операторы vec и vech

[32/100%]
Показать
LaTeX
Задача 16.1
?
(a)

Проверьте, что для любых матриц A\mathbf{A} и B\mathbf{B} размера m×nm \times n и любых матриц C\mathbf{C} и D\mathbf{D} размера p×qp \times q

(A+B)⊗(C+D)=(A⊗C)+(A⊗D)+(B⊗C)+(B⊗D). (\mathbf{A} + \mathbf{B}) \otimes (\mathbf{C} + \mathbf{D}) = (\mathbf{A} \otimes \mathbf{C}) + (\mathbf{A} \otimes \mathbf{D}) + (\mathbf{B} \otimes \mathbf{C}) + (\mathbf{B} \otimes \mathbf{D}).
(b)

Проверьте, что для любых матриц A1,A2,…,Ar\mathbf{A}_{1}, \mathbf{A}_{2}, \ldots , \mathbf{A}_{r} размера m×nm \times n и матриц B1,B2,…,Bs\mathbf{B}_{1}, \mathbf{B}_{2}, \ldots , \mathbf{B}_{s} размера p×qp \times q

(∑i=1rAi)⊗(∑j=1sBj)=∑i=1r∑j=1s(Ai⊗Bj). \left(\sum _{i=1}^{r} \mathbf{A}_{i}\right) \otimes \left(\sum _{j=1}^{s} \mathbf{B}_{j}\right) = \sum _{i=1}^{r} \sum _{j=1}^{s} (\mathbf{A}_{i} \otimes \mathbf{B}_{j}).
Задача 16.2

Покажите, что для любой матрицы A\mathbf{A} размера m×nm \times n и матрицы B\mathbf{B} размера p×qp \times q

A⊗B=(A⊗Ip)diag⁡(B,B,…,B). \mathbf{A} \otimes \mathbf{B} = (\mathbf{A} \otimes \mathbf{I}_{p}) \operatorname {diag}(\mathbf{B}, \mathbf{B}, \ldots , \mathbf{B}).
?
Задача 16.3

Покажите, что для любого вектора a\mathbf{a} размера m×1m \times 1 и любого вектора b\mathbf{b} размера p×1p \times 1 (1) a⊗b=(a⊗Ip)b\mathbf{a} \otimes \mathbf{b} = (\mathbf{a} \otimes \mathbf{I}_{p})\mathbf{b} и (2) a′⊗b′=b′(a′⊗Ip)\mathbf{a}' \otimes \mathbf{b}' = \mathbf{b}'(\mathbf{a}' \otimes \mathbf{I}_{p}).

?
Задача 16.4

Пусть A\mathbf{A} и B\mathbf{B} — квадратные матрицы.

?
(a)

Покажите, что если A\mathbf{A} и B\mathbf{B} ортогональны, то A⊗B\mathbf{A} \otimes \mathbf{B} ортогональна.

(b)

Покажите, что если A\mathbf{A} и B\mathbf{B} идемпотентны, то A⊗B\mathbf{A} \otimes \mathbf{B} идемпотентна.

Задача 16.5

Пусть mm, nn, pp и qq — произвольные положительные целые числа. Покажите

?
(a)

что для любой матрицы B\mathbf{B} размера p×qp \times q (с p>1p > 1 или q>1q > 1) существует такая матрица A\mathbf{A} размера m×nm \times n, что обобщённые обратные матрицы для A⊗B\mathbf{A} \otimes \mathbf{B} не все представимы в виде A−⊗B−\mathbf{A}^{-} \otimes \mathbf{B}^{-};

(b)

что для любой матрицы A\mathbf{A} размера m×nm \times n (с m>1m > 1 или n>1n > 1) существует такая матрица B\mathbf{B} размера p×qp \times q, что обобщённые обратные матрицы для A⊗B\mathbf{A} \otimes \mathbf{B} не все представимы в виде A−⊗B−\mathbf{A}^{-} \otimes \mathbf{B}^{-}.

Задача 16.6

Пусть X=A⊗B\mathbf{X} = \mathbf{A} \otimes \mathbf{B}, где A\mathbf{A} — матрица размера m×nm \times n, а B\mathbf{B} — матрица размера p×qp \times q. Покажите, что PX=PA⊗PB\mathbf{P}_{\mathbf{X}} = \mathbf{P}_{\mathbf{A}} \otimes \mathbf{P}_{\mathbf{B}}.

?
Задача 16.7

Покажите, что произведение Кронекера A⊗B\mathbf{A} \otimes \mathbf{B} матрицы A\mathbf{A} размера m×mm \times m и матрицы B\mathbf{B} размера n×nn \times n является

?
(a)

симметрической неотрицательно определённой, если A\mathbf{A} и B\mathbf{B} обе симметрические неотрицательно определённые или обе симметрические неположительно определённые, и

(b)

симметрической положительно определённой, если A\mathbf{A} и B\mathbf{B} обе симметрические положительно определённые или обе симметрические отрицательно определённые.

Задача 16.8

Пусть A\mathbf{A} и B\mathbf{B} — симметрические матрицы размера m×mm \times m, а C\mathbf{C} и D\mathbf{D} — симметрические матрицы размера n×nn \times n. Используя результат упражнения 7 (или иным способом), покажите, что если A−B\mathbf{A} - \mathbf{B}, C−D\mathbf{C} - \mathbf{D}, B\mathbf{B} и C\mathbf{C} неотрицательно определены, то (A⊗C)−(B⊗D)(\mathbf{A} \otimes \mathbf{C}) - (\mathbf{B} \otimes \mathbf{D}) симметрическая неотрицательно определённая.

?
Задача 16.9

Пусть A\mathbf{A} — матрица размера m×nm \times n, а B\mathbf{B} — матрица размера p×qp \times q. Покажите, что в случае обычной нормы

∥A⊗B∥=∥A∥⋅∥B∥. \left\| \mathbf{A} \otimes \mathbf{B}\right\| = \left\| \mathbf{A}\right\| \cdot \left\| \mathbf{B}\right\| .
?
Задача 16.10

Проверьте, что для разбитой на блоки матрицы размера m×nm \times n

A=(A11A12…A1cA21A22…A2c⋮⋮⋮Ar1Ar2…Arc) \mathbf{A} = \begin{pmatrix} \mathbf{A}_{11} & \mathbf{A}_{12} & \ldots & \mathbf{A}_{1c} \\ \mathbf{A}_{21} & \mathbf{A}_{22} & \ldots & \mathbf{A}_{2c} \\ \vdots & \vdots & & \vdots \\ \mathbf{A}_{r1} & \mathbf{A}_{r2} & \ldots & \mathbf{A}_{rc} \end{pmatrix}

и матрицы B\mathbf{B} размера p×qp \times q

A⊗B=(A11⊗BA12⊗B…A1c⊗BA21⊗BA22⊗B…A2c⊗B⋮⋮⋮Ar1⊗BAr2⊗B…Arc⊗B), \mathbf{A} \otimes \mathbf{B} = \begin{pmatrix} \mathbf{A}_{11} \otimes \mathbf{B} & \mathbf{A}_{12} \otimes \mathbf{B} & \ldots & \mathbf{A}_{1c} \otimes \mathbf{B} \\ \mathbf{A}_{21} \otimes \mathbf{B} & \mathbf{A}_{22} \otimes \mathbf{B} & \ldots & \mathbf{A}_{2c} \otimes \mathbf{B} \\ \vdots & \vdots & & \vdots \\ \mathbf{A}_{r1} \otimes \mathbf{B} & \mathbf{A}_{r2} \otimes \mathbf{B} & \ldots & \mathbf{A}_{rc} \otimes \mathbf{B} \end{pmatrix},

то есть A⊗B\mathbf{A} \otimes \mathbf{B} равна матрице размера mp×nqmp \times nq, полученной заменой каждого блока Aij\mathbf{A}_{ij} матрицы A\mathbf{A} произведением Кронекера Aij\mathbf{A}_{ij} и B\mathbf{B}.

?
Задача 16.11

Покажите, что

?
(a)

если T\mathbf{T} и U\mathbf{U} — обе верхнетреугольные матрицы, то T⊗U\mathbf{T} \otimes \mathbf{U} — верхнетреугольная матрица, и

(b)

если T\mathbf{T} и L\mathbf{L} — обе нижнетреугольные матрицы, то T⊗L\mathbf{T} \otimes \mathbf{L} — нижнетреугольная матрица.

Задача 16.12

Пусть A\mathbf{A} — матрица размера m×mm \times m, а B\mathbf{B} — матрица размера n×nn \times n. Предположим, что A\mathbf{A} и B\mathbf{B} имеют LDU-разложения, скажем A=L1D1U1\mathbf{A} = \mathbf{L}_{1}\mathbf{D}_{1}\mathbf{U}_{1} и B=L2D2U2\mathbf{B} = \mathbf{L}_{2}\mathbf{D}_{2}\mathbf{U}_{2}. Используя результаты упражнения 11, покажите, что A⊗B\mathbf{A} \otimes \mathbf{B} имеет LDU-разложение A⊗B=LDU\mathbf{A} \otimes \mathbf{B} = \mathbf{L}\mathbf{D}\mathbf{U}, где L=L1⊗L2\mathbf{L} = \mathbf{L}_{1} \otimes \mathbf{L}_{2}, D=D1⊗D2\mathbf{D} = \mathbf{D}_{1} \otimes \mathbf{D}_{2} и U=U1⊗U2\mathbf{U} = \mathbf{U}_{1} \otimes \mathbf{U}_{2}.

?
Задача 16.13

Пусть A1,A2,…,Ak\mathbf{A}_{1}, \mathbf{A}_{2}, \ldots , \mathbf{A}_{k} — kk матриц (одинаковых размеров). Покажите, что A1,A2,…,Ak\mathbf{A}_{1}, \mathbf{A}_{2}, \ldots , \mathbf{A}_{k} линейно независимы тогда и только тогда, когда vec⁡(A1),vec⁡(A2),…,vec⁡(Ak)\operatorname {vec}(\mathbf{A}_{1}), \operatorname {vec}(\mathbf{A}_{2}), \ldots , \operatorname {vec}(\mathbf{A}_{k}) линейно независимы.

?
Задача 16.14

Пусть mm — положительное целое число, ei\mathbf{e}_{i} обозначает ii-й столбец матрицы Im\mathbf{I}_{m} (i=1,…,mi = 1, \ldots , m), и (для i,j=1,…,mi, j = 1, \ldots , m) пусть Uij=eiej′\mathbf{U}_{ij} = \mathbf{e}_{i}\mathbf{e}_{j}' (в этом случае Uij\mathbf{U}_{ij} — матрица размера m×mm \times m, у которой элемент на месте ijij равен 11, а остальные m2−1m^{2}-1 элементов равны 00).

?
(a)

Покажите, что

vec⁡(Im)=∑i=1mei⊗ei. \operatorname {vec}(\mathbf{I}_{m}) = \sum _{i=1}^{m} \mathbf{e}_{i} \otimes \mathbf{e}_{i}.
(b)

Покажите, что (для i,j,r,s=1,…,mi, j, r, s = 1, \ldots , m)

vec⁡(Uri)[vec⁡(Usj)]′=Uij⊗Urs. \operatorname {vec}(\mathbf{U}_{ri})[\operatorname {vec}(\mathbf{U}_{sj})]' = \mathbf{U}_{ij} \otimes \mathbf{U}_{rs}.
(c)

Покажите, что

∑i=1m∑j=1mUij⊗Uij=vec⁡(Im)[vec⁡(Im)]′. \sum _{i=1}^{m} \sum _{j=1}^{m} \mathbf{U}_{ij} \otimes \mathbf{U}_{ij} = \operatorname {vec}(\mathbf{I}_{m})[\operatorname {vec}(\mathbf{I}_{m})]'.
Задача 16.15

Пусть A\mathbf{A} — матрица размера n×nn \times n.

?
(a)

Покажите, что если A\mathbf{A} ортогональна, то (vec⁡A)′vec⁡A=n(\operatorname {vec} \mathbf{A})'\operatorname {vec} \mathbf{A} = n.

(b)

Покажите, что если A\mathbf{A} идемпотентна, то [vec⁡(A′)]′vec⁡A=rank⁡(A)[\operatorname {vec}(\mathbf{A}')]'\operatorname {vec} \mathbf{A} = \operatorname {rank}(\mathbf{A}).

Задача 16.16

Покажите, что для любой матрицы A\mathbf{A} размера m×nm \times n, матрицы X\mathbf{X} размера p×np \times n, матрицы B\mathbf{B} размера p×pp \times p и матрицы C\mathbf{C} размера n×mn \times m

tr⁡(AX′BXC)=(vec⁡X)′[(A′C′)⊗B]vec⁡X=(vec⁡X)′[(CA)⊗B′]vec⁡X. \operatorname {tr}(\mathbf{A}\mathbf{X}'\mathbf{B}\mathbf{X}\mathbf{C}) = (\operatorname {vec} \mathbf{X})'[(\mathbf{A}'\mathbf{C}') \otimes \mathbf{B}]\operatorname {vec} \mathbf{X} = (\operatorname {vec} \mathbf{X})'[(\mathbf{C}\mathbf{A}) \otimes \mathbf{B}']\operatorname {vec} \mathbf{X}.
?
Задача 16.17
?
(a)

Пусть V\mathcal{V} — линейное пространство матриц размера m×nm \times n, и пусть gg — функция, сопоставляющая значение A∗B\mathbf{A} * \mathbf{B} каждой паре матриц A\mathbf{A} и B\mathbf{B} из V\mathcal{V}. Обозначим через U\mathcal{U} линейное пространство векторов размера mn×1mn \times 1, определённое как

U={x∈Rmn×1:x=vec⁡(A) для некоторой A∈V}, \mathcal{U} = \left\{ \mathbf{x} \in \mathbb {R}^{mn \times 1} : \mathbf{x} = \operatorname {vec}(\mathbf{A}) \text{ для некоторой } \mathbf{A} \in \mathcal{V}\right\} ,

и пусть x⋅y\mathbf{x} \cdot \mathbf{y} — значение, сопоставляемое каждой паре векторов x\mathbf{x} и y\mathbf{y} из U\mathcal{U} некоторым произвольным скалярным произведением ff. Покажите, что gg является скалярным произведением (для V\mathcal{V}) тогда и только тогда, когда существует такое ff, что (для всех A\mathbf{A} и B\mathbf{B} из V\mathcal{V})

A∗B=vec⁡(A)⋅vec⁡(B). \mathbf{A} * \mathbf{B} = \operatorname {vec}(\mathbf{A}) \cdot \operatorname {vec}(\mathbf{B}).
(b)

Пусть gg — функция, сопоставляющая значение A∗B\mathbf{A} * \mathbf{B} произвольной паре матриц A\mathbf{A} и B\mathbf{B} из Rm×n\mathbb {R}^{m \times n}. Покажите, что gg является скалярным произведением (для Rm×n\mathbb {R}^{m \times n}) тогда и только тогда, когда существует разбитая на блоки симметрическая положительно определённая матрица W\mathbf{W} размера mn×mnmn \times mn

W=(W11W12…W1nW21W22…W2n⋮⋮⋱⋮Wn1Wn2…Wnn) \mathbf{W} = \begin{pmatrix} \mathbf{W}_{11} & \mathbf{W}_{12} & \ldots & \mathbf{W}_{1n} \\ \mathbf{W}_{21} & \mathbf{W}_{22} & \ldots & \mathbf{W}_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ \mathbf{W}_{n1} & \mathbf{W}_{n2} & \ldots & \mathbf{W}_{nn} \end{pmatrix}

(где каждая подматрица имеет размеры m×mm \times m), такая что (для всех A\mathbf{A} и B\mathbf{B} из Rm×n\mathbb {R}^{m \times n})

A∗B=∑i,jai′Wijbj, \mathbf{A} * \mathbf{B} = \sum _{i,j} \mathbf{a}_{i}' \mathbf{W}_{ij} \mathbf{b}_{j},

где a1,a2,…,an\mathbf{a}_{1}, \mathbf{a}_{2}, \ldots , \mathbf{a}_{n} и b1,b2,…,bn\mathbf{b}_{1}, \mathbf{b}_{2}, \ldots , \mathbf{b}_{n} — соответственно первый, второй, , nn-й столбцы A\mathbf{A} и B\mathbf{B}.

(c)

Пусть gg — функция, сопоставляющая значение x′∗y′\mathbf{x}' * \mathbf{y}' произвольной паре (строчных) векторов из R1×n\mathbb {R}^{1 \times n}. Покажите, что gg является скалярным произведением (для R1×n\mathbb {R}^{1 \times n}) тогда и только тогда, когда существует такая симметрическая положительно определённая матрица W\mathbf{W} размера n×nn \times n, что (для любой пары nn-мерных строчных векторов x′\mathbf{x}' и y′\mathbf{y}')

x′∗y′=x′Wy. \mathbf{x}' * \mathbf{y}' = \mathbf{x}'\mathbf{W}\mathbf{y}.
Задача 16.18
?
(a)

Определим (для m≥2m \geq 2) P\mathbf{P} как mn×mnmn \times mn матрицу перестановки такую, что для любой матрицы A\mathbf{A} размера m×nm \times n

(vec⁡A∗r)=Pvec⁡A, \begin{pmatrix} \operatorname {vec} \mathbf{A}_{*} \\ \mathbf{r} \end{pmatrix} = \mathbf{P} \operatorname {vec} \mathbf{A},

где A∗\mathbf{A}_{*} — матрица размера (m−1)×n(m-1) \times n, строки которой являются соответственно первой, , (m−1)(m-1)-й строками A\mathbf{A}, а r′\mathbf{r}' — mm-й строкой A\mathbf{A} [и, следовательно, A=(A∗r′)\mathbf{A} = \begin{pmatrix} \mathbf{A}_{*} \\ \mathbf{r}' \end{pmatrix} и A′=(A∗′,r)\mathbf{A}' = (\mathbf{A}_{*}', \mathbf{r})].

  1. Покажите, что Kmn=(Km−1,n00In)P\mathbf{K}_{mn} = \begin{pmatrix} \mathbf{K}_{m-1,n} & \mathbf{0} \\ \mathbf{0} & \mathbf{I}_{n} \end{pmatrix} \mathbf{P}.

  2. Покажите, что ∣P∣=(−1)(m−1)n(n−1)/2\left|\mathbf{P}\right| = (-1)^{(m-1)n(n-1)/2}.

  3. Покажите, что ∣Kmn∣=(−1)(m−1)n(n−1)/2∣Km−1,n∣\left|\mathbf{K}_{mn}\right| = (-1)^{(m-1)n(n-1)/2} \left|\mathbf{K}_{m-1,n}\right|.

(b)

Покажите, что ∣Kmn∣=(−1)m(m−1)n(n−1)/4\left|\mathbf{K}_{mn}\right| = (-1)^{m(m-1)n(n-1)/4}.

(c)

Покажите, что ∣Kmm∣=(−1)m(m−1)/2\left|\mathbf{K}_{mm}\right| = (-1)^{m(m-1)/2}.

Задача 16.19

Покажите, что для любой матрицы A\mathbf{A} размера m×nm \times n, вектора a\mathbf{a} размера p×1p \times 1 и вектора b\mathbf{b} размера q×1q \times 1

  1. b′⊗A⊗a=Kmp[(ab′)⊗A]\mathbf{b}' \otimes \mathbf{A} \otimes \mathbf{a} = \mathbf{K}_{mp}[(\mathbf{a}\mathbf{b}') \otimes \mathbf{A}];

  2. a⊗A⊗b′=Kpm[A⊗(ab′)]\mathbf{a} \otimes \mathbf{A} \otimes \mathbf{b}' = \mathbf{K}_{pm}[\mathbf{A} \otimes (\mathbf{a}\mathbf{b}')].

?
Задача 16.20

Пусть mm и nn — положительные целые числа, ei\mathbf{e}_{i} обозначает ii-й столбец матрицы Im\mathbf{I}_{m} (i=1,…,mi = 1, \ldots , m), а uj\mathbf{u}_{j} обозначает jj-й столбец матрицы In\mathbf{I}_{n} (j=1,…,nj = 1, \ldots , n). Покажите, что

Kmn=∑j=1nuj′⊗Im⊗uj=∑i=1mei⊗In⊗ei′. \mathbf{K}_{mn} = \sum _{j=1}^{n} \mathbf{u}_{j}' \otimes \mathbf{I}_{m} \otimes \mathbf{u}_{j} = \sum _{i=1}^{m} \mathbf{e}_{i} \otimes \mathbf{I}_{n} \otimes \mathbf{e}_{i}'.
?
Задача 16.21

Пусть mm, nn и pp — положительные целые числа. Используя результат упражнения 20, покажите, что

?
(a)

Kmp,n=Kp,mnKm,np\mathbf{K}_{mp,n} = \mathbf{K}_{p,mn}\mathbf{K}_{m,np};

(b)

Kmp,nKnp,mKmn,p=I\mathbf{K}_{mp,n}\mathbf{K}_{np,m}\mathbf{K}_{mn,p} = \mathbf{I};

(c)

Kn,mp=Knp,mKmn,p\mathbf{K}_{n,mp} = \mathbf{K}_{np,m}\mathbf{K}_{mn,p};

(d)

Kp,mnKm,np=Km,npKp,mn\mathbf{K}_{p,mn}\mathbf{K}_{m,np} = \mathbf{K}_{m,np}\mathbf{K}_{p,mn};

(e)

Knp,mKmn,p=Kmn,pKnp,m\mathbf{K}_{np,m}\mathbf{K}_{mn,p} = \mathbf{K}_{mn,p}\mathbf{K}_{np,m};

(f)

Km,npKmp,n=Kmp,nKm,np\mathbf{K}_{m,np}\mathbf{K}_{mp,n} = \mathbf{K}_{mp,n}\mathbf{K}_{m,np}.

Примечание.
?

Начните с того, что обозначьте через uj\mathbf{u}_{j} jj-й столбец матрицы In\mathbf{I}_{n} и покажите, что Kmp,n=∑j(uj′⊗Ip)⊗(Im⊗uj)\mathbf{K}_{mp,n} = \sum_{j} (\mathbf{u}_{j}' \otimes \mathbf{I}_{p}) \otimes (\mathbf{I}_{m} \otimes \mathbf{u}_{j}), а затем используйте результат о том, что для любой матрицы A\mathbf{A} размера m×nm \times n и матрицы B\mathbf{B} размера p×qp \times q B⊗A=Kpm(A⊗B)Knq\mathbf{B} \otimes \mathbf{A} = \mathbf{K}_{pm}(\mathbf{A} \otimes \mathbf{B})\mathbf{K}_{nq}.

Задача 16.22

Пусть A\mathbf{A} — матрица размера m×nm \times n, и определим B=Kmn(A′⊗A)\mathbf{B} = \mathbf{K}_{mn}(\mathbf{A}' \otimes \mathbf{A}). Покажите

?
(a)

что B\mathbf{B} симметрична,

(b)

что rank⁡(B)=[rank⁡(A)]2\operatorname {rank}(\mathbf{B}) = [\operatorname {rank}(\mathbf{A})]^{2},

(c)

что B2=(AA′)⊗(A′A)\mathbf{B}^{2} = (\mathbf{A}\mathbf{A}') \otimes (\mathbf{A}'\mathbf{A}), и

(d)

что tr⁡(B)=tr⁡(A′A)\operatorname {tr}(\mathbf{B}) = \operatorname {tr}(\mathbf{A}'\mathbf{A}).

Задача 16.23

Покажите, что для любой матрицы A\mathbf{A} размера m×nm \times n и любой матрицы B\mathbf{B} размера p×qp \times q

vec⁡(A⊗B)=(In⊗G)vec⁡A=(H⊗Ip)vec⁡B, \operatorname {vec}(\mathbf{A} \otimes \mathbf{B}) = (\mathbf{I}_{n} \otimes \mathbf{G})\operatorname {vec} \mathbf{A} = (\mathbf{H} \otimes \mathbf{I}_{p})\operatorname {vec} \mathbf{B},

где G=(Kqm⊗Ip)(Im⊗vec⁡B)\mathbf{G} = (\mathbf{K}_{qm} \otimes \mathbf{I}_{p})(\mathbf{I}_{m} \otimes \operatorname {vec} \mathbf{B}) и H=(In⊗Kqm)[vec⁡(A)⊗Iq]\mathbf{H} = (\mathbf{I}_{n} \otimes \mathbf{K}_{qm})[\operatorname {vec}(\mathbf{A}) \otimes \mathbf{I}_{q}].

?
Задача 16.24

Покажите, что для Hn=(Gn′Gn)−1Gn′\mathbf{H}_{n} = (\mathbf{G}_{n}'\mathbf{G}_{n})^{-1}\mathbf{G}_{n}'

GnHnHn′=Hn′. \mathbf{G}_{n}\mathbf{H}_{n}\mathbf{H}_{n}' = \mathbf{H}_{n}'.
?
Задача 16.25

Существует единственная матрица Ln\mathbf{L}_{n} такая, что

vech⁡A=Lnvec⁡A \operatorname {vech} \mathbf{A} = \mathbf{L}_{n} \operatorname {vec} \mathbf{A}

для любой матрицы A\mathbf{A} размера n×nn \times n (симметрической или нет). [Матрица Ln\mathbf{L}_{n} — один из вариантов выбора для матрицы Hn\mathbf{H}_{n}, то есть для левой обратной матрицы Gn\mathbf{G}_{n}. Магнус и Нойдекер (1980) называют её матрицей исключения — эффект умножения слева на Ln\mathbf{L}_{n} вектора vec⁡\operatorname {vec} матрицы A\mathbf{A} размера n×nn \times n состоит в исключении (из vec⁡A\operatorname {vec} \mathbf{A}) «наддиагональных» элементов A\mathbf{A}.]

?
(a)

Выпишите элементы L1\mathbf{L}_{1}, L2\mathbf{L}_{2} и L3\mathbf{L}_{3}.

(b)

Для произвольного положительного целого числа nn опишите Ln\mathbf{L}_{n} через её строки.

Задача 16.26

Пусть A\mathbf{A} — матрица размера n×nn \times n, а b\mathbf{b} — вектор размера n×1n \times 1.

?
(a)

Покажите, что (1/2)[(A⊗b′)+(b′⊗A)]Gn=(A⊗b′)Gn(1/2)[(\mathbf{A} \otimes \mathbf{b}') + (\mathbf{b}' \otimes \mathbf{A})]\mathbf{G}_{n} = (\mathbf{A} \otimes \mathbf{b}')\mathbf{G}_{n}.

(b)

Покажите, что для Hn=(Gn′Gn)−1Gn′\mathbf{H}_{n} = (\mathbf{G}_{n}'\mathbf{G}_{n})^{-1}\mathbf{G}_{n}'

  1. (1/2)Hn[(b⊗A)+(A⊗b)]=Hn(b⊗A)(1/2)\mathbf{H}_{n}[(\mathbf{b} \otimes \mathbf{A}) + (\mathbf{A} \otimes \mathbf{b})] = \mathbf{H}_{n}(\mathbf{b} \otimes \mathbf{A});

  2. (A⊗b′)GnHn=(1/2)[(A⊗b′)+(b′⊗A)](\mathbf{A} \otimes \mathbf{b}')\mathbf{G}_{n}\mathbf{H}_{n} = (1/2)[(\mathbf{A} \otimes \mathbf{b}') + (\mathbf{b}' \otimes \mathbf{A})];

  3. GnHn(b⊗A)=(1/2)[(b⊗A)+(A⊗b)]\mathbf{G}_{n}\mathbf{H}_{n}(\mathbf{b} \otimes \mathbf{A}) = (1/2)[(\mathbf{b} \otimes \mathbf{A}) + (\mathbf{A} \otimes \mathbf{b})].

Задача 16.27

Пусть A={aij}\mathbf{A} = \left\{ a_{ij}\right\} — (возможно, несимметрическая) матрица размера n×nn \times n.

?
(a)

Покажите, что для Hn=(Gn′Gn)−1Gn′\mathbf{H}_{n} = (\mathbf{G}_{n}'\mathbf{G}_{n})^{-1}\mathbf{G}_{n}'

Hnvec⁡A=(1/2)vech⁡(A+A′). \mathbf{H}_{n} \operatorname {vec} \mathbf{A} = (1/2) \operatorname {vech}(\mathbf{A} + \mathbf{A}').
(b)

Покажите, что

Gn′Gnvech⁡A=vech⁡[2A−diag⁡(a11,a22,…,ann)]. \mathbf{G}_{n}'\mathbf{G}_{n} \operatorname {vech} \mathbf{A} = \operatorname {vech}[2\mathbf{A} - \operatorname {diag}(a_{11}, a_{22}, \ldots , a_{nn})].
(c)

Покажите, что

Gn′vec⁡A=vech⁡[A+A′−diag⁡(a11,a22,…,ann)]. \mathbf{G}_{n}' \operatorname {vec} \mathbf{A} = \operatorname {vech}[\mathbf{A} + \mathbf{A}' - \operatorname {diag}(a_{11}, a_{22}, \ldots , a_{nn})].
Задача 16.28

Пусть A\mathbf{A} — квадратная матрица порядка nn. Покажите, что для Hn=(Gn′Gn)−1Gn′\mathbf{H}_{n} = (\mathbf{G}_{n}'\mathbf{G}_{n})^{-1}\mathbf{G}_{n}'

GnHn(A⊗A)Hn′=(A⊗A)Hn′. \mathbf{G}_{n}\mathbf{H}_{n}(\mathbf{A} \otimes \mathbf{A})\mathbf{H}_{n}' = (\mathbf{A} \otimes \mathbf{A})\mathbf{H}_{n}'.
?
Задача 16.29

Покажите, что если матрица A={aij}\mathbf{A} = \left\{ a_{ij}\right\} размера n×nn \times n верхнетреугольная, нижнетреугольная или диагональная, то Hn(A⊗A)Gn\mathbf{H}_{n}(\mathbf{A} \otimes \mathbf{A})\mathbf{G}_{n} соответственно верхнетреугольная, нижнетреугольная или диагональная с диагональными элементами aiiajja_{ii}a_{jj} (i=1,…,ni = 1, \ldots , n; j=i,…,nj = i, \ldots , n).

?
Задача 16.30

Пусть A1,…,Ak\mathbf{A}_{1}, \ldots , \mathbf{A}_{k} и B\mathbf{B} — матрицы размера m×nm \times n, и пусть b=vec⁡B\mathbf{b} = \operatorname {vec} \mathbf{B}.

?
(a)

Покажите, что матричное уравнение ∑i=1kxiAi=B\sum_{i=1}^{k} x_{i}\mathbf{A}_{i} = \mathbf{B} (относительно неизвестных x1,…,xkx_{1}, \ldots , x_{k}) равносильно линейной системе вида Ax=b\mathbf{A}\mathbf{x} = \mathbf{b}, где x=(x1,…,xk)′\mathbf{x} = (x_{1}, \ldots , x_{k})' — вектор неизвестных.

(b)

Покажите, что если A1,…,Ak\mathbf{A}_{1}, \ldots , \mathbf{A}_{k} и B\mathbf{B} симметричны, то матричное уравнение ∑i=1kxiAi=B\sum_{i=1}^{k} x_{i}\mathbf{A}_{i} = \mathbf{B} (относительно неизвестных x1,…,xkx_{1}, \ldots , x_{k}) равносильно линейной системе вида A∗x=b∗\mathbf{A}^{*}\mathbf{x} = \mathbf{b}^{*}, где b∗=vech⁡B\mathbf{b}^{*} = \operatorname {vech} \mathbf{B} и x=(x1,…,xk)′\mathbf{x} = (x_{1}, \ldots , x_{k})' — вектор неизвестных.

Задача 16.31

Пусть F\mathbf{F} — матрица размера p×pp \times p, элементы которой — функции, определённые на множестве SS, от вектора x=(x1,…,xm)′\mathbf{x} = (x_{1}, \ldots , x_{m})' из mm переменных. Покажите, что при k=2,3,…k = 2, 3, \ldots

∂vec⁡(Fk)∂x′=∑s=1k[(Fs−1)′⊗Fk−s]∂vec⁡F∂x′ \frac{\partial \operatorname {vec}(\mathbf{F}^{k})}{\partial \mathbf{x}'} = \sum _{s=1}^{k} [(\mathbf{F}^{s-1})' \otimes \mathbf{F}^{k-s}]\frac{\partial \operatorname {vec} \mathbf{F}}{\partial \mathbf{x}'}

(где F0=Ip\mathbf{F}^{0} = \mathbf{I}_{p}).

?
Задача 16.32

Пусть F={fis}\mathbf{F} = \left\{ f_{is}\right\} и G\mathbf{G} — матрицы размера p×qp \times q и r×sr \times s соответственно, элементы которых — функции, определённые на множестве SS, от вектора x=(x1,…,xm)′\mathbf{x} = (x_{1}, \ldots , x_{m})' из mm переменных.

?
(a)

Покажите, что (для j=1,…,mj = 1, \ldots , m)

∂(F⊗G)∂xj=(F⊗∂G∂xj)+(∂F∂xj⊗G). \frac{\partial (\mathbf{F} \otimes \mathbf{G})}{\partial x_{j}} = \left(\mathbf{F} \otimes \frac{\partial \mathbf{G}}{\partial x_{j}}\right) + \left(\frac{\partial \mathbf{F}}{\partial x_{j}} \otimes \mathbf{G}\right).
(b)

Покажите, что (для j=1,…,mj = 1, \ldots , m)

∂vec⁡(F⊗G)∂xj=(Iq⊗Ksp⊗Ir)[(vec⁡F)⊗∂vec⁡G∂xj+∂vec⁡F∂xj⊗(vec⁡G)]. \frac{\partial \operatorname {vec}(\mathbf{F} \otimes \mathbf{G})}{\partial x_{j}} = (\mathbf{I}_{q} \otimes \mathbf{K}_{sp} \otimes \mathbf{I}_{r})\left[(\operatorname {vec} \mathbf{F}) \otimes \frac{\partial \operatorname {vec} \mathbf{G}}{\partial x_{j}} + \frac{\partial \operatorname {vec} \mathbf{F}}{\partial x_{j}} \otimes (\operatorname {vec} \mathbf{G})\right].
(c)

Покажите, что

∂vec⁡(F⊗G)∂x′=(Iq⊗Ksp⊗Ir)[(vec⁡F)⊗∂vec⁡G∂x′+∂vec⁡F∂x′⊗(vec⁡G)]. \frac{\partial \operatorname {vec}(\mathbf{F} \otimes \mathbf{G})}{\partial \mathbf{x}'} = (\mathbf{I}_{q} \otimes \mathbf{K}_{sp} \otimes \mathbf{I}_{r})\left[(\operatorname {vec} \mathbf{F}) \otimes \frac{\partial \operatorname {vec} \mathbf{G}}{\partial \mathbf{x}'} + \frac{\partial \operatorname {vec} \mathbf{F}}{\partial \mathbf{x}'} \otimes (\operatorname {vec} \mathbf{G})\right].
(d)

Покажите, что в частном случае, когда x′=[(vec⁡X)′,(vec⁡Y)′]\mathbf{x}' = [(\operatorname {vec} \mathbf{X})', (\operatorname {vec} \mathbf{Y})'], F(x)=X\mathbf{F}(\mathbf{x}) = \mathbf{X} и G(x)=Y\mathbf{G}(\mathbf{x}) = \mathbf{Y} для некоторых матриц переменных X\mathbf{X} и Y\mathbf{Y} размера p×qp \times q и r×sr \times s, формула из пункта (c) упрощается до

∂vec⁡(X⊗Y)∂x′=(Iq⊗Ksp⊗Ir)[Ipq⊗(vec⁡Y), (vec⁡X)⊗Irs]. \frac{\partial \operatorname {vec}(\mathbf{X} \otimes \mathbf{Y})}{\partial \mathbf{x}'} = (\mathbf{I}_{q} \otimes \mathbf{K}_{sp} \otimes \mathbf{I}_{r})[\mathbf{I}_{pq} \otimes (\operatorname {vec} \mathbf{Y}),\ (\operatorname {vec} \mathbf{X}) \otimes \mathbf{I}_{rs}].