基于SSA-Sn法的GNSS坐标时间序列粗差探测

针对GNSS高程坐标时间序列中含有粗差、高程坐标时间序列不平稳、非线性变化导致异常值探测困难等问题,提出一种奇异谱分析SSA与S_n估计量相结合的粗差探测方法。利用实验模拟数据和实测数据分别对该方法的探测效果同常规3σ法和MAD法进行对比,结果表明,SSA-S_n粗差探测方法的总体探测效果更好,能够满足GNSS坐标时间序列高可靠性的需求。     

基于DIA的GNSS坐标时间序列异常值探测

为解决GNSS坐标时间序列中存在异常值导致地壳形变分析受到影响的问题,首先将经典数据探测法DIA(detection, identification, adaptation)引入GNSS坐标时间序列分析中,基于谐波模型构造相应的统计量;然后采用计算效率更高的改进DIA算法进行异常值探测;最后采用模拟数据和中国大陆地壳运动观测网络(陆态网络)CMONOC实测数据对DIA算法和传统3σ法进行比较。结果表明,DIA算法能够有效探测出GNSS坐标时间序列中的异常值,准确率高于传统3σ法。相较于传统DIA算法,改进后的DIA算法计算效率显著提高。     

城市CORS坐标时间序列噪声模型建立与分析

以GAMIT/GLOBK处理得到高精度常州CZCORS坐标时间序列,用奇异谱分析方法对其进行空间滤波,使用最大似然估计方法对剔除共性误差的坐标时间序列进行噪声特征分析,得出最优白噪声模型为“WH”,N、U分量最优有色噪声模型为“WH+FN”,E分量最优有色噪声模型为“WH+FN+RWN”。根据最优有色噪声模型估计出的CZCORS站速度场表现为东南方向整体运动的趋势,其整体位移变化值及运动方向与周围IGS站基本一致。     

环境负载对GNSS坐标时间序列噪声特性的影响研究

基于中国大陆构造环境监测网络的连续GNSS观测数据,选取云南省25个测站的坐标时间序列研究地表质量变化对建立噪声模型的影响。计算陆地水负载、大气压负载和非潮汐海洋负载引起的测站非线性变化,采用极大似然估计对不同噪声组合进行分析,分别建立负载改正前后各GNSS测站的最优噪声模型,同时估计速度的不确定度。结果发现,云南省GNSS测站各分量的噪声特性主要表现为闪烁噪声+白噪声(FN+WN);环境负载改正能够改变坐标时间序列的最优噪声模型;用加入环境负载改正后建立的最优噪声模型估计各测站的速度不确定度,对于E、N、U分量,分别有72%、80%、68%的测站的速度不确定度变小。     

顾及GNSS坐标时间序列中季节信号的CEEMD降噪方法

针对互补集合经验模态分解（complementary ensemble empirical mode decomposition,CEEMD）在GNSS坐标时间序列降噪过程中需要选取筛选准则的问题,提出一种顾及GNSS坐标时间序列中季节信号的CEEMD降噪方法。首先利用CEEMD方法对GNSS坐标时间序列进行分解,然后计算各本征模态函数（intrinsic mode function,IMF）的平均周期,将平均周期低于120 d的IMF分量作为噪声分量扣除,并重构剩余分量为信号分量。利用该方法对中国大陆227个GNSS垂向坐标时间序列进行降噪,并与以连续均方误差和相关系数为筛选准则的CEEMD降噪方法的结果进行对比。结果表明,本文方法未出现过度降噪,而其他2种方法均导致部分测站的季节信号被扣除;在未过度降噪站点,本文方法GNSS坐标时间序列的RMS、幂律噪声和速度不确定度的平均改正率分别为19.13%、88.29%和86.46%,优于其他2种方法。     

GNSS形变监测时间序列的粗差数据处理方法

针对GNSS形变监测的传统粗差探测算法性能受限于数据长度的问题,提出基于小波分析的一阶导数粗差剔除法。首先使用新算法对原始信号进行粗差剔除,再利用广义延拓插值补充残缺点。实测算例结果表明,新算法执行时间是3倍中误差法（3σ法）的0.01倍,插值后与真实值偏离仅0.03 mm,远优于其他传统粗差探测手段。     

