数学实验第3次作业线性方程组的数值解法.docxVIP

线性方程组的数值解法 一 实验目的 1 学会用MATLAB软件求解线性代数方程组，对迭代法的收敛性和解的稳定性做初步分析； 2 通过实例学习线性代数方程组解决简化的实际问题。 二 实验内容 1 已知方程组Ax=b，其中A∈R20×20，定义为： 试通过迭代法求解此方程组，认识迭代法收敛的含义以及迭代初值和方程组系数矩阵性质对收敛速度的影响。实验要求： 选取不同的初始向量x(0)和不同的方程组右端项向量b，给定迭代误差要求，用雅可比迭代法和高斯-赛德尔迭代法计算，观测得到的迭代向量序列是否均收敛？若收敛，记录迭代次数，分析计算结果并得出你的结论。 取定右端项向量b和初始向量x(0)，将A的主对角线元素成倍增长若干次，非对角线元素不变，每次用雅可比迭代法计算，要求迭代误差满足x(k+1)-x(k)∞10-5，比较收敛速度，分析现象并得出你的结论。 初步解决： 首先建立利用雅克比迭代法和高斯-赛德尔迭代法计算的函数M文件。 雅克比迭代法： 高斯-赛德尔迭代法： 关于迭代法的收敛性： 由原式Ax=b化简可得x=Bx+f，也就是说，研究此处迭代法是否收敛实际上就是研究有原矩阵A经变换之后所得矩阵B的性质。而由于该题第一问中，A保持不变，所以对于收敛性只需研究一次。在命令栏中输入以下命令： 构造题目中的矩阵A： 提取对角矩阵、上三角矩阵、下三角矩阵： 对于两种迭代方法分别求B： 分别计算两个矩阵的三种范数： n1和m1、n2和m2、n3和m3分别是矩阵B1、B2的1-范数、2-范数和∞范数。由书上定理，矩阵的谱半径不超过任何一个范数，即ρ(B)≤B，而由图中可以看出，两个矩阵的六个范数没有一个大于1，所以两个矩阵的谱半径一定都小于1，所以此时两种迭代方法均收敛。 利用两个函数研究原方程组： 对于相同的初始值x(0)，选取不同的右端项向量b。 令x(0)=(0…0)20T，b分别取(1 2…19 20)20T、(20 19…2 1)20T、(10 10…10 10)20T、(-10 -9…8 9)20T，分别观察两种迭代方法最后的迭代次数。(迭代误差拟定为0.001) k(1 2…19 20)20T(20 19…2 1)20T(10 10…10 10)20T(-10 -9…8 9)20TJacobi13131312Gauss-S10998相关截图如下： 从表中的数据可以看出，当初始值不变，右端项向量变化时，使用雅克比迭代法和高斯-赛德尔迭代法都可以得到最后的结果，但是明显高斯-赛德尔迭代法比雅克比迭代法的收敛速度要快。从表中还可以看出，b的变化是会影响迭代次数的，因为b变化以后方程组的解发生了变化，使得计算得过程也发生了变化，进而影响了迭代次数。 对于相同的右端项向量b，选取不同的初始值x(0)。 令b=(1 2…19 20)20T，x(0)分别取(0…0)20T、(1…1)20T、(-1…-1)20T、(1 2…19 20)20T，分别观察两种迭代方法最后的迭代次数。(迭代误差拟定为0.001) k(0…0)20T(1…1)20T(-1…-1)20T(1 2…19 20)20TJacobi13131413Gauss-S109109相关截图如下： 从表中可以看出，初始值的变化会影响迭代次数，但是当迭代误差不大的时候，初始值的变化对迭代次数的变化的影响也不会非常大。但是通过实验可以发现，将最后一组数据的迭代误差改为0.0001是，用雅克比迭代法的迭代次数变为16，而高斯-赛德尔迭代法的迭代次数变为11。但是总体而言，高斯-赛德尔迭代法的收敛速度还是比雅克比迭代法的收敛速度快。 下面研究第二问。将初始值x(0)定为(1…1)20T，右端向量b定为(1 2…19 20)20T。改变的A的主对角元素，分别变为原来的2倍、3倍、5倍、10倍，观察迭代次数变化。(用雅克比迭代法，迭代误差定为10-5) 主对角元素691530k10865相关截图如下： 从表格中的数据可以看出，当改变主对角元素的值时，对迭代次数的影响很大，而且在此时给出的条件下，当主对角元素的值增大时，迭代次数同步的减小。 2 对于5.1.2节假设c=10，a1=0.5，a2=0.25，b=0.20，第一年有50颗植物，且第50年后有600颗植物。试分别用追赶法、稀疏系数矩阵和满矩阵求解；若b有10%的误差，估计对结果的影响。 初步解决：用三种不同的方法解决问题。首先对数据做简单处理，由题目中c=10，a1=0.5，a2=0.25，b=0.20，第一年有50颗植物，且第50年后有600颗植物，这些条件可知，p=-a1bc=-1，q=-a2b1-a1bc=-0.05，n=50，x0=50，xn=