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

news2024/11/30 0:50:47

1.算法运行效果图预览

2.算法运行软件版本

3.部分核心程序

4.算法理论概述

5.算法完整程序工程

1.算法运行效果图预览

2.算法运行软件版本

MATLAB2022a

3.部分核心程序

....................................................................................
% 将所有PRN对应的ECEF位置信息合并成一个表格  
Info2_tab = cat(1,Info2{:});
% 将ECEF位置转换为地理坐标系下的位置  
Geo_tab   = func_calpos(Info2_tab);
% 按照时间对地理坐标位置进行排序  
Geo_tab   = sortrows(Geo_tab);
% 获取排序后的唯一时间点  
Times_tab = unique(Geo_tab.time);
% 读取卫星图像文件'sat.png'，获取图像数据和透明度信息  
[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');

%创建63个卫星图标，初始状态下可能不可见或不在地图上显示（具体取决于func_addsat函数的实现） 
for ii = 1:63
    sat_pos(ii) = func_addsat(ax,0,0,cdata,alpha);
end

% 定义一个匿名函数funPixelPos，用于将地理坐标转换为像素坐标（具体转换方式取决于func_Pixel2Pos函数的实现）
funPixelPos = @(lat,lon) func_Pixel2Pos(ax.UserData,lat,lon,cdata);
% 对于每个唯一的时间点，更新卫星图标的位置和可见性，并更新文本标签的时间显示，然后保存地图图像到变量AA中（不保存到文件）  

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
    end
    pause(0.01);
    title([num2str(ii)]);
end
91