2014 dxdy logo

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

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


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


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



Начать новую тему Ответить на тему
 
 Линейная алгебра, минимизация нормы
Сообщение16.11.2016, 11:25 


27/09/16
8
Есть задача:
Матрицы $A$ и $B$ одинаковой размерности. Нужно найти такую ортогональную матрицу $Q$, что $||QA-B||_F$ минимальна. $F$ - норма Фробениуса.

Мое решение: найдем вместо минимума $||QA-B||_F$ минимум $||QA-B||^2_F$. Минимум будем находить с помощью метода множителей Лагранжа с условием $E - QQ^T = 0$

В итоге я получил $Q(2AA^T+L^T+L) = 2BA^T$. как отсюда выразить $Q$?пробовал просто транспонировать это равенство и перемножить их, получил $(2AA^T+L^T+L)^T(2AA^T+L^T+L) = (BA^T)^T(BA^T)$. Отсюда $Q = 2E$. А это какая-то ерунда! Помогите разобраться, пожалуйста.

 Профиль  
                  
 
 Re: Линейная алгебра, минимизация нормы
Сообщение16.11.2016, 14:56 
Заслуженный участник
Аватара пользователя


01/03/06
13626
Москва
stanger в сообщении #1169390 писал(а):
В итоге я получил $Q(2AA^T+L^T+L) = 2BA^T$

Давайте сначала предложим форумчанам догадаться, что такое $L$! Экстрасенсы, вперед! Форум в вас верит! :shock:

 Профиль  
                  
 
 Re: Линейная алгебра, минимизация нормы
Сообщение16.11.2016, 16:12 


27/09/16
8
$L$ - это собственно матрица, образованная из множителей Лагранжа.

То есть, вообще говоря, функция Лагранжа выглядит следующим образом:

$L(Q, \lambda) = ||QA-B||_F^2 - \sum_{ij}\lambda_{ij}(E-QQ^T)_{ij}$

Перепишем это в матричном виде $L = (\lambda)_{ij}$:

$L(Q, \lambda) = \operatorname{tr} ((QA-B)^T(QA-B))- \operatorname{tr} (L^T(E-QQ^T))$

Отсюда мы должны взять производную по $Q$ и приравнять ее к 0.
У меня получилось: $QAA^T + QAA^T - BA^T - BA^T + (L^T+L)Q = 0$
Что я делаю не так?

 Профиль  
                  
 
 Re: Линейная алгебра, минимизация нормы
Сообщение16.11.2016, 20:33 
Заслуженный участник
Аватара пользователя


01/03/06
13626
Москва
stanger в сообщении #1169466 писал(а):
Что я делаю не так?

Вы не так излагаете свои попытки решения. Чтобы вас проверить, кому-то нужно восстановить пропущенные вами куски. Это имеет смысл делать, когда я читаю научную статью. Но зачем мне так корячиться, чтобы найти вашу ошибку? Предложите мне хотя бы одну разумную причину для этого, и я сразу же кинусь к столу с ручкой и бумагой!

 Профиль  
                  
 
 Re: Линейная алгебра, минимизация нормы
Сообщение16.11.2016, 21:13 
Заслуженный участник
Аватара пользователя


06/10/08
6422
stanger в сообщении #1169466 писал(а):
У меня получилось: $QAA^T + QAA^T - BA^T - BA^T + (L^T+L)Q = 0$
Уравнение правильное. Непонятно, как Вы сделали вот этот вывод:
stanger в сообщении #1169390 писал(а):
получил $(2AA^T+L^T+L)^T(2AA^T+L^T+L) = (BA^T)^T(BA^T)$. Отсюда $Q = 2E$.

Из написанного Вами соотношения надо выражать $L$ (а точнее, какое-нибудь удобное подвыражение с участием $L$), а из него уже $Q$.

 Профиль  
                  
 
 Re: Линейная алгебра, минимизация нормы
Сообщение16.11.2016, 23:36 


27/09/16
8
я уже разобрался, спасибо.

данное равенство $QAA^T + QAA^T - BA^T - BA^T + (L^T+L)Q = 0$ я домножил на $Q^T$ справа. Получил $2AA^T + Q^T(L^T+L)Q = 2Q^TBA^T$. Затем транспонируем полученное равенство: $2AA^T+Q^T(L+L^T)Q = 2(BA^T)^TQ$. Вычтем из одного второе, получим, $Q^TBA^T = (BA^T)^TQ$. Ну отсюда уже и достанем $Q$

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

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



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

Сейчас этот форум просматривают: Geen


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

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