Skip to content

Commit 759311a

Browse files
thiagotpsthiagotps
committed
The scripts used for getting some of the published results were added.
1 parent 356d2c5 commit 759311a

File tree

4 files changed

+179
-0
lines changed

4 files changed

+179
-0
lines changed

‎experiments/mse-evolution/extra.sh

+52
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,52 @@
1+
#!/usr/bin/env bash
2+
3+
set -euo pipefail
4+
5+
export NITER=175
6+
export N=2
7+
export M=3
8+
export SIGMAV2=0.001
9+
export EXP=6
10+
export LANG=en_US
11+
export DIST="gauss"
12+
export DIR_PREFIX="data_${DIST}"
13+
14+
foo()
15+
{
16+
local beta=$1
17+
local mse_eea_cache="mse_eea_cache_${N}_${M}.bin"
18+
local DIR="${DIR_PREFIX}_${N}_${M}/$beta"
19+
local mse_ia_file="${DIR}/mse_ia_${N}_${M}.txt"
20+
local mse_st_file="${DIR}/mse_st_${N}_${M}.txt"
21+
local mse_eea_file="${DIR}/mse_eea_${N}_${M}.txt"
22+
local mse_png_file="${DIR}/mse_${N}_${M}.png"
23+
local mse_eps_file="${DIR}/mse_${N}_${M}.eps"
24+
25+
mkdir -p "$DIR"
26+
27+
# Por questões de desempenho, guardamos todas as recursões em um arquivo.
28+
# Se esse arquivo não existir, criamos ele.
29+
[[ ! -f "$mse_eea_cache" ]] && skewness --writecache "$mse_eea_cache" -N $N -M $M --indmode eea --outmode mse
30+
31+
# Obtemos o MSE usando EEA
32+
skewness --readcache "$mse_eea_cache" -N "$N" -M "$M" --indmode eea -n "$NITER" -b "$beta" --sv2 \
33+
"$SIGMAV2" --outmode mse -d "$DIST" -o "$mse_eea_file"
34+
35+
# Obtemos o MSE usando IA
36+
skewness -N "$N" -M "$M" --indmode ia -n "$NITER" -b "$beta" --sv2 "$SIGMAV2" --outmode mse -d \
37+
"$DIST" -o "$mse_ia_file"
38+
39+
# Obtemos o MSE usando ensaios de Monte Carlo
40+
skemp --filter-length $N --data-length $M --niter $NITER \
41+
--exp "$EXP" --beta "$beta" --sigmav2 "$SIGMAV2" --dist "$DIST" \
42+
--mse-file "$mse_st_file" \
43+
44+
echo "Creating the graphs"
45+
46+
# Finalmente criamos os gráficos
47+
sgg --yt db --xmin 0 --xmax "$NITER" -y "MSE (dB)" -s ":" " -" " --" -l "IA" "EEA" "Empírico" \
48+
-f "$mse_ia_file" "$mse_eea_file" "$mse_st_file" -d "$mse_png_file" --eps-file "$mse_eps_file" \
49+
-c "lightseagreen" "darkorange" "royalblue"
50+
}
51+
52+
foo "0.04"

‎experiments/skewness-evolution/experiment.sh

