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

Merge branch 'feat-fst280' of bitbucket.org:k1jt/wsjtx into feat-fst280

Browse Source
This commit is contained in:
Bill Somerville 2020-07-15 01:48:01 +01:00
commit c5ab12d1f4
No known key found for this signature in database
GPG Key ID: D864B06D1E81618F
1 changed files with 25 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

30
lib/fst240_decode.f90
View File

@ -64,7 +64,7 @@ contains
integer mcq(29),mrrr(19),m73(19),mrr73(19) integer mcq(29),mrrr(19),m73(19),mrr73(19)
logical badsync,unpk77_success,single_decode logical badsync,unpk77_success,single_decode
logical first,nohiscall,lwspr logical first,nohiscall,lwspr,ex
integer*2 iwave(300*12000) integer*2 iwave(300*12000)
@ -547,7 +547,8 @@ contains
else else
call get_fst240_tones_from_bits(message74,itone,1) call get_fst240_tones_from_bits(message74,itone,1)
endif endif
if(.false.) then inquire(file='plotspec',exist=ex)
if(ex) then
call write_ref(itone,iwave,nsps,nmax,ndown,hmod, & call write_ref(itone,iwave,nsps,nmax,ndown,hmod, &
isbest,fc_synced) isbest,fc_synced)
endif endif
@ -808,6 +809,7 @@ contains
subroutine write_ref(itone,iwave,nsps,nmax,ndown,hmod,i0,fc) subroutine write_ref(itone,iwave,nsps,nmax,ndown,hmod,i0,fc)
complex cwave(nmax) complex cwave(nmax)
complex c(0:1440000-1)
integer itone(160) integer itone(160)
integer*2 iwave(nmax) integer*2 iwave(nmax)
integer hmod integer hmod
@ -817,10 +819,28 @@ contains
nsym=160 nsym=160
call gen_fst240wave(itone,nsym,nsps,nmax,fsample,hmod,fc, & call gen_fst240wave(itone,nsym,nsps,nmax,fsample,hmod,fc, &
1,cwave,wave) 1,cwave,wave)
cwave=cshift(cwave,-i0*ndown) cwave=cshift(cwave,-i0*ndown)
do i=1,nmax ! do i=1,nmax
write(51,*) i,iwave(i),real(cwave(i)),imag(cwave(i)) ! write(51,1000) i,iwave(i),cwave(i)
!1000 format(2i10,f12.6)
! enddo
fac=1.0/32768
c=fac*float(iwave)*conjg(cwave)
call four2a(c,nmax,1,-1,1) !Forward c2c FFT
df=12000.0/nmax
ia=-10.0/df
ib=10.0/df
do i=ia,ib
j=i
if(j.lt.0) j=i+nmax
s=real(c(j))**2 + aimag(c(j))**2
f=i*df
write(52,1010) f,s,db(s)
1010 format(f10.3,e12.3,f10.3)
enddo enddo
return
end subroutine write_ref end subroutine write_ref
end module fst240_decode end module fst240_decode