rpitx/csdr
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Very detailed, but correct looking phase variance diagram

Browse source
This commit is contained in:
ha7ilm 2017-04-29 23:56:29 +02:00
parent ec38e3be7d
commit 3f527311d6
1 changed files with 4 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
grc_tests/bpsk31_tedvar.m
View file

@ -5,7 +5,7 @@
%pkg install -forge parallel %pkg install -forge parallel
pkg load parallel pkg load parallel
system('cat /dev/urandom | csdr pack_bits_8to1_u8_u8 | csdr psk_modulator_u8_c 2 | csdr gain_ff 0.25 | csdr psk31_interpolate_sine_cc 256 | csdr add_n_zero_samples_at_beginning_f 170 | dd iflag=fullblock bs=4M count=1 of=/tmp/psk31-raw-data'); system('cat /dev/urandom | csdr pack_bits_8to1_u8_u8 | csdr psk_modulator_u8_c 2 | csdr gain_ff 0.25 | csdr psk31_interpolate_sine_cc 256 | csdr add_n_zero_samples_at_beginning_f 170 | dd iflag=fullblock bs=8M count=1 of=/tmp/psk31-raw-data');
function output=shrun(cmd, type, minsize) function output=shrun(cmd, type, minsize)
SIGTERM=15; SIGTERM=15;
@ -30,10 +30,10 @@ end
function variance=run_var(snr, which_ted) function variance=run_var(snr, which_ted)
disp('ran a command') disp('ran a command')
out_vect=shrun(sprintf('cat /tmp/psk31-raw-data | csdr awgn_cc %d --snrshow-no| csdr timing_recovery_cc %s 256 --add_q --output_indexes | CSDR_FIXED_BUFSIZE=131072 csdr normalized_timing_variance_u32_f 256 85', snr, which_ted), 'float32', 1); out_vect=shrun(sprintf('cat /tmp/psk31-raw-data | csdr awgn_cc %d --snrshow-no| csdr timing_recovery_cc %s 256 --add_q --output_indexes | CSDR_FIXED_BUFSIZE=2048 csdr normalized_timing_variance_u32_f 256 85', snr, which_ted), 'float32', 2);
disp('run_var output:'); disp('run_var output:');
out_vect' out_vect'
variance=out_vect(1); variance=out_vect(2);
end end
function variances=mkvarplot(which_ted, snrs) function variances=mkvarplot(which_ted, snrs)
@ -58,7 +58,7 @@ function fmtplot(h)
ylabel('Phase error variance [rad^2]'); ylabel('Phase error variance [rad^2]');
end end
snrs_gardner=-70:5:40 snrs_gardner=-70:2:100
%snrs_gardner=[10] %snrs_gardner=[10]
error_values_gardner=mkvarplot('GARDNER',snrs_gardner); error_values_gardner=mkvarplot('GARDNER',snrs_gardner);
%{ %{