Bir matris kümesinin tüm ürünlerinin sonunda sıfıra eşit olup olmadığını kontrol etme


19

Aşağıdaki sorunla ilgileniyorum: tamsayı matrisleri verildiğinde A1,A2,,Ak, bu matrislerin her sonsuz ürününün sonunda sıfır matrise eşit olup olmadığına karar verir.

Bu araçlar tam olarak öyle düşündüğünü: Biz matrislerin kümesi diyecekler {A1,,Ak} ürünlerinin sonunda eşit sıfır tüm orada eğer o özelliğine sahiptir değil sonsuz bir dizi var i1,i2,i3 , hepsi {1,,k} , öyle ki

Ai1Ai2Ail0
herkes için l .

Her ürünün nihayetinde sıfıra eşit olup olmadığına karar verme sorunu daha önce incelenmiş midir? Karar verilebilir mi?

Kararsız olan matris mortalitesi ile ilgili gibi görünebilir, ancak net bir bağlantı görmüyorum.


Sonsuz ürünün tanımlandığından emin olmak için matris setinde bir tür yakınsama özelliğine ihtiyacınız vardır.
András Salamon

Sınırlı bir alanda mı yoksa sınırsız büyüme gösteren tamsayılarda mı çalışıyorsunuz? = 1 vaka vardır kendi başlarına ilginç. 5x5 matrisinde -100'den 100'e kadar tamsayılar kullanarak sıfırlamadan önce alabileceğiniz en yüksek güç nedir? k
Chad Brewbaker

2
@YuvalFilmus - Ölümden farklı olduğuna inanıyorum. Matrislerin boyutları , sadece sayılarımız olsun ve A 0 = 0 , A 1 = 1 olduğunu varsayalım . Mortal? Evet çünkü A 0 = 0 . Her ürün sıfıra eşit mi? Hayır: ürün değil 1 1 1 . Öte yandan, A 0 = 0 , A 1 = 0 ise, hem ölümlü hem de her ürün sıfır olan bir diziniz vardır. 1A0=0,A1=1A0=0111A0=0,A1=0
Robinson

1
@ChadBrewbaker - Matrislerin girişlerinin sadece tamsayı olduğunu düşünüyordum. Sanırım , bakış açısından ilginç: matrisin nilpotent olduğunu kontrol etmek için kaç işlem yapmanız gerekiyor? Eğer Not A nilpotenttir, o zaman kolay olduğunu görmek için bir n = 0 , n boyutudur A böylece muhtemelen sen matris karesi alınarak çözebilir günlüğünü n kere. Yapabileceğin en iyi şey olup olmadığı hakkında hiçbir fikrim yok. k=1AAn=0nAlogn
Robinson

1
İlginçtir, bu sadece: arxiv.org/abs/1306.0729 . Tüm ürünlerin sonunda sıfır olup olmadığını sormak yerine, bazı ürünlerin sonunda olumlu olup olmadığını soruyorlar. Sorunun NP-zor olduğunu gösteriyorlar (ya da en azından özetten topladığım şey).
Joshua Grochow

Yanıtlar:


17

Sorunuz olmadığını eşdeğerdir bir nilpotentlik cebir üretmek sırayla her birine denktir, A i nilpotenttir olmak . Bu nedenle sadece karar verilebilir değil, aynı zamanda ω matris çarpımının üssü olan ˜ O ( n 2 ω ) zamanında .A1,,AkAiO~(n2ω)ω

Let tarafından oluşturulan birleştirici cebri olarak A ı tüm doğrusal kombinasyonlarını almak olup: A i ve bunların tüm sonlu ürünler. Bir denir nilpotentlik bazı varsa , N , her ürün, öyle ki , N unsurları A sıfırdır.AAiAiANNA

