数值分析与实验

1、数值分析与实验作者：wlf5775196时间：2010-3-12 10:08:00 论文关键词： 列主元高斯消去法 雅可比法 高斯赛德尔迭代法 幂法论文摘要：本文通过实例对线性方程组数值解法和矩阵的特征值及特向量的计算进行了探讨。在对线性方程组数值解法的讨论下用到了列主元高斯消去法、雅可比法和高斯赛德尔迭代法。正是高斯消去法在消元时存在一些必须的条件，才启发我们通过列主元高斯消去法来对线性方程组数值解法作进一步的研究，达到了很好的的效果。同时用雅可比法和高斯赛德尔迭代法对相类似的问题的探讨来比较它们的优劣，使我们在分析问题时能更好的把握方法。在求矩阵按模最大的特征值及对应特征向量时，本文用到了

2、幂法，可以使现实中很多复杂的计算简单。第一章：线性方程组数值解法实验目的熟悉求解线性方程组的有关理论和方法 ；会编制列主元消去法，雅可比及高斯-赛德尔迭代法的程序 ；通过实际计算，进一步了解各种方法的优缺点，选择合适的数值方法。实验内容列主元高斯消去法求解线形方程组；雅可比法和高斯赛德尔迭代法解方程组；1.1 题目：列主元高斯消去法求解线形方程组 方程组为： 将方程用增广矩阵 表示1） 消元过程对k=1，2，.,n-1 1 选主元，找 使得 2 如果 则矩阵A奇异，程序结束；否则执行33 如果 则交换第k行与第 行对应元素位置， j=k,n+14 消元，对i=k+1,n计算 对j=k+1,n+

3、1计算 2） 回代过程1 若 则矩阵A奇异，程序结束；否则执行22 ；对i=n-1,2,1计算 #include #include void ColPivot(float *c,int n,float x) int i,j,t,k; float p; for(i=0;i=n-2;i+)k=i;for(j=i+1;j(fabs(*(c+k*(n+1)+i)k=j; if(k!=i) for(j=i;j=n;j+) p=*(c+i*(n+1)+j); *(c+i*(n+1)+j)=*(c+k*(n+1)+j); *(c+k*(n+1)+j)=p; for(j=i+1;j=n-1;j+) p=(*(

