实验五_线性代数方程组数值解法.pdfVIP

实验五 线性代数方程组的数值解法 材料系 材93 2009011976 邓陟 一． 实验目的 1. 学会用MATLAB 软件数值求解线性代数方程组，对迭代法的收敛性和解的稳定性作 初步分析。 2. 通过实例学习用线性代数方程组解决简化的实际问题。 二． 实验内容 1. 已知方程组Ax b ；其中A  2020 ，定义为  3 1 2 1 4    1 2 3 1 2 1 4  1 4 1 2 3 1 2  A    1 4  1 4 1 2 3 1 2    1 4 1 2 3  试通过迭代法求解此方程组，认识迭代法收敛的含义以及迭代初值和方程组系数矩 阵性质对收敛速度的影响。实验要求： (0) (1) 选取不同的初始向量 和不同的方程组右端项向量 ，给定迭代误差要求， x b 用雅可比迭代法和高斯-赛德尔迭代法计算，观测得到的迭代向量序列是否均收 敛？若收敛，记录迭代次数，分析计算结果并得出你的结论。 (2) 取定右端向量 和初始向量x (0) ，将 的主对角线元素成倍增长若干次，非主 b A 对角线元素不变，每次用雅可比迭代法计算，要求迭代误差满足 (k 1) (k ) 5 x x 10 ，比较收敛速度，分析现象并得出你的结论。  根据雅可比迭代法 k 1 k     1 1 x B x  f ,B D L U , f D b J J J J 以及高斯-赛德尔迭代法 k 1 k 1 1 x   BGS x    f GS ,BGS D L  U , f GS D L  b 编写相应M 文件jacobi .m 和gau