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

Add optional code to q65sim to confirm calkibration of SNR_2500.

Browse Source
This commit is contained in:
Joe Taylor 2024-02-19 13:52:37 -05:00
parent fb19d27602
commit d234986165
1 changed files with 12 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
lib/qra/q65/q65sim.f90
View File

@ -222,14 +222,20 @@ program q65sim
fac=sqrt(1.0/avep) fac=sqrt(1.0/avep)
cspread=fac*cspread !Normalize to constant avg power cspread=fac*cspread !Normalize to constant avg power
cdat=cspread*cdat !Apply Rayleigh fading cdat=cspread*cdat !Apply Rayleigh fading
! do i=0,nfft-1
! p=real(cspread(i))**2 + aimag(cspread(i))**2
! write(14,3010) i,p,cspread(i)
!3010 format(i8,3f12.6)
! enddo
endif endif
! psig=0.
! pnoise=0.
! do i=1,npts
! if(i.gt.12000 .and. i.lt.624000) then
! psig=psig + aimag(cdat(i))**2
! pnoise=pnoise + xnoise(i)*xnoise(i)
! endif
! enddo
! pnoise=pnoise*bandwidth_ratio
! snr_2500=db(psig/pnoise) !Calibration confirmation!
! print*,'SNR_2500:',snr_2500
dat=aimag(cdat) + xnoise !Add generated AWGN noise dat=aimag(cdat) + xnoise !Add generated AWGN noise
fac=32767.0 fac=32767.0
if(snrdb.ge.90.0) iwave(1:npts)=nint(fac*dat(1:npts)) if(snrdb.ge.90.0) iwave(1:npts)=nint(fac*dat(1:npts))