双共轭梯度算法

双共轭梯度算法的基本流程如下所示，具体如下所示： 1）、初步给定 x 0 2）、计算 r0=b-Ax 0 ； 令 r0=r0 ； 3）、循环体计算： for    i=1,2,… if   i=1 pi =ri -1  ;  pi =ri -1   ; else pi =ri -1+βi-1 pi-1  ;  pi =ri -1+βi-1 pi-1   ; end αi =ri -1Tri -1p

>>阅读原文<<