2014 dxdy logo

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

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




 
 Помогите решить матричное уравнение...
Сообщение02.04.2007, 03:12 
Уважаемые математики, в процессе работы возникла необходимость решить следующую задачку. Дано: матрица А (n x m, n>m) - известная; матрица C, такой же размерности - неизвестная. Необходимо решить уравнение:

$A^T C + C^T A + C^T C = \alpha E$

где $\alpha $ - скаляр; $E$ - единичная матрица, размерности m x m. Требуется, чтобы при решении или не использовались итерационные методы (типа Ньютона, релаксации и т.п.), или их сходимость была бы 100%. Понятно, что однозначно можно определить только m*m элементов матрицы C.

Заранее благодарен...

 
 
 
 
Сообщение03.04.2007, 02:45 
Аватара пользователя
:evil:
Может, это несколько упростит задачу: искать $X^T X = A^T A + \alpha E$, $C = X - A$.

Правильно ли я понимаю, что $n \geqslant m$?

 
 
 
 
Сообщение03.04.2007, 05:44 
Спасибо, идея мне понятна. Правда, как бы потом еще из $X^T X$
вычленить саму Х (точнее любую, удовлетворяющую уравнению)? Число строк действительно больше числа столбцов.

 
 
 
 
Сообщение03.04.2007, 19:18 
Аватара пользователя
:evil:
У Вас слева всегда положительно определенная симметричная матрица, справа симметричная, но вот положительно ли определенная? — это зависит от $\alpha$. При достаточно большом отрицательном $\alpha$ вещественного решения просто не существует.

Пусть все хорошо. Более того, пусть все собственные числа $A^T A$ различны (я не уверен, что это нужно). Тогда существует ортогональная $\Phi$, приводящая $A^T A$ к диагональному виду: $\Phi^T A^T A \, \Phi = \Lambda$. Тогда легко извлечь корень из $\Lambda + \alpha E$, и $X = \Phi \sqrt{\Lambda + \alpha E} \, \Phi^T$ должно быть одним из решений. (Хвост матрицы можно всегда добить нулями...)

P.S. Фактически, нам подходит любой алгорифм извлечения симметричного квадратного корня. Указанный метод лишь доказывает существование оного.

 
 
 
 
Сообщение04.04.2007, 00:01 
Суммарная матрица справа - всегда положительно определенная. А Вы не могли бы намекнуть, как найти Ф (хотя бы подсказать направление поисков)? Или вычисление Ф - составная часть алгоритмов извлечения корня? Т.е. нам фактически нужно просто извлечь корень из правой части?
Огромное Вам спасибо за помощь!

 
 
 
 
Сообщение04.04.2007, 00:58 
Аватара пользователя
:evil:
Например, система собственных векторов образует нужную нам матрицу (после нормирования). Так что, литературы довольно много. (Можно искать преобразование матрицы к диагональному виду, Жорданову нормальную форму, Jordan.) Я надеюсь, что меня кто-нибудь поправит, но, по-моему, правая часть является нормальной матрицей, и, соответственно, всегда диагонализируема согласно спектральной теореме. Обратите внимание на симметричность Вашей матрицы: для таких есть специализированные алгорифмы.

 
 
 
 
Сообщение04.04.2007, 01:05 
Направление понятно. Пошел искать :)

 
 
 [ Сообщений: 7 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group