14

Линейные, билинейные и квадратичные формы

[51/100%]
Показать
LaTeX
Задача 14.1

Покажите, что симметрическую билинейную форму x′Ayx'Ay (от nn-мерных векторов xx и yy) можно выразить через соответствующую квадратичную форму, то есть квадратичную форму, матрицей которой является AA. Сделайте это, проверив, что

x′Ay=(1/2)[(x+y)′A(x+y)−x′Ax−y′Ay]. x'Ay = (1/2)\left[(x+y)'A(x+y) - x'Ax - y'Ay\right] .
?
Задача 14.2

Покажите, что для любой квадратичной формы x′Axx'Ax (от nn-мерного вектора xx) существует единственная верхняя треугольная матрица BB, такая что x′Axx'Ax и x′Bxx'Bx тождественно равны, и выразите элементы BB через элементы AA.

?
Задача 14.3

Покажите на примере, что сумма двух положительно полуопределённых матриц может быть положительно определённой.

?
Задача 14.4

Покажите на примере, что существуют (несимметрические) невырожденные положительно полуопределённые матрицы.

?
Задача 14.5

Покажите на примере, что существуют положительно полуопределённая матрица AA порядка n×nn \times n и матрица PP размера n×mn \times m (где m<nm < n), такие что P′APP'AP положительно определена.

?
Задача 14.6

Для матрицы AA порядка n×nn \times n и матрицы PP размера n×mn \times m имеет место следующее: (1) если AA неотрицательно определена, то P′APP'AP неотрицательно определена; (2) если AA неотрицательно определена и rk⁡(P)<m\operatorname {rk}\left(P\right) < m, то P′APP'AP положительно полуопределена; и (3) если AA положительно определена и rk⁡(P)=m\operatorname {rk}\left(P\right) = m, то P′APP'AP положительно определена. Переформулируйте эти результаты, относящиеся к неотрицательно определённым (положительно определённым или положительно полуопределённым) матрицам, в эквивалентные результаты для неположительно определённых матриц.

?
Задача 14.7

Пусть {X1,…,Xr}\{ X_1, \ldots , X_r\} — набор матриц из линейного пространства V\mathcal{V}. И пусть A={aij}A = \{ a_{ij}\} — матрица порядка r×rr \times r, ijij-й элемент которой равен ⟨Xi,Xj⟩\left\langle X_i, X_j \right\rangle -- эта матрица называется матрицей Грама (или грамианом) набора {X1,…,Xr}\{ X_1, \ldots , X_r\}, а её определитель называется грамианом (или определителем Грама) набора {X1,…,Xr}\{ X_1, \ldots , X_r\}.

?
(a)

Покажите, что AA симметрична и неотрицательно определена.

(b)

Покажите, что X1,…,XrX_1, \ldots , X_r линейно независимы тогда и только тогда, когда AA невырождена.

Задача 14.8

Пусть A={aij}A = \{ a_{ij}\} — симметрическая положительно определённая матрица порядка n×nn \times n, и пусть B={bij}=A−1B = \{ b_{ij}\} = A^{-1}. Покажите, что для i=1,…,ni = 1, \ldots , n

bii≥1/aii, b_{ii} \geq 1/a_{ii} ,

причём равенство имеет место тогда и только тогда, когда для всех j≠ij \neq i выполняется aij=0a_{ij} = 0.

?
Задача 14.9

Пусть AA — матрица размера m×nm \times n, а DD — диагональная матрица, такие что A=PDQA = PDQ для некоторой матрицы PP полного столбцового ранга и некоторой матрицы QQ полного строчного ранга. Покажите, что rk⁡(A)\operatorname {rk}\left(A\right) равен числу ненулевых диагональных элементов DD.

?
Задача 14.10

Пусть AA — симметрическая идемпотентная матрица порядка n×nn \times n, а VV — симметрическая положительно определённая матрица порядка n×nn \times n. Покажите, что rk⁡(AVA)=tr⁡(A)\operatorname {rk}\left(AVA\right) = \operatorname {tr}(A).

?
Задача 14.11

Покажите, что если матрица AA порядка n×nn \times n такова, что x′Ax≠0x'Ax \neq 0 для всякого ненулевого вектора xx размера n×1n \times 1, то AA либо положительно определена, либо отрицательно определена.

?
Задача 14.12
?
(a)

