Invers Cosinus Probleme

[Weazy]

Hi,

ich habe da noch so ne blöde frage:


wenn ich in folgender Funktion als z eine real Zahl einsetze, ist das Resultat "NAN". Wenn ich hingegen -1 einsetze ist das Ergebniss korrekterweise 180.

radtodeg(arccos((z)));

Wisst ihr woran das liegt?

[dizzy]

Was für Werte hast du denn eingegeben? Definiert ist der Arcuscosinus nur über -1..1. Alles darüber und darunter sollte jedoch eher eine Exception statt NaN liefern :gruebel:

[himitsu]

Es kommt auf die Fehlerbehandlung an, wenn diese abgestellt ist, wird korrekter Weise NaN (Not a Number) zurückgegeben, wenn sich ein Fehler ereignet hat.

Und wwie dizzy schon sagte, was gibts du für Werte dort ein?
Diese dürfen ja schließlich nur in einem bestimmten Bereich liegen.

-1 <= Z <= 1

[Weazy]

Ich hab es gefunden! Ich runde den eingegbenen Wert einfach mit Roud()!

Vielen Dank trotzdem!

[dizzy]

Äääähm, damit sorgst du dafür dass arccos() ganze 3 mögliche Ergebnisse liefert - und das ist gewiss nicht akkurat :roll:. Eine wirkliche Lösung des Problems kann das imho nicht sein ;)

[Weazy]

Wie denn das?

[3_of_8]

:shock:

Ganz einfach: Dann hast du entweder:
-1,0, oder 1. Und damit auch nur drei verschiedene Ergebnisse. Das bringts auch voll.

[Weazy]

Wie kann man dann z.b auf 5 Stellen genau runden?

[himitsu]

Ganz genau Runden geht eh nicht ... liegt halt an der internen Struktur der Typen (Single, Double...).

Aber du kannst es so machen:
X := Runden(X * 10^5) / 10^5;

Delphi-Quellcode:

X := Round(X * 100000) / 100000;

Und hast du mal nachgesehn, was du für Werte in z übergibst?
Notfalls mal vorher testen:

Delphi-Quellcode:

If (z < -1) or (z > 1) Then 
  MessageBox('Fehler: ' + FloatToStr(z), ...

[3_of_8]

Ich hab das ganze mal getestet:

Delphi-Quellcode:

with memo1.lines do
begin
add(floattostr(radtodeg(arccos(-1))));
add(floattostr(radtodeg(arccos(-0.75))));
add(floattostr(radtodeg(arccos(0))));
add(floattostr(radtodeg(arccos(0.75))));
add(floattostr(radtodeg(arccos(1))));
add(floattostr(radtodeg(arccos(1.0001))));
end;

Alle Befehle haben ein korrektes Ergebnis gebracht, der letzte Befehl (quod erat exspectandum) eine Exception (Ungültige Gleitkommaoperation).

[himitsu]

@3_of_8:
Versuch es mal mit deaktivierten FPU-Exceptions :zwinker:

Delphi-Quellcode:

var
  Saved8087CW: Word;
begin
  Saved8087CW := Default8087CW;
  Set8087CW($133f);
  with memo1.lines do
  begin
    add(floattostr(radtodeg(arccos(-1.0001)))); 
    add(floattostr(radtodeg(arccos(-1))));
    add(floattostr(radtodeg(arccos(-0.75))));
    add(floattostr(radtodeg(arccos(0))));
    add(floattostr(radtodeg(arccos(0.75))));
    add(floattostr(radtodeg(arccos(1))));
    add(floattostr(radtodeg(arccos(1.0001))));
  end;
  Set8087CW(Saved8087CW);
end;

@Weazy:
Wie du siehst gibt es "nur" Probleme, wenn z außerhalb des zulässigen Bereichs liegt.

[3_of_8]

Was genau ist FPU? Und NaN heißt Not a Number, nicht Natriumnitrat, oder?

[Daniel G]

Zitat:

Zitat von 3_of_8

Was genau ist FPU?

Floating-Point-Unit. Der Teil deiner CPU, der für die Gleitkommaoperationen zuständig ist. (Grob gesagt)

[3_of_8]

Achso.

[himitsu]

Zitat:

Zitat von 3_of_8 #12

Was genau ist FPU? Und NaN heißt Not a Number, nicht Natriumnitrat, oder?

Zitat:

Zitat von himitsu #3

Es kommt auf die Fehlerbehandlung an, wenn diese abgestellt ist, wird korrekter Weise NaN (Not a Number) zurückgegeben, wenn sich ein Fehler ereignet hat.