线性方程组程序.docVIP

线性方程组程序 //数值计算实验，解线性方程组 #include iostream #include cmath #include cstdlib #include ctime #include vector #include mylib/InitData.h #include mylib/MyClock.h #include gsl/gsl_linalg.h using namespace std; //初始化系数矩阵 void Init(vectorvectordouble A) { A.clear(); int n; cout输入方程组的阶数:; cinn; if (n = 1) { cout方程组的阶数必须大于1endl; exit(1); } for (int i = 0; i n; i++) { vectordouble temp(n+1); for (int j = 0; j n+1; j++) { cout输入第i个方程的第j个系数:; double a = 0.0; cina; temp[j] = a; } A.push_back(temp); temp.clear(); } } //随机初始化系数矩阵 void RandInit(vectorvectordouble A, const int n) { A.clear(); for (int i = 0; i n; i++) { vectordouble temp; InitData(temp, n+1, 1, 100); A.push_back(temp); temp.clear(); } } //输出系数矩阵 void DisplayData(const vectorvectordouble A) { int n = A.size(); int m = A[0].size(); for (int i = 0; i n; i++) { ShowData(A[i], cout); } } //输出解向量 void DisplaySolver(const vectordouble solve) { int n = solve.size(); for (int i = 0; i n; i++) { coutxi = solve[i]endl; } } //高斯消去法 //A为系数矩阵，x为解向量，若成功，返回true，否则返回false，并将x清空。 bool Gauss(const vectorvectordouble A, vectordouble x) //由于输入函数已经保证了系数矩阵是对的，就不检查A了 int n = A.size(); int m = A[0].size(); x.resize(n); //复制系数矩阵，防止修改原矩阵 vectorvectordouble Atemp(n); for (int i = 0; i n; i++) { vectordouble temp(m); for (int j = 0; j m; j++) { temp[j] = A[i][j]; } Atemp[i] = temp; temp.clear(); } //第一次消元 for (int i = 1; i n; i++) { double l = Atemp[i][0]/Atemp[0][0]; for (int j = 0; j m; j++) { Atemp[i][j] = Atemp[i][j] - l*Atemp[0][j]; } } //剩余的消元 for (int k = 1; k n-1; k++) { for (int i = k+1; i n; i+