matheraum.de
Raum für Mathematik
Offene Informations- und Nachhilfegemeinschaft

Für Schüler, Studenten, Lehrer, Mathematik-Interessierte.
Hallo Gast!einloggen | registrieren ]
Startseite · Forum · Wissen · Kurse · Mitglieder · Team · Impressum
Forenbaum
^ Forenbaum
Status Mathe
  Status Schulmathe
    Status Primarstufe
    Status Mathe Klassen 5-7
    Status Mathe Klassen 8-10
    Status Oberstufenmathe
    Status Mathe-Wettbewerbe
    Status Sonstiges
  Status Hochschulmathe
    Status Uni-Analysis
    Status Uni-Lin. Algebra
    Status Algebra+Zahlentheo.
    Status Diskrete Mathematik
    Status Fachdidaktik
    Status Finanz+Versicherung
    Status Logik+Mengenlehre
    Status Numerik
    Status Uni-Stochastik
    Status Topologie+Geometrie
    Status Uni-Sonstiges
  Status Mathe-Vorkurse
    Status Organisatorisches
    Status Schule
    Status Universität
  Status Mathe-Software
    Status Derive
    Status DynaGeo
    Status FunkyPlot
    Status GeoGebra
    Status LaTeX
    Status Maple
    Status MathCad
    Status Mathematica
    Status Matlab
    Status Maxima
    Status MuPad
    Status Taschenrechner

Gezeigt werden alle Foren bis zur Tiefe 2

Navigation
 Startseite...
 Neuerdings beta neu
 Forum...
 vorwissen...
 vorkurse...
 Werkzeuge...
 Nachhilfevermittlung beta...
 Online-Spiele beta
 Suchen
 Verein...
 Impressum
Das Projekt
Server und Internetanbindung werden durch Spenden finanziert.
Organisiert wird das Projekt von unserem Koordinatorenteam.
Hunderte Mitglieder helfen ehrenamtlich in unseren moderierten Foren.
Anbieter der Seite ist der gemeinnützige Verein "Vorhilfe.de e.V.".
Partnerseiten
Dt. Schulen im Ausland: Mathe-Seiten:Weitere Fächer:

Open Source FunktionenplotterFunkyPlot: Kostenloser und quelloffener Funktionenplotter für Linux und andere Betriebssysteme
StartseiteMatheForenLineare Algebra - MatrizenHessenberg-Matrix-Transf.
Foren für weitere Studienfächer findest Du auf www.vorhilfe.de z.B. Astronomie • Medizin • Elektrotechnik • Maschinenbau • Bauingenieurwesen • Jura • Psychologie • Geowissenschaften
Forum "Lineare Algebra - Matrizen" - Hessenberg-Matrix-Transf.
Hessenberg-Matrix-Transf. < Matrizen < Lineare Algebra < Hochschule < Mathe < Vorhilfe
Ansicht: [ geschachtelt ] | ^ Forum "Lineare Algebra - Matrizen"  | ^^ Alle Foren  | ^ Forenbaum  | Materialien

Hessenberg-Matrix-Transf.: Frage (überfällig)
Status: (Frage) überfällig Status 
Datum: 17:15 Mi 20.06.2007
Autor: Burdy

Aufgabe
Entwickeln sie einen Algorithmus, der eine Matrix mittels Ähnlichkeitstransformation durch Givensrotation auf Hessenbergform bringt.

Ich soll einen Algorithmus aufschreiben, der eine Matrix mittels Givensrotation auf Hessenberg-Form [mm] (h_{ij}=0 [/mm]  für alle i>j+1) bringt und gleichzeitig eine Ähnlichkeitstransformation ist, damit die Eigenwerte erhalten bleiben.
Für die gewöhnliche Givensrotation auf obere Dreiecksgestalt ohne Ähnlichkeitstransformation  hab ich das schon mal aufgeschrieben:

H = A
for i = 1:n-1
    for j = i+1:n        
        a = H(i,i);
        b = H(j,i);
        if ( b ~= 0)
            c = a / [mm] sqrt(a^2 [/mm] + [mm] b^2); [/mm]
            s = b / [mm] sqrt(a^2 [/mm] + [mm] b^2); [/mm]
            G = eye(n); % Einheitsmatrix vom Rang n
            G(i,i) = c;
            G(j,j) = c;
            G(i,j) = s;
            G(j,i) = -s;            
            H = G*H;
        end          
    end
end

Wenn ich da jetzt einfach zur Erzeugung einer Hessenberg-Matrix "j = i+1" auf "j = i+2" ändere, damit nicht zuerst der erste Eintrag unter der Diagonale eliminiert wird, sondern der 2. geht das zwar korrekt für die erste Spalte, aber bei der zweiten Spalte werden dann die in der ersten Spalte erzeugten Nullen teilweise wieder zerstört.
Eine Idee, wie man das umgehen kann wäre super.

Und wenn das gelöst ist hab ich noch das Problem, das, damit das ganze eine Ähnlichkeitstransformation ist, ich ja eigentlich jedes mal [mm] H=G^{T}*H*G [/mm] rechnen müsste und nicht H = G*H. Nur so, in der Form, wie das da jetzt steht, sind dann gleich alle Nullen weg.
Was muss ich denn dafür ändern?

Also irgendwie steh ich da grade aufm Schlauch und mir fällt nichts ein, ich wär für jede Hilfe dankbar.

Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.

        
Bezug
Hessenberg-Matrix-Transf.: Fälligkeit abgelaufen
Status: (Mitteilung) Reaktion unnötig Status 
Datum: 17:20 Sa 23.06.2007
Autor: matux

$MATUXTEXT(ueberfaellige_frage)
Bezug
Ansicht: [ geschachtelt ] | ^ Forum "Lineare Algebra - Matrizen"  | ^^ Alle Foren  | ^ Forenbaum  | Materialien


^ Seitenanfang ^
www.matheraum.de
[ Startseite | Forum | Wissen | Kurse | Mitglieder | Team | Impressum ]