İlk olarak, durumunuzun neden nilpotent olduğunu ima ettiğini görelim . Bu Konig Lemma'sı (kompakt) dan şöyle: uzunluk her dize n alfabenin üzerine { 1 , ... , k } bir ürüne karşılık gelir A 1 , ... , A k uzunluğunun n bariz bir şekilde. Düğümleri doğal olarak { 1 , , k } üzerindeki dizelerle iki yönlü yazışma içinde olan sonsuz k -kök köklü ağacı düşünün . T alt ağacını düşününAn{1,,k}A1,,Aknk{1,,k}Tkarşılık gelen ürün bu düğümden oluşan sıfırdan farklıdır. Konig'den Lemma, eğer T sonsuzsa, o zaman sonsuz bir yola sahip olduğunu (mülkünüzü tam olarak ihlal ettiğini), bu nedenle T'nin sonlu olduğunu söylüyor . Daha sonra N'yi T'deki herhangi bir dizenin maksimum uzunluğu olarak alabiliriz . Yani mülkünüz A'nın nilpotent olduğunu ima ediyor .AiTTNTA

Her bir elemanı yana tersi de doğrudur ürünlerinin lineer bir kombinasyonu olan A i .AAi

Daha sonra, n × n matrislerinin bir alt cebiri ve dolayısıyla sonlu boyutlu olduğuna dikkat edin .An×n

Son olarak: karakteristik sıfırdaki sonlu boyutlu çağrışımsal bir cebir nilpotent elementlerin temeline sahiptir (gidip gelme - bu Yuval'ın cevabı ile çelişen kısımdır) nilpotent ise (bakınız, örneğin, burada ).

Böylece, probleminizi çözmek için, (en genişlik ilk aramanın lineer-cebir versiyonu ile) tarafından oluşturulan ilişkisel cebir için bir temel bulun ve temeldeki her matrisin nilpotent olduğunu kontrol edin. Üst sınır ˜ O ( n 2 ω ) , ilk genişlik araştırmasında n 2 değişkeninde doğrusal denklemler sisteminin çözülmesinden gelir . Olarak loş An 2 , bu çünkü BFS değil, son derece uzun ve olabilir , n x n bir matris kontrol etmek amacıyla matrisler bir olan nilpotentlik bir ihtiyacı sadece kontrol A , n =AiO~(n2ω)n2dimAn2n×nA .An=0


2
Bunu herhangi bir seçim prensibi kullanmadan göstermenin bir yolu olduğunu düşünüyor musunuz ( eşdeğer König'in Lemma'sı kadar zayıf olsa bile )? ACω
András Salamon

2
@Andras: Bunun Chris Conidis için bir soru olduğunu söyleyebilirim. (Hesaplanabilir) ters matematikte böyle sorular okudu. Ona soracağım ve onu burada göstereceğim.
Joshua Grochow

1
@robinson: 1) Evet sorun, aslında, Karar verilebilen bir zaman ω matris çarpım üssüdür. Bu , lineer cebir genişliği ilk araması yapılırken Q üzerinden lineer denklem sistemlerinin çözülmesinden gelir . 2) Evet, taban rotasyonunu vektör olarak matrisler görüntülerken Q , n , 2 (ya da üzerinde R ve C ). O(n2ω)ωQQn2RC
Joshua Grochow

1
Bir baz ile başlayan içinde A . Şimdi A A ve B B matrislerini bulmaya çalışıyorsunuz, böylece A B veya B A B aralığında olmayacak . Başarılı olursa ürünü B'ye ekleyin ve devam edin. Aksi takdirde, B'deki herhangi bir matrisin A'daki herhangi bir matrisin sonlu ürünü ile çarpılması her zaman B'nin aralığında olur . Cebirin boyutu sınırlandırıldığı için işlem sona erer (en fazla n 2 adımda). BAAABBABBABBBABn2
Yuval Filmus

1
@robinson: Hayır. Cebir nilpotent ise, cebirin her elemanı nilpotenttir. Eğer herhangi bir nilpotent olmayan bir eleman bulursanız, cebir nilpotent değildir (ve sonra matrislerinizin asla sıfır olmayan sonsuz ürünleri vardır).
Joshua Grochow

