CN109508510A - 一种基于改进的卡尔曼滤波的铷原子钟参数估计算法 - Google Patents
一种基于改进的卡尔曼滤波的铷原子钟参数估计算法 Download PDFInfo
- Publication number
- CN109508510A CN109508510A CN201811561770.8A CN201811561770A CN109508510A CN 109508510 A CN109508510 A CN 109508510A CN 201811561770 A CN201811561770 A CN 201811561770A CN 109508510 A CN109508510 A CN 109508510A
- Authority
- CN
- China
- Prior art keywords
- parameter
- atomic clock
- value
- parameter estimation
- matrix
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于改进的卡尔曼滤波的铷原子钟参数估计算法,本发明包括参与计算的数学模型的建立,并运用建立的数学模型按照一定的步骤与流程进行计算。本发明对调整参数a进一步处理,通过加入惩罚措施和幅值限制,解决参数预测可能出现较大跳变的问题。本发明具有以下有益效果:本发明通过对调整参数加入惩罚措施和幅值限制,从而在卡尔曼滤波算法的基础上提高了钟差预测精度,一定程度上抑制了滤波发散的问题。
Description
技术领域
本发明属时间频率、信号处理技术领域,特别是涉及一种基于改进的卡尔曼滤波的铷原子钟参数估计算法。
背景技术
目前,时间频率标准以铯钟和氢钟为主,由于他们昂贵的价格和对环境的高要求,大多数用户无法接受。但在实际应用中,很多地方都需要高精度的时频标准。如在电力传输网要求时刻准确度小于±1us,频率准确度优于±1×10-11。因此价格相对较低的铷原子钟便成了首选。
但是,铷原子存在频率偏移和老化的问题,导致输出结果逐渐偏离标准时间,不能满足时间频率精度的要求。因此,现有技术中通过卡尔曼滤波算法对铷钟进行控制,分析GPS接收机输出的信号的噪声的特性,提出适应控制铷原子钟的卡尔曼滤波器初始化参数。但现有的卡尔曼滤波算法中的预测调整参数a可能出现较大的误差,最终导致卡尔曼滤波发散,造成预测可能出现较大跳变的问题。
发明内容
本发明通过对调整参数加入惩罚措施和幅值限制,从而在卡尔曼滤波算法的基础上提高了钟差预测精度,一定程度上抑制了现有技术中的卡尔曼滤波发散的问题。
本发明的基于改进的卡尔曼滤波的铷原子钟参数估计算法采用以下技术方案进行:
一种基于改进的卡尔曼滤波的铷原子钟参数估计算法,所述铷原子钟参数估计方法包括以下模型:
S1,钟差特性
钟差数据的噪声特性会随数据采样间隔δ的不同而发生变化;原子钟的噪声有白色相位噪声、闪变相位噪声、白色频率噪声、闪变频率噪声、随机游走频率噪声5种;铷原子钟与标准参考源的差表现为如式(1)所示规律:
其中ξ(t)是噪声,受各种不确定因素的影响;a是初始时刻铷原子钟与参考源的钟差;b是初始时刻原子钟的频率与参考频率的偏差,即频率偏移;c是老化系数,它是由于铷原子钟本身参数老化引起的,在频率上表现为线性的偏离参考频率,在时间上呈二次方的趋势偏离参考时间源;
S2,传统卡尔曼滤波算法
传统卡尔曼滤波算法通过上一时刻的估计值和当前时刻的测量值对当前时刻的调整参数a,b,c进行估计;根据估计值得到铷原子钟的频率调整量,最终实现对原子钟系统的驾驭;
调整参数a,b,c构成了状态向量xi=[a(i) b(i) c(i)]T,系统从第i-1次采样到第i次采样的递推公式为:
xi=Axi-1+qi-1
(2)
其中,A表示状态转移矩阵,qi-1为输入到系统的三维噪声,零均值正态分布且时间上互不相关;
测量向量yi的矩阵方程表示如下:
yi=Hxi+ri
(3)
其中H为量测矩阵;ri为一维白噪声;
在这里定义R为测量噪声的方差阵;Q为状态噪声的方差阵;
在i次测量后,状态向量的误差方差阵为:
其中为误差协方差阵的先验估计,Pi-1为前一时刻最优估计;
状态向量的先验估计值为:
其中为前一时刻的最优估计;
新状态参量的估计依赖于前一时刻的估计值和当前的测量值:
其中yi为测量值,为估计值,ki为增益矩阵,如下式
最后,更新误差协方差阵
这里,I为单位矩阵;
S3,改进的卡尔曼滤波
由于预测调整参数a可能出现较大的误差,最终导致卡尔曼滤波发散,因此加入惩罚措施用以提高预测精度;
S3.1,调整参数预测误差的计算
预测误差:
其中调整参数a的预测值,yi为实测值;
S3.2,对预测值的修正
其中,为最优的参数估计,μ为惩罚因子;
通过引入惩罚因子,将上一步的预测误差对当前预测带来的影响考虑在内;即用上一步调整参数a的预测值与实测值yi的差δ(i)来修正调整参数a;
S3.3,预测值变化幅度的判断及修正
当则最优参数估计可调整为:
否则,保持不变;其中λ为限制幅度;
通过加入幅值限制措施,解决参数预测可能出现较大跳变的问题,从而可以提高预测精度。
优选的,所述铷原子钟参数估计算法包括以下步骤:
步骤一,系统第一次测量i=1时,状态向量初始值为x0,误差方差阵初始值为P0;
步骤二,计算状态向量的先验估计值为
步骤三,计算调整参数a的预测值为
步骤四,对系统第i次测量进行判断;
(1)若i>1,则最优的参数估计
(2)若i≤1,计算预测误差:随后计算最优的参数估计
对值进行判断;
若则最优参数估计调整为
若则优的参数估计按计算;
步骤五,计算误差方差阵
步骤六,计算增益矩阵
步骤七,计算新状态参量的估计
步骤八,计算更新误差协方差阵
步骤九,根据上述步骤计算,最终得到状态向量xi与误差方差阵Pi。
优选的,所述步骤八中,若计算i=i+1,则继续重复步骤二到步骤八。
优选的,惩罚因子μ的取值范围在0.5~0.7之间;所述限制幅度λ的取值范围在10~15之间。
本发明的作用原理如下:
现有的卡尔曼滤波算法存在滤波发散的问题,本发明通过对调整参数加入惩罚措施和幅值限制,从而在卡尔曼滤波算法的基础上提高了钟差预测精度,一定程度上抑制了滤波发散的问题。传统卡尔曼模型是对调整参数a,b,c进行估计。而本专利是在卡尔曼算法的基础上,对调整参数a进一步处理,通过加入惩罚措施和幅值限制,解决参数预测可能出现较大跳变的问题。
本发明的改进点,首先在于引入了惩罚因子将上一步的预测误差对当前预测带来的影响考虑在内。即用上一步调整参数a的预测值与实测值yi的差δ(i)来修正调整参数a。
其次,通过加入幅值限制措施,解决参数预测可能出现较大跳变的问题,从而可以提高预测精度。对值进行判断,来决定选择哪个模型进行下一步的运算。
在本发明中引入的μ为惩罚因子:在最优化函数中,对离群点的惩罚因子,也是对离群点的重视程度体现。因为算法预测可能出现误差,所以通过加入惩罚因子,将前一时刻产生的误差引入到当前时刻的目标函数的预测中,以此提高预测精度。
在本发明中引入的λ为限制幅度。通过参考时间源对铷原子钟进行控制,当参考时间源发生切换时,钟差数据可能产生大的跳变,此时传统卡尔曼滤波算法可能产生滤波发散的问题。通过加入限制幅度,一定程度上改善了对于钟差出现跳变时产生滤波发散的问题。
以上两个因子的引入都是基于支持向量机(SVM)的惩罚思想。根据一系列的计算与分析,最终确定在本发明中的惩罚因子μ取值就在在0.5~0.7之间,因为此时平均相对预测误差相对较小且相对比较平稳;本发明中的限制幅度λ取值就在在10~15之间,此时平均相对预测误差达到最小。
本发明具有以下有益效果:
本发明通过对调整参数加入惩罚措施和幅值限制,从而在卡尔曼滤波算法的基础上提高了钟差预测精度,一定程度上抑制了滤波发散的问题。
附图说明
图1为本发明改进的卡尔曼滤波的铷原子钟参数估计算法的流程图。
具体实施方式
下面以具体实施例对本发明作进一步描述,在此发明的示意性实施例以及说明用来解释本发明,但并不作为对本发明的限定。
一种基于改进的卡尔曼滤波的铷原子钟参数估计算法,包括以下模型:
S1,钟差特性
钟差数据的噪声特性会随数据采样间隔δ的不同而发生变化;原子钟的噪声有白色相位噪声、闪变相位噪声、白色频率噪声、闪变频率噪声、随机游走频率噪声5种;
铷原子钟与标准参考源的差表现为如式(1)所示规律:
其中ξ(t)是噪声,受各种不确定因素的影响;a是初始时刻铷原子钟与参考源的钟差;b是初始时刻原子钟的频率与参考频率的偏差,即频率偏移;c是老化系数,它是由于铷原子钟本身参数老化引起的,在频率上表现为线性的偏离参考频率,在时间上呈二次方的趋势偏离参考时间源;
S2,传统卡尔曼滤波算法
传统卡尔曼滤波算法通过上一时刻的估计值和当前时刻的测量值对当前时刻的调整参数a,b,c进行估计;根据估计值得到铷原子钟的频率调整量,最终实现对原子钟系统的驾驭;
调整参数a,b,c构成了状态向量xi=[a(i) b(i) c(i)]T,系统从第i-1次采样到第i次采样的递推公式为:
xi=Axi-1+qi-1
(2)
其中,A表示状态转移矩阵,qi-1为输入到系统的三维噪声,零均值正态分布且时间上互不相关;
测量向量yi的矩阵方程表示如下:
yi=Hxi+ri
(3)
其中H为量测矩阵;ri为一维白噪声;
在这里定义R为测量噪声的方差阵;Q为状态噪声的方差阵;在i次测量后,状态向量的误差方差阵为:
其中为误差协方差阵的先验估计,Pi-1为前一时刻最优估计;状态向量的先验估计值为:
其中为前一时刻的最优估计;
新状态参量的估计依赖于前一时刻的估计值和当前的测量值:
其中yi为测量值,为估计值,ki为增益矩阵,如下式
最后,更新误差协方差阵
这里,I为单位矩阵;
S3,改进的卡尔曼滤波
由于预测调整参数a可能出现较大的误差,最终导致卡尔曼滤波发散,因此加入惩罚措施用以提高预测精度;
S3.1,调整参数预测误差的计算
预测误差:
其中调整参数a的预测值,yi为实测值;
S3.2,对预测值的修正
其中,为最优的参数估计,μ为惩罚因子;
通过引入惩罚因子,将上一步的预测误差对当前预测带来的影响考虑在内;即用上一步调整参数a的预测值与实测值yi的差δ(i)来修正调整参数a;
S3.3,预测值变化幅度的判断及修正
当则最优参数估计可调整为:
否则,保持不变;其中λ为限制幅度;
通过加入幅值限制措施,解决参数预测可能出现较大跳变的问题,从而可以提高预测精度。
运用上述模型,来进行铷原子钟参数估计,铷原子钟参数估计算法包括以下步骤:
步骤一,系统第一次测量i=1时,状态向量初始值为x0,误差方差阵初始值为P0;
步骤二,计算状态向量的先验估计值为
步骤三,计算调整参数a的预测值为
步骤四,对系统第i次测量进行判断;
(1)若i>1,则最优的参数估计
(2)若i≤1,计算预测误差:随后计算最优的参数估计
对值进行判断;
若则最优参数估计调整为
若则优的参数估计按计算;
步骤五,计算误差方差阵
步骤六,计算增益矩阵
步骤七,计算新状态参量的估计
步骤八,计算更新误差协方差阵
步骤九,根据上述步骤计算,最终得到状态向量xi与误差方差阵Pi。
其中,若i值为i=i+1,则继续重复步骤二到步骤八。在本次算法中惩罚因子μ=0.65;限制幅度λ=13。
本发明的技术方案不限于上述具体实施例的限制,凡是根据本发明的技术方案做出的技术变形,均落入本发明的保护范围之内。
Claims (4)
1.一种基于改进的卡尔曼滤波的铷原子钟参数估计算法,其特征在于,所述铷原子钟参数估计方法包括以下模型:
S1,钟差特性
钟差数据的噪声特性会随数据采样间隔δ的不同而发生变化;铷原子钟与标准参考源的差表现为如式(1)所示规律:
其中ξ(t)是噪声,受各种不确定因素的影响;a是初始时刻铷原子钟与参考源的钟差;b是初始时刻原子钟的频率与参考频率的偏差,即频率偏移;c是老化系数;
S2,传统卡尔曼滤波算法
传统卡尔曼滤波算法通过上一时刻的估计值和当前时刻的测量值对当前时刻的调整参数a,b,c进行估计;根据估计值得到铷原子钟的频率调整量,最终实现对原子钟系统的驾驭;
调整参数a,b,c构成了状态向量xi=[a(i)b(i)c(i)]T,系统从第i-1次采样到第i次采样的递推公式为:
xi=Axi-1+qi-1
(2)
其中,A表示状态转移矩阵,qi-1为输入到系统的三维噪声,零均值正态分布且时间上互不相关;
测量向量yi的矩阵方程表示如下:
yi=Hxi+ri
(3)
其中H为量测矩阵;ri为一维白噪声;
在这里定义R为测量噪声的方差阵;Q为状态噪声的方差阵;在i次测量后,状态向量的误差方差阵为:
Pi -=APi-1AT+Q
(4)
其中Pi -为误差协方差阵的先验估计,Pi-1为前一时刻最优估计;
状态向量的先验估计值为:
其中为前一时刻的最优估计;
新状态参量的估计依赖于前一时刻的估计值和当前的测量值:
其中yi为测量值,为估计值,ki为增益矩阵,如下式
ki=Pi -HT(HPi -HT+R)-1
(7)
最后,更新误差协方差阵
Pi=(I-kiH)Pi -
(8)
这里,I为单位矩阵;
S3,改进的卡尔曼滤波
由于预测调整参数a可能出现较大的误差,最终导致卡尔曼滤波发散,因此加入惩罚措施用以提高预测精度;
S3.1,调整参数预测误差的计算
预测误差:
其中调整参数a的预测值,yi为实测值;
S3.2,对预测值的修正
其中,为最优的参数估计,μ为惩罚因子;
通过引入惩罚因子,将上一步的预测误差对当前预测带来的影响考虑在内;即用上一步调整参数a的预测值与实测值yi的差δ(i)来修正调整参数a;
S3.3,预测值变化幅度的判断及修正
当则最优参数估计可调整为:
否则,保持不变;其中λ为限制幅度。
2.根据权利要求1所述的一种基于改进的卡尔曼滤波的铷原子钟参数估计算法,其特征在于:所述铷原子钟参数估计算法包括以下步骤:
步骤一,系统第一次测量i=1时,状态向量初始值为x0,误差方差阵初始值为P0;
步骤二,计算状态向量的先验估计值为
步骤三,计算调整参数a的预测值为
步骤四,对系统第i次测量进行判断;
(1)若i>1,则最优的参数估计
(2)若i≤1,计算预测误差:随后计算最优的参数估计
对值进行判断;
若则最优参数估计调整为
若则优的参数估计按计算;
步骤五,计算误差方差阵Pi -=APi-1AT+Q;
步骤六,计算增益矩阵ki=Pi -HT(HPi -HT+R)-1;
步骤七,计算新状态参量的估计
步骤八,计算更新误差协方差阵Pi=(I-kiH)Pi -;
步骤九,根据上述步骤计算,最终得到状态向量xi与误差方差阵Pi。
3.根据权利要求2所述的一种基于改进的卡尔曼滤波的铷原子钟参数估计算法,其特征在于:所述步骤八中,若计算i=i+1,则继续重复步骤二到步骤八。
4.根据权利要求1-3任一项所述的一种基于改进的卡尔曼滤波的铷原子钟参数估计算法,其特征在于:惩罚因子μ的取值范围在0.5~0.7之间;所述限制幅度λ的取值范围在10~15之间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811561770.8A CN109508510B (zh) | 2018-12-20 | 2018-12-20 | 一种基于改进的卡尔曼滤波的铷原子钟参数估计算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811561770.8A CN109508510B (zh) | 2018-12-20 | 2018-12-20 | 一种基于改进的卡尔曼滤波的铷原子钟参数估计算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109508510A true CN109508510A (zh) | 2019-03-22 |
CN109508510B CN109508510B (zh) | 2022-10-28 |
Family
ID=65753903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811561770.8A Active CN109508510B (zh) | 2018-12-20 | 2018-12-20 | 一种基于改进的卡尔曼滤波的铷原子钟参数估计算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109508510B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110837219A (zh) * | 2019-10-06 | 2020-02-25 | 中国计量科学研究院 | 一种用于监测实体原子钟的虚拟原子钟系统及工作方法 |
CN110850450A (zh) * | 2019-12-03 | 2020-02-28 | 航天恒星科技有限公司 | 一种卫星钟差参数的自适应估计方法 |
CN111641471A (zh) * | 2020-05-28 | 2020-09-08 | 中国计量科学研究院 | 一种原子钟信号组合控制中预测的权重设计策略 |
CN112597622A (zh) * | 2020-10-12 | 2021-04-02 | 北京卫星导航中心 | 一种用于检测铯原子钟频率异常的方法、系统和介质 |
CN112697215A (zh) * | 2020-11-19 | 2021-04-23 | 浙江工业大学 | 一种用于超声波水表数据滤波的卡尔曼滤波参数调试方法 |
CN113176725A (zh) * | 2021-03-05 | 2021-07-27 | 北京大学 | 基于卡尔曼滤波和/或dfb的激光芯片原子钟及实现方法 |
CN113758503A (zh) * | 2021-08-12 | 2021-12-07 | 清华大学 | 一种过程参数估计方法、装置、电子设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6108555A (en) * | 1996-05-17 | 2000-08-22 | Ksi, Inc. | Enchanced time difference localization system |
WO2011085585A1 (zh) * | 2010-07-23 | 2011-07-21 | 华为技术有限公司 | 时间同步的方法和设备 |
CN202102264U (zh) * | 2011-06-24 | 2012-01-04 | 成都可为科技发展有限公司 | 铷钟驯服系统 |
CN107976700A (zh) * | 2017-10-17 | 2018-05-01 | 武汉华中天纬测控有限公司 | 一种卫星导航接收机秒脉冲的稳定输出方法 |
-
2018
- 2018-12-20 CN CN201811561770.8A patent/CN109508510B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6108555A (en) * | 1996-05-17 | 2000-08-22 | Ksi, Inc. | Enchanced time difference localization system |
WO2011085585A1 (zh) * | 2010-07-23 | 2011-07-21 | 华为技术有限公司 | 时间同步的方法和设备 |
CN202102264U (zh) * | 2011-06-24 | 2012-01-04 | 成都可为科技发展有限公司 | 铷钟驯服系统 |
CN107976700A (zh) * | 2017-10-17 | 2018-05-01 | 武汉华中天纬测控有限公司 | 一种卫星导航接收机秒脉冲的稳定输出方法 |
Non-Patent Citations (2)
Title |
---|
T. IAKYMCHUK等: "Versatile Direct and Transpose Matrix Multiplication with Chained Operations: An Optimized Architecture Using Circulant Matrices", 《IEEE TRANSACTIONS ON COMPUTERS》 * |
孟茁等: "卡尔曼滤波在原子钟驯服中的应用", 《现代导航》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110837219A (zh) * | 2019-10-06 | 2020-02-25 | 中国计量科学研究院 | 一种用于监测实体原子钟的虚拟原子钟系统及工作方法 |
CN110837219B (zh) * | 2019-10-06 | 2022-04-05 | 中国计量科学研究院 | 一种用于监测实体原子钟的虚拟原子钟系统及工作方法 |
CN110850450A (zh) * | 2019-12-03 | 2020-02-28 | 航天恒星科技有限公司 | 一种卫星钟差参数的自适应估计方法 |
CN111641471A (zh) * | 2020-05-28 | 2020-09-08 | 中国计量科学研究院 | 一种原子钟信号组合控制中预测的权重设计策略 |
CN111641471B (zh) * | 2020-05-28 | 2022-11-04 | 中国计量科学研究院 | 一种原子钟信号组合控制中预测的权重设计方法 |
CN112597622A (zh) * | 2020-10-12 | 2021-04-02 | 北京卫星导航中心 | 一种用于检测铯原子钟频率异常的方法、系统和介质 |
CN112597622B (zh) * | 2020-10-12 | 2024-01-19 | 北京卫星导航中心 | 一种用于检测铯原子钟频率异常的方法、系统和介质 |
CN112697215A (zh) * | 2020-11-19 | 2021-04-23 | 浙江工业大学 | 一种用于超声波水表数据滤波的卡尔曼滤波参数调试方法 |
CN113176725A (zh) * | 2021-03-05 | 2021-07-27 | 北京大学 | 基于卡尔曼滤波和/或dfb的激光芯片原子钟及实现方法 |
CN113176725B (zh) * | 2021-03-05 | 2023-01-24 | 北京大学 | 基于卡尔曼滤波和/或dfb的激光芯片原子钟及实现方法 |
CN113758503A (zh) * | 2021-08-12 | 2021-12-07 | 清华大学 | 一种过程参数估计方法、装置、电子设备及存储介质 |
CN113758503B (zh) * | 2021-08-12 | 2022-03-18 | 清华大学 | 一种过程参数估计方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109508510B (zh) | 2022-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109508510A (zh) | 一种基于改进的卡尔曼滤波的铷原子钟参数估计算法 | |
CN107765347B (zh) | 一种高斯过程回归和粒子滤波的短期风速预测方法 | |
CN107045125B (zh) | 一种基于预测值量测转换的交互多模型雷达目标跟踪方法 | |
CN105718642B (zh) | 一种基于门限自回归模型的参考时间尺度产生方法 | |
Van Overloop et al. | Identification of resonance waves in open water channels | |
CN110837219B (zh) | 一种用于监测实体原子钟的虚拟原子钟系统及工作方法 | |
US20130204403A1 (en) | Apparatus and Methods for Non-Invasive Closed Loop Step Testing Using a Tunable Trade-Off Factor | |
CN101957598A (zh) | 一种大时滞系统的灰色无模型控制方法 | |
WO2007001252A1 (en) | Apparatuses, systems, and methods utilizing adaptive control | |
Truong et al. | Design of an advanced time delay measurement and a smart adaptive unequal interval grey predictor for real-time nonlinear control systems | |
CN111650617A (zh) | 基于新息加权自适应不敏卡尔曼滤波的晶振频率驯服方法、系统及介质 | |
CN108628288B (zh) | 一种用于一阶时滞组合积分控制系统的性能评价方法 | |
Cheng et al. | Highly accurate real-time GPS carrier phase-disciplined oscillator | |
CN110850450A (zh) | 一种卫星钟差参数的自适应估计方法 | |
Rangegowda et al. | Simultaneous state and parameter estimation using receding-horizon nonlinear kalman filter | |
CN101231508A (zh) | 采用时间序列分析预测-校正集成电路制造工艺控制方法 | |
CN108958027B (zh) | 一种区间系统的动态矩阵控制方法 | |
Goffaux et al. | Bioprocess state estimation: some classical and less classical approaches | |
CN111259340B (zh) | 一种基于logistic回归的饱和负荷预测方法 | |
CN117035001A (zh) | 一种基于Kriging插值的转换力传感器温度补偿方法 | |
Yang et al. | Time-delay system identification using genetic algorithm–part one: Precise fopdt model estimation | |
CN108733622A (zh) | 基于自相关不等间隔迭代法的电离层tec单站预报方法 | |
CN111711446B (zh) | 一种利用gps信号驯服晶振频率的方法、系统及介质 | |
CN112182864B (zh) | 一种基于氢原子钟漂移情况选择钟差预报的方法 | |
van Overloop et al. | Identification of properties of open water channels for controller design |
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 |