mirrors/gnss-sdr
1
0
Fork 0
mirror of https://github.com/gnss-sdr/gnss-sdr synced 2024-12-16 13:10:35 +00:00
Code Issues Releases Wiki Activity

Minor changes in matlab acquisition grid plot script

Browse Source
This commit is contained in:
Antonio Ramos 2018-01-24 11:17:33 +01:00
parent 2153e24169
commit 2e2c3474be
1 changed files with 14 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
src/utils/matlab/plot_acq_grid.m
View File

@ -31,7 +31,7 @@
% */ % */
%%%%%%%%% ¡¡¡ CONFIGURE !!! %%%%%%%%%%%%% %%%%%%%%% ¡¡¡ CONFIGURE !!! %%%%%%%%%%%%%
sat = 27; sat = 30;
n_chips = 1023; n_chips = 1023;
system = 'G'; % GPS = 'G', Galileo = 'E' system = 'G'; % GPS = 'G', Galileo = 'E'
%%% True for light grid representation %%% True for light grid representation
@ -40,23 +40,33 @@ lite_view = true;
n_samples_per_chip = 4; n_samples_per_chip = 4;
path='/home/aramos/signals/'; path='/home/aramos/signals/GNSS-IN-THE-SPACE/CAPTURES SPIRENT/acq/';
file=['acq_' system '_sat_' num2str(sat) '.mat']; file=['acq_' system '_sat_' num2str(sat) '.mat'];
load([path file]); load([path file]);
[n_fft n_dop_bins] = size(grid); [n_fft n_dop_bins] = size(grid);
[d_max f_max] = find(grid == max(max(grid)));
freq = (0 : n_dop_bins - 1) * doppler_step - doppler_max; freq = (0 : n_dop_bins - 1) * doppler_step - doppler_max;
delay = (0 : n_fft - 1) / n_fft * n_chips; delay = (0 : n_fft - 1) / n_fft * n_chips;
figure(1) figure(1)
if(lite_view == false) if(lite_view == false)
mesh(freq, delay, grid) surf(freq, delay, grid)
else else
delay_interp = (0 : n_samples_per_chip * n_chips - 1) / n_samples_per_chip; delay_interp = (0 : n_samples_per_chip * n_chips - 1) / n_samples_per_chip;
grid_interp = spline(delay, grid', delay_interp)'; grid_interp = spline(delay, grid', delay_interp)';
mesh(freq, delay_interp, grid_interp) surf(freq, delay_interp, grid_interp)
end end
xlabel('Doppler shift / Hz') xlabel('Doppler shift / Hz')
ylabel('Code delay / chips') ylabel('Code delay / chips')
zlabel('Test statistics') zlabel('Test statistics')
figure(2)
subplot(2,1,1)
plot(freq, grid(d_max, :))
xlabel('Doppler shift / Hz')
ylabel('Test statistics (fixed delay)')
subplot(2,1,2)
plot(delay, grid(:, f_max))
xlabel('Code delay / chips')
ylabel('Test statistics (fixed Doppler shift)')