本文基于MATLAB，首先计算卫星轨道在地球上的垂直投影x,y，然后设置miller地图投影方式，最后用m_line( )函数将星下点轨迹绘制出来。中间过程演示了不同相位差时，卫星轨道在地球上的垂直投影（星下点轨迹）的变化。工具/原料moreMATLABm_mapm_line方法/步骤1/6分步阅读

第一，启动MATLAB，新建脚本（Ctrl+N），输入如下代码，用于绘制卫星的星下点轨迹。

close all; clear all; clc

x=;

y1=atan(tan(60*pi/180)*cos((x-30)*pi/180))*180/pi;

m_proj(miller,lat,82);

m_coast(color,);

m_line(x,y1,linewi,3,color,r);hold on

m_grid(linestyle,none,box,fancy,ticklength,0.01);

其中x,y1是计算的卫星轨道在地球上的垂直投影（星下点轨迹）。m_proj设置地图投影为miller方式。m_coast绘制海岸线。m_line绘制x,y1星下点轨迹。m_grid设置网格线为none和设置坐标轴边框属性等。

需要注意的是，使用上述代码时，请确保MATLAB已经安装m_map工具箱，具体可以安装方法百度“MATLAB学习与使用：安装M_Map工具箱和测试使用”。

2/6

第二，保存和运行上述脚本，得到如下图形：地图投影为miller的地图投影方式，红线代表卫星运行轨道在地球上的垂直投影（星下点轨迹）。

3/6

第三，增加y2=atan(tan(60*pi/180)*cos((x-120)*pi/180))*180/pi，y2相比y1相位差90度，因为y1是x-30，y2是x-120。同时绘制y1和y2，观测相位差90情况下，星下点轨迹的变化情况。

4/6

第四，保存和运行上述脚本，得到如下图形：红色和绿色两条星下点轨迹，彼此相位差90度。

5/6

第五，增加y3=atan(tan(60*pi/180)*cos((x-210)*pi/180))*180/pi，y3是x-210，与y2相比相位差90，与y1相比相位差180度。同时绘制y1，y2，y3，观测不同相位差下的星下点轨迹。

6/6

第六，保存和运行上述脚本，最终得到如下图形：三条星下点轨迹图形，它们彼此存在相位差。

MATLABM_LINE卫星轨道星下点轨迹