数值线性的代数课设课件资料

《数值线性的代数课设课件资料》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案数值线性代数课程设计报告姓名:陶英学号:081410124任课教师:杨熙南京航空航天大学2016年6月22日精彩文档实用标准文案求解线性方程组的三种迭代法及其结果比较摘要当今的环境下，数值计算越来越依赖于计算机。大规模科学计算和工程技术中许多问题的解决，最终归结为大型稀疏线性方程组的求解，其求解时间在整个问题求解时间中占有很大的比重，有的甚至达到80%。由于现今科学研究和大型项目中各种复杂的可以对计算精度和计算速度的要求越来越高。因此，作为大规模科学计算基础的线性代数方程组的高效数值求解引起了人们的普遍关注。这种方程组的求解一般采用迭代法。关于迭代法，是有很多

2、种解决公式的：Jacobi，G-S和超松弛迭代法。这三种方法的原理大致相同，Jacobi需要给定初向量，G-S则需要给定初值，超松弛法是对Guass-Seidel迭代法的加权平均改造。而本文则是对大型稀疏线性方程组迭代求解与三种迭代法（Jacobi，Gauss-Seidel和超松弛迭代法）的收敛速度与精确解的误差比较做出研究。关键词：Jacobi迭代法；Gauss-Seidel迭代法；SOR迭代法；线性方程组精彩文档实用标准文案1方法与理论的叙述1.1迭代法简介1.Jacobi迭代法：对于非奇异线性方程组Ax=b，令A=D-L-U，其中则原方程组可改写为：（2.2）其中给

3、定初始向量：由（2.2）可以构造迭代公式：其分量形式为：精彩文档实用标准文案2.Guass-Seidel迭代法：类似于Jacobi迭代法，给定初值：令则得到Guass-Seidel公式：其分量形式为：3.超松弛迭代法（SOR迭代法）：SOR迭代法是对Guass-Seidel迭代法的加权平均改造，即精彩文档实用标准文案为Guass-Seidel迭代解，即它的分量形式为：其中ω称为松弛因子，当ω>1时称为超松弛；当ω<1时叫低松弛；ω=1时就是Guass-Seidel迭代。上述三种经典迭代法收敛的充分必要条件是迭代矩阵谱半径小于1。谱半径不易求解，而在一定条件下，通过系数矩阵

4、A的性质可判断迭代法的收敛性。定理1：若系数矩阵A是严格对角占优或不可约对角占优，则Jacobi迭代法和Gauss-Seidel迭代法均收敛。定理2：（1）SOR迭代法收敛的必要条件是0

5、S，超松弛迭代法分别求线性方程组的解，要求4位有效数字，然后比较与精确解的误差。对=0.1，=0.01，=0.001，考虑同样问题。1.方程的表示及存储由于本题中线性方程组的系数矩阵为三对角矩阵，所以可以采用紧缩方法存储，即然后在矩阵乘法时对下标处理一下即可。但是考虑到三种迭代方法的一般性，且本题中n=200并不是很大，所以实验中并没有采用紧缩存储，而是采用了直接存储。2.边值条件的处理由于差分得到的方程组的第一行和最后一行中分别出现了边值y(0)与y(1)作为常数项，因此要在常向量的第一项和最后一项作一些修改：精彩文档实用标准文案3.迭代终止条件首先确定要求的精度tol

6、，我们希望当则停止迭代。对于迭代格式，若且，则迭代序列的第k次近似解和精确解之间有估计式。由题目要求知我们需要有，而由上面的迭代估计，只要，即即可。而本题中q可近似取为，因此最后令迭代终止条件为4.SOR迭代中最佳松弛因子的选取由于SOR迭代法的效果和其松弛因子w的选取有关，所以有必要选取合适的松弛因子。当选择最佳松弛因子时，SOR方法的迭代速度最快。Matlab实现：迭代矩阵是n-1阶的，不是n阶；等号右端向量b的最后一项，不是ah^2，而是ah^2-eps-h精彩文档实用标准文案2.2精确解带入a=1/2,=1代码：>>clear>>x=linspace(0,1);t

7、ruy=(1-0.5)/(1-exp(-1/1))*(1-exp(-x./1))+x.*0.5;figure;plot(x,truy,'g','LineWidth',1.5);holdon;Grid图：精彩文档实用标准文案2.3三种迭代法Jacobi法：代码见附录Eps=1结果：迭代次数k：22273结果与精确解的比较图（绿色粗线是精确解，黑色细线是迭代结果）Eps=0.1结果：迭代次数k：8753结果与精确解的比较图（绿色粗线是精确解，黑色细线是迭代结果）精彩文档实用标准文案Eps=0.01结果：迭代次数k：661结果与精确