Unendliche Matrizen mit dem Pascalschen Dreieck als Elementen
In der Mathematik , insbesondere der Matrixtheorie und Kombinatorik , ist eine Pascal-Matrix eine (möglicherweise unendliche ) Matrix, die die Binomialkoeffizienten als ihre Elemente enthält. Es ist somit eine Kodierung des Pacalschen Dreiecks in Matrixform. Es gibt drei natürliche Wege, dies zu erreichen: als untere Dreiecksmatrix , als obere Dreiecksmatrix oder als symmetrische Matrix . Zum Beispiel sind die 5 × 5-Matrizen:
Es gibt noch andere Möglichkeiten, das Pascalsche Dreieck in Matrixform zu bringen, aber diese lassen sich nicht ohne weiteres ins Unendliche erweitern.
Definition
Die von Null verschiedenen Elemente einer Pascal-Matrix sind durch die Binomialkoeffizienten gegeben :
wobei die Indizes i , j bei 0 beginnen und ! bezeichnet die Fakultät .
Eigenschaften
Die Matrizen haben die angenehme Beziehung S n = L n U n . Daraus ist leicht ersichtlich, dass alle drei Matrizen die Determinante 1 haben, da die Determinante einer Dreiecksmatrix einfach das Produkt ihrer diagonalen Elemente ist, die sowohl für L n als auch für U n alle 1 sind . Mit anderen Worten, die Matrizen S n , L n und U n sind unimodular , wobei L n und U n die Spur n haben .
Die Spur von S n ist gegeben durch
![{\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)
mit den ersten Termen der Sequenz 1, 3, 9, 29, 99, 351, 1275, … (Sequenz A006134 im OEIS ).
Konstruktion
Die Pascal - Matrix kann tatsächlich , indem man die konstruiert werden Exponentialmatrix einer speziellen subdiagonal oder superdiagonal Matrix. Das folgende Beispiel konstruiert eine 7 × 7 Pascal-Matrix, aber das Verfahren funktioniert für beliebige n × n Pascal-Matrizen. Die Punkte in den folgenden Matrizen repräsentieren Nullelemente.
![{\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];\\\\\daher &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]\ rechts)=\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&9&9&126&210&1&6&21&56&126&25&sub2;](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/82e99eb0250c3ea2302d5649dfad45d19d1d034b)
Es ist wichtig zu beachten, dass man nicht einfach exp( A ) exp( B ) = exp( A + B ) für n × n Matrizen A und B annehmen kann ; diese Gleichheit gilt nur, wenn AB = BA (dh wenn die Matrizen A und B kommutieren ). Bei der Konstruktion symmetrischer Pascal-Matrizen wie der obigen kommutieren die sub- und superdiagonalen Matrizen nicht, so dass die (vielleicht) verlockende Vereinfachung mit der Addition der Matrizen nicht gemacht werden kann.
Eine nützliche Eigenschaft der bei der Konstruktion verwendeten sub- und superdiagonalen Matrizen besteht darin, dass beide nilpotent sind ; das heißt, wenn sie auf eine ausreichend hohe ganzzahlige Potenz angehoben werden , degenerieren sie in die Nullmatrix . (Siehe Schiebematrix für weitere Details.) Da die n × n verallgemeinerten Schiebematrizen, die wir verwenden, null werden, wenn sie mit n potenziert werden, brauchen wir bei der Berechnung der Matrixexponentialfunktion nur die ersten n + 1 Terme der unendlichen Reihe zu berücksichtigen , um an . zu erhalten genaues Ergebnis.
Varianten
Interessante Varianten erhält man durch offensichtliche Modifikation des Matrix-Logarithmus PL 7 und anschließende Anwendung der Matrix-Exponentialfunktion.
Das erste Beispiel unten verwendet die Quadrate der Werte der Log-Matrix und konstruiert eine 7 × 7 "Laguerre"-Matrix (oder Koeffizientenmatrix von Laguerre-Polynomen
![{\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\end{ smallmatrix}}\right];\quad \end{array}}](/criselda-https-wikimedia.org/api/rest_v1/media/math/render/svg/e04ba5c1436107e3df27418ec52690307031066f)
Die Laguerre-Matrix wird tatsächlich mit einer anderen Skalierung und/oder dem Schema alternierender Vorzeichen verwendet. (Literatur über Verallgemeinerungen auf höhere Potenzen ist noch nicht vorhanden)
Das zweite Beispiel unten verwendet die Produkte v ( v + 1) der Werte der Log-Matrix und konstruiert eine 7 × 7 "Lah"-Matrix (oder Matrix von Koeffizienten von Lah-Zahlen )
![{\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)
Die Verwendung von v ( v − 1) bietet stattdessen eine diagonale Verschiebung nach unten rechts.
Das dritte Beispiel unten verwendet das Quadrat der ursprünglichen PL 7 -Matrix, dividiert durch 2, mit anderen Worten: die Binomiale erster Ordnung (binomial( k , 2)) in der zweiten Subdiagonalen und konstruiert eine Matrix, die im Kontext von die Ableitungen und Integrale der Gaußschen Fehlerfunktion :
![{\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)
Wird diese Matrix invertiert (zB mit dem negativen Matrix-Logarithmus), dann hat diese Matrix alternierende Vorzeichen und gibt die Koeffizienten der Ableitungen (und damit die Integrale) der Gaußschen Fehlerfunktion an. (Literatur über Verallgemeinerungen auf höhere Potenzen ist noch nicht vorhanden.)
Eine andere Variante kann erhalten werden, indem die ursprüngliche Matrix auf negative Werte erweitert wird :
![{\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)
Siehe auch
Verweise
- GS Call und DJ Velleman, "Pascal's matrices", American Mathematical Monthly , Band 100, (April 1993) Seiten 372–376
-
Edelmann, Alan; Strang, Gilbert (März 2004), "Pascal Matrices" (PDF) , American Mathematical Monthly , 111 (3): 361–385, doi : 10.2307/4145127 , archiviert vom Original (PDF) am 2010-07-04
Externe Links