From 29599830ebacd2f1b2f45fcf9efb80c2d3fa85ca Mon Sep 17 00:00:00 2001 From: Joe Taylor Date: Wed, 5 Sep 2012 18:02:50 +0000 Subject: [PATCH] Remove DC component of symbol powers before computing FFT. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@2559 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 --- libm65/ccf65.f90 | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/libm65/ccf65.f90 b/libm65/ccf65.f90 index 75a552f3c..f69cbd249 100644 --- a/libm65/ccf65.f90 +++ b/libm65/ccf65.f90 @@ -60,6 +60,8 @@ subroutine ccf65(ss,nhsym,nfast,ssmax,sync1,ipol1,jpz,dt1,flipk, & ! s(i)=ss(ip,i)+ss(ip,i+1) s(i)=min(ssmax,ss(ip,i)+ss(ip,i+1)) enddo + call pctile(s,tmp1,nhsym-1,50,base) + s(1:nhsym-1)=s(1:nhsym-1)-base s(nhsym:NFFT)=0. call four2a(s,NFFT,1,-1,0) !Real-to-complex FFT do i=0,NH @@ -68,7 +70,7 @@ subroutine ccf65(ss,nhsym,nfast,ssmax,sync1,ipol1,jpz,dt1,flipk, & enddo call four2a(cs,NFFT,1,1,-1) !Complex-to-real inv-FFT call four2a(cs2,NFFT,1,1,-1) !Complex-to-real inv-FFT - + do lag=-27,27 !Check for best JT65 sync ccf(lag,ip)=s(lag+28) if(abs(ccf(lag,ip)).gt.ccfbest) then