%%?提取觀測值文件（O文件）里的偽距，北斗和GPS
%%?格式為每顆衛星的數據按列以時間順序排列
%%?得到的周內時是GPS時間




function?[??Obs??]?=?Extra3_02O（?filename??mark?）

%------------定義常量和字符-----------
line_num???????=?0;???????????????????%行數
markEnd?=?‘END?OF?HEADER‘?;
%------------------------------------
fid??=?fopen（filename‘r‘）;
%-----------------------
while?feof（fid）~=1?
????%讀取頭文件
????tline????????=?fgetl（fid）;??
????chartline?=?char（tline）?;???
????END?=?strfind（?chartline??markEnd?）?;
%?????schar??????=?chartline（3:12）;???????????????????????????
????if??isempty（?END?）?
????????continue?;
????else
????????break?;
????end
end
while?feof（fid）~=1
????????line_num?=?line_num+1;???????????????????????????%?歷元計數
????????tline????????=?fgetl（fid）;??
????????chartline?=?char（tline）?;
????????%------------------數據觀測時間----------------
????????time?????=??chartline（3:21）;?????????????????????????????????????%存儲時間
????????UTC??=?sscanf（?time??‘%d%d%d%d%d%d‘??[16]?）?;????????????%轉換為雙精度時間信息，UTC時間
?????????[Toc_WN?Toc_SOW]?=?utc2gpstime（UTC?）?;??????????????????????%將toc從UTC格式轉為GPS時格式（周數和周內時）
?????????Obs.ObsTime.Toc_WN（line_num）?=?Toc_WN;????????????????????????????????%周數
?????????Obs.ObsTime.Toc_TOW?（line_num）?=?Toc_SOW;?????????????????????????%周內時
?????????
????????search_SatNum?=?str2num（chartline（34:35））;???????????????????????????%搜索到的衛星總數
%?????????Rec_clock（?line_num1?）?=?str2num（chartline（42:56））;?????????%接收機時鐘校正?
????????%-----------------------------------
????????NumPRN1?=?0?;
????????NumPRN2?=?0?;
????????for?Nr?=?1:search_SatNum
????????????chartline?=?char（fgetl（fid））;?
????????????if??chartline（1）?==?mark.gnss1
????????????????NumPRN1?=?NumPRN1+1?;
????????????????satNr?=?str2num（chartline（2:3））?;?
????????????????Obs.BD（line_num?）??.PRN（NumPRN1）?=???satNr??;
????????????????Obs.BD（line_num?）?.PseudoRange（NumPRN1）?=?str2num（chartline（6:17））;
????????????????Obs.BD（line_num?）?.phase（NumPRN1）?=?str2num（chartline（21:33））;
????????????elseif?chartline（1）?==?mark.gnss2
????????????????NumPRN2?=?NumPRN2+1?;
????????????????satNr?=?str2num（chartline（2:3））?;?
????????????????Obs.GPS（line_num?）??.PRN（NumPRN2）?=???satNr??;
????????????????Obs.GPS（line_num?）?.PseudoRange（NumPRN2）?=?str2num（chartline（6:17））;
????????????????Obs.GPS（line_num?）?.phase（NumPRN2）?=?str2num（chartline（21:33））;
????????????end
????????end
end
Obs.epoch?=?line_num?;
fclose（fid）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2565??2016-02-22?17:05??Extra3_02O.m