Пусть AA — симметрическая матрица ранга rr порядка n×nn \times n. Возьмём невырожденную матрицу PP порядка n×nn \times n и диагональную матрицу DD порядка n×nn \times n, такие что A=P′DPA = P'DP -- существование таких матриц гарантировано. Число, скажем mm, положительных диагональных элементов DD называется индексом инерции матрицы AA (или квадратичной формы x′Axx'Ax, матрицей которой является AA). Покажите, что индекс инерции определён корректно в том смысле, что mm не зависит от выбора PP или DD. То есть покажите, что если P1P_1 и P2P_2 — невырожденные матрицы, а D1D_1 и D2D_2 — диагональные матрицы, такие что A=P1′D1P1=P2′D2P2A = P_1'D_1P_1 = P_2'D_2P_2, то D2D_2 содержит столько же положительных диагональных элементов, сколько и D1D_1. Покажите также, что число отрицательных диагональных элементов DD равно r−mr - m.

(b)

Пусть AA — симметрическая матрица порядка n×nn \times n. Покажите, что A=P′diag⁡(Im,−Ir−m,0)PA = P'\operatorname {diag}(I_m, -I_{r-m}, 0)P для некоторой невырожденной матрицы PP порядка n×nn \times n и некоторых неотрицательных целых чисел mm и rr. Покажите далее, что mm равно индексу инерции матрицы AA, а r=rk⁡(A)r = \operatorname {rk}\left(A\right).

(c)

Симметричная матрица BB порядка n×nn \times n называется конгруэнтной симметрической матрице AA порядка n×nn \times n, если существует невырожденная матрица PP порядка n×nn \times n, такая что B=P′APB = P'AP. (Если BB конгруэнтна AA, то, очевидно, AA конгруэнтна BB.) Покажите, что BB конгруэнтна AA тогда и только тогда, когда BB имеет тот же ранг и тот же индекс инерции, что и AA. Этот результат называется законом инерции Сильвестра, в честь Джеймса Джозефа Сильвестра (1814--1897).

(d)

Пусть AA — симметрическая матрица ранга rr и с индексом инерции mm порядка n×nn \times n. Покажите, что AA неотрицательно определена тогда и только тогда, когда m=rm = r, и положительно определена тогда и только тогда, когда m=r=nm = r = n.

Задача 14.13

Пусть AA — симметрическая неотрицательно определённая матрица ранга rr порядка n×nn \times n. Тогда существует матрица BB размера n×rn \times r (ранга rr), такая что A=BB′A = BB'. Пусть XX — произвольная матрица размера n×mn \times m (где m≥rm \geq r), такая что A=XX′A = XX'.

?
(a)

Покажите, что X=PBXX = P_BX.

(b)

Покажите, что X=(B,0)QX = (B, 0)Q для некоторой ортогональной матрицы QQ.

Задача 14.14

Покажите, что если симметрическая матрица AA имеет неотрицательно определённую обобщённую обратную, то AA неотрицательно определена.

?
Задача 14.15

Предположим, что матрица AA порядка n×nn \times n имеет LDU-разложение, скажем A=LDUA = LDU, и пусть d1,d2,…,dnd_1, d_2, \ldots , d_n — диагональные элементы диагональной матрицы DD. Покажите, что

∣A∣=d1d2⋯dn. \left|A\right| = d_1d_2 \cdots d_n .
?
Задача 14.16
?
(a)

Предположим, что матрица AA порядка n×nn \times n (где n≥2n \geq 2) имеет единственное LDU-разложение, скажем A=LDUA = LDU, и пусть d1,d2,…,dnd_1, d_2, \ldots , d_n — первый, второй, , nn-й диагональные элементы DD. Покажите, что di≠0d_i \neq 0 (i=1,2,…,n−1)(i = 1, 2, \ldots , n-1) и что dn≠0d_n \neq 0 тогда и только тогда, когда AA невырождена.

(b)

Предположим, что (симметрическая) матрица AA порядка n×nn \times n (где n≥2n \geq 2) имеет единственное разложение U′DUU'DU, скажем A=U′DUA = U'DU, и пусть d1,d2,…,dnd_1, d_2, \ldots , d_n — первый, второй, , nn-й диагональные элементы DD. Покажите, что di≠0d_i \neq 0 (i=1,2,…,n−1)(i = 1, 2, \ldots , n-1) и что dn≠0d_n \neq 0 тогда и только тогда, когда AA невырождена.

Задача 14.17

Предположим, что (симметрическая) матрица AA порядка n×nn \times n имеет единственное разложение U′DUU'DU, скажем A=U′DUA = U'DU. Используя результат части (b) упражнения 14.16, покажите, что AA не имеет LDU-разложений, отличных от A=U′DUA = U'DU.

?
Задача 14.18

Покажите, что если невырожденная матрица имеет LDU-разложение, то это разложение единственно.

?
Задача 14.19

