线性方程组的直接解法迭代解法.docx

PAGE \* MERGEFORMAT13 广东金融学院实验报告 课程名称：数值分析 实验编号 及实验名称 线性方程组的直接解法/迭代解法系 别应用数学系姓 名学 号班 级实验地点实验日期实验时数2指导教师同组其他成员无成 绩实验目的及要求 实验目的： 题一：通过数值实验，从中体会解线性方程组选主元的必要性和LU分解法的优点，以及方程组系数矩阵和右端向量的微小变化对解向量的影响。比较各种直接接法在解线性方程组中的效果； 题二：认识各种迭代法收敛的含义、影响各迭代法收敛速度的因素。 实验要求： 题一： （1）在MATLAB中编写程序用列主元高斯消去法和LU分解求解上述方程组，输出Ax=b中矩阵A及向量b和 A=LU分解中的L及U，detA及解向量x. （2）将方程组中的2.099999改为2.1，5.900001改为5.9，用列主元高斯消去法求解变换后的方程组，输出解向量x及detA，并与（1）中的结果比较。 （3）用MATLAB的内部函数inv求出系数矩阵的逆矩阵，再输入命令x=inv(A)*b,即可求出方程组的解。请与列主元高斯消去法和LU分解法求出的解进行比较，体会选主元的方法具有良好的数值稳定性。用MATLAB的内部函数det求出系数行列式的值，并与（1）、（2）中输出的系数行列式的值进行比较。 （4）比较以上各种直接解法在解线性方程组中的效果。 题二： 选取不同的初始向量X及右端向量b，给定迭代误差要求，用Jacobi迭代法和 Gauss-Seidel迭代法求解，观察得到的序列是否收敛？若收敛，记录迭代次数，分析计算结果并得出你的结论。列出算法清单。 （2）用SOR迭代法求上述方程组的解，松弛系数取的不同的三个值，在时停止迭代，记录迭代次数，分析计算结果与松弛系数的关系并得出你的结论。 （3）用MATLAB的内部函数inv求出系数矩阵的逆矩阵，再输入命令x=inv(A)*b,即可求 出上述各个方程组的解，并与上述三种方法求出的解进行比较。请将比较结果列入下表。  方程组的解迭代次数误 差精 确 解Jacibi解法Gause-seidel解法SOR 解法 实验环境及相关情况（包含使用软件、实验设备、主要仪器及材料等） 1.Win7 2.Matlab 7.0 实验内容及步骤（包含简要的实验步骤流程） 实验内容： 题一： 解下列线性方程组 题二： 研究解线性方程组迭代法的收敛性、收敛速度以及SOR方法中最佳松弛因子的选取问题，用迭代法求解 ，其中 ， 实验结果（包括程序或图表、结论陈述、数据记录及分析等，可附页） 题一：直接解法解线性方程组 列主元高斯消去法与LU分解求解 列主元高斯消去法: 编写matalab程序（见附录gaosi.m），输出矩阵 A=10.000-7.0000.0002.5000 0.0001-5.000-1.50.0000.0000.0000.000 6.0000022.3000.0005.080 向量 b=818.3005.0800 解向量： X=(0，-1，1，1)T 其中系数行列式的值det(A)=762.00009 LU分解求解： 编写matalab程序（见附录zhjLU.m和LU.m），执行输出: L=1.00000.0000-3.0001.00000 0.00000.00000.00000.00000.50002000-2400000 1.00000.00000.96001.0000 U=10.0000-7.00000.0000-0.000001 0.00001.00000.00002.30000.00000.00000.00000.0000 1500000057500000.00005.0800 在matlab命令窗口输入L*U ，可以得到A= L*U ,即分解结果正确。 解向量: X=(0，-1，1，1)T 其中系数行列式的值det(A)=-762.00009 将方程组中的2.099999改为2.1，5.900001改为5.9,使用列主元法求解后与（1）比较结果 修改程序（见附录gaosi2.m），输出解向量： X=(0，-1，1，1)T 其中系数行列式的值det(A)=762，与(1)中结果对比： 方法解向量|det(A)|列主元高斯消去法X=(0，-1，1，1)T762.00009微小变化后列主元法X=(0，-1，1，1)T762 从上表可以