matlab编程求卫星轨道长度,GPS卫星轨道计算及其MATLAB仿真.pdf

GPS 卫星轨道计算及其MATLAB 仿真

黎奇，白征东，李帅，陈波波

(清华大学地球空间信息研究所，北京100084)

一、程序设计思路

1. 读取RINEX 文件(注意：文件路径)

2. 计算测量日周积秒(测量日的格里历→GPST )

3. 按卫星轨道计算步骤计算WGS-84 坐标系坐标 (内插)

4. 按需要将WGS-84 坐标系下坐标转换为所需坐标系坐标

5. 画图输出

主程序

function tracePlot() % 主程序

format long

%读取RINEX格式n文件

data = RinexNreader('C:\Users\Richie\Desktop\卫星轨迹\brdc2750.16n',1); % 注意读取路径和卫星编号

%计算测量日的周积秒

[JD,FOD,GPSW,SOW,DOY,DOW] =

GCtoGPS(data(1,1),data(2,1),data(3,1),data(4,1),data(5,1),data(6,1));

t0 = SOW; % SOW-周积秒

%计算卫星坐标

i = 1;

for t = t0:10*60:(t0+86400) % 从 t0 到 t0+86400秒(1天) 间隔 600秒(10min)，决定轨迹疏密

satPosECEF(:,i) = orbitDetermine(data,t);

% 输入：卫星PRN编号、时刻t

% 输出：卫星WGS-84坐标x、y、z

[B_ECEF(i),L_ECEF(i)] =