矩陣迭代：這段代碼的功能是什麼？

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);