CN113507279B - 一种高精度星间距离观测数据的降速率滤波方法 - Google Patents

一种高精度星间距离观测数据的降速率滤波方法 Download PDF

Info

Publication number
CN113507279B
CN113507279B CN202110653549.0A CN202110653549A CN113507279B CN 113507279 B CN113507279 B CN 113507279B CN 202110653549 A CN202110653549 A CN 202110653549A CN 113507279 B CN113507279 B CN 113507279B
Authority
CN
China
Prior art keywords
measurement value
filter
satellite
calculating
inter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110653549.0A
Other languages
English (en)
Other versions
CN113507279A (zh
Inventor
刘玄
钟兴旺
杨珊珊
王登峰
薛大雷
张攀
任帅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN202110653549.0A priority Critical patent/CN113507279B/zh
Publication of CN113507279A publication Critical patent/CN113507279A/zh
Application granted granted Critical
Publication of CN113507279B publication Critical patent/CN113507279B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/06Non-recursive filters
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H2017/0072Theoretical filter design
    • H03H2017/0081Theoretical filter design of FIR filters
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明涉及一种高精度星间距离观测数据的降速率滤波方法,属于卫星星间测量技术领域;步骤一、设置采样率为10Hz的原始距离测量值输入数据为步骤二、计算归一化因子FNorm;步骤三、计算用于产生距离测量值的滤波器系数Fn;步骤四、计算用于产生速度测量值的滤波器系数步骤五、计算用于产生加速度测量值的滤波器系数步骤六、获得数字低通滤波器的表达式;本发明具备降速和差分运算两种功能,具有良好的频率截止特性,最大程度保留了重力场信息的同时减少了高频噪声的干扰,有利于星间距离数据在高精度重力场反演计算中的应用。

Description

一种高精度星间距离观测数据的降速率滤波方法
技术领域
本发明属于卫星星间测量技术领域,涉及一种高精度星间距离观测数据的降速率滤波方法。
背景技术
在以重力测量卫星微米级精度测距系统(KBR)、导航卫星高精度信号处理以及其他科学测量领域,观测数据的高精度滤波处理十分重要,在保留有效观测信息的同时也要降低高频噪声干扰,需要针对性的算法来实现。
本项目来源于我国低低跟踪重力测量卫星的工程研制。地球重力场信息分布在0.1Hz~0.0001Hz之间,属于极低频的信息。有效载荷输出的KBR测量数据速率为10Hz,为了降低高频噪声、满足重力场测量频带需求,需要对10Hz数据进行降速率滤波,同时实现一次和二次变率的提取。
针对上述应用背景,现有典型的方法可分为两类:
(1)基于不同窗函数的FIR滤波器,可实现典型的数字低通滤波功能,根据不同的阶数,可以实现一定的幅频截止特性,适应性比较广泛。比如专利CN201710858438.7《低通滤波器的滤波方法、低通滤波器及伺服驱动器》公开了一种低通滤波器的滤波方法、低通滤波器及伺服驱动器。本发明中,一种低通滤波器的滤波方法,包括:周期性采样数据,得到低通滤波器的滤波输入值;根据采样到的数据进行滤波计算,输出滤波后的数据。
(2)基于CIC抽取的滤波器,一般应用在信号处理、降速率运算,常见于基于FPGA的软件无线电设备,比如导航接收机。该方法结构简单,设计灵活,适应性比较广泛。典型的比如北京宇航系统工程研究所田园园等发表的论文《积分清零与积分梳状滤波器在抽取滤波中的应用研究》主要论述了选择合适的抽取滤波器可以使得效率和资源达到最佳的平衡。积分清零和积分梳状滤波器是两种实现简单、滤波性能较好的数字抽取滤波器.这两种滤波器在中频数字接收机中都有重要的应用。
分析现有方法可以总结出如下几个不足之处:
(1)没有一种能够同时融合降速率采样和差分运算两种功能的算法结构,要么仅涉及低通滤波,要么仅涉及降速率抽取算法;
(2)相关方法在应用于极低频信号处理中,在有限的资源和运算量下,很难获得良好的幅频截止特性,抗干扰和抗混叠效能略差。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提出一种高精度星间距离观测数据的降速率滤波方法,具备降速和差分运算两种功能,具有良好的频率截止特性,最大程度保留了重力场信息的同时减少了高频噪声的干扰,有利于星间距离数据在高精度重力场反演计算中的应用。
本发明解决技术的方案是:
一种高精度星间距离观测数据的降速率滤波方法,包括如下步骤:
步骤一、设置采样率为10Hz的原始距离测量值输入数据为
步骤二、计算归一化因子FNorm
步骤三、计算用于产生距离测量值的滤波器系数Fn
步骤四、计算用于产生速度测量值的滤波器系数
步骤五、计算用于产生加速度测量值的滤波器系数
步骤六、获得数字低通滤波器的表达式,通过数字低通滤波器实现对星间距离观测数据的降速率滤波。
在上述的一种高精度星间距离观测数据的降速率滤波方法,所述步骤二中,归一化因子FNorm的计算方法为:
式中,Nh为积分点数;
f0为重力场2次谐波分量;
fs为原始数据采样率;
Hk为非归一化频率响应;
k为采样点,k取值范围为(-Nh,+Nh);
n为频率因子;
Nf为滤波器时域点数,且Nf为奇整数。
在上述的一种高精度星间距离观测数据的降速率滤波方法,所述积分点数Nh的计算方法为:
Nh=(Nf-1)/2。
在上述的一种高精度星间距离观测数据的降速率滤波方法,所述非归一化频率响应Hk的计算方法为:
式中,NB为通带内频点数量;
m为循环计数值;
Nc为卷积次数。
在上述的一种高精度星间距离观测数据的降速率滤波方法,所述通带内频点数量NB的计算方法为:
NB=BTf
式中,B为重力卫星数据信号带宽;
Tf为滤波器时域宽度。
在上述的一种高精度星间距离观测数据的降速率滤波方法,所述步骤三中,用于产生距离测量值的滤波器系数Fn的计算方法为:
在上述的一种高精度星间距离观测数据的降速率滤波方法,所述步骤四中,用于产生速度测量值的滤波器系数的计算方法为:
在上述的一种高精度星间距离观测数据的降速率滤波方法,所述步骤五中,用于产生加速度测量值的滤波器系数的计算方法为:
在上述的一种高精度星间距离观测数据的降速率滤波方法,所述步骤六中,数字低通滤波器的表达式为:
式中,为采样率为0.2Hz的距离输出数据;
为采样率为0.2Hz的速度输出数据;
为采样率为0.2Hz的加速度输出数据。
本发明与现有技术相比的有益效果是:
(1)本发明面向重力场测量的高精度测距载荷原始观测量的降速滤波方法,在充分的需求分析基础上,基于FIR滤波器设计了改进的卷积滤波器,具备降速和差分运算两种功能,具有良好的频率截止特性,最大程度保留了重力场信息的同时减少了高频噪声的干扰,有利于星间距离数据在高精度重力场反演计算中的应用;
(2)本发明滤波算法的架构可以适应并行计算,可并行获得降速率的距离、速度和加速度三种测量值。
附图说明
图1为本发明降速滤波流程示意图。
具体实施方式
下面结合实施例对本发明作进一步阐述。
对于重力卫星KBR的双单向有偏距离测量值而言,其标称数据速率是10Hz。本发明提供一种高精度星间距离观测数据的降速率滤波方法,能够降速率提取输出有偏距离、速度以及加速度数据,在降低高频噪声的同时,按需保留地球重力场信息。因此,需要严格设计滤波器在频率域的性能表现,尤其是抗混叠效能。根据地球重力场球谐系数的估计理论,球谐系数的估计方法实质上是从相关直接观测量(比如星间距离观测量)提取各个单音频率幅度值的总和,这些频率数值代表特定的地球重力场谐波信息。因此,所得谐波系数的估计误差将由对应单音频率上面的噪声决定。据此可知,球谐系数的估计过程是有深度的频率关联特性,所以数字滤波算法的性能应该在频率域上着重分析,需仔细考虑所有的噪声混叠进入重力场的信号频带内。根据重力场反演相关理论,对于100x100阶的地球重力场模型系数,KBR观测数据的有效带宽约0.1mHz~18mHz;对于200x200阶模型系数,带宽约0.1mHz~36mHz。由于重力场模型球谐系数是幅度量,因此对于KBR观测量滤波算法的一个重要需求就是在有效信号带宽内尽可能减小幅频失真。
高精度星间距离观测数据的降速率滤波方法,基于经典FIR滤波器,通过优化设计阶数、截止特性、窗函数以及运算结构,提出了一种结合抽取降速率、差分运算的卷积滤波器,称作CEDF滤波器(convolution extract differential filter),如图1所示,具体包括如下步骤:
步骤一、设置采样率为10Hz的原始距离测量值输入数据为
步骤二、计算归一化因子FNorm;归一化因子FNorm的计算方法为:
式中,Nh为积分点数;积分点数Nh的计算方法为:
Nh=(Nf-1)/2。
f0为重力场2次谐波分量;f0=0.37mHz;
fs为原始数据采样率;fs=10Hz;
Hk为非归一化频率响应;非归一化频率响应Hk的计算方法为:
式中,NB为通带内频点数量;通带内频点数量NB的计算方法为:
NB=BTf
式中,B为重力卫星数据信号带宽;B=0.1Hz。
Tf为滤波器时域宽度,滤波器时域宽度Tf=70.7s。
m为循环计数值;
Nc为卷积次数。
k为采样点,k取值范围为(-Nh,+Nh);
n为频率因子;
Nf为滤波器时域点数,且Nf为奇整数,优选707。
步骤三、计算用于产生距离测量值的滤波器系数Fn;用于产生距离测量值的滤波器系数Fn的计算方法为:
步骤四、计算用于产生速度测量值的滤波器系数用于产生速度测量值的滤波器系数/>的计算方法为:
步骤五、计算用于产生加速度测量值的滤波器系数用于产生加速度测量值的滤波器系数/>的计算方法为:
步骤六、获得数字低通滤波器的表达式,数字低通滤波器的表达式为:
式中,为采样率为0.2Hz的距离输出数据;
为采样率为0.2Hz的速度输出数据;
为采样率为0.2Hz的加速度输出数据。。
通过数字低通滤波器实现对星间距离观测数据的降速率滤波。
在频域内,滤波器旁瓣需要能够抑制带外噪声,使噪声不能进入重力场信号带内,因此滤波器频率响应必须有足够小的旁瓣。对于重力场模型,最大的谐波分量为J2项(重力信号的谐波分量,频率近似为0.37mHz),其旁瓣抑制度需要≤-70dB。
对于KBR测距数据处理来说,幅频响应为矩形的滤波器是最理想的。在理论上,只要Sa函数(sinx/s)是无限长,这样可以通过在时域内将一个信号与Sa函数卷积运算来实现所需的滤波处理。而实际情况往往是Sa函数被限定在一个可实现、易于处理的时域长度内,这种截断一般是通过加窗函数来实现的。
Fn可以形象的看成是从矩形窗口看到的一段Fd(n)序列,即Fn=Fd(n)RN(n)。用有限长序列代替无限长序列肯定会引起误差,产生通带和阻带内的波动,为了缓和这种矛盾,CEDF算法采用在时域内逐个乘以矩形窗函数。从时域分析,该方法的优点就是滤波器的主瓣宽度和旁瓣幅度计算简单。另外,利用傅立叶变换的微分性质,在时域内对函数求一阶导数(产生速度测量值)和二阶导数(产生加速度测量值),在频域内就是分别乘以-2πf和-(2πf)2。从频域分析,CEDF滤波器可以按需选择矩形窗函数自身卷积的次数,定义为Nc。由于矩形窗函数在频域内有sinx/x的傅立叶变换形式,则由卷积定理可知
据此,一个时域内自身卷积N次的矩形窗函数的频域响应是(sinx/x)N的形式。随着自身卷积次数Nc的增加,窗函数傅立叶变换的旁瓣幅度迅速降低,当Nc=7时,第一旁瓣和第二旁瓣幅度衰减量相对于主瓣来说分别为(约-94dB)和/>(约-125dB),即可以满足保留重力场有效信息的需求。
本发明解决了基于卫星的地球重力场反演应用中的微米级精度距离测量值降速滤波及并行差分运算的工程技术难题。面向重力场测量的高精度测距载荷原始观测量的降速滤波方法,在充分的需求分析基础上,基于FIR滤波器设计了改进的卷积滤波器,具备降速和差分运算两种功能,具有良好的频率截止特性,最大程度保留了重力场信息的同时减少了高频噪声的干扰,有利于星间距离数据在高精度重力场反演计算中的应用。同时,滤波算法的架构可以适应并行计算,可并行获得降速率的距离、速度和加速度三种测量值。
本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案的保护范围。

Claims (4)

1.一种高精度星间距离观测数据的降速率滤波方法,其特征在于:包括如下步骤:
步骤一、设置采样率为10Hz的原始距离测量值输入数据为
步骤二、计算归一化因子FNorm
所述步骤二中,归一化因子FNorm的计算方法为:
式中,Nh为积分点数;
f0为重力场2次谐波分量;
fs为原始数据采样率;
Hk为非归一化频率响应;
k为采样点,k取值范围为(-Nh,+Nh);
n为频率因子;
Nf为滤波器时域点数,且Nf为奇整数;
步骤三、计算用于产生距离测量值的滤波器系数Fn
所述步骤三中,用于产生距离测量值的滤波器系数Fn的计算方法为:
步骤四、计算用于产生速度测量值的滤波器系数
所述步骤四中,用于产生速度测量值的滤波器系数的计算方法为:
步骤五、计算用于产生加速度测量值的滤波器系数
所述步骤五中,用于产生加速度测量值的滤波器系数的计算方法为:
步骤六、获得数字低通滤波器的表达式,通过数字低通滤波器实现对星间距离观测数据的降速率滤波;
所述步骤六中,数字低通滤波器的表达式为:
式中,为采样率为0.2Hz的距离输出数据;
为采样率为0.2Hz的速度输出数据;
为采样率为0.2Hz的加速度输出数据。
2.根据权利要求1所述的一种高精度星间距离观测数据的降速率滤波方法,其特征在于:所述积分点数Nh的计算方法为:
Nh=(Nf-1)/2。
3.根据权利要求2所述的一种高精度星间距离观测数据的降速率滤波方法,其特征在于:所述非归一化频率响应Hk的计算方法为:
式中,NB为通带内频点数量;
m为循环计数值;
Nc为卷积次数。
4.根据权利要求3所述的一种高精度星间距离观测数据的降速率滤波方法,其特征在于:所述通带内频点数量NB的计算方法为:
NB=BTf
式中,B为重力卫星数据信号带宽;
Tf为滤波器时域宽度。
CN202110653549.0A 2021-06-11 2021-06-11 一种高精度星间距离观测数据的降速率滤波方法 Active CN113507279B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110653549.0A CN113507279B (zh) 2021-06-11 2021-06-11 一种高精度星间距离观测数据的降速率滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110653549.0A CN113507279B (zh) 2021-06-11 2021-06-11 一种高精度星间距离观测数据的降速率滤波方法

Publications (2)

Publication Number Publication Date
CN113507279A CN113507279A (zh) 2021-10-15
CN113507279B true CN113507279B (zh) 2024-05-03

Family

ID=78009901

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110653549.0A Active CN113507279B (zh) 2021-06-11 2021-06-11 一种高精度星间距离观测数据的降速率滤波方法

Country Status (1)

Country Link
CN (1) CN113507279B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5537368A (en) * 1993-09-22 1996-07-16 The United States Of America As Represented By The Secretary Of The Navy Enhanced adaptive statistical filter providing improved performance for target motion analysis noise discrimination
EP0930725A2 (de) * 1998-01-16 1999-07-21 Alcatel Einrichtung und Verfahren zur Echounterdrückung mit adaptiven FIR-Filtern
DE10138720C1 (de) * 2001-08-07 2003-04-24 Siemens Ag Verfahren zur Bestimmung der optimalen Filterkoeffizienten eines adaptiven Filters durch Lösung der Wiener-Hopfgleichung
CN101211177A (zh) * 2006-12-29 2008-07-02 中国科学院沈阳计算技术研究所有限公司 一种基于滤波技术的数控系统加减速控制方法
CN101277103A (zh) * 2007-03-31 2008-10-01 索尼德国有限责任公司 自适应滤波器装置和用于确定滤波器系数的方法
CN102735207A (zh) * 2012-06-18 2012-10-17 西安空间无线电技术研究所 一种基于距离和载波相差的二维角度解算方法
KR20180094571A (ko) * 2017-02-16 2018-08-24 한국전자통신연구원 주파수 영역 nlms 알고리즘을 사용하는 적응 필터의 수렴 인자 최적화 방법 및 그 장치
KR20210030864A (ko) * 2019-09-10 2021-03-18 고려대학교 산학협력단 주파수 필터 회로

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7050491B2 (en) * 2001-10-15 2006-05-23 Mcdonald James Douglas Adaptive equalization of digital modulating signal recovered from amplitude-modulated signal subject to multipath
KR101225556B1 (ko) * 2008-08-25 2013-01-23 돌비 레버러토리즈 라이쎈싱 코오포레이션 프리-화이트닝을 사용한 lms 알고리즘에 의해 적응된 적응형 필터의 갱신된 필터 계수를 결정하기 위한 방법
US9391591B2 (en) * 2014-02-19 2016-07-12 Pentomics, Inc. Optimal factoring of FIR filters

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5537368A (en) * 1993-09-22 1996-07-16 The United States Of America As Represented By The Secretary Of The Navy Enhanced adaptive statistical filter providing improved performance for target motion analysis noise discrimination
EP0930725A2 (de) * 1998-01-16 1999-07-21 Alcatel Einrichtung und Verfahren zur Echounterdrückung mit adaptiven FIR-Filtern
DE10138720C1 (de) * 2001-08-07 2003-04-24 Siemens Ag Verfahren zur Bestimmung der optimalen Filterkoeffizienten eines adaptiven Filters durch Lösung der Wiener-Hopfgleichung
CN101211177A (zh) * 2006-12-29 2008-07-02 中国科学院沈阳计算技术研究所有限公司 一种基于滤波技术的数控系统加减速控制方法
CN101277103A (zh) * 2007-03-31 2008-10-01 索尼德国有限责任公司 自适应滤波器装置和用于确定滤波器系数的方法
CN102735207A (zh) * 2012-06-18 2012-10-17 西安空间无线电技术研究所 一种基于距离和载波相差的二维角度解算方法
KR20180094571A (ko) * 2017-02-16 2018-08-24 한국전자통신연구원 주파수 영역 nlms 알고리즘을 사용하는 적응 필터의 수렴 인자 최적화 방법 및 그 장치
KR20210030864A (ko) * 2019-09-10 2021-03-18 고려대학교 산학협력단 주파수 필터 회로

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
多普勒积分时长对星间测速精度影响分析;尚湘安 等;《航天器工程》;第29卷(第4期);34-39 *

Also Published As

Publication number Publication date
CN113507279A (zh) 2021-10-15

Similar Documents

Publication Publication Date Title
CN101304404B (zh) 一种宽带时域数字波束形成方法
CN110210101B (zh) 一种基于cem frm的动态非均匀窄过渡带滤波器组及设计方法
CN112350691A (zh) 一种基于多级滤波抽取和两级wola结构的地磁弱信号检测算法
Lim et al. On the synthesis of very sharp decimators and interpolators using the frequency-response masking technique
CN107294512B (zh) 一种基于树型结构的非均匀滤波器组滤波方法
CN113507279B (zh) 一种高精度星间距离观测数据的降速率滤波方法
US7061992B2 (en) Parallel correlator architecture
CN107748354A (zh) 基于分析与综合的宽带数字波束形成装置
CN106603036A (zh) 一种基于低阶内插滤波器的自适应时延估计方法
CN102129073A (zh) 一种合成孔径雷达方位预处理的两级滤波的方法
CN104410388A (zh) 基于非线性变换的自适应时间延迟估计方法
CN114252899B (zh) 一种核信号的级联冲激卷积成形方法和装置
Khademi et al. Reducing the computational complexity of narrowband 2D fan filters using shaped 2D window functions
CN111010144B (zh) 改进的两通道iir的qmfb设计方法
CN110703198B (zh) 基于频率挑选的四元十字阵包络谱估计方法
Zhang et al. The design of digital down converter based on FPGA
CN110932697B (zh) 基于全通滤波器的两通道iir的qmf组设计方法
Li et al. Compensation method for the CIC filter in digital down converter
JP3984131B2 (ja) 波形復元装置及び波形復元方法
CN106338743B (zh) 基于ccd解算的空频联合自适应调零算法
Gautam et al. Application of linear optimization technique to design a fixed coefficient multiplierless sample rate conversion filter for multiple-of-three decimation factor
CN104360354A (zh) 空频抗干扰处理方法和装置
CN108259024A (zh) 一种过d采样多相dft滤波器组设计方法
CN104880721A (zh) 一种卫星导航抗干扰天线数字信号正交分解方法
CN116545412A (zh) 一种基于frm的可配置完美重构滤波器组低复杂度实现方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant