Cholesky 分解的课件挺不错,来自中科院研究生院,讲得清楚,还配了详细的 MATLAB 例子。适合你平时对称正定矩阵的时候用,分解速度快,占内存也少,尤其在大规模计算里优势。课件里还顺带提了 QR 分解和 LU 分解,整体是打包讲的,干货挺多。
Cholesky 分解的核心点就是把一个对称正定矩阵拆成一个上三角矩阵和它的转置,方便后面解线性方程组。比如你在解 A * X = b
,可以变成 R' * R * X = b
,算起来轻松多。
顺带说一句,QR 分解也讲到了,任何矩阵都能用,而且适合最小二乘法场景。你只要记住公式 A = QR
,后面解 AX = b
的时候就能变成 RX = Q'b
,逻辑也清楚。
如果你平时用 MATLAB 写科学计算或者搞点建模的活儿,这课件挺值得一看。还贴心附了不少例题和链接,比如:
建议你用 MATLAB 自带的 chol()
函数试试,输入对称正定矩阵,直接就能出上三角结果。用的时候注意矩阵是不是对称正定,否则直接报错。
如果你最近刚好在调线性方程组、优化算法或者数据拟合那块的代码,可以抽空看看这套课件,思路清晰,上手快。