w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2025-08-03 14:42:25 -04:00
Code Issues Packages Projects Releases Wiki Activity

Align displayed S/N with S/N in qra64sim.

Browse Source 
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6851 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Joe Taylor 2016-07-02 15:35:31 +00:00
parent eb77d3d65a
commit e195b3ce15
2 changed files with 11 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

29
lib/qra/qra64/qra64sim.f90
View File

@ -9,21 +9,19 @@ program qra64sim
parameter (NFFT=10*65536,NH=NFFT/2) parameter (NFFT=10*65536,NH=NFFT/2)
type(hdr) h !Header for .wav file type(hdr) h !Header for .wav file
integer*2 iwave(NMAX) !Generated waveform integer*2 iwave(NMAX) !Generated waveform
! integer*4 itone(126) !Channel symbols (values 0-65)
integer*4 itone(84) !Channel symbols (values 0-63) integer*4 itone(84) !Channel symbols (values 0-63)
integer dgen(12) !Twelve 6-bit data symbols
integer sent(63) !RS(63,12) codeword
real*4 xnoise(NMAX) !Generated random noise real*4 xnoise(NMAX) !Generated random noise
real*4 dat(NMAX) !Generated real data real*4 dat(NMAX) !Generated real data
complex cdat(NMAX) !Generated complex waveform complex cdat(NMAX) !Generated complex waveform
complex cspread(0:NFFT-1) !Complex amplitude for Rayleigh fading complex cspread(0:NFFT-1) !Complex amplitude for Rayleigh fading
complex z complex z
real*8 f0,dt,twopi,phi,dphi,baud,fsample,freq,sps real*8 f0,dt,twopi,phi,dphi,baud,fsample,freq
character msg*22,fname*11,csubmode*1,c,optarg*500,numbuf*32 character msg*22,fname*11,csubmode*1,c,optarg*500,numbuf*32
character msgsent*22 character msgsent*22
logical :: display_help=.false.,seed_prngs=.true. logical :: display_help=.false.,seed_prngs=.true.
type (option) :: long_options(8) = [ & type (option) :: long_options(9) = [ &
option ('help',.false.,'h','Display this help message',''), & option ('help',.false.,'h','Display this help message',''), &
option ('sub-mode',.true.,'M','message','MSG'), &
option ('sub-mode',.true.,'m','sub mode, default MODE=A','MODE'), & option ('sub-mode',.true.,'m','sub mode, default MODE=A','MODE'), &
option ('num-sigs',.true.,'n','number of signals per file, default SIGNALS=10','SIGNALS'), & option ('num-sigs',.true.,'n','number of signals per file, default SIGNALS=10','SIGNALS'), &
option ('doppler-spread',.true.,'d','Doppler spread, default SPREAD=0.0','SPREAD'), & option ('doppler-spread',.true.,'d','Doppler spread, default SPREAD=0.0','SPREAD'), &
@ -42,13 +40,16 @@ program qra64sim
nfiles=1 nfiles=1
do do
call getopt('hm:n:d:t:f:ps:',long_options,c,optarg,narglen,nstat,noffset,nremain,.true.) call getopt('hM:m:n:d:t:f:ps:',long_options,c,optarg,narglen,nstat,noffset,nremain,.true.)
if( nstat .ne. 0 ) then if( nstat .ne. 0 ) then
exit exit
end if end if
select case (c) select case (c)
case ('h') case ('h')
display_help = .true. display_help = .true.
case ('M')
read (optarg(:narglen), *) msg
print*,msg
case ('m') case ('m')
read (optarg(:narglen), *) csubmode read (optarg(:narglen), *) csubmode
if(csubmode.eq.'A') mode65=1 if(csubmode.eq.'A') mode65=1
@ -109,13 +110,11 @@ program qra64sim
twopi=8.d0*atan(1.d0) twopi=8.d0*atan(1.d0)
npts=54*12000 !Total samples in .wav file npts=54*12000 !Total samples in .wav file
nsps=6912 nsps=6912
! baud=11025.d0/4096.d0 !Keying rate baud=12000.d0/nsps !Keying rate = 1.7361111111
baud=12000.d0/6912 !Keying rate = 1.7361111111
! sps=12000.d0/baud !Samples per symbol, at fsample=12000 Hz
! nsym=126 !Number of channel symbols
nsym=84 !Number of channel symbols nsym=84 !Number of channel symbols
h=default_header(12000,npts) h=default_header(12000,npts)
dfsig=2000.0/nsigs !Freq spacing between sigs in file (Hz) dfsig=2000.0/nsigs !Freq spacing between sigs in file (Hz)
ichk=0
do ifile=1,nfiles !Loop over requested number of files do ifile=1,nfiles !Loop over requested number of files
write(fname,1002) ifile !Output filename write(fname,1002) ifile !Output filename
@ -130,7 +129,7 @@ program qra64sim
enddo enddo
endif endif
msg="K1ABC W9XYZ EN37" ! msg="K1ABC W9XYZ EN37"
do isig=1,nsigs !Generate requested number of sigs do isig=1,nsigs !Generate requested number of sigs
if(mod(nsigs,2).eq.0) f0=1500.0 + dfsig*(isig-0.5-nsigs/2) if(mod(nsigs,2).eq.0) f0=1500.0 + dfsig*(isig-0.5-nsigs/2)
if(mod(nsigs,2).eq.1) f0=1500.0 + dfsig*(isig-(nsigs+1)/2) if(mod(nsigs,2).eq.1) f0=1500.0 + dfsig*(isig-(nsigs+1)/2)
@ -139,16 +138,8 @@ program qra64sim
if(csubmode.eq.'B' .and. snrdb.eq.0.0) xsnr=-21 - isig if(csubmode.eq.'B' .and. snrdb.eq.0.0) xsnr=-21 - isig
if(csubmode.eq.'C' .and. snrdb.eq.0.0) xsnr=-21 - isig if(csubmode.eq.'C' .and. snrdb.eq.0.0) xsnr=-21 - isig
xsnr=xsnr+5 !### TEMPORARY ###
call genqra64(msg,ichk,msgsent,itone,itype) call genqra64(msg,ichk,msgsent,itone,itype)
! call packmsg(msg,dgen,itype) !Pack message into 12 six-bit bytes
! call qra64_enc(dgen,sent) !Encode using QRA64
!! call qra64_dec(sent,dgen,ierr) !Decode (### for test only ###)
!! write(*,3001) sent
!!3001 format(21i3)
bandwidth_ratio=2500.0/6000.0 bandwidth_ratio=2500.0/6000.0
sig=sqrt(2*bandwidth_ratio)*10.0**(0.05*xsnr) sig=sqrt(2*bandwidth_ratio)*10.0**(0.05*xsnr)
if(xsnr.gt.90.0) sig=1.0 if(xsnr.gt.90.0) sig=1.0

2
lib/qra64a.f90
View File

@ -118,7 +118,7 @@ subroutine qra64a(dd,nf1,nf2,nfqso,ntol,mycall_12,sync,nsnr,dtx,nfreq, &
enddo enddo
enddo enddo
if(sync.gt.1.0) nsnr=nint(10.0*log10(sync) - 39.0) if(sync.gt.1.0) nsnr=nint(10.0*log10(sync) - 38.0)
! if(sync.lt.12.8) go to 900 !### Temporary ### ! if(sync.lt.12.8) go to 900 !### Temporary ###
mycall=mycall_12(1:6) !### May need fixing ### mycall=mycall_12(1:6) !### May need fixing ###