資源簡(jiǎn)介
matlab代碼,Jacobi迭代算法,親測(cè)可用,直接調(diào)用該函數(shù)即可。
代碼片段和文件信息
%雅克比迭代法:建立函數(shù)文件Jacobi.m
function?Jacobi(AbXnmw_error)
%用雅克比迭代法求解方程組Ax=b
%輸入:A為方程組的系數(shù)矩陣,b為方程組右端的列向量,X為迭代初值構(gòu)成的列向量,nm為最大迭代次數(shù),w_error為誤差精度
%輸出:x為求得的方程組的解構(gòu)成的列向量,n為迭代次數(shù)
n=1;
m=length(A);
D=diag(diag(A));????????%令A(yù)=D-L-U計(jì)算矩陣D
L=tril(-A)+D;???????????%令A(yù)=D-L-U計(jì)算矩陣L
U=triu(-A)+D;??????????%令A(yù)=D-L-U計(jì)算矩陣U
M=inv(D)*(L+U);???????%計(jì)算迭代矩陣
g=inv(D)*b;????????????%計(jì)算迭代格式中的常數(shù)項(xiàng)
%
評(píng)論
共有 條評(píng)論