贾哥高等数值分析第一次实验.pdfVIP

学习资料收集于网络，仅供学习和参考 ,如有侵权，请联系网站删除 高等数值分析第一次实验 第一题： 构造例子说明 CG 的数值形态。当步数 = 阶数时 CG 的解如何？当 A 的最大特征值远大于 第二个最大特征值，最小特征值远小于第二个最小特征值时，方法的收敛性如何？ 解： 用 Housholder 变换和对角阵构造 1000 阶正定对称矩阵 A ： 1）构造对角阵 D = diag( linspace(1, 1000, 1000) ) ； 2 ）构造 Householder 阵 H 。取单位向量 w=[1,0,0,0] T ，I 为 1000 阶单位矩阵， H = I – wTw 。 T 3）构造对称正定矩阵 A 。A = H DH 。由于 D 是对角阵， H 是对称的，所以 A 对称；且 A 与 D 具有相同的特征值 linspace(1, 1000, 1000) 0 ，因此 A 对阵正定。 4 ）b=rand(1000,1); 取初始解 x0=zer