Пусть AA — матрица порядка n×nn \times n (где n≥2n \geq 2). Используя, например, результаты упражнений 14.16, 14.17 и 14.18, покажите, что если AA имеет единственное LDU-разложение либо (в частном случае, когда AA симметрична) единственное разложение U′DUU'DU, то главные угловые подматрицы (матрицы AA) порядков 1,2,…,n−11, 2, \ldots , n-1 невырождены и имеют единственные LDU-разложения.

?
Задача 14.20
?
(a)

Пусть A={aij}A = \{ a_{ij}\} — ненулевая матрица ранга rr размера m×nm \times n. Покажите, что существуют матрица перестановки PP порядка m×mm \times m и матрица перестановки QQ порядка n×nn \times n, такие что

PAQ=(B11B12B21B22), PAQ = \left(\begin{smallmatrix} B_{11} & B_{12} \\ B_{21} & B_{22} \end{smallmatrix}\right) ,

где B11B_{11} -- невырожденная матрица порядка r×rr \times r, главные угловые подматрицы которой (порядков 1,2,…,r−11, 2, \ldots , r-1) невырождены.

(b)

Пусть B=(B11B12B21B22)B = \left(\begin{smallmatrix} B_{11} & B_{12} \\ B_{21} & B_{22} \end{smallmatrix}\right) -- произвольная ненулевая матрица ранга rr размера m×nm \times n, такая что B11B_{11} -- невырожденная матрица порядка r×rr \times r, главные угловые подматрицы которой (порядков 1,2,…,r−11, 2, \ldots , r-1) невырождены. Покажите, что существует единственное разложение BB вида

B=(L1L2)D(U1,U2), B = \left(\begin{smallmatrix} L_1 \\ L_2 \end{smallmatrix}\right) D (U_1, U_2) ,

где L1L_1 -- унитреугольная нижняя матрица порядка r×rr \times r, U1U_1 -- унитреугольная верхняя матрица порядка r×rr \times r, а DD -- диагональная матрица порядка r×rr \times r. Покажите далее, что это разложение таково, что B11=L1DU1B_{11} = L_1DU_1 -- единственное LDU-разложение B11B_{11}, DD невырождена, L2=B21U1−1D−1L_2 = B_{21}U_1^{-1}D^{-1}, а U2=D−1L1−1B12U_2 = D^{-1}L_1^{-1}B_{12}.

Задача 14.21

Покажите на примере, что существуют (несимметрические) положительно полуопределённые матрицы порядка n×nn \times n, не имеющие LDU-разложений.

?
Задача 14.22

Пусть AA — неотрицательно определённая (возможно, несимметрическая) матрица порядка n×nn \times n, имеющая LDU-разложение, скажем A=LDUA = LDU. Покажите, что диагональные элементы диагональной матрицы DD неотрицательны.

?
Задача 14.23

