对两种GPS解算方案的残差统计检验评估

文[1]和[2]中提出根据拟合后剩余残差的统计检验来全面评估解算,其想法是除了一般采用的ＲＭＳ和协方差外还试图以更多的统计量来反映解算更细致的差别.将其应用于Ｌａｇｅｏｓ激光资料的处理,文[3]鉴别了具有相同ＲＭＳ水平而差异较大的不同解算.但是除了理论上的依据外,文[3]事实上未能表明其中残差分布较合理的解确实优于其余的情况.本文的目的之一是试图改进这一不足之处.采用由美国喷气推进实验室(ＪＰＬ)发展的ＧＰＳ数据处理软件ＧＩＰＳＹ,本文比较仅固定精密星历(ＦｉｘｅｄＰｒｅｃｉｓｅＯｒｂｉｔｓ,简称为ＦＰＯ)和同时固定精…     

多种资料综合处理的精度分析和相对权因子

在统计轨道理论的框架下分析多种技术资料的综合处理方法的精度．提出除了拟合剩余残差的ＲＭＳ外,皮尔逊（Ｐｅａｒｓｏｎ） χ２ 统计量和高阶矩也应作为判断解算好坏的重要指标．初步探讨了利用这些统计量评价不同技术间的相对权因子给定问题．作为对系统差来源的初步探索,利用考察协方差分析理论分析了考察参数误差对解算的影响,讨论了考察协方差与相对权因子的关系     

月球天平动模型对激光测月解算UTO的影响

本文利用1987年1月－1989年3月的LLR资料解算了UTO－UTC序列，在归算中分别应用了Eckhardt-500和LLB200两个天平动模型，以便分析和比较它们对解算UTO的影响。结果表明两个天平动模型间存在着明显的系统差。它反映在两者的事前残差中，最大差值可达100ns，这就使得用Eckhardt-500天平动模型的事前残差解算得到的UTO－UTC序列对IERS的平均偏离比采用LLB200天平动模型时要大1ms，但对两结果序列的外部精度没有影响。用事后残差解算，系统差的影响减小了，确定的UT0－UTC序列更靠近IERS系统，用两种天平动模型解算的UTO－UTC结果也更接近，外部精度也几乎一样且双用事前残差解的结果略好。     

IGS‘92联测期间地极坐标的高频波动

本文分析了ＩＧＳ＇９２联测期间七个ＧＰＳ数据处理中心提供的极坐标序列。通过谱分析、最小二乘拟合和Ｆ检验，表明在这些序列中存在一些共同的高频波动：在Ｘ方向上具有２７．０，１６．５，１３．４和１０．４天的周期，在Ｙ方向上的波动周期约为２０．５，１５．８和１０．０天，并且每个序列与ＥＯＰ（ＩＥＲＳ）９２Ｃ０４之间都存在一个系统差。计算与分析表明，这些系统偏离的主要原因是由于在用ＧＰＳ资料解算Ｘ、Ｙ时，不     

IGS′92联测期间地极坐标的高频波动

本文分析了ＩＧＳ＇９２联测期间七个ＧＰＳ数据处理中心提供的极坐标序列。通过谱分析、最小二乘拟合和Ｆ检验，表明在这些序列中存在一些共同的高频波动：在Ｘ方向上具有２７．０，１６．５，１３．４和１０．４天的周期，在Ｙ方向上的波动周期约为２０．５，１５．８和１０．０天。并且每个序列与ＥＯＰ（ＩＥＲＳ）９２Ｃ０４之间都存在一个系统差。计算与分析表明，这些系统偏离的主要原因是由于在用ＧＰＳ资料解算Ｘ、Ｙ时，不同分析中心采用了不同系统的台站坐标（或者说只有部分台站采用了固定的台站坐标），从而造成这些序列所在的参考架与ＩＴＲＦ９１之间存在一个平移和旋转。最后，计算了该期间的大气角动量激发函数，可部分地解释该期间的Ｘ、Ｙ高频波动的原因。     

利用非差资料的精密点定位方案解算区域GPS网

