《计算方法》课内实验报告（实验三）.docVIP

PAGE 《计算方法》课内实验报告 学生姓名: 及 学 号： 学 院: 班 级: 课程名称： 计算方法 实验题目: 解线性方程组的直接方法和迭代法 指导教师 姓名及职称： 宋云飞 讲 师 朱秀丽 讲 师 尚宝欣 讲 师 2012年12月10日 目 录 一、实验题目 1 二、实验目的 1 三、实验内容 1 四、实验结果 2 五、实验体会或遇到问题 9 9 - 实验题目 解线性方程组的直接方法和迭代法 二、实验目的 1．熟悉Matlab编写及运行数值计算程序的方法。 2．进一步理解求解线性方程组的直接方法和迭代法基础理论。 3．进一步掌握应用不同的方法求解线性方程组的收敛速度及误差分析。 三、实验内容 1．用LU分解及列主元高斯消去法解线性方程组 . 输出中系数分解法的矩阵及，解向量及；列主元法的行交换次序，解向量及；比较两种方法所得的结果。 2．线性方程组的及为 则解.用Matlab内部函数求及的所有特征值和若令 ， 求解，输出向量和，从理论结果和实际计算两方面分析线性方程组解的相对误差及的相对误差的关系。 3．给出线性方程组，其中系数矩阵为希尔伯特矩阵： ， 假设，若取，分别用雅可比迭代法及SOR迭代法求解，比较计算结果。 四、实验结果 1. A=[10 -7 0 1;-3 2.099999 6 2;5 -1 5 -1;2 1 0 2] [L,U,y,x]=zhjLU(A,b) 结果： 请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解.A的秩RA和各阶顺序主子式值hl依次如下： 请注意：因为A的各阶主子式都不等于零，所以A能进行LU分解. L = 1.0e+006 * 0.0000 0 0 0 -0.0000 0.0000 0 0 0.0000 -2.5000 0.0000 0 0.0000 -2.4000 0.0000 0.0000 U = 1.0e+007 * 0.0000 -0.0000 0 0.0000 0 -0.0000 0.0000 0.0000 0 0 1.5000 0.5750 0 -0.0000 0 0.0000 y = 1.0e+007 * 0.0000 0.0000 2.0750 0.0000 x = 0.0000 -1.0000 1.0000 1.0000 列主元消去法： b=[8 5.9 5 1] b=b [RA,RB,n,X]=liezhu(A,b) 结果： 请注意：因为RA=RB=n，所以此方程组有唯一解. RA =4 RB =4 n =4 X = 0.0000 -1.0000 1.0000 1.0000 2. A=[10 7 8 7;7 5 6 5;8 6 10 9; 7 5 9 10]; B=det(A) B =1.0000 [V,D]=eig(A) V = 0.5016 -0.3017 0.6149 0.5286 -0.8304 0.0933 0.3963 0.3803 0.2086 0.7603 -0.2716 0.5520 -0.1237 -0.5676 -0.6254 0.5209 D = 0.0102 0 0 0 0 0.8431 0 0 0 0 3.8581 0 0 0 0 30.2887 C=cond(A,2) C =2.9841e+003 D=[10 7 8.1 7.2;7.08 5.04 6 5;8 5.98 9.89 9;6.99 5 9 9.98]; b=[32;23;33;31]; x=[1;1;1;1]; E=inv(D)*b; F=E-x; F = -10.5863 17.3741 -4.2258 2.5240