WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Addition nicht korrekt

 

Schokohoernle
Forum-Newbie

Forum-Newbie


Beiträge: 3
Anmeldedatum: 13.01.11
Wohnort: AC
Version: R2010b
     Beitrag Verfasst am: 13.01.2011, 16:46     Titel: Addition nicht korrekt
  Antworten mit Zitat      
Hallo!

Ich bin da auf ein Problem in MATLAB gestoßen.
Mein Code sieht folgender Maßen aus:

Code:

a = 4.982435678234596e+05
b = 5.237647823678964e+04
c = -b
d = -a
 


Nun sollte man ja meinen dass a+b+c+d = 0 ist, ist es aber nicht!

Code:

>> a+b+c+d

ans =

 -5.8208e-011
 


Allerdings ist wenn man b und c verstauscht:
Code:

>> a+c+b+d

ans =

     0
 


und zudem:

Code:

>> (a+b)+(c+d)

ans =

     0
 


Meine Vermutung ist, dass es daran liegt, dass die beiden Zahlen a und b eine Größenordnung unterschied haben. und Beim Addieren der beiden unterschiedlichen Größenordnungen Ungenauigkeiten auftreten.

Kennt einer das Problem bzw. weiß einer wie ich das lösen kann?

Vielen Dank!

Niko
Private Nachricht senden Benutzer-Profile anzeigen


Maddy
Ehrenmitglied

Ehrenmitglied



Beiträge: 494
Anmeldedatum: 02.10.08
Wohnort: Greifswald
Version: ---
     Beitrag Verfasst am: 13.01.2011, 16:54     Titel:
  Antworten mit Zitat      
Das sind halt die Tücken der numerischen Genauigkeit. Wenn du es überschlägst, ist die Rechnung auf 16 Größenordnungen genau.
Also vollkommen im numerischen Limit.

Was hast du denn vor mit dem Ergebnis bzw. warum benötigst du diese Rechnung?
_________________

>> why
The computer did it.
Private Nachricht senden Benutzer-Profile anzeigen
 
Schokohoernle
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 3
Anmeldedatum: 13.01.11
Wohnort: AC
Version: R2010b
     Beitrag Verfasst am: 14.01.2011, 13:23     Titel:
  Antworten mit Zitat      
Das ist für die Berechnung von Faserverbundwerkstoffen.
Da ist bei einem symmetrischen Schichtaufbau die Koppelmatrix genau 0. Die Koppelmatrix ist aber eine Summe und ergibt in manchen Fällen wegen o. g. Fehler eben nicht genau 0.
Das ist nicht nur unschön sondern führt im weiteren Berechnungsverlauf auch zu kleinen Fehlern und genau diese kleinen Fehler würde ich gerne beseitigen.

MfG

Niko
Private Nachricht senden Benutzer-Profile anzeigen
 
Jan S
Moderator

Moderator


Beiträge: 11.057
Anmeldedatum: 08.07.10
Wohnort: Heidelberg
Version: 2009a, 2016b
     Beitrag Verfasst am: 14.01.2011, 13:45     Titel: Re: Addition nicht korrekt
  Antworten mit Zitat      
Hallo Schokohoernle,

Herzlich willkommen in der Welt der numerischen Mathematik.
Rundungsfehler sind allgegenwärtig, wenn man mit begrenzter Genauigkeit rechnet. Versuche dies:
Code:
1e9 + 1e-9 - 1e9

Es kommt beim Summieren also auf die Reihenfolge an, weil intern die 64-Bit-Repräsentation der DOUBLEs nicht ausreichen, um einen Unterschied zwischen 1e9 und 1e9+1e-9 zu speichern.

Gruß, Jan
Private Nachricht senden Benutzer-Profile anzeigen
 
Andreas Goser
Forum-Meister

Forum-Meister


Beiträge: 3.654
Anmeldedatum: 04.12.08
Wohnort: Ismaning
Version: 1.0
     Beitrag Verfasst am: 14.01.2011, 13:49     Titel:
  Antworten mit Zitat      
Ich kann auch noch diesen Link empfehlen:

http://www.mathworks.com/support/so.....1-16FOQ/?solution=1-16FOQ

Andreas
Private Nachricht senden Benutzer-Profile anzeigen E-Mail senden
 
Maddy
Ehrenmitglied

Ehrenmitglied



Beiträge: 494
Anmeldedatum: 02.10.08
Wohnort: Greifswald
Version: ---
     Beitrag Verfasst am: 14.01.2011, 14:41     Titel:
  Antworten mit Zitat      
Schokohoernle hat Folgendes geschrieben:
Das ist für die Berechnung von Faserverbundwerkstoffen.
Da ist bei einem symmetrischen Schichtaufbau die Koppelmatrix genau 0. Die Koppelmatrix ist aber eine Summe und ergibt in manchen Fällen wegen o. g. Fehler eben nicht genau 0.
Das ist nicht nur unschön sondern führt im weiteren Berechnungsverlauf auch zu kleinen Fehlern und genau diese kleinen Fehler würde ich gerne beseitigen.

MfG

Niko


Was für Genauigkeiten bzw. Größenordnungen hast du denn sonst so in der Koppelmatrix? Im günstigsten Fall würde es wohl reichen einfach ein bisschen zu runden. Wink
_________________

>> why
The computer did it.
Private Nachricht senden Benutzer-Profile anzeigen
 
Schokohoernle
Themenstarter

Forum-Newbie

Forum-Newbie


Beiträge: 3
Anmeldedatum: 13.01.11
Wohnort: AC
Version: R2010b
     Beitrag Verfasst am: 18.01.2011, 11:21     Titel:
  Antworten mit Zitat      
Ja das habe ich mir auch schon überlegt. Es reicht wenn man die Letzte Kommastelle bei den Zahlen mit niedrigerer Größenordnung weglässt, bzw. rundet.
Was ich mal ausprobieren muss ist folgendes:
Ich speichere die einzelnen Summanden in einer 9xn Matrix (die Koppelmatrix ist eine 3x3 Matrix), sortiere dann jede Zeile dieser 9xn Matrix nach der Größe und summiere sie dann auf. Müsste ja eigentlich funktionieren...

MfG

Niko
Private Nachricht senden Benutzer-Profile anzeigen
 
Jan S
Moderator

Moderator


Beiträge: 11.057
Anmeldedatum: 08.07.10
Wohnort: Heidelberg
Version: 2009a, 2016b
     Beitrag Verfasst am: 18.01.2011, 15:56     Titel:
  Antworten mit Zitat      
Hallo Schokohoernle,

Zitat:
Ich speichere die einzelnen Summanden in einer 9xn Matrix (die Koppelmatrix ist eine 3x3 Matrix), sortiere dann jede Zeile dieser 9xn Matrix nach der Größe und summiere sie dann auf. Müsste ja eigentlich funktionieren...

Nein, leider gibt es kein allgemeines Konzept, um die Berechnung einer Summe stabil zu machen. Ich habe in der FEX die Funktion XSum veröffentlicht (http://www.mathworks.com/matlabcentral/fileexchange/26800-xsum), die intern mit 30 oder 45 Ziffern arbeitet. Damit kann man zwar die Summen schon viel genauer berechnen, aber trotzdem gibt es für jede Anzahl von Stellen (oder jede Stelle, auf die man runden möchte) eine Folge von Zahlen, die einen Verlust von Bits provoziert.

Ich würde eher normal summieren und das Endergebnis mit einer gewissen Entfernung von Null akzeptieren, z.B. Anzahl der Summanden * MAX(Summanden) * EPS.

Gruß, Jan
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

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 | goMatlab RSS Button 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.