![]() |
|
![]() ![]() ![]() |
|
mrgloom |
|
|||
Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 829 Регистрация: 8.6.2011 Репутация: нет Всего: нет |
есть матричное уравнение, все матрицы 3х3, надо найти M, А и B известны.
пните в нужную сторону |
|||
|
||||
Pavia |
|
|||
Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 418 Регистрация: 6.12.2008 Репутация: 11 Всего: 12 |
A*M-M*B=0
Далее решаем систему уравнений C*D=0 D=(m00,m01,m02,m10,m11,m12,m20,m21,m22); |
|||
|
||||
mrgloom |
|
|||
Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 829 Регистрация: 8.6.2011 Репутация: нет Всего: нет |
а С чему равно?
|
|||
|
||||
W4FhLF |
|
|||
![]() found myself ![]() ![]() ![]() ![]() Профиль Группа: Участник Клуба Сообщений: 2831 Регистрация: 2.12.2006 Репутация: 5 Всего: 121 |
Pavia, далеко не всё так просто. Вообще по виду уравнение не имеет решений в общем случае, кроме тривиального М = 0.
Может это поможет: http://en.wikipedia.org/wiki/Sylvester_equation Хотя форма другая всё-таки. Это сообщение отредактировал(а) W4FhLF - 5.8.2012, 20:20 -------------------- "Бог умер" © Ницше "Ницше умер" © Бог |
|||
|
||||
entrix |
|
|||
Новичок Профиль Группа: Участник Сообщений: 3 Регистрация: 25.4.2010 Репутация: нет Всего: нет |
Смотрите GNU Octave, там есть соответствующая функция syl:
Solve the Sylvester equation A X + X B + C = 0 using standard lapack subroutines. For example: syl ([1, 2; 3, 4], [5, 6; 7, 8], [9, 10; 11, 12]) ⇒ [ -0.50000, -0.66667; -0.66667, -0.50000 ] http://www.gnu.org/software/octave/doc/int...#index-syl-2110 |
|||
|
||||
mrgloom |
|
|||
Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 829 Регистрация: 8.6.2011 Репутация: нет Всего: нет |
похоже, то что нужно.
попробую в Mathematica. |
|||
|
||||
mrgloom |
|
|||
Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 829 Регистрация: 8.6.2011 Репутация: нет Всего: нет |
попробовал в Mathematica выдало все нули, ну т.е. тривиальное решение.
может что то сделал не правильно? мне всё таки надо наверно сперва доказать, что существует решение отличное от тривиального. |
|||
|
||||
maxim1000 |
|
|||
![]() Эксперт ![]() ![]() ![]() ![]() Профиль Группа: Участник Сообщений: 3334 Регистрация: 11.1.2003 Где: Киев Репутация: 33 Всего: 110 |
это обычное однородное линейное уравнение, как и сказал Pavia
так что можно просто построить матрицу (9*9) и посмотреть на её определитель -------------------- qqq |
|||
|
||||
mrgloom |
|
|||
Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 829 Регистрация: 8.6.2011 Репутация: нет Всего: нет |
вы походу перепутали
A*X=B*X и A*X=X*B тут Х так просто не выразить Это сообщение отредактировал(а) mrgloom - 6.8.2012, 12:52 |
|||
|
||||
W4FhLF |
|
|||
![]() found myself ![]() ![]() ![]() ![]() Профиль Группа: Участник Клуба Сообщений: 2831 Регистрация: 2.12.2006 Репутация: 5 Всего: 121 |
maxim1000, умножение матриц операция не коммутативная.
-------------------- "Бог умер" © Ницше "Ницше умер" © Бог |
|||
|
||||
mrgloom |
|
||||||||||||||
Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 829 Регистрация: 8.6.2011 Репутация: нет Всего: нет |
http://math.stackexchange.com/questions/17...vester-equation
проверяю на существование не тривиального решения через первый критерий
выдает {} матрицы задавал как
так же пробовал
хотелось бы узнать как по матрице A
можно найти хотя бы 1 матрицу B которая будет удовлетворять вышеописанному условию? так же пробовал через другое эквивалентное условие.
хотя возможно оно выглядит так
вопрос опять же остается, хотелось бы узнать как по матрице A
можно найти хотя бы 1 матрицу B которая будет удовлетворять вышеописанному условию? |
||||||||||||||
|
|||||||||||||||
maxim1000 |
|
|||
![]() Эксперт ![]() ![]() ![]() ![]() Профиль Группа: Участник Сообщений: 3334 Регистрация: 11.1.2003 Где: Киев Репутация: 33 Всего: 110 |
а этого и не нужно просто записать каждый коеэффициент левого произведения и каждый коэффициент правого произведения, ну и приравнять 9 уравнений 9 неизвестных все уравнения линейные тот факт, что произведение некоммутативно, приведёт разве что к тому, что формулы будут разные и разные наборы неизвестных будут слева и справа, но для СЛАУ это не так важно -------------------- qqq |
|||
|
||||
mrgloom |
|
|||
Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 829 Регистрация: 8.6.2011 Репутация: нет Всего: нет |
дело в том, что эти 9 уравнений не решить как обычную СЛАУ.
т.к. хотя бы для примера AX-XB первая строчка зависит от x11,x12,x13,x21,x31 AX первая строчка зависит от x11,x21,x31 возможно AX-XB=0 как то и переписывается в виде типа P*X=C но как построить матрицы P,C ? |
|||
|
||||
maxim1000 |
|
|||
![]() Эксперт ![]() ![]() ![]() ![]() Профиль Группа: Участник Сообщений: 3334 Регистрация: 11.1.2003 Где: Киев Репутация: 33 Всего: 110 |
да почему??? даже если в каждом уравнении будут участвовать все 9 неизвестных, менее линейной система не станет на всякий случай уточню, систему будет представлена уже матрицей 9*9 (т.к. 9 неищвестных и 9 уравнений) Это сообщение отредактировал(а) maxim1000 - 7.8.2012, 10:06 -------------------- qqq |
|||
|
||||
mrgloom |
|
|||
Опытный ![]() ![]() Профиль Группа: Участник Сообщений: 829 Регистрация: 8.6.2011 Репутация: нет Всего: нет |
да, я кажется понял.
вместо системы 3х3*3х3 получаем 9х9*1х9 A*X-X*B=new_A*new_Xr=0 я это вручную сделал, а вот как записать формулу для матрицы 9х9 в общем виде? потом загнал в LeastSquares и получил опять тривиальное решение =0. |
|||
|
||||
![]() ![]() ![]() |
Правила форума "Алгоритмы" | |
|
Форум "Алгоритмы" предназначен для обсуждения вопросов, связанных только с алгоритмами и структурами данных, без привязки к конкретному языку программирования и/или программному продукту.
Если Вам понравилась атмосфера форума, заходите к нам чаще! С уважением, maxim1000. |
1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей) | |
0 Пользователей: | |
« Предыдущая тема | Алгоритмы | Следующая тема » |
|
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности Powered by Invision Power Board(R) 1.3 © 2003 IPS, Inc. |