2017-11-30 15:34:05 -05:00
|
|
|
subroutine foxgen()
|
2017-11-16 16:32:07 -05:00
|
|
|
|
2017-12-08 12:03:11 -05:00
|
|
|
! Called from MainWindow::foxTxSequencer() to generate the Tx waveform in
|
|
|
|
|
! FT8 Fox mode. The Tx message can contain up to 5 "slots", each carrying
|
|
|
|
|
! its own FT8 signal.
|
|
|
|
|
|
|
|
|
|
! Encoded messages can be of the form "HoundCall FoxCall rpt" (a standard FT8
|
|
|
|
|
! message with i3bit=0) or "HoundCall_1 RR73; HoundCall_2 <FoxCall> rpt",
|
|
|
|
|
! a new message type with i3bit=1. The waveform is generated with
|
|
|
|
|
! fsample=48000 Hz; it is compressed to reduce the PEP-to-average power ratio,
|
|
|
|
|
! with (currently disabled) filtering afterware to reduce spectral growth.
|
|
|
|
|
|
|
|
|
|
! Input message information is provided in character array cmsg(5), in
|
|
|
|
|
! common/foxcom/. The generated wave(NWAVE) is passed back in the same
|
|
|
|
|
! common block.
|
|
|
|
|
|
2018-09-30 13:02:14 -04:00
|
|
|
parameter (NN=79,ND=58,NSPS=4*1920)
|
2020-07-23 14:13:05 -04:00
|
|
|
parameter (NWAVE=(160+2)*134400*4) !the biggest waveform we generate (FST4-1800 at 48kHz)
|
2020-07-17 18:44:14 -04:00
|
|
|
parameter (NFFT=614400,NH=NFFT/2)
|
2017-12-21 15:55:38 -05:00
|
|
|
character*40 cmsg
|
2018-09-13 12:39:08 -04:00
|
|
|
character*37 msg,msgsent
|
2018-09-30 13:02:14 -04:00
|
|
|
integer itone(79)
|
|
|
|
|
integer*1 msgbits(77),msgbits2
|
2017-12-07 16:36:25 -05:00
|
|
|
integer*1, target:: mycall
|
2017-12-12 12:59:09 -05:00
|
|
|
real x(NFFT)
|
2017-11-16 16:32:07 -05:00
|
|
|
real*8 dt,twopi,f0,fstep,dfreq,phi,dphi
|
2017-12-12 12:59:09 -05:00
|
|
|
complex cx(0:NH)
|
2017-12-27 16:06:57 -05:00
|
|
|
common/foxcom/wave(NWAVE),nslots,nfreq,i3bit(5),cmsg(5),mycall(12)
|
2018-09-30 13:02:14 -04:00
|
|
|
common/foxcom2/itone2(NN),msgbits2(77)
|
2017-11-16 16:32:07 -05:00
|
|
|
equivalence (x,cx),(y,cy)
|
|
|
|
|
|
|
|
|
|
fstep=60.d0
|
|
|
|
|
dfreq=6.25d0
|
|
|
|
|
dt=1.d0/48000.d0
|
|
|
|
|
twopi=8.d0*atan(1.d0)
|
2017-12-08 12:03:11 -05:00
|
|
|
irpt=0
|
2017-12-12 12:59:09 -05:00
|
|
|
nplot=0
|
|
|
|
|
wave=0.
|
2017-11-16 16:32:07 -05:00
|
|
|
|
2017-11-27 15:22:44 -05:00
|
|
|
do n=1,nslots
|
2018-09-30 13:02:14 -04:00
|
|
|
msg=cmsg(n)(1:37)
|
2018-11-05 11:59:48 -05:00
|
|
|
call genft8(msg,i3,n3,msgsent,msgbits,itone)
|
2017-12-08 12:03:11 -05:00
|
|
|
! Make copies of itone() and msgbits() for ft8sim
|
|
|
|
|
itone2=itone
|
|
|
|
|
msgbits2=msgbits
|
2017-12-27 16:06:57 -05:00
|
|
|
f0=nfreq + fstep*(n-1)
|
2017-11-16 16:32:07 -05:00
|
|
|
phi=0.d0
|
|
|
|
|
k=0
|
|
|
|
|
do j=1,NN
|
|
|
|
|
f=f0 + dfreq*itone(j)
|
|
|
|
|
dphi=twopi*f*dt
|
|
|
|
|
do ii=1,NSPS
|
|
|
|
|
k=k+1
|
|
|
|
|
phi=phi+dphi
|
|
|
|
|
xphi=phi
|
|
|
|
|
wave(k)=wave(k)+sin(xphi)
|
|
|
|
|
enddo
|
|
|
|
|
enddo
|
|
|
|
|
enddo
|
2017-12-12 12:59:09 -05:00
|
|
|
kz=k
|
2017-11-21 16:00:54 -05:00
|
|
|
|
2017-12-12 12:59:09 -05:00
|
|
|
peak1=maxval(abs(wave))
|
|
|
|
|
wave=wave/peak1
|
2017-12-12 15:15:36 -05:00
|
|
|
width=50.0
|
2017-12-27 16:06:57 -05:00
|
|
|
call foxfilt(nslots,nfreq,width,wave)
|
2017-12-12 12:59:09 -05:00
|
|
|
peak3=maxval(abs(wave))
|
|
|
|
|
wave=wave/peak3
|
|
|
|
|
|
2017-11-16 16:32:07 -05:00
|
|
|
return
|
|
|
|
|
end subroutine foxgen
|
2017-12-12 12:59:09 -05:00
|
|
|
|
2017-12-13 12:18:34 -05:00
|
|
|
! include 'plotspec.f90'
|
