Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
![]() |
#1 |
初级会员
注册日期: 2011-10-15
年龄: 38
帖子: 1
声望力: 0 ![]() |
![]()
%这个程序是利用高斯消元法求解线性方程
%这个程序用于求解线性方程组,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, ![]() ![]() ![]() 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开始)没有运行 |
![]() |
![]() |