Поиск:

Ответ в темуСоздание новой темы Создание опроса
> матричное уравнение 
:(
    Опции темы
mrgloom
Дата 24.7.2012, 14:18 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 829
Регистрация: 8.6.2011

Репутация: нет
Всего: нет



есть матричное уравнение, все матрицы 3х3, надо найти M, А и B известны. 
Цитата

A*M=M*B

пните в нужную сторону 
PM MAIL   Вверх
Pavia
Дата 24.7.2012, 18:17 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 418
Регистрация: 6.12.2008

Репутация: 11
Всего: 12



A*M-M*B=0
Далее решаем систему уравнений

C*D=0
D=(m00,m01,m02,m10,m11,m12,m20,m21,m22);
PM MAIL   Вверх
mrgloom
Дата 25.7.2012, 08:50 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 829
Регистрация: 8.6.2011

Репутация: нет
Всего: нет



а С чему равно?
PM MAIL   Вверх
W4FhLF
Дата 5.8.2012, 20:13 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


found myself
****


Профиль
Группа: Участник Клуба
Сообщений: 2831
Регистрация: 2.12.2006

Репутация: 5
Всего: 121



Pavia, далеко не всё так просто. Вообще по виду уравнение не имеет решений в общем случае, кроме тривиального М = 0. 

Может это поможет: http://en.wikipedia.org/wiki/Sylvester_equation

Хотя форма другая всё-таки. 

Это сообщение отредактировал(а) W4FhLF - 5.8.2012, 20:20


--------------------
"Бог умер" © Ницше
"Ницше умер" © Бог
PM ICQ   Вверх
entrix
Дата 6.8.2012, 02:35 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 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
PM MAIL   Вверх
mrgloom
Дата 6.8.2012, 11:15 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 829
Регистрация: 8.6.2011

Репутация: нет
Всего: нет



похоже, то что нужно.
попробую в Mathematica.
PM MAIL   Вверх
mrgloom
Дата 6.8.2012, 12:12 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 829
Регистрация: 8.6.2011

Репутация: нет
Всего: нет



попробовал в Mathematica выдало все нули, ну т.е. тривиальное решение.

может что то сделал не правильно?

мне всё таки надо наверно сперва доказать, что существует решение отличное от тривиального.


PM MAIL   Вверх
maxim1000
Дата 6.8.2012, 12:21 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
****


Профиль
Группа: Участник
Сообщений: 3334
Регистрация: 11.1.2003
Где: Киев

Репутация: 33
Всего: 110



это обычное однородное линейное уравнение, как и сказал Pavia

так что можно просто построить матрицу (9*9) и посмотреть на её определитель



--------------------
qqq
PM WWW   Вверх
mrgloom
Дата 6.8.2012, 12:51 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 829
Регистрация: 8.6.2011

Репутация: нет
Всего: нет



вы походу перепутали 
A*X=B*X
и 
A*X=X*B

тут Х так просто не выразить

Это сообщение отредактировал(а) mrgloom - 6.8.2012, 12:52
PM MAIL   Вверх
W4FhLF
Дата 6.8.2012, 13:34 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


found myself
****


Профиль
Группа: Участник Клуба
Сообщений: 2831
Регистрация: 2.12.2006

Репутация: 5
Всего: 121



maxim1000, умножение матриц операция не коммутативная.


--------------------
"Бог умер" © Ницше
"Ницше умер" © Бог
PM ICQ   Вверх
mrgloom
Дата 6.8.2012, 15:31 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 829
Регистрация: 8.6.2011

Репутация: нет
Всего: нет



http://math.stackexchange.com/questions/17...vester-equation


проверяю на существование не тривиального решения
через первый критерий
Цитата
NullSpace[KroneckerProduct[IdentityMatrix[3],A]-KronekerProduct[B,IdentityMatrix[3]]]

выдает {}
матрицы задавал как 
Цитата
A={{a11,a12,a13},{a21,a22,a23},{a31,a32,a33}}
B={{b11,b12,b13},{b21,b22,b23},{b31,b32,b33}}

так же пробовал
Цитата
A={{a11,a12,a13},{a21,a22,a23},{0,0,1}}
B={{b11,b12,b13},{b21,b22,b23},{b31,b32,1}}


хотелось бы узнать как по матрице A
Цитата
A={{a11,a12,a13},{a21,a22,a23},{0,0,1}}

можно найти хотя бы 1 матрицу B  которая будет удовлетворять вышеописанному условию?
 
так же пробовал через другое эквивалентное условие.
Цитата
Resultant[Det[A-x IdentityMatrix[3]],Det[-B-x IdentityMatrix[3]],x] =0

хотя возможно оно выглядит так 
Цитата
Resultant[Det[A-x IdentityMatrix[3]],Det[B-x IdentityMatrix[3]],x]=0

вопрос опять же остается, 
хотелось бы узнать как по матрице A
Цитата
A={{a11,a12,a13},{a21,a22,a23},{0,0,1}}

можно найти хотя бы 1 матрицу B которая будет удовлетворять вышеописанному условию?
PM MAIL   Вверх
maxim1000
Дата 6.8.2012, 18:18 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
****


Профиль
Группа: Участник
Сообщений: 3334
Регистрация: 11.1.2003
Где: Киев

Репутация: 33
Всего: 110



Цитата(W4FhLF @  6.8.2012,  13:34 Найти цитируемый пост)
maxim1000, умножение матриц операция не коммутативная. 

а этого и не нужно
просто записать каждый коеэффициент левого произведения и каждый коэффициент правого произведения, ну и приравнять

9 уравнений
9 неизвестных
все уравнения линейные

тот факт, что произведение некоммутативно, приведёт разве что к тому, что формулы будут разные и разные наборы неизвестных будут слева и справа, но для СЛАУ это не так важно


--------------------
qqq
PM WWW   Вверх
mrgloom
Дата 7.8.2012, 09:34 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 829
Регистрация: 8.6.2011

Репутация: нет
Всего: нет



дело в том, что эти 9 уравнений не решить как обычную СЛАУ.
т.к.  хотя бы для примера
AX-XB
первая строчка зависит от x11,x12,x13,x21,x31
AX
первая строчка зависит от x11,x21,x31


возможно 
AX-XB=0
как то и переписывается в виде типа  
P*X=C
но как построить матрицы P,C ?
PM MAIL   Вверх
maxim1000
Дата 7.8.2012, 10:01 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
****


Профиль
Группа: Участник
Сообщений: 3334
Регистрация: 11.1.2003
Где: Киев

Репутация: 33
Всего: 110



Цитата(mrgloom @  7.8.2012,  09:34 Найти цитируемый пост)
дело в том, что эти 9 уравнений не решить как обычную СЛАУ.
т.к.  хотя бы для примера
AX-XB
первая строчка зависит от x11,x12,x13,x21,x31
AX
первая строчка зависит от x11,x21,x31

да почему???
даже если в каждом уравнении будут участвовать все 9 неизвестных, менее линейной система не станет

на всякий случай уточню, систему будет представлена уже матрицей 9*9 (т.к. 9 неищвестных и 9 уравнений)

Это сообщение отредактировал(а) maxim1000 - 7.8.2012, 10:06


--------------------
qqq
PM WWW   Вверх
mrgloom
Дата 7.8.2012, 10:53 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 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.
PM MAIL   Вверх
Ответ в темуСоздание новой темы Создание опроса
Правила форума "Алгоритмы"

maxim1000

Форум "Алгоритмы" предназначен для обсуждения вопросов, связанных только с алгоритмами и структурами данных, без привязки к конкретному языку программирования и/или программному продукту.


Если Вам понравилась атмосфера форума, заходите к нам чаще! С уважением, maxim1000.

 
1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Алгоритмы | Следующая тема »


 




[ Время генерации скрипта: 0.0819 ]   [ Использовано запросов: 20 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.