+84
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,84 @@
1+
#!/usr/bin/env bash
2+
3+
set -euo pipefail
4+
5+
print_options()
6+
{
7+
echo "DIST=$DIST"
8+
echo "NITER=$NITER"
9+
echo "N=$N"
10+
echo "M=$M"
11+
echo "SIGMAV2=$SIGMAV2"
12+
echo "EXP=$EXP"
13+
echo "LANG=$LANG"
14+
echo "BETA=$BETA"
15+
echo "K=$K"
16+
echo "PDF_START=$PDF_START"
17+
echo "PDF_END=$PDF_END"
18+
echo "PDF_SAMPLES=$PDF_SAMPLES"
19+
echo "KERNEL_EXP=$KERNEL_EXP"
20+
}
21+
22+
experiment()
23+
{
24+
DIR="data/${DIST}"
25+
sk_ia_file="${DIR}/sk_ia.txt"
26+
sk_eea_file="${DIR}/sk_eea.txt"
27+
sk_st_file="${DIR}/sk_st.txt"
28+
sk_eps_file="${DIR}/sk.eps"
29+
30+
mse_ia_file="${DIR}/mse_ia.txt"
31+
mse_eea_file="${DIR}/mse_eea.txt"
32+
mse_st_file="${DIR}/mse_st.txt"
33+
mse_eps_file="${DIR}/mse.eps"
34+
35+
pdf_file="${DIR}/pdf_st.txt"
36+
pdf_gauss_file="${DIR}/pdf_gauss.txt"
37+
pdf_eps_file="${DIR}/pdf.eps"
38+
39+
sk_eea_cache="sk_eea_cache_${N}_${M}.bin"
40+
mse_eea_cache="mse_eea_cache_${N}_${M}.bin"
41+
42+
mkdir -p "$DIR"
43+
44+
print_options > "${DIR}/README"
45+
46+
[[ ! -f "$sk_eea_cache" ]] && skewness --writecache "$sk_eea_cache" -N $N -M $M --indmode eea --outmode sk
47+
skewness --readcache "$sk_eea_cache" -N "$N" -M "$M" --indmode eea -n "$NITER" -b "$BETA" --sv2 "$SIGMAV2" --outmode sk -d "$DIST" -o "$sk_eea_file"
48+
skewness -N "$N" -M "$M" --indmode ia -n "$NITER" -b "$BETA" --sv2 "$SIGMAV2" --outmode sk -d "$DIST" -o "$sk_ia_file"
49+
50+
51+
[[ ! -f "$mse_eea_cache" ]] && skewness --writecache "$mse_eea_cache" -N $N -M $M --indmode eea --outmode mse
52+
skewness --readcache "$mse_eea_cache" -N "$N" -M "$M" --indmode eea -n "$NITER" -b "$BETA" --sv2 "$SIGMAV2" --outmode mse -d "$DIST" -o "$mse_eea_file"
53+
skewness -N "$N" -M "$M" --indmode ia -n "$NITER" -b "$BETA" --sv2 "$SIGMAV2" --outmode mse -d "$DIST" -o "$mse_ia_file"
54+
55+
56+
skemp --filter-length "$N" --data-length "$M" --niter "$NITER" \
57+
--exp "$EXP" --beta "$BETA" --sigmav2 "$SIGMAV2" --skewness-file "$sk_st_file" \
58+
--mse-file "$mse_st_file" \
59+
--pdf-file "$pdf_file" --pdf-instant "$K" --pdf-start "$PDF_START" --pdf-end "$PDF_END" --kernel-exp "$KERNEL_EXP" \
60+
--pdf-samples "$PDF_SAMPLES" \
61+
--dist "$DIST"
62+
63+
mean_sd=$(awk -F '#' 'FNR == 2 {print $2}' $pdf_file)
64+
mean=$(echo "$mean_sd" | cut -d ',' -f1)
65+
sd=$(echo "$mean_sd" | cut -d ',' -f2)
66+
67+
printf "mean = %s, σ = %s\n" "$mean" "$sd"
68+
69+
gauss_emp --mean "$mean" --sd "$sd" --start "$PDF_START" --end "$PDF_END" --samples "$PDF_SAMPLES" \
70+
--exp 6 --kernel-exp "$KERNEL_EXP" > "$pdf_gauss_file"
71+
72+
sgg -y "Skewness" -s ":" " -" " --" -l "IA" "EEA" "Empírico" \
73+
-f "$sk_ia_file" "$sk_eea_file" "$sk_st_file" --eps-file "$sk_eps_file" \
74+
-c "lightseagreen" "darkorange" "royalblue"
75+
76+
sgg -y "MSE" -s ":" " -" " --" -l "IA" "EEA" "Empírico" \
77+
-f "$mse_ia_file" "$mse_eea_file" "$mse_st_file" --eps-file "$mse_eps_file" \
78+
-c "lightseagreen" "darkorange" "royalblue" --yt "db" --yt-files "$mse_ia_file" "$mse_eea_file" "$mse_st_file"
79+
80+
sgg -y "Função Densidade de Probabilidade" -s " -" " --" -l "Gaussiano" "Empírico" \
81+
-f $pdf_gauss_file $pdf_file --eps-file $pdf_eps_file \
82+
-c "darkorange" "royalblue"
83+
84+
}

‎experiments/skewness-evolution/experiments_gaussian.sh

+22
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,22 @@
1+
#!/usr/bin/env bash
2+
# GAUSSIAN !!!
3+
4+
set -euo pipefail
5+
6+
export NITER=120
7+
export N=2
8+
export M=2
9+
export SIGMAV2=0.001
10+
export EXP=11
11+
export LANG=en_US
12+
export BETA=0.09
13+
export K=16
14+
export PDF_START=-0.5
15+
export PDF_END=0.5
16+
export PDF_SAMPLES=1000
17+
export KERNEL_EXP=-3
18+
export DIST="gauss"
19+
20+
source ./experiment.sh
21+
22+
experiment

‎experiments/skewness-evolution/experiments_laplacian.sh

+21
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,21 @@
1+
#!/usr/bin/env bash
2+
# LAPLACIAN !!!
3+
4+
set -euo pipefail
5+
6+
export NITER=1000
7+
export N=3
8+
export M=2
9+
export SIGMAV2=0.00000000001
10+
export EXP=6
11+
export LANG=en_US
12+
export BETA=0.005
13+
export K=256
14+
export PDF_START=0.0
15+
export PDF_END=0.20
16+
export PDF_SAMPLES=10000
17+
export KERNEL_EXP=-3
18+
export DIST="lap"
19+
20+
source ./experiment.sh
21+
experiment

0 commit comments

Comments
 (0)