Estou interessado no seguinte problema: dadas as matrizes inteiras decidem se todo produto infinito dessas matrizes eventualmente é igual à matriz zero.
Isso significa exatamente o que você pensa: diremos que o conjunto de matrizes tem a propriedade de que todos os seus produtos acabam sendo iguais a zero se não existir uma sequência infinita , todos em , de modo que para todos os .i 1 , i 2 , i 3 … { 1 , … , k } A i 1 A i 2 ⋯ A i l ≠ 0 l
O problema de decidir se todo produto eventualmente é igual a zero foi estudado antes? É decidível?
Parece que isso pode estar relacionado à mortalidade da matriz, que é indecidível, mas não vejo uma conexão clara.
linear-algebra
decidability
Robinson
fonte
fonte
Respostas:
Sua pergunta é equivalente a gerar uma álgebra nilpotenteUMA1, … , Ak UMAEu O~( n2 ω) ω
, que por sua vez é equivalente a cada um dosnilpotentes. Portanto, não é apenas decidível, mas em tempo em que é o expoente da multiplicação de matrizes.Ai˜ O ( n 2 ω ) ωSeja a álgebra associativa gerada pelo : ou seja, tome todas as combinações lineares do e todos os seus produtos finitos. é chamado nilpotent se houver algum tal que todo produto de elementos de seja zero.Um i A i A N N UmaUMA UMAEu UMAEu UMA N N UMA
Primeiro, vamos ver por que sua condição implica que é nilpotente. Isto decorre de Konig Lema (compacidade): cada cadeia de comprimento n sobre o alfabeto { 1 , ... , k } corresponde a um produto de A 1 , ... , A k de comprimento n de uma forma óbvia. Considere a infinita árvore k- arraigada, cujos nós estão naturalmente em correspondência bijetiva com cadeias acima de { 1 , … , k } . Considere a subárvore TUMA n { 1 , … , k } UMA1, … , Ak n k { 1 , … , k } T que consiste naqueles nós em que o produto correspondente da é diferente de zero. O Lema de Konig diz que se T é infinito, então ele tem um caminho infinito (violando exatamente sua propriedade), portanto, T é finito. Podemos, então, tomar N ser o comprimento máximo de qualquer cadeia em T . Portanto, sua propriedade implica que A é nilpotente.UMAEu T T N T UMA
O inverso também é verdadeiro, pois cada elemento de é uma combinação linear de produtos de A i .UMA UMAEu
Em seguida, observe que é uma subalgebra de n × n matrizes e, portanto, é de dimensão finita.UMA n × n
Finalmente: uma álgebra associativa de dimensão finita na característica zero tem uma base de elementos nilpotentes (pendulares ou não - esta é a parte que contradiz a resposta de Yuval) se for nilpotente (veja, por exemplo, aqui ).
Assim, para resolver o problema, encontrar uma base para a álgebra associativa gerado pelo (pela versão linear-álgebra de procura em largura) e verificar que cada matriz na base é nilpotent. O limite superior ˜ O ( n 2 ω ) vem da solução de um sistema de equações lineares em n 2 variáveis na busca pela primeira vez. Como dim A ≤ n 2, o BFS não pode durar muito, e como essas são n × n matrizes para verificar se uma matriz A é nilpotente, é necessário apenas verificar se A n =UMAEu O~( n2 ω) n2 escuroUMA≤ n2 n × n UMA .UMAn= 0
fonte
Eu obtive um algoritmo poli-tempo para esse problema (problema bastante trivial), ou seja, para verificar se o raio espectral da junta (JSR) é zero ou não, em 1995: http://en.wikipedia.org/wiki/Joint_spectral_radius
A história por trás do algoritmo é mais ou menos a seguinte: Blondel e Tsitsiklis afirmaram erroneamente que para matrizes booleanas verificando se JSR <1 é NP-HARD. Para qualquer conjunto de matrizes inteiras, o JSR é igual a zero ou maior ou igual a 1. Portanto, o exemplo de contador para a declaração deles foi o meu algoritmo (veja as erratas no artigo). A principal moral: consulte a Wikipedia primeiro!
fonte
A pergunta que você está fazendo é exatamente equivalente a decidir se o raio espectral da junta (JSR) do conjunto de matrizes é estritamente menor que um. A decidibilidade dessa questão permanece em aberto há algum tempo. (Na teoria do controle, isso é equivalente à decidibilidade da estabilidade de sistemas lineares comutados sob comutação arbitrária.)
Sabe-se que a seguinte variante de sua pergunta é indecidível: Dado um conjunto finito de matrizes quadradas, decida se todos os produtos permanecem delimitados; veja aqui .
A indecisão do acima exposto permanece válida mesmo se você tiver apenas 2 matrizes de tamanho 47x47: veja aqui .
Na linguagem JSR, a questão do teste "é JSR ?" é indecidível (consulte as referências acima), mas a decidibilidade do teste "é JSR < 1 ?" está aberto. A última pergunta está relacionada à chamada "conjectura racional da finitude": se a conjectura racional da finitude é verdadeira, então a pergunta que você está fazendo é decidível.≤ 1 < 1
Finalmente, a menos que P = NP, o JSR não seja aproximado em tempo polinomial (no sentido preciso definido neste artigo ).
Como resultado, uma das respostas acima, que afirma que um algoritmo eficiente deve ser falso.
No lado positivo, existem vários algoritmos (por exemplo, baseados em programação semidefinida) para aproximar o JSR. Os diferentes algoritmos vêm com diferentes garantias de desempenho. Veja, por exemplo, o seguinte (sem vergonha, por mim e pelos meus colegas - mas veja também as referências aqui ).
Em vários casos especiais, a pergunta que você está fazendo é um tempo polinomial decidível. Por exemplo, quando as matrizes são simétricas, ou classificam uma, ou se comutam.
Finalmente, um ótimo livro sobre o assunto é o seguinte .
fonte
Editar: infelizmente esta resposta está incorreta. O erro é destacado abaixo. O argumento funciona se nos for permitido transpor as matrizes.
Começamos provando um lema.
Lema. Deixe ser um n × n matriz e deixar que N seja o n x n matriz com aqueles na diagonal secundário. Se A n t e N t A são nilpotentes para todos os t ≥ 0, então A = 0 . Conclusão correta: A é triangular superior com zeros na diagonal. (A conclusão original é recuperada se também nos é permitido multiplicar por potências de transposição de N. )A n×n N n×n ANt NtA t≥0 A=0 A N
Prova. Suponha, por exemplo, que e escreva A = ( a b c d e f g h i ) ,n=3
Começamos calculando A N 2 :
A N 2 = ( 0 0 a 0 0 d 0 0 g ) .
Essa matriz está na forma triangular e, portanto, se A N 2 é nilpotente, então g = 0 . Continue com A N 1 :
A N 1 = ( 0
Se considerarmos agora , concluímos que A é triangular inferior com zeros na diagonal. Na verdade, nós não temos nada de novo a partir considerando N t A . Portanto A = 0 . ◻N2A,N1A,N0A A NtA A=0 □
Resumindo, a propriedade P mantém-se se todas as matrizes forem nilpotentes e todas elas comutarem.
fonte