mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2024-11-26 14:18:38 -05:00
a51c5c4251
As at ^/branches/wsjtx_exp@5271 git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@5272 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
29 lines
669 B
Fortran
29 lines
669 B
Fortran
subroutine wav11(d2,npts,dd)
|
|
|
|
! Convert i*2 data sampled at 12000 Hz to r*4 sampled at 11025 Hz.
|
|
|
|
parameter (NZ11=60*11025,NZ12=60*12000)
|
|
parameter (NFFT1=64*12000,NFFT2=64*11025)
|
|
integer*2 d2(NZ12)
|
|
real*4 dd(NZ11)
|
|
real x(NFFT2)
|
|
complex cx(0:NFFT1/2)
|
|
equivalence (x,cx)
|
|
save x,cx
|
|
|
|
jz=min(NZ12,npts)
|
|
x(1:jz)=d2(1:jz)
|
|
x(jz+1:)=0.0
|
|
call four2a(x,nfft1,1,-1,0) !Forwarxd FFT, r2c
|
|
df=12000.0/NFFT1
|
|
ia=5000.0/df
|
|
cx(ia:)=0.0
|
|
call four2a(cx,nfft2,1,1,-1) !Inverse FFT, c2r
|
|
npts=jz*11025.0/12000.0
|
|
fac=1.e-6
|
|
dd(1:npts)=fac*x(1:npts)
|
|
if(npts.lt.NZ11) dd(npts+1:NZ11)=0.0
|
|
|
|
return
|
|
end subroutine wav11
|