贾哥高等数值第一次实验规范.doc

高等数值分析第一次实验 第一题： 构造例子说明CG的数值形态。当步数 = 阶数时CG的解如何？当A的最大特征值远大于第二个最大特征值，最小特征值远小于第二个最小特征值时，方法的收敛性如何？ 解： 用Housholder变换和对角阵构造1000阶正定对称矩阵A： 构造对角阵D = diag( linspace(1, 1000, 1000) )； 构造Householder阵H。取单位向量w=[1,0,0,.....0]T，I为1000阶单位矩阵，H = I – wTw。 构造对称正定矩阵A。A = HTDH。由于D是对角阵，H是对称的，所以A对称；且A与D具有相同的特征值linspace(1, 1000, 1000) 0，因此A对阵正定。 b=rand(1000,1);取初始解x0=zeros(1000,1); 1.计算Ax = b 利用matlab编程实现CG算法。由于实际计算存在机器误差，因此迭代1000步后的残差不等于0，因此不能用rk=0作为停机准则，否则matlab会无休止地计算下去。本例采用停机准则为：迭代步数=1000步。 当D = diag( linspace(1, 1000, 1000) )时，条件数k=1000； 当D = diag( linspace(1, 100, 1000) )时，条件数k=100； 当D = diag( linspace(1, 10,