矩阵迭代：这段代码的功能是什么？

Question

它是解决矩阵的gauss seidel迭代方法的一部分，但我不确定这段代码是做什么的。矩阵“A”是预定义的。

[columns,rows] = size(A); 
L = zeros(size(A)); 
U = zeros(size(A)); 

for i = 1:columns 
for j = 1:i 
    L(i,j) = A(i,j); 
end 
for j = i+1:columns 
    U(i,j) = A(i,j); 
end 
end 

Answer 1

该函数只是发现矩阵的lower and upper triangular portions。你得到那个“未定义”的错误，因为你还没有定义A。限定A，然后运行该代码：

A = rand(6,6); % For example 
% Rest of your code goes here.... 

L和U将包含升奥尔和ù矩阵A的PPER部分。

小注：您正在使用size错误。第一个输出参数是行，第二个输出参数是列。交换变量：

[rows, columns] = size(A);