цос. Туралиев лаба2. Преобразование координат
 Скачать 2.13 Mb.
|
|
Министерство образования и науки Кыргызской Республики Кыргызский Государственный Технический Университет им. И.Раззакова Кыргызско-Германский Технический Институт Кафедра: ” Телематика” О Т Ч Е Т По дисциплине: Глобальные сети На тему: Преобразование координат Выполнил: Туралиев Ж Проверил: Мамбеталиев З.С. БИШКЕК 2022 Цель лабораторной работы: Изучение и практическое освоение систем координат, применяемых в спутниковых радионавигационных системах Задание 1. По географической карте определите широту (градусы/минуты/секунды) и долготу (градусы/минуты/секунды) любого города Европы. Преобразуйте выбранные значения в радианы.      Задание 2. Используя результаты, полученные в п. 12, в качестве входных данных в соответствии с п. п. 6, 7 сформируйте m-файл и выполните.  Задание 3. Задайте координаты двух объектов, находящихся в прямой видимости широта (градусы, минуты, секунды), долгота (градусы, минуты, секунды), высота (метры). Преобразуйте заданные координаты широты и долготы в градусы.   Задание 4. Задайте координаты и скорости объекта satpos_eci, истинное звездное время s0, текущее время ti. Используя эти исходные данные и п. п. 16, 17 сформируйте m-файл и выполните его.   Задание 5. Выберите эллипсоид и задайте его полуоси a, b; задайте текущее время ti, истинное звездное время time_s0, координаты потребителя llh_loc. Используя эти исходные данные и п. п. 19, 20 сформируйте m-файл и выполните его.  function [XYZ] = ECEFLLH_N(llh,ab) %Имя функции: ECEFLLH_N %Назначение функции: преобразование координат из географической системы в прямоугольную %Входные данные: %llh.lon-долгота; %llh.lat-широта; %llh.h-высота; %ab.a-большая полуось эллипсоида; %ab.b- малая полуось эллипсоида в WGS-84; %Выходные данные: %XYZ.x,XYZ.y,XYZ.z- координаты X, Y, Z соответственно в ECEF % Справочные данные: %ECEF- прямоугольная геоцентрическая система координат %a=6378137.0 (м)- большая полуось эллипсоида для WGS-84; %b=6356752.314 (м)- малая полуось эллипсоида для WGS-84; %A_PZ90_M =6 378 136 (м)- большая полуось эллипсоида для ПЗ 90; %B_PZ90_M = 6356751.36174 (м)- малая полуось эллипсоида для ПЗ 90; a2=ab.a*ab.a; b2=ab.b*ab.b; r=a2/sqrt(a2*cos(llh.lat)*cos(llh.lat)+b2*sin(llh.lat)*sin(llh.lat)); XYZ.x=(r+llh.h)*cos(llh.lat)*cos(llh.lon); XYZ.y=(r+llh.h)*cos(llh.lat)*sin(llh.lon); XYZ.z=(b2/a2*r+llh.h)*sin(llh.lat); function [XYZ] = ECEFLLH_N(llh,ab) %Имя функции: ECEFLLH_N %Назначение функции: преобразование координат из географической системы в прямоугольную %Входные данные: %llh.lon-долгота; %llh.lat-широта; %llh.h-высота; %ab.a-большая полуось эллипсоида; %ab.b- малая полуось эллипсоида в WGS-84; %Выходные данные: %XYZ.x,XYZ.y,XYZ.z- координаты X, Y, Z соответственно в ECEF % Справочные данные: %ECEF- прямоугольная геоцентрическая система координат %a=6378137.0 (м)- большая полуось эллипсоида для WGS-84; %b=6356752.314 (м)- малая полуось эллипсоида для WGS-84; %A_PZ90_M =6 378 136 (м)- большая полуось эллипсоида для ПЗ 90; %B_PZ90_M = 6356751.36174 (м)- малая полуось эллипсоида для ПЗ 90; a2=ab.a*ab.a; b2=ab.b*ab.b; r=a2/sqrt(a2*cos(llh.lat)*cos(llh.lat)+b2*sin(llh.lat)*sin(llh.lat)); XYZ.x=(r+llh.h)*cos(llh.lat)*cos(llh.lon); XYZ.y=(r+llh.h)*cos(llh.lat)*sin(llh.lon); XYZ.z=(b2/a2*r+llh.h)*sin(llh.lat); |