w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2025-10-31 13:10:19 -04:00
Code Issues Packages Projects Releases Wiki Activity

Build synctest from top-level CMakeLists.txt. Add 'nfiles' capability.

Browse Source
This commit is contained in:
Joe Taylor 2024-01-31 15:13:42 -05:00
parent 0531f339e8
commit b0bc73e3f3
2 changed files with 59 additions and 46 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
CMakeLists.txt
View File

@ -586,6 +586,10 @@ set (wsjt_FSRCS
lib/fst4/get_crc24.f90 lib/fst4/get_crc24.f90
lib/fst4/fst4_baseline.f90 lib/fst4/fst4_baseline.f90
lib/77bit/hash22calc.f90 lib/77bit/hash22calc.f90
lib/superfox/gen_sfox.f90
lib/superfox/sync_sf.f90
lib/superfox/hard_symbols.f90
lib/superfox/ran1.f90
) )
# temporary workaround for a gfortran v7.3 ICE on Fedora 27 64-bit # temporary workaround for a gfortran v7.3 ICE on Fedora 27 64-bit

101
lib/superfox/synctest.f90
View File

@ -18,9 +18,9 @@ program synctest
character fname*17,arg*12 character fname*17,arg*12
nargs=iargc() nargs=iargc()
if(nargs.ne.6) then if(nargs.ne.8) then
print*,'Usage: synctest f0 DT fspread delay width snr' print*,'Usage: synctest f0 DT fspread delay width nran nfiles snr'
print*,'Example: synctest 1500.0 2.5 0.0 0.0 100 -20' print*,'Example: synctest 1500.0 2.5 0.0 0.0 100 0 10 -20'
go to 999 go to 999
endif endif
call getarg(1,arg) call getarg(1,arg)
@ -34,71 +34,80 @@ program synctest
call getarg(5,arg) call getarg(5,arg)
read(arg,*) syncwidth read(arg,*) syncwidth
call getarg(6,arg) call getarg(6,arg)
read(arg,*) nran
call getarg(7,arg)
read(arg,*) nfiles
call getarg(8,arg)
read(arg,*) snrdb read(arg,*) snrdb
rms=100. rms=100.
fsample=12000.0 !Sample rate (Hz) fsample=12000.0 !Sample rate (Hz)
baud=12000.0/nsps !Keying rate, 11.719 baud for nsps=1024 baud=12000.0/nsps !Keying rate, 11.719 baud for nsps=1024
idummy=0
do i=1,ND
idat(i)=128*ran1(idummy)
enddo
h=default_header(12000,NMAX) h=default_header(12000,NMAX)
fname='000000_000001.wav' idummy=0
open(10,file=trim(fname),access='stream',status='unknown')
xnoise=0.
cnoise=0.
if(snrdb.lt.90) then
do i=1,NMAX !Generate Gaussian noise
x=gran()
y=gran()
xnoise(i)=x
cnoise(i)=cmplx(x,y)
enddo
endif
bandwidth_ratio=2500.0/6000.0 bandwidth_ratio=2500.0/6000.0
sig=sqrt(2*bandwidth_ratio)*10.0**(0.05*snrdb) sig=sqrt(2*bandwidth_ratio)*10.0**(0.05*snrdb)
if(snrdb.gt.90.0) sig=1.0 if(snrdb.gt.90.0) sig=1.0
ngood=0
do ifile=1,nfiles
do i=1,ND
call random_number(r)
if(nran.eq.1) r=ran1(idummy)
idat(i)=128*r
enddo
fname='000000_000001.wav'
write(fname(8:13),'(i6.6)') ifile
open(10,file=trim(fname),access='stream',status='unknown')
xnoise=0.
cnoise=0.
if(snrdb.lt.90) then
do i=1,NMAX !Generate Gaussian noise
x=gran()
y=gran()
xnoise(i)=x
cnoise(i)=cmplx(x,y)
enddo
endif
!Generate cdat (SuperFox waveform) and clo (LO needed for sync detection) !Generate cdat (SuperFox waveform) and clo (LO needed for sync detection)
call gen_sfox(idat,f0,fsample,syncwidth,cdat,clo) call gen_sfox(idat,f0,fsample,syncwidth,cdat,clo)
crcvd=0. crcvd=0.
crcvd(1:NMAX)=cshift(sig*cdat(1:NMAX),-nint(xdt*fsample)) + cnoise crcvd(1:NMAX)=cshift(sig*cdat(1:NMAX),-nint(xdt*fsample)) + cnoise
dat=aimag(sig*cdat(1:NMAX)) + xnoise !Add generated AWGN noise dat=aimag(sig*cdat(1:NMAX)) + xnoise !Add generated AWGN noise
fac=32767.0 fac=32767.0
if(snrdb.ge.90.0) iwave(1:NMAX)=nint(fac*dat(1:NMAX)) if(snrdb.ge.90.0) iwave(1:NMAX)=nint(fac*dat(1:NMAX))
if(snrdb.lt.90.0) iwave(1:NMAX)=nint(rms*dat(1:NMAX)) if(snrdb.lt.90.0) iwave(1:NMAX)=nint(rms*dat(1:NMAX))
write(10) h,iwave(1:NMAX) !Save the .wav file write(10) h,iwave(1:NMAX) !Save the .wav file
close(10) close(10)
if(fspread.ne.0 .or. delay.ne.0) call watterson(crcvd,NMAX,NZ,fsample, & if(fspread.ne.0 .or. delay.ne.0) call watterson(crcvd,NMAX,NZ,fsample, &
delay,fspread) delay,fspread)
! Find signal freq and DT ! Find signal freq and DT
call sync_sf(crcvd,clo,f,t) call sync_sf(crcvd,clo,f,t)
call hard_symbols(crcvd,f,t,jdat) call hard_symbols(crcvd,f,t,jdat)
nharderr=count(jdat.ne.idat) nharderr=count(jdat.ne.idat)
write(*,1100) f0,xdt write(*,1100) f0,xdt
1100 format(/'f0:',f7.1,' xdt:',f6.2) 1100 format(/'f0:',f7.1,' xdt:',f6.2)
write(*,1112) f,t write(*,1112) f,t
1112 format('f: ',f7.1,' DT:',f6.2) 1112 format('f: ',f7.1,' DT:',f6.2)
write(*,1110) f-f0,t-xdt write(*,1110) f-f0,t-xdt
1110 format('err:',f6.1,f12.2) 1110 format('err:',f6.1,f12.2)
write(*,1120) nharderr write(*,1120) nharderr
1120 format('Hard errors:',i4) 1120 format('Hard errors:',i4)
if(nharderr.le.38) ngood=ngood+1
! write(13,1200) snrdb,nharderr
!1200 format(f7.2,i5)
enddo
write(*,1300) snrdb,nfiles,ngood,float(ngood)/nfiles
1300 format(f7.2,2i5,f7.2)
999 end program synctest 999 end program synctest
include 'gen_sfox.f90'
include 'sync_sf.f90'
include 'hard_symbols.f90'
include 'ran1.f90'