有點(diǎn)線代的知識(shí)：

const double EPS = 1e-8; typedef vector<double> vec; typedef vector<vec> mat; //Ax = b vec gauss_jordan(const mat &A, const vec &b) {int n = A.size();mat B(n,vec(n+1)); //定義大小for(int i=0;i < n;i++)for(int j=0;j < n;j++)B[i][j] = A[i][j];for(int i=0;i < n;i++) B[i][n] = b[i]; //把 b 放在 B 的右邊一起運(yùn)算for(int i=0;i < n;i++) //把正在處理的未知數(shù)系數(shù)的絕對值最大的式子換到第i行（為什么是第i行）{int pivot = i;for(int j=i;j < n;j++)if(abs(B[j][i]) > abs(B[pivot][i]))pivot = j;swap(B[pivot],B[i]);if(abs(B[i][i]) < EPS) return vec(); //這種情況下，無解或者是無窮多解。for(int j=i+1;j <= n;j++) B[i][j] = B[i][j]/B[i][i]; //把正在處理的未知數(shù)的系數(shù)變?yōu)?。for(int j=0;j < n;j++){if(i != j)for(int k=i+1;k <= n;k++)B[j][k] -= B[j][i]*B[i][k]; //消參}}vec x(n);for(int i=0;i < n;i++) x[i] = B[i][n];return x; }

先了解一下，以后拿例題練。

轉(zhuǎn)載于:https://www.cnblogs.com/cunyusup/p/8452904.html

總結(jié)

以上是生活随笔為你收集整理的高斯消元法的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。