w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2025-07-12 22:55:18 -04:00
Code Issues Packages Projects Releases Wiki Activity

Do polyfit to fil vector in refspectrum.f90 and write the coeffs to the first line of refspec.dat.

Browse Source 
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@7576 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Steven Franke 2017-02-22 20:43:40 +00:00
parent 7383ba46f5
commit a659b59c7e
1 changed files with 16 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
lib/refspectrum.f90
View File

@ -7,7 +7,7 @@ subroutine refspectrum(id2,bclear,brefspec,buseref,fname)
parameter (NFFT=6912,NH=NFFT/2) parameter (NFFT=6912,NH=NFFT/2)
integer*2 id2(NFFT) integer*2 id2(NFFT)
logical*1 bclear,brefspec,buseref logical*1 bclear,brefspec,buseref
real x0(0:NH-1) !Input samples real x0(0:NH-1) !Input samples
real x1(0:NH-1) !Output samples (delayed by one block) real x1(0:NH-1) !Output samples (delayed by one block)
real x0s(0:NH-1) !Saved upper half of input samples real x0s(0:NH-1) !Saved upper half of input samples
@ -16,6 +16,7 @@ subroutine refspectrum(id2,bclear,brefspec,buseref,fname)
real*4 w(0:NFFT-1) !Window function real*4 w(0:NFFT-1) !Window function
real*4 s(0:NH) !Average spectrum real*4 s(0:NH) !Average spectrum
real*4 fil(0:NH) real*4 fil(0:NH)
real*8 xfit(1500),yfit(1500),sigmay(1500),a(5),chisqr !Polyfit arrays
logical first,firstuse logical first,firstuse
complex cx(0:NH) !Complex frequency-domain work array complex cx(0:NH) !Complex frequency-domain work array
character*(*) fname character*(*) fname
@ -92,12 +93,25 @@ subroutine refspectrum(id2,bclear,brefspec,buseref,fname)
call smo121(fil,NH) call smo121(fil,NH)
enddo enddo
do i=0,NH do i=0,NH
filter(i)=-60.0 filter(i)=-60.0
if(s(i).gt.0.0) filter(i)=20.0*log10(fil(i)) if(s(i).gt.0.0) filter(i)=20.0*log10(fil(i))
enddo enddo
il=nint(400.0/df)
ih=nint(2600.0/df)
nfit=ih-il+1
mode=0
nterms=5
do i=1,nfit
xfit(i)=((i+il-1)*df-1500.0)/1000.0
yfit(i)=fil(i+il-1)
sigmay(i)=1.0
enddo
call polyfit(xfit,yfit,sigmay,nfit,nterms,mode,a,chisqr)
open(16,file=fname,status='unknown') open(16,file=fname,status='unknown')
write(16,'(5f10.4)') a
do i=1,NH do i=1,NH
freq=i*df freq=i*df
ref(i)=db(s(i)/avemid) ref(i)=db(s(i)/avemid)