|
|
LGS lösen mit Messunsicherheit |
|
Weihnachtsmann |
Forum-Newbie
|
|
Beiträge: 2
|
|
|
|
Anmeldedatum: 31.07.12
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 31.07.2012, 23:26
Titel: LGS lösen mit Messunsicherheit
|
|
Guten Abend,
ich habe folgendes Problem:
Im Rahmen eines Physik Protokols möchte ich ein LGS Lösen. Sagen wir der Form
Dabei ist x ein unbekannter Vektor der Dimension n, A eine (nicht diagonale) nxn-Matrix und B ein bekannter Vektor der Dimension n.
Ich möchte den Vektor x berechnen. Soweit ein Standardbeispiel. Nur hat die Matrix B messunsicherheiten statistischer Art, die als gaußverteilt angenommen werden.
Wie das Prinzipiell geht ist mir klar, nach der Fehlerbehafteten Größe Ableiten, mit Fehler multiplizieren und quadratisch summieren.
Wie kann ich Fortpflanzung der Fehler in Matlab berechnen?
Gibts da ne eingebaute Funktion oder irgendwas vorgefertigtes? Das dauert ja sonst Jahre bei einer großen Matrix.
Vielen Dank und viele Grüße.
|
|
|
|
|
MaFam |
Forum-Meister
|
|
Beiträge: 799
|
|
|
|
Anmeldedatum: 02.05.12
|
|
|
|
Wohnort: ---
|
|
|
|
Version: R2009b
|
|
|
|
|
|
Verfasst am: 01.08.2012, 08:58
Titel:
|
|
Hallo,
ja, sowas gibt es. Das ist die Kondition der Matrix und diese wird mit cond(A) berechnet. Dabei gilt die Gleichung:
Grüße, Marc
|
|
|
Weihnachtsmann |
Themenstarter
Forum-Newbie
|
|
Beiträge: 2
|
|
|
|
Anmeldedatum: 31.07.12
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 01.08.2012, 12:58
Titel: Kondition auf Einträge anwenden
|
|
Vielen Dank für den Beitrag bezüglich der Kondition!
Ich könnte also jetzt beispielsweise mit einer euklidischen Norm meinen Fehler ausrechnen. Dann habe ich aber einen Fehler für die Norm von x, nicht aber für jeden Eintrag x_i von x. Ginge das auch? Kann ich also auch jeden Eintrag mit der Kondition der Matrix A multiplizieren?
Viele Grüße.
|
|
|
MaFam |
Forum-Meister
|
|
Beiträge: 799
|
|
|
|
Anmeldedatum: 02.05.12
|
|
|
|
Wohnort: ---
|
|
|
|
Version: R2009b
|
|
|
|
|
|
Verfasst am: 01.08.2012, 13:48
Titel:
|
|
Das ginge auch. Ersetze einfach x durch eine Komponente von x. Dabei sollte dir bewusst sein, dass eine Störung in b sich auf alle Komponenten von x auswirkt, es sei denn, A wäre eine Diagonalmatrix.
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
| RSS
Hosted by:
Copyright © 2007 - 2024
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|