Algoritmus I-DFT pre delphi

[Ivo Strasil]

>> P.S. pochlubite se pak s tim algoritmem?

>Ja v nem mel nejakou chybu?

no, ted sem se na to dival, aby to odpovidalo tomu matlabu,
tak by to melo byt takhle:
(zbezne jsem to zkousel a zda se to byt OK)

for k := 0 to N-1 do
 begin
  y_re[k] := 0.0; y_im[k] := 0.0;
  for t := 0 to N-1 do
   begin
    u := k*2*pi*t/N;
    y_re[k] := y_re[k] + cos(u)*x_re[t] - sin(u)*x_im[t];
    y_im[k] := y_im[k] + cos(u)*x_im[t] + sin(u)*x_re[t];
   end;
  y_re[k]:=y_re[k]/N;
  y_im[k]:=y_im[k]/N;
 end;

Ivos
_______________________________________________
HW-list mailing list  -  sponsored by www.HW.cz
Hw-list@list.hw.cz
http://list.hw.cz/mailman/listinfo/hw-list