'프로그래밍/SciLab'에 해당되는 글 1건

  1. 2012.12.04 SciLab 으로 행렬의 LU분해 구하기

-->A = [1, 2, 3; 4, 5, 6]; C = A'*A;
 
-->[L, U, E] = lu(C)
 E  =
 
    0.    0.    1. 
    1.    0.    0. 
    0.    1.    0. 
 U  =
 
    27.    36.          45.       
    0.   - 0.6666667  - 1.3333333 
    0.     0.           0.        
 L  =
 
    1.           0.     0. 
    0.6296296    1.     0. 
    0.8148148    0.5    1. 
 
-->L*U
 ans  =
 
    27.    36.    45. 
    17.    22.    27. 
    22.    29.    36. 
 
-->E*C
 ans  =
 
    27.    36.    45. 
    17.    22.    27. 
    22.    29.    36. 
 
-->E'*L*U
 ans  =
 
    17.    22.    27. 
    22.    29.    36. 
    27.    36.    45. 
 
-->C
 C  =
 
    17.    22.    27. 
    22.    29.    36. 
    27.    36.    45. 
 
-->E'*L*U - C
 ans  =

   1.0D-14 *
 
  - 0.3552714    0.    0. 
    0.           0.    0. 
    0.           0.    0. 
 
-->E'*L*U - C == zeros(2, 3)
 ans  =
 
  F 
 
-->E'*L*U - C == zeros(3, 3)
 ans  =
 
  F T T 
  T T T 
  T T T 
 

-->E'*L*U == C
 ans  =
 
  F T T 
  T T T 
  T T T 
 

Posted by Scripter
,