Fattorizzazione di cholesky
Ciao a tutti:)!!
Ho dei grossi problemi con la fattorizzazione di cholesky.
Non riesco a trovare l'algoritmo per matlab.
Qualcuno può aiutarmi?..grazie:)...
Ho dei grossi problemi con la fattorizzazione di cholesky.
Non riesco a trovare l'algoritmo per matlab.
Qualcuno può aiutarmi?..grazie:)...
Risposte
function [L] = cholesky(A) %fattorizza A simm e def pos in A=L*L'
n=size(A,1);
L=ones(n,n);
for k=1:n
for j=1:k-1 %calcola tutti gli elementi della k-esima riga
somma=0; %che servono a completare L_k
for i=1:j-1
somma=somma+L(j,i)*L(k,i);
end
L(k,j)=(A(j,k)-somma)/L(j,j);
end
somma2=0;
for j=1:k-1
somma2=somma2+(L(k,j)^2);
end
L(k,k)=sqrt(A(k,k)-somma2); %ultimo elemento della riga,quello diagonale non ancora calcolato
end
end.
a me funziona .
dove studi?
n=size(A,1);
L=ones(n,n);
for k=1:n
for j=1:k-1 %calcola tutti gli elementi della k-esima riga
somma=0; %che servono a completare L_k
for i=1:j-1
somma=somma+L(j,i)*L(k,i);
end
L(k,j)=(A(j,k)-somma)/L(j,j);
end
somma2=0;
for j=1:k-1
somma2=somma2+(L(k,j)^2);
end
L(k,k)=sqrt(A(k,k)-somma2); %ultimo elemento della riga,quello diagonale non ancora calcolato
end
end.
a me funziona .
dove studi?
grazie mille!!!ora provo subito!!:)
Studio a bologna e ho finito il secondo anno di matematica.
E mi manca solo sto esame del cavolo..mi porta via troppo tempo..
Studio a bologna e ho finito il secondo anno di matematica.
E mi manca solo sto esame del cavolo..mi porta via troppo tempo..