矩阵计算上机实验告(哈工程).docVIP

2012-2013学年第二学期《矩阵计算》上机实验报告 哈尔滨工程大学 矩阵计算上机实验报告 姓名： 学号： 矩阵计算实验报告 实验名称 矩阵计算上机实验 指导教师 姓名 年级 学号 成绩 实验一、用共轭梯度法求实对称线性方程组的解 一、实验目的： 　　１、学习使用matlab编写数值计算程序。 　　２、了解共轭梯度法原理和编程方法。 二、实验要求： 　　１、阐述共轭梯度法原理和编程方法。 ２、根据共轭梯度法算法原理编写matlab程序，并运行的出结果。 （提示：下文是Jacobi迭代求解线性方程组的实验过程，仅供参考。提交作业时，请删除！） 三、实验原理： 设n阶非奇异矩阵Ａ的主对角元全不为０，记Ｄ＝diag（ ）是非奇异对角阵，做Ａ的一个分裂： 。 记。迭代过程式为。 这种迭代方法称为Jacobi迭代法。其迭代矩阵记为：右端向量记为： = - 记，Jacobi迭代过程为： 对照等号两边，得到计算(i=1,2,…,n)的公式为 ，或，i=1,2,… 四、实验内容与步骤： 　　１、实验内容：依照实验原理编写共轭梯度法的程序。 　　２、实验步骤：首先，在电脑上安装matlab，然后，启动matlab，新建一个M文件。 实验程序如下： 程序代码 程序代码说明 function [x,k]=jacobi(,b); x=