資源簡介
運(yùn)用matlab求解三對(duì)角方程組的追趕法原始函數(shù)
代碼片段和文件信息
function?[k]=zhuigan(abcd)
%對(duì)角線下方的元素,個(gè)數(shù)比A少一個(gè)
%?%對(duì)角線元素
%對(duì)角線上方的元素,個(gè)數(shù)比A少一個(gè)
%d為方程常數(shù)項(xiàng)
%用追趕法解三對(duì)角矩陣方程
r=size(a);
m=r(2);
r=size(b);
n=r(2);
if?size(a)~=size(c)|m~=n-1|size(b)~=size(d)
????error(‘變量不匹配檢查變量輸入情況!‘);
end
%%
%LU分解
u(1)=b(1);
for?i=2:n
????l(i-1)=a(i-1)/u(i-1);
????u(i)=b(i)-l(
評(píng)論
共有 條評(píng)論