单站精密定位（Precise Point Positioning，以下简称PPP）是在同时固定GPS精密星历和卫星钟的前提下，利用载波相位和伪距资料进行单台站的精密点定位，采用该方法时不同台站之间不存在共同的待估参数，即各台站互不相关，这一特点大大降低了计算量，采用美国喷气推进实验室JPL发展的数据处理软件GIPSY处理APSG联测资料，计算表明PPP的重复率相当于目前国内普遍采用的双差解算结果，采用较好保持地面网构型的无基准算法，计算表明通过Helmert参考系转换后，PPP的解算结果与双差算法的外符精度大致相当，解算表明，采用PPP处理100个台站约需3．5小时，而处理同样的资料采用双差算法则需18－20小时，对于我国即将建成的大科学工程或地震监测的多达2000个接收机的GPS网而言，在保持精度前提下的节省计算资源和计算时间的PPP解算方案值得广泛的应用。     

对EOP（GSFC）96R 01中极偏移序列的分析

分析了１９７９年８月至１９９６年１月之间由ＭａｒｋⅢ ＶＬＢＩ时延资料得到的ＩＡＵ１９８０间动模型在经度（δψ）和交角（δψ）方向的偏移序列,得到了对ＩＡＵ１９８０章动模型中主项系数的改正,并检测了自由地核章动。计算步骤依照常规资料处理方法设计,并采用了一些分析技巧以便充分利用资料的整个时段和减小端部效应对参数解的影响。对拟合残差序列的进一步分析表明,自由地核章动是统计显著的,但其参数尤其是其振幅     

利用非差资料的精密点定位方案

单站精密定位 (PrecisePointPositioning ,以下简称PPP)是在同时固定GPS精密星历和卫星钟的前提下 ,利用载波相位和伪距资料进行单台站的精密点定位 .采用该方法时不同台站之间不存在共同的待估参数 ,即各台站互不相关 ,这一特点大大降低了计算量 .采用美国喷气推进实验室JPL发展的数据处理软件GIPSY处理APSG联测资料 ,计算表明PPP的重复率相当于目前国内普遍采用的双差解算结果 .采用较好保持地面网构型的无基准算法 ,计算表明通过Helmert参考系转换后 ,PPP的解算结果与双差算法的外符精度大致相当 .解算表明 ,采用PPP处理 1 0 0个台站约需 3 .5小时 ,而处理同样的资料采用双差算法则需 1 8～2 0小时 .对于我国即将建成的大科学工程或地震监测的多达 2 0 0 0个接收机的GPS网而言 ,在保持精度前提下的节省计算资源和计算时间的PPP解算方案值得广泛的应用     

1979-1999全球天测与测地VLBI观测资料处理结果与讨论

采用国际地民服务（ＩＥＲＳ）１９９６年规范和上海天文台甚长基线干涉测量（ＶＬＢＩ）资料分析软件,对自１９７９年８月以来全球约２９００次天测与测地ＶＩＢＬ实验进行了综合归算,得到６０１颗河外射电源在Ｊ２０００．０的坐标、１０１台站在Ｊ１９９７．０的坐标和运动速度、以及约２０年地地球定向参数序列。与国际相关结果的比较表明,归算结果不存在显著系统差。通过对归算结果的分析,指出目前国际上评估天球参考架局部     

恒星统计视差研究方法的新进展

本文对按经典方法确定恒星统计视差这一重要恒星天文学问题作了简单的回顾,并指出这种方法在统计学上是不严格的。太阳运动和速度椭球的确定严格来说是不可分离的,长期视差和统计视差不能看作为彼此独立的量。在上述讨论的基础上,对于按最大似然原理解算包括统计视差在内的运动学参数问题的提出和发展过程作了介绍,重点介绍了Clube等和Murray在七十年代和八十年代所提出的方法。Murray所用的数学模型认为星群内各个恒星的残差v服从三维正态分布,因而自然引入本动协方差阵。就目前来说,它可算是最为严格的方法,其中共需通过迭代算法解算11个未知数。最后,对这两种方法作了简单的比较。 统计视差的确定是恒星天文学的重要内容之一,在有关宇宙距离尺度问题的研究中它有着不可忽视的地位。比如,按统计视差法求得的天琴座RR型变星的绝对星等,可用来确定球状星团以至较近河外星系的距离,从而被誉为“量天尺”。因此,关于统计视差确定方法的研究长期来受到人们的重视,并一直处于不断改进之中。本文试图对这一工作的主要方面作一综合性的评述,并着重介绍最近一些年代提出的新方法。     

