12

Проекции и матрицы проектирования

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

Пусть Y\mathbf{Y} — матрица в линейном пространстве V\mathcal{V}, пусть U\mathcal{U} и W\mathcal{W} — подпространства V\mathcal{V}, и пусть {X1,…,Xs}\{ \mathbf{X}_{1}, \ldots , \mathbf{X}_{s}\} — множество матриц, порождающее U\mathcal{U}, а {Z1,…,Zt}\{ \mathbf{Z}_{1}, \ldots , \mathbf{Z}_{t}\} — множество, порождающее W\mathcal{W}. Проверьте, что Y⊥U\mathbf{Y} \perp \mathcal{U} тогда и только тогда, когда Y⋅Xi=0\mathbf{Y} \cdot \mathbf{X}_{i} = 0 для i=1,…,si = 1, \ldots , s (то есть что Y\mathbf{Y} ортогональна U\mathcal{U} тогда и только тогда, когда Y\mathbf{Y} ортогональна каждой из матриц X1,…,Xs\mathbf{X}_{1}, \ldots , \mathbf{X}_{s}); и, аналогично, что U⊥W\mathcal{U} \perp \mathcal{W} тогда и только тогда, когда Xi⋅Zj=0\mathbf{X}_{i} \cdot \mathbf{Z}_{j} = 0 для i=1,…,si = 1, \ldots , s и j=1,…,tj = 1, \ldots , t (то есть что U\mathcal{U} ортогонально W\mathcal{W} тогда и только тогда, когда каждая из матриц X1,…,Xs\mathbf{X}_{1}, \ldots , \mathbf{X}_{s} ортогональна каждой из матриц Z1,…,Zt\mathbf{Z}_{1}, \ldots , \mathbf{Z}_{t}).

?
Задача 12.2

Пусть U\mathcal{U} и V\mathcal{V} — подпространства Rm×n\mathbb {R}^{m \times n}. Покажите, что если dim⁡(V)>dim⁡(U)\dim (\mathcal{V}) > \dim (\mathcal{U}), то в V\mathcal{V} найдётся ненулевая матрица, ортогональная U\mathcal{U}.

?
Задача 12.3

Пусть U\mathcal{U} — подпространство линейного пространства Rm\mathbb {R}^{m} всех mm-мерных вектор-столбцов. Пусть M\mathcal{M} — подпространство Rm×n\mathbb {R}^{m \times n}, определяемое условием: W∈M\mathbf{W} \in \mathcal{M} тогда и только тогда, когда W=(w1,…,wn)\mathbf{W} = (\mathbf{w}_{1}, \ldots , \mathbf{w}_{n}) для некоторых векторов w1,…,wn\mathbf{w}_{1}, \ldots , \mathbf{w}_{n} из U\mathcal{U}. Пусть Z\mathbf{Z} — проекция (относительно обычного скалярного произведения) m×nm \times n матрицы Y\mathbf{Y} на M\mathcal{M}, а X\mathbf{X} — произвольная m×pm \times p матрица, столбцы которой порождают U\mathcal{U}. Покажите, что Z=XB∗\mathbf{Z} = \mathbf{X}\mathbf{B}^{*} для любого решения B∗\mathbf{B}^{*} линейной системы

X′XB=X′Y(относительно B). \mathbf{X}'\mathbf{X}\mathbf{B} = \mathbf{X}'\mathbf{Y} \quad \text{(относительно }\mathbf{B}\text{)}.
?
Задача 12.4

Проекция (относительно обычного скалярного произведения) nn-мерного вектор-столбца y\mathbf{y} на подпространство U\mathcal{U} пространства Rn\mathbb {R}^{n} в частном случае n=3n = 3, y=(3,−38/5,74/5)′\mathbf{y} = (3, -38/5, 74/5)' и U=sp⁡{x1,x2,x3}\mathcal{U} = \operatorname {sp}\{ \mathbf{x}_{1}, \mathbf{x}_{2}, \mathbf{x}_{3}\}, где

x1=(036),x2=(−224),x3=(−212), \mathbf{x}_{1} = \left(\begin{smallmatrix} 0 \\ 3 \\ 6 \end{smallmatrix}\right), \quad \mathbf{x}_{2} = \left(\begin{smallmatrix} -2 \\ 2 \\ 4 \end{smallmatrix}\right), \quad \mathbf{x}_{3} = \left(\begin{smallmatrix} -2 \\ 1 \\ 2 \end{smallmatrix}\right),

была определена как вектор (3,22/5,44/5)′(3, 22/5, 44/5)' --- при этом было отмечено, что x1\mathbf{x}_{1} и x2\mathbf{x}_{2} линейно независимы и что x3=x2−(1/3)x1\mathbf{x}_{3} = \mathbf{x}_{2} - (1/3)\mathbf{x}_{1}, откуда dim⁡(U)=2\dim (\mathcal{U}) = 2. Заново вычислите проекцию y\mathbf{y} на U\mathcal{U} (в этом частном случае), взяв в качестве X\mathbf{X} матрицу размера 3×23 \times 2

(0−23264) \left(\begin{smallmatrix} 0 & -2 \\ 3 & 2 \\ 6 & 4 \end{smallmatrix}\right)

