Как использовать «решить» при сравнении двух матриц в MAXIMA? - PullRequest
1 голос
/ 04 апреля 2019

Я хочу сделать простую вещь:

M: matrix([a,b], [c,d]);
N: matrix([1,2], [3,4]);
solve(M = N, [a,b,c,d]);

но я получаю результат []. Если я разбью вышеупомянутое на много ограничений, как это:

M: matrix([a,b], [c,d]);
N: matrix([1,2],[3,4]);
c1: M[1,1] = N[1,1];
c2: M[1,2] = N[1,2];
c3: M[2,1] = N[2,1];
c4: M[2,2] = N[2,2];
solve([c1, c2, c3, c4], [M[1,1], M[1,2], M[2,1], M[2,2]]);

Я получаю хороший результат [[a=1,b=2,c=3,d=4]]. Но этот метод будет болезненным для больших матриц. Есть ли какой-нибудь более простой способ сделать это?

1 Ответ

2 голосов
/ 06 апреля 2019

Преобразование матриц в список может помочь

(%i1) M: matrix([a,b], [c,d])$
(%i2) N: matrix([1,2] [c,d])$
(%i3) f(M):= apply(join, args(M))$
(%i4) e: map("=",  f(M), f(N));
(%o4)                    [a = 1, c = 3, b = 2, d = 4]
(%i5) solve(e, [a, b, c, d]);
(%o4)                    [a = 1, c = 3, b = 2, d = 4]
...