Lesen Sie auch:
- [Matlab] Gauss-Seidel et Cholesky
- Intranet et extranet - Themen
- Matlab funktion definieren und plotten - Forum - Programmierung
- Classe et instance d'objet - Themen
- Et zeichen mac - Tipps - macOS
- Back office et front office - Themen
1 Antwort
Al-Jabr
Hallo,
versuche dies:
function B=cholesky(A)
[n,m]=size(A);
if m ~=n
Error('A ist keine quadratische Matrix')
end
B=zeros(size(A));
for j=1:n
S(j,j)=A(j,j)
for k=1:j-1
S(j,j)=S(j,j)-B(j,k)^2;
end
B(j,j)=sqrt(S(j,j));
for i=j+1:n
B(i,j)=A(i,j)
for k=1:j-1
B(i,j)=A(i,j)-B(i,k)*B(j,k);
end
B(i,j)=B(i,j)/B(j,j)
end
end
versuche dies:
function B=cholesky(A)
[n,m]=size(A);
if m ~=n
Error('A ist keine quadratische Matrix')
end
B=zeros(size(A));
for j=1:n
S(j,j)=A(j,j)
for k=1:j-1
S(j,j)=S(j,j)-B(j,k)^2;
end
B(j,j)=sqrt(S(j,j));
for i=j+1:n
B(i,j)=A(i,j)
for k=1:j-1
B(i,j)=A(i,j)-B(i,k)*B(j,k);
end
B(i,j)=B(i,j)/B(j,j)
end
end
Melden
thoerdim
vielen Dank