LU分解CROUT减少

-3

这是我的LU分解CROUT方法代码：LU分解CROUT减少

function [L, U] = croutreduction(A) 
    [row,column]=size(A); 
    L=eye(row,column); 

    //A = 3x3 
    if row==3 then 
     U(1,1)=A(1,1); U(1,2)=A(1,2); U(1,3)=A(1,3); 
     L(2,1)=A(2,1)/U(1,1); L(3,1)=A(3,1)/U(1,1); 

     U(2,2)=A(2,2)-L(2,1)*U(1,2); 
     U(2,3)=A(2,3)-L(2,1)*U(1,3); 
     L(3,2)=A(3,2)/U(2,2); 

     U(3,3)=A(3,3)-L(3,2)*U(2,3);  
    end 

    //A = 4x4 
    if column==4 then 
     U(1,1)=A(1,1); U(1,2)=A(1,2); U(1,3)=A(1,3); U(1,4)=A(1,4); 
     L(2,1)=A(2,1)/U(1,1); L(3,1)=A(3,1)/U(1,1); L(4,1)=A(4,1)/U(1,1); 

     U(2,2)=A(2,2)-L(2,1)*U(1,2); 
     U(2,3)=A(2,3)-L(2,1)*U(1,3); 
     U(2,4)=A(2,4)-L(2,1)*U(1,4); 
     L(3,2)=(A(3,2)-L(3,1)*U(1,2))/U(2,2); 
     L(4,2)=(A(4,2)-L(4,1)*U(1,2))/U(2,2); 

     U(3,3)=A(3,3)-(L(3,1)*U(1,3)+L(3,2)*U(2,3)); 
     U(3,4)=A(3,4)-(L(3,1)*U(1,4)+L(3,2)*U(2,4)); 
     L(4,3)=(A(4,3)-(L(4,1)*U(1,3)+L(4,2)*U(2,3)))/U(3,3); 

     U(4,4)=A(4,4)-(L(4,1)*U(1,4)+L(4,2)*U(2,4)+L(4,3)*U(3,4)); 
    end 
endfunction

如何修改我的代码具有不同尺寸的矩阵来工作？如您所见，上面的代码仅适用于3x3和4x4矩阵。

来源

2015-04-04 Elra Ghifary

Crout缩减代码是通过谷歌avaialbe和关于第一个问题，请尝试做到这一点，并在您有任何问题和问题后，我们将更乐意帮助您 – madbitloman 2015-04-04 21:47:59

您可以检查我的新代码？ – 2015-04-06 16:33:15

如果您有N维矩阵，我建议您将其重新整形为正方形2D矩阵，然后执行反演，然后重新整形。 – madbitloman 2015-04-06 17:49:51

您应该使用for循环代替硬编码索引。在此基础上例如：http://en.wikipedia.org/wiki/Crout_matrix_decomposition我修改的Scilab代码（原代码是对C和Matlab /八度）：

function [L,U]=LUdecompCrout(A) 
    [r,c]=size(A); 
    for i=1:r 
    L(i,1)=A(i,1); 
    U(i,i)=1; 
    end 
    for j=2:r 
    U(1,j)=A(1,j)/L(1,1); 
    end 
    for i=2:r 
    for j=2:i 
     L(i,j)=A(i,j)-L(i,1:j-1)*U(1:j -1,j); 
    end 

    for j=i+1:r 
     U(i,j)=(A(i,j)-L(i,1:i-1)*U(1:i-1,j))/L(i,i); 
    end 
    end 
endfunction

Hovewer这给了不同的结果，你的代码，我没有检查哪一个是错误和其中...

来源

2015-04-06 21:26:05 Attila

@Elra Ghifary此代码适合您吗？ – Attila 2015-04-08 08:40:16

LU分解CROUT减少

回答

相关问题