Fehler beim erstellen einer Sinus-Tabelle

[Trigiton]

Hallo zusammen,

ich hab ein kleines Problem, ich will in einem Progamm eine ganz simple Tabelle mit Sinuswerten erzeugen. Hier mal der Code:

Delphi-Quellcode:

procedure TForm1.Button1Click(Sender: TObject);
var
  i : Integer;
  z, x : single;

begin
  x := 0; z := 0;
  for i:= 0 to 200 do begin
    x := sin(z);
    z := z + ((2 * Pi) / 200);
    Memo1.Lines.Add(IntToStr(i)+' :  '+FloatToStr(x));
  end;
end;

Das Problem ist das bei i=100 und i=200, für x unsinnige werte raus kommen: -1,1... und -2,2...
Ich weis das es schon mal einen Thread zu diesem Thema gab, finde ihn aber nicht mehr.

Kann mir einer dabei helfen?

Gruß
Trigiton

[glkgereon]

ich finds irgendwie komisch, denn der wert von i hat doch keinen einfluss auf das x, oder? :gruebel:

[Christian Seehase]

Moin Trigiton,

also bei mir kommt folgendes heraus:

Code:

100 :  -1,04109710719058[color=red][b]E-6[/b][/color]

Code:

200 :  2,08219421438116[color=red][b]E-6[/b][/color]

kann es sein, dass Du das E-6 übersehen hast?

[Trigiton]

Hallo...

@christian: Ja hab vergessen die mit reinzuschreiben, bei mir kommen natürlich auch die E-6'er vor, aber das ändert ja nix dran das es falsch ist, schliesslich müsste normalerweise an den stellen Null stehen, wenn ich in mathe aufgepasst habe :-D

@glkbkk: Es hat auch nicht direkt etwas mit dem I zu tun, sondern an der tatsache das es sich dort um die Nulldurchgänge der Kurve handelt, und die rechnet Delphi nicht korrekt aus.

Gruß
Trigiton

[Luckie]

Was rechnest du eigentlich hier:

Delphi-Quellcode:

z := z + ((2 * Pi) / 200);

aus? Und wo kommt die 200 her? Deg nach Rad ist es nicht.

[Christian Seehase]

Moin Trigiton,

dann nimm' mal als Datentyp Extended und nicht single.
Da es ja auch Rundungsdifferenzen gibt, lass Dir dass dann mal mit

Delphi-Quellcode:

    Memo1.Lines.Add(Format(IntToStr(i)+' : %.16f',[x]));

ins Memo schreiben.

[Trigiton]

Hallo...

die 200 kommen aus der Schleife -> 200 werte, 2*pi ist ja ein voller sinus.
Also errechne ich die schrittweite für Z in dem ich einfach immer 2*pi/200 addiere.

Damit bekomme ich aber dann das problem das bei den nulldurchgängen z nicht korrekt ist damit sin(z) = 0 ergibt.

So stehts zumindest in der Formelsammlung die hier rumliegt :-)


Gruß
Trigiton

[Christian Seehase]

Moin Trigiton,

Da die Schleife von 0 bis 200 läuft, hat sie 201 Durchläufe ;-)

[Trigiton]

Hallo...

@Christian:
Dachte ich mir schon das es sich um rundungsfehler handelt, an extended hab ich nicht gedacht.
Stimmt die schleife läuft 201mal aber dann kommt halt trotzdem bei 101 noch ein fehler.

Aber mit dem extended-typen ist es jetzt für mich genau genug. Danke.

Gruß
Trigiton

[Dagon]

Noch eine winzig kleine Bemerkung (nur für Mathematiker :wink: ):

Delphi-Quellcode:

z := z + ((2 * Pi) / 200);

ist das selbe wie

Delphi-Quellcode:

z := z + Pi/100;