Создайте модель датчика GNSS как gnssSensor (Navigation Toolbox) система Object™.

gnss = gnssSensor; 

Задайте положение и скорость датчика. Симулируйте показания датчика и получите состояние от видимых спутников. Сохраните азимут и углы возвышения как векторы.

pos = [0 0 0]; 
vel = [0 0 0]; 
[~, ~, status] = gnss(pos, vel); 
satAz = status.SatelliteAzimuth; 
satEl = status.SatelliteElevation; 

Постройте спутниковые постионы.

skyplot(satAz,satEl) 

Анимируйте траекторию спутниковых положений в зависимости от времени от датчика GNSS.

Инициализируйте фигуру графика неба. Укажите соответствующую продвигающуюся во время информацию.

skyplotHandle = skyplot(0,0); 

numHours = 12; 
dt = 100; 
numSeconds = numHours * 60 * 60; 
numSimSteps = numSeconds/dt;

Создайте модель датчика GNSS как gnssSensor (Navigation Toolbox) система Object™.

gnss = gnssSensor('SampleRate', 1/dt); 

Выполните итерации через временные шаги и сделайте следующее:

for i = 1:numSimSteps 

    [~, ~, status] = gnss([0 0 0],[0 0 0]); 

    satAz = status.SatelliteAzimuth; 
    satEl = status.SatelliteElevation; 

    set(skyplotHandle,'AzimuthData',satAz,'ElevationData',satEl);
    
    drawnow 
end 

Загрузите азимут и данные о вертикальном изменении из файла журнала, сгенерированного датчиком спутника Adafruit® GPS. Данные, обеспеченные в этом примере, содержат азимут и вертикальное изменение каждого спутника и псевдослучайного шума (PRN) коды. Сохраните эти значения как векторы.

load('gpsHWInfo','hwInfo') 
satAz = hwInfo.SatelliteAzimuths; 
satEl = hwInfo.SatelliteElevations; 
prn = hwInfo.SatellitePRNs; 

Разделите спутники на основе кодов PRN. Чтобы коррелировать каждое положение с группой, создайте categorical массив. Для этого набора спутников только те с PRNs меньше чем 32 используются в решении для расположения.

isUnused = (prn > 32);
group = categorical(isUnused,[false true],["Used in Positioning Solution" "Unused"]);

Визуализируйте спутники и задайте категориальные группы в GroupData аргумент значения имени. Задайте PRN как метку для каждой точки. Покажите легенду.

skyplot(satAz,satEl,prn,GroupData=group) 
legend('Used','Unused')