From 8772a101dcd2e0ea0319910c92e57e96d680386a Mon Sep 17 00:00:00 2001 From: Steven Franke Date: Sun, 8 Apr 2018 17:07:08 +0000 Subject: [PATCH] Add some development routines. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@8604 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 --- lib/fsk4hf/genwsprdpsk.f90 | 57 ++++++++++++ lib/fsk4hf/ldpcsim204.f90 | 2 +- lib/fsk4hf/wsprdpsksim.f90 | 175 +++++++++++++++++++++++++++++++++++++ 3 files changed, 233 insertions(+), 1 deletion(-) create mode 100644 lib/fsk4hf/genwsprdpsk.f90 create mode 100644 lib/fsk4hf/wsprdpsksim.f90 diff --git a/lib/fsk4hf/genwsprdpsk.f90 b/lib/fsk4hf/genwsprdpsk.f90 new file mode 100644 index 000000000..353ef1339 --- /dev/null +++ b/lib/fsk4hf/genwsprdpsk.f90 @@ -0,0 +1,57 @@ +subroutine genwsprdpsk(msg,msgsent,imessage) + +! Encode a WSPRDPSK message, producing array txwave(). +! + use crc + include 'wsprdpsk_params.f90' + + character*22 msg,msgsent + character*64 cbits + character*32 sbits + integer iuniqueword0 + integer*1,target :: idat(9) + integer*1 msgbits(68),codeword(ND) + logical first + integer ipreamble(16) !Freq estimation preamble + integer isync(32) !Long sync vector + integer imessage(NN) + data ipreamble/1,1,1,1,0,0,0,0,0,0,0,0,1,1,1,1/ + data first/.true./ + data iuniqueword0/z'30C9E8AD'/ + save first,isync,ipreamble + + if(first) then + write(sbits,'(b32.32)') iuniqueword0 + read(sbits,'(32i1)') isync(1:32) + first=.false. + endif + + idat=0 + call wqencode(msg,ntype0,idat) !Source encoding + id7=idat(7) + if(id7.lt.0) id7=id7+256 + id7=id7/64 +write(*,*) 'idat ',idat + icrc=crc14(c_loc(idat),9) +write(*,*) 'icrc: ',icrc +write(*,'(a6,b16.16)') 'icrc: ',icrc + call wqdecode(idat,msgsent,itype) + print*,msgsent,itype + write(cbits,1004) idat(1:6),id7,iand(icrc,z'3FFF') +1004 format(6b8.8,b2.2,b14.14) + msgbits=0 + read(cbits,1006) msgbits(1:64) +1006 format(64i1) + +write(*,'(50i1,1x,14i1,1x,4i1)') msgbits + + call encode204(msgbits,codeword) !Encode the test message + +! Message structure: +! d100 p16 d100 + imessage(1:100)=codeword(1:100) + imessage(101:132)=isync + imessage(133:232)=codeword(101:200) + + return +end subroutine genwsprdpsk diff --git a/lib/fsk4hf/ldpcsim204.f90 b/lib/fsk4hf/ldpcsim204.f90 index 224b2624f..794b1476b 100644 --- a/lib/fsk4hf/ldpcsim204.f90 +++ b/lib/fsk4hf/ldpcsim204.f90 @@ -227,7 +227,7 @@ do idb = 20,-18,-1 endif endif enddo - snr2500=db+10*log10(1.389/2500.0) + snr2500=db+10*log10(200.0/116.0/2500.0) pberr=real(nberr)/(real(ntrials*N)) write(*,"(f4.1,4x,f5.1,1x,i8,1x,i8,1x,i8,8x,f5.2,8x,e10.3)") db,snr2500,ngood,nue,nbadcrc,ss,pberr diff --git a/lib/fsk4hf/wsprdpsksim.f90 b/lib/fsk4hf/wsprdpsksim.f90 new file mode 100644 index 000000000..80d294279 --- /dev/null +++ b/lib/fsk4hf/wsprdpsksim.f90 @@ -0,0 +1,175 @@ +program wsprdpsksim + +! Generate simulated data for a 2-minute "WSPR-DPSK" mode. Output is saved +! to a *.c2 or *.wav file. + + use wavhdr + include 'wsprdpsk_params.f90' !Set various constants + parameter (NMAX=120*12000) + type(hdr) hwav !Header for .wav file + character arg*12,fname*16 + character msg*22,msgsent*22 + complex c0(0:NMAX/NDOWN-1) + complex c(0:NMAX/NDOWN-1) + complex c0wav(0:NMAX-1) + complex cwav(0:NMAX-1) + real*8 fMHz + integer imessage(NN) + integer idemess(NN) !Differentially encoded message + integer*2 iwave(NMAX) !Generated full-length waveform + +! Get command-line argument(s) + nargs=iargc() + if(nargs.ne.8) then + print*,'Usage: wsprdpsksim "message" f0 DT fsp del nwav nfiles snr' + print*,'Example: wsprdpsksim "K1ABC FN42 30" 50 0.0 0.1 1.0 1 10 -33' + go to 999 + endif + call getarg(1,msg) !Message to be transmitted + call getarg(2,arg) + read(arg,*) f0 !Freq relative to WSPR-band center (Hz) + call getarg(3,arg) + read(arg,*) xdt !Time offset from nominal (s) + call getarg(4,arg) + read(arg,*) fspread !Watterson frequency spread (Hz) + call getarg(5,arg) + read(arg,*) delay !Watterson delay (ms) + call getarg(6,arg) + read(arg,*) nwav !1 for *.wav file, 0 for *.c2 file + call getarg(7,arg) + read(arg,*) nfiles !Number of files + call getarg(8,arg) + read(arg,*) snrdb !SNR_2500 + + twopi=8.0*atan(1.0) + pi=twopi/2.0 + fs=12000.0/NDOWN + dt=1.0/fs !Sample interval (s) + tt=NSPS*dt !Duration of "itone" symbols (s) + baud=1.0/tt !Keying rate for "itone" symbols (baud) + txt=NZ*dt !Transmission length (s) + bandwidth_ratio=2500.0/(fs/2.0) + sig=sqrt(bandwidth_ratio) * 10.0**(0.05*snrdb) + if(snrdb.gt.90.0) sig=1.0 + txt=NN*NSPS0/12000.0 + + call genwsprdpsk(msg,msgsent,imessage) !Encode the message, get itone + write(*,1000) f0,xdt,txt,snrdb,fspread,delay,nfiles,msgsent +1000 format('f0:',f9.3,' DT:',f6.2,' txt:',f6.1,' SNR:',f6.1, & + ' fspread:',f6.1,' delay:',f6.1,' nfiles:',i3,2x,a22) + +! Assume 0'th message bit is 0 + imessage=2*imessage-1 + idemess(1)=-1*imessage(1) + do i=2,NN + idemess(i)=imessage(i-1)*imessage(i) + enddo + + beta=1.0 ! excess bandwidth + if(nwav.eq.0) then + df=fs/(NMAX/NDOWN) ! + c=0 + bw=(1+beta)*baud/2.0 + bf=(1-beta)*baud/2.0 + iw=bw/df + if=bf/df + c(0:if-1)=1.0 + if(iw.gt.if) then + do i=if,iw + c(i)=((1.0+cos(pi*(i-if)/(iw-if)))/2.0)**0.5 + enddo + endif + c(NMAX/NDOWN-1:NMAX/NDOWN-iw:-1)=c(1:iw) + + istart=xdt/dt + c0=0.0 + do i=1,NN + c0(istart+(i-1)*200)=idemess(i) + enddo + call four2a(c0,NMAX/NDOWN,1,1,1) + c0=c0*conjg(c) + ic=f0/df + c0=cshift(c0,-ic) + call four2a(c0,NMAX/NDOWN,1,-1,1) + xx=sum(abs(c0(istart:istart+(NN-1)*200)**2))/((NN-1)*200) + c0=c0/sqrt(xx) + + call sgran() + do ifile=1,nfiles + c=c0 + if( fspread .ne. 0.0 .or. delay .ne. 0.0 ) then + call watterson(c,NMAX/NDOWN,fs,delay,fspread) + endif + c=c*sig + if(snrdb.lt.90) then + do i=0,NMAX/NDOWN-1 !Add gaussian noise at specified SNR + xnoise=gran() + ynoise=gran() + c(i)=c(i) + cmplx(xnoise,ynoise) + enddo + endif + write(fname,1100) ifile +1100 format('000000_',i4.4,'.c2') + open(10,file=fname,status='unknown',access='stream') + fMHz=10.1387d0 + nmin=2 + write(10) fname,nmin,fMHz,c !Save to *.c2 file + close(10) + enddo + else + fs=12000.0 + df=fs/NMAX + dt=1/fs + bandwidth_ratio=2500.0/(fs/2.0) + sig=sqrt(bandwidth_ratio) * 10.0**(0.05*snrdb) + if(snrdb.gt.90.0) sig=1.0 + cwav=0 + bw=(1+beta)*baud/2.0 + bf=(1-beta)*baud/2.0 + iw=bw/df + if=bf/df + cwav(0:if-1)=1.0 + if(iw.gt.if) then + do i=if,iw + cwav(i)=((1.0+cos(pi*(i-if)/(iw-if)))/2.0)**0.5 + enddo + endif + cwav(NMAX-1:NMAX-iw:-1)=cwav(1:iw) + + istart=xdt/dt + c0wav=0.0 + do i=1,NN + c0wav(istart+(i-1)*200*NDOWN)=idemess(i) + enddo + call four2a(c0wav,NMAX,1,1,1) + c0wav=c0wav*conjg(cwav) + ic=f0/df + c0wav=cshift(c0wav,-ic) + call four2a(c0wav,NMAX,1,-1,1) + xx=sum(abs(c0wav(istart:istart+(NN-1)*200*NDOWN))**2)/((NN-1)*200*NDOWN) + c0wav=c0wav/sqrt(xx) + call sgran() + do ifile=1,nfiles + cwav=c0wav + if( fspread .ne. 0.0 .or. delay .ne. 0.0 ) then + call watterson(cwav,NMAX,fs,delay,fspread) + endif + cwav=cwav*sig + if(snrdb.lt.90) then + do i=1,NMAX !Add gaussian noise at specified SNR + xnoise=gran() + iwave(i)=100*(real(cwav(i-1)) + xnoise) + enddo + endif + hwav=default_header(12000,NMAX) + write(fname,1102) ifile +1102 format('000000_',i4.4,'.wav') + open(10,file=fname,status='unknown',access='stream') + write(10) hwav,iwave !Save to *.wav file + close(10) + enddo + endif + write(*,1110) ifile,xdt,f0,snrdb,fname +1110 format(i4,f7.2,f8.2,f7.1,2x,a16) + +999 end program wsprdpsksim