5.6

Унитарные и ортогональные матрицы

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

Определите, какие из следующих матриц являются изометриями.

?
(a)

(1/2−1/201/61/6−2/61/31/31/3)\begin{pmatrix} 1/\sqrt{2} & -1/\sqrt{2} & 0 \\ 1/\sqrt{6} & 1/\sqrt{6} & -2/\sqrt{6} \\ 1/\sqrt{3} & 1/\sqrt{3} & 1/\sqrt{3} \end{pmatrix}.

(b)

(10110−1010)\begin{pmatrix} 1 & 0 & 1 \\ 1 & 0 & -1 \\ 0 & 1 & 0 \end{pmatrix}.

(c)

(0010100000010100)\begin{pmatrix} 0 & 0 & 1 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 1 & 0 & 0 \end{pmatrix}.

(d)

(eiθ10⋯00eiθ2⋯0⋮⋮⋱⋮00⋯eiθn)\begin{pmatrix} e^{i\theta_{1}} & 0 & \cdots & 0 \\ 0 & e^{i\theta_{2}} & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & e^{i\theta_{n}} \end{pmatrix}.

Задача 5.6.2

Является ли (1+i31+i6i3−2i6)\begin{pmatrix} \dfrac {1+i}{\sqrt{3}} & \dfrac {1+i}{\sqrt{6}} \\ \dfrac {i}{\sqrt{3}} & \dfrac {-2i}{\sqrt{6}} \end{pmatrix} унитарной матрицей?

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

Сколько существует матриц 3×33 \times 3, которые одновременно диагональны и ортогональны?

(b)

Сколько существует матриц n×nn \times n, которые одновременно диагональны и ортогональны?

(c)

Сколько существует матриц n×nn \times n, которые одновременно диагональны и унитарны?

Задача 5.6.4
?
(a)

При каких условиях на вещественные числа α\alpha и β\beta матрица P=(α+ββ−αα−ββ+α)P = \begin{pmatrix} \alpha +\beta & \beta -\alpha \\ \alpha -\beta & \beta +\alpha \end{pmatrix} будет ортогональной?

(b)

При каких условиях на вещественные числа α\alpha и β\beta матрица U=(0α0iβα0iβ00iβ0αiβ0α0)U = \begin{pmatrix} 0 & \alpha & 0 & i\beta \\ \alpha & 0 & i\beta & 0 \\ 0 & i\beta & 0 & \alpha \\ i\beta & 0 & \alpha & 0 \end{pmatrix} будет унитарной?

Задача 5.6.5

Пусть UU и VV — две унитарные (ортогональные) матрицы n×nn \times n.

?
(a)

Объясните, почему произведение UVUV обязательно унитарно (ортогонально).

(b)

Объясните, почему сумма U+VU+V не обязана быть унитарной (ортогональной).

(c)

Объясните, почему (Un×n00Vm×m)\begin{pmatrix} U_{n \times n} & 0 \\ 0 & V_{m \times m} \end{pmatrix} обязательно унитарна (ортогональна).

Задача 5.6.6

Докажите, как это сделал Кэли в 1846 году, что если AA косоэрмитова (или вещественная кососимметричная), то U=(I−A)(I+A)−1=(I+A)−1(I−A)U = (I-A)(I+A)^{-1} = (I+A)^{-1}(I-A) унитарна (ортогональна); для этого сначала покажите, что (I+A)−1(I+A)^{-1} существует для косоэрмитовых матриц, а также что (I−A)(I+A)−1=(I+A)−1(I−A)(I-A)(I+A)^{-1} = (I+A)^{-1}(I-A).

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

Упражнение 3.7.6 устанавливает, что A(I+A)−1=(I+A)−1AA(I+A)^{-1} = (I+A)^{-1}A всякий раз, когда I+AI+A невырождена. Существует более прямой подход к этому упражнению, но он требует теоремы о диагонализации нормальных матриц (Упражнение 7.5.5).

Задача 5.6.7

Предположим, что RR и SS — элементарные отражения.

?
(a)

Является ли (I00R)\begin{pmatrix} I & 0 \\ 0 & R \end{pmatrix} элементарным отражением?

(b)

Является ли (R00S)\begin{pmatrix} R & 0 \\ 0 & S \end{pmatrix} элементарным отражением?

Задача 5.6.8
?
(a)

