Matrices infinies avec le triangle de Pascal comme éléments
En mathématiques , en particulier la théorie de la matrice et combinatoire , une matrice Pascal est un (éventuellement infinie ) la matrice contenant les coefficients binomiaux que ses éléments. Il s'agit donc d'un codage du triangle de Pacal sous forme matricielle. Il existe trois manières naturelles d'y parvenir : en tant que matrice triangulaire inférieure , matrice triangulaire supérieure ou matrice symétrique . Par exemple, les matrices 5 × 5 sont :
Il existe d'autres manières de mettre le triangle de Pascal sous forme matricielle, mais celles-ci ne sont pas facilement étendues à l'infini.
Définition
Les éléments non nuls d'une matrice de Pascal sont donnés par les coefficients binomiaux :
où les indices i , j commencent à 0, et ! désigne la factorielle .
Propriétés
Les matrices ont la relation agréable S n = L n U n . De là, on voit facilement que les trois matrices ont le déterminant 1, car le déterminant d'une matrice triangulaire est simplement le produit de ses éléments diagonaux, qui sont tous 1 pour L n et U n . En d'autres termes, les matrices S n , L n et U n sont unimodulaires , L n et U n ayant la trace n .
La trace de S n est donnée par
![{\text{tr}}(S_{n})=\sum _{i=1}^{n}{\frac {[2(i-1)]!}{[(i-1)!]^ {2}}}=\somme _{k=0}^{n-1}{\frac {(2k)!}{(k!)^{2}}}](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/807ef023bf273f4d6285d938df8948ab45df0332)
avec les premiers termes donnés par la séquence 1, 3, 9, 29, 99, 351, 1275, … (séquence A006134 dans l' OEIS ).
Construction
La matrice de Pascal peut en fait être construite en prenant l' exponentielle matricielle d'une matrice spéciale sous- diagonale ou superdiagonale . L'exemple ci-dessous construit une matrice Pascal 7 × 7, mais la méthode fonctionne pour toutes les matrices Pascal n × n souhaitées . Les points dans les matrices suivantes représentent zéro élément.
![{\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&{&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];\\\\\donc &S_{7}=\exp \left(\left[{\begin{smallmatrix}.&.&.&.&.&.&.\ \1&.&.&.&.&.&.\\.&2&.&.&.&.&.\\.&.&3&.&.&.&.\\.&.&.&4&.&. &.\\.&.&.&.&5&.&.\\.&.&.&.&.&6&.\end{smallmatrix}}\right]\right)\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&7\\1&3&6&10&15&21&28\\1&4&10&20&35&56&84\\1&5&15&35&70&126&210\&&&1&626&21&15&21&28\\1&4&10&20&35&56&84\\1&5&15&35&70&126&210\&&&&126&21}arra\&1&126&21}](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/82e99eb0250c3ea2302d5649dfad45d19d1d034b)
Il est important de noter qu'on ne peut pas simplement supposer exp( A ) exp( B ) = exp( A + B ), pour n × n matrices A et B ; cette égalité n'est vraie que lorsque AB = BA (c'est-à-dire lorsque les matrices A et B commutent ). Dans la construction de matrices Pascal symétriques comme celle ci-dessus, les matrices sous-diagonales et superdiagonales ne commutent pas, de sorte que la simplification (peut-être) tentante impliquant l'ajout des matrices ne peut pas être effectuée.
Une propriété utile des matrices sous-diagonales et superdiagonales utilisées dans la construction est qu'elles sont toutes deux nilpotentes ; c'est-à-dire que lorsqu'ils sont élevés à une puissance entière suffisamment élevée , ils dégénèrent en la matrice zéro . (Voir la matrice de décalage pour plus de détails.) Comme les matrices de décalage généralisées n × n que nous utilisons deviennent nulles lorsqu'elles sont élevées à la puissance n , lors du calcul de la matrice exponentielle, nous n'avons qu'à considérer les n + 1 premiers termes de la série infinie pour obtenir un résultat exact.
Variantes
Des variantes intéressantes peuvent être obtenues par modification évidente de la matrice-logarithme PL 7 puis application de l'exponentielle matricielle.
Le premier exemple ci-dessous utilise les carrés des valeurs de la log-matrice et construit une matrice 7 × 7 "Laguerre" (ou matrice de coefficients de polynômes de 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 petite matrice}}\right];\quad \end{array}}](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/e04ba5c1436107e3df27418ec52690307031066f)
La matrice de Laguerre est en fait utilisée avec une autre mise à l'échelle et/ou le schéma de signes alternés. (La littérature sur les généralisations à des puissances supérieures n'est pas encore trouvée)
Le deuxième exemple ci-dessous utilise les produits v ( v + 1) des valeurs de la log-matrice et construit une matrice 7 × 7 "Lah" (ou matrice de coefficients de nombres de 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)
L'utilisation de v ( v − 1) fournit à la place un décalage diagonal vers le bas à droite.
Le troisième exemple ci-dessous utilise le carré de la matrice PL 7 originale , divisé par 2, en d'autres termes : les binômes du premier ordre (binôme( k , 2)) dans la deuxième sous-diagonale et construit une matrice, qui se produit dans le contexte de les dérivées et intégrales de la fonction d'erreur gaussienne :
![{\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)
Si cette matrice est inversée (en utilisant, par exemple, la matrice-logarithme négative), alors cette matrice a des signes alternés et donne les coefficients des dérivées (et par extension les intégrales) de la fonction d'erreur de Gauss. (La littérature sur les généralisations à des puissances supérieures n'a pas encore été trouvée.)
Une autre variante peut être obtenue en étendant la matrice d'origine aux valeurs négatives :
![{\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&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)
Voir également
Les références
- GS Call et DJ Velleman, "Pascal's matrices", American Mathematical Monthly , volume 100, (avril 1993) pages 372-376
-
Edelman, Alain ; Strang, Gilbert (mars 2004), "Pascal Matrices" (PDF) , American Mathematical Monthly , 111 (3) : 361–385, doi : 10.2307/4145127 , archivé à partir de l'original (PDF) le 04/07/2010
Liens externes