Red de conocimiento informático - Aprendizaje de código fuente - El análisis numérico de la descomposición de LU requiere programación en Matlab para encontrar L y U detA

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

detA =

0

gt;gt;