Fix the sign of one of the Nuttal window coefficients. Change BT to 2.0.

This commit is contained in:
Steve Franke 2019-02-25 12:57:34 -06:00
parent 17ac99344c
commit ec2d5b6233
2 changed files with 4 additions and 4 deletions

View File

@ -17,7 +17,7 @@ subroutine gen_ft8wave(itone,nsym,nsps,fsample,f0,cwave,wave,icmplx,nwave)
twopi=8.0*atan(1.0)
dt=1.0/fsample
hmod=1.0
bt=1.5
bt=2.0
! Compute the frequency-smoothing pulse
do i=1,3*nsps
tt=(i-1.5*nsps)/real(nsps)

View File

@ -58,11 +58,11 @@ subroutine symspec(shared_data,k,ntrperiod,nsps,ingain,nminw,pxdb,s, &
! discrete Hilbert Transform, by Albert Nuttall"
!
a0=0.3635819
a1=0.4891775;
a1=-0.4891775;
a2=0.1365995;
a3=0.0106411;
a3=-0.0106411;
do i=1,nfft3
w3(i)=a0-a1*cos(2*pi*(i-1)/(nfft3))+ &
w3(i)=a0+a1*cos(2*pi*(i-1)/(nfft3))+ &
a2*cos(4*pi*(i-1)/(nfft3))+ &
a3*cos(6*pi*(i-1)/(nfft3))
enddo