Объясните, почему стандартное скалярное произведение инвариантно относительно унитарного преобразования. То есть, если UU — произвольная унитарная матрица и если u=Uxu = Ux и v=Uyv = Uy, то u∗v=x∗yu^{*}v = x^{*}y.

(b)

Для произвольных двух векторов x,y∈Rnx, y \in \mathbb {R}^{n} объясните, почему угол между ними инвариантен относительно ортогонального преобразования. То есть, если u=Pxu = Px и v=Pyv = Py, где PP — ортогональная матрица, то cos⁡θu,v=cos⁡θx,y\cos \theta_{u,v} = \cos \theta_{x,y}.

Задача 5.6.9

Пусть Um×rU_{m \times r} — матрица с ортонормированными столбцами, а Vk×nV_{k \times n} — матрица с ортонормированными строками. Для произвольной A∈Cr×kA \in \mathbb {C}^{r \times k} решите следующие задачи, используя матричную 2-норму и норму Фробениуса.

?
(a)

Определите значения ∥U∥2\left\| U\right\|_{2}, ∥V∥2\left\| V\right\|_{2}, ∥U∥F\left\| U\right\|_{F} и ∥V∥F\left\| V\right\|_{F}.

(b)

Покажите, что ∥UAV∥2=∥A∥2\left\| UAV\right\|_{2} = \left\| A\right\|_{2}.

(c)

Покажите, что ∥UAV∥F=∥A∥F\left\| UAV\right\|_{F} = \left\| A\right\|_{F}.

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

В частности, эти свойства верны, когда UU и VV — унитарные матрицы. Из-за пунктов (б) и (в) 2-норму и норму Фробениуса называют унитарно инвариантными нормами.

Задача 5.6.10

Пусть u=(−213−1)u = \begin{pmatrix} -2 \\ 1 \\ 3 \\ -1 \end{pmatrix} и v=(140−1)v = \begin{pmatrix} 1 \\ 4 \\ 0 \\ -1 \end{pmatrix}.

?
(a)

Определите ортогональную проекцию uu на ⟨v⟩\left\langle v\right\rangle.

(b)

Определите ортогональную проекцию vv на ⟨u⟩\left\langle u\right\rangle.

(c)

Определите ортогональную проекцию uu на v⊥v^{\perp }.

(d)

Определите ортогональную проекцию vv на u⊥u^{\perp }.

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

Формула (5.6.6): для ненулевого вектора u∈Cnu \in \mathbb {C}^{n} ортогональный проектор на ⟨u⟩\left\langle u\right\rangle равен uu∗u∗u\dfrac {uu^{*}}{u^{*}u}, а ортогональный проектор на u⊥u^{\perp } равен I−uu∗u∗uI - \dfrac {uu^{*}}{u^{*}u}.

Задача 5.6.11

Рассмотрим элементарные ортогональные проекторы Q=I−uu∗Q = I - uu^{*} (где ∥u∥=1\left\| u\right\| =1).

?
(a)

Докажите, что QQ вырождена.

(b)

Теперь докажите, что если QQ имеет размер n×nn \times n, то rk⁡(Q)=n−1\operatorname {rk}\left(Q\right) = n-1.

Задача 5.6.12

Для векторов u,x∈Cnu, x \in \mathbb {C}^{n} таких, что ∥u∥=1\left\| u\right\| =1, пусть pp — ортогональная проекция xx на ⟨u⟩\left\langle u\right\rangle. Объясните, почему ∥p∥≤∥x∥\left\| p\right\| \leq \left\| x\right\|, причём равенство выполняется тогда и только тогда, когда xx является скалярным кратным uu.

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

Формула (5.6.5) утверждает, что ∥p∥=∣u∗x∣\left\| p\right\| = \left|u^{*}x\right|. Неравенство Коши—Буняковского—Шварца (5.1.3) утверждает, что ∣u∗x∣≤∥u∥∥x∥\left|u^{*}x\right| \leq \left\| u\right\| \left\| x\right\|.

Задача 5.6.13

Пусть x=13(1−2−2)x = \dfrac {1}{3}\begin{pmatrix} 1 \\ -2 \\ -2 \end{pmatrix}.

?
(a)

Определите элементарное отражение RR такое, что RxRx лежит на оси xx.

(b)

Проверьте прямым вычислением, что ваше отражение RR симметрично, ортогонально и инволютивно.

(c)

Дополните xx до ортонормированного базиса пространства R3\mathbb {R}^{3}, используя элементарное отражение.

