Cómo resolver este tipo de ecuación matricial relacionada con el rango
Vi una publicación aquí para resolver la ecuación de relación de rango de la matriz:
symbolicMatrixRank[mat_, assumptions_] := Assuming[assumptions,
Simplify @ Total @ Map[
Boole @ Simplify @ Reduce[ConditionalExpression[#, $Assumptions] != 0]&,
SingularValueList[mat]
]
];
adj[m_](*adjoint matrix*):=
Map[Reverse, Minors[Transpose[m], Length[m] - 1], {0, 1}]*
Table[(-1)^(i + j), {i, Length[m]}, {j, Length[m]}]
Reduce[symbolicMatrixRank[adj[{{a, b, b}, {b, a, b}, {b, b, a}}],
Element[a, Reals] && Element[b, Reals]] == 1]
La función personalizada symbolicMatrixRank
de Sjoerd Smit es excelente. Pero encuentro que la siguiente ecuación matricial no se puede resolver con el método anterior:
Reduce[symbolicMatrixRank[{{a, 2*b}, {b, 2*c}, {c, 2*a}},
Element[a, Reals] && Element[b, Reals] && Element[c, Reals]] ==
symbolicMatrixRank[{{a, 2*b, -3*c}, {b, 2*c, -3*a}, {c,
2*a, -3*b}},
Element[a, Reals] && Element[b, Reals] && Element[c, Reals]] ==
2]
¿Cómo puedo mejorar el código anterior para poder resolver el tipo de ecuación matricial?
Nota: Las siguientes preguntas pertenecen a la décima pregunta del Examen de ingreso para graduados en matemáticas de China de 2003 (primer conjunto).
La respuesta de referencia de esta pregunta es a + b + c = 0
.
Respuestas
Solo logré obtener este simplemente jugando con los valores singulares hasta que alguna combinación de Reduce
y Simplify
produjo algo legible:
mat1 = {{a, 2*b}, {b, 2*c}, {c, 2*a}}
mat2 = {{a, 2*b, -3*c}, {b, 2*c, -3*a}, {c, 2*a, -3*b}}
assumptions = {a, b, c} \[Element] Reals
rank1 = Assuming[assumptions,
Total[
Boole @ Simplify @ Reduce @ Reduce[# != 0 && assumptions] & /@ SingularValueList[mat1]
]
]
rank2 = Assuming[assumptions,
Total[
Boole @ Simplify @ Reduce @ Reduce[# != 0 && assumptions] & /@ SingularValueList[mat2]
]
]
Simplify[Reduce[assumptions && rank1 == rank2], Assumptions -> assumptions]
(c! = 0 && ((a == c && b == c) || a + b + c == 0)) || (a + b == 0 && c == 0)