w4kek
/
WSJT-X
mirror of https://github.com/saitohirga/WSJT-X.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Revised figures for ftrsd paper. 

git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6267 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Steven Franke 2015-12-14 02:09:44 +00:00
parent 81e11c06ad
commit 01b65be2f1
10 changed files with 52 additions and 51 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
lib/ftrsd/ftrsd_paper/bmdata-2.dat
View File

@ -1,8 +0,0 @@
fspread0 0.2, BM (ntrials=0)
-21.5 0.05
-21.0 0.17
-20.5 0.34
-20.0 0.60
-19.5 0.79
-19.0 0.94
-18.5 0.984

15
lib/ftrsd/ftrsd_paper/fig_psuccess.gnuplot
View File

@ -2,7 +2,7 @@
# run: gnuplot fig_psuccess.gnuplot
# then: pdflatex fig_psuccess.tex
#
set term epslatex standalone size 12cm,8cm
set term epslatex standalone size 20cm,10cm
set output "fig_psuccess.tex"
set xlabel "SNR in 2500 Hz BW (dB)"
set ylabel "Percent copy"
@ -12,7 +12,12 @@ set tics in
set mxtics 2
set mytics 2
set grid
plot "ftdata-10000-tight.dat" using 1:2 every ::1 with linespoints pt 5 title 'FT-1', \
"ftdata-100000-loose.dat" using 1:2 every ::1 with linespoints pt 5 title 'FT-2', \
"bmdata.dat" using 1:2 with linespoints pt 7 title 'BM'
# "mfsk64.dat" using ($1-29.7):(1-$2) every ::2 with linespoints pt 9 title 'symbols'
plot "ftdata-10000.dat" using 1:2 every ::1 with linespoints pt 4 title 'FT-10K', \
"ftdata-100000.dat" using 1:2 every ::1 with linespoints pt 5 title 'FT-100K', \
"ftdata-1000-rf.dat" using 1:2 every ::1 with linespoints pt 7 title 'FT-1K-RF', \
"ftdata-100-rf.dat" using 1:2 every ::1 with linespoints pt 8 title 'FT-100-RF', \
"bmdata.dat" using 1:2 with linespoints pt 9 title 'BM', \
"bmdata-rf.dat" using 1:2 with linespoints pt 10 title 'BM-RF', \
"kvasd-7999-rf.dat" using 1:2 every ::1 with linespoints pt 12 title 'KV-8-RF', \
"kvasd-11999.dat" using 1:2 every ::1 with linespoints pt 13 title 'KV-12', \
"kvasd-11999-rf.dat" using 1:2 every ::1 with linespoints pt 14 title 'KV-12-RF'

BIN
lib/ftrsd/ftrsd_paper/fig_psuccess.pdf
View File

Binary file not shown.

5
lib/ftrsd/ftrsd_paper/fig_wer.gnuplot
View File

@ -13,6 +13,7 @@ set mxtics 2
set mytics 10
set grid
set logscale y
plot "ftdata-10000-tight.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 5 title 'FT-1', \
"ftdata-100000-loose.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 5 title 'FT-2', \
plot "ftdata-10000.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 5 title 'FT-10K', \
"ftdata-100000.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 7 title 'FT-100K', \
"kvasd-11999.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 8 title 'KV-11.999', \
"bmdata.dat" using ($1+29.7):(1-$2) with linespoints pt 7 title 'BM'

22
lib/ftrsd/ftrsd_paper/fig_wer2.gnuplot
View File

@ -2,16 +2,30 @@
# run: gnuplot fig_wer2.gnuplot
# then: pdflatex fig_wer2.tex
#
set term epslatex standalone size 12cm,8cm
set term epslatex standalone size 16cm,10cm
set output "fig_wer2.tex"
set xlabel "$E_s/N_o$ (dB)"
set x2label "SNR in 2500 Hz (dB)"
set ylabel "WER"
#set autoscale xfix
#set autoscale x2fix
set style func linespoints
set key on top outside nobox
set tics in
set mxtics 2
set mytics 10
set grid
set grid ytics
set logscale y
plot "ftdata-2.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 7 title "FT", \
"bmdata-2.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 5 title 'BM'
set x2tics out
set xtics nomirror
set mx2tics 2
#set xrange [3:13]
#set x2range [(3-29.7):(13-29.7)]
set xrange [4:7.5]
set x2range [(4-29.7):(7.5-29.7)]
set yrange [0.001:1.0]
plot "ftdata-1000-rf.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 4 title "FT-1K-RF", \
"ftdata-100-rf.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 5 title 'FT-100-RF', \
"kvasd-7999-rf.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 6 title 'KV-8-RF', \
"kvasd-11999-rf.dat" using 1:(1-$2) every ::1 with linespoints pt 7 title 'KV-12-RF' axes x2y1
#"bmdata-rf.dat" using 1:(1-$2) every ::1 with linespoints pt 8 title 'BM-RF' axes x2y1

18
lib/ftrsd/ftrsd_paper/fig_wer3.gnuplot Normal file
View File

@ -0,0 +1,18 @@
# gnuplot script for AWGN vs Rayleigh figure
#
set term epslatex standalone size 16cm,8cm
set output "fig_wer3.tex"
set xlabel "$E_s/N_o$ (dB)"
set ylabel "WER"
set style func linespoints
set key on top outside nobox
set tics in
set mxtics 2
set mytics 10
set grid
set logscale y
#set format y "10^{%L}"
plot "ftdata-1000-rf.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 7 title "FT-1K-RF", \
"bmdata-rf.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 5 title 'BM-RF', \
"ftdata-10000.dat" using ($1+29.7):(1-$2) every ::1 with linespoints pt 7 title 'FT-10K-AWGN', \
"bmdata.dat" using ($1+29.7):(1-$2) with linespoints pt 5 title 'BM-AWGN'

10
lib/ftrsd/ftrsd_paper/ftdata-10000-tight.dat
View File

@ -1,10 +0,0 @@
snr psuccess ntrials 10000, ntotal<72, nhard<42
-26.5 0.0
-26.0 0.016
-25.5 0.090
-25.0 0.281
-24.5 0.581
-24.0 0.844
-23.5 0.977
-23.0 0.998

11
lib/ftrsd/ftrsd_paper/ftdata-100000-loose.dat
View File

@ -1,11 +0,0 @@
snr psuccess 100000 trials, ntotal<76, nhard<44
-27.0 0.0
-26.5 0.008
-26.0 0.04
-25.5 0.2
-25.0 0.41
-24.5 0.75
-24.0 0.93700
-23.5 0.99000
-23.0 0.99954

8
lib/ftrsd/ftrsd_paper/ftdata-2.dat
View File

@ -1,8 +0,0 @@
fspread0 0.2, FT (ntrials=1000)
-26.0 0.020
-25.5 0.09
-25.0 0.23
-24.5 0.58
-24.0 0.83
-23.5 0.94
-23.0 0.983

6
lib/ftrsd/ftrsd_paper/mfsk.f90
View File

@ -9,8 +9,8 @@ program mfsk
nargs=iargc()
if(nargs.ne.1) then
print*,'Probability of symbol error for noncoherent MFSK'
print*,'Usage: prob M'
print*,'Example: prob 64'
print*,'Usage: mfsk M'
print*,'Example: mfsk 64'
go to 999
endif
call getarg(1,arg)
@ -19,7 +19,7 @@ program mfsk
write(*,1012)
1012 format('Es/No P(symbol error)'/ &
'----------------------')
do isnr=0,20
do isnr=0,40
esno=10**(isnr/2.0/10.0)
hsum=0.d0
do k=1,M-1