24小时热门版块排行榜    

查看: 1659  |  回复: 1

米尔蜜儿

铜虫 (初入文坛)

[求助] 求解三对角线矩阵方程组的追赶法matlab源程序 已有1人参与

Ax=b
A 系数矩阵 主对角线元素B1...Bj.....Bn  上三角对角线C1....Cj.....Cn-1
下三角对角线A2到An
x 为xi1.....xij.... xi,n
b  为D1....Dj....Dn

并且讲解一下

[ Last edited by 米尔蜜儿 on 2011-4-24 at 16:39 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

2011301157

金虫 (正式写手)

【答案】应助回帖

function x=Trid(a,b,c,d)
%  追赶法求解三对角的线性方程组   Ax=d
%  b为主对角线元素,a,c分别为次对角线元素,d为右端项
%  A=[ b1 c1
%      a2 b2 c2
%         ......
%            a_(n-1) b_(n-1) c_(n-1)
%                    a_(n)    b_(n) ]
%  b=[b1...b_(n)]
%  a=[0 a2...a_(n)]
%  c=[c1...c_(n-1)]

n=length(b);

u(1)=b(1);
for i=2:n
    l(i)=a(i)/u(i-1);
    u(i)=b(i)-l(i)*c(i-1);
end

y(1)=d(1);
for i=2:n
    y(i)=d(i)-l(i)*y(i-1);
end

x(n)=y(n)/u(n);
for i=n-1:-1:1
    x(i)=(y(i)-c(i)*x(i+1))/u(i);
end
2楼2015-08-10 14:48:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 米尔蜜儿 的主题更新
信息提示
请填处理意见