Matroids Matheplanet Forum Index
Moderiert von Fabi Dune ligning
Mathematik » Lineare Algebra » Algorithmus für Gauss mittels Cholesky
Druckversion
Druckversion
Autor
Universität/Hochschule J Algorithmus für Gauss mittels Cholesky
Ehemaliges_Mitglied
Zum letzten BeitragZum nächsten BeitragZum vorigen BeitragZum erstem Beitrag  Themenstart: 2014-10-29


fed-Code einblenden
MATLAB
function [x,y] = choleskysol(L,b)
%CHOLESKYSOL Loest LGS A*x=b mittels Cholesky-Zerlegung.
%   L: invertierbare untere nxn-Dreiecksmatrix.
 
n=length(L(:,1));
y=zeros(n,1);
x=zeros(n,1);
T=L';
 
%VORWAERTSEINSETZEN
for i=1:n
    summe=0;
    for k=1:i-1
        summe=summe+L(i,k)*y(k);
    end;
    y(i)=1/L(i,i)*(b(i)-summe);
end;
 
%RUECKWAERTSEINSETZEN
for j=n:-1:1
    summe2=0;
    for q=j+1:n
        summe2=summe2+T(j,q)*x(q);
    end;
    x(j)=1/T(j,j)*(y(j)-summe2);
end;
end

Dieser Code funktioniert nicht. Ich habe mir als Test mal das y mit ausgeben lassen. Daran scheitert es schon.
Kann mir jemand sagen, wo das Problem liegt?

Danke! Freundliche Grüße,
Pauline



Eine Notiz zu diese Forumbeitrag schreiben Notiz   Profil  Quote  Link auf diesen Beitrag Link
Ehemaliges_Mitglied hat selbst das Ok-Häkchen gesetzt.
Neues Thema [Neues Thema]  Druckversion [Druckversion]

 


Wechsel in ein anderes Forum:
 Suchen    
 
All logos and trademarks in this site are property of their respective owner. The comments are property of their posters, all the rest © 2001-2020 by Matroids Matheplanet
This web site was originally made with PHP-Nuke, a former web portal system written in PHP that seems no longer to be maintained nor supported. PHP-Nuke is Free Software released under the GNU/GPL license.
Ich distanziere mich von rechtswidrigen oder anstößigen Inhalten, die sich trotz aufmerksamer Prüfung hinter hier verwendeten Links verbergen mögen.
Lesen Sie die Nutzungsbedingungen, die Distanzierung, die Datenschutzerklärung und das Impressum.
[Seitenanfang]