24小时热门版块排行榜

返回列表

米尔蜜儿

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 75.2
散金: 3
帖子: 47
在线: 12.2小时
虫号: 1232678
注册: 2011-03-14
专业: 分离过程

[求助] 求助详细讲解下述程序三对角矩阵方程追赶法

function x=followup(A,b)
%采用追赶法求线性方程组Ax=b的解
%线性方程组的系数矩阵：A
%线性方程组中的常数向量:b
%线性方程组的解：x
n=rank(A);
for(i=1:n)
if(A(i,i)==0)
disp('Error:对角所有元素为0!');
return;
end
end;
d=ones(n,1);
a=ones(n-1,1);
c=ones(n-1);
for(i=1:n-1)
a(i,1)=A(i+1,i);
c(i,1)=A(i,i+1);
d(i,1)=A(i,i);
end
d(n,1)=A(n,n);
%求解Ly=b的解y,解保存在b中
for(i=2:n)
d(i,1)=d(i,1)-(a(i-1,1)/d(i-1,1))*c(i-1,1);
b(i,1)=b(i,1)-(a(i-1,1)/d(i-1,1))*b(i-1,1);
end
%求解Ux=y的解x
x(n,1)=b(n,1)/d(n,1);
for(i=(n-1):-1:1)
x(i,1)=(b(i,1)-c(i,1)*x(i+1,1))/d(i,1);
end

回复此楼

» 猜你喜欢

1楼 2011-04-24 17:05:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主米尔蜜儿的主题更新

返回列表

24小时热门版块排行榜

米尔蜜儿

[求助] 求助详细讲解下述程序 三对角矩阵方程追赶法

» 猜你喜欢

[求助] 求助详细讲解下述程序三对角矩阵方程追赶法