引力常数变化对地球自转长期变化的影响

探讨和估计了各种引力常数变化理论对地球角速度和日长变化的影响。各种引力常数变化理论包括了引力常数G随时间、空间以及速度变化等几个方面的影响。另外也估计了对地球自转角速度和日长变化产生的效应。其中有些研究对探讨地球自转变化也有启发意义。     

近地小行星轨道演化的研究

本文主要阐述了近年来我们在近地小行星轨道演化研究所获得的一些基本结果,即合理的力学模型和相应的有效算法,并以实际预报算例与有关权威性的结果作了比较,证实这些研究结果确实是可信的。     

表面边界位置对太阳ρ模本征振荡频率的影响

本文研究了表面边界位置对太阳ρ模绝热本征振荡频率的影响。数值计算结果表明，对于ｖ４０００μＨｚ的中低阶ｐ模，表面边界置于温度极小点所引起的本征振荡频率的计算误差随着ｖ和ｌ的增大，表面边界点位置对太阳理论振荡频率的影响增大，色球结构对太阳ｐ模振荡频率的影响已变得不可忽略。     

关于土卫八运动的轨道解

为了适应星际探测的需求,本文建立了在新的精度要求下土星卫星运动对应的力学模型,具体讨论了土卫八的运动,并针对主要摄动源土卫六的引力作用,建立了轨道变化的分析解,以此表明建立了土卫运动理论该采取的途径和精密定轨宜采用以轨道根数作为状态量的数值定轨方法。     

对钱德勒摆动椭率及长轴方向的检测

证明地球的三轴性对钱德勒摆动椭率的影响只有１．８Ｅ５,,可以被忽略．假定极潮为平衡潮,求出由于海陆分布的影响,将使钱德勒摆动有０．０１３的椭率,主轴方向大约指向东经６°．采用傅里叶反卷积方法对ＳＰＡＣＥ９２和ＢＩＨ的极移资料进行分析,求出逐年的长期极移、周年极移和钱德勒摆动的参数,然后将钱德勒摆动归算到同一历元,证明钱德摆动椭率约为０．０１１,长轴方向约为东经９°,与理论结果十分符合．     

埃尔塔宁小行星质量的下限估计

埃尔塔宁小行星是一颗２１５ 万年前在南太平洋撞击地球的小行星,它也是目前唯一的深海区与地球碰撞的小行星。根据埃尔塔宁号和Ｐｏｌａｒｓｔｅｒｎ 号海洋调查船在陨击海区勘查的小行星残骸分布密度,我们估计埃尔塔宁小行星质量为２ ．３ ×１０１３ｇ ,对应于它的直径为０ ．２３ 公里。这比Ｋｙｔｅ 等人（１９８８） 估计的质量小了１０ 倍,然而它更接近埃尔塔宁小行星撞击事件的实际情况,可以作为埃尔塔宁小行星质量下限和半径下限更合理的估计值。     

日长变化预报中BP神经网络拓扑结构的选择

日长变化的预报具有重要的科学意义和实际应用价值。非线性的人工神经网络技术中的反向传播模型（BP网络）可用于预报日长变化。BP网络的拓扑结构决定了神经网络解决问题的能力，针对不同的问题需要采用不同的网络结构。该文分析了神经网络的拓扑结构算法，选用最小均方误差法确定网络的拓扑结构，并将此应用于日长变化预报。结果表明，该方法是可靠和有效的。     

银河系运动学参数的确定

利用 2 0 2个太阳附近疏散星团的视向速度和自行观测资料 ,对太阳的运动和银河系的运动学参数进行了研究。其中 ,距离在 0 .5kpc到 2kpc之间的 12 8个疏散星团对平均太阳运动分量的解算结果是 (u0 ,v0 ,w0 ) =(- 13.8± 1.4 ,- 5 .0± 1.6 ,- 11.6± 2 .9)km/s ;Oort常数和银河系径向运动参数的解算结果分别为 (A ,B) =(16 .9± 1.1,- 11.6± 2 .6 )km·s- 1·kpc- 1及 (C ,D) =(2 .5± 1.1,- 2 .1± 0 .9)km·s- 1·kpc- 1。