seltsames Problem < Matlab < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) beantwortet  | | Datum: | 11:23 Do 25.10.2007 | | Autor: | JayKayLE |
allo,
ich habe in die Matrix InclinosMat Messdaten mit Zeitstempel eingelesen. diese sollen nun in eine Matrix CompMat kopiert werden. Da nicht zu jeden Zeitpunkt Messdaten vorhanden sind, ich aber über jeden Zeitpunkt plotten will (gibt noch andere Messwerte), will ich in die Zeilen von CompMat, zu denen in InclinosMat kein Zeitstempel existiert, ein NaN reinschreiben.
Die Zeitstempel stehen jeweils in der ersten Spalte der Matrizen.
Code:
for i=1:size(CompMat,1)
for k=1:size(InclinosMat,1)
if (CompMat(i,1) == InclinosMat(k,1))
CompMat(i,3) = InclinosMat(k,10); % Aus Spalte 10 in Spalte 3
InclinosMat(k,:) = []; % betreffende Zeile löschen, dann gehts beim nächsten Suchen schneller
break;
end;
if (k == size(InclinosMat,1)) % Kein Treffer, Abbruch
CompMat(i,3) = NaN;
break;
end;
end;
end;
Problem: Auch wenn in beiden Matrizen die gleichen Zeitstempel stehen, werden die Daten manchmal nicht kopiert - und zwar genau bei x.1 x.2 x.6 und x.7 (wobei x für die Zeit in Sekunden steht). Das heißt: Wird zB der Zeitstempel 500.1 verglichen, der in beiden Matrizen steht, werden beide if-Anweisungen einfach übersprungen! Wenn ich vor der ersten if Anweisung einen Haltepunkt setze und die werte der beiden Matrizen über die Matlab-Console abfrage, stimmen sie überein, die if Anweisung müsste true werden - wird sie aber nicht. Habe ich einen Denkfehler?
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 14:41 Do 25.10.2007 | | Autor: | JayKayLE |
Habe es selbst herausgefunden. Ich erzeuge die Zeitstempel mit
x = StartTime:0.1:EndTime;
x = x';
Matlab erzeugt dabei aber nicht wirklich die Werte, die ich haben will - aus 518,6 wird zum Beispiel immer 518,5999999999... Dasselbe bei x,1 x,2 und x,7. Alles andere passt bzw wird "genau" erzeugt. Darauf soll man erstmal kommen.
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig | | Datum: | 16:12 Do 25.10.2007 | | Autor: | rainerS |
Hallo!
> Habe es selbst herausgefunden. Ich erzeuge die Zeitstempel
> mit
>
> x = StartTime:0.1:EndTime;
> x = x';
>
> Matlab erzeugt dabei aber nicht wirklich die Werte, die ich
> haben will - aus 518,6 wird zum Beispiel immer
> 518,5999999999... Dasselbe bei x,1 x,2 und x,7. Alles
> andere passt bzw wird "genau" erzeugt. Darauf soll man
> erstmal kommen.
Ich möchte wetten, dass auch andere ungenau dargestellt werden. Du addierst ja immer 0.1, was sich nicht als als endliche Binärbruch darstellen lässt und daher im Rahmen der Rechengenauigkeit auf die nächste Binärzahl gerundet wird. Die Rundungsfehler akkumulieren sich durch die fortgesetzte Addition. Was du siehst, sind nur die Ausgaben als Dezimalzahl, die durch Runden auf die nächstgelegene Dezimalzahl im Rahmen der Ausgabegenauigkeit entsteht. Sieh ![[]](/images/popup.gif) hier.
Faustregel: bei Gleitkommazahlen nie auf Gleichheit testen.
Viele Grüße
Rainer
|
|
|
| |
|
Hallo,
ich weiß ja nicht genau, wo die Zeitstempel herkommen, aber könnte es sein, dass sie sich (vielleicht Matlabs Schuld?) um 0.00001 (oder so ähnlich) unterscheiden? Ändert sich das Verhalten, wenn du den Vergleich auf Gleichheit durch if abs(a-b)<0.001 ersetzt?
Sonst ist deine Logik natürlich richtig, auch wenn ich den sonstigen Zusammenhang nicht ganz durchschaue.
Gruß
Martin
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig | | Datum: | 15:22 Do 25.10.2007 | | Autor: | JayKayLE |
*g* Jo das wäre schon denkbar ;)
Danke für deine Hilfe!
|
|
|
|
