高丹计算方法高斯.docVIP

对于高斯消去法、高斯列主元消去法、高斯全主元消去法的分析与实现 ——计算方法结业设计 高丹 高斯（Gauss）消元法 1.高斯消元法的基本思想 高斯（Gauss）消元法的基本思想是：按照线性方程组给定的方程和变量顺序进行消元，将线性方程组化为等价的上三角形方程组，然后再通过回代法求出此线性方程组的解。 例 用高斯消元法解方程组解 第一次消元： ，得 第二次消元： ，得 最后，原方程组的解为 2高斯消元法的计算公式 用高斯消元法解阶线性方程组 设，则以乘以第1个方程加到第个方程，则方程组化为 第2次消元：设，则以乘以第个方程加到第个方程，则方程组化为 第n-2次消元后化为同解方程组： 第n-1次消元：设，则以乘以第个方程加到第个方程，则方程组化为 从原方程组化为上三角形方程组的过程，称为消元过程。求解上三角形方程组的过程，称为回代过程。消的计算为 回代过程的计算公式为：　消元过程和回代过程一起，组成了解阶线性方程组的高斯消元法 高斯消元法可实行的条件是元素全不为零。高斯消元法的总运算量为 消元的计算为 第k次消元设 以乘以第个方程加到第个方程 即：对的每个i 对,执行 注:用代替 回代过程如下: 根据最后一个方程直接可得 如此再解出,一般有 1.4高斯消去法程序: #include iostream #include vector #include fstream #include cmath #include iomanip #includestring using namespace std; class GS{ public: GS(string fn=string()) { if(fn.size()==0){ cout 请输入方程组系数增广矩阵文件名：; cinfn; } //fn文件名 ifstream fin(fn.c_str()); if(!fin) { throw fn+文件不存在!; } finN; A.resize(N+1); for(int i=0;iN+1;i++) A[i].resize(N+2); x.resize(N+1); for(int ii=1;ii=N;ii++) for(int j=1;j=N+1;j++) finA[ii][j]; } void print() { coutN阶endl; for(int i=1;i=N;i++){ for(int j=1;j=N+1;j++) coutfixedsetw(12) A[i][j] ; coutendl; } cout*************************endl; } void go() { for(int k=1;k=N-1;k++) { if((fabs(A[k][k])-1e-6)0) throw string(主元为0，失败！); for(int i=k+1;i=N;i++) { for(int j=k+1;j=N+1;j++) A[i][j]+=(-1)*A[i][k]/A[k][k] * A[k][j]; A[i][k]=0; } print(); } //回代 if(fabs(A[N][N])1e-6) throw string(主元为0，失败！); x[N]=A[N][N+1]/A[N][N]; for(int kk=N-1;kk=1;kk--) { double s=0.0; for(int j=k+1;j=N;j++) s=s+A[kk][j]*x[j]; x[k]=(A[kk][N+1]-s)/A[kk][kk]; } for(int i=1;i=N;i++) coutX[i]=x[i]endl; } private: int N; vectorvectordouble A; vectordouble x; }; int main() { try { GS f(); f.print(); f.go(); } catch(string s) { cout s; } catch(int n) { coutn; } return 0; } 2.高斯列主元消元法