Matlab中的矩阵索引和除法

Question

我想编码高斯消除，这是我正在使用的代码段。矩阵已经是三角形的，而是有这表达了一个错误：

x(n) = b(n)/A(n, n); 

下面是完整的一段代码：

function [x] = gauss(A, b) 
    n = size(A); 
    for k=1 : n-1 
     for i=k+1 : n 
      m = A(i, k)/A(k, k); 
      for j=1 : n 
       A(i, j) = A(i, j) - m * A(k, j); 
      end 
      b(i) = b(i) - m * b(k); 
     end 
    end 
    x(n) = b(n)/A(n, n); 
    disp(double(b(n))); 
    for k=n-1 : -1 : 1 
     s=0; 
     for i=k+1 : n 
      s = s + A(k, i) * x(i); 
     end 
     x(k) = (b(k) - s)/A(k, k); 
    end 
    disp(x); 
end 

感谢您的帮助。顺便说一句IM在Matlab新手...

编辑：添加一些更多的信息

林调用这样此功能： A = [6 3 2; 9 -1 4; 10 5 3] B = [12 37 21]

sol = gauss(A, b); 

Answer 1

此行

n = size(A); 

将有一个矢量结果作为A是一个矩阵。看来，你期待此行

x(n) = b(n)/A(n, n); 

表现得像一个标量划分，但n矢量你实际上是企图分裂不同维度的矩阵。根据您的示例代码检查这些语句的输出

n = size(A); 
A(n,n) 
b(n) 

并且看到您没有处理标量。如果你想n等于的A行数标尝试

[n n1] = size(A); 

更换

n = size(A); 

现在n应该真正代表的A