//基本廣播星歷塊
struct?EPHEMERISBLOCK
//每小時一個衛星對應一個基本星歷塊
{
//PRN號?
int?i1PRN;
double?f8a0f8a1f8a2;//時間改正數
//六個軌道參數
double?f8IODEf8Crsf8Deltanf8M0;//?ORBIT?-?1
double?f8Cucf8ef8Cusf8SqrtA;//?ORBIT?-?2
double?f8Toef8Cicf8OMEGAf8Cis;//?ORBIT?-?3
double?f8i0f8Crcf8omegaf8OMEGAdot;//?ORBIT?-?4
double?f8IDOTf8GpsWeekNumberf8L2Cf8L2P;//?ORBIT?-?5
double?f8SatAccuracyf8SatHealthf8TGDf8IODC;//?ORBIT?-?6
};

struct?XYZ
{
int?prn;
double?xyz;
};



////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////*參數定義如下，在DOC類中，EPHEMERISBLOCK?*m_pGpsEphemeris;
XYZ?*m_pxyz;
int?m_i4EphemerisBlockNum;


CGpsTime為一個類，部分代碼如下



int?CGpsTime::InitgGpsTime?（int?nYear?int?nMounth?int?nDay?int?nHour?int?nMinute?double?dSecond）
{
int?DayofMonth?=?0;
int?DayofYear?=?0;
int?weekno?=?0;
int?dayofweek;
int?m;
if?（nYear??12?||?nDay??31）?return?-1;
//計算從1980年到當前的前一年的天數
for（?m?=?1980?;?m?{
???if?（?（m%4?==?0?&&?m%100?!=?0）?||?（m%400?==?0）?）?
???{
????DayofYear?+=?366;
???}
???else
????DayofYear?+=?365;
}
//計算當前一年內從元月到當前前一月的天數
for（?m?=?1;m?{
???if（m==1?||?m==3?||?m==5?||?m==7?||?m==8?||?m==10?||?m==12）
????DayofMonth?+=?31;
???else?if?（m==4?||?m==6?||?m==9?||?m==11）?
????DayofMonth?+=?30;
???else?if?（m?==2）
???{
????if?（?（nYear%4?==?0?&&?nYear%100?!=?0）?||?（nYear%400?==?0）?）
?????DayofMonth?+=?29;
????else?
?????DayofMonth?+=?28;
????
???}
}
DayofMonth?=?DayofMonth?+?nDay?-?6;//加上當月天數/減去1980年元月的6日??
weekno?=?（DayofYear?+?DayofMonth）?/?7;//計算GPS周
dayofweek?=?（DayofYear?+?DayofMonth）?%?7;
//計算GPS?周秒時間
m_dWeekSecond?=?dayofweek*86400?+?nHour*3600?+?nMinute*60?+?dSecond;
m_nWeek?=?weekno;

return?weekno;
}


*////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

////strEpheNam為廣博星歷文件，如：c:\15372370.06n
int?CGpsPosDoc::ReadBrodcastEphemeris（CString?strEpheNam）
{
int?i4HeadLineNum?=?0;
//打開文件
CStdioFile?pfEph;
BOOL?IsEn?=?pfEph.Open（strEpheNam?CFile::modeRead）;
if（!IsEn）?return?0;
//讀入頭文件
CString?strLine;

while（IsEn）
{
???IsEn?=?pfEph.ReadString（strLine）;
???i4HeadLineNum++;
???int?index?=?strLine.Find（“END?OF?HEADER“）;
???if（?-1?!=?index?）
????break;
}
//計算星歷塊數
int?AllNum?=?0;
while（IsEn）
{
???IsEn?=?pfEph.ReadString（strLine）;
???AllNum++;
}
//臨時讀入星歷塊
int?i4AllBlockNum?=?（AllNum?+?1）?/?8;
EPHEMERISBLOCK?*pTemEpheme?=?new?EPHEMERISBLOCK[i4AllBlockNum];
CGpsTime?*pGpsTime?=?new?CGpsTime[i4AllBlockNum];

if（!pTemEpheme?||?!pGpsTime）?return?0;
//將文件指針調整到數據位置
pfEph.SeekToBegin（）;
for（int?i=0;?i???IsEn?=?pfEph.ReadString（strLine）;
????//定義讀取的參數
int?mPrn;//衛星號PRNo?
int?yearmonthdayhourminute;//衛星鐘參考時刻
double???msecond;?
double???a0a1a2;//衛星鐘飄參數
double???IO

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????9411??2009-08-19?09:56??通過GPS廣播星歷計算衛星坐標\代碼.c

?????目錄??????????0??2009-08-19?09:53??通過GPS廣播星歷計算衛星坐標

-----------?---------??----------?-----??----

?????????????????9411????????????????????2