我试图分解以下拉普拉斯矩阵$D^TD$, D 是一阶导数矩阵。使用 Neumann 边界条件的二阶导数矩阵的三对角形式由下式给出，

$\begin{bmatrix} -1 & 1 & 0 & 0 & 0 \\ 1 & -2 & 1 & 0 &0\\ 0 & 1 & -2 &1&0\\0&0&1&-2&1\\0&0&0&1&-1\end{bmatrix}$.

为了用梯度（一阶导数）算子来写上面的内容，我尝试使用 Cholesky 分解来分解拉普拉斯算子。但是，我无法成功分解。

关于如何找到表格分解的任何建议$D^TD$?