LOESS用于GNSS垂向坐标时间序列噪声识别与提取

针对GNSS垂向坐标时间序列噪声复杂、精度较差等特点,采用局部加权回归模型(locally weighted regression, LOESS)对中国大陆构造环境监测网络中289个GNSS站的垂向坐标时间序列进行降噪分析。首先利用LOESS方法对预处理后的时间序列进行降噪处理,得到降噪后的时间序列及噪声序列;然后采用Durbin-Watson(DW)检验对降噪后的噪声序列进行自相关性检验,同时采用Wilcoxon秩和检验方法对降噪前后序列的标准差、噪声项、速度不确定度等指标进行显著性检验;最后采用降噪前后序列的信噪比及前3个指标来定量评价降噪效果。结果表明,各测站降噪后的噪声序列不存在自相关性,采用LOESS方法降噪处理后各评价指标均有显著改正,表明LOESS方法能够有效减少GNSS垂向坐标时间序列中的噪声,进一步提高GNSS垂向坐标时间序列的精度。     

2015年尼泊尔MW7.8地震近场GNSS时间序列分析

收集尼泊尔地震近场密集的GNSS测站数据,采用GAMIT/GLOBK软件包解算2011～2017年尼泊尔近场58个GNSS连续观测站和41个震后GNSS流动观测站数据,对获得的数据进行周期性、共模误差以及噪声的时间序列分析。结果显示：1)GRACE改正GNSS时间序列周期项,可使E、N、U三个分量的WRMS平均降低15.52%、26.41%和45.06%;2)针对时间序列中存在大量缺失数据,改进传统的PCA方法,利用第一主成分计算共模误差,3个分量有效率分别达到61.91%、53.91%和53.88%;3)噪声分析表明,尼泊尔地震近场GNSS站点最佳噪声模型为WN+FN,相比单纯的WN模型,WN+FN模型估计震间速度、同震和震后位移误差分别约为其6～9倍、4～7倍和2～4倍。该研究结果对区域坐标框架维护、尼泊尔地震的同震破裂和震后余滑的地球物理模型反演具有重要意义。     

一种基于MAD改进的GNSS高程时间序列粗差探测方法

针对GNSS高程时间序列中不可避免地含有粗差,以及在非线性、不平稳性的高程时间序列中粗差难以探测的问题,在传统MAD方法基础上,构建一种引入小波分析的WT-MAD粗差探测方法.利用模拟数据和LHAZ、BJFS、TWTF三个IGS站的实测高程数据进行实验,将WT-MAD法与基于最小二乘的3σ法、IQR法和MAD法进行对比...     

顾及三维方向的GNSS坐标时间序列粗差探测方法

针对GNSS坐标时间序列中粗差为三维位置粗差的特点,提出综合三维坐标的迭代最小二乘粗差探测方法。通过实验模拟数据,分别利用3倍中误差（3σ）准则和四分位距（interquartile range, IQR）准则对单方向和综合N、E、U三方向的粗差进行探测。结果表明,2种判别准则中,综合考虑3个方向的探测效果远好于单方向的探测效果,且探测效率均超过95%;3σ准则的探测效率稍高于IQR准则,但误判率也较高。在实测数据中的应用表明,SOPAC(Scrips Orbit and Permanent Array Center)提供的部分raw类型全球IGS站时间序列中粗差占0.1%~4.7%;运用IQR准则综合3个方向的粗差探测得到的时间序列要比cleaned trend类型更加干净。     

����GA�Ľ������񶯶�Ƶ�ʳɷֱ�ʶ

????????????????????????????????????????????????????????????????????????????·???????GA?????????????GA??????????????????????????????????????????????????????????????????????? ???????????????????????????????齨??????????????????????????????????????????????????????б?????????????GA??????????????????     

一种基于M估计的Robust-ELM滑坡位移预测方法

