ALS Invité
|
Posté le: 05 Juin 2006 8:08 Sujet du message: |
|
|
Bonjour, désolé, j'avais lu le sujet un peu trop rapidement. Il semble que ce problème n'ait pas de solution pour les matrices 2x2, car les solutions obtenues qui commutent toutes 2 avec AB-BA commutent aussi entre elles.
Code: |
> restart: with(linalg):A:=matrix(2,2): B:=matrix(2,2):
Warning, the protected names norm and trace have been redefined and unprotected
> C:=evalm(A&*B-B&*A):
> Y:=evalm(A&*C-C&*A): Z:=evalm(B&*C-C&*B):
> S:={seq(seq(Y[i,j],j=1..2),i=1..2),seq(seq(Z[i,j],j=1..2),i=1..2)}:
> solutions:={solve(S)}:
> for k to nops(solutions) do
> A:=matrix(2,2): B:=matrix(2,2): assign(solutions[k]): print(evalm(A),evalm(B),simplify(evalm(A&*B-B&*A)));
> end do:
[A[1, 1] A[1, 2]] [B[2, 2] 0 ] [0 0]
[ ], [ ], [ ]
[A[2, 1] A[2, 2]] [ 0 B[2, 2]] [0 0]
[A[2, 1] B[1, 1] + A[2, 2] B[2, 1] - B[2, 2] A[2, 1] ]
[--------------------------------------------------- , 0]
[ B[2, 1] ]
[A[2, 1] , A[2, 2]] [B[1, 1] 0 ] [0 0]
, [ ], [ ]
[B[2, 1] B[2, 2]] [0 0]
[A[1, 1] 0 ] [B[1, 1] 0 ] [0 0]
[ ], [ ], [ ]
[ 0 A[2, 2]] [ 0 B[2, 2]] [0 0]
[B[1, 2] A[2, 2] - A[1, 2] B[2, 2] + B[1, 1] A[1, 2]
[--------------------------------------------------- ,
[ B[1, 2]
]
A[1, 2]]
]
[A[1, 2] B[2, 1] ] [B[1, 1] B[1, 2]] [0 0]
[--------------- , A[2, 2]], [ ], [ ]
[ B[1, 2] ] [B[2, 1] B[2, 2]] [0 0]
|
Essayez de la même manière avec des matrices 3x3.
A plus tard. |
|