Algoritmus I-DFT pre delphi

[Daniel Valuch]

Tak som sa s tym cez vikend trochu pohral. Trocha matematiky na papieri 
a potom uz ani nebolo treba ziadne pocitanie s komplexnymi cislami a su 
to dva for cykly :-) Problem co som mal predtym bol v usporiadani 
vstupnych vzoriek :-)
Este aby ma niekto naucil ako sa deklaruju a pouzivaju polia s 
dynamickou velkostou a uz budem king :-P
b.



Marek Peca wrote:

>(v matlabu i s overenim, napsano za 7 minut)
>
>N=length(x);
>s=[];
>for k=1:N,
>    s(k)=0;
>    for t=1:N,
>        s(k)=s(k)+exp(j*(k-1)*2*pi*(t-1)/N)*x(t);
>    end
>    s(k)=s(k)/N;
>end
>
>Takze v Delfach, jak by to mohlo bejt...? Pozuivas nejakou komplexni
>tridu? Kdyby ne, tak asi
>
>{ vstup x_re,x_im[0..N-1] }
>
>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];
>    x_re[k] := y_re[k] + cos(u)*x_im[t] + sin(u)*y_re[t];
>   end
> end
>
>{ vystup v y_re,y_im[k], k=0..N-1, priemz zaporne frekvence lezi
>  v horni pulce pole (zkratka vraci to stejne jako matlabovska ifft }
>
>ZdraviM.P.
>
>  
>