GNSS精密星历读取（matlab）

该代码用于从SP3文件中提取第一颗卫星的X、Y、Z坐标以及2到9颗卫星的X坐标。代码首先读取文件头信息，然后循环读取数据，将数据存储到数组中。目前代码仅提取了部分卫星数据，若需其他数据，需进行相应改进。

该代码仅获取了精密星历SP3文件中第一颗卫星在不同时间点的X、Y、Z坐标，以及卫星2至9号的X坐标位置。若需要其他数据，建议对代码进行优化。

clc
clear
使用fopen函数打开SP3文件...，其中文件路径为'ESA0OPSRAP_20222770000_01D_15M_ORB.SP3'。
% 读取文件头
读取文件头的第一行，赋值给变量line。
% 判断文件头是否为'NONE'
使用strcmp函数判断当前行是否为'NONE'，若不为'NONE'则继续读取下一行。
% 读取文件内容
初始化变量j为1，变量str为'EOF'，然后使用fgetl函数读取每一行，直到当前行等于'EOF'为止。

if the line from positions 1 to 3 is a string,
continue to the next iteration,
%Check if the line is a string type; if not, exit the loop
ephx(1,j) is assigned the value of X(1,1),
and similarly for the other variables.
for k ranging from 1 to 51,
assign X(1,k) to the double value of the substring from position 6 to 18,
Y(1,k) to the substring from 20 to 32,
Z(1,k) from 34 to 46,
and G(1,k) from 50 to 60.
then,
read the next line of data,
%First, test with 10 satellites
assign the first 10 values of X, Y, and Z to the ephx and ephy arrays.

j=j+1;
end
fclose(fid);