5线性代数方程组数值解法.docxVIP

PAGE / NUMPAGES 实验5 线性代数方程组的数值解法 化工系毕啸天 2018011811 【实验目的】 1. 学会用MATLAB 软件数值求解线性代数方程组，对迭代法的收敛性和解的稳定性作初步分析； 2. 通过实例学习用线性代数方程组解决简化的实际问题。 【实验内容】 题目3 已知方程组Ax=b，其中A，定义为 试通过迭代法求解此方程组，认识迭代法收敛的含义以及迭代初值和方程组系数矩阵性质对收敛速度的影响。实验要求： 1）选取不同的初始向量x(0和不同的方程组右端项向量b ，给定迭代误差要求，用雅可比迭代法和高斯-赛德尔迭代法计算，观测得到的迭代向量序列是否均收敛？若收敛，记录迭代次数，分析计算结果并得出你的结论； 2）取定右端向量b和初始向量x(0，将A的主对角线元素成倍增长若干次，非主对角线元素不变，每次用雅可比迭代法计算，要求迭代误差满足，比较收敛速度，分析现象并得出你的结论。 3.1 模型分析 选取初始向量x(0=(1,1,…,1T，b=(1,1,…,1T，迭代要求为误差满足，编写雅各比、高斯-赛德尔迭代法的函数，迭代求解。 3.2 程序代码 function x = Jacobi( x0,A,b,m D=diag(diag(A。 U=-triu(A,1。 L=-tril(A,-1。 B1=D\(L+U。 f1=D\b。 x(:,1=x0。 x(:,2