program jt9sim ! Generate simulated data for testing of WSJT-X parameter (NMAX=1800*12000) integer ihdr(11) integer*2 iwave !Generated waveform (no noise) real*8 f0,f,dt,twopi,phi,dphi,baud,fspan character msg*22,msg0*22,message*22,msgsent*22,arg*8,fname*11 integer*4 d6(85) integer*4 mettab(0:255,0:1) integer*4 t0(13) !72-bit message as 6-bit words integer*1 t1(13) !72 bits and zero tail as 8-bit bytes integer*1 t2(207) !Encoded information-carrying bits integer*1 t3(207) !Bits from t4, after interleaving integer*4 t4(69) !Symbols from t5, values 0-7 integer*4 t5(69) !Gray-coded symbols, values 0-7 integer*4 t6(85) !Channel symbols including sync, values 0-8 integer*1 tmp(72) integer*1 i1 equivalence (i1,i4) integer isync(85) !Sync vector data isync/ & 1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,0, & 1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,0,0,0,1,0, & 0,0,0,0,1,0,0,0,0,0,1,0,0,0,0,0,1,0,0,0, & 0,0,1,0,0,0,0,0,1,0,0,0,0,0,1,0,0,0,0,0, & 1,0,0,0,1/ common/acom/dat(NMAX),iwave(NMAX) nargs=iargc() if(nargs.ne.6) then print*,'Usage: jt9sim "message" fspan nsigs minutes SNR nfiles' print*,'Example: "CQ K1ABC FN42" 200 20 2 -28 1' print*,' ' print*,'Enter message = "" to use entries in msgs.txt.' print*,'Enter SNR = 0 to generate a range of SNRs.' go to 999 endif call getarg(1,msg0) message=msg0 !Transmitted message call getarg(2,arg) read(arg,*) fspan !Total freq range (Hz) call getarg(3,arg) read(arg,*) nsigs !Number of signals in each file call getarg(4,arg) read(arg,*) minutes !Length of file (1 2 5 10 30 minutes) call getarg(5,arg) read(arg,*) snrdb !S/N in dB (2500 hz reference BW) call getarg(6,arg) read(arg,*) nfiles !Number of files rmsdb=25. rms=10.0**(0.05*rmsdb) f0=1500.d0 !Center frequency (MHz) fsample=12000.d0 !Sample rate (Hz) dt=1.d0/fsample !Sample interval (s) twopi=8.d0*atan(1.d0) npts=12000*(60*minutes-6) nsps=0 if(minutes.eq.1) nsps=6912 if(minutes.eq.2) nsps=15360 if(minutes.eq.5) nsps=40960 if(minutes.eq.10) nsps=82944 if(minutes.eq.30) nsps=250880 if(nsps.eq.0) stop 'Bad value for minutes.' ihdr=0 !Temporary ### open(12,file='msgs.txt',status='old') ! Get the metric table bias=0.37 scale=10 !Optimize? open(19,file='met8.21',status='old') do i=0,255 read(19,*) x00,x0,x1 mettab(i,0)=nint(scale*(x0-bias)) mettab(i,1)=nint(scale*(x1-bias)) !### Check range, etc. ### enddo close(19) write(*,1000) 1000 format('File N freq S/N Message'/ & '---------------------------------------------------') do ifile=1,nfiles nmin=(ifile-1)*minutes ihr=nmin/60 imin=mod(nmin,60) write(fname,1002) ihr,imin !Create the output filenames 1002 format('000000_',2i2.2) open(10,file=fname//'.wav',access='stream',status='unknown') if(snrdb.lt.90) then do i=1,npts dat(i)=gran() enddo else dat(1:npts)=0. endif if(msg0.ne.' ') then call genjt9(message,minutes,msgsent,d6) endif rewind 12 do isig=1,nsigs if(msg0.eq.' ') then read(12,1004) message 1004 format(a22) call genjt9(message,minutes,msgsent,d6) endif f=f0 if(nsigs.gt.1) f=f0 - 0.5d0*fspan + fspan*(isig-1.d0)/(nsigs-1.d0) snrdbx=snrdb ! if(snrdb.ge.-1.0) snrdbx=-15.0 - 15.0*(isig-1.0)/nsigs sig=sqrt(2500.0/12000.0) * 10.0**(0.05*snrdbx) write(*,1020) ifile,isig,f,snrdbx,msgsent 1020 format(i3,i4,f10.3,f7.1,2x,a22) phi=0. baud=12000.0/nsps k=12000 !Start at t = 1 s do isym=1,85 freq=f + d6(isym)*baud dphi=twopi*freq*dt do i=1,nsps phi=phi + dphi if(phi.lt.-twopi) phi=phi+twopi if(phi.gt.twopi) phi=phi-twopi xphi=phi k=k+1 dat(k)=dat(k) + sig*sin(xphi) enddo enddo enddo do i=1,npts iwave(i)=nint(rms*dat(i)) enddo write(10) ihdr,iwave(1:npts) close(10) ! We're done! Now decode the data in d6, as a test. j=0 do i=1,85 if(isync(i).eq.1) cycle j=j+1 t5(j)=d6(i)-1 enddo call graycode(t5,69,-1,t4) call unpackbits(t4,69,3,t3) call interleave9(t3,-1,t2) nbits=72 ndelta=17 limit=1000 do i=1,206 i4=-10 if(t2(i).eq.1) i4=10 i4=i4+128 t2(i)=i1 enddo call fano232(t2,nbits+31,mettab,ndelta,limit,t1,ncycles,metric,ierr, & maxmetric,maxnp) nbytes=(nbits+7)/8 do i=1,nbytes n=t1(i) t4(i)=iand(n,255) enddo call unpackbits(t4,nbytes,8,tmp) call packbits(tmp,12,6,t4) do i=1,12 if(t4(i).lt.128) t1(i)=t4(i) if(t4(i).ge.128) t1(i)=t4(i)-256 enddo do i=1,12 t4(i)=t1(i) enddo call unpackmsg(t4,msg) !Unpack decoded msg if(msg.ne.msg0) print*,'Decode error: ',msg0,' ',msg enddo 999 end program jt9sim