El análisis numérico de la descomposición de LU requiere programación en Matlab para encontrar L y U detA
Hay muchos programas escritos por expertos en Internet, sólo búscalos.
Encontré uno, échale un vistazo.
Programa de descomposición 1.lu
función [l, u, detA] = mylu( a )
Este programa realiza principalmente la descomposición LU de la matriz de coeficientes a , el algoritmo específico se puede encontrar en el análisis numérico.
n=size(a, 2);
u=zeros(size(a));
l= ojo(tamaño(a));
u(1,:)=a(1,:);
l(2:end,1)=a( 2:fin, 1)/a(1, 1);
para r=2: n
para j=r: n
u(r , j) =a(r,j)-l(r,1:r-1)*u(1:r-1,j);
fin
para i =r 1 :n
l(i,r)=(a(i,r)-l(i,1:r-1)*u(1:r-1,r))/ u(r , r);
fin
fin
detA=det(a);
2. >
gt; gt; A=magia(4);
gt; [L, U, detA]=mylu(A)
L =
1.0000 0 0 0
0.3125 1.0000 0 0
0.5625 0.5663 1.0000 0
0.2500 1.3012 -3.0000 1.0000
U =
p>16.0000 2.0000 3.0000 13.0000
0 10.3750 9.0625 3.9375
0 0 -0.8193 2.4578
0 0 0 0 p>
detA =
0
gt;gt;