2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки


Правила форума


Посмотреть правила форума



Начать новую тему Ответить на тему
 
 Доказать эквивалентность минимизации матричных выражений
Сообщение18.11.2023, 18:34 
Аватара пользователя


18/11/23
3
Здравствуйте,
помогите пожалуйста найти путь решения

$H \in \mathbb{R}^{n \times k}$ матрица с ортонормированными столбцами, $H^\top H=I$, $S(X)=X X^\top$

Доказать, что

$$\min_{\Theta \in \mathbb{R}^{k \times 3}} \|Z-S(H \Theta)\|_F^2
$$
эквивалентно
$$
\min_{\Theta \in \mathbb{R}^{k \times 3}} \left\|H^T Z H-\Theta \Theta^T\right\|_F^2.
$$

Я пробовал использовать свойство ортогональности матрицы $H$ и доказал эквивалентность для случая, когда $H$ квадратная, но для произвольной формы $H$ не знаю что делать, пробовал преобразовать выражение, но в таком случае множители $H$ и $H^T$ не сокращаются...

Натолкните на мысль пожалуйста :roll:

 Профиль  
                  
 
 Re: Доказать эквивалентность минимизации матричных выражений
Сообщение19.11.2023, 02:12 
Заслуженный участник
Аватара пользователя


23/07/08
10908
Crna Gora
Пусть
$E_m$ — единичная квадратная матрица порядка $m$,
$O_{m,p}$ — нулевая матрица размера $m\times p$.
Ещё определим матричные функции
$A(\Theta)=Z-H \Theta\Theta^\top H^\top$
$B(\Theta)=H^\top Z H-\Theta \Theta^\top$
Ясно, что $k\leqslant n$, иначе равенство $H^\top H=E_k$ невозможно. Раз у Вас всё получилось для $k=n$, давайте теперь считать, что $k<n.$

Рассмотрим простой случай, когда $Z=E_n$, а $H$ — блочная матрица
$H=\begin{bmatrix}E_k\\O_{n-k,k}\end{bmatrix}$
(условие $H^\top H=E_k$ при этом выполняется). Тогда
$A(\Theta)=\begin{bmatrix}E_k-\Theta \Theta^\top&O_{k,n-k}\\O_{n-k,k}&E_{n-k} \end{bmatrix}$
$B(\Theta)=E_k-\Theta \Theta^\top$

Но ведь это означает, что
$\|A(\Theta)\|_F^2=\|B(\Theta)\|_F^2+(n-k)$
:shock:

 Профиль  
                  
 
 Re: Доказать эквивалентность минимизации матричных выражений
Сообщение19.11.2023, 15:23 
Аватара пользователя


18/11/23
3
svv в сообщении #1618654 писал(а):
Пусть
$E_m$ — единичная квадратная матрица порядка $m$,
$O_{m,p}$ — нулевая матрица размера $m\times p$.
Ещё определим матричные функции
$A(\Theta)=Z-H \Theta\Theta^\top H^\top$
$B(\Theta)=H^\top Z H-\Theta \Theta^\top$
Ясно, что $k\leqslant n$, иначе равенство $H^\top H=E_k$ невозможно. Раз у Вас всё получилось для $k=n$, давайте теперь считать, что $k<n.$

Рассмотрим простой случай, когда $Z=E_n$, а $H$ — блочная матрица
$H=\begin{bmatrix}E_k\\O_{n-k,k}\end{bmatrix}$
(условие $H^\top H=E_k$ при этом выполняется). Тогда
$A(\Theta)=\begin{bmatrix}E_k-\Theta \Theta^\top&O_{k,n-k}\\O_{n-k,k}&E_{n-k} \end{bmatrix}$
$B(\Theta)=E_k-\Theta \Theta^\top$

Но ведь это означает, что
$\|A(\Theta)\|_F^2=\|B(\Theta)\|_F^2+(n-k)$
:shock:


Спасибо большое!
Еще в процессе ночных блужданий родилась идея воспользоваться ортогональностью оператора $H^T \cdot H$, что это сохраняет норму Фробениуса.

 Профиль  
                  
 
 Re: Доказать эквивалентность минимизации матричных выражений
