gnssconstellation

Спутниковые местоположения в требуемое время

Описание

пример

[satPos,satVel] = gnssconstellation(t) возвращает спутниковые положения и скорости в datetime t. Функция возвращает положения и скорости в системе координат Сосредоточенного землей зафиксированного землей (ECEF) в метрах и метрах в секунду, соответственно.

Примеры

свернуть все

Получите текущие спутниковые положения и скорости от спутников GNSS. Доступ к орбитальным параметрам из Спецификации Интерфейса IS-GPS-200K и вычисляет положение и скорости в координатах ECEF на данное время. Отобразите спутниковые положения.

t = datetime('now','TimeZone','Local');
[satPos,satVel] = gnssconstellation(t);
disp(satPos)
   1.0e+07 *

   -0.7187   -1.5277   -2.0503
   -2.2398    0.6413    1.2752
    2.6423    0.1437   -0.2272
    2.3368    0.8957    0.8896
   -1.1300    1.1290    2.1220
   -0.3523    2.0639    1.6341
   -1.5386   -1.9706    0.8965
    1.6941    0.3642   -2.0128
   -0.8861   -2.4909   -0.2541
   -1.1043   -1.0615    2.1698
   -1.4485    2.2150   -0.2227
   -0.4932    2.2398   -1.3396
    1.9885   -0.9085   -1.5082
    1.3749   -1.2014    1.9288
    0.0563    1.5341   -2.1674
    1.1559    1.4650   -1.8899
   -2.4644   -0.7441    0.6537
    2.3657   -0.5020    1.0981
    1.0717    2.4300    0.0237
   -1.6227   -0.0421   -2.1022
    0.3436   -2.0711    1.6268
   -1.7308   -1.3325   -1.5109
   -1.9963    1.0069   -1.4335
    1.1491   -2.3064   -0.6439
    0.3000   -2.1699   -1.5019
    1.2905    0.8149    2.1736
    0.3488    1.7945    1.9267

Используйте lookangles функция, чтобы получить азимут и углы возвышения спутников для данных положений спутника и приемника. Задайте угол маски 5 градусов. Получите спутниковые положения с помощью gnssconstellation функция.

Задайте положение приемника в геодезических координатах (широта, долгота, высота).

recPos = [42 -71 50];

Получите спутниковые положения в течение текущего времени.

t = datetime('now');
gpsSatPos = gnssconstellation(t);

Задайте угол маски 5 градусов.

maskAngle = 5;

Получите азимут и углы взгляда вертикального изменения для спутниковых положений. vis выведите указывает, какие спутники отображаются. Получите общее количество с помощью nnz.

[az,el,vis] = lookangles(recPos,gpsSatPos,maskAngle);
fprintf('%d satellites visible at %s.\n',nnz(vis),t);
8 satellites visible at 25-Aug-2021 19:44:43.

Входные параметры

свернуть все

Текущее время для спутниковой симуляции в виде скалярного datetime массив.

Часовой пояс по умолчанию для datetime массивом является UTC. Для получения информации об определении различного часового пояса смотрите datetime.

Время начала GPS является 6 января 1980 полуночью, UTC. Определение любого datetime до этого времени будет использовать время начала GPS.

Пример: datetime('now','TimeZone','Local');

Типы данных: datetime

Выходные аргументы

свернуть все

Спутниковые положения в системе координат Сосредоточенного землей зафиксированного землей (ECEF) в метрах, возвращенных как N-by-3 матрица скаляров. N является количеством спутников в созвездии.

Типы данных: single | double

Спутниковые скорости в системе координат Сосредоточенного землей зафиксированного землей (ECEF) в метрах в секунду, возвращенный как N-by-3 матрица скаляров. N является количеством спутников в созвездии.

Типы данных: single | double

Больше о

свернуть все

Орбитальные параметры

Начальные спутниковые положения и скорости заданы орбитальными параметрами в Таблице A.2-2 в GPS Стандарт Эффективности SPS и даны в координатах Сосредоточенного землей зафиксированного землей (ECEF).

Вычисления положения используют уравнения из Таблицы, 30-II в той же Спецификации Интерфейса IS-GPS-200K.

Скоростные расчеты используют уравнения 8.21-8.27 в [1].

Ссылки

[1] Рощи, Пол Д. Принципы GNSS, инерционного, и мультидатчик интегрированные системы навигации. Бостон: дом Artech, 2013.

Расширенные возможности

Введенный в R2021a