19

Минимизация полинома второй степени (от n переменных) при линейных ограничениях

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

Пусть a\mathbf{a} — n×1n \times 1-вектор (не ограниченных) переменных, и пусть f(a)=a′Va−2b′af(\mathbf{a}) = \mathbf{a}'\mathbf{V}\mathbf{a} - 2\mathbf{b}'\mathbf{a}, где V\mathbf{V} — n×nn \times n-матрица, а b\mathbf{b} — n×1n \times 1-вектор. Покажите, что если V\mathbf{V} не является неотрицательно определённой или b∉C(V)\mathbf{b} \notin C(\mathbf{V}), то f(a)f(\mathbf{a}) не ограничен снизу, то есть для любого скаляра cc найдётся вектор a∗\mathbf{a}_{*} такой, что f(a∗)<cf(\mathbf{a}_{*}) < c.

?
Задача 19.2

Пусть V\mathbf{V} — симметрическая n×nn \times n-матрица, а X\mathbf{X} — n×pn \times p-матрица. Покажите, что для любой p×pp \times p-матрицы U\mathbf{U} такой, что C(X)⊂C(V+XUX′)C(\mathbf{X}) \subset C(\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}'),

  1. (V+XUX′)(V+XUX′)−V=V(\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}')(\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}')^{-}\mathbf{V} = \mathbf{V};

  2. V(V+XUX′)−(V+XUX′)=V\mathbf{V}(\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}')^{-}(\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}') = \mathbf{V}.

?
Задача 19.3

Пусть V\mathbf{V} — симметрическая неотрицательно определённая n×nn \times n-матрица, X\mathbf{X} — n×pn \times p-матрица, B\mathbf{B} — n×sn \times s-матрица такая, что C(B)⊂C(V,X)C(\mathbf{B}) \subset C(\mathbf{V}, \mathbf{X}), а D\mathbf{D} — p×sp \times s-матрица такая, что C(D)⊂C(X′)C(\mathbf{D}) \subset C(\mathbf{X}'). Далее, пусть U\mathbf{U} — произвольная p×pp \times p-матрица такая, что C(X)⊂C(V+XUX′)C(\mathbf{X}) \subset C(\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}'), а W\mathbf{W} — произвольная обобщённая обратная матрица к V+XUX′\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}'. Дайте короткое доказательство того, что A∗\mathbf{A}_{*} и R∗\mathbf{R}_{*} являются соответственно первой (n×sn \times s) и второй (p×sp \times s) частями решения (совместной) линейной системы

(VXX′0)(AR)=(BD)(∗) \begin{pmatrix} \mathbf{V} & \mathbf{X} \\ \mathbf{X}' & \mathbf{0} \end{pmatrix} \begin{pmatrix} \mathbf{A} \\ \mathbf{R} \end{pmatrix} = \begin{pmatrix} \mathbf{B} \\ \mathbf{D} \end{pmatrix} \tag {$*$}

(относительно n×sn \times s-матрицы A\mathbf{A} и p×sp \times s-матрицы R\mathbf{R}) тогда и только тогда, когда

R∗=T∗+UD \mathbf{R}_{*} = \mathbf{T}_{*} + \mathbf{U}\mathbf{D}

и

A∗=WB−WXT∗+[I−W(V+XUX′)]L \mathbf{A}_{*} = \mathbf{W}\mathbf{B} - \mathbf{W}\mathbf{X}\mathbf{T}_{*} + [\mathbf{I} - \mathbf{W}(\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}')]\mathbf{L}

для некоторого решения T∗\mathbf{T}_{*} (совместной) линейной системы

X′WXT=X′WB−D(∗∗) \mathbf{X}'\mathbf{W}\mathbf{X}\mathbf{T} = \mathbf{X}'\mathbf{W}\mathbf{B} - \mathbf{D} \tag {$**$}

(относительно p×sp \times s-матрицы T\mathbf{T}) и некоторой n×sn \times s-матрицы L\mathbf{L}. Сделайте это, воспользовавшись результатом о том, что если коэффициентная матрица, правая часть и матрица неизвестных линейной системы HY=S\mathbf{H}\mathbf{Y} = \mathbf{S} (относительно Y\mathbf{Y}) разбиты (согласованно) на блоки

H=(H11H12H21H22),S=(S1S2),иY=(Y1Y2), \mathbf{H} = \begin{pmatrix} \mathbf{H}_{11} & \mathbf{H}_{12} \\ \mathbf{H}_{21} & \mathbf{H}_{22} \end{pmatrix}, \quad \mathbf{S} = \begin{pmatrix} \mathbf{S}_{1} \\ \mathbf{S}_{2} \end{pmatrix}, \quad \text{и} \quad \mathbf{Y} = \begin{pmatrix} \mathbf{Y}_{1} \\ \mathbf{Y}_{2} \end{pmatrix},

и если

C(H12)⊂C(H11),C(S1)⊂C(H11),иR(H21)⊂R(H11), C(\mathbf{H}_{12}) \subset C(\mathbf{H}_{11}), \quad C(\mathbf{S}_{1}) \subset C(\mathbf{H}_{11}), \quad \text{и} \quad R(\mathbf{H}_{21}) \subset R(\mathbf{H}_{11}),

то матрица Y∗=(Y1∗Y2∗)\mathbf{Y}^{*} = \begin{pmatrix} \mathbf{Y}_{1}^{*} \\ \mathbf{Y}_{2}^{*} \end{pmatrix} является решением линейной системы HY=S\mathbf{H}\mathbf{Y} = \mathbf{S} тогда и только тогда, когда Y2∗\mathbf{Y}_{2}^{*} является решением линейной системы

(H22−H21H11−H12)Y2=S2−H21H11−S1(относительно Y2) (\mathbf{H}_{22} - \mathbf{H}_{21}\mathbf{H}_{11}^{-}\mathbf{H}_{12})\mathbf{Y}_{2} = \mathbf{S}_{2} - \mathbf{H}_{21}\mathbf{H}_{11}^{-}\mathbf{S}_{1} \quad \text{(относительно }\mathbf{Y}_{2}\text{)}

и Y1∗\mathbf{Y}_{1}^{*} и Y2∗\mathbf{Y}_{2}^{*} являются решением линейной системы

H11Y1+H12Y2=S1(относительно Y1 и Y2). \mathbf{H}_{11}\mathbf{Y}_{1} + \mathbf{H}_{12}\mathbf{Y}_{2} = \mathbf{S}_{1} \quad \text{(относительно }\mathbf{Y}_{1}\text{ и }\mathbf{Y}_{2}\text{)}.
?
Задача 19.4

Пусть V\mathbf{V} — симметрическая n×nn \times n-матрица, а X\mathbf{X} — n×pn \times p-матрица. Далее, пусть U=X′TT′X\mathbf{U} = \mathbf{X}'\mathbf{T}\mathbf{T}'\mathbf{X}, где T\mathbf{T} — произвольная матрица, столбцы которой порождают нулевое пространство матрицы V\mathbf{V}. Покажите, что C(X)⊂C(V+XUX′)C(\mathbf{X}) \subset C(\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}') и что C(V)C(\mathbf{V}) и C(XUX′)C(\mathbf{X}\mathbf{U}\mathbf{X}') существенно не пересекаются, а R(V)R(\mathbf{V}) и R(XUX′)R(\mathbf{X}\mathbf{U}\mathbf{X}') существенно не пересекаются (даже без предположения о неотрицательной определённости V\mathbf{V}).

?
Задача 19.5

Пусть V\mathbf{V} — симметрическая неотрицательно определённая n×nn \times n-матрица, а X\mathbf{X} — n×pn \times p-матрица. Далее, пусть Z\mathbf{Z} — произвольная матрица, столбцы которой порождают N(X′)N(\mathbf{X}') или, что то же самое, C⊥(X)C^{\perp }(\mathbf{X}). И примем ту же терминологию, что и в задаче 17.20.

?
(a)

Используя результат части (a) задачи 17.20 (или иным способом), покажите, что n×nn \times n-матрица H\mathbf{H} является матрицей проектирования на C(X)C(\mathbf{X}) вдоль C(VZ)C(\mathbf{V}\mathbf{Z}) тогда и только тогда, когда H′\mathbf{H}' является первой (n×nn \times n) частью решения совместной линейной системы

(VXX′0)(AR)=(0X′)(E.1) \begin{pmatrix} \mathbf{V} & \mathbf{X} \\ \mathbf{X}' & \mathbf{0} \end{pmatrix} \begin{pmatrix} \mathbf{A} \\ \mathbf{R} \end{pmatrix} = \begin{pmatrix} \mathbf{0} \\ \mathbf{X}' \end{pmatrix} \tag {E.1}

(относительно n×nn \times n-матрицы A\mathbf{A} и p×np \times n-матрицы R\mathbf{R}).

(b)

Пусть U\mathbf{U} — произвольная p×pp \times p-матрица такая, что C(X)⊂C(V+XUX′)C(\mathbf{X}) \subset C(\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}'), а W\mathbf{W} — произвольная обобщённая обратная матрица к V+XUX′\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}'; покажите, что n×nn \times n-матрица H\mathbf{H} является матрицей проектирования на C(X)C(\mathbf{X}) вдоль C(VZ)C(\mathbf{V}\mathbf{Z}) тогда и только тогда, когда

H=PX,W+K[I−(V+XUX′)W] \mathbf{H} = \mathbf{P}_{\mathbf{X},\mathbf{W}} + \mathbf{K}[\mathbf{I} - (\mathbf{V} + \mathbf{X}\mathbf{U}\mathbf{X}')\mathbf{W}]

для некоторой n×nn \times n-матрицы K\mathbf{K}.

Задача 19.6

Пусть V\mathbf{V} — симметрическая неотрицательно определённая n×nn \times n-матрица, W\mathbf{W} — n×nn \times n-матрица, а X\mathbf{X} — n×pn \times p-матрица. Покажите, что для того, чтобы матрица WX(X′WX)−X′\mathbf{W}\mathbf{X}(\mathbf{X}'\mathbf{W}\mathbf{X})^{-}\mathbf{X}' была первой (n×nn \times n) частью некоторого решения (совместной) линейной системы

(VXX′0)(AR)=(0X′)(E.2) \begin{pmatrix} \mathbf{V} & \mathbf{X} \\ \mathbf{X}' & \mathbf{0} \end{pmatrix} \begin{pmatrix} \mathbf{A} \\ \mathbf{R} \end{pmatrix} = \begin{pmatrix} \mathbf{0} \\ \mathbf{X}' \end{pmatrix} \tag {E.2}

(относительно n×nn \times n-матрицы A\mathbf{A} и p×np \times n-матрицы R\mathbf{R}), необходимо и достаточно, чтобы C(VWX)⊂C(X)C(\mathbf{V}\mathbf{W}\mathbf{X}) \subset C(\mathbf{X}) и rk⁡(X′WX)=rk⁡(X)\operatorname {rk}\left(\mathbf{X}'\mathbf{W}\mathbf{X}\right) = \operatorname {rk}\left(\mathbf{X}\right).

?
Задача 19.7

Пусть a\mathbf{a} — n×1n \times 1-вектор переменных, на который наложено ограничение X′a=d\mathbf{X}'\mathbf{a} = \mathbf{d}, где X\mathbf{X} — n×pn \times p-матрица, а d\mathbf{d} — p×1p \times 1-вектор такой, что d∈C(X′)\mathbf{d} \in C(\mathbf{X}'). Пусть f(a)=a′Va−2b′af(\mathbf{a}) = \mathbf{a}'\mathbf{V}\mathbf{a} - 2\mathbf{b}'\mathbf{a}, где V\mathbf{V} — симметрическая неотрицательно определённая n×nn \times n-матрица, а b\mathbf{b} — n×1n \times 1-вектор такой, что b∈C(V,X)\mathbf{b} \in C(\mathbf{V}, \mathbf{X}). Далее, пусть g(a)=a′(V+W)a−2(b+c)′ag(\mathbf{a}) = \mathbf{a}'(\mathbf{V} + \mathbf{W})\mathbf{a} - 2(\mathbf{b} + \mathbf{c})'\mathbf{a}, где W\mathbf{W} — произвольная n×nn \times n-матрица такая, что C(W)⊂C(X)C(\mathbf{W}) \subset C(\mathbf{X}) и R(W)⊂R(X′)R(\mathbf{W}) \subset R(\mathbf{X}'), а c\mathbf{c} — произвольный n×1n \times 1-вектор из C(X)C(\mathbf{X}). Покажите, что минимизация g(a)g(\mathbf{a}) при ограничении X′a=d\mathbf{X}'\mathbf{a} = \mathbf{d} эквивалентна минимизации f(a)f(\mathbf{a}) при том же ограничении [в том смысле, что g(a)g(\mathbf{a}) и f(a)f(\mathbf{a}) достигают своих минимальных значений в одних и тех же точках].

?
Задача 19.8

Пусть V\mathbf{V} — симметрическая неотрицательно определённая n×nn \times n-матрица, W\mathbf{W} — n×nn \times n-матрица, X\mathbf{X} — n×pn \times p-матрица, f\mathbf{f} — n×1n \times 1-вектор, а d\mathbf{d} — p×1p \times 1-вектор. Далее, пусть b\mathbf{b} — n×1n \times 1-вектор такой, что b∈C(V,X)\mathbf{b} \in C(\mathbf{V}, \mathbf{X}). Покажите, что для того, чтобы вектор W(I−PX,W)f+WX(X′WX)−d\mathbf{W}(\mathbf{I} - \mathbf{P}_{\mathbf{X},\mathbf{W}})\mathbf{f} + \mathbf{W}\mathbf{X}(\mathbf{X}'\mathbf{W}\mathbf{X})^{-}\mathbf{d} был решением, для каждого d∈C(X′)\mathbf{d} \in C(\mathbf{X}'), задачи минимизации полинома второй степени a′Va−2b′a\mathbf{a}'\mathbf{V}\mathbf{a} - 2\mathbf{b}'\mathbf{a} (относительно a\mathbf{a}) при ограничении X′a=d\mathbf{X}'\mathbf{a} = \mathbf{d}, необходимо и достаточно, чтобы

VWf−b∈C(X),(E.3) \mathbf{V}\mathbf{W}\mathbf{f} - \mathbf{b} \in C(\mathbf{X}), \tag {E.3} C(VWX)⊂C(X),(E.4) C(\mathbf{V}\mathbf{W}\mathbf{X}) \subset C(\mathbf{X}), \tag {E.4}

и

rk⁡(X′WX)=rk⁡(X).(E.5) \operatorname {rk}\left(\mathbf{X}'\mathbf{W}\mathbf{X}\right) = \operatorname {rk}\left(\mathbf{X}\right). \tag {E.5}
?
Задача 19.9

Пусть V\mathbf{V} и W\mathbf{W} — n×nn \times n-матрицы, а X\mathbf{X} — n×pn \times p-матрица. Покажите, что если V\mathbf{V} и W\mathbf{W} невырожденны, то условие C(VWX)⊂C(X)C(\mathbf{V}\mathbf{W}\mathbf{X}) \subset C(\mathbf{X}) эквивалентно условию C(V−1X)⊂C(WX)C(\mathbf{V}^{-1}\mathbf{X}) \subset C(\mathbf{W}\mathbf{X}), а также эквивалентно условию C(W−1V−1X)⊂C(X)C(\mathbf{W}^{-1}\mathbf{V}^{-1}\mathbf{X}) \subset C(\mathbf{X}).

?
Задача 19.10

Пусть V\mathbf{V} — симметрическая положительно определённая n×nn \times n-матрица, W\mathbf{W} — n×nn \times n-матрица, X\mathbf{X} — n×pn \times p-матрица, а d\mathbf{d} — p×1p \times 1-вектор. Покажите, что для того, чтобы вектор WX(X′WX)−d\mathbf{W}\mathbf{X}(\mathbf{X}'\mathbf{W}\mathbf{X})^{-}\mathbf{d} был решением, для каждого d∈C(X′)\mathbf{d} \in C(\mathbf{X}'), задачи минимизации квадратичной формы a′Va\mathbf{a}'\mathbf{V}\mathbf{a} (относительно a\mathbf{a}) при ограничении X′a=d\mathbf{X}'\mathbf{a} = \mathbf{d}, необходимо и достаточно, чтобы V−1PX,W′\mathbf{V}^{-1}\mathbf{P}_{\mathbf{X},\mathbf{W}'} была симметрична и rk⁡(X′WX)=rk⁡(X)\operatorname {rk}\left(\mathbf{X}'\mathbf{W}\mathbf{X}\right) = \operatorname {rk}\left(\mathbf{X}\right). Покажите, что также необходимо и достаточно, чтобы (I−PX,W′′)V−1PX,W′=0(\mathbf{I} - \mathbf{P}_{\mathbf{X},\mathbf{W}'}')\mathbf{V}^{-1}\mathbf{P}_{\mathbf{X},\mathbf{W}'} = \mathbf{0} и rk⁡(X′WX)=rk⁡(X)\operatorname {rk}\left(\mathbf{X}'\mathbf{W}\mathbf{X}\right) = \operatorname {rk}\left(\mathbf{X}\right).

?
Задача 19.11

Пусть V\mathbf{V} — симметрическая неотрицательно определённая n×nn \times n-матрица, X\mathbf{X} — n×pn \times p-матрица, а d\mathbf{d} — p×1p \times 1-вектор. Покажите, что каждое из следующих шести условий необходимо и достаточно для того, чтобы вектор X(X′X)−d\mathbf{X}(\mathbf{X}'\mathbf{X})^{-}\mathbf{d} был решением, для каждого d∈C(X′)\mathbf{d} \in C(\mathbf{X}'), задачи минимизации квадратичной формы a′Va\mathbf{a}'\mathbf{V}\mathbf{a} (относительно a\mathbf{a}) при ограничении X′a=d\mathbf{X}'\mathbf{a} = \mathbf{d}:

?
(a)

C(VX)⊂C(X)C(\mathbf{V}\mathbf{X}) \subset C(\mathbf{X}) (или, что то же самое, VX=XQ\mathbf{V}\mathbf{X} = \mathbf{X}\mathbf{Q} для некоторой матрицы Q\mathbf{Q});

(b)

PXV(I−PX)=0\mathbf{P}_{\mathbf{X}}\mathbf{V}(\mathbf{I} - \mathbf{P}_{\mathbf{X}}) = \mathbf{0} (или, что то же самое, PXV=PXVPX\mathbf{P}_{\mathbf{X}}\mathbf{V} = \mathbf{P}_{\mathbf{X}}\mathbf{V}\mathbf{P}_{\mathbf{X}});

(c)

PXV=VPX\mathbf{P}_{\mathbf{X}}\mathbf{V} = \mathbf{V}\mathbf{P}_{\mathbf{X}} (или, что то же самое, PXV\mathbf{P}_{\mathbf{X}}\mathbf{V} симметрична);

(d)

C(VPX)⊂C(PX)C(\mathbf{V}\mathbf{P}_{\mathbf{X}}) \subset C(\mathbf{P}_{\mathbf{X}});

(e)

C(VPX)=C(V)∩C(PX)C(\mathbf{V}\mathbf{P}_{\mathbf{X}}) = C(\mathbf{V}) \cap C(\mathbf{P}_{\mathbf{X}});

(f)

C(VX)=C(V)∩C(X)C(\mathbf{V}\mathbf{X}) = C(\mathbf{V}) \cap C(\mathbf{X}).

Задача 19.12

Пусть V\mathbf{V} — симметрическая неотрицательно определённая n×nn \times n-матрица, W\mathbf{W} — n×nn \times n-матрица, X\mathbf{X} — n×pn \times p-матрица, а d\mathbf{d} — p×1p \times 1-вектор. Далее, пусть K\mathbf{K} — произвольная n×qn \times q-матрица такая, что C(K)=C(I−PX,W′)C(\mathbf{K}) = C(\mathbf{I} - \mathbf{P}_{\mathbf{X},\mathbf{W}'}). Покажите, что если rk⁡(X′WX)=rk⁡(X)\operatorname {rk}\left(\mathbf{X}'\mathbf{W}\mathbf{X}\right) = \operatorname {rk}\left(\mathbf{X}\right), то каждое из следующих двух условий необходимо и достаточно для того, чтобы вектор WX(X′WX)−d\mathbf{W}\mathbf{X}(\mathbf{X}'\mathbf{W}\mathbf{X})^{-}\mathbf{d} был решением, для каждого d∈C(X′)\mathbf{d} \in C(\mathbf{X}'), задачи минимизации квадратичной формы a′Va\mathbf{a}'\mathbf{V}\mathbf{a} (относительно a\mathbf{a}) при ограничении X′a=d\mathbf{X}'\mathbf{a} = \mathbf{d}:

?
(a)

V=XR1X′+(I−PX,W′)R2(I−PX,W′)′\mathbf{V} = \mathbf{X}\mathbf{R}_{1}\mathbf{X}' + (\mathbf{I} - \mathbf{P}_{\mathbf{X},\mathbf{W}'})\mathbf{R}_{2}(\mathbf{I} - \mathbf{P}_{\mathbf{X},\mathbf{W}'})'

для некоторой p×pp \times p-матрицы R1\mathbf{R}_{1} и некоторой n×nn \times n-матрицы R2\mathbf{R}_{2};

(b)

V=XS1X′+KS2K′\mathbf{V} = \mathbf{X}\mathbf{S}_{1}\mathbf{X}' + \mathbf{K}\mathbf{S}_{2}\mathbf{K}'

для некоторой p×pp \times p-матрицы S1\mathbf{S}_{1} и некоторой q×qq \times q-матрицы S2\mathbf{S}_{2}.

И покажите, что если rk⁡(X′WX)=rk⁡(X)\operatorname {rk}\left(\mathbf{X}'\mathbf{W}\mathbf{X}\right) = \operatorname {rk}\left(\mathbf{X}\right) и W\mathbf{W} невырожденна, то ещё одним необходимым и достаточным условием является: [(c)] V=tW−1+XT1X′+KT2K′\mathbf{V} = t\mathbf{W}^{-1} + \mathbf{X}\mathbf{T}_{1}\mathbf{X}' + \mathbf{K}\mathbf{T}_{2}\mathbf{K}'

для некоторого скаляра tt, некоторой p×pp \times p-матрицы T1\mathbf{T}_{1} и некоторой q×qq \times q-матрицы T2\mathbf{T}_{2}.

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

Чтобы установить необходимость условия (a), начните с того, что V=CC′\mathbf{V} = \mathbf{C}\mathbf{C}' для некоторой матрицы C\mathbf{C}, и представьте C\mathbf{C} в виде C=PX,W′C+(I−PX,W′)C\mathbf{C} = \mathbf{P}_{\mathbf{X},\mathbf{W}'}\mathbf{C} + (\mathbf{I} - \mathbf{P}_{\mathbf{X},\mathbf{W}'})\mathbf{C}.