musicdsp.org source code archive - WSInf

egin
skn:=0;
skn1:=0;
skn2:=0;
temp1:=2*PI*frequency/samplerate;
temp2:=Cos(temp1);
for i:=0 to Length(Buffer) do
begin
Skn2 = Skn1;
Skn1 = Skn;
Skn = 2*temp2*Skn1 - Skn2 + Buffer[i];
end;
Result:=(Skn - exp(-temp1)*Skn1);
end;
Maybe someone can use it...
Christian
from : GPOLAT51@YAHOO.COM
comment : HI . I NEED DSP C62X CODE FOR DTMF DETECTÝON WÝTH GOERTZEL. ÝF YOU HAVE SEND ÝT GPOLAT51@YAHOO.COM
TAHANKSSS..
from : kevin@subatomicglue.com
comment : why's the 2*cos(2*PI*frequency/samplerate) in that inner loop? It's all constants... I'd precompute this outside the for loop.
also. use the float version of cos and exp... cosf expf... that way you're not wasting time converting and using the more expensive functions.
from : net_day2004@hotmail.com
comment : if anyone please can help me, i need a c++ code about the fft.the code should input and array of N rows and 1 column, all set to 1, and
also it should input the length of fft(interval from f1 to f2) and output the results.
thank you,
i would appreciate your help,
bassem