%EASY2 Convert observation time into sow.
% We read the corresponding RINEX navigation file
% and reformat the data into the Matlab matrix Eph.
% For given SV we find the corresponding column in Eph
% and call the basic satpos function
% Compute sow for first epoch in observation file xbgd173A.16o
% 2016 6 21 6 44 16.0000000
clear all;
jd = julday(2016,6,21,6+44/60);
[week,sow] = gps_time(jd);
% Read RINEX ephemerides星历表 file and convert to
% internal Matlab format
Eph =rinexe('brdm1730.16p','eph.dat');%星历表,输出文件
% We identify the observed satellites in line 29 of RINEX file brdm1730.16p
svs = [1 2 3 4 5 6 7 8 9 10 11 12 14 15]; %卫星号
for t = 1:length(svs)
col_Eph(t) = find_eph(Eph,svs(t),sow);
sat(1:3,t) = satpos(sow,Eph(:,col_Eph(t))); %由星历计算ECEF卫星位置
end
sat % position of svs in ECEF system
%%%%%%%%%%%%%%%%%%%%% end easy2.m %%%%%%%%%%%%%%%
matlab学习.rar_GPS_GPS 定位 MATLAB_gps 原理 matlab_gps定位matlab
版权申诉
63 浏览量
2022-09-14
23:42:05
上传
评论
收藏 1KB RAR 举报
weixin_42651887
- 粉丝: 79
- 资源: 1万+