mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2024-11-03 07:51:16 -05:00
612f13a0e7
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@357 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
72 lines
2.1 KiB
Fortran
72 lines
2.1 KiB
Fortran
subroutine symspec(id,nz,ss,savg)
|
|
|
|
C Compute spectra at four polarizations, using half-symbol steps.
|
|
|
|
parameter (NFFT=32768)
|
|
integer*2 id(4,nz)
|
|
real ss(5,322,NFFT)
|
|
real savg(4,NFFT)
|
|
complex cx(NFFT),cy(NFFT) ! pad to 32k with zeros
|
|
complex z
|
|
real*8 ts,hsym
|
|
|
|
fac=1.e-4
|
|
hsym=2048.d0*96000.d0/11025.d0 !Samples per half symbol
|
|
npts=hsym !Integral samples per half symbol
|
|
nsteps=322 !Half symbols per transmission
|
|
|
|
do ip=1,4
|
|
do i=1,NFFT
|
|
savg(ip,i)=0.
|
|
enddo
|
|
enddo
|
|
|
|
ts=1.d0 - hsym
|
|
do n=1,nsteps
|
|
ts=ts+hsym !Update exact sample pointer
|
|
i0=ts !Starting sample pointer
|
|
do i=1,npts !Copy data to FFT arrays
|
|
xr=fac*id(1,i0+i)
|
|
xi=fac*id(2,i0+i)
|
|
cx(i)=cmplx(xr,xi)
|
|
yr=fac*id(3,i0+i)
|
|
yi=fac*id(4,i0+i)
|
|
cy(i)=cmplx(yr,yi)
|
|
enddo
|
|
|
|
do i=npts+1,NFFT !Pad to 32k with zeros
|
|
cx(i)=0.
|
|
cy(i)=0.
|
|
enddo
|
|
|
|
call four2a(cx,NFFT,1,1,1) !Do the FFTs
|
|
call four2a(cy,NFFT,1,1,1)
|
|
|
|
do i=1,NFFT !Save and accumulate power spectra
|
|
s1=real(cx(i))**2 + aimag(cx(i))**2
|
|
ss(1,n,i)=s1 ! Pol = 0
|
|
savg(1,i)=savg(1,i) + s1
|
|
|
|
z=cx(i) + cy(i)
|
|
s2=0.5*(real(z)**2 + aimag(z)**2)
|
|
ss(2,n,i)=s2 ! Pol = 45
|
|
savg(2,i)=savg(2,i) + s2
|
|
|
|
s3=real(cy(i))**2 + aimag(cy(i))**2
|
|
ss(3,n,i)=s3 ! Pol = 90
|
|
savg(3,i)=savg(3,i) + s3
|
|
|
|
z=cx(i) - cy(i)
|
|
s4=0.5*(real(z)**2 + aimag(z)**2)
|
|
ss(4,n,i)=s4 ! Pol = 135
|
|
savg(4,i)=savg(4,i) + s4
|
|
|
|
z=cx(i) * conjg(cy(i))
|
|
ss(5,n,i)=0.5*(s1+s3) + (real(z)**2 + aimag(z)**2 -
|
|
+ s1*s3)/(s1+s3)
|
|
enddo
|
|
enddo
|
|
|
|
return
|
|
end
|