jacobi迭代法求解線性方程組

2021-08-07 20:39:39 字數 674 閱讀 8614

首先將

方程組中的

係數矩陣

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 ...