Ich interessiere mich für folgendes Problem: Bei gegebenen ganzzahligen Matrizen
Dieses Mittel genau das, was Sie denken , es tut: Wir haben die Menge von Matrizen sagen
Wurde das Problem der Entscheidung, ob jedes Produkt irgendwann gleich Null ist, bereits untersucht? Ist es entscheidbar?
Es scheint, als könnte es mit der Matrixmortalität zusammenhängen, was nicht zu entscheiden ist, aber ich sehe keinen klaren Zusammenhang.
linear-algebra
decidability
Robinson
quelle
quelle
Antworten:
Ihre Frage ist äquivalent dazu, ob eine nilpotente Algebra erzeugenA1,…,Ak Ai O~(n2ω) ω
, die wiederum äquivalent dazu ist, dass jedes derAinilpotent ist. Daher ist es nicht nur entscheidbar, aber in ~ O ( n 2 ω ) die Zeit , wo ω der Exponent der Matrixmultiplikation ist.LetA be the associative algebra generated by the Ai : that is, take all linear combinations of the Ai and all finite products thereof. A is called nilpotent if there is some N such that every product of N elements of A is zero.
First, let's see why your condition implies thatA is nilpotent. This follows from Konig's Lemma (compactness): every string of length n over the alphabet {1,…,k} corresponds to a product of A1,…,Ak of length n in an obvious manner. Consider the infinite k -ary rooted tree, whose nodes are naturally in bijective correspondence with strings over {1,…,k} . Consider the sub-tree T consisting of those nodes where the corresponding product of the Ai is nonzero. Konig's Lemma says that if T is infinite, then it has an infinite path (exactly violating your property), hence T is finite. We can then take N to be the maximum length of any string in T . So your property implies that A is nilpotent.
The converse is also true, since every element ofA is a linear combination of products of the Ai .
Next, note thatA is a subalgebra of n×n matrices, and hence is finite-dimensional.
Finally: a finite-dimensional associative algebra in characteristic zero has a basis of nilpotent elements (commuting or not - this is the part that contradicts Yuval's answer) iff it is nilpotent (see, e.g., here).
Thus, to solve your problem, find a basis for the associative algebra generated by theAi (by the linear-algebra version of breadth-first search) and check that every matrix in the basis is nilpotent. The upper bound O~(n2ω) comes from solving a system of linear equations in n2 variables in the breadth-first search. As dimA≤n2 the BFS can't last very long, and because these are n×n matrices to check if a matrix A is nilpotent one needs only check that An=0 .
quelle
I got a poly-time algorithm for this (rather trivial problem)problem, i.e. for checking whether the joint spectral radius(JSR) is zero or not, in 1995: http://en.wikipedia.org/wiki/Joint_spectral_radius
The story behind the algorithm is roughly as follows: Blondel and Tsitsiklis wrongly stated that for boolean matrices checking whether JSR < 1 is NP-HARD. For any set of integer matrices JSR is ether zero or greater or equal 1. So the counter example to their statement was my algorithm(see the errata to their paper). The main moral: consult the Wikipedia first!
quelle
The question you are asking is exactly equivalent to deciding whether the joint spectral radius (JSR) of the set of matrices is strictly less than one. Decidability of this question has remained open for quite some time now. (In control theory, this is equivalent to decidability of stability of switched linear systems under arbitrary switching.)
The following variant of your question is known to be undecidable: Given a finite set of square matrices, decide whether all products remain bounded; see here.
The undecidability of the above remains valid even if you have only 2 matrices of size 47x47: see here.
In the JSR language, the question of testing "is JSR≤1 ?" is undecidable (see references above), but decidability of testing "is JSR <1 ?" is open.
The latter question is related to the so-called "rational finiteness conjecture":
If the rational finiteness conjecture is true, then the question you are asking is decidable.
Finally, unless P=NP, the JSR is not approximable in polynomial time (in the precise sense defined in this paper).
As a result, one of the answers above which claims an efficient algorithm must be false.
On the positive side, there are several algorithms (e.g. based on semidefinite programming) for approximating the JSR. The different algorithms come with different performance guarantees. See e.g. the following (shamelessly by myself and my colleagues - but see also references therein).
In several special cases, the question you are asking is polynomial time decidable. For example, when the matrices are symmetric, or rank one, or if they commute.
Finally, a great book on the subject is the following.
quelle
Edit: This answer is unfortunately incorrect. The error is highlighted below. The argument does work if we are allowed to transpose the matrices.
We start by proving a lemma.
Lemma. LetA be an n×n matrix and let N be the n×n matrix with ones on the secondary diagonal. If ANt and NtA are nilpotent for all t≥0 then A=0 . Correct conclusion: A is upper triangular with zeroes on the diagonal. (The original conclusion is recovered if we are also allowed to multiply by powers of the transpose of N .)
Proof. Suppose for example thatn=3 , and write
If we now considerN2A,N1A,N0A instead, then we conclude that A is lower triangular with zeroes on the diagonal. In fact, we don't get anything new from considering NtA . Therefore A=0 . □
This is how the proof would go if the original version of the lemma were correct. Now back to the problem at hand. Say that the matricesA1,…,Ak satisfy property P if for every infinite sequence i1,…∈[k] , we have Ai1⋯Aim=0 for some m . If one of the matrices Ai is not nilpotent then property P clearly fails, so suppose that all the matrices are nilpotent. If all matrices commute then property P clearly holds, so suppose that A1A2≠A2A1 . Change basis so that A1 is in Jordan normal form, and let the corresponding decomposition of the vector space be V1⊕⋯⊕Vt . Let Vi be a vector space on which A1A2≠A2A1 ; note that dimVi>1 since 0 commutes with everything. Restricted to Vi , A1=N and A2≠0 . Therefore the lemma implies that for some t≥0 , either A2At1 or At1A2 is not nilpotent, and therefore property P clearly fails.
Summarizing, property P holds iff all matrices are nilpotent and all of them commute.
quelle