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

Fixed the frequency problem.

Browse Source 
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@2532 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Joe Taylor 2012-08-28 14:36:29 +00:00
parent 6bfe2ed709
commit 37d91076fc
1 changed files with 11 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
libm65/mapsim.f90
View File

@ -34,23 +34,23 @@ program mapsim
b=sin(pol/rad) b=sin(pol/rad)
dt=1.d0/fsample dt=1.d0/fsample
do isig=10,10 do isig=1,20
f=-23 + 3*isig f=-23000 + 3000*(isig-1)
dt=0.05d0*(isig-1) dphi=twopi*f*dt + 0.5*twopi
! snrdb=-(19.0 + (isig-1)/2.0) ! snrdb=-(19.0 + (isig-1)/2.0)
snrdb=-20.0 snrdb=-20.0
sig=sqrt(2500.0/96000.0) * 10.0**(0.05*snrdb) sig=sqrt(2.0*2500.0/96000.0) * 10.0**(0.05*snrdb)
sig=1.6*sig
dphi=twopi*f*dt
phi=0. phi=0.
i0=fsample*(3.5d0+dt) i0=fsample*(3.5d0+0.05d0*(isig-1))
print*,f,dt,dphi,i0,sig
do i=1,nwave do i=1,nwave
phi=phi + dphi phi=phi + dphi
! if(phi.gt.twopi) phi=phi-twopi if(phi.lt.-twopi) phi=phi+twopi
! xphi=phi if(phi.gt.twopi) phi=phi-twopi
z=sig*cwave(i)*cmplx(cos(phi),sin(phi)) xphi=phi
z=sig*cwave(i)*cmplx(cos(xphi),-sin(xphi))
zx=a*z zx=a*z
zy=b*z zy=b*z
j=i+i0 j=i+i0