数学问题：利用Cholesky分解法解决线性方程问题Ax=b,其中(A,b)分别为:解决代码：​#include"windows.h"#include<stdio.h>#include<string.h>#include<stdlib.h>#include<iostream>#include<iomanip>#include<math.h>usingnamespacest

（搬运外网的代码，非原创；原网址）（其实是专业课作业，但感觉国内博客没有合适的代码实现，所以就搬运到自己博客了）背景-Cholesky分解：若A为n阶实对称正定矩阵，则存在非奇异下三角矩阵L，使得A=LL^T。是特殊的LU分解（下三角上三角分解）。若限定L的对角元素为正，则这种分解

对于正定对称矩阵\(\mathbf{H}\)，可以分解为\(\mathbf{H}=\mathbf{XX}^T\)，其中\(\mathbf{X}\)是下三角矩阵。这个分解方法就是cholesky分解，pytorch对应的函数是torch.linalg.cholesky使用\(\mathbf{X}\)可以求出\(\mathbf{H}^{-1}\)，pytorch对应的函数是torch.cholesky_inverse计

求L的步骤：首先求出l11然后求第一列其他元素 

CeresSolver中的线性求解器类型（linear_solver_type）有多个选项，包括：DENSE_QR：使用稠密QR分解方法求解线性方程组。适用于内存足够的小规模问题，求解速度较快。DENSE_SCHUR：

  伟大的算法：  Newtonandquasi-Newtonmethods  Matrixfactorizations(LU,Cholesky,QR)  Singularvaluedecomposition,QRandQZalgorithms  Mon

（17）对称矩阵的cholesky分解补发笔记