线性的方程组求解高高质量C语言程序.docVIP

实用标准文案 精彩文档 课 题： 线性方程组求解 目录 课题描述及要求··········································2 项目分析················································2 算法流程················································2 方法说明················································3 源代码··················································3 程序说明················································6 运行结果················································7 总结····················································7 参考文献················································8 线性方程组求解 陈龙 一．课题描述和功能要求 1.描述：求解线性方程组Ax=b，写成函数。其中，A为n乘n阶矩阵，x为n元未知向量，b为n个常数组成的矩阵。 2.要求：采用高斯先列主元消元法（也可采用其他方法）求解线性方程组AX=b。 二．项目分析 数学上，高斯消去法或称高斯－约当消去法，由高斯和约当得名（很多人将高斯消去作为完整的高斯－约当消去的前半部分），它是线性代数中的一个算法，用于决定线性方程组的解，决定矩阵的秩，以及决定可逆方矩阵的逆。当用于一个矩阵时，高斯消去产生“行消去梯形形式”。　　例如：一个二元一次方程组，设法对每个等式进行变形，使两个等式中的同一个未知数的系数相等，这两个等式相减，得到一个新的等式，在这个新的等式中，细数相等的未知数就被除去了（系数为0）。　　同样的也适合多元多次方程组。 我们知道m*n矩阵（用大写字母表示）是一个m行n列的数阵，n维向量（用加粗的小写字母表示）是n个数的数组，也就是一个n*1矩阵（列向量。我们不考虑行向量）。另外，大家也都知道矩阵乘法。因此一个m*n线性方程组可以表示为　　Ax=b，其中A是由系数aij组成的m*n矩阵即系数矩阵，x是n维的未知数向量，b是m维的结果向量。如果把向量b写到A的右边得到m*(n+1)的矩阵，得到的新矩阵称为这个方程组的增广矩阵。每一个方程组均对应于一个增广矩阵。 三．算法流程图 四．方法说明 （1）第1步消元——在增广矩阵（A，b）第一列中找到绝对值最大的元素，将其所在行与第一行交换 （2）第2步消元——在增广矩阵（A，b）中的第二列中（从第二行开始）找到绝对值最大的元素，将其所在行与第二行交换 （3）第3步消元——在增广矩阵（A，b）中的第三列中（从第三行开始）找到绝对值最大的元素，将其所在行与第二行交换 （4）按x4?- x3- x2- x1?的顺序回代求解出方程组的解。 五．源代码 #include stdio.h #include stdlib.h #include malloc.h #include math.h int GS(int,double**,double *,double); double **TwoArrayAlloc(int,int); void TwoArrayFree(double **); void main() { int i,j,n; double ep,**a,*b; ep = 1e-4; printf(有几个未知数); scanf(%d,n); a = TwoArrayAlloc(n,n); b = (double *)calloc(n,sizeof(double)); if(b == NULL) { printf(内存分布失败\n); exit(1); } for(i=0;in;i++) { printf(请输入第%d行相应的系数：\n,i+1); for(j=0;jn;j++) { printf(a[%d][%d]: ,i,j); scanf(%lf,a[i]+j); fflush(stdin); } printf(请输入第%d行相应的常数：\n,i+1); printf(b[%d]: ,i); scanf(%lf,b+i); fflush(stdin); } if(!GS(n,a,b,ep)) { p