Berechnung mit cos führt zum Fehler

CG schrieb:

>

1

 Zeitdifferenz = ( 12*acos( ( sin(h[i]) - sin(B)*sin(DK) ) / 

2

> (cos(B)*cos(DK)) ) ) / pi;

>
> inludiert ist die math.h
> für
> h[0] = -0,833
> B = 0,91583333
> DK = -0,1747629
>
> wenn ich mich nicht verrechnet habe müsste für die Zeitdiferenz
> 5,21131484 errechnet werden.
> Leider rechnet mein Controller nur Unsinn. Kann es sein das er mit den
> cos(DK) nicht zurecht kommt?

kann ich mir nicht vorstellen.
Was in solchen Sachen immer eine gute Hilfe ist, das ist wenn man den 
Monsterausdruck einfach mal in Teilausdrücke aufdröselt

1

  float cosDK = cos(DK);

2

  ....

und sich die Teilergebnisse ansieht bzw. mit dem Taschenrechner 
vergleicht.

Aber Vorsicht: Im Gegensatz zu der Normaleinstellung eines 
Taschenrechners, rechnet die C Funktion den Winkel in Radianten und 
nicht in Grad! Das ist sogar ein recht häufiger Fehler, dass dieses 
Faktum übersehen wird. Sei es im C-Programm oder sei es bei der 
Kontrolle im Taschenrechner.

Rechnest DU ihn Grad oder Bogenmaß?

CG schrieb:

> @ Karl Heinz
> Das ausrechnen in einzelnen Schritten habe ich auch schon ausprobiert
> und leider kommt auch nicht gescheites dabei raus.

Was muss ich mir unter 'nichts gescheites' vorstellen?

> Im Debuger steht unter Zeitdifferenz bei Value n. def.

Na ja.
Da ist eine komplexe Berechnung und am Ende kommt 'nicht definiert' 
raus. Ich schätze mal: Division durch 0.0

D.h. nächste Frage: wird tatsächlcih durch 0 dividiert?

Also aufteilen:

1

 float Dividend = sin(h[i]) - sin(B)*sin(DK);

2

 float Divisor = cos(B)*cos(DK);

3

4

 Zeitdifferenz = ( 12*acos( Dividend / Divisor  ) ) / pi;

und die Einzelwerte ansehen bzw. mit dem Taschenrechner verlgeichen.

Es tut mir leid, wenn ich dir nicht so helfen kann, dass ich mit dem 
Finger auf eine Stelle zeige und sage: hier ist der Fehler.
Ich kann dir nur zeigen, wie ich das angehen würde. Alle die mich hier 
kennen wissen sehr wahrscheinlich, dass ich mit der Methode "wenn etwas 
zu komplex ist, dann muss man es eben in einfachere Einheiten zerlegen" 
ganz gut fahre.
Und genau das würde ich auch hier machen. Die Berechnung ist in ihrer 
Gesamtheit komplex und ganz zum Schluss kommt Unsinn raus. Ergo zerleg 
ich die eben in Teile und sieh mir die Teile an. Wenn feststeht, welcher 
dieser Teile für den Unsinn verantwortlich ist, dann zerleg ich den 
wieder weiter. Solange bis ich bei der eigentlichen Ursache angelangt 
bin.

So würde ich das angehen. Du magst das gerne anders halten. Aber sich 
hinstellen und sagen "zum Schluss kommt n.def raus - was könnte die 
Ursache sein", bringt dich nicht weiter.

Karl Heinz schrieb:

> Da ist eine komplexe Berechnung und am Ende kommt 'nicht definiert'
> raus. Ich schätze mal: Division durch 0.0

Der acos könnte auch die Ursache sein.
Dividend / Divisor muss im Bereich -1.0 bis +1.0 sein, damit der acos 
definiert ist. Schon ein ganz klein wenig drüber .... und das Ergebnis 
ist nicht mehr definiert.

Also: Was ergibt

1

  float Dividend = sin(h[i]) - sin(B)*sin(DK);

2

  float Divisor = cos(B)*cos(DK);

3

  float Quotient = Dividend / Divisor;

4

5

  Zeitdifferenz = ( 12*acos( Quotient  ) ) / pi;

Welchen Wert hat Quotient?

PS. die math.lib hast du aber schon eingebunden?

Josef schrieb:
> = -1.00368

Seien wir mal kühn und sagen -1.

Dann kommt bei acos(-1) pi raus.
Und als Endergebnis 12. Ist doch etwas anderes als 5,21131484.



Karl Heinz schrieb:
> PS. die math.lib hast du aber schon eingebunden?

Müßte nicht der Linker meckern, wenn die fehlt?

CG schrieb:
> @Karl Heinz
> du meinst die math.h oder?

nein, ich meine die Linkeroption -lm

> Mich wundert es das es in Excel so ohne weiteres berechnet werden
> kann.....hmmm

du hast gerade entdeckt, dass float mit seinen 5 bis 6 signifikanten 
Stellen auch nicht das gelbe vom Ei und die Lösung aller Probleme ist.

>
> habe dann vorhin auch versucht dies zu machen:
>

1

> if(x5<-1)

2

>   x5 = x5 + 2;

3

>

>
> aber geht auch nicht so wie ich es mir vorgestellt habe

wieso + 2 ?

1

  if( x5 < -1.0 )

2

    x5 = 1.0;

3

4

  if( x5 > 1.0 )

5

    x5 = 1.0;

CG schrieb:
> @Karl Heinz
> du meinst die math.h oder?

Nein, meint er nicht.

In der math.h stehen nur die Deklarationen, damit der Compiler die 
Funktionen kennt.
In der math.lib steht der (schon compilierte) Code dazu.

Wenn du die math.lib nutzen möchtest musst du sie manchmal extra mit 
dazu linken. Die Option ist meist -lm  (minus el em)

Die Cosinusfunktion hat einen Wertebereich von einschließlich -1 bis 
einschließlich 1. Entsprechen sieht der Defintionsbereich der 
Acosfunktion aus. Beträge größer 1 sind in der Tat nicht Bestandteil des 
Definitonsbereiches dieser Funktion. Mathestoff Klasse 7 in einer POS


Namaste

Karl Heinz schrieb:

>

1

>   if( x5 < -1.0 )

2

>     x5 = 1.0;

3

>

Minus! Natürlich Minus Eins!
Blöder Tippfehler.

1

  if( x5 < -1.0 )

2

    x5 = -1.0;

3

4

  if( x5 > 1.0 )

5

    x5 = 1.0;