Пусть AA — матрица полного столбцового ранга размера m×km \times k. И пусть A=QRA = QR -- QR-разложение AA; то есть пусть QQ -- единственная матрица размера m×km \times k, столбцы которой ортонормированы относительно обычного скалярного произведения, а RR -- единственная верхняя треугольная матрица порядка k×kk \times k с положительными диагональными элементами, такие что A=QRA = QR. Покажите, что A′A=R′RA'A = R'R (так что A′A=R′RA'A = R'R -- разложение Холецкого матрицы A′AA'A).

?
Задача 14.24

Пусть AA — матрица ранга rr размера m×km \times k (где rr, возможно, меньше kk). Рассмотрим разложение A=QR1A = QR_1, где QQ -- матрица размера m×rm \times r с ортонормированными столбцами, а R1R_1 -- подматрица размера r×kr \times k, строками которой служат rr ненулевых строк верхней треугольной матрицы RR порядка k×kk \times k, имеющей rr положительных диагональных элементов и n−rn - r нулевых строк. (Такое разложение можно получить, используя результаты упражнения 6.4 -- см. упражнение 6.5.) Обобщите результат упражнения 14.23, показав, что если скалярное произведение, относительно которого столбцы QQ ортонормированы, -- это обычное скалярное произведение, то A′A=R′RA'A = R'R (так что A′A=R′RA'A = R'R -- разложение Холецкого матрицы A′AA'A).

?
Задача 14.25

Пусть A={aij}A = \{ a_{ij}\} — матрица порядка n×nn \times n, имеющая LDU-разложение, скажем A=LDUA = LDU. И определим G=U−1D−L−1G = U^{-1}D^{-}L^{-1} (которая является обобщённой обратной матрицы AA).

?
(a)

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

G=D−L−1+(I−U)G=U−1D−+G(I−L). G = D^{-}L^{-1} + (I-U)G = U^{-1}D^{-} + G(I-L) .
(b)

Для i=1,…,ni = 1, \ldots , n пусть did_i -- ii-й диагональный элемент диагональной матрицы DD; и для i,j=1,…,ni, j = 1, \ldots , n пусть ℓij,uij\ell_{ij}, u_{ij} и gijg_{ij} -- ijij-е элементы L,UL, U и GG соответственно. Положим D−=diag⁡(d1∗,…,dn∗)D^{-} = \operatorname {diag}(d_1^{*}, \ldots , d_n^{*}), где di∗=1/did_i^{*} = 1/d_i, если di≠0d_i \neq 0, и di∗d_i^{*} -- произвольный скаляр, если di=0d_i = 0. Покажите, что

gii=di∗−∑k=i+1nuikgki=di∗−∑k=i+1ngikℓki(E.1) g_{ii} = d_i^{*} - \sum _{k=i+1}^{n} u_{ik}g_{ki} = d_i^{*} - \sum _{k=i+1}^{n} g_{ik}\ell _{ki} \tag {E.1}

и что

g_{ij} = \begin{cases} -\sum _{k=j+1}^{n} g_{ik}\ell _{kj}, & \text{при } j < i , \tag {E.2a} \\ -\sum _{k=i+1}^{n} u_{ik}g_{kj}, & \text{при } j > i \tag {E.2b} \end{cases}

(где вырожденные суммы ∑k=n+1ngikℓki\sum_{k=n+1}^{n} g_{ik}\ell_{ki} и ∑k=n+1nuikgki\sum_{k=n+1}^{n} u_{ik}g_{ki} следует понимать как 00).

(c)

Разработайте рекурсивную процедуру, использующую формулы из части (b) для построения обобщённой обратной матрицы AA.

Задача 14.26

Проверьте, что главная подматрица кососимметрической матрицы кососимметрична.

?
Задача 14.27
?
(a)

Покажите, что сумма кососимметрических матриц кососимметрична.

(b)

Покажите, что сумма A1+A2+⋯+AkA_1+A_2+\cdots +A_k неотрицательно определённых матриц A1,A2,…,AkA_1, A_2, \ldots , A_k порядка n×nn \times n кососимметрична тогда и только тогда, когда A1,A2,…,AkA_1, A_2, \ldots , A_k кососимметричны.

(c)

Покажите, что сумма A1+A2+⋯+AkA_1+A_2+\cdots +A_k симметрических неотрицательно определённых матриц A1,A2,…,AkA_1, A_2, \ldots , A_k порядка n×nn \times n является нулевой матрицей тогда и только тогда, когда A1,A2,…,AkA_1, A_2, \ldots , A_k — нулевые матрицы.

Задача 14.28
?
(a)

Пусть A1,A2,…,AkA_1, A_2, \ldots , A_k — неотрицательно определённые матрицы порядка n×nn \times n. Покажите, что tr⁡(∑i=1kAi)≥0\operatorname {tr}(\sum_{i=1}^{k} A_i) \geq 0, причём равенство имеет место тогда и только тогда, когда ∑i=1kAi\sum_{i=1}^{k} A_i кососимметрична, или, что эквивалентно, тогда и только тогда, когда A1,A2,…,AkA_1, A_2, \ldots , A_k кососимметричны.

(b)

Пусть A1,A2,…,AkA_1, A_2, \ldots , A_k — симметрические неотрицательно определённые матрицы порядка n×nn \times n. Покажите, что tr⁡(∑i=1kAi)≥0\operatorname {tr}(\sum_{i=1}^{k} A_i) \geq 0, причём равенство имеет место тогда и только тогда, когда ∑i=1kAi=0\sum_{i=1}^{k} A_i = 0, или, что эквивалентно, тогда и только тогда, когда A1,A2,…,AkA_1, A_2, \ldots , A_k -- нулевые матрицы.

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

То, что кососимметричность ∑i=1kAi\sum_{i=1}^{k} A_i равносильна кососимметричности A1,A2,…,AkA_1, A_2, \ldots , A_k, -- это результат части (b) упражнения 14.27.

Задача 14.29

Покажите на примере, что (при n>1n > 1) существуют (несимметрические) положительно определённые матрицы AA и BB порядка n×nn \times n, такие что tr⁡(AB)<0\operatorname {tr}(AB) < 0.

?
Задача 14.30
?
(a)

Покажите на примере, что (при n>1n > 1) существуют симметрические положительно определённые матрицы AA и BB порядка n×nn \times n, такие что произведение ABAB имеет один или несколько отрицательных диагональных элементов (а значит, ABAB не является неотрицательно определённой).

(b)

Покажите, однако, что произведение двух симметрических положительно определённых матриц порядка n×nn \times n не может быть неположительно определённым.

Задача 14.31

Пусть A={aij}A = \{ a_{ij}\} и B={bij}B = \{ b_{ij}\} — матрицы порядка n×nn \times n, и пусть CC — матрица порядка n×nn \times n, ijij-й элемент которой cij=aijbijc_{ij} = a_{ij}b_{ij} является произведением ijij-х элементов AA и BB -- CC называется произведением Адамара матриц AA и BB. Покажите, что если AA неотрицательно определена, а BB симметрична и неотрицательно определена, то CC неотрицательно определена. Покажите далее, что если AA положительно определена, а BB симметрична и положительно определена, то CC положительно определена.

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

Взяв x=(x1,…,xn)′x = (x_1, \ldots , x_n)' -- произвольный вектор размера n×1n \times 1 -- и F=(f1,…,fn)F = (f_1, \ldots , f_n) -- матрицу, такую что B=F′FB = F'F, начните с того, чтобы показать, что x′Cx=tr⁡(AH)x'Cx = \operatorname {tr}(AH), где H=G′GH = G'G при G=(x1f1,…,xnfn)G = (x_1f_1, \ldots , x_nf_n).

Задача 14.32

Пусть A1,A2,…,AkA_1, A_2, \ldots , A_k и B1,B2,…,BkB_1, B_2, \ldots , B_k — симметрические неотрицательно определённые матрицы порядка n×nn \times n. Покажите, что tr⁡(∑i=1kAiBi)≥0\operatorname {tr}(\sum_{i=1}^{k} A_iB_i) \geq 0, причём равенство имеет место тогда и только тогда, когда для i=1,2,…,ki = 1, 2, \ldots , k выполняется AiBi=0A_iB_i = 0, тем самым обобщив результаты части (b) упражнения 14.28.

?
Задача 14.33

Пусть AA — симметрическая неотрицательно определённая матрица, разбитая на блоки как

A=(TUVW) A = \left(\begin{smallmatrix} T & U \\ V & W \end{smallmatrix}\right)

где TT (а значит, и WW) квадратная. Покажите, что VT−UVT^{-}U и UW−VUW^{-}V симметричны и неотрицательно определены.

?
Задача 14.34

Покажите на примере, что существует (несимметрическая) положительно полуопределённая матрица AA порядка (m+n)×(m+n)(m+n) \times (m+n) вида A=(TUVW)A = \left(\begin{smallmatrix} T & U \\ V & W \end{smallmatrix}\right), где TT имеет размер m×mm \times m, WW -- размер n×nn \times n, UU -- размер m×nm \times n, а VV -- размер n×mn \times m, для которой C(U)⊄C(T)C(U) \not\subset C(T) и/или R(V)⊄R(T)R(V) \not\subset R(T), выражение rk⁡(T)+rk⁡(W−VT−U)\operatorname {rk}\left(T\right) + \operatorname {rk}\left(W-VT^{-}U\right) не обязательно равно rk⁡(A)\operatorname {rk}\left(A\right), а формула

(T−+T−UQ−VT−−T−UQ−−Q−VT−Q−),(*) \left(\begin{smallmatrix} T^{-} + T^{-}UQ^{-}VT^{-} & -T^{-}UQ^{-} \\ -Q^{-}VT^{-} & Q^{-} \end{smallmatrix}\right) , \tag {*}

где Q=W−VT−UQ = W - VT^{-}U, не обязательно даёт обобщённую обратную AA.

?
Задача 14.35

Покажите на примере, что существует симметрическая блочная матрица AA порядка (m+n)×(m+n)(m+n) \times (m+n) вида A=(TUU′W)A = \left(\begin{smallmatrix} T & U \\ U' & W \end{smallmatrix}\right), где TT имеет размер m×mm \times m, UU -- размер m×nm \times n, а WW -- размер n×nn \times n, такая что TT неотрицательно определена и (в зависимости от выбора T−T^{-}) дополнение Шура W−U′T−UW - U'T^{-}U матрицы TT относительно T−T^{-} неотрицательно определено, но AA не является неотрицательно определённой.

?
Задача 14.36

Матрица A={aij}A = \{ a_{ij}\} порядка n×nn \times n называется матрицей с диагональным преобладанием, если для i=1,2,…,ni = 1, 2, \ldots , n выполняется ∣aii∣>∑j=1(j≠i)n∣aij∣\left|a_{ii}\right| > \sum_{j=1(j \neq i)}^{n} \left|a_{ij}\right|. (В вырожденном частном случае n=1n = 1 матрица AA называется матрицей с диагональным преобладанием, если она ненулевая.)

?
(a)

Покажите, что главная подматрица матрицы с диагональным преобладанием сама обладает диагональным преобладанием.

(b)

Пусть A={aij}A = \{ a_{ij}\} — матрица порядка n×nn \times n с диагональным преобладанием, разобьём AA на блоки как A=(A11ab′ann)A = \left(\begin{smallmatrix} A_{11} & a \\ b' & a_{nn} \end{smallmatrix}\right) [так что A11A_{11} имеет размер (n−1)×(n−1)(n-1) \times (n-1)], и пусть C=A11−(1/ann)ab′C = A_{11} - (1/a_{nn})ab' -- дополнение Шура элемента anna_{nn}. Покажите, что CC обладает диагональным преобладанием.

(c)

Покажите, что матрица с диагональным преобладанием невырождена.

(d)

Покажите, что матрица с диагональным преобладанием имеет единственное LDU-разложение.

(e)

Пусть A={aij}A = \{ a_{ij}\} — симметрическая матрица порядка n×nn \times n. Покажите, что если AA обладает диагональным преобладанием и если диагональные элементы a11,a22,…,anna_{11}, a_{22}, \ldots , a_{nn} матрицы AA все положительны, то AA положительно определена.

Задача 14.37

Пусть A={aij}A = \{ a_{ij}\} — симметрическая положительно определённая матрица порядка n×nn \times n. Покажите, что det⁡(A)≤∏i=1naii\det (A) \leq \prod_{i=1}^{n} a_{ii}, причём равенство имеет место тогда и только тогда, когда AA диагональна.

?
Задача 14.38

Пусть A=(abcd)A = \left(\begin{smallmatrix} a & b \\ c & d \end{smallmatrix}\right), где a,b,ca, b, c и dd — скаляры.

?
(a)

Покажите, что AA положительно определена тогда и только тогда, когда a>0a > 0, d>0d > 0 и ∣b+c∣/2<ad\left|b+c\right|/2 < \sqrt{ad}.

(b)

Покажите, что в частном случае, когда AA симметрична (то есть c=bc = b), AA положительно определена тогда и только тогда, когда a>0a > 0, d>0d > 0 и ∣b∣<ad\left|b\right| < \sqrt{ad}.

Задача 14.39

Используя, например, результат упражнения 14.38, покажите, что если матрица A={aij}A = \{ a_{ij}\} порядка n×nn \times n симметрична и положительно определена, то для j≠i=1,…,nj \neq i = 1, \ldots , n

∣aij∣<aiiajj≤max⁡(aii,ajj). \left|a_{ij}\right| < \sqrt{a_{ii}a_{jj}} \leq \max (a_{ii}, a_{jj}) .
?
Задача 14.40

Покажите на примере, что определители обеих главных угловых подматриц симметрической матрицы 2×22 \times 2 могут быть неотрицательны без того, чтобы сама матрица была неотрицательно определена, и что при n≥3n \geq 3 определители всех nn главных угловых подматриц симметрической матрицы порядка n×nn \times n могут быть неотрицательны, а сама матрица -- невырожденной, без того, чтобы матрица была неотрицательно определена.

?
Задача 14.41

Пусть V\mathcal{V} — подпространство Rn×1\mathbb {R}^{n \times 1} размерности rr (где r≥1r \geq 1). Возьмём B=(b1,b2,…,br)B = (b_1, b_2, \ldots , b_r) -- произвольную матрицу размера n×rn \times r, столбцы которой b1,b2,…,brb_1, b_2, \ldots , b_r образуют базис V\mathcal{V}, и пусть LL -- произвольная левая обратная BB. Пусть gg -- функция, сопоставляющая значение x∗yx * y произвольной паре векторов xx и yy из V\mathcal{V}.

?
(a)

Пусть ff — произвольное скалярное произведение на Rr×1\mathbb {R}^{r \times 1}, и обозначим через s⋅ts \cdot t значение, сопоставляемое ff произвольной паре rr-мерных векторов ss и tt. Покажите, что gg является скалярным произведением (на V\mathcal{V}) тогда и только тогда, когда существует такое ff, что (для всех xx и yy из V\mathcal{V})

x∗y=(Lx)⋅(Ly). x * y = (Lx) \cdot (Ly) .
(b)

Покажите, что gg является скалярным произведением (на V\mathcal{V}) тогда и только тогда, когда существует симметрическая положительно определённая матрица WW порядка r×rr \times r, такая что (для всех xx и yy из V\mathcal{V})

x∗y=x′L′WLy. x * y = x'L'WLy .
(c)

Покажите, что gg является скалярным произведением (на V\mathcal{V}) тогда и только тогда, когда существует симметрическая положительно определённая матрица WW порядка n×nn \times n, такая что (для всех xx и yy из V\mathcal{V})

x∗y=x′Wy. x * y = x'Wy .
Задача 14.42

Пусть V\mathcal{V} — линейное пространство матриц размера m×nm \times n, и пусть A⋅BA \cdot B -- значение, сопоставляемое квазискалярным произведением произвольной паре матриц AA и BB (из V\mathcal{V}). Покажите, что множество

U={A∈V:A⋅A=0}, \mathcal{U} = \{ A \in \mathcal{V} : A \cdot A = 0\} ,

состоящее из всех матриц в V\mathcal{V} с нулевой квазинормой, является линейным пространством.

?
Задача 14.43

Пусть WW — симметрическая положительно определённая матрица порядка m×mm \times m, а VV -- симметрическая положительно определённая матрица порядка n×nn \times n.

?
(a)

Покажите, что функция, сопоставляющая значение tr⁡(A′WBV)\operatorname {tr}(A'WBV) произвольной паре матриц AA и BB размера m×nm \times n, является скалярным произведением для линейного пространства Rm×n\mathbb {R}^{m \times n}.

(b)

Покажите, что функция, сопоставляющая значение tr⁡(A′WB)\operatorname {tr}(A'WB) произвольной паре матриц AA и BB размера m×nm \times n, является скалярным произведением на Rm×n\mathbb {R}^{m \times n}.

(c)

Покажите, что функция, сопоставляющая значение tr⁡(A′WBW)\operatorname {tr}(A'WBW) произвольной паре матриц AA и BB размера m×mm \times m, является скалярным произведением на Rm×m\mathbb {R}^{m \times m}.

Задача 14.44

Пусть AA — матрица размера q×pq \times p, BB -- матрица размера p×np \times n, а CC -- матрица размера m×qm \times q. Покажите, что

?
(a)

CAB(CAB)−C=CCAB(CAB)^{-}C = C тогда и только тогда, когда rk⁡(CAB)=rk⁡(C)\operatorname {rk}\left(CAB\right) = \operatorname {rk}\left(C\right), и

(b)

B(CAB)−CAB=BB(CAB)^{-}CAB = B тогда и только тогда, когда rk⁡(CAB)=rk⁡(B)\operatorname {rk}\left(CAB\right) = \operatorname {rk}\left(B\right).

Задача 14.45

Пусть U\mathcal{U} — подпространство Rn×1\mathbb {R}^{n \times 1}, пусть XX -- матрица размера n×pn \times p, столбцы которой порождают U\mathcal{U}, и пусть WW и VV -- симметрические положительно определённые матрицы порядка n×nn \times n. Покажите, что каждое из следующих двух условий необходимо и достаточно для того, чтобы проекция PX,WyP_{X,W}y вектора yy на U\mathcal{U} относительно WW совпадала (для всякого yy из Rn\mathbb {R}^{n}) с проекцией PX,VyP_{X,V}y вектора yy на U\mathcal{U} относительно VV:

?
(a)

V=PX,W′VPX,W+(I−PX,W)′V(I−PX,W)V = P_{X,W}'VP_{X,W} + (I-P_{X,W})'V(I-P_{X,W});

(b)

существуют скаляр cc, матрица KK порядка p×pp \times p и матрица HH порядка n×nn \times n, такие что

V=cW+WXKX′W+(I−PX,W)′H(I−PX,W). V = cW + WXKX'W + (I-P_{X,W})'H(I-P_{X,W}) .
Задача 14.46

Пусть yy — nn-мерный вектор-столбец, пусть U\mathcal{U} -- подпространство Rn×1\mathbb {R}^{n \times 1}, и пусть XX -- матрица размера n×pn \times p, столбцы которой порождают U\mathcal{U}. ``Напомним'', что для любой матрицы LL размера m×nm \times n, для любого подпространства W\mathcal{W} пространства Rn×1\mathbb {R}^{n \times 1} и для V={v∈Rm:v=Lx при некотором x∈W}\mathcal{V} = \{ v \in \mathbb {R}^{m} : v = Lx \text{ при некотором } x \in \mathcal{W}\}

x⊥HW  ⟺  Lx⊥IV,(*) x \perp _H \mathcal{W} \iff Lx \perp _I \mathcal{V} , \tag {*}

где H=L′LH = L'L (и где xx -- произвольный вектор размера n×1n \times 1), и

C(LZ)=V,(**) C(LZ) = \mathcal{V} , \tag {**}

где ZZ -- произвольная матрица размера n×qn \times q, столбцы которой порождают W\mathcal{W}. Используя результаты (*) и (**), покажите, что для любой симметрической неотрицательно определённой матрицы WW порядка n×nn \times n y⊥WUy \perp_W \mathcal{U} тогда и только тогда, когда X′Wy=0X'Wy = 0.

?
Задача 14.47

Пусть WW — симметрическая неотрицательно определённая матрица порядка n×nn \times n.

?
(a)

Покажите, что для любой матрицы XX размера n×pn \times p и любой матрицы UU размера n×qn \times q, такой что C(U)⊂C(X)C(U) \subset C(X),

  1. WPX,WU=WUWP_{X,W}U = WU, и U′WPX,W=U′WU'WP_{X,W} = U'W;

  2. PU,WPX,W=PU,WP_{U,W}P_{X,W} = P_{U,W}, и PX,W′WPU,W=WPX,WPU,W=WPU,WP_{X,W}'WP_{U,W} = WP_{X,W}P_{U,W} = WP_{U,W}.

(b)

Покажите, что для любой матрицы XX размера n×pn \times p и любой матрицы UU размера n×qn \times q, такой что C(U)=C(X)C(U) = C(X),

WPU,W=WPX,W. WP_{U,W} = WP_{X,W} .
Задача 14.48

Пусть U\mathcal{U} — подпространство Rn×1\mathbb {R}^{n \times 1}, пусть AA -- матрица порядка n×nn \times n, а WW -- симметрическая неотрицательно определённая матрица порядка n×nn \times n, и пусть XX -- произвольная матрица размера n×pn \times p, столбцы которой порождают U\mathcal{U}.

?
(a)

Покажите, что AA является матрицей проектирования для U\mathcal{U} относительно WW тогда и только тогда, когда

A=PX,W+(I−PX,W)XK A = P_{X,W} + (I-P_{X,W})XK

для некоторой матрицы KK размера p×np \times n.

(b)

Покажите, что если AA является матрицей проектирования для U\mathcal{U} относительно WW, то WA=WPX,WWA = WP_{X,W}.

Задача 14.49

Пусть AA — матрица порядка n×nn \times n, а WW -- симметрическая неотрицательно определённая матрица порядка n×nn \times n.

?
(a)

Покажите (например, используя результаты упражнения 14.48), что если A′WA=WAA'WA = WA [или, что эквивалентно, если (I−A)′WA=0(I-A)'WA = 0], то AA является матрицей проектирования относительно WW, и в частности AA является матрицей проектирования для C(A)C(A) относительно WW, и, обратно, покажите, что если AA является матрицей проектирования относительно WW, то A′WA=WAA'WA = WA.

(b)

Покажите, что если AA является матрицей проектирования относительно WW, то, в частности, AA является матрицей проектирования для C(A)C(A) относительно WW.

(c)

Покажите, что AA является матрицей проектирования относительно WW тогда и только тогда, когда WAWA симметрична и WA2=WAWA^2 = WA.

Задача 14.50

Пусть U\mathcal{U} — подпространство Rn×1\mathbb {R}^{n \times 1}, пусть XX -- матрица размера n×pn \times p, столбцы которой порождают U\mathcal{U}, и пусть WW и VV -- симметрические неотрицательно определённые матрицы порядка n×nn \times n. Покажите (например, используя результат упражнения 14.46), что каждое из следующих двух условий необходимо и достаточно для того, чтобы всякая проекция yy на U\mathcal{U} относительно WW была проекцией (для всякого yy из Rn\mathbb {R}^{n}) yy на U\mathcal{U} относительно VV:

?
(a)

X′VPX,W=X′VX'VP_{X,W} = X'V, или, что эквивалентно, X′V(I−PX,W)=0X'V(I-P_{X,W}) = 0;

(b)

существует матрица QQ порядка p×pp \times p, такая что VX=WXQVX = WXQ, или, что эквивалентно, C(VX)⊂C(WX)C(VX) \subset C(WX).

Задача 14.51

Пусть XX — матрица размера n×pn \times p, а WW -- симметрическая неотрицательно определённая матрица порядка n×nn \times n. Как и в частном случае, когда WW положительно определена, положим

CW⊥(X)={y∈Rn×1:y⊥WC(X)}. C_W^{\perp }(X) = \{ y \in \mathbb {R}^{n \times 1} : y \perp _W C(X)\} .
?
(a)

Используя, например, результат упражнения 14.46, покажите, что

CW⊥(X)=N(X′W)=C(I−PX,W). C_W^{\perp }(X) = N(X'W) = C(I-P_{X,W}) .
(b)

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

dim⁡[CW⊥(X)]=n−rk⁡(WX)≥n−rk⁡(X)=n−dim⁡[C(X)]. \dim [C_W^{\perp }(X)] = n - \operatorname {rk}\left(WX\right) \geq n - \operatorname {rk}\left(X\right) = n - \dim [C(X)] .
(c)

Используя, например, результат упражнения 14.46, покажите, что для любого решения b∗b^{*} линейной системы X′WXb=X′WyX'WXb = X'Wy (относительно bb) вектор y−Xb∗y - Xb^{*} является проекцией yy на CW⊥(X)C_W^{\perp }(X) относительно WW.