Matlab读入观测数据和导航电文文件的代码

以下是读取观测数据（以RINEX格式为例）和导航电文文件（以SP3格式为例）的Matlab代码示例：

读取RINEX观测数据文件：

% 打开文件
obsfile = 'obs.19o';
fid = fopen(obsfile);

% 读取文件头
line = fgetl(fid);
while ~contains(line, 'END OF HEADER')
    % 处理文件头信息
    line = fgetl(fid);
end

% 读取数据
data = [];
while ~feof(fid)
    line = fgetl(fid);
    if ~isempty(line)
        % 处理数据行
        obs = strsplit(line);
        % 将数据转换为数值型
        obs = cellfun(@str2double, obs, 'UniformOutput', false);
        data = [data; cell2mat(obs)];
    end
end

% 关闭文件
fclose(fid);

读取SP3导航电文文件：

% 打开文件
navfile = 'igs20694.sp3';
fid = fopen(navfile);

% 读取文件头
line = fgetl(fid);
while ~contains(line, 'END OF HEADER')
    % 处理文件头信息
    line = fgetl(fid);
end

% 读取数据
data = [];
while ~feof(fid)
    line = fgetl(fid);
    if ~isempty(line)
        % 处理数据行
        if startsWith(line, '*')
            % 处理卫星编号
            prn = strsplit(line);
            prn = str2double(prn{2});
        else
            % 处理导航电文数据
            nav = strsplit(line);
            % 将数据转换为数值型
            nav = cellfun(@str2double, nav, 'UniformOutput', false);
            nav = cell2mat(nav);
            nav(1) = nav(1) + 2000; % 将年份转换为4位数字
            % 将数据保存到数组中
            data(prn).time(end+1,:) = nav(1:6);
            data(prn).pos(end+1,:) = nav(7:9);
            data(prn).vel(end+1,:) = nav(10:12);
            data(prn).acc(end+1,:) = nav(13:15);
        end
    end
end

% 关闭文件
fclose(fid);

上述代码中，假设RINEX观测数据文件名为`obs.19o`，SP3导航电文文件名为`igs20694.sp3`。读取RINEX文件时，首先读取文件头信息，然后逐行读取数据，将数据转换为数值型后保存到`data`数组中。读取SP3文件时，也先读取文件头信息，然后逐行读取数据，根据行首的卫星编号来区分不同的卫星，将导航电文数据保存到`data`数组中相应的卫星结构体中。