График магнитного поля орбиты спутника в матлабе - PullRequest
Новая
       30

График магнитного поля орбиты спутника в матлабе

0 голосов
/ 15 мая 2018

Я сделал графический интерфейс в MATLAB, и в один момент я получил трехмерный график Земли, где я показываю орбиту спутника. Что я хочу, это показать в каждой точке (например, я увеличиваю координаты [широта,долготы] на 1 градус и сохраняют ту же высоту) этой орбиты, что величина магнитного поля Земли с векторами (вектор на север, вектор на восток, вектор вниз). Это часть кода, но она не работает 

magnetic_field_ = igrf(time,T(i,1),T(i,2),T(i,3),'geodetic');
[hh(i,1) hh(i,2) hh(i,3)]=igrf(time,T(i,1),T(i,2),T(i,3),'geodetic');

.,. 

% earth
axesm('globe','geoid',wgs84Ellipsoid)
load topo
meshm(topo, topolegend, size(topo)); 
demcmap(topo);
hold on;

.. 

 vv=lla2ecef([T(i,1),T(i,2),T(i,3)]);
 spheroid = referenceEllipsoid('GRS 80');

[Bh(1) Bh(2) Bh(3)]=ned2ecefv(hh(i,1),hh(i,2),hh(i,3),T(i,1),T(i,2),T(i,3));%,spheroid,'degrees');
h=quiver3(vv(1),vv(2),vv(3),Bh(1),Bh(2),Bh(3),0.1,'r');
h.MaxHeadSize=3;
hold on
...