2021-10-08

共轭梯度算法matlab程序实现

首先要把矩阵弄成方阵，然后再进行迭代。
程序如下：

初始化：S_CG=S’S;
(处理为方阵）
r=b_CG-S_CGsigma_CG;
p=r;

之后迭代：
sigma_CG=sigma_CG+alphap;
r=r-alphaS_CGp;
p=r-betap;

注意迭代时数据改变的先后顺序。