Nekonečné matice s Pascalovým trojúhelníkem jako prvky
V matematice , zejména maticové teorii a kombinatorice , je Pascalova matice (možná nekonečná ) matice obsahující jako prvky binomické koeficienty . Jedná se tedy o kódování Pacalova trojúhelníku v maticové formě. Existují tři přirozené způsoby, jak toho dosáhnout: jako nižší trojúhelníková matice , horní trojúhelníková matice nebo symetrická matice . Například matice 5 × 5 jsou:
Existují i jiné způsoby, kterými lze Pascalův trojúhelník vložit do maticové formy, ale ty nelze snadno rozšířit do nekonečna.
Definice
Nenulové prvky Pascalovy matice jsou dány binomickými koeficienty :
kde indexy i , j začínají na 0, a! označuje faktoriál .
Vlastnosti
Matice mají potěšující vztah S n = L n U n . Z toho je snadno vidět, že všechny tři matice mají determinant 1, protože determinant trojúhelníkové matice je jednoduše součin jejích diagonálních prvků, které jsou všechny 1 pro L n i U n . Jinými slovy, matice S n , L n a U n jsou unimodulární , přičemž L n a U n mají stopu n .
Stopa S n je dána vztahem
![{\ text {tr}} (S_ {n}) = \ sum _ {i = 1}^{n} {\ frac {[2 (i-1)]!} {[(i-1)!]^ {2}}} = \ sum _ {k = 0}^{n-1} {\ frac {(2k)!} {(K!)^{2}}}](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/807ef023bf273f4d6285d938df8948ab45df0332)
s několika prvními termíny danými sekvencí 1, 3, 9, 29, 99, 351, 1275, ... (sekvence A006134 v OEIS ).
Konstrukce
Pascalovu matici lze ve skutečnosti sestrojit tak, že vezmeme exponenciální matici speciální subdiagonální nebo superdiagonální matice. Níže uvedený příklad konstruuje matici Pascal 7 × 7, ale metoda funguje pro libovolné požadované n × n Pascalových matic. Tečky v následujících maticích představují nulové prvky.
![{\ begin {array} {lll} & L_ {7} = \ exp \ left (\ left [{\ begin {smallmatrix}. &. &. &. &. &. &. &. &. \\ 1 &. &. &. & . &. &. \\. & 2 &. &. &. &. &. \\. &. & 3 &. &. &. &. \.. . &. & 5 &. &. \\. &. &. &. &. &. & 6 &. \ End {smallmatrix}} \ right] \ right) = \ left [{\ begin {smallmatrix} 1 &. &. &. &. &. &. \\ 1 & 1 &. &. &. &. &. \. 1 & 2 & 1 &. &. &. &. \. 1 & 3 & 3 & 1 &. &. &. \. 1 & 4 & 6 & 4 & 1 &. &. \. 1 & 5 & 10 & 10 & 5 & 1 &. \\ 1 & 6 & 15 & 20 & 15 & 6 & 1 \ end { smallmatrix}} \ right]; \ quad \\\\ & U_ {7} = \ exp \ left (\ left [{\ begin {smallmatrix}. & 1 &. &. &. &. &. &. \\. &. & 2 &. &. &. &. \\. &. &. & 3 &. &. &. \\. &. &. &. &. 4 &. &. \\. &. &. &. &. &. & 5 &. \\. &. &. &. &. &. &. 6. \.. &. &. &. &. & 1 & 2 & 3 & 4 & 5 & 6 \\. &. & 1 & 3 & 6 & 10 & 15 \\. &. &. & 1 & 4 & 10 & 20 \\. &. &. &. & 1 & 5 & 15 \\. &. &. &. &. &. & 1 & 6 \\. &. &. &. &. &. & 1 \ end {smallmatrix}} \ right]; \\\\\ proto & S_ {7} = \ exp \ left (\ left [{\ begin {smallmatrix}. &. &. &. &. &. &. &. \ \ 1 &. &. &. &. &. &. \\. & 2 &. &. &. &. &. &. \\. &. & 3 &. &. &. &. \.. &. &. &. \\. &. &. &. & 5 &. &. \\. &. &. &. &. &. & 6 &. \ End {smallmatrix}} \ right] \ right) \ exp \ left (\ left [{\ začít {smallmatrix}. & 1 &. &. &. &. &. \\. &. & 2 &. &. &. &. \\. &. &. &. & 3 &. &. \. & 4 &. &. \\ . &. &. &. &. & 5 &. \\. &. &. &. &. &. &. & 6 \\. &. &. &. &. &. &. \. \ End {smallmatrix}} \ right] \ vpravo) \ \ left [smallmatrix} 1 & 1 & 1 & 1 & 1 & 1 \\ 1 & 2 & 3 & 4 & 5 \ 6](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/82e99eb0250c3ea2302d5649dfad45d19d1d034b)
Je důležité si uvědomit, že nelze jednoduše předpokládat exp ( A ) exp ( B ) = exp ( A + B ), pro n × n matic A a B ; tato rovnost platí pouze tehdy, když AB = BA (tj. když dojíždějí matice A a B ). Při konstrukci symetrických Pascalových matic, jako je výše, sub- a superdiagonální matice nepojíždějí, takže (snad) lákavé zjednodušení zahrnující sčítání matic nelze provést.
Užitečnou vlastností sub- a superdiagonálních matic použitých při konstrukci je, že obě jsou nilpotentní ; to znamená, že když se zvýší na dostatečně vysokou celočíselnou mocnost, degenerují do nulové matice . ( Další podrobnosti viz matice posunu .) Protože n × n generalizovaných matic posunu, které používáme, se stanou nulovými, když se zvýší na výkon n , při výpočtu exponenciální matice musíme vzít v úvahu pouze prvních n + 1 členů nekonečné řady, abychom získali přesný výsledek.
Varianty
Zajímavé varianty lze získat zřejmou modifikací maticového logaritmu PL 7 a poté aplikací maticového exponenciálu.
První příklad níže používá čtverce hodnot log-matice a sestrojuje matici „Laguerre“ 7 × 7 (nebo matici koeficientů Laguerrových polynomů
![{\ begin {array} {lll} & LAG_ {7} = \ exp \ left (\ left [{\ begin {smallmatrix}. &. &. &. &. &. &. &. \. 1 &. &. &. & . &. &. \\. & 4 &. &. &. &. &. \\. &. &. & 9 &. &. &. &. \\. &. &. &. & 16 &. &. \.. . &. & 25 &. &. \\. &. &. &. &. &. & 36 &. \ End {smallmatrix}} \ right] \ right) = \ left [{\ begin {smallmatrix} 1 &. &. &. &. &. &. \\ 1 & 1 &. &. &. &. &. \\ 2 & 4 & 1 &. &. &. &. \\ 6 & 18 & 9 & 1 &. &. &. &. \\ 24 & 96 & 72 & 16 & 1 &. &. \\ 120 & 600 & 600 & 200 & 25 & 1 &. \\ 720 & 4320 & 5400 & 2400 & 450 & 36 & 1 & smallmatrix}} \ right]; \ quad \ end {array}}](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/e04ba5c1436107e3df27418ec52690307031066f)
Laguerreova matice se ve skutečnosti používá s jiným měřítkem a/nebo schématem střídajících se znaků. (Literatura o zobecnění na vyšší síly zatím nebyla nalezena)
Druhý příklad níže používá součiny v ( v + 1) hodnot log-matice a sestrojí matici 7 × 7 „Lah“- (nebo matici koeficientů čísel Lah )
![{\ Displaystyle {\ begin {array} {lll} & LAH_ {7} = \ exp \ left (\ left [{\ begin {smallmatrix}. &. &. &. &. &. &. &. \. 2 \. &. &. &. &. &. &. \\. & 6 &. &. &. &. &. \\. &. & 12 &. &. &. &. \\. &. &. &. & 20 &. &. &. \\. &. &. &. & 30 &. &. \\. &. &. &. &. &. & 42 &. \ End {smallmatrix}} \ right] \ right) = \ left [{\ begin {smallmatrix} 1 &. &. & . &. &. &. &. \. 2 & 1 &. &. &. &. &. &. &. \\ 6 & 6 & 1 &. &. &. &. &. \. 24 & 36 & 12 & 1 &. &. &. &. \\ 120 & 240 & 120 & 20 & 1 &. & . &. \\ 720 & 1800 & 1200 & 300 & 30 & 1 &. &. \\ 5040 & 15120 & 12600 & 4200 & 630 & 42 & 1 &. \\ 40320 & 141120 & 141120 & 58800 & 11760 & 1176 & 56 & 1 \ end {smallmatrix}} \ right]; \ quad \ end {array}}}](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/a5f5fa329b5cef35f151bd8d7a04c7021bbe48ba)
Použití v ( v -1) místo toho poskytuje diagonální posun vpravo dole.
Třetí příklad níže používá druhou mocninu původní matice PL 7 , dělenou 2, jinými slovy: binomiky prvního řádu (binomické ( k , 2)) ve druhém subdiagonálu a konstruuje matici, která se vyskytuje v kontextu že deriváty a integrály Gaussova chybové funkce :
![{\ Displaystyle {\ begin {array} {lll} & GS_ {7} = \ exp \ left (\ left [{\ begin {smallmatrix}. &. &. &. &. &. &. &. \\. &. & . &. &. &. &. \\ 1 &. &. &. &. &. &. \.. & 3 &. &. &. &. &. \\. &. & 6 &. &. &. &. \ \. &. &. & 10 &. &. &. \\. &. &. &. &. & 15 &. &. \ End {smallmatrix}} \ right] \ right) = \ left [{\ begin {smallmatrix} 1 &. & . &. &. &. &. \\. & 1 &. &. &. &. &. \. 1 &. & 1 &. &. &. &. \\. & 3 &. & 1 &. &. &. \\ 3 &. & 6 & . & 1 &. &. \\. & 15 &. & 10 &. & 1 &. \\ 15 &. & 45 &. & 15 &. & 1 \ end {smallmatrix}} \ right]; \ quad \ end {array}}}](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/1258220d051b41ff146243babe85ae98c066b71d)
Pokud je tato matice invertována (například pomocí záporného maticového logaritmu), pak má tato matice střídavá znaménka a udává koeficienty derivací (a potažmo integrálů) Gaussovy chybové funkce. (Literatura o zobecnění na vyšší moc zatím nebyla nalezena.)
Další variantu lze získat rozšířením původní matice na záporné hodnoty :
![{\ Displaystyle {\ begin {array} {lll} & \ exp \ left (\ left [{\ begin {smallmatrix}. &. &. &. &. &. &. &. &. &. &. &. &. &. \\-5 &. &. &. &. &. &. &. &. &. &. &. &. \\. &-4 &. &. &. &. &. &. &. &. &. &. \\. &. &-3 &. &. &. &. &. &. &. &. &. &. \\. &. &. &.-2 &. &. &. &. &. &. &. &. & \\. &. &. &. &-1 &. &. &. &. &. &. &. \. \\. &. &. &. &. & 0 &. &. &. &. &. &. &. \\ . &. &. &. &. &. & 1 &. &. &. &. &. \\. &. &. &. &. &. &. & 2 &. &. &. &. \.. . &. &. &. &. &. & 3 &. &. &. \\. &. &. &. &. &. &. &. &. &. & 4 &. &. \\. &. &. &. & . &. &. &. &. &. &. & 5 &. \ End {smallmatrix}} \ right] \ right) = \ left [{\ begin {smallmatrix} 1 &. &. &. &. &. &. &. &. & . &. &. &. \\-5 & 1 &. &. &. &. &. &. &. &. &. &. &. \\ 10 & -4 & 1 &. &. &. &. &. &. &. &. & &. \\-10 & 6 & -3 & 1 &. &. &. &. &. &. &. &. &. \\ 5 & -4 & 3 & -2 & 1 &. &. &. &. &. &. &. \.-1 & 1 & -1 & 1 &- 1 & 1 &. &. &. &. &. &. \\. &. &. &. &. &. & 0 & 1 &. &. &. &. &. \\. &. &. &. &. &. \\. &. &. &. &. &. &. & 1 & 2 & 1 &. &. &. \\. &. &. &. &. &. & 1 & 3 & 3 & 1 &. &. \\. &. &. &. & . &. & 1 & 4 & 6 & 4 & 1 &. \\. &. &. &. &. &. &. & 1 & 5 & 10 & 10 & 5 & 1 \ end {smallmatrix}} \ right]. \ End {array}}}](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/fcc8c5d05a829836beceb789cb8275b2c86b5c2f)
Viz také
Reference
- GS Call a DJ Velleman, „Pascalovy matice“, American Mathematical Monthly , svazek 100, (duben 1993), strany 372–376
-
Edelman, Alan; Strang, Gilbert (březen 2004), „Pascal Matrices“ (PDF) , American Mathematical Monthly , 111 (3): 361–385, doi : 10,2307/4145127 , archivováno z originálu (PDF) dne 2010-07-04
externí odkazy