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 ] 

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



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

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


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

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