6

1995'te bu (oldukça önemsiz bir problem) problemi için bir poli-zaman algoritması aldım, yani eklem spektral yarıçapının (JSR) sıfır olup olmadığını kontrol etmek için: http://en.wikipedia.org/wiki/Joint_spectral_radius

Algoritmanın arkasındaki hikaye kabaca şu şekildedir: Blondel ve Tsitsiklis yanlış bir şekilde boole matrisleri için JSR <1'in NP-HARD olup olmadığını kontrol ettiğini belirtmiştir. Herhangi bir tamsayı matrisi kümesi için JSR eter sıfır ya da daha büyük ya da 1'e eşittir. Bu yüzden ifadelerine karşı olan örnek algoritmamdı (kağıtlarına bakınız). Ana ahlak: önce Wikipedia'ya danışın!


5

Sorduğunuz soru, matris kümesinin eklem spektral yarıçapının (JSR) kesinlikle birinden daha az olup olmadığına karar vermekle tamamen eşdeğerdir. Bu sorunun karar verilebilirliği bir süredir açık kalmıştır. (Kontrol teorisinde, bu, anahtarlamalı doğrusal sistemlerin keyfi anahtarlama altında kararlılığının karar verilebilirliğine eşdeğerdir.)

Sorunuzun aşağıdaki varyantının kararsız olduğu biliniyor: Sonlu bir kare matris seti verildiğinde, tüm ürünlerin sınırlı kalıp kalmayacağına karar verin; buraya bakın .

Yukarıdakilerin kararsızlığı, 47x47 boyutunda sadece 2 matrisiniz olsa bile geçerli kalır: buraya bakın .

JSR dilinde, "JSR nedir?" karar verilemez (yukarıdaki referanslara bakın), ancak "JSR < 1 ?" açık. Son soru "rasyonel sonluluk varsayımı" ile ilgilidir: Rasyonel sonluluk varsayımı doğruysa, sorduğunuz soru karar verilebilir.1<1

Son olarak, P = NP olmadıkça, JSR polinom zamanında ( bu makalede tanımlanan tam anlamıyla) yakın olamaz .

Sonuç olarak, etkili bir algoritma iddia eden yukarıdaki cevaplardan biri yanlış olmalıdır.

Olumlu tarafta, JSR'ye yaklaşmak için çeşitli algoritmalar (örn. Yarı-yan programlamaya dayalı) vardır. Farklı algoritmalar farklı performans garantilerine sahiptir. Örneğin aşağıdaki Bkz (- ama aynı zamanda referanslara bakın utanmadan kendim ve meslektaşları tarafından oradaki ).

Bazı özel durumlarda, sorduğunuz soru polinom zamanına karar verilebilir. Örneğin, matrisler simetrik olduğunda veya birinci sırada olduğunda veya işe gidip gelmediklerinde.

Son olarak, konuyla ilgili harika bir kitap şu .


Lütfen sorduğum sorunun resmi ifadesini okuyun - JSR'nin kesinlikle birinden daha az olup olmadığına karar vermekle eşdeğer değildir . Belki de sorunun başlığıyla yanıltılmışsınızdır. Kısacası, asimtotik anlamda değil , sonlu zamanda sıfıra eşit olan her ürünü soruyorum .
Robinson

2
O zaman sorduğunuz soru çok daha basit. Aşağıdakiler eşdeğerdir: (i) Tanımladığınız koşul (ii) Tüm sonlu ürünler sıfır noktadır (iii) JSR = 0 (iv) n uzunluğundaki tüm ürünler sıfırdır (n boyuttur, bu sayıdan bağımsızdır matrisler k). Son durum açıkça karar verilebilirliği ifade eder ve eğer gerçekte durumu polinom zamanında kontrol edebilirsiniz. Yazımın sonunda bağlantılı olan Jungers'ın kitabının 2.3.1 bölümüne bakın. Asimptotik versiyonu kastettiğini düşündüğüm için özür dilerim. ("Tüm ürünler sonunda sıfıra eşit" ifadesiyle yanıltılmıştım.)
Amir Ali Ahmadi

