mirrors/gnss-sdr
1
0
Fork 0
mirror of https://github.com/gnss-sdr/gnss-sdr synced 2024-11-16 06:44:57 +00:00
Code Issues Releases Wiki Activity

ADD: matlab skyplot implemented

Browse Source
This commit is contained in:
miguekf 2023-01-03 13:40:58 +01:00
parent aa07ba17fd
commit 1b0a3f6b37
1 changed files with 30 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
src/utils/matlab/vtl/pvt_raw_plotting_SPIRENT_GnssSDR.m
View File

@ -30,7 +30,7 @@ channels=6;
TTFF_sec=41.48; TTFF_sec=41.48;
%% ============================ PARSER TO STRUCT ============================ %% ============================ PARSER TO STRUCT ============================
plot_skyplot=0; % not yet implemented plot_skyplot=1;
plot_reference=1; plot_reference=1;
load_observables=1; load_observables=1;
advance_vtl_data_available=1; advance_vtl_data_available=1;
@ -371,4 +371,32 @@ if(load_observables)
grid on grid on
end end
%% %%
dopp_filtered_plotting dopp_filtered_plotting
%%
labels = ["G1", "G22", "E17", "E11"]; %inizialitation neccesary
if(plot_skyplot)
for i=1:length(refSatData.GPS.SIM_time)
for j=1:length(refSatData.GPS.series(i).Sat_PRN)
labels{i,j} = ['G' num2str(refSatData.GPS.series(i).Sat_PRN(j))];
if(refSatData.GPS.series(i).azimuth(j)<0)
az(i,j) = rad2deg(refSatData.GPS.series(i).azimuth(j))+360;
else
az(i,j) = rad2deg(refSatData.GPS.series(i).azimuth(j));
end
el(i,j) = rad2deg(refSatData.GPS.series(i).elevation(j));
end
groups(i,:) = categorical([0 0 1 1], [0 1], ["GPS", "Galileo"]);
end
% skyplot(az(i,:), el(i,:), labels, "GroupData", groups)
figure;
%skyplot(az(500,:), el(500,:), labels(500,:))
sp = skyplot(az(1,:),el(1,:),labels(1,:));
for idx = 1:length(az)
set(sp,AzimuthData=az(1:idx,:),ElevationData=el(1:idx,:));
drawnow limitrate
end
legend('GPS')
end