www.pudn.com > jaguar2s.zip > Fftmult.m

function y = fftmult(x,w,fwd); 
 
A=real(x) * cos(w); 
B=imag(x) * sin(w); 
C=imag(x) * cos(w); 
D=real(x) * sin(w); 
 
if (fwd == 1)  
   %                        -jw 
   % multiply a complex by e 
   E=A; 
   F=B; 
   G=C; 
   H=D; 
else 
   %                        jw 
   % multiply a complex by e 
   E=C; 
   F=D; 
   G=A; 
   H=B; 
end; 
 
SUM = E + F; 
DIFF = G - H; 
 
if (fwd == 1)  
    y = SUM + j * DIFF; 
else 
    y = DIFF + j * SUM; 
end;