Hangi durumda, Joshua Grochow'un cevabı @AmirAliAhmadi kapsamıyor mu?
Suresh Venkat

2
Bana öyle geliyor ki, aklımdakinden farklı bir algoritma ile. (Yine, sorunun "tüm ürünler sıfıra yakınlaşıyor mu" (yani JSR <1?) Olduğunu düşündüğüm için özür dilerim. Karar verilebilirliği açık olan birkaç farklılık var. (1) Önceki yorumumda (i) - (iv) 'nin eşdeğerinde, Konig'in Lemma'sının kullanılması gerektiğini düşünmüyorum. (2) Matrislerin doğrusal kombinasyonlarını neden aldığını anlamıyorum. (3) Jungers'ın kitabının 2.3.1 Bölümünden Leonid Gurvits'e atfedilen basit bir alternatif algoritmanın altına kopyalarım.
Amir Ali Ahmadi

4
Elimizde çek gerek [yukarıda ... den devam] uzunluğunun tüm ürünler olup olmadığıdır sıfır, ama orada k n böyle matrisler. Bundan kaçınmak için, aşağıdaki matrisleri yinelemeli olarak tanımlayın: X 0 = I , X j = k i = 1 A T i X j - 1 A i . Daha sonra, X n =  n A T A uzunluğunda bir ürün vardır . Bu matris k n ile hesaplanabilirnknX0=I, Xj=i=1kAiTXj1AiXn=A product of length nATAknmatris çarpımlarıdır ve yalnızca uzunluğundaki tüm ürünler sıfırsa sıfırdır. n
Amir Ali Ahmadi

0

Düzenleme: Bu cevap ne yazık ki yanlış. Hata aşağıda vurgulanmıştır. Matrisleri transpoze etmemize izin verilirse argüman işe yarar.

Bir lemmayı kanıtlayarak başlıyoruz.

Lemma. Let A 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 t0 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 that n=3, and write

A=(abcdefghi),N=(010001000).
We start by calculating AN2:
AN2=(00a00d00g).
This matrix is in triangular form, and so if AN2 is nilpotent then g=0. Continue with AN1:
AN1=(0ab0de0gh)=(0ab0de00h).
Again the matrix is in triangular form, and so if AN1 is nilpotent then d=h=0. Continuing,
AN0=(abc0ef00i).
As before, we conclude that a=e=i=0, and so A is upper triangular with zeroes on the diagonal.

If we now consider N2A,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 matrices A1,,Ak satisfy property P if for every infinite sequence i1,[k], we have Ai1Aim=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 A1A2A2A1. Change basis so that A1 is in Jordan normal form, and let the corresponding decomposition of the vector space be V1Vt. Let Vi be a vector space on which A1A2A2A1; note that dimVi>1 since 0 commutes with everything. Restricted to Vi, A1=N and A20. Therefore the lemma implies that for some t0, either A2A1t or A1tA2 is not nilpotent, and therefore property P clearly fails.

Summarizing, property P holds iff all matrices are nilpotent and all of them commute.


4
The last sentence of your lemma's proof is not correct. N2A nilpotent implies g=0, N1A nilpotent gives d=h=0, and N0A nilpotent gives a=e=i=0. So we only conclude that A is upper triangular with zeros on the diagonal, not that A is diagonal (and hence zero).
Joshua Grochow

Indeed , this answer is not correct. If no one else does, I'll post a counter example to both the lemma and the final assertion when I get home later today.
robinson

5
As usual, it is when something is claimed but not proved that the proof fails. Oh well...
Yuval Filmus

1
So the example I had in mind was:
A0=(010001000),A1=(011000000)
One can verify that every product of sufficient length of these two matrices is zero but they don't commute, and the second one is not zero.
robinson
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.