15bf364e5e
commit dbca224ee15f598ba86b860c19e902b594f6b888
Author: Jordan Sherer <jordan@widefido.com>
Date: Wed Aug 14 21:10:17 2019 -0400
Cleanup decoder with easy parameters
commit 427a6ac6bf02a7de80e76b921bbd48ec23dab0c4
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 22:56:32 2019 -0400
Working on downsampling
commit ef5e1fb9cdafb6069b6f67377c3bc9e073d2785e
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 22:01:34 2019 -0400
Fixed NN parameters
commit 216dbc83a931127540fb028d3365a205b4ef1d28
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 21:40:20 2019 -0400
Fixed dupe during sync
commit 9505d50ba6c42b05e5fff1a059864a643d9e7558
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 16:54:37 2019 -0400
Tweaking downsampling
commit 86cd07e1ac8ce8d64cdd2d4321cf295cd867f80e
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 10:11:07 2019 -0400
Working with logging
commit 0de3849a442d20081f545c6fef5fadf2def6e0d2
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 03:46:51 2019 -0400
250 millisecond startup delay
commit a35cc5302b6f8532a424d504a4e2feaacf74ef99
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 03:39:20 2019 -0400
1200 sps, 10 seconds, 80Hz
commit 6f36db20f90c510336b497605d9d7667d6f0f44f
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 03:27:12 2019 -0400
Working word count
commit 741c55ac0840512d68e1fbdb2f1fff6adfcd4612
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 02:40:40 2019 -0400
Updated LDPC simulator
commit ce0992644315402503f44bdae4131bacf2cc3e21
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 02:34:01 2019 -0400
Weird. Apparently symbol stop is the only way to get decodes to work in app...need to figure out why
commit 8b330fbff6574f2ee6d7a471dd93a2cb57b492dd
Author: Jordan Sherer <jordan@widefido.com>
Date: Tue Aug 13 01:50:37 2019 -0400
More parameter fixups (downsampling) and first decode
commit 37c780b45e6579294b5a4692b8e8d0602d9bd8b1
Author: Jordan Sherer <jordan@widefido.com>
Date: Mon Aug 12 23:20:58 2019 -0400
More params fixing
commit bc8547e1e8f55bfb4bce5f95890f3cf6798e36d7
Author: Jordan Sherer <jordan@widefido.com>
Date: Mon Aug 12 22:48:00 2019 -0400
Fixed params in fortran code
commit 10fbf9f4a78c6244e0082a31819ee26aeae39325
Author: Jordan Sherer <jordan@widefido.com>
Date: Mon Aug 12 16:08:32 2019 -0400
Using constants in place of literals in a few places to make it easier to change
52 lines
1.3 KiB
Fortran
52 lines
1.3 KiB
Fortran
subroutine baseline(s,nfa,nfb,sbase)
|
|
|
|
! Fit baseline to spectrum (for FT8)
|
|
! Input: s(npts) Linear scale in power
|
|
! Output: sbase(npts) Baseline
|
|
|
|
implicit real*8 (a-h,o-z)
|
|
|
|
include 'ft8_params.f90'
|
|
|
|
real*4 s(NSPS)
|
|
real*4 sbase(NSPS)
|
|
real*4 base
|
|
real*8 x(1000),y(1000),a(5)
|
|
data nseg/10/,npct/10/
|
|
|
|
df=12000.0/(NSPS*2.0d0) !3.125 Hz
|
|
ia=max(1,nint(nfa/df))
|
|
ib=nint(nfb/df)
|
|
do i=ia,ib
|
|
s(i)=10.0*log10(s(i)) !Convert to dB scale
|
|
enddo
|
|
|
|
nterms=5
|
|
nlen=(ib-ia+1)/nseg !Length of test segment
|
|
i0=(ib-ia+1)/2 !Midpoint
|
|
k=0
|
|
do n=1,nseg !Loop over all segments
|
|
ja=ia + (n-1)*nlen
|
|
jb=ja+nlen-1
|
|
call pctile(s(ja),nlen,npct,base) !Find lowest npct of points
|
|
do i=ja,jb
|
|
if(s(i).le.base) then
|
|
if (k.lt.1000) k=k+1 !Save all "lower envelope" points
|
|
x(k)=i-i0
|
|
y(k)=s(i)
|
|
endif
|
|
enddo
|
|
enddo
|
|
kz=k
|
|
a=0.
|
|
call polyfit(x,y,y,kz,nterms,0,a,chisqr) !Fit a low-order polynomial
|
|
do i=ia,ib
|
|
t=i-i0
|
|
sbase(i)=a(1)+t*(a(2)+t*(a(3)+t*(a(4)+t*(a(5))))) + 0.65
|
|
! write(51,3051) i*df,s(i),sbase(i)
|
|
!3051 format(3f12.3)
|
|
enddo
|
|
|
|
return
|
|
end subroutine baseline
|