数值计算方法实验5.doc

实验报告 学院（系）名称： 计算机与通信工程 姓名 学号专业 计算机科学与技术 班级 2012级1班 实验项目 实验五 线性方程组的解法 课程名称 数值计算方法 课程代码 0665026 实验时间 2015年6月23日第3-4节 2015年6月30日第3-4节 实验地点 7-219 批改意见： 成绩 教师签字： 实验目的： 编写列选主元高斯消去法、雅可比和高斯－塞德尔迭代法的算法程序，上机调试通过。 实验环境： 硬件环境：IBM－PC或兼容机 软件环境：Windows操作系统 编程语言：C语言 实验内容： 用列主元高斯消去法求解下列线性方程组 用高斯－塞德尔迭代方法求解线性方程组 要求误差为 实验步骤： 列主元高斯消去法 高斯消去法分为消元过程与回代过程，消元过程将所给方程组加工成上三角形方程组，再经回代过程求解 列选主元法 1. 列选主元是当消元到第k步时，从k列的akk以下（包括akk）的各元素中选出绝对值最大的，然后通过行交换将其交换到akk的位置上。交换系数矩阵中的两行（包括常数项）,只相当于两个方程的位置交换了，因此列选主元不影响求解的结果 2. 列选主元比全选主元的运算量小，但一般可以满足精度要求，所以列选主元更常被采用 3. 在列主元消去法中，未知数仍然是顺序地消去的，但是把各方程中要消去的那个未知数的系数按绝对值最大值作为主元素，然后用顺序消去法的公式求解 列主元高斯消去法的流程图： 主程序部分 列选主元部分 高斯-塞德尔迭代法 高斯－塞德尔迭代法的思想 迭代收敛时，新值xi (k+1) 比老值xi (k) 更准确 算出新值xi (k+1)后，用新值xi (k+1)代替用于后面计算的老值xi (k) ，使每次迭代计算，都是利用“最新求解信息” 这样，必然会使迭代求解的速度加快 高斯-塞德尔(Seidel)迭代公式 高斯－塞德尔迭代的流程图 实验结果： 列主元消去法 输入各个数据，最终使用列选主元法，得到结果为：x1=x2=x3=1 高斯-赛德尔迭代法 输入各个数据，输出每一步迭代数据，最终结果为：x1=0.285716，x2=2.142855，x3=-3.285712 附录（源程序及运行结果） 列主元高斯消去法 #includestdio.h #includemath.h void print(double a[3][3],int n,double b[3]){ printf(输出矩阵：\n); for(int i=0;in;i++){ for(int j=0;jn;j++) printf(%lf ,a[i][j]); printf(%lf,b[i]); printf(\n); } } void liezhuyuan(double a[3][3],double b[3],int k,int n){ double d,t; int l=k,i=k+1; d=a[k][k]; while(i!=n){ if(fabs(a[i][k])fabs(d)){ d=a[i][k]; l=i; } i++; } printf(选出主元：%lf\n,d); if(d==0) printf(矩阵奇异！\n); else if(l!=k){ for(int j=k;jn;j++){ t=a[l][j]; a[l][j]=a[k][j]; a[k][j]=t; } t=b[l]; b[l]=b[k]; b[k]=t; } } void main(){ double a[3][3]={2,1,1,1,3,2,1,2,2}, b[3]={4,6,5}; int k=0,i,j,n=3; print(a,n,b); while(k!=n-1){ liezhuyuan(a,b,k,n); for(i=k+1;in;i++) a[i][k]=a[i][k]/a[k][k]; for(i=k+1;in;i++){ for(j=k+1;jn;j++) a[i][j]=a[i][j]-a[i][k]*a[k][j]; b[i]=b[i]-a[i][k]*b[k]; } k++; } b[n-1]=b[n-1]/a[n-1][n-1]; for(i=n-2;i=0;i--){ for(j=i+1;jn;j++) b[