In mathematics, particularly matrix theory and combinatorics, the Pascal matrix is an infinite matrixcontaining the binomial coefficients as its elements. There are three ways to achieve this: as either an upper-triangular matrix, a lower-triangular matrix, or a symmetric matrix. The 5×5 truncations of these are shown below. Lower triangular: Symmetric: Upper triangular: These matrices have the pleasing relationship Sn = LnUn. From this it is easily seen that all three matrices have determinant 1, as the determinant of a triangular matrix is simply the product of its diagonal elements, which are all 1 for both Ln and Un. In other words, matrices Sn, Ln, and Un are unimodular, with Ln and Un having tracen. The elements of the symmetric Pascal matrix are the binomial coefficients, i.e. In other words, Thus the trace of Sn is given by with the first few terms given by the sequence 1, 3, 9, 29, 99, 351, 1275, ….
Construction
The Pascal matrix can actually be constructed by taking the matrix exponential of a special subdiagonal or superdiagonal matrix. The example below constructs a 7-by-7 Pascal matrix, but the method works for any desired n×n Pascal matrices. It is important to note that one cannot simply assume expexp = exp, for A and Bn×n matrices. Such an identity only holds when AB = BA. In the construction of symmetric Pascal matrices like that above, the sub- and superdiagonal matrices do not commute, so the tempting simplification involving the addition of the matrices cannot be made. A useful property of the sub- and superdiagonal matrices used in the construction is that both are nilpotent; that is, when raised to a sufficiently high integer power, they degenerate into the zero matrix. As the n×n generalised shift matrices we are using become zero when raised to power n, when calculating the matrix exponential we need only consider the first n + 1 terms of the infinite series to obtain an exact result.
Variants
Interesting variants can be obtained by obvious modification of the matrix-logarithm PL7 and then application of the matrix exponential. The first example below uses the squares of the values of the log-matrix and constructs a 7-by-7 "Laguerre"- matrix The second example below uses the products v of the values of the log-matrix and constructs a 7-by-7 "Lah"- matrix Using v instead provides a diagonal shifting to bottom-right. The third example below uses the square of the original PL7-matrix, divided by 2, in other words: the first-order binomials in the second subdiagonal and constructs a matrix, which occurs in context of the derivatives and integrals of the Gaussian error function: If this matrix is inverted, then this matrix has alternating signs and gives the coefficients of the derivatives the integrals of the Gauss' error-function. Another variant can be obtained by extending the original matrix to negative values: