From addd52f14851cfe2dcccb49a7a13f9b13a192620 Mon Sep 17 00:00:00 2001 From: ha7ilm Date: Sun, 30 Apr 2017 10:06:31 +0200 Subject: [PATCH] Now showing E_b/N_0 in the plot --- grc_tests/bpsk31_tedvar.m | 16 ++++++++++------ 1 file changed, 10 insertions(+), 6 deletions(-) diff --git a/grc_tests/bpsk31_tedvar.m b/grc_tests/bpsk31_tedvar.m index 9416ed8..32729de 100755 --- a/grc_tests/bpsk31_tedvar.m +++ b/grc_tests/bpsk31_tedvar.m @@ -54,13 +54,13 @@ function fmtplot(h) set(FN,'FontName','times'); FS = findall(h,'-property','FontSize'); set(FS,'FontSize',18); - xlabel('SNR [dB]'); + xlabel('E_b/N_0 [dB]'); ylabel('Phase error variance [rad^2]'); end -snrs_gardner=0:5:30 -%snrs_gardner=[10] -error_values_gardner=mkvarplot('EARLYLATE',snrs_gardner); +snrs=-5:5:30 +%snrs=[10] +error_values_gardner=mkvarplot('EARLYLATE',snrs); %{ snrs_earlylate=0:256 error_values_earlylate=mkvarplot('EARLYLATE',snrs_earlylate); @@ -69,8 +69,12 @@ error_values_earlylate=mkvarplot('EARLYLATE',snrs_earlylate); %graphics_toolkit("gnuplot") h=figure(1); -semilogy(snrs_gardner, error_values_gardner, 'linewidth', 2); -title('Variance curve'); + +ebn0=snrs-13.26-10*log10(1/256.) +%13.56 dB is the difference between the real (measured) SNR and the number input to awgn_cc + +semilogy(ebn0, error_values_gardner, 'linewidth', 2); +title('Estimation variance'); fmtplot(h) pause