以darw站2014年数据为例
clear all;
clc;day='缺失天数'; %定义变量名day,记录缺失的文件年积日
m=[]; %定义空矩阵,存储数据
for n=1:366 %m文件天数1天到366天
ns=num2str(n,'%03d'); % 数值转字符串
name=['darw',ns,'0.14m']; % 字符串拼接得到文件名称
ffid=fopen(name,'rt'); % 打开文本文件,从第一行开始读取,数据存储在cell里
if ffid>0
i = 1;
while feof(ffid) == 0
tline{i,1} = fgetl(ffid);
i = i+1;
end
% 提取间隔2h的数据,从0时至22时,RINEX m文件采样间隔为30s,2h间隔为240行
for i=18:240:2658
a((i+222)/240,:)=str2num(tline{i,1});
end
% 提取压强p 和温度 t的值
for i=1:12
b(i,1)=a(i,7);
b(i,2)=a(i,8);
end
m=cat(1,m,b);
else
day=strcat(day,'-',ns);
end
end
xlswrite('14m.xlsx',m)