资源都在这里了所需卫星轨道计算

GPS种与无关的RINEX(TheReceiverIndependentExchangeFormat)格式，目前已使用2号版本。下面分别介绍RINEX2格式的广播星历文件、观测数据文件、和地面气象数据文件。ssssdddf.其中：ssss～以4f～该站该日收到的某类文件的顺序号，0t～文件种类：为了便于交流，RINEX2格式的GPS数据文件均以①无带标；②ASCII码；③每个记录长度为80个字符，块大小为8000；录制在磁带上，磁带上的第一个文件是全部文件的。但目前国际上的IGS等组织是通过通讯方式表2.1.1RINEX2 NAVIGATIONDATA RINEXVERSION/TYPE Version1.09 29-NOV-9521:05 PGM/RUNBY/DATE .2235D-07-.1192D-06-.1192D- ION .1310D-07-.1310D+06- ION.133179128170D-06.107469588780D- LEAPENDOFHEADER9941021800.0-0.103851780295D-04-0.909494701773D-120.000000000000D+000.720000000000D+020.106062500000D+030.476841277575D-080.548548996449D-050.312971079256D-020.747293233871D-050.233938313166D-090.100000000000D+010.771000000000D+030.700000000000D+010.000000000000D+000.139698386192D-08.0-0.635907053947D-04-0.909494701773D-120.000000000000D+000.228000000000D+030.167187500000D+020.424946255961D-080.566244125366D-060.790101150051D-020.883266329765D-050.404659722397D-090.100000000000D+010.771000000000D+030.700000000000D+010.000000000000D+000.139698386192D-08上表中的前几行为表头，表头中的第60～80个字符是相应行的说明，如第一行的“2”是RINEX版本号，COMMENT中的第八行至第十五行是PRN9的广播星历，对应量的含义如表2.1.2。表2.1.2RINEX2PRN号，年，月，日，时，分，秒，af0，af1，afaode，Crs，M0Cuc，e，CusAtoe，Cic，0i0，Crc，svacc，svhlth，tgdaf0～钟差常数项，时秒af1～钟差漂移项，时秒/时秒 2 fn～平近点角的长期变化(近地点参数)，弧度/M0～e～A～长半轴的平方根，米1/20～i0～～～升交点赤经在赤道平面中的长期变化(主要是由地球引力场的2阶带谐系数C20J2)引起)，弧度/Cuc～在星历参考时刻toe在轨道延迹方向上周期改正余弦项的振幅，弧度Cus～在星历参考时刻toeCrc～在星历参考时刻toe在轨道径向方向上周期改正余弦项的振幅，米；Crs～在星历参考时刻toeCic～在星历参考时刻toe轨道倾角(近似于法向)周期改正余弦项的振幅，弧度Cis～在星历参考时刻toe轨道倾角(近似于法向)周期改正正弦项的振幅，弧度toe～星历参考时刻(星期中的秒数)，秒；aode～星历数据的cflgl2～l2上存在哪些码指示weekno～G星期数；pflgl2～l2上Psvacc～svhlth～是否健康指标tgd～aodc～钟数据ttm～信息传送时间(与对接收到的有关)，时秒表2.1.3RINEX2 OBSERVATIONDATA RINEXVERSION/95-07-2022:20:20PGM/RUNBY/TurboSIIrinexMODE:1MARKERTURBOProduction REC#/TYPE/TURBOANT#/ APPROXPOSITION ANTENNA:DELTAOriginalslant #/TYPESOFENDOF95719021 0527281531------021 0527281531----后一个观测历元时刻、测得的号码，说明为“COMMENT”的注解行等。 数 数 号码 号码，...，钟当质量标记(EventFlag)为：0～1～一历元和本历元之间停电3～在新的台站重新开始观测(紧接着会出现新的点号5～其它外部6～②号码计为s：系统G=-钟差伪距=测得伪距-钟差×光速相位=测得相位-钟差×频率当失周指示(取值范围是0～7)为0或空格时数据正常。0或空格～1～5～9～表2.1.4RINEX2 METEOROLOGICALDATA RINEXVERSION/195-07-20PGM/RUNBY/MARKER3#/TYPESOFENDOF701071107210二．由广播星历计算位A

