function?[OBSXYZ]=?ReadObsfile（）
%此文件用來讀取RINEX格式的觀測數據文件
%讀取前先設定每顆衛星所提供的觀測參數個數
fprintf（‘開始讀取O文件\n‘）;
[FileNamePathName]=uigetfile（‘*.*o‘‘select?observation?rinex?file‘）;
fido=fopen（FileName）;
if?fido<1
????fprintf（‘打開o文件失敗\n‘）;
end

fseek（fido0‘eof‘）;
f_size=ftell（fido）;
fseek（fido0‘bof‘）;%把指針放到文件開頭
head_find=0;
while（~feof（fido）&head_find==0）
????line=fgetl（fido）;
????k=findstr（line‘APPROX?POSITION?XYZ‘）;
????if~isempty（k）
????????XYZ=[str2num（line（2:14））str2num（line（16:28））str2num（line（31:42））];
????end
????k=findstr（line‘TYPES?OF?OBSERV‘）;
????if?~isempty（k）
????????[numtypetype]=parse_type（line）;?????????%函數在最下面，用于得到觀測值的類型數量和類型
????end
????
????k=findstr（line‘INTERVAL‘）;
????if~isempty（k）
????????interval=str2num（line（1:11））;
????end
????
????k=findstr（line‘END?OF?HEADER‘）;
????if~isempty（k）
????????head_find=1;
????end
end
epoch=0;
while（feof（fido）~=1）
????epoch=epoch+1;
????clear?sat_num?num_sats;
????line=fgetl（fido）;
????
????long=length（line）;
????if（long<80）
????????line（long+1:80）=‘0‘;
????end
????
????utc_time（1）=str2num（line（1:3））;
????utc_time（2）=str2num（line（4:6））;
????utc_time（3）=str2num（line（7:9））;
????utc_time（4）=str2num（line（10:12））;
????utc_time（5）=str2num（line（13:15））;
????utc_time（6）=str2num（line（16:26））;
????utc_time;
????OBS（epoch）.t=cal2gps（utc_time）;???%%%%%%
????
????num_sats=str2num（line（30:32））;
????if（num_sats>12）
????????fprintf（‘程序不能讀取大于12課衛星的o文件‘）;
????????return
????end
????for?i=1:num_sats
????????start_dx=33+（i-1）*3;
????????sat=sscanf（line（start_dx:start_dx+2）‘%c‘）;
????????if?~isempty（deblank（sat（2:3）））
????????????sat_num（i）=str2num（sat（2:3））;???
???