CN106126961B - 灰色粒子群卫星钟差自适应预测方法及系统 - Google Patents
灰色粒子群卫星钟差自适应预测方法及系统 Download PDFInfo
- Publication number
- CN106126961B CN106126961B CN201610605339.3A CN201610605339A CN106126961B CN 106126961 B CN106126961 B CN 106126961B CN 201610605339 A CN201610605339 A CN 201610605339A CN 106126961 B CN106126961 B CN 106126961B
- Authority
- CN
- China
- Prior art keywords
- clock correction
- grey
- model
- adaptive
- factor
- 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
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
一种灰色粒子群卫星钟差自适应预测方法及系统,通过在GM(1,1)模型中引入优化因子以及遗忘因子最小二乘法建立灰色粒子群模型,然后采用粒子群优化算法求得最佳优化因子,最终将最佳优化因子代入PGM(1,1)模型并进行钟差预测得到钟差的预测序列,本发明能够提高卫星钟差的预报精度,进而辅助精密单点定位接收机进行高精度实时定位解算,PGM(1,1)预报模型所需训练样本小,模型参数通过自适应方法调节,计算量较小,便于工程运用,且卫星钟差预报精度较高,误差控制在1ns以内。
Description
技术领域
本发明涉及的是一种卫星导航定位领域的技术,具体是一种灰色粒子群卫星钟差自适应预测方法及系统。
背景技术
卫星钟差预测的准确性直接影响卫星导航定位的精度。现有的灰色理论方法,即灰色系统GM(1,1),实质是将无规律的原始数据进行累加,得到规律性较强的生成数列后,再用微分拟合法从新建模,由GM(1,1)模型得到的数据通过累减反向得到预测值。但这种方法模型适应性差,初值选取有缺陷,不适用于近指数序列拟合,背景值计算方式过于简单。
发明内容
本发明针对现有技术难以适用于近指数序列拟合,背景值计算方式过于简单以及预测精度不高等缺陷,提出一种灰色粒子群卫星钟差自适应预测方法及系统,能够显著提高卫星钟差的预报精度,进而辅助精密单点定位接收机进行高精度实时定位解算,PGM(1,1)预报模型所需训练样本小,模型参数通过自适应方法调节,计算量较小,便于工程运用,且卫星钟差预报精度较高,误差控制在1ns以内。
本发明是通过以下技术方案实现的:
本发明涉及一种灰色粒子群卫星钟差自适应预测方法,首先根据钟差的原始序列,通过在GM(1,1)模型中引入优化因子以及遗忘因子最小二乘法建立灰色粒子群模型,即PGM(1,1)模型,然后采用粒子群优化(PSO)算法求得最佳优化因子,最终将最佳优化因子代入PGM(1,1)模型并进行钟差预测得到钟差的预测序列。
所述的优化因子包括:自适应权值ω、发展系数λ和灰作用量θ。
本发明具体包括以下步骤:
1)对钟差的原始序列进行处理得到非负的等维序列x(0)={x(0)(1),x(0)(2),..x(0)(k).,x(0)(n)}其中:k=1,2,...,n;
2)对等维序列进行累加得到一次累加序列(1‐AGO序列)x(1)={x(1)(1),x(1)(2),..x(1)(k).,x(1)(n)},其中:
3)建立包含优化因子和遗忘因子最小二乘法的PGM(1,1)模型;
4)采用PSO算法以历史钟差作为样本求解出最佳优化因子;
5)将最佳优化因子代入PGM(1,1)模型;
6)计算钟差的预测值。
所述的步骤3)具体包括以下步骤:
3.1)通过对1‐AGO序列进行一次紧邻数据自适应寻优生成背景值z(1)={z(1)(2),z(1)(3),...,z(1)(n)},其中:z(1)(k)=ωx(1)(k)+(1-ω)x(1)(k-1),k=2,...,n,0≤ω≤1。
3.2)建立灰微分方程x(0)(k)=az(1)(k)+b后,得到白化方程其中:a,b为模型待辨识参数;
3.3)得到离散化后的初始解其中:k=1,2,...,n,1≤l≤k;
3.4)引入发展系数λ和灰作用量θ后获得预测值
所述的发展系数λ∈[0.5,1.5],灰作用量θ∈(0.9,1.1],且发展系数θb∈[0,1)。
所述的模型待辨识参数[a b]T=(LTWL)-1LTWY,其中:遗忘因子W=diag(e-(n-1),e-(n-2),...,e-(n-n)),Y=[x(0)(2),x(0)(3),...,x(0)(n)]T,
所述的步骤4)具体以下步骤:
4.1)生成初始群体(φ1,φ2,...,φn),其中:φi=[ωi,λi,θi]T;
4.2)以归一化平均相对误差最小为准则,比较每个粒子当前适应值和其个体历史最好适应值,获得全局最好适应值;
4.3)更新各微粒的速度和位置;
4.3)当适应值误差达到设定误差限或到达最大迭代次数则停止搜索并输出最佳优化因子,否则回到步骤4.2)。
所述的归一化平均相对误差
本发明涉及一种实现上述方法的系统,包括:网络数据更新模块、原始数据预处理模块、卫星钟差预报模块以及数据预警模块,其中:网络数据更新模块与原始数据预处理模块相连并传输快速精密星历数据;原始数据预处理模块提取钟差数据截取跳跃点后的钟差数据作为原始钟差序列;卫星钟差预报模块与原始数据预处理模块相连,并接收预处理后原始钟差序列;卫星钟差预报模块与数据预警模块相连,并传输钟差预测值,数据预警模块通过比对钟差预测值和实际钟差值,评估系统预报结果的可用性和准确性。
附图说明
图1为本发明流程示意图;
图2为本发明系统结构示意图;
图3为实施例中G01预测结果及误差示意图;
图4为实施例中G11预测结果及误差示意图;
图5为实施例中G27预测结果及误差示意图;
图6为实施例中G30预测结果及误差示意图;
图7为实施例中G31预测结果及误差示意图。
具体实施方式
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例1
如图1所示,本实施例中,首先获得钟差的原始序列,通过在GM(1,1)模型中引入优化因子以及遗忘因子最小二乘法建立PGM(1,1)模型,然后采用PSO算法求得最佳优化因子,最终将最佳优化因子代入PGM(1,1)模型并进行钟差预测得到钟差的预测序列。
如图2所示,实现本发明方法的系统,包括:网络数据更新模块、原始数据预处理模块、卫星钟差预报模块以及数据预警模块,其中:网络数据更新模块与原始数据预处理模块相连并传输快速精密星历数据;原始数据预处理模块提取钟差数据并处理为非负原始钟差序列,截取跳跃点后的钟差数据作为原始钟差序列;卫星钟差预报模块和数据预警模块分别与原始数据预处理模块相连,并接收预处理后非负原始钟差序列;卫星钟差预报模块与数据预警模块相连,并传输钟差预测值,数据预警模块通过比对钟差预测值和实际钟差值,评估系统预报结果的可用性和准确性。
本实施例中的方法具体包括以下步骤:
1)对钟差的原始序列进行处理得到非负的等维序列x(0)={x(0)(1),x(0)(2),..x(0)(k).,x(0)(n)}其中:k=1,2,...,n。该等维序列作为建模序列,x(0)(1)作为最旧值,x(0)(n)作为最新值。
2)对等维序列进行累加得到一次累加序列x(1)={x(1)(1),x(1)(2),..x(1)(k).,x(1)(n)},其中:
3)建立包含优化因子和遗忘因子最小二乘法的PGM(1,1)模型。
3.1)通过对1‐AGO序列进行一次紧邻数据自适应寻优生成背景值z(1)={z(1)(2),z(1)(3),...,z(1)(n)},其中:z(1)(k)=ωx(1)(k)+(1-ω)x(1)(k-1),k=2,...,n,0≤ω≤1。
3.2)建立灰微分方程x(0)(k)=az(1)(k)+b后,得到白化方程其中:a,b为模型待辨识参数。令A=[a b]T,由包括遗忘因子的最小二乘法得到:其中:W=diag(e-(n-1),e-(n-2),...,e-(n-n)),Y=[x(0)(2),x(0)(3),...,x(0)(n)]T,
3.3)方程初始解为:则得到离散化后的初始解其中:k=1,2,...,n,1≤l≤n,x(1)(l)表示1‐AGO序列中的第l个值。若取l=1,就是直接利用初值即最旧值进行预测。若取l=n则为利用距离预测点最近的建模数据点即最新值进行预测。
3.4)引入发展系数λ和灰作用量θ后获得预测值所述的发展系数λ∈[0.5,1.5],灰作用量θ∈(0.9,1.1],且发展系数θb∈[0,1)。优化因子的组合φ=[ω,λ,θ]T。
4)采用PSO算法以历史钟差作为样本求解出最佳优化因子。
选用归一化平均相对误差作为精度检验标准,以采样点精度为权值对平均相对误差进行归一化平均,归一化平均相对误差其中:q(k)采样点精度权值,若q(k)=1则采样点精度一致。在本实施例中以精密星历中的钟差数据点精度作为精度权值进行归一化。
假设在D维搜索空间中,有m个粒子组成一个微粒群,第i个粒子的飞行速度和位置分别为vi=[vi1,vi2,...,viD]和si=[si1,si2,...,siD],其中:i=1,2,...,m。
根据对环境的适应度将群体中的微粒移动到更好的区域,每个粒子的位置就是一个潜在的解。在每次迭代中,粒子通过跟踪两个“极值”更新自己。第一个极值是粒子本身所找到的最优解,即个体极值pibest=[pi1,pi2,...,piD];另一个极值是整个种群中所找到的最优解,即全局极值Pgbest=[pg1,pg2,...,pg3]。对每一代粒子,在找到这两个最优值时,粒子根据公式来更新自己的速度和位置:其中:w为惯性权值,为改善PSO算法的全局和局部搜索能力,采用惯性权值非线性递减策略,即w=(wstart-wend)(t/tmax)2+(wstart-wend)(2t/tmax)+wstart,wstart和wend为惯性权值上下限,t和tmax为迭代次数和迭代上限。c1、c2为正常数,称为加速系数,通常取值范围在[2,4]之间。r1、r2为在范围[0,1]内变化的随机数。
所述的PSO算法求解最佳优化因子φbest=[ωbest,λbest,θbest]步骤具体包括:
4.1)生成初始群体(φ1,φ2,...,φn),惯性权值上下限wmax和wmin,加速系数c1和c2,最大允许迭代次数tmax以及各个微粒的初始位置和初始速度,其中:φi=[ωi,λi,θi]T;
4.2)以归一化平均相对误差σ最小为各个适应度评判准则,比较每个粒子当前适应值和其个体历史最好适应值,获得全局最好适应值,即:其中:为优化因子组合为φi时计算得到的预测值;
4.3)更新各微粒的速度和位置;
4.3)当适应值误差达到设定误差限或到达最大迭代次数则停止搜索并输出最佳优化因子φbest=[ωbest,λbest,θbest],否则回到步骤4.2)。
5)将最佳优化因子代入PGM(1,1)模型,根据最佳优化因子计算参数a、b,得到预测值的公式。
6)计算钟差的预测值,依据步骤5)得到的预测公式计算预测值。
采用2015年12月1日至2015年12月2日的快速精密星历实测数据进行验证,通过12月1日的钟差数据进行PGM(1,1)模型建模自适应优化,形成最终预测函数后,再对12月2日的钟差数据进行预报,采样间隔与IGS精密星历保持一致为15min,同时与传统的GM(1,1)模型预报精度进行比较。目前在轨运行的4个不同时期发射的GPS卫星所搭载的原子钟由于自身物理结构和运行环境的不同,所表现出的钟差变化也各有其特点。因此选取5颗钟差变化典型的GPS卫星进行预报如表1所示,每颗卫星代表了相应时期的原子钟特性。
表1精密星历数据选取表
如图3~图7所示,并结合表2,PGM(1,1)模型对4颗搭载RB钟卫星钟差的平均预报误差都在0.1ns以内,标准差平均在0.5ns左右,而IGS快速精密星历的平均精度也在0.5ns左右,证明了本发明预测方法的有效性。G01卫星属于BLOCK IIF最新一代的GPS导航卫星,其星载RB钟稳定性较好,钟差测量数据精度较高,相应的钟差预报精度也最高。G31卫星这一时段的钟差测量精度不高,钟差数据有所起伏,因此预报精度相对于其他RB钟卫星稍差。CS钟比RB钟稳定性差,其钟差测量精度最低,因此相应的G27卫星钟差预报精度也低一些,但残差也稳定在1ns以内,证明PGM(1,1)模型有较好的适应性。同时也可说明预测序列的变化趋势不影响钟差预报的精度,而序列的平稳程度对精度有较大的影响。
表2PGM(1,1)模型预测精度统计表
如表2和表3所示,比较两种模型的预报结果可知,PGM(1,1)模型比GM(1,1)模型平均相对误差降低80%至98%,平均预报残差也从纳秒量级降到了亚纳秒级。GM(1,1)模型在预报G01卫星钟差时出现了误差发散的情况,说明传统GM(1,1)模型对不同序列的敏感性有很大差异。而PGM(1,1)模型对G01卫星钟差预报效果极佳,证明PSO算法的引入改善了传统模型的适应性和敛散性,在各数据情况下提供较高精度的预测结果。
表3GM(1,1)模型预报精度统计
与现有技术相比,本发明能够提高卫星钟差的预报精度,进而辅助精密单点定位接收机进行高精度实时定位解算,PGM(1,1)预报模型所需训练样本小,模型参数通过自适应方法调节,计算量较小,便于工程运用,且卫星钟差预报精度较高,误差控制在1ns以内。
Claims (6)
1.一种灰色粒子群卫星钟差自适应预测方法,其特征在于,通过在GM(1,1)模型中引入优化因子以及遗忘因子最小二乘法建立灰色粒子群模型,然后采用粒子群优化算法求得最佳优化因子,最终将最佳优化因子代入PGM(1,1)模型并进行钟差预测得到钟差的预测序列;
所述的方法具体包括以下步骤:
1)对钟差的原始序列进行处理得到非负的等维序列x(0)={x(0)(1),x(0)(2),..x(0)(k).,x(0)(n)}其中:k=1,2,...,n;
2)对等维序列进行累加得到一次累加序列x(1) = {x(1)(1),x(1)(2),..x(1)(k).,x(1)(n)} , 其中:
3)建立包含优化因子和遗忘因子最小二乘法的PGM(1,1)模型,具体包括以下步骤:
3.1)通过对1-AGO序列进行一次紧邻数据自适应寻优生成背景值z(1)={z(1)(2),z(1)(3),...,z(1)(n)},其中:z(1)(k)=ωx(1)(k)+(1-ω)x(1)(k-1),k=2,...,n,0≤ω≤1;
3.2)建立灰微分方程x(0)(k)=az(1)(k)+b后,得到白化方程其中:a,b为模型待辨识参数;
3.3)得到离散化后的初始解 其中:k=1,2,...,n,1≤l≤k;
3.4)引入发展系数λ和灰作用量θ后获得预测值k=1,2,...,n;
所述的优化因子包括:自适应权值ω、发展系数λ和灰作用量θ;
4)采用PSO算法以历史钟差作为样本求解出最佳优化因子;
5)将最佳优化因子代入PGM(1,1)模型;
6)计算钟差的预测值。
2.根据权利要求1所述的灰色粒子群卫星钟差自适应预测方法,其特征是,所述的发展系数λ∈[0.5,1.5],灰作用量θ∈(0.9,1.1],且发展系数θb∈[0,1)。
3.根据权利要求2所述的灰色粒子群卫星钟差自适应预测方法,其特征是,所述的模型待辨识参数[a b]T=(LTWL)-1LTWY,其中:遗忘因子W=diag(e-(n-1),e-(n-2),...,e-(n-n)),Y=[x(0)(2),x(0)(3),...,x(0)(n)]T,
4.根据权利要求3所述的灰色粒子群卫星钟差自适应预测方法,其特征是,所述的步骤4)具体以下步骤:
4.1)生成初始群体(φ1,φ2,...,φn),其中:φi=[ωi,λi,θi]T;
4.2)以归一化平均相对误差最小为准则,比较每个粒子当前适应值和其个体历史最好适应值,获得全局最好适应值;
4.3)更新各微粒的速度和位置;
4.3)当适应值误差达到设定误差限或到达最大迭代次数则停止搜索并输出最佳优化因子,否则回到步骤4.2)。
5.根据权利要求4所述的灰色粒子群卫星钟差自适应预测方法,其特征是,所述的归一化平均相对误差k=1,2,...,n。
6.一种实现上述任一权利要求所述方法的系统,其特征在于,包括:网络数据更新模块、原始数据预处理模块、卫星钟差预报模块以及数据预警模块,其中:网络数据更新模块与原始数据预处理模块相连并传输快速精密星历数据;
所述的原始数据预处理模块提取钟差数据截取跳跃点后的钟差数据作为原始钟差序列;
所述的卫星钟差预报模块与原始数据预处理模块相连,并接收预处理后原始钟差序列;
所述的卫星钟差预报模块与数据预警模块相连,并传输钟差预测值,数据预警模块通过比对钟差预测值和实际钟差值,评估系统预报结果的可用性和准确性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610605339.3A CN106126961B (zh) | 2016-07-28 | 2016-07-28 | 灰色粒子群卫星钟差自适应预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610605339.3A CN106126961B (zh) | 2016-07-28 | 2016-07-28 | 灰色粒子群卫星钟差自适应预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106126961A CN106126961A (zh) | 2016-11-16 |
CN106126961B true CN106126961B (zh) | 2018-07-27 |
Family
ID=57254279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610605339.3A Active CN106126961B (zh) | 2016-07-28 | 2016-07-28 | 灰色粒子群卫星钟差自适应预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106126961B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108387915B (zh) * | 2018-03-05 | 2021-06-01 | 中国人民解放军空军工程大学 | 一种融合钟差预报方法 |
CN109034491A (zh) * | 2018-08-15 | 2018-12-18 | 长沙学院 | 一种卫星钟差预测方法 |
CN110187630B (zh) * | 2019-05-16 | 2020-11-03 | 中国人民解放军空军工程大学 | 预报时长不确定条件下的短期组合钟差预报方法 |
CN112630805B (zh) * | 2019-09-24 | 2022-06-28 | 千寻位置网络有限公司 | 卫星钟差预报方法及其系统 |
CN111626343B (zh) * | 2020-05-13 | 2022-05-03 | 哈尔滨工程大学 | 一种基于pgm与pso聚类的船舰数据关系抽取方法 |
CN112329197A (zh) * | 2020-09-23 | 2021-02-05 | 北京无线电计量测试研究所 | 一种基于灰色模型的综合原子时建立方法 |
CN112579580B (zh) * | 2020-11-26 | 2022-07-05 | 广东工业大学 | 一种基于工业大数据预测的预报警方法 |
CN115529606A (zh) * | 2021-06-25 | 2022-12-27 | 中国移动通信集团吉林有限公司 | 参数更新方法、系统以及电子设备 |
CN113570120A (zh) * | 2021-07-07 | 2021-10-29 | 上海电机学院 | 基于改进粒子群算法的风电功率预测方法及装置 |
CN114265094B (zh) * | 2021-12-06 | 2023-07-18 | 中国人民解放军海军工程大学 | 基于初始条件优化gm模型的gps卫星钟差预报方法及系统 |
CN114265093B (zh) * | 2021-12-06 | 2023-07-18 | 中国人民解放军海军工程大学 | 一种卫星时钟钟差预报方法、系统、介质、设备及终端 |
CN116050483B (zh) * | 2023-01-10 | 2024-02-13 | 辽宁工程技术大学 | 一种ssa-bp神经网络与灰色模型卫星钟差预报方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105117593A (zh) * | 2015-08-18 | 2015-12-02 | 河海大学 | 基于小波变换和粒子群改进灰色模型的短期风速预测方法 |
-
2016
- 2016-07-28 CN CN201610605339.3A patent/CN106126961B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105117593A (zh) * | 2015-08-18 | 2015-12-02 | 河海大学 | 基于小波变换和粒子群改进灰色模型的短期风速预测方法 |
Non-Patent Citations (5)
Title |
---|
GM(1,1)模型的背景值构造方法和应用(I);谭冠军;《系统工程理论与实践》;20000430;第20卷(第4期);98-103 * |
Grey particle swarm optimization;Min-Shyang Leu等;《Applied Soft Computing》;20120930;第12卷(第9期);2985-2996 * |
The parameter estimation of time-varying GM(1,1);Shih N.Y.等;《Journal of Grey System》;20060131;第9卷(第1期);51-56 * |
基于灰色系统模型的IGS精密钟差预报;潘绍林等;《测绘工程》;20150430;第24卷(第4期);20-24 * |
灰色预测技术研究进展综述;党耀国等;《上海电机学院学报》;20150131;第18卷(第1期);1-8 * |
Also Published As
Publication number | Publication date |
---|---|
CN106126961A (zh) | 2016-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106126961B (zh) | 灰色粒子群卫星钟差自适应预测方法及系统 | |
Chen et al. | Novel hybrid of strong tracking Kalman filter and wavelet neural network for GPS/INS during GPS outages | |
CN110058236A (zh) | 一种面向三维地表形变估计的InSAR和GNSS定权方法 | |
CN110232471B (zh) | 一种降水传感网节点布局优化方法及装置 | |
CN111680870B (zh) | 目标运动轨迹质量综合评估方法 | |
Mo et al. | Estimating the extreme wind speed for regions in China using surface wind observations and reanalysis data | |
CN104048676B (zh) | 基于改进粒子滤波的mems陀螺随机误差补偿方法 | |
CN108508457B (zh) | 基于精密星历的阻力系数自适应调制方法 | |
CN111649763A (zh) | 一种基于重力灯塔建立的潜艇导航方法及系统 | |
CN104915534A (zh) | 基于序列学习的电力铁塔变形分析与决策方法 | |
CN110990505A (zh) | 一种基于神经网络的Loran-C ASF修正方法 | |
CN106770620B (zh) | 用溅射深度剖析技术确定薄膜中元素成分深度分布的方法 | |
CN112686481A (zh) | 一种径流量预报方法、处理器 | |
CN110532621A (zh) | 一种飞行器气动参数在线辨识方法 | |
CN110018501A (zh) | 一种基于系统间随机模型在线估计调整的多模精密单点定位方法 | |
CN110489879A (zh) | 一种适用于空间环境扰动情况下的空间目标陨落预报方法 | |
CN114970341B (zh) | 基于机器学习的低轨卫星轨道预报精度提升模型建立方法 | |
CN112131752B (zh) | 一种基于拟准检定的超强崩溃污染率抗差估计算法 | |
CN103793614B (zh) | 一种突变滤波方法 | |
CN105447228A (zh) | 一种地球扰动引力赋值模型适用性能评估方法 | |
CN110598171A (zh) | 基于Shannon熵-能量比的小波神经网络卫星钟差预报方法 | |
CN113642785B (zh) | 基于先验信息的空间碎片轨道长期预报方法、系统及设备 | |
CN112434261B (zh) | 基于标校卫星的测控装备精度鉴定方法 | |
Cheng et al. | An OW-FCE model based on MDE algorithm for evaluating integrated navigation system | |
Ma et al. | A Method for Establishing Tropospheric Atmospheric Refractivity Profile Model Based on Multiquadric RBF and k-means Clustering |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |