数值分析与实验.docVIP

数值分析与实验 论文关键词：列主元高斯消去法雅可比法高斯－赛德尔迭代法幂法 论文本文实例对线性方程组数值解法和矩阵的特征值及特向量的计算了探讨。在对线性方程组数值解法的下用到了列主元高斯消去法、雅可比法和高斯－赛… 论文关键词： 列主元高斯消去法 雅可比法 高斯－赛德尔迭代法 幂法 论文本文实例对线性方程组数值解法和矩阵的特征值及特向量的计算了探讨。在对线性方程组数值解法的下用到了列主元高斯消去法、雅可比法和高斯－赛德尔迭代法。正是高斯消去法在消元时的条件，才启发列主元高斯消去法来对线性方程组数值解法作的，了的的。用雅可比法和高斯－赛德尔迭代法对相类似的问题的探讨来它们的优劣，使在分析问题时能的把握方法。在求矩阵按模最大的特征值及对应特征向量时，本文用到了幂法，可以使现实中的计算简单。 章：线性方程组数值解法 实验目的 熟悉求解线性方程组的理论和方法 ；会编制列主元消去法，雅可比及高斯-赛德尔迭代法的程序 ；计算，方法的优缺点，选择合适的数值方法。 实验内容 列主元高斯消去法求解线形方程组； 雅可比法和高斯－赛德尔迭代法解方程组； 1.1 题目：列主元高斯消去法求解线形方程组 方程组为： 1.1.1 列主元高斯消去法算法 将方程用增广矩阵 表示 1） 消元过程 对k=1，2，….,n-1 1 选主元，找 使得 2 则矩阵A奇异，程序结束；否则3 3 则交换第k行与第 行对应元素位置， j=k,…,n+1 4 消元，对i=k+1,…,n计算 对j=k+1,…,n+1计算 2） 回代过程 1 若 则矩阵A奇异，程序结束；否则2 2 ；对i=n-1,…2,1计算 1.1.2 程 序 #inlude stdi.h #inlude ath.h vid lPivt(flat *,int n,flat x[]) { int i,j,t,k; flat p; fr(i=0;i =n-2;i++) {k=i; fr(j=i+1;j =n-1;j++) if(fabs(*(+j*(n+1)+i)) (fabs(*(+k*(n+1)+i))))k=j; if(k!=i) fr(j=i;j j++) { p=*(+i*(n+1)+j); *(+i*(n+1)+j)=*(+k*(n+1)+j); *(+k*(n+1)+j)=p; } fr(j=i+1;j =n-1;j++) { p=(*(+j*(n+1)+i))/(*(+i*(n+1)+i)); fr(t=i;t t++)*(+j*(n+1)+t)-=p*(*(+i*(n+1)+t)); } } fr(i=n-1;i i--) { fr(j=n-1;j =i+1;j--) (*(+i*(n+1)+n))-=x[j]*(*(+i*(n+1)+j)); x[i]=*(+i*(n+1)+n)/(*(+i*(n+1)+i)); } } vid ain() { vid lPivt(flat*,int,flat[]); int i; flat x[4]; flat [4][5]={1,-1,2,-1,-8,2,-2,3,-3,-20,1,1,1,0,-2,1,-1,4,3,4,}; lPivt([0],4,x); fr(i=0;i i++)printf(x[%d]=%f\n,i,x[i]); } 1.1.3 输出结果 1.1.4结果分析 从输出结果可以 =-6.999999， =3.000000， =2.000000， =2.000000 从结果和过程可以知道方法能舍入误差不扩散，方法上是稳定的。[1][2][3][4][5]下一页 【 1.2 题目 雅可比法解方程组 方程组为： 1.2.1 雅可比迭代法算法 设方程组Ax=b的系数矩阵的对角线元素 （i=1,2,…,n），为迭代次数容许的最大值 为容许误差。 1 取初始向量 令k=0. 2 对i=1,2,…,n计算 3 则输出 结果；否则4 4 则不收敛，终止程序；否则 ，转2 1.2.2 程 序 #inlude stdi.h #inlude ath.h #define eps 1e-6 #define ax 100 vid Jabi(flat *a,int n,flat x[]) { int i,j,k=0; flat epsiln,s; flat *y= ne flat [n]; fr(i=0;i i++) x[i]=0; hile(1) { epsiln=0; k++; fr(i=0;i i++) { s=0; fr(