24小时热门版块排行榜    

查看: 161  |  回复: 0

米尔蜜儿

铜虫 (初入文坛)

[求助] 求助详细讲解以下程序

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
回复此楼

» 猜你喜欢

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

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 米尔蜜儿 的主题更新
信息提示
请填处理意见