基于奇异谱分析的BDS卫星钟差周期项提取.pdf
《基于奇异谱分析的BDS卫星钟差周期项提取.pdf》由会员分享,可在线阅读,更多相关《基于奇异谱分析的BDS卫星钟差周期项提取.pdf(9页珍藏版)》请在咨信网上搜索。
1、第 31 卷第 9 期 中国惯性技术学报 Vol.31 No.9 2023年09月 Journal of Chinese Inertial Technology Sep.2023 收稿日期:收稿日期:2023-04-19;修回日期:修回日期:2023-07-30 基金项目:基金项目:国家自然科学基金(11503031);陕西省自然科学基础研究计划(2023-JC-YB-057,2022-JM031)作者简介:作者简介:雷雨(1983),男,讲师,博士,从事时间保持方法研究。文章编号:文章编号:1005-6734(2023)09-0909-09 doi.10.13695/ki.12-1222/o
2、3.2023.09.008 基于奇异谱分析的基于奇异谱分析的 BDS 卫星钟差周期项提取卫星钟差周期项提取 雷 雨1,赵丹宁2(1.西安邮电大学 计算机学院,西安 710121;2.宝鸡文理学院 电子电气工程学院,宝鸡 721016)摘要:摘要:在建立星座自主时间基准时,必须扣除卫星钟差的周期波动,以避免将其引入系统时间。为准确地扣除周期波动,提出一种基于奇异谱分析(SSA)的 BDS 卫星钟差周期项提取方法。首先,经SSA 分解获得钟差信号的多个重构成分;然后,引入重标极差分析方法计算各重构成分的 Hurst 指数,根据 Hurst 指数辨识钟差的低频信号主导分量和高频信号随机分量,从而实现
3、周期项的重建提取。对不同类型的 BDS 卫星钟进行了分析和研究,结果表明:所提方法能够准确提取钟差周期项,所提取周期项的频谱特征比钟差多项式拟合残差的频谱更为清晰;利用多项式拟合方法扣除周期项后,万秒稳和日稳平均分别提升 4.1%和 2.4%,而利用所提方法扣除周期项后,万秒稳和日稳平均分别提升 20.9%和 34.1%,提高了卫星钟频率稳定度,为星座自主时间基准的建立提供了基础。关 键 词:关 键 词:北斗卫星导航系统;卫星钟;周期项;奇异谱分析;频率稳定度 中图分类号:中图分类号:P228.1 文献标志码:文献标志码:A Periodic term extraction of BDS sa
4、tellite clock offset based on singular spectrum analysis LEI Yu1,ZHAO Danning2(1.School of Computer Science&Technology,Xian University of Posts&Telecommunications,Xian 710121,China;2.School of Electrical&Electronic Engineering,Baoji University of Arts and Sciences,Baoji 721016,China)Abstract:In esta
5、blishing the autonomous time reference of the navigation constellation,the periodic fluctuations of satellite clock offset must be deducted to avoid introducing them into the system time.In order to accurately remove the periodic fluctuations,a decomposition,identification and rebuilt method based o
6、n singular spectrum analysis(SSA)is proposed to extract the periodic terms of BeiDou satellite clock offset.Firstly,multiple reconstructed components of the clock offset signal are obtained by SSA decomposition.Then,rescaled range analysis is introduced to calculate the Hurst index of each reconstru
7、cted component.The dominant low-frequency components and random high-frequency components of clock offset subsequently are identified according to the Hurst index,so as to realize the reconstruction and extraction of the period term.The analysis and study of different types of BDS satellite clocks s
8、how that the proposed method can accurately extract the period terms of satellite clock offset,and the spectral features of the extracted periodic terms are clearer than that of the polynomial fitting residuals.While the mean improvements are 4.1%and 2.4%respectively for the frequency stability at t
9、he ten-thousand and daily sampling intervals via the polynomial fitting method,through the proposed method the improvements can reach up to 20.9%and 34.1%,respectively,which improves the frequency stability of the satellite clock and provides the reference for the establishment of an autonomous time
10、 reference for the navigation constellation.Key words:Beidou navigation satellite system;satellite clock;periodic terms;singular spectrum analysis;frequency stability 910 中国惯性技术学报 第 31 卷 卫星导航系统的定位、导航和授时服务依赖于准确的时间基准,目前北斗卫星导航系统(Beidou Navigation Satellite System,BDS)时间基准的建立和保持由地面监测站原子钟通过综合原子时算法计算产生,这种
11、模式过度依赖地面站,一旦地面监测站出现异常,系统时间基准就会不连续。为保证系统时间基准的连续性和稳定性,必须降低对地面监测站的依赖。随着 BDS 卫星数量与星载钟性能的不断提高,利用星载钟建立和保持星座自主时间基准,对提高系统的生存能力具有重要的现实意义1,2。受星载钟自身因素和太空环境的影响,卫星钟差序列一般含有趋势项、周期项和随机项3,其中,星载钟受太空环境影响的主要表现为钟差的周期波动,导致星载钟在轨性能逊于地面钟。卫星钟差的周期波动不仅和其轨道周期相关,而且和日食周期、测量性能等有关4,5。在未扣除卫星钟差的周期波动之前,不应将星载钟纳入星座自主时间基准的建立和保持,以避免将周期波动引
12、入系统时间。现有的卫星钟差特征分量提取方法有多项式拟合、傅里叶变换、小波变换和经验模态分解等,Huang 等构建多项式模型拟合提取钟差周期项,提升了 BDS 卫星钟差预报精度3;李骁逸等基于傅里叶变换的频谱分析方法校正 BDS 卫星钟的周期波动,校正后 BDS 各类星载钟的万秒稳均得到提高6;雷雨等利用小波变换将钟差序列分解成具有不同频率特征的信号分量,并根据各分量的特征构建预报模型,改进了钟差预报精度7;梁益丰等提出以完备集合经验模态分解为基础的钟差信噪分离方法,并将其应用于 BDS 卫星钟差周期项识别,扣除周期项后卫星钟万秒稳获得提升8。上述方法虽然在钟差特征分量分离和提取方面取得了较好的
13、效果,但均存在一定的局限性。传统的多项式拟合方法无法提取钟差序列中的时变信号,提取的钟差趋势项、周期项和随机项不够准确;傅立叶变换将信号从时间域转化到频率域,以分离信号中的周期分量,但该方法不适用于非线性、非平稳信号;经验模态分解能将非线性、非平稳信号分解为一系列固有模态分量,并根据不同模态分量的特征识别周期成分,但这类方法存在端部效应和模态混叠的问题。奇异谱分析(Singular Spectrum Analysis,SSA)是20 世纪 90 年代兴起的一种研究非线性、非平稳信号的有效方法。根据所观测到的时间序列构造出轨迹矩阵,并对轨迹矩阵进行分解和重构,从而提取出代表原时间序列不同成分的信
14、号,如长期趋势信号、周期信号和噪声信号,从而对时间序列的波动特征进行分析9,10。SSA 方法提取信号不需要先验信息,不受正弦波假定的约束,能够较好地从含噪声的时间序列中提取时变信号,目前已在 GPS 坐标时间序列分析、多路径效应消除、钟差分析和预报11,12等方面取得了成功应用。其中,肖胜红等11提出 SSA 和傅立叶带通滤波器相结合的卫星钟差周期项提取方法,有效抑制了重建周期项时频率混叠和边界效应;Xue 等12利用SSA 分离卫星钟差的趋势项和随机项,提高了钟差短期预报精度。上述研究表明,SSA 可以较好地从复杂卫星钟差序列中提取趋势项和周期项信息。本文在分析卫星钟差特性的基础上,提出一
15、种基于 SSA 的“分解辨识重建”钟差周期项提取方法。首先应用 SSA 对钟差序列进行分解,获得时间序列的多个重构成分,然后利用重标极差分析(Rescaled Range Analysis,R/S)方法计算各重构成分的 Hurst 指数,最后根据 Hurst 指数对各重构成分的波动特征进行辨识,进而完成钟差趋势项、周期项和随机项的分离和提取。利用 BDS 卫星钟差数据分析和检验了本文方法的有效性与实用性。1 星载原子钟物理特性 1 星载原子钟物理特性 1.1 星载原子钟物理模型 受星载原子钟自身因素及空间外界环境的影响,星载原子钟信号不仅存在趋势性变化,而且表现出复杂的周期性波动。星载原子钟模
16、型通常可表示为13:20001122301()sin(2)2()()d()tx txy tdtAf tW tW sst (1)其中,()x t为t时刻原子钟相位;0 x为原子钟初始相位;0y为原子钟钟速;d为原子钟频漂;A、0f和分别为原子钟相位周期分量的振幅、频率和初始相位;11()W t和220()dtW ss表示两种起主导作用的原子钟 噪声,1()W t和2()W s分别为两个独立的维纳过程,1和2分别为两个维纳过程的扩散系数,用来表示两种原子钟噪声强度;3()t 为原子钟测量噪声,3表示测量噪声强度。与趋势分量相比,周期分量和随机分量的数量级较小,但对频率稳定度的影响会逐渐累积。1.2
17、 各分量对频率稳定度的影响 阿伦方差与原子钟频漂、相位白噪声、维纳过程的扩散系数、周期分量之间的关系可表示为13:2422222230122203sin()11()23()yfdAf (2)其中,2()y表示平滑时间为时的阿伦方差;21/、22 和2233/分别表示频率白噪声、频率随机游走噪声和相位白噪声对频率稳定度的影响。根据式(2)可知,当星载原子钟在轨运行时钟差周期特征显著时,第9期 雷雨等:基于奇异谱分析的BDS卫星钟差周期项提取 911 频率稳定度也将出现周期性波动,表现为时域频率稳定度的异常突变,从而影响原子钟噪声系数的拟合估计。因此,在建立星座自主时间基准时,必须消除卫星钟差的的
18、周期性效应,以免将其引入系统时间。2 基于2 基于 SSA 的卫星钟差周期项分析方法的卫星钟差周期项分析方法 基于SSA的卫星钟差周期项分析方法主要包括钟差数据预处理、钟差分解与重构以及钟差重构成分辨识等过程,图1给出卫星钟差周期项分析流程。图 1 卫星钟差周期项分析流程 Fig.1 Analysis process of satellite clock offset periodic terms 已知卫星钟差序列为12,Nx xx,钟差周期项分析的具体步骤为:1)钟差数据预处理 运用中位数法检测和剔除钟差序列中的粗差,若钟差的一次差分数据ix(1,21iN)满足式(3)即认为ix为粗差14:
19、MAD=medianMAD=median/0.6745iiixmnmxxm (3)其中,m表示钟差一次差分序列的中位数;n表示检测阈值,为保证粗差检测的准确性和可靠性,本文取=6n。对于粗差点,利用线性内插法对其进行插补。2)钟差序列分解 利用SSA对数据预处理后的钟差序列进行分解,主要计算过程为9,10:奇异值分解。根据钟差序列构建时滞矩阵121,NMXY YY,矩阵X的每行、每列为钟差序列的子序列,即T121,ii Mx xxY,11i NM,其中M表示时滞窗口长度。对时滞矩阵的协方差矩阵TWXX进行奇异值分解,获得特征值120M及矩阵W对应的特征向量是12M、UUU,其中jU称为时间正交
20、函数。令T(1,2)jjjjMVX U,则时滞矩阵的奇异值分解可写为:12+IXXX+X (4)其中,max|0jIj;初等矩阵TjjjjXU V。对角平均。计算时滞矩阵X在jU上的投影:T,i jijaX U (5)其中,,i ja称为时间主成分。根据时间正交函数与时间主成分进行对角平均,得到钟差序列的重构成分(Reconstructed Component,RC):*,1*,*11*,11 11 1 1ii k jj kkMi ji k jj kkNKi k jj kk i KiMizMiKMKi NNi aUaUaU (6)其中,*min(,1)MM NM,*max(,1)KM NM。分
21、组。卫星钟差序列中通常包含多种周期性变化、趋势性变化和随机性变化,利用w-correlation分析RC成分之间的相关性,将信号特征相似的RC分组。将不同RC分别用()iZ、()jZ表示,则重构分量之间的w-correlation可表示为15:()()()(),()()()()(),1,ijiji jwwiijwNijkk ik jkw z zZZZZZZZZZ 1,i j M(7)其 中,kw为 权 重 系 数,其 定 义 为min(,)kwk M Nk;,i j的绝对值越接近于1,说明()iZ、()jZ两者之间的相关性越强。根据经验,,0.6i j即认为两者之间存在相关性,因此,将RC成分
22、之间相关性大于0.6的两者视为是同一组信号。3)RC成分辨识 一个时间序列在一段时间内的波动如何随时间跨度大小而变化往往可以揭示该时间序列的特性,对于卫星钟差的趋势项和周期项,其时间序列当前或过去的取值以远超随机扰动所能达到的程度影响该序列在未来的取值,统计学上称为时间序列存在长期记忆性(Long-term Memory)和长时间相关效应,而随机项的未来值和当前或过去的取值相关性不强,长期记忆性较弱,属于均值回复过程。由英国水利学家赫尔斯特提出Hurst指数是用来衡量时间序列是否有长期记忆性的一个指标16,17,体现了时间序列的自相关性,尤其能够反映时间序列中隐藏的长期趋势。本文采用Hurst
23、指数H来辨识RC成分序列特征,以定性分析钟差序列中的信号主导分量和噪声主导分量。若RC成分的Hurst指数0.51H,则RC成分属于趋势分量或周期信号主导分量,且H越大,规律性越强;若00.5H,则判定RC成分为随机噪声主导分量,912 中国惯性技术学报 第 31 卷 且H越小,随机性越强。基于R/S分析法计算Hurst指数的原理为16,17:对于时间序列12,Nx xx,将其分为C个长度为L的等长子区间,/LN C,则累积离差为:,1()kc kc icidxe,1,2kL (8)其中,ce为第c个子区间序列的平均值,1,2cC。令第c个子区间序列的标准差为cS,则极差cR为:,max()m
24、in()cc kc kRdd (9)定义重标极差L为:11CHcLccROLCS (10)其中,H为Hurst指数。对L和L进行双对数线性回归拟合,回归方程的截距就是式(11)中的常数O,而斜率就是H。3 算例分析 3 算例分析 采 用 德 国 地 球 科 学 研 究 中 心(Deutsches Geoforschungs Zentrum,GFZ)提供的BDS精密卫星钟差数据(ftp:/ftp.gfz-potsdam.de/GNSS/products/mgex),选取钟差序列连续性较好的C05、C09、C11、C43星载钟,包括GEO(C05)、IGSO(C09)与MEO(C11与C43)轨道
25、类型,铷原子钟(C05、C09、C11)与氢原子钟(C43)类型,选择2022年7月29日至8月11日共14天4032个数据点,数据采样间隔为5 min。利用SSA对时间序列进行分析,重点是确定窗口长度M,一般1/2MN且为周期的最小公倍数。BDS卫星钟差序列中通常存在6 h、12 h和24 h的周期18-20,因此,选择窗口长度M为卫星钟差周期的最小公倍数,即288M。3.1 钟差序列分解和重建试验 由于篇幅限制,以C05卫星钟差为例展示钟差序列分解和重建的过程及效果。卫星钟差数据经SSA分解获得288个RC序列,频率按照RC阶次依次从低到高排列。原始钟差序列与所有RC之和的残差均方根为1.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 奇异 谱分析 BDS 卫星 周期 提取
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。