и выполнив следующие два шага: (1) вычислите решение нормальных уравнений X′Xb=X′y\mathbf{X}'\mathbf{X}\mathbf{b} = \mathbf{X}'\mathbf{y}; и (2) умножьте X\mathbf{X} справа на решение, вычисленное на шаге (1).

?
Задача 12.5

Пусть X\mathbf{X} — произвольная матрица размера n×pn \times p. Если p×np \times n матрица B∗\mathbf{B}^{*} является решением линейной системы X′XB=X′\mathbf{X}'\mathbf{X}\mathbf{B} = \mathbf{X}' (относительно B\mathbf{B}), то B∗\mathbf{B}^{*} является обобщённой обратной матрицей для X\mathbf{X}, а XB∗\mathbf{X}\mathbf{B}^{*} симметрична. Покажите, что, обратно, если p×np \times n матрица G\mathbf{G} является обобщённой обратной матрицей для X\mathbf{X} и если XG\mathbf{X}\mathbf{G} симметрична, то X′XG=X′\mathbf{X}'\mathbf{X}\mathbf{G} = \mathbf{X}' (то есть G\mathbf{G} — решение X′XB=X′\mathbf{X}'\mathbf{X}\mathbf{B} = \mathbf{X}').

?
Задача 12.6

Используя результат части (b) упражнения 9.3 (или иным способом), покажите, что для любой ненулевой симметрической матрицы A\mathbf{A}

PA=B(TB)−1T, \mathbf{P}_{\mathbf{A}} = \mathbf{B}(\mathbf{T}\mathbf{B})^{-1}\mathbf{T} ,

где B\mathbf{B} — произвольная матрица полного ранга по столбцам, а T\mathbf{T} — произвольная матрица полного ранга по строкам, такие что A=BT\mathbf{A} = \mathbf{B}\mathbf{T}. (Невырожденность TB\mathbf{T}\mathbf{B} следует из результата упражнения 8.3.)

?
Задача 12.7

Пусть V\mathcal{V} — kk-мерное подпространство линейного пространства Rn\mathbb {R}^{n} всех nn-мерных вектор-столбцов. Пусть X\mathbf{X} — произвольная n×pn \times p матрица, столбцы которой порождают V\mathcal{V}, пусть U\mathcal{U} — подпространство V\mathcal{V}, и пусть A\mathbf{A} — матрица проектирования для U\mathcal{U}. Покажите (1), что матрица B\mathbf{B} (размера n×nn \times n) такова, что By\mathbf{B}\mathbf{y} является проекцией y\mathbf{y} на U\mathcal{U} для каждого y∈V\mathbf{y} \in \mathcal{V}, тогда и только тогда, когда B=A+Z∗\mathbf{B} = \mathbf{A} + \mathbf{Z}_{*} для некоторого решения Z∗\mathbf{Z}_{*} однородной линейной системы X′Z=0\mathbf{X}'\mathbf{Z} = \mathbf{0} (относительно n×nn \times n матрицы Z\mathbf{Z}), и (2), что, если только не k=nk = n, существует более одной матрицы B\mathbf{B}, такой что By\mathbf{B}\mathbf{y} является проекцией y\mathbf{y} на U\mathcal{U} для каждого y∈V\mathbf{y} \in \mathcal{V}.

?
Задача 12.8

Пусть {A1,…,Ak}\{ \mathbf{A}_{1}, \ldots , \mathbf{A}_{k}\} — непустое линейно независимое множество матриц в линейном пространстве V\mathcal{V}. И определим (как в процессе ортогонализации Грама--Шмидта) kk ненулевых ортогональных линейных комбинаций, скажем, B1,…,Bk\mathbf{B}_{1}, \ldots , \mathbf{B}_{k}, матриц A1,…,Ak\mathbf{A}_{1}, \ldots , \mathbf{A}_{k} следующим образом:

B1=A1,B2=A2−x12B1,⋮Bj=Aj−xj−1,jBj−1−⋯−x1jB1,⋮Bk=Ak−xk−1,kBk−1−⋯−x1kB1, \begin{align} \mathbf{B}_{1} & = \mathbf{A}_{1}, \\ \mathbf{B}_{2} & = \mathbf{A}_{2} - x_{12}\mathbf{B}_{1}, \\ & \vdots \\ \mathbf{B}_{j} & = \mathbf{A}_{j} - x_{j-1,j}\mathbf{B}_{j-1} - \cdots - x_{1j}\mathbf{B}_{1}, \\ & \vdots \\ \mathbf{B}_{k} & = \mathbf{A}_{k} - x_{k-1,k}\mathbf{B}_{k-1} - \cdots - x_{1k}\mathbf{B}_{1}, \end{align}

где (для i<j=1,…,ki < j = 1, \ldots , k)

xij=Aj⋅BiBi⋅Bi. x_{ij} = \frac{\mathbf{A}_{j} \cdot \mathbf{B}_{i}}{\mathbf{B}_{i} \cdot \mathbf{B}_{i}} .

Покажите, что Bj\mathbf{B}_{j} является (ортогональной) проекцией Aj\mathbf{A}_{j} на некоторое подпространство Uj\mathcal{U}_{j} (пространства V\mathcal{V}), и опишите Uj\mathcal{U}_{j} (j=2,…,kj = 2, \ldots , k).

?