第八章 计算方法——共轭梯度法求解线性方程组.pdf

计算方法上机报告 计算方法上机报告 1 共轭梯度法求解线性方程组 1.1 算法原理及程序框图 当线性方程组Ax = b 的系数矩阵A 是对称正定矩阵是，可以采用共轭梯度法对该 方程组进行求解，可以证明，式(1)所示的n 元二次函数 1 T T f (x ) x Ax b x (1) 2 * 取得极小值点x 是方程Ax = b 的解。共轭梯度法是把求解线性方程组的问题转化为求 解一个与之等价的二次函数极小化的问题。从任意给定的初始点出发，沿一组关于矩 阵A 的共轭方向进行线性搜索，在无舍入误差的假定下，最多迭代n 次（其中n 为矩 阵A 的阶数），就可求得二次函数的极小点，也就求得线性方程组Ax = b 的解。其迭 代格式为公式(2) 。 x(k 1) x(k )  d (k ) (2) k 共轭梯度法中关键的两点是迭代格式(2) 中最佳步长 和搜索方向d (k) 的确定。其 k 中 可以通过一元函数f (x(k)+d(k)) 的极小化来求得，其表达式为公式(3) ；取d (0) k = r(0) = b-Ax (0) ，则d(k+ 1) = r(k+1) + d(k) ，要求d(k+ 1)满足 (d(k+ 1) , Ad(k)) = 0 ，可得 的表达 k k 式(4) 。 k T k r  d   k k T k (3) d   Ad   r(k 1)T Ad (k )   (4) k d (k )T Ad (k ) 经过一系列的证明和简化，最终可得共轭梯度法的计算过程如下，计算程序框图 如图1。 (1) 给定初始计算向量x(0) 即精度0 ； (2) 计算r(0) = b -Ax (0) ，取d(0) = r(0) ； (3) for k =0 to n- 1 do