function?readGPSdata
%觀測文件讀取程序
clear?all;
clc;
[filenamepathname]=uigetfile（‘*.??O‘‘打開O文件‘）;
fid=fopen（strcat（pathnamefilename）‘rt‘）;
if?fid==-1
????msgbox（‘文件選取出錯，讀取程序無法運行‘‘warning‘‘warn‘）;
????return;
end
fid2=fopen（‘觀測文件備份.txt‘‘wt‘）;
fprintf（fid2‘觀測文件數據備份\n‘）;
%讀取文件頭
while?（1）
????line=fgets（fid）;
????%讀取一行數據
????if?（line==-1）
????????break;
????end
?????if?（findstr（line‘APPROX?POSITION?XYZ‘）~=0）%測站近似坐標
????????approxx=str2num（line（1:14））;
????????approxy=str2num（line（15:28））;
????????approxz=str2num（line（29:42））;
?????end
????if?（findstr（line‘DELTA?H/E/N‘）~=0）%獲取天線參數:高，東向、北向的偏心
????????anth=str2num（line（1:14））;
????????ante=str2num（line（15:28））;
????????anth=str2num（line（29:42））;
????end
????if?（findstr（line‘TYPES?OF?OBSERV‘）~=0）%觀測值種類
????????signalnum=str2num（line（5:6））;
???????typeNum=1;
????????for?i=1:signalnum%使編號對應為0-L1;1-L2;2-C1;3-P1;4-P2;5-D1;6-D2
????????????ii=11+（i-1）*6;
????????????if?（line（ii:ii+1）==‘L1‘）
????????????????type（typeNum）=0;
????????????????typeNum=typeNum+1;
????????????elseif?（line（ii:ii+1）==‘L2‘）
????????????????type（typeNum）=1;
????????????????typeNum=typeNum+1;
????????????elseif?（line（ii:ii+1）==‘C1‘）
????????????????type（typeNum）=2;
????????????????typeNum=typeNum+1;
????????????elseif?（line（ii:ii+1）==‘C2‘）
????????????????type（typeNum）=3;
????????????????typeNum=typeNum+1;?
????????????????
????????????elseif?（line（ii:ii+1）==‘P1‘）
????????????????type（typeNum）=4;
????????????????typeNum=typeNum+1;
????????????elseif?（line（ii:ii+1）==‘P2‘）
????????????????type（typeNum）=5;
????????????????typeNum=typeNum+1;
????????????elseif?（line（ii:ii+1）==‘D1‘）%多普勒頻率
????????????????type（typeNum）=6;
????????????????typeNum=typeNum+1;
????????????elseif?（line（ii:ii+1）==‘D2‘）
????????????????type（typeNum）=7;
???????????????break;
????????????end
????????end
????end
????if?（findstr（line‘INTERVAL‘）~=0）%觀測歷元的間隔
????????inter=str2num（line（1:11））;
????end
????if?（findstr（line‘TIME?OF?FIRST?OBS‘）~=0）%數據文件中第一個記錄的時刻
????????nian（11）=str2num（line（1:6））;
????????yue（11）=str2num（line（7:12））;
????????ri（11）=str2num（line（13:18））;
????????shi（11）=str2num（line（19:24））;
????????fen（11）=str2num（line（25:30））;
????????miao（11）=str2double（line（31:43））;
????end
????if?（findstr（line‘TIME?OF?LAST?OBS‘）~=0）%數據文件中最后一個記錄的時刻
????????nian（21）=str2num（line（1:6））;
????????yue（21）=str2num（line（7:12））;
????????ri（21）=str2num（line（13:18））;
????????shi（21）=str2num（line（19:24））;
????????fen（21）=str2num（line（25:30））;
????????miao（21）=str2double（line（31:43））;
????end
????if?（findstr（line‘#?OF?SATELLITES‘）~=0）
????????satenum=str2num（line（1:6））;
????end
???%?if?（findstr（line‘PRN?/?#?OF?OBS‘）~=0）
????????%讀取個衛星觀測數據的相應數據的個數
????%end
????if?（findstr（line‘END?OF?HEADER‘）~=0）
????????break;
????end
end
?fprintf（fid2‘測站近似坐標X0=%14.4f‘approxx）;
?fprintf（fid2‘?????Y0=%14.4f‘approxy）;
?fprintf（fid2‘?????Z0=%14.4f\n‘approxz）;
?fprintf（fid2‘測站觀測間隔??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????13222818??2017-06-18?13:51??程序與測試數據\abmf0800.15o
?????文件???????10321??2017-06-20?08:13??程序與測試數據\readGPSdata.m
?????目錄???????????0??2017-06-20?08:16??程序與測試數據\