matlab绘制星下点轨迹,MATLAB绘制GPS星下点轨迹图

MATLAB绘制GPS星下点轨迹图

轨迹计算部分参考链接:

https://wenku.baidu.com/view/45bd098d4a7302768e9939cf.html

本文对上述matlab代码进行了整理与修改:

增加了底图并进行了相关的图形美化。

将轨道六参数设置为GPS相关参数。

从原有的1颗卫星轨迹绘制增加至24颗。

可进一步修改的部分:

通过修改底图绘制部分相关参数可以进行图片的基础修改。

通过修改轨道六根数的大小以及循环次数可以实现其他卫星导航系统的星下点轨迹绘制。

代码如下:

clc

close all

clear all

%% % % % % % % % % % % % % % % % % 底图绘制 % % % % % % % % % % % % % % % % % % % % % %

h = geoshow('landareas.shp', 'FaceColor', 'c');

grid on

hold on

xlabel('Longitude');

ylabel('Latitude');%坐标轴标题

set(gca,'Ylim',[-90,90],'ytick',[-90:30:90]);

set(gca,'yticklabel',{'90°S','60°','30°','0°','30°','60°','90°N'});

set(gca,'Xlim',[-180,1

  • 2
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
您好,要绘制星下点轨迹,需要知道卫星的轨道参数和地球的形状参数。一般来说,可以使用开源的轨道计算工具如Orekit或者SPG4来计算卫星的位置信息,然后再将其投影到地球表面上得到星下点的位置。接下来使用Matlab中的绘函数即可绘制星下点轨迹。 以下是一个简单的Matlab代码示例,用于绘制一个卫星在地球表面上的星下点轨迹: ```matlab % 定义地球形状参数 R = 6378.137; % 地球半径,单位:km f = 1/298.257223563; % 地球扁率 % 定义卫星轨道参数 a = 7000; % 半长轴,单位:km e = 0.001; % 偏心率 i = 45; % 轨道倾角,单位:度 Omega = 0; % 升交点赤经,单位:度 omega = 0; % 近地点幅角,单位:度 M0 = 0; % 平近点角,单位:度 n = sqrt(398600.4418/(a^3)); % 平均运动角速度,单位:rad/s % 计算卫星位置信息 t = linspace(0, 2*pi/n, 1000); % 时间序列 E = kepler(M0, e, t); % 求解偏近点角 r = a*(1 - e*cos(E)); % 距离 v = sqrt(398600.4418*(2/r - 1/a)); % 速度 theta = atan2(sqrt(1 - e^2)*sin(E), cos(E) - e); % 真近点角 u = omega + theta; % 真近点经度 i = i*pi/180; Omega = Omega*pi/180; u = u*pi/180; X = r.*(cos(Omega).*cos(u) - sin(Omega).*sin(u).*cos(i)); Y = r.*(sin(Omega).*cos(u) + cos(Omega).*sin(u).*cos(i)); Z = r.*sin(u).*sin(i); % 投影到地球表面上得到星下点位置 lat = asin(Z./sqrt(X.^2 + Y.^2 + Z.^2))*180/pi; % 纬度 lon = atan2(Y, X)*180/pi; % 经度 % 绘制星下点轨迹 figure; worldmap('World'); load coastlines; plotm(coastlat, coastlon); plotm(lat, lon, 'r'); title('Satellite Ground Track'); ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值