Logarithmus Log zu Delphi

[EWeiss]

Zitat von Sherlock:

Stimmt, und das ist nach dem ln.
und ich habe vom ln an sich zu wenig Ahnung.

Sherlock

Mir geht's genauso..
Fand das Original nicht schlecht und wollte es halt mal nach Delphi umsetzen (Nebenbei Vaporware ist das nicht) zu viel Aufwand wie man sehen kann.
Und man kann von lernen.

Mein Hauptproblem zur zeit ist das der Linker mich einfach rausschmeißt.

Delphi-Quellcode:

procedure TSpectrum.FFT(Dat: array of TComplex);
var
  i, j, n, K , io, ie, in_, nn: Integer;
  u, tp, tq, w: TComplex;
  sr, mGain: Single;

 function LimitedSingleValue(doubleVal: double): single; stdcall;
 begin
   if doubleVal > MaxSingle then
     result := MaxSingle
   else if doubleVal < MinSingle then
     result := MinSingle
   else
     result := doubleVal;
 end;

begin

Der linker entfernt variablen.

Delphi-Quellcode:

  mGain :=  (FGain / 100);
  if FView = 0 then
    sr := (4096 * mGain) / FFFTSize
  else
  sr := 1 / FFFTSize;

hier FFFTSize.

EurekaLog springt dann mit einem Access Violation auf das letzte end im oberen Quelltext. Mit der Meldung vom Speicher kann nicht gelesen werden.
Da ist einiges im argen.

Zitat:

.. wie groß können denn Spectrum[i1].r und Spectrum[i1].i maximal werden?

Kann ich leider mit Bestimmtheit nicht sagen.

Wenn Interesse besteht kann ich den Quelltext zum debuggen aber gerne mal hochladen.
Oder ich schicke einen Link per PN.

gruss

[Delphi.Narium]

Nur 'ne blöde Idee, ein Schuss ins Blaue:

Mach' bitte mal aus

Delphi-Quellcode:

  mGain := (FGain / 100);
  if FView = 0 then
    sr := (4096 * mGain) / FFFTSize
  else
  sr := 1 / FFFTSize;

dashier:

Delphi-Quellcode:

  mGain := (FGain / 100);
  sr := IfThen(FView = 0, 4096 * mGain, 1) / FFFTSize;

Geht das und ändert sich dann was?

[EWeiss]

IfThen(FView = 0, (4096 * mGain) / FFFTSize, 1);

Ist am ende das gleiche oder?

Das Ergebnis ist das gleiche.
Habe es versucht. Danke.

gruss

[Delphi.Narium]

nicht ganz:
entweder IfThen(FView = 0, (4096 * mGain) / FFFTSize, 1 / FFFTSize);
oder IfThen(FView = 0, (4096 * mGain), 1) / FFFTSize;

Im ursprünglichen Quelltext findet für sr ja auch im Else noch eine Division durch FFFTSize statt, muss also im IfThen entweder in beiden Fällen durchgeführt werden oder auf das Ergebnis von IfThen.

[EWeiss]

IfThen(FView = 0, (4096 * mGain), 1) / FFFTSize;
Sorry aber das funktioniert nicht weil IFThen so nicht ausgelegt ist.

Zitat:

[DCC Error] uSpectrum.pas(512): E2014 Statement expected, but expression of type 'Extended' found

Deshalb hatte ich meine Zeile angepasst.

gruss

[Delphi.Narium]

Ok, dann gibt es keine überladene Version von IfThen für den passenden Datentypen. Schade.

[Michael II]

Was passiert, wenn du deine Formel

Delphi-Quellcode:

  c := 0.0001;
  d := -ln(c);
...
  b2 := (ln(Sqrt(r * r + i * i) + c ) + d ) * q1;

Schritt für Schritt auswertest?

Also so:

Delphi-Quellcode:

  c := 0.0001;
  d := -ln(c);
...
  b := Sqrt(r * r + i * i);
  b := ln( b + c );
  b := b + d;
  b := b*q1;

Wo haut's dich da raus?


Nebenbei:
Wenn du statt d := 9.21034037197618; (als single Wert speichert Delphi 9.21034049987793) d := -ln(0.0001); verwendest, dann wird b oben für r=i=0 wie gewünscht 0 [auch wenn du dich später für "double" entscheiden solltest].

Man darf es Resultatkosmetik nennen...