Задача 5.6.14

Пусть R=I−2uu∗R = I - 2uu^{*}, где ∥un×1∥=1\left\| u_{n \times 1}\right\| = 1. Если xx — неподвижная точка RR в том смысле, что Rx=xRx = x, и если n>1n>1, докажите, что xx должен быть ортогонален uu, а затем изобразите эту ситуацию в R3\mathbb {R}^{3}.

Плоскость неподвижных точек u^{\perp } отражения R = I - 2uu^{*}, с изображённым u и неподвижной точкой x.Плоскость неподвижных точек u^{\perp } отражения R = I - 2uu^{*}, с изображённым u и неподвижной точкой x.

?
Задача 5.6.15

Пусть x,y∈Rnx, y \in \mathbb {R}^{n} — векторы такие, что ∥x∥=∥y∥\left\| x\right\| = \left\| y\right\|, но x≠yx \neq y. Объясните, как построить элементарное отражение RR такое, что Rx=yRx = y.

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

Вектор uu, задающий RR, можно определить визуально в R3\mathbb {R}^{3}, рассматривая картину отражения относительно плоскости из построения элементарного отражения (рисунок 5.6.2 источника).

Задача 5.6.16

Пусть xn×1x_{n \times 1} — вектор такой, что ∥x∥=1\left\| x\right\| = 1, и разобьём xx на блоки x=(x1x~)x = \begin{pmatrix} x_{1} \\ \tilde{x} \end{pmatrix}, где x~\tilde{x} имеет размер (n−1)×1(n-1) \times 1.

?
(a)

Если компоненты xx вещественны и x1≠1x_{1} \neq 1, покажите, что P=(x1x~Tx~I−αx~x~T)P = \begin{pmatrix} x_{1} & \tilde{x}^{T} \\ \tilde{x} & I - \alpha \tilde{x}\tilde{x}^{T} \end{pmatrix}, где α=11−x1\alpha = \dfrac {1}{1-x_{1}}, является ортогональной матрицей.

(b)

Предположим, что компоненты xx комплексны. Если ∣x1∣≠1\left|x_{1}\right| \neq 1 и μ\mu — число, определённое как μ=1\mu = 1, если x1x_{1} вещественно, и μ=x1/∣x1∣\mu = x_{1}/\left|x_{1}\right|, если x1x_{1} не вещественно, покажите, что матрица U=(x1μ2x~∗x~μ(I−αx~x~∗))U = \begin{pmatrix} x_{1} & \mu^{2}\tilde{x}^{*} \\ \tilde{x} & \mu (I-\alpha \tilde{x}\tilde{x}^{*}) \end{pmatrix}, где α=11−∣x1∣\alpha = \dfrac {1}{1-\left|x_{1}\right|}, унитарна.

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

Эти результаты дают простой способ дополнить заданный вектор до ортонормированного базиса всего пространства Rn\mathbb {R}^{n} или Cn\mathbb {C}^{n}.

Задача 5.6.17

Выполните следующую последовательность поворотов в R3\mathbb {R}^{3}, начиная с v0=(11−1)v_{0} = \begin{pmatrix} 1 \\ 1 \\ -1 \end{pmatrix}.

  1. Поверните v0v_{0} против часовой стрелки на 45°45° вокруг оси xx, получив v1v_{1}.

  2. Поверните v1v_{1} по часовой стрелке на 90°90° вокруг оси yy, получив v2v_{2}.

  3. Поверните v2v_{2} против часовой стрелки на 30°30° вокруг оси zz, получив v3v_{3}. Найдите координаты v3v_{3}, а также ортогональную матрицу QQ такую, что Qv0=v3Qv_{0} = v_{3}.

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

Матрицы поворота против часовой стрелки вокруг каждой координатной оси на угол θ\theta имеют вид Px(θ)=(1000cos⁡θ−sin⁡θ0sin⁡θcos⁡θ)P_{x}(\theta ) = \begin{pmatrix} 1 & 0 & 0 \\ 0 & \cos \theta & -\sin \theta \\ 0 & \sin \theta & \cos \theta \end{pmatrix}, Py(θ)=(cos⁡θ0sin⁡θ010−sin⁡θ0cos⁡θ)P_{y}(\theta ) = \begin{pmatrix} \cos \theta & 0 & \sin \theta \\ 0 & 1 & 0 \\ -\sin \theta & 0 & \cos \theta \end{pmatrix}, Pz(θ)=(cos⁡θ−sin⁡θ0sin⁡θcos⁡θ0001)P_{z}(\theta ) = \begin{pmatrix} \cos \theta & -\sin \theta & 0 \\ \sin \theta & \cos \theta & 0 \\ 0 & 0 & 1 \end{pmatrix}; поворот по часовой стрелке на угол θ\theta — это P⋆(−θ)P_{\star }(-\theta ).

