%本程序用于將subfram 1中得到的導(dǎo)航信號(hào)轉(zhuǎn)化為星歷表信號(hào)，
ephermeris data
clear angmat navd rymkmat
load d:\gps\big data\srvy8210;
angmat=ang;
ptnumat=transpt;
rymkmat=fintime;
[st navd]=navdat(angmat); %*find naviga data (20ms) & first
angle transition
[points, navdata]=matsubf(navd); %尋找subfram 2的電文頭
prod 0 1s
nsat=length(points);
for m=1:nsat;
iode2(m,:)=navdata(m,points(m)-1+61:points(m)-1+68);%獲得星歷數(shù)據(jù)期號(hào)域，IODE由IODC的8個(gè)最低有效位(LSB)組成，可以被用戶設(shè)備用于檢測(cè)廣播軌道根數(shù)
crs(m)=comp2dec(navdata(m,points(m)-1+84:-1:points(m)-
1+69),-5);%軌道向徑正弦調(diào)和改正的振幅
deln(m)=comp2dec(navdata(m,points(m)-1+106:-1:points(m)-
1+91),-43)*pi;%獲得平均運(yùn)動(dòng)角速度的校正量，乘以PI轉(zhuǎn)化為弧度制表示
m01=navdata(m,points(m)-1+114:-1:points(m)-1+107);
m02=navdata(m,points(m)-1+144:-1:points(m)-1+121);
m03=[m02 m01];
m0(m)=comp2dec(m03,-31)*pi;%參考時(shí)刻的平近點(diǎn)角,乘以PI轉(zhuǎn)化為弧度制表示
clear m01 m02 m03;
cuc(m)=comp2dec(navdata(m,points(m)-1+166:-1:points(m)-
1+151),-29);%升交角距余弦調(diào)和校正的振幅
e1=navdata(m,points(m)-1+174:-1:points(m)-1+167);
e2=navdata(m,points(m)-1+204:-1:points(m)-1+181);
e3=[e2 e1];
e(m)=bi2de(e3)*2^(-33);
clear e1 e2 e3;%衛(wèi)星軌道的離心率
cus(m)=comp2dec(navdata(m,points(m)-1+226:-1:points(m)-
1+211),-29);%升交角距正弦調(diào)和校正的振幅
sa1=navdata(m,points(m)-1+234:-1:points(m)-1+227);
sa2=navdata(m,points(m)-1+264:-1:points(m)-1+241);
sa3=[sa2 sa1];
sa(m)=bi2de(sa3)*2^(-19);%獲得軌道長(zhǎng)半軸的平方根
clear sa1 sa2 sa3;%清除中間變量
toe(m)=bi2de(navdata(m,points(m)-1+286:-1:points(m)-
1+271))*2^4;%星歷參考時(shí)刻，從GPS時(shí)每星期六/星期日子夜零時(shí)起算，變化于0～604800秒
end