Average over 3 frames on the second pass.

git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6735 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Steven Franke 2016-06-05 23:11:51 +00:00
parent 7501185723
commit 6e54009797

View File

@ -161,23 +161,19 @@ subroutine syncmsk144(cdat,npts,msgreceived,fest)
! we want ic to be the index of the first sample of the message ! we want ic to be the index of the first sample of the message
ic=ipeaks(ipk) ic=ipeaks(ipk)
! This needs to be improved - it's used to protect the edges of the array from
! overruns.
if( ic .lt. 12 .or. ic .gt. 2*864-12 ) then
! write(*,*) "Peak not in central section: ",ipk,is,ic
cycle
endif
! bb is used to place the sampling index at the center of the eye ! bb is used to place the sampling index at the center of the eye
do i=1,6 do i=1,6
io=i-3 io=i-3
bb(i) = sum( ( cdat(ic+io:ic+io+864:6) * conjg( cdat(ic+io+6:ic+io+6+864:6) ) )*2 ) ul=min(ic+io+6+864,npts)
bb(i) = sum( ( cdat(ic+io:ul-6:6) * conjg( cdat(ic+io+6:ul:6) ) )*2 )
enddo enddo
iloc=maxloc(abs(bb)) iloc=maxloc(abs(bb))
ibb=iloc(1) ibb=iloc(1)
! write(*,*) 'ic0: ',ic,'bb peak is at : ',ibb ! write(*,*) 'ic0: ',ic,'bb peak is at : ',ibb
! Adjust frame index to place peak of bb at desired lag ! Adjust frame index to place peak of bb at desired lag
ic=ic + ibb-2+is ic=ic + ibb-2+is
if( ic .le. 864 ) ic=ic+864
if( ic .gt. 2*864 ) ic=ic-864
! Sanity check - recompute bb and verify that peak is now at designated lag. ! Sanity check - recompute bb and verify that peak is now at designated lag.
! do i=1,6 ! do i=1,6
@ -188,17 +184,17 @@ subroutine syncmsk144(cdat,npts,msgreceived,fest)
! ibb=iloc(1) ! ibb=iloc(1)
! write(*,*) 'ic1: ',ic,'bb peak is at : ',ibb ! write(*,*) 'ic1: ',ic,'bb peak is at : ',ibb
! Average two frames on the second pass only if its incredibly easy to do ! Average three frames on the second pass - truncation of the high frame
! better than nothing. Should be improved. ! is not properly accounted for.
c=cdat(ic:ic+864-1) c=cdat(ic:ic+864-1)
if( iav .eq. 1 ) then if( iav .eq. 1 ) then
id0=ic+864 id0=ic+864
id1=ic+864+863 id1=min(ic+864+863,npts)
if( id1 .le. npts ) then np=id1-id0+1
c=c+cdat(id0:id1) c(1:np)=c(1:np)+cdat(id0:id1)
else id0=ic-864
cycle id1=ic-1
endif c=c+cdat(id0:id1)
endif endif
! Estimate fine frequency error. ! Estimate fine frequency error.