4、c+j*(n+1)+i)/(*(c+i*(n+1)+i);for(t=i;t=0;i-) for(j=n-1;j=i+1;j-) (*(c+i*(n+1)+n)-=xj*(*(c+i*(n+1)+j); xi=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i); void main()void ColPivot(float*,int,float);int i;float x4;float c45=1,-1,2,-1,-8,2,-2,3,-3,-20,1,1,1,0,-2,1,-1,4,3,4,;ColPivot(c0,4,x);for(i=0;i=3;i+)printf(x%d=%f

5、n,i,xi);从输出结果可以得到 =-6.999999， =3.000000， =2.000000， =2.000000从结果和过程可以知道这种方法一般能保证舍入误差不扩散，这个方法基本上是稳定的。1.2 题目 雅可比法解方程组方程组为： 设方程组Ax=b的系数矩阵的对角线元素 （i=1,2,n），M为迭代次数容许的最大值 为容许误差。1 取初始向量 令k=0.2 对i=1,2,n计算 3 如果 则输出 结果；否则执行44 如果 则不收敛，终止程序；否则 ，转2#include #include #define eps 1e-6#define max 100void Jacobi(float

6、 *a,int n,float x) int i,j,k=0; float epsilon,s; float *y= new float n; for(i=0;in;i+) xi=0; while(1) epsilon=0; k+; for(i=0;in;i+) s=0; for(j=0;jn;j+) if(j=i)continue; s+=*(a+i*(n+1)+j)*xj; yi=(*(a+i*(n+1)+n)-s)/(*(a+i*(n+1)+i); epsilon+=fabs(yi-xi);for(i=0;in;i+)xi=yi;if(epsilon=max)printf(die dai

7、 fa san);return;delete y;void main()sint i;float a45=10,-1,2,0,-11,0,8,-1,3,-11,2,-1,10,0,6,-1,3,-1,11,25;float x4;Jacobi(a0,4,x);for(i=0;i4;i+)printf(x%d=%fn,i,xi);迭代次数增加时，精度越高。从输出结果可以看出此方程组的迭代次数为17，迭代结果越来越接近精确解了，于是 =-1.467391, =-2.358696, =0.657609, =2.8423911.3 题目 高斯-赛德尔迭代法解方程组方程组为： 设方程组Ax=b的系数矩阵

8、的对角线元素 （i=1,2,n），M为迭代次数容许的最大值 为容许误差。1 取初始向量 令k=0.2 对i=1,2,n计算 3 如果 则输出 结束；否则执行44 如果 则不收敛，终止程序；否则 ，转2#include #include #define N 600void main()int i;float x4;float c45=10,-1,2,0,-11,0,8,-1,3,-11,2,-1,10,0,6,-1,3,-1,11,25;void GaussSeidel(float *,int,float);GaussSeidel(c0,4,x);for(i=0;i=3;i+)printf(x%

9、d=%fn,i,xi);void GaussSeidel(float *a,int n,float x) int i,j,k=1; float d,dx,eps; for(i=0;in-1;i+)xi=0.0; while(1)eps=0; for(i=0;in-1;i+) d=0; for(j=0;jn-1;j+) if(j=i)continue; d+=*(a+i*(n+1)+j)*xj; dx=(*(a+i*(n+1)+n)-d)/(*(a+i*(n+1)+i);eps+=fabs(dx-xi);xi=dx;if(epsN)printf(迭代发散nn);return;k+;从输出结果可以

10、看出此方程组的迭代次数为7，此时能得到精确结果是 =-1.377632, =-1.281579, =0.747368, =-107374176从结果和原有知识可以知道其系数矩阵是严格对角占优的。所以此迭代解法有很好的收敛性。1.4 方法比较雅可比法和高斯赛德尔迭代法解方程组两种方法的比较 。由于此题的系数矩阵是严格对角占优的，所以雅克比迭代法和高斯-赛德尔迭代法都是收敛的，这两种迭代法没迭代一步均是作一次矩阵和向量的乘法，但前者需要2组工作单元分别存放 和 ，而后者只需要1组工作单元。对于同一个线性方程组，这两种方法可能同时收敛，也可能同时发散，也可能其一收敛，而另一发散。但当两者皆收敛时，一

11、般来说高斯-赛德尔迭代法比雅克比迭代法收敛快。实际中更多的是使用逐次超松弛迭代法。第二章 矩阵的特征值及特征向量的计算实验目的在数学和物理中，很多问题都需要计算矩阵的特征值及特征向量，它们是线性代数中的一个重要课题，而在实际问题中，这样的计算是很复杂的，有的要求矩阵按模最大特征值及相应的特征向量，有些则要求全部特征值及特征向量，根据不同的要求计算方法大体上可分为2种类型。本实验用的是幂法求矩阵按模最大的特征值及对应特征向量，要求领会求矩阵特征值及特征向量的幂法的方法，并要求会编制幂法的计算程序，来计算有关问题。实验内容利用幂法求矩阵按模最大的特征值及对应特征向量。2.1 幂法求矩阵按模最大的特

12、征值及对应特征向量用幂法求矩阵按模最大的特征值 及其相应的特征向量 ，使 ， 幂法是求矩阵主特征值的一种迭代方法。设 有n个线性无关的特征向量 ，而相应的特征值满足 ，则对任意非零初始向量 按下述公式构造向量序列： 其中 表示 中最大的分量，并且有 ， 。 用幂法计算实对称矩阵的特征值时，可用Rayleigh商作加速。设 的Rayleigh商为 则， 当 时，将比 更快趋于 。主特征值为：98.521690；相应的特征向量为 。幂法是求矩阵主特征值的一种有效方法，特别当矩阵为大型稀疏（即矩阵元素中0元素较多）时，更显得如此。但由于特征值的分布无法事先预测，因此不能控制收敛速度，往往需要利用某些加速技巧。所以计算时我们要根据需要选择计算方法来