Choleskey Decomposition

Choleskey Decomposition is a decomposition of a Hermitian, positive-definite matrix into the product of a lower triangular matrix and its conjugate transpose, which is useful for efficient numerical solutions

A matrix Input

B matrix Input

System of linear equations

[61555155522555225979]\begin{bmatrix}6 & 15 & 55 \\ 15 & 55 & 225 \\ 55 & 225 & 979\end{bmatrix}
\cdot
[x1x2x3]\begin{bmatrix} x_{1} \\ x_{2} \\ x_{3} \end{bmatrix}
==
[762951259]\begin{bmatrix}76 \\ 295 \\ 1259\end{bmatrix}

Idea is to solve 2 less complex systems

Ax=b(UTU)x=bUTy=bUx=y\begin{align*} Ax &= b \\ (U^T U) x &= b \\ U^T y &= b \\ Ux &= y \end{align*}

Finding U matrix

U

==
[0.000.000.000.000.000.000.000.000.00]\begin{bmatrix}0.00 & 0.00 & 0.00 \\ 0.00 & 0.00 & 0.00 \\ 0.00 & 0.00 & 0.00\end{bmatrix}

Solution vector

[x1x2x3]\begin{bmatrix} x_{1} \\ x_{2} \\ x_{3} \end{bmatrix}
==
[000]\begin{bmatrix}0 \\ 0 \\ 0\end{bmatrix}

Checking if

AxB0\begin{align*} Ax - B &\approx 0 \\ \end{align*}
[000]\begin{bmatrix}0 \\ 0 \\ 0\end{bmatrix}