计算方法 实验 线性方程组的数值解法.docVIP

姓名： 学号： 班级： 线性方程组的数值解法 实习目的 （1）通过实习进一步掌握高斯消去法、列主元高斯消去法、柯朗分解法、追赶法以及雅克比迭代法和高斯-赛德尔迭代法的基本思想； （2）通过实习进一步掌握高斯消去法、列主元高斯消去法、柯朗分解法、追赶法以及雅克比迭代法和高斯-赛德尔迭代法的计算步骤，并能灵活应用； （3）通过对高斯消去法、列主元高斯消去法、柯朗分解法、追赶法以及雅克比迭代法和高斯-赛德尔迭代法的调试练习，进一步体会各种算法的特点； （4）通过对上机调试运行，组不培养解决实际问题的编程能力。 实习要求 （1）熟悉Turbo C的编译环境； （2）实习前复习高斯消去法、列主元高斯消去法、直接三角分解法、雅克比迭代法、高斯-赛德尔迭代法以及追赶法的基本思想和过程； （3）实习前复习高斯消去法、列主元高斯消去法、直接三角分解法、雅克比迭代法、高斯-赛德尔迭代法以及追赶法的计算步骤。 实习设备 （1）硬件设备：单机或网络环境下的微型计算机一台； （2）软件设备：DOS 3.3以上操作系统，Turbo C 2.0编译器。 实习内容 实习一 高斯消去法 （1）用高斯消去法求解线性方程组： （2）要求： ①将线性方程组写成用矩阵表示的形式，即Ax=b的形式。 ②输出系数矩阵的原始元素和经高斯消去法消去后的矩阵元素。 ③经高斯消去法消去后的矩阵是一个什么形式的矩阵？ ④请写出程序的运行结果。 #include stdio.h #define N 4 int Gauss(float a[N][N],float b[N]) { int i,j,k,flag=1; float t; for(i=0;iN-1;i++) { if(a[i][i]==0) { flag=0; break; } else { for(j=i+1;jN;j++) { t=-a[j][i]/a[i][i]; b[j]=b[j]+t*b[i]; for(k=i;kN;k++) a[j][k]=a[j][k]+t*a[i][k]; } } } return(flag); } void zg_matric(float a[N][N],float b[N]) { int i,j; for(i=0;iN;i++) { for(j=0;jN;j++) printf(%10f,a[i][j]); printf(%10f,b[i]); printf(\n); } printf(\n); } void main() { static float a[N][N]={{1,-4,6,-5},{-5,21,-33,32},{6,-26,43,-48},{5,-24,45,-64}}; float b[N]={8,-32,25,-10}; float x[N]={0,0,0,0}; int i,j,flag; zg_matric(a,b); flag=Gauss(a,b); if(flag==0) printf(Gauss method does not run.); else { x[N-1]=b[N-1]/a[N-1][N-1]; for(i=N-2;i=0;i--) { x[i]=b[i]; for(j=i+1;jN;j++) x[i]=x[i]-a[i][j]*x[i]; x[i]=x[i]/a[i][i]; } for(i=0;iN;i++) printf(x%d=%11.7f\n,i+1,x[i]); } } （3）思考题： ①高斯消去法的基本思想是什么？ 答：高斯消去法是最古老的求解线性代数方程组的方法之一，是消去法的一种特殊形式。 ②高斯消去法由哪两个过程组成？ 答：高斯消去法包括消元与回代两个过程。 实习二 列主元高斯消去法求解线性方程组 （1）用列主元高斯消去法求解线性方程组： （2）要求： ①输出系数矩阵的原始元素和经列主高斯消去法消去后的矩阵元素。 ②经列主高斯消去法消去后的矩阵是一个什么形式的矩阵？ ③请写出输出后的运行结果。 ④若线性方程组的准确解为：，请通过计算绝对误差证明，在一般情况下用经列主高斯消去法求解线性方程组是稳定的。 #includestdio.h #includemath.h #define N 3 void ColGauss(float a[N][N],float b[N]) { float t,max_fab; int i,j,k,l; for(i=0;iN-1;i++) { l=i; max_fab=fabs(