0 Daumen
329 Aufrufe

Hi. Ich habe keine Aufgabe, sondern eine Frage bzgl. einer Umformung. Es geht um folgende Gleichung:

(JT J + λE)-1 JT  =   JT (J JT + λE)-1

J ist eine m x n (Jacobi-)Matrix, JT die Transponierte. Also ist JT J ist eine n x n symmetrische Matrix und E die n x n Einheitsmatrix. λ ist ein kleiner Skalar > 0. Beachte das rechts in der Klammer J JT steht. Es ist vorgegeben, dass (JT J + λ2 E) regulär ist, darum existiert die Inverse.


Problem: Über welche Gedankengänge/Rechenregeln gelangt man von der linken Seite auf die rechte Seite der Gleichung? Die Literatur sagt nur "Es ist einfach zu zeigen". Ich schätze ich übersehe irgendetwas offensichtliches, kann mir jemand auf die Sprünge helfen? Der Kontext ist die Methode der kleinsten gedämpften Quadrate, aber das ist vermutlich egal.

Avatar von

Ich weiß nicht, ob ein Ausmultiplizieren erlaubt wäre. Die Inverse macht mich stutzig. Kann man etwa das JT von rechts in die Klammer ziehen?

Falls ja, könnte man λ2E⋅JT anhand von AE=EA umstellen. Dann hätten beide Summanden JT auf der linken Seite und man könnte es von links ausklammern.

Falls das die Lösung ist bitte ich um Bestätigung mit einem Hinweis, warum das erlaubt ist, so banal dieser auch sein mag.

Vielen Dank :)

Mit ist nicht ganz klar, was Du mit Hineinziehen in die Inverse meinst, aber es hört sich für mich nicht schlüssig an. Alternativ siehe meine Lösung

Damit meine ich, ob evtl. folgendes geht:

(JT J + λ2 E)-1 JT = (JT J JT + λ2 E JT)-1

Ich schau mir deine Lösung gerade an, hab sie übersehen.

Es wäre besser die Einheitsmatrizen hier zu unterscheiden. Auf der linken Seite \(E_n\), die nxn Einheitsmatrix, auf der rechten dagegen \(E_m\), die mxm-Einheitsmatrix. Das hilft auch bei den Überlegungen, welche Umformungen überhaupt möglich sind.

1 Antwort

0 Daumen
 
Beste Antwort

Ich würde äquivalent umformen. Zur Abkürzung schreibe ich q für lambda^2:

$$(J^TJ+qE)^{-1}J^T=J^T(JJ^T+qE)^{-1} \iff (J^TJ+qE)^{-1}J^T(JJ^T+qE)=J^T\\ \quad \iff J^T(JJ^T+qE)=(J^TJ+qE)J^T \iff J^TJJ^T+qJ^T=J^TJJ^T+qJ^T$$

Avatar von 14 k

Vielen Dank, für die sehr schlüssige Antwort, nun bin ich beruhigt. Äquivalenzumformung muss ich definitiv meinem mathematischen Werkzeugkasten hinzufügen.

Ein anderes Problem?

Stell deine Frage

Willkommen bei der Mathelounge! Stell deine Frage einfach und kostenlos

x
Made by a lovely community