參考：
// jacobi.c ://雅克比迭代演算法
#define _crt_secure_no_warnings
#include#include#include#define dim 10
//jacobi迭代法為整體代換法，n為係數矩陣a的階數，b為右端常向量，x0為迭代初始值，
// a(n×n)*x0(n×1)=b(n×1)
//maxnum為最大迭代次數，方程組的解返回在陣列x1。
void jacobi(int n, long double a[dim][dim], long double b[dim], long double x0[dim], int maxnum, long double *x1)
fprintf(fp, "\n");
for (k = 1; k <= maxnum; k++)
}for (i = 1; i <= n; i++)
//換行
fprintf(fp, "\n");
} //關閉檔案
fclose(fp);
}//主程式
int main()
//最大迭代次數
maxnum = 100;
//呼叫jacobi迭代法求解過程
jacobi(n, a, b, x0, maxnum, x);
//輸出方程組的解
結果：
2020 03 31 思維永遠只有一種
08 30 乙個人沒有愛情的時候可以滿足於虛榮，但一旦有了愛情，虛榮就變得庸俗不堪了。小仲馬 何太看到了我出門前後那魂不守舍的狀態，各種的出門的遺漏，便問我是不是因為這些資訊讓自己焦慮了。當何太說出焦慮二字時，我的第一反應是在反省我自己到底是不是在焦慮。我嘗試跳出當前的狀態，重新思考了一番早上出門前...
數值分析Jacobi迭代法c c 實現
matlab還沒學精通，正好好久沒寫 了，藉此溫習溫習。include jacobi.h include include include using namespace std using std vector vector double matrix a vector double vector ...
C 解方程組之Jacobi迭代法
迭代過程 首先將 方程組中的 係數矩陣 a分解成三部分，即 a l d u，如圖1所示，其中 d為對角陣，l為下三角矩陣，u為上三角矩陣。之後確定迭代格式，x k 1 b x k f 這裡 表示的是上標，括號內數字即迭代次數 如圖2所示，其中 b稱為迭代矩陣，雅克比迭代法中一般記為 j。k 0,1，...