Make refspec filter causal for overlap and add.

git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@7895 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Steven Franke 2017-07-16 00:10:37 +00:00
parent f4eaca5cf0
commit d50f3a009d
2 changed files with 15 additions and 6 deletions

View File

@ -1,11 +1,11 @@
subroutine refspectrum(id2,id2b,kk,bclear,brefspec,buseref,fname) subroutine refspectrum(id2,bclear,brefspec,buseref,fname)
! Input: ! Input:
! id2 i*2 Raw 16-bit integer data, 12000 Hz sample rate ! id2 i*2 Raw 16-bit integer data, 12000 Hz sample rate
! brefspec logical True when accumulating a reference spectrum ! brefspec logical True when accumulating a reference spectrum
parameter (NFFT=6912,NH=NFFT/2,NPOLYLOW=400,NPOLYHIGH=2600) parameter (NFFT=6912,NH=NFFT/2,NPOLYLOW=400,NPOLYHIGH=2600)
integer*2 id2(NFFT),id2b(120*12000) integer*2 id2(NFFT)
logical*1 bclear,brefspec,buseref,blastuse logical*1 bclear,brefspec,buseref,blastuse
real xs(0:NH-1) !Saved upper half of input chunk convolved with h(t) real xs(0:NH-1) !Saved upper half of input chunk convolved with h(t)
@ -16,6 +16,7 @@ subroutine refspectrum(id2,id2b,kk,bclear,brefspec,buseref,fname)
real*8 xfit(1500),yfit(1500),sigmay(1500),a(5),chisqr !Polyfit arrays real*8 xfit(1500),yfit(1500),sigmay(1500),a(5),chisqr !Polyfit arrays
logical first logical first
complex cx(0:NH) !Complex frequency-domain work array complex cx(0:NH) !Complex frequency-domain work array
complex cfil(0:NH)
character*(*) fname character*(*) fname
common/spectra/syellow(6827),ref(0:NH),filter(0:NH) common/spectra/syellow(6827),ref(0:NH),filter(0:NH)
equivalence(x,cx) equivalence(x,cx)
@ -130,15 +131,23 @@ subroutine refspectrum(id2,id2b,kk,bclear,brefspec,buseref,fname)
30 do i=1,NH 30 do i=1,NH
read(16,1005,err=100,end=100) freq,s(i),ref(i),fil(i),filter(i) read(16,1005,err=100,end=100) freq,s(i),ref(i),fil(i),filter(i)
enddo enddo
! Make the filter causal for overlap and add.
cx(0)=0.0
cx(1:NH)=fil(1:NH)/NFFT
call four2a(x,NFFT,1,1,-1)
x=cshift(x,-400)
x(800:NH)=0.0
call four2a(cx,NFFT,1,-1,0)
cfil=cx
100 close(16) 100 close(16)
110 continue 110 continue
endif endif
! Use overlap and add method to apply reference filter ! Use overlap and add method to apply causal reference filter.
x(0:NH-1)=id2(1:NH) x(0:NH-1)=id2(1:NH)
x(NH:NFFT-1)=0.0 x(NH:NFFT-1)=0.0
x=x/NFFT x=x/NFFT
call four2a(x,NFFT,1,-1,0) call four2a(x,NFFT,1,-1,0)
cx=fil*cx cx=cfil*cx
call four2a(cx,NFFT,1,1,-1) call four2a(cx,NFFT,1,1,-1)
x(0:NH-1)=x(0:NH-1)+xs x(0:NH-1)=x(0:NH-1)+xs
xs=x(NH:NFFT-1) xs=x(NH:NFFT-1)

View File

@ -124,7 +124,7 @@ extern "C" {
void wav12_(short d2[], short d1[], int* nbytes, short* nbitsam2); void wav12_(short d2[], short d1[], int* nbytes, short* nbitsam2);
void refspectrum_(short int d2[], short int d2b[], int* k, bool* bclearrefspec, void refspectrum_(short int d2[], bool* bclearrefspec,
bool* brefspec, bool* buseref, const char* c_fname, int len); bool* brefspec, bool* buseref, const char* c_fname, int len);
void freqcal_(short d2[], int* k, int* nkhz,int* noffset, int* ntol, void freqcal_(short d2[], int* k, int* nkhz,int* noffset, int* ntol,
@ -1154,7 +1154,7 @@ void MainWindow::dataSink(qint64 frames)
} }
m_bUseRef=m_wideGraph->useRef(); m_bUseRef=m_wideGraph->useRef();
refspectrum_(&dec_data.d2[k-m_nsps/2],&dec_data.d2[0],&k,&m_bClearRefSpec,&m_bRefSpec, refspectrum_(&dec_data.d2[k-m_nsps/2],&m_bClearRefSpec,&m_bRefSpec,
&m_bUseRef,c_fname,len); &m_bUseRef,c_fname,len);
m_bClearRefSpec=false; m_bClearRefSpec=false;