A2计算平角速度n0式中的GMtkt改正平角速度nn0MkM0n按下式迭代计算偏近点角EkMkEkesin由下两式计算真近点角

cosEk1ecosEk1e2sin

kvkr

Ccos2

ukk计算改正后的向径ikioiki*计 在轨道平面内的坐标x y' x'rcosu y'rsinu e ek e ek式中的e最后计

Zk x'cosy'cosisin Yx'siny'cosi Zy'sini 根据以上过程计算位置时，应注意，求出坐标表示在地固系中，计算时刻t必须是以GPST计的信号发射时刻，它由信号接收时刻减去时间延迟(需迭代求得，见后)后得到。当tkttoe大于302400秒时，因在tk中减去一个GPS星期相应的秒数604800秒，当tkttoe小于0秒时，应在tk中加上604800秒。CTS中。GPS观测值是和振荡器频率及时间延迟的函数时间延迟包括几何时间延迟g和大气时间延迟a，即：g 我们进行GPS测量需要求定的参数均包含在g严格来说，g的求定应该在广义相对论框架下的以系质心为原点的坐标系中计算。但实用时，在地心惯如果在地固坐标系中计算在tj时刻台站j收到 i的信号的几 时间延迟为g，按电磁 定律g

RStiRRSti是在信号发射时刻tiR是台站在地固坐标系中的坐标；C是电磁波在真空中的速度(光速，米/秒)。由于包含在2.2.16式中的信号发射时刻ti是未知的，要在计算出之后才能求出，因此需要一个迭代过程才能求得g，具体步骤是：i①按第k次迭代求得的几何时间延迟值i

gi 信号发射时刻tkgikk ijtktkij式中的大 时间延迟a可按大气改正模型计算(详见后)② 坐标Rk gg kSg kSgC 以上 中的(NGS格式的精密星历)，而又不进行定轨，则选择在地固系中计算较为方便；如果要求定轨据量非常大，我们一般只采用高度角大于15的观测数据，所以大部分的现有天顶改正和映射函数模型均可采用。经常采用的SAASTAMOINEN1973模型如下：dzddmapwzd dzd

f wzd

f 70erh6.1110t02370由于4.1式中的湿项模型，一般可引进一个待定参数xatm，使4.1式变为 dzddmapwzdxatm

mz

cosz

actgz cosz00cp0为地面气压(毫巴)，t0为地面温度(摄氏)，e0为水气压(毫巴)，rh为相对湿度，为台站的地心纬度，h为以TA表示水汽辐射计测得的任意方向上的天空亮温度，则大气湿项引起的电磁波路径增长L(即前面a与光速的乘积)可由以下一 LC0C1T'Af1C2T'Af C

1

T' C

bgf22k

1f21

0W1 f2W1 C2

2f2W2

T TLn1T

T

af

v1v1f2f2 f2f2 a afT'

T

0101ff ff

上式中的

和f

为两个波段的频率，TAf

TAf

可近似地取为2.8 T和T

ke

(tK0.92ke0.95)。水汽辐射计测得的L的精度可达到1.5电磁波通过电离层产生的时间延迟f为：2f

式中的f为观测频率，C为光速，Ne为路径上的电子总含量，r0为经典电子半径由于4.1.5式中的Ne与辐射压有关，昼夜可相差一个数量级，很难用模型来精确表示。但可以看出它产生如假设1、2、1、2为某 在同一历元测得的L1和L2两个波段上的伪距和相位，则电离层对

f f1

f f2

、

，这里的CONST对一个对同一在同一线性组合后的观测量LC作为观测值。

1式中的g

f2，即L1与L2基准站R测得至GPSJ的伪距j'j dj)djj 其中j为基准站R至第J颗之间的真实距离；dj是GPS星历误差所引起的距离偏差；d为 收机时钟相对于GPS时间系统的偏差；dj是J颗时钟相对于GPS时间系统的偏差；j电离层时延所引 的距离偏差；j是对流层时延所引起的距离偏差；c是为电磁波的速度

jj'

dj)dj

在基准进行伪距测量的同时，流动站K也对第J颗进行了伪距测量，流动站所测得j'jc(ddj)djj djdj,(jj,j j'jjc(dd