Algoritmus I-DFT pre delphi

[Marek Peca]

Nojo, ja tam v zbrklosti pomotal x a y. Diky za opravu (snad Dr. Balu
neobetoval vikend kvuli tomu :-D).

mp

On Mon, 6 Jun 2005, Ivo Strasil wrote:

> >> 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
> 
> _______________________________________________
> HW-list mailing list  -  sponsored by www.HW.cz
> Hw-list@list.hw.cz
> http://list.hw.cz/mailman/listinfo/hw-list
>