Matrici infinite cu triunghiul lui Pascal ca elemente
În matematică , în special teoria matrice și combinatorică , o matrice Pascal este un (posibil infinit ) matrice care conține coeficienții binomiali ca elementele sale. Este astfel o codificare a triunghiului lui Pacal sub formă de matrice. Există trei moduri naturale de a realiza acest lucru: ca o matrice triunghiulară inferioară , o matrice triunghiulară superioară sau o matrice simetrică . De exemplu, matricile 5 × 5 sunt:
Există și alte moduri în care triunghiul lui Pascal poate fi pus sub formă de matrice, dar acestea nu sunt ușor extinse până la infinit.
Definiție
Elementele diferite de zero ale unei matrice Pascal sunt date de coeficienții binomiali :
unde indicii i , j încep de la 0 și! denotă factorialul .
Proprietăți
Matricile au relația plăcută S n = L n U n . Din aceasta se vede ușor că toate cele trei matrice au determinantul 1, deoarece determinantul unei matrice triunghiulare este pur și simplu produsul elementelor sale diagonale, care sunt toate 1 atât pentru L n cât și pentru U n . Cu alte cuvinte, matricile S n , L n și U n sunt unimodulare , cu L n și U n având urme n .
Urma lui S n este dată de
![{\ 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)
cu primii termeni dați de secvența 1, 3, 9, 29, 99, 351, 1275,… (secvența A006134 din OEIS ).
Constructie
Matricea Pascal poate fi de fapt construită luând matricea exponențială a unei matrice speciale subdiagonale sau superdiagonale . Exemplul de mai jos construiește o matrice Pascal 7 × 7, dar metoda funcționează pentru orice matrice n × n Pascal dorită . Punctele din următoarele matrice reprezintă zero elemente.
![{\ begin {array} {lll} & L_ {7} = \ exp \ left (\ left [{\ begin {smallmatrix}. &. &. &. &. &. &. \\ 1 &. &. &. & . &. &. \\. & 2 &. &. &. &. &. \\. &. & 3 &. &. &. &. \\. &. &. & 4 &. &. &. \\. &. & . &. & 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 \\. &. &. &. &. &. &. \ End {smallmatrix}} \ right] \ right) = \ left [{\ begin {smallmatrix} 1 & 1 & 1 & 1 & 1 & 1 & 1 \\. & 1 & 2 & 3 & 4 & 5 & 6 \\. &. & 1 & 3 & 6 & 10 & 15 \\. &. &. & 1 & 4 & 10 & 20 \\. &. &. &. & 1 & 5 & 15 \\. &. &. &. & 1 & 6 \\. &. &. &. &. &. & 1 \ end {smallmatrix}} \ right]; \\\\\ așadar & S_ {7} = \ exp \ left (\ left [{\ begin {smallmatrix}. &. &. &. &. &. &. \ \ 1 &. &. &. &. &. &. \\. & 2 &. &. &. &. &. \\. &. & 3 &. &. &. &. \\. &. &. & 4 &. &. &. \\. &. &. &. & 5 &. &. \\. &. &. &. &. & 6 &. \ End {smallmatrix}} \ right] \ right) \ exp \ left (\ left [{\ începe {smallmatrix}. & 1 &. &. &. &. &. \\. &. & 2 &. &. &. &. \\. &. &. & 3 &. &. &. \\. &. &. &. & 4 &. &. \\ . &. &. &. &. & 5 &. \\. &. &. &. &. &. & 6 \\. &. &. &. &. &. \. {End {smallmatrix}} \ right] \ dreapta) = \ left [{\ begin {smallmatrix} 1 & 1 & 1 & 1 & 1 & 1 & 1 \\ 1 & 2 & 3 & 4 & 5 & 6 & 7 \\ 1 & 3 & 6 & 10 & 15 & 21 & 28 \\ 1 & 4 & 10 & 20 & 35 & 56 & 84 \\ 1 & 5 & 15 & 35 & 70 & 126 & 210 & 1 & 6 & 21 & 56 & 46 & 56 & 46 & 56](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/82e99eb0250c3ea2302d5649dfad45d19d1d034b)
Este important de reținut că nu se poate presupune pur și simplu exp ( A ) exp ( B ) = exp ( A + B ), pentru n × n matricile A și B ; această egalitate se menține numai atunci când AB = BA (adică atunci când matricile A și B fac naveta ). În construcția matricilor Pascal simetrice ca cea de mai sus, matricile sub și superdiagonale nu fac naveta, astfel încât simplificarea (poate) tentantă care implică adăugarea matricilor nu poate fi realizată.
O proprietate utilă a matricilor sub- și superdiagonale utilizate în construcție este că ambele sunt nilpotente ; adică, când sunt ridicați la o putere întreagă suficient de mare , degenerează în matricea zero . (A se vedea matricea de schimbare pentru detalii suplimentare.) Deoarece n × n matricele de schimbare generalizate pe care le folosim devin zero atunci când sunt ridicate la puterea n , atunci când calculăm exponențialul matricei trebuie să luăm în considerare doar primii termeni n + 1 din seria infinită pentru a obține un rezultatul exact.
Variante
Variante interesante pot fi obținute prin modificarea evidentă a matricei-logaritmului PL 7 și apoi prin aplicarea exponențialei matricei.
Primul exemplu de mai jos folosește pătratele valorilor matricei log și construiește o matrice "Laguerre" de 7 × 7 (sau matricea coeficienților polinomilor Laguerre
![{\ 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)
Matricea Laguerre este de fapt folosită cu alte scalări și / sau schema de semne alternative. (Literatura despre generalizări la puteri superioare nu este încă găsită)
Al doilea exemplu de mai jos folosește produsele v ( v + 1) ale valorilor matricei log și construiește o matrice de 7 × 7 "Lah" - (sau matricea coeficienților numerelor 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)
Utilizarea v ( v - 1) oferă în schimb o deplasare în diagonală în partea dreaptă jos.
Al treilea exemplu de mai jos folosește pătratul matricei PL 7 originale , împărțit la 2, cu alte cuvinte: binomii de ordinul întâi (binom ( k , 2)) în al doilea subdiagonal și construiește o matrice, care apare în contextul a derivatelor și integralelor ale Gaussian funcției de eroare :
![{\ 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)
Dacă această matrice este inversată (folosind, de exemplu, matricea-logaritm negativ), atunci această matrice are semne alternante și dă coeficienții derivaților (și prin extensie integralele) funcției de eroare a lui Gauss. (Literatura despre generalizări la puteri superioare nu este încă găsită.)
O altă variantă poate fi obținută prin extinderea matricei originale la valori negative :
![{\ displaystyle {\ begin {array} {lll} & \ exp \ left (\ left [{\ begin {smallmatrix}. &. &. &. &. &. &. &. &. &. &. \\ - 5 &. &. &. &. &. &. &. &. &. &. &. \\. & - 4 &. &. &. &. &. &. &. &. &. \\. &. & - 3 &. &. &. &. &. &. &. &. &. \\. &. &. & - 2 &. &. &. &. &. &. &. &. \\. &. &. &. & - 1 &. &. &. &. &. &. &. \\. &. &. &. &. & 0 &. &. &. &. &. &. \\ . &. &. &. &. &. & 1 &. &. &. &. &. \\. &. &. &. &. . &. &. &. &. &. & 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)
Vezi si
Referințe
- GS Call și DJ Velleman, „Matricile lui Pascal”, American Mathematical Monthly , volumul 100, (aprilie 1993) paginile 372–376
-
Edelman, Alan; Strang, Gilbert (martie 2004), „Pascal Matrices” (PDF) , American Mathematical Monthly , 111 (3): 361–385, doi : 10.2307 / 4145127 , arhivat din original (PDF) în data de 04.07.2010
linkuri externe