file='D:\matlab\par&res\sh0310_G.sp3';
% gc_eqwgt=load('D:\matlab\par&res\par_2022031_gc_eqwgt');
file1='D:\matlab\par&res\tian21951.grc.sp3';
fid=fopen(file,'rt');
fid1=fopen(file1,'rt');
% year=[];
flag1=1;
flag2=1;
flag3=1;
flag4=1;
while ~feof(fid)
line=fgetl(fid);
if (length(line)<33)
year_sh(flag1,1)=str2num(line(4:7));
month_sh(flag1,1)=str2num(line(10:11));
day_sh(flag1,1)=str2num(line(12:13));
hour_sh(flag1,1)=str2num(line(15:16));
min_sh(flag1,1)=str2num(line(18:19));
s_sh(flag1,1)=str2num(line(21:31));
flag1=flag1+1;
elseif findstr(line,'PL01')
x_sh(flag2,1)=str2num(line(7:18));
y_sh(flag2,1)=str2num(line(21:32));
z_sh(flag2,1)=str2num(line(35:46));
flag2=flag2+1;
end
end
sh_out=[year_sh,month_sh,day_sh,hour_sh,min_sh,s_sh,x_sh,y_sh,z_sh] ;
% while ~feof(fid1)
% line2=fgetl(fid1);
% if (length(line2)<33)
% year_gc(flag3,1)=str2num(line2(4:7));
% month_gc(flag3,1)=str2num(line2(10:11));
% day_gc(flag3,1)=str2num(line2(12:13));
% hour_gc(flag3,1)=str2num(line2(15:16));
% min_gc(flag3,1)=str2num(line2(18:19));
% s_gc(flag3,1)=str2num(line2(21:31));
% flag3=flag3+1;
% elseif findstr(line2,'PT01')
% x_gc(flag4,1)=str2num(line2(7:18));
% y_gc(flag4,1)=str2num(line2(21:32));
% z_gc(flag4,1)=str2num(line2(35:46));
% flag4=flag4+1;
% end
%
% end
% gc_out=[year_gc,month_gc,day_gc,hour_gc,min_gc,s_gc,x_gc,y_gc,z_gc];
%做差
% time_gc=(gc_eqwgt(:,1)-59610)*86400;
% time_sh=hour*3600+min*60+s;
% if ((abs(time_sh-time_gc))<0.1)
% dx=x-gc_eqwgt(:,3);
% dy=y-gc_eqwgt(:,4);
% dz=z_gc_eqwgt(:,5);
% rms_x=rms(dx);
% rms_y=rms(dy);
% rms_z=rms(dz);
% end
%
Matlab读取SP3文件_matlab
版权申诉
5星 · 超过95%的资源 119 浏览量
2022-09-20
13:22:40
上传
评论 3
收藏 737B RAR 举报
weixin_42653672
- 粉丝: 93
- 资源: 1万+
评论13