登录论坛

查看完整版本 : [MATLAB基础] 脚本文件中的编程问题


shaoshuaithe
2011-10-24, 11:38
%这个程序是利用高斯消元法求解线性方程
%这个程序用于求解线性方程组,A的前三列是系数,B是根矩阵初始值。最后一列是方程等号右边的数
%第一个嵌套循环是将A前三列转换成上三角型。
%第二循环是求根,然后赋值给B矩阵
A=[4 -1 1 6;-1 4.25 2.75 -0.5;1 2.75 3.5 1.25];
B=[0 0 0];
for i=1:2
for j=i:2
m=A(j+1,i)/A(i,i);
A(j+1,:)=A(j+1,:)-A(i,:)*m;
end
end
A;
for k=3:1
for h=3:k
A(k,4)=A(k,4)-A(k,h)*B(h);
B(k)=A(k,4)/A(k,k);
end
end
B
这是我在脚本文件中编写的程序,问题是:矩阵B的值一直没变,我怀疑第二个嵌套循环(从for k=3:1开始)没有运行