首先將
方程組中的
係數矩陣
a分解成三部分,即:
a = l+d+u,如圖1所示,其中
d為對角陣,
l為下三角矩陣,
u為上三角矩陣。
之後確定迭代格式,x^(k+1) =
b*x^(k) +
f ,(這裡^表示的是上標,括號內數字即迭代次數),如圖2所示,其中
b稱為迭代矩陣,雅克比迭代法中一般記為
j。(k = 0,1,......)
再選取初始迭代
向量x^(0),開始逐次迭代。
設ax= b
,其中a=d+l+u
為非奇異矩陣
,且對角陣
d也非奇異,則當迭代矩陣j的
譜半徑ρ(j)
<1時,雅克比迭代法收斂。
#include
#include
#include
main(),b[3]=,x[3]=,y[3];
int n=3,j,i,k=1;
while(1)
i=0;
while(i<3)
if(i!=3)
else if(i==3)
break; }
printf("%f\n%f\n%f\n",y[0],y[1],y[2]); }
求解線性方程組 Jacobi迭代法
求解線性方程組 jacobi迭代法 求解線性方程組 jacobi迭代法 include include using namespace std 二維陣列動態分配模板 template t allocation2d int m,int n return a 一維陣列動態分配模板 template t ...
迭代法求解非線性方程的根
1 求解的是非線性方程放f x 0的根 1 判斷使用全域性收斂性定理還是區域性收斂性定理 全域性收斂性 未知根,已知根的取值區間 定理1.1 1.2 區域性收斂性 已知根,即已知根的鄰域 定理1.3 1.4 2 解題步驟 根的存在性 零點存在定理 唯一性 單調性 反證法 迭代格式的收斂性 一般用 a...
牛頓迭代法求解非線性方程 C
一 牛頓迭代法演算法 設 已知 第1步計算函式 第2步計算雅可比矩陣 第3步求線性方程組的解 第4步計算下乙個點 重複上述過程。二 牛頓迭代法c 程式 牛頓迭代解非線性方程組 include include using namespace std double f1 double x double ...