Сообщение19.11.2023, 16:51 
Заслуженный участник
Аватара пользователя


23/07/08
10908
Crna Gora
student0431, чтобы я не беспокоился: Вы поняли, что я привёл контрпример к утверждению, которое надо доказать? В моём примере минимумы разные получаются. Отличаются на $n-k$. И ничего с этим не сделаешь. А если $n=k$, то всё в порядке.

 Профиль  
                  
 
 Re: Доказать эквивалентность минимизации матричных выражений
Сообщение21.11.2023, 16:36 
Аватара пользователя


18/11/23
3
svv в сообщении #1618771 писал(а):
student0431, чтобы я не беспокоился: Вы поняли, что я привёл контрпример к утверждению, которое надо доказать? В моём примере минимумы разные получаются. Отличаются на $n-k$. И ничего с этим не сделаешь. А если $n=k$, то всё в порядке.


Уважаемый svv,

Да, про разные минимумы довольно наглядно получилось показать, но все же, они же могут достигаться при одном и том же $\Theta=\Theta_{\min}$

 Профиль  
                  
 
 Re: Доказать эквивалентность минимизации матричных выражений
Сообщение22.11.2023, 04:53 
Заслуженный участник
Аватара пользователя


23/07/08
10908
Crna Gora
Если нужен не сам минимум, а аргумент (в нашем случае $\Theta$), при котором он достигается, это обозначается $\operatorname{arg\;min}$. См. здесь.
Значит, Вам достаточно, если минимумы этих функций достигаются на одной и той же матрице $\Theta$.

Пусть $m=n-k$. Рассмотрим столбцы $H \in \mathbb{R}^{n \times k}$ как векторы в $\mathbb{R}^n$. Используя процесс ортогонализации Грама-Шмидта, дополним их набор до ортонормированного базиса в $\mathbb{R}^n$. Получим ортогональную матрицу $P\in \mathbb{R}^{n \times n}$:
$P=\begin{bmatrix}P_1&P_2\end{bmatrix},$
где $P_1\in \mathbb{R}^{n \times k}$ — просто другое обозначение для $H$, а $P_2\in \mathbb{R}^{n \times m}$ — матрица из добавленных векторов.
Из $P^\top P=E_n$ следует
$\begin{array}{ll}P_1^\top P_1=E_k&P_1^\top P_2=O_{km} \\P_2^\top P_1=O_{mk}&P_2^\top P_2=E_m\end{array}$

Пусть $C(\Theta)=P^\top A(\Theta) P$. Поскольку $P$ ортогональна, $\|C(\Theta)\|=\|A(\Theta)\|$. Имеем:
$C(\Theta)=P^\top (Z-H \Theta\Theta^\top H^\top) P=\begin{bmatrix}P_1^\top\\P_2^\top\end{bmatrix}(Z-P_1 \Theta\Theta^\top P_1^\top)\begin{bmatrix}P_1&P_2\end{bmatrix}=$
$=\begin{bmatrix}P_1^\top Z P_1-\Theta\Theta^\top&P_1^\top Z P_2\\P_2^\top Z P_1&P_2^\top Z P_2\end{bmatrix}=\begin{bmatrix}B(\Theta)&P_1^\top Z P_2\\P_2^\top Z P_1&P_2^\top Z P_2\end{bmatrix}$

Поскольку квадрат нормы Фробениуса матрицы — это сумма квадратов её элементов, то
$\|A(\Theta)\|^2=\|C(\Theta)\|^2=\|B(\Theta)\|^2+\|P_1^\top Z P_2\|^2+\|P_2^\top Z P_1\|^2+\|P_2^\top Z P_2\|^2$
Но последние три слагаемых не зависят от $\Theta$. Поэтому минимум $\|A(\Theta)\|$ и $\|B(\Theta)\|$ достигается при одной и той же матрице $\Theta$ (точнее, на одном и том же множестве матриц).

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 6 ] 

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group