采用传统ELM算法进行滑坡位移预测时,其网络输出权值由最小二乘估计得出,导致ELM抗差能力较差,从而造成网络训练参数不准确。为此,将M估计与ELM相结合,提出一种基于M估计的RobustELM滑坡变形预测方法。该方法利用加权最小二乘方法来取代最小二乘法计算ELM输出权值,以减少滑坡监测数据中粗差对ELM预测的干扰。分别以链子崖、古树屋滑坡体为例,将Robust-ELM进行了单维、多维粗差的抵御性验证。结果表明,该方法能够有效降低粗差对预测的影响,具有良好的抗差能力。     

基于改进Cao算法的奇异谱分析法及其在北斗多路径去噪中的应用

研究了基于奇异谱分析的北斗恒星日滤波算法,采用相空间重构Cao算法来确定奇异谱嵌入维度,并针对Cao算法的不足进行改进,提高奇异谱分析法的准确性和计算效率。分析北斗系统不同星座卫星的轨道重复周期特性,通过计算确定北斗系统多路径误差的周期约为86 160s。利用奇异谱分析法和传统小波分析法对北斗短基线解算结果进行恒星日滤波处理。结果表明,本文提出的奇异谱分析法多路径滤波效果优于小波滤波法,能较好地消除原始坐标序列中的多路径误差。     

基于GPS高程时间序列粗差的抗差探测与插补研究

研究抗差估计和稳健估计粗差探测IQR准则的粗差剔除算法以及缺值点插值的正交多项式最小二乘拟合方法。以中国3个CORS基准站从1999-03-01—2009-12-31日的单天解GPS高程时间序列为对象进行粗差探测与插补的试验。结果表明：抗差估计比IQR准则能更有效地探测和剔除粗差;三次样条插值当连续缺失历元超过5~10个时插值结果易发生扭曲;低次正交多项式拟合的结果稳定,一般选择2~3次拟合即可获得理想的效果,当连续缺失历元达3个月时,正交多项式拟合后的曲线与序列趋势仍基本一致。     

一种基于小波分析的卫星钟差数据粗差处理方法

针对传统MAD方法在探测钟差粗差方面存在的不足,提出一种基于小波分析的钟差数据粗差探测与处理方法。该方法利用小波变换的多尺度分析能力,将含有粗差的数据分解为低频小波系数和各层的高频小波系数,并在不同时间尺度下进行粗差探测和消除。使用CODE提供的BDS精密钟差数据进行实验,分析不同小波函数及不同分解尺度对预处理效果的影响,并与传统MAD方法比较,发现利用小波分解预处理后的数据在钟差预报方面有明显优势,预报精度平均提高10%。     

����Kullback-Leiber����ֲ�̽���Bayes����

????????????????Kullback-Leiber???????????????Bayes?????????????????????Σ?????????????￡??????????????????????????????????????????????????3??????????Kullback-Leiber????????????????????????б???????÷???????????????????????????Ч????     

基于改进的经验模态分解和SSA联合算法的GPS坐标时间序列重构

基于ICEEMDAN算法无需先验信息即可准确分离和提取低频信号与趋势信息的特性,以及SSA具有较好的信号重构优势,提出基于ICEEMDAN和SSA的联合重构方法。该方法将弱周期信号利用ICEEMDAN方法进行提取与重构,可弥补SSA方法中当弱周期信号对应的Hankel矩阵的奇异值和噪声Hankel矩阵的奇异值接近时容易被噪声掩盖而难以提取的不足。通过模拟实验和真实站点数据验证该算法分解和重构精度,并与奇异谱分析法、小波分解法、滑动最小二乘法进行比较。实验结果表明,ICEEMDAN-SSA联合算法相对于已有方法具有更好的重构精度。     

HVF������GPS��·��ЧӦ�о��е�Ӧ��

????????Helmert????????????Vondrak????????????HVF??????????????GPS??·??Ч????о????????????????????????????÷???????????Vondrak?????????????????????????????????Ч????????????????е????????????????÷????????·????????????GPS??·??Ч?????????????????????Ч????????·??Ч??????????????????GPS??λ?????