Задача 5.6.18

Имеет ли значение порядок, в котором выполняются повороты в R3\mathbb {R}^{3}? Например, предположим, что вектор v∈R3v \in \mathbb {R}^{3} сначала поворачивается против часовой стрелки вокруг оси xx на угол θ\theta, а затем этот вектор поворачивается против часовой стрелки вокруг оси yy на угол φ\varphi. Будет ли результат таким же, как если сначала повернуть vv против часовой стрелки вокруг оси yy на угол φ\varphi, а затем выполнить поворот против часовой стрелки вокруг оси xx на угол θ\theta?

?
Задача 5.6.19

Для каждого ненулевого вектора u∈Cnu \in \mathbb {C}^{n} докажите, что dim⁡u⊥=n−1\dim u^{\perp } = n-1.

?
Задача 5.6.20

Матрица, удовлетворяющая условию A2=IA^{2} = I, называется инволюцией или инволютивной матрицей, а матрица PP, удовлетворяющая условию P2=PP^{2} = P, называется проектором или идемпотентной матрицей. Покажите, что между множеством инволюций и множеством проекторов в Cn×n\mathbb {C}^{n \times n} существует взаимно однозначное соответствие.

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

Формулы (5.6.6)-(5.6.7) показывают, что элементарный отражатель R=I−2uu∗/u∗uR = I - 2uu^{*}/u^{*}u (являющийся инволюцией) строится из проектора P=uu∗/u∗uP = uu^{*}/u^{*}u по формуле R=I−2PR = I - 2P.

Задача 5.6.21

При использовании компьютера для построения и отображения трёхмерного выпуклого многогранника желательно не отрисовывать те грани, которые должны быть скрыты от наблюдателя. Векторное произведение в R3\mathbb {R}^{3} можно использовать для того, чтобы определить, какие грани видимы, а какие нет: если u=(u1u2u3)u = \begin{pmatrix} u_{1} \\ u_{2} \\ u_{3} \end{pmatrix} и v=(v1v2v3)v = \begin{pmatrix} v_{1} \\ v_{2} \\ v_{3} \end{pmatrix}, то u×v=(u2v3−u3v2u3v1−u1v3u1v2−u2v1)u \times v = \begin{pmatrix} u_{2}v_{3}-u_{3}v_{2} \\ u_{3}v_{1}-u_{1}v_{3} \\ u_{1}v_{2}-u_{2}v_{1} \end{pmatrix} — это вектор, ортогональный как uu, так и vv, направление которого задаётся правилом правой руки.

Слева: наблюдатель, расположенный на положительной полуоси x и смотрящий в сторону начала координат (источник — рисунок 5.6.6). В центре: правило правой руки для u \times v (источник — рисунок 5.6.8). Справа: грань с вершинами p_0,p_1,p_2 и внешней нормалью n (источник — рисунок 5.6.9).Слева: наблюдатель, расположенный на положительной полуоси x и смотрящий в сторону начала координат (источник — рисунок 5.6.6). В центре: правило правой руки для u \times v (источник — рисунок 5.6.8). Справа: грань с вершинами p_0,p_1,p_2 и внешней нормалью n (источник — рисунок 5.6.9).

Предположим, что начало координат находится внутри многогранника, и рассмотрим отдельную грань и три вершины p0,p1,p2p_{0}, p_{1}, p_{2} на этой грани, перечисленные против часовой стрелки, если смотреть снаружи грани. Вектор n=(p1−p0)×(p2−p1)n = (p_{1}-p_{0}) \times (p_{2}-p_{1}) ортогонален грани и направлен наружу. Объясните, почему внешняя сторона грани видна наблюдателю, расположенному на положительной полуоси xx и смотрящему в сторону начала координат, тогда и только тогда, когда первая компонента внешней нормали nn положительна. Другими словами, грань отрисовывается тогда и только тогда, когда n1>0n_{1} > 0.

?