%?函數(shù)功能說明:讀取IGS精密星歷數(shù)據(jù)子函數(shù)freadnfile.m

function?[igsdata]=freadnfile（nfilename）;

%?打開IGS精密星歷文件
[fidnmessage]=fopen（nfilename‘rt‘）;
fidn;

%?開始讀取IGS精密星歷文件的內(nèi)容
disp（‘開始讀取衛(wèi)星星歷表頭文件的內(nèi)容‘）;
linen=fgetl（fidn）;

block=str2num（linen（33:39））;%歷元個(gè)數(shù)

linen=fgetl（fidn）;
linen=fgetl（fidn）;

sumsv=str2num（linen（5:6））;%衛(wèi)星個(gè)數(shù)
????
while?（feof（fidn）~=1）????
????linen=fgetl（fidn）;????
????[rowcolumn]=size（linen）;
????
????%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
????%?獲取表頭文件的結(jié)尾的指針位置sitefidn
????if?（linen（1）==‘*‘）
%??????????if?（（column<=3）&（linen（1:2）==‘/*‘））
????????sitefidn=ftell（fidn）;????%?sitefidn為文件表頭結(jié)束指針位置??
????????disp（‘讀取衛(wèi)星星歷表頭文件完畢‘）;
????????break;
????end?
????%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
end

igsdata=zeros（block5sumsv）;??%%block為衛(wèi)星在此數(shù)據(jù)時(shí)段中出現(xiàn)的歷元數(shù)，sumsv為衛(wèi)星總數(shù)，5為衛(wèi)星PRN編號（1個(gè)）、衛(wèi)星坐標(biāo)（3個(gè)）和衛(wèi)星鐘差（1個(gè)）共5個(gè)參數(shù)

for?b=1:block?????
????for?cen=1:sumsv
????????linen=fgetl（fidn）;
????????igsdata（b1cen）=str2num（linen（3:4））;?%%衛(wèi)星PRN編號
????????igsdata（b2cen）=1000*str2num（linen（6:18））;?
????????igsdata（b3cen）=1000*str2num（linen（20:32））;
????????igsdata（b4cen）=1000*str2num（linen（34:46））;
????????igsdata（b5cen）=str2num（linen（48:60））;?
????end
????linen=fgetl（fidn）;
end

igsdata;
????????
%?關(guān)閉文件
fclose（fidn）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件???????31704??2015-05-14?15:13??單點(diǎn)定位\X.xls
?????文件????????1501??2015-05-14?11:06??單點(diǎn)定位\freadnfile.m
?????文件?????????399??2015-05-13?15:49??單點(diǎn)定位\fsvposion.m
?????文件??????251250??2014-05-05?14:16??單點(diǎn)定位\igs17761.sp3
?????文件????????6482??2015-05-14?15:12??單點(diǎn)定位\kalman.m
?????文件??????146618??2015-05-13?18:48??單點(diǎn)定位\monp0200.14o
?????目錄???????????0??2015-05-18?09:15??單點(diǎn)定位\