| 查看: 807 | 回复: 0 | ||
[求助]
求助详细讲解下述程序 三对角矩阵方程追赶法
|
|
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 |
回复此楼» 猜你喜欢
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？
已经有5人回复
-
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有6人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有16人回复
-
版面费该交吗
已经有13人回复
-
面上可以超过30页吧？
已经有6人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有10人回复
-
什么是人一生最重要的？
已经有4人回复