From 35fd23c39278b196f3dc194ba5e30777a1c6e484 Mon Sep 17 00:00:00 2001 From: Steven Franke Date: Mon, 1 Feb 2016 03:01:04 +0000 Subject: [PATCH] Some more work on code comparisons. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6454 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 --- lib/JTMSKfer.f90 | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/lib/JTMSKfer.f90 b/lib/JTMSKfer.f90 index 04c638435..6c8642021 100644 --- a/lib/JTMSKfer.f90 +++ b/lib/JTMSKfer.f90 @@ -33,9 +33,10 @@ program JTMSKfer bias=0.0 scale=20.0 xln2=log(2.0) + mettab=0 do i=128,156 - x0=log(max(0.0001,2.0*xp(i-127)))/xln2 - x1=log(max(0.001,2.0*(1.0-xp(i-127))))/xln2 + x0=log(max(0.001,2.0*xp(i-127)))/xln2 + x1=log(max(0.001,2.0*(1-xp(i-127))))/xln2 mettab(i,0)=nint(scale*(x0-bias)) mettab(i,1)=nint(scale*(x1-bias)) mettab(256-i,0)=mettab(i,1) @@ -48,7 +49,7 @@ program JTMSKfer mettab(256-i,1)=mettab(i,0) enddo - rdscale=2.0 ! empirically optimized + rdscale=2.0 ntrials=1000000 rate=72.0/198.0 msg="123" @@ -67,7 +68,7 @@ program JTMSKfer ! call sgran() - do idb=0,11 + do idb=6,11 db=idb/2.0-0.5 ! Eb/N0=1/(2*R*sigma^2), so sigma= sqrt( 1/(2*R*Eb/N0) ) sigma=1/sqrt( 2*rate*(10**(db/10.0)) )