mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2024-11-03 00:51:13 -04:00
2c17544f3f
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/WSJT/trunk@1 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
30 lines
627 B
Fortran
30 lines
627 B
Fortran
subroutine detect(data,npts,f,y)
|
|
|
|
C Compute powers at the tone frequencies using 1-sample steps.
|
|
|
|
parameter (NZ=11025,NSPD=25)
|
|
real data(npts)
|
|
real y(npts)
|
|
complex c(NZ)
|
|
complex csum
|
|
data twopi/6.283185307/
|
|
|
|
dpha=twopi*f/11025.0
|
|
do i=1,npts
|
|
c(i)=data(i)*cmplx(cos(dpha*i),-sin(dpha*i))
|
|
enddo
|
|
|
|
csum=0.
|
|
do i=1,NSPD
|
|
csum=csum+c(i)
|
|
enddo
|
|
|
|
y(1)=real(csum)**2 + imag(csum)**2
|
|
do i=2,npts-(NSPD-1)
|
|
csum=csum-c(i-1)+c(i+NSPD-1)
|
|
y(i)=real(csum)**2 + imag(csum)**2
|
|
enddo
|
|
|
|
return
|
|
end
|