Wilson matrix

1. ${\displaystyle W}$ is a symmetric matrix.
2. ${\displaystyle W}$ is positive definite.
3. The determinant of ${\displaystyle W}$ is ${\displaystyle 1}$.
4. The inverse of ${\displaystyle W}$ is ${\displaystyle W^{-1}={\begin{bmatrix}68&-41&-17&10\\-41&25&10&-6\\-17&10&5&-3\\10&-6&-3&2\end{bmatrix}}}$
5. The characteristic polynomial of ${\displaystyle W}$ is ${\displaystyle \lambda ^{4}-35\lambda ^{3}+146\lambda ^{2}-100\lambda +1}$.
6. The eigenvalues of ${\displaystyle W}$ are ${\displaystyle \quad 0.01015004839789187,\quad 0.8431071498550294,\quad 3.858057455944953,\quad 30.28868534580213}$.
7. Since ${\displaystyle W}$ is symmetric, the 2-norm condition number of ${\displaystyle W}$ is ${\displaystyle \kappa _{2}(W)=({\text{max eigen value}})/({\text{min eigen value}})=30.28868534580213/0.01015004839789187=2984.09270167549}$.
8. The solution of the system of equations ${\displaystyle (S1)}$ is ${\displaystyle x=y=z=u=1}$.
9. The Cholesky factorisation of ${\displaystyle W}$ is ${\displaystyle W=R^{T}R}$ where ${\displaystyle R={\begin{bmatrix}{\sqrt {5}}&{\frac {7}{\sqrt {5}}}&{\frac {6}{\sqrt {5}}}&{\sqrt {5}}\\0&{\frac {1}{\sqrt {5}}}&-{\frac {2}{\sqrt {5}}}&0\\0&0&{\sqrt {2}}&{\frac {3}{\sqrt {2}}}\\0&0&0&{\frac {1}{\sqrt {2}}}\end{bmatrix}}}$.
10. ${\displaystyle W}$ has the factorisation ${\displaystyle W=LDL^{T}}$ where ${\displaystyle L={\begin{bmatrix}1&0&0&0\\{\frac {7}{5}}&1&0&0\\{\frac {6}{5}}&-2&1&0\\1&0&{\frac {3}{2}}&1\end{bmatrix}},\quad D={\begin{bmatrix}5&0&0&0\\0&{\frac {1}{5}}&0&0\\0&0&2&0\\0&0&0&{\frac {1}{2}}\end{bmatrix}}}$.
11. ${\displaystyle W}$ has the factorisation ${\displaystyle W=Z^{T}Z}$ with ${\displaystyle Z}$ as the integer matrix[7] ${\displaystyle Z={\begin{bmatrix}2&3&2&2\\1&1&2&1\\0&0&1&2\\0&0&1&1\end{bmatrix}}}$.

A consideration of the condition number of the Wilson matrix has spawned several interesting research problems relating to condition numbers of matrices in certain special classes of matrices having some or all the special features of the Wilson matrix. In particular, the following special classes of matrices have been studied:[1]

1. ${\displaystyle S=}$ the set of ${\displaystyle 4\times 4}$ nonsingular, symmetric matrices with integer entries between 1 and 10.
2. ${\displaystyle P=}$ the set of ${\displaystyle 4\times 4}$ positive definite, symmetric matrices with integer entries between 1 and 10.

An exhaustive computation of the condition numbers of the matrices in the above sets has yielded the following results:

1. Among the elements of ${\displaystyle S}$, the maximum condition number is ${\displaystyle 7.6119\times 10^{4}}$ and this maximum is attained by the matrix ${\displaystyle {\begin{bmatrix}2&7&10&10\\7&10&10&9\\10&10&10&1\\10&9&1&10\end{bmatrix}}}$.
2. Among the elements of ${\displaystyle P}$, the maximum condition number is ${\displaystyle 3.5529\times 10^{4}}$ and this maximum is attained by the matrix ${\displaystyle {\begin{bmatrix}9&1&1&5\\1&10&1&9\\1&1&10&1\\5&9&1&10\end{bmatrix}}}$.