根据星历文件实现卫星的动态运行模拟matlab仿真

发布时间：2025-05-21 19:34

Info2_tab = cat(1,Info2{:});

Geo_tab = func_calpos(Info2_tab);

Geo_tab = sortrows(Geo_tab);

Times_tab = unique(Geo_tab.time);

[cdata,~,alpha] = imread('dat\sat.png');

cdata=cdata/1.5;

figure;

ax = worldmap('World');

setm(ax, 'Origin', [0 0 0]);

ax.Position =[0 0 1 1];

land = shaperead('landareas', 'UseGeoCoords', true);

geoshow(ax, land, 'FaceColor', [0.5 0.7 0.5]);

lakes = shaperead('worldlakes', 'UseGeoCoords', true);

geoshow(lakes, 'FaceColor', 'blue');

rivers = shaperead('worldrivers', 'UseGeoCoords', true);

geoshow(rivers, 'Color', 'blue');

for ii = 1:63

sat_pos(ii) = func_addsat(ax,0,0,cdata,alpha);

end

funPixelPos = @(lat,lon) func_Pixel2Pos(ax.UserData,lat,lon,cdata);

for ii = 1:numel(Times_tab)

t = Times_tab(ii);

data = rowfun(funPixelPos,Geo_tab(t,1:2),'NumOutputs',2,'OutputFormat','cell');

for jj = 11:63

if jj<=size(data,1)

sat_pos(jj).Visible = 'on';

sat_pos(jj).XData = data{jj,1};

sat_pos(jj).YData = data{jj,2};

else

sat_pos(jj).Visible = 'off';

end

pause(0.01);

title([num2str(ii)]);

end