CN109581524A - 一种旋转加速度计式重力梯度敏感器动态测量解调方法 - Google Patents

一种旋转加速度计式重力梯度敏感器动态测量解调方法 Download PDF

Info

Publication number
CN109581524A
CN109581524A CN201811406506.7A CN201811406506A CN109581524A CN 109581524 A CN109581524 A CN 109581524A CN 201811406506 A CN201811406506 A CN 201811406506A CN 109581524 A CN109581524 A CN 109581524A
Authority
CN
China
Prior art keywords
signal
gravity gradient
gradient sensor
dynamic measurement
demodulation method
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
Application number
CN201811406506.7A
Other languages
English (en)
Other versions
CN109581524B (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.)
707th Research Institute of CSIC
Original Assignee
707th Research Institute of CSIC
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 707th Research Institute of CSIC filed Critical 707th Research Institute of CSIC
Priority to CN201811406506.7A priority Critical patent/CN109581524B/zh
Publication of CN109581524A publication Critical patent/CN109581524A/zh
Application granted granted Critical
Publication of CN109581524B publication Critical patent/CN109581524B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种旋转加速度计式重力梯度敏感器动态测量解调方法，其技术特点在于：包括以下步骤：步骤1、对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号进行离散采样；步骤2、通过解调参考信号sin2ωt和cos2ωt将敏感器原始输出信号进行频谱搬移；步骤3、对步骤2得到的信号进行降采样；步骤4、将m组不同相位的信号(Γxx‑Γyy)gk和(Γxy)gk分别进行低通滤波；步骤5、将m组信号(Γxx‑Γyy)fk和(Γxy)fk的每点数据进行等权重融合。本发明可有效抑制系统输出信号中与位置无关的噪声，提高系统动态测量精度。

Description

一种旋转加速度计式重力梯度敏感器动态测量解调方法
技术领域
本发明属于重力梯度敏感器技术领域，涉及重力梯度敏感器动态测量解调方法，尤其是一种旋转加速度计式重力梯度敏感器动态测量解调方法。
背景技术
重力梯度仪是用于测量重力场分布梯度的精密设备，基于Bell Aerospace公司提出的旋转加速度计测量原理的重力梯度仪是迄今唯一实用的近地表动态测量重力梯度仪。碍于重力梯度信号的微弱和加速度计性能的限制，梯度仪输出信号的信噪比极低，需通过解调、滤波等数据处理手段进一步提高测量精度。
如图2所示，旋转加速度计式重力梯度仪基于加速度计位置差分测量原理，通过机械旋转的方式将重力梯度张量分量调制到系统旋转频率的二倍频处，加速度计输出与重力梯度张量分量之间的关系可表示为：
(a1+a3)-(a2+a4)＝2R[(Γxx-Γyy)sin(2at+θ)-2Γxycos(2at+θ)]
式中，ω为圆盘旋转速率；R为加速度计中心到圆盘中心的距离；θ为初始相位。最后，加速度计组合信号通过2ω频率解调得到最终的重力梯度张量测量信号Γxx-Γyy和Γxy。
传统的重力梯度敏感器解调方法基于整周期的正弦解调方法，在数字系统中其解调原理公式为：
式中Aout(i)是经位置同步的加速度计输出采样值，N是重力梯度仪圆盘旋转半周时采样数据个数，ts是梯度仪的采样时间间隔，2是将解调值转化为幅值的比例系数。
上述方法虽可准确解调出重力梯度信号的零频成分，但其滚动平滑滤波的特性对重力梯度动态测量有效频段内会随着频率的增加会不同程度的衰减，滚动平滑滤波幅频响应如图3所示，在航空重力梯度测量关心的0.01-0.1Hz有效频带内，可存在0.5dB的衰减，这在高精度动态重力梯度测量中是不可接受的。同时在有效测量频带外，该解调方式又无法满足重力梯度测量对动态噪声抑制比的要求，需在数据处理环节中增加低通滤波以满足系统要求。
因此，传统的解调方法并不适合重力梯度动态测量，需设计一种适应动态测量的旋转加速度计式重力梯度敏感器解调方法。
发明内容
本发明的目的在于克服现有技术的不足，提出一种旋转加速度计式重力梯度敏感器动态测量的解调方法，在保证重力梯度有效测量频带内的信号不受衰减的前提下最大程度抑制系统原始输出信号重力梯度调制频率附近的不相关噪声，提高系统动态测量精度。
本发明解决其现实问题是采取以下技术方案实现的：
一种旋转加速度计式重力梯度敏感器动态测量解调方法，包括以下步骤：
步骤1、对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号进行离散采样；
步骤2、通过解调参考信号sin2ωt和cos2ωt将敏感器原始输出信号进行频谱搬移，得到两组相位正交的信号(Γxx-Γyy)ci和(Γxy)ci；
步骤3、对步骤2得到的信号(Γxx-Γyy)ci和(Γxy)ci进行降采样；
步骤4、将m组不同相位的信号(Γxx-Γyy)gk和(Γxy)gk分别进行低通滤波；
步骤5、将m组信号(Γxx-Γyy)fk和(Γxy)fk的每点数据进行等权重融合，得到解调后的重力梯度敏感器输出重力梯度信号(Γxx-Γyy)j和(Γxy)j，j＝1,2,3……，N/m；
而且，所述步骤1的具体方法为：对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号从初始相位角为0的时刻开始进行离散采样，采样时间间隔为Δt，得到敏感器原始输出信号Aout(i)和与之对应的时间ti，其中i＝1,2,3，……N；
而且，所述步骤2的两组相位正交的信号(Γxx-Γyy)ci和(Γxy)ci的计算公式为：
(Γxx-Γyy)ci＝Aout(i)·sin(2ω·ti)
(Γxy)ci＝Aout(i)·cos(2ω·ti)
而且，所述步骤3的具体步骤包括：
(1)降采样比例为m(应保证降采样后的信号采样频率大于1即m≤fs)，所用低通滤波器参数为：通带上限频率阻带下限频率通带容限ap＝0.001dB，阻带容限as＝60dB，滤波器算法为雷米兹交替算法；其中fs为原始信号采样频率；
(2)将两组相位正交的信号(Γxx-Γyy)ci和(Γxy)ci经过上述低通滤波器，得到信号(Γxx-Γyy)di和(Γxy)di，对两路信号等间隔抽样，抽样间隔为m，得到m组不同的信号(Γxx-Γyy)gk和(Γxy)gk，其中抽样相位k的取值为[0,m-1]；
而且，所述步骤4的具体方法为：
所用低通滤波器参数为：通带上限频率wp＝0.1Hz，阻带下限频率ws＝0.12Hz，通带容限ap＝0.001dB，阻带容限as＝60dB，滤波器算法为雷米兹交替算法，得到m组低通滤波后的信号(Γxx-Γyy)fk和(Γxy)fk；
而且，所述步骤5的计算公式为：
本发明的优点和有益效果：
1、本发明在低通滤波器设计环节中采用雷米兹交替算法，通过对通带容限的参数设定可有效保证在动态重力梯度有效测量频带(0.01-0.1Hz)内，重力梯度信号衰减程度小于0.005dB；
2、本发明在降采样环节中采用等间隔抽样的方式建立多组数据，并将其等权重融合的方法可有效抑制系统输出信号中与位置无关的噪声，提高系统动态测量精度。
附图说明
图1为本发明的处理流程图；
图2为本发明的旋转加速度计式重力梯度敏感器原理示意图；
图3为本发明的滚动平滑滤波幅频响应结果图。
具体实施方式
以下结合附图对本发明实施例作进一步详述：
一种旋转加速度计式重力梯度敏感器动态测量解调方法，如图1所示，包括以下步骤：
步骤1、对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号进行离散采样；
所述步骤1的具体方法为：对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号从初始相位角为0的时刻开始进行离散采样，采样时间间隔为Δt，得到敏感器原始输出信号Aout(i)和与之对应的时间ti，其中i＝1,2,3，……N；
步骤2、通过解调参考信号sin2ωt和cos2ωt将敏感器原始输出信号进行频谱搬移，得到两组相位正交的信号(Γxx-Γyy)ci和(Γxy)ci；
所述步骤2的两组相位正交的信号(Γxx-Γyy)ci和(Γxy)ci的计算公式为：
(Γxx-Γyy)ci＝Aout(i)·sin(2ω·ti)
(Γxy)ci＝Aout(i)·cos(2ω·ti)
步骤3、对步骤2得到的信号(Γxx-Γyy)ci和(Γxy)ci进行降采样；
所述步骤3的具体步骤包括：
(1)降采样比例为m(应保证降采样后的信号采样频率大于1即m≤fs)，所用低通滤波器参数为：通带上限频率阻带下限频率通带容限ap＝0.001dB，阻带容限as＝60dB，滤波器算法为雷米兹交替算法；其中fs为原始信号采样频率；
(2)将两组相位正交的信号(Γxx-Γyy)ci和(Γxy)ci经过上述低通滤波器，得到信号(Γxx-Γyy)di和(Γxy)di，对两路信号等间隔抽样，抽样间隔为m，得到m组不同的信号(Γxx-Γyy)gk和(Γxy)gk，其中抽样相位k的取值为[0,m-1]；
步骤4、将m组不同相位的信号(Γxx-Γyy)gk和(Γxy)gk分别进行低通滤波；
所述步骤4的具体方法为：
所用低通滤波器参数为：通带上限频率wp＝0.1Hz，阻带下限频率ws＝0.12Hz，通带容限ap＝0.001dB，阻带容限as＝60dB，滤波器算法为雷米兹交替算法，得到m组低通滤波后的信号(Γxx-Γyy)fk和(Γxy)fk；
步骤5、将m组信号(Γxx-Γyy)fk和(Γxy)fk的每点数据进行等权重融合，得到解调后的重力梯度敏感器输出重力梯度信号(Γxx-Γyy)j和(Γxy)j，j＝1,2,3……，N/m；
所述步骤5的计算公式为：
本发明的工作原理是：
基于同步解调的原理，将传统的时域平滑改为低通滤波，在滤波器算法采用雷米兹交替方法，通过对通带容限的约束保证重力梯度动态测量下信号衰减小于0.005dB。并在低通滤波前增加降采样环节，并在抽取过程中依据相位的不同抽取成多组只与位置有关的数据，同时低通滤波成多组重力梯度量测信号，将上述多组量测信号等权重融合成最后输出结果以抑制与位置无关的噪声，提高系统动态测量精度。
需要强调的是，本发明所述实施例是说明性的，而不是限定性的，因此本发明包括并不限于具体实施方式中所述实施例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，同样属于本发明保护的范围。

Claims (6)

1.一种旋转加速度计式重力梯度敏感器动态测量解调方法，其特征在于：包括以下步骤：
步骤1、对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号进行离散采样；
步骤2、通过解调参考信号sin2ωt和cos2ωt将敏感器原始输出信号进行频谱搬移，得到两组相位正交的信号(Γxx-Γyy)ci和(Γxy)ci；
步骤3、对步骤2得到的信号(Γxx-Γyy)ci和(Γxy)ci进行降采样；
步骤4、将m组不同相位的信号(Γxx-Γyy)gk和(Γxy)gk分别进行低通滤波；
步骤5、将m组信号(Γxx-Γyy)fk和(Γxy)fk的每点数据进行等权重融合，得到解调后的重力梯度敏感器输出重力梯度信号(Γxx-Γyy)j和(Γxy)j，j＝1,2,3……，N/m。
2.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器动态测量解调方法，其特征在于：所述步骤1的具体方法为：对工作中的旋转加速度计式重力梯度敏感器中加速度计组合输出模拟信号从初始相位角为0的时刻开始进行离散采样，采样时间间隔为Δt，得到敏感器原始输出信号Aout(i)和与之对应的时间ti，其中i＝1,2,3，……N。
3.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器动态测量解调方法，其特征在于：所述步骤2的两组相位正交的信号(Γxx-Γyy)ci和(Γxy)ci的计算公式为：
(Γxx-Γyy)ci＝Aout(i)·sin(2ω·ti)
(Γxy)ci＝Aout(i)·cos(2ω·ti)
4.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器动态测量解调方法，其特征在于：所述步骤3的具体步骤包括：
(1)降采样比例为m，应保证降采样后的信号采样频率大于1即m≤fs，所用低通滤波器参数为：通带上限频率阻带下限频率通带容限ap＝0.001dB，阻带容限as＝60dB，滤波器算法为雷米兹交替算法；其中fs为原始信号采样频率；
(2)将两组相位正交的信号(Γxx-Γyy)ci和(Γxy)ci经过上述低通滤波器，得到信号(Γxx-Γyy)di和(Γxy)di，对两路信号等间隔抽样，抽样间隔为m，得到m组不同的信号(Γxx-Γyy)gk和(Γxy)gk，其中抽样相位k的取值为[0,m-1]。
5.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器动态测量解调方法，其特征在于：所述步骤4的具体方法为：
所用低通滤波器参数为：通带上限频率wp＝0.1Hz，阻带下限频率ws＝0.12Hz，通带容限ap＝0.001dB，阻带容限as＝60dB，滤波器算法为雷米兹交替算法，得到m组低通滤波后的信号(Γxx-Γyy)fk和(Γxy)fk。
6.根据权利要求1所述的一种旋转加速度计式重力梯度敏感器动态测量解调方法，其特征在于：所述步骤5的计算公式为：
CN201811406506.7A 2018-11-23 2018-11-23 一种旋转加速度计式重力梯度敏感器动态测量解调方法 Active CN109581524B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811406506.7A CN109581524B (zh) 2018-11-23 2018-11-23 一种旋转加速度计式重力梯度敏感器动态测量解调方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811406506.7A CN109581524B (zh) 2018-11-23 2018-11-23 一种旋转加速度计式重力梯度敏感器动态测量解调方法

Publications (2)

Publication Number Publication Date
CN109581524A true CN109581524A (zh) 2019-04-05
CN109581524B CN109581524B (zh) 2020-06-19

Family

ID=65924220

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811406506.7A Active CN109581524B (zh) 2018-11-23 2018-11-23 一种旋转加速度计式重力梯度敏感器动态测量解调方法

Country Status (1)

Country Link
CN (1) CN109581524B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111624671A (zh) * 2020-06-30 2020-09-04 东南大学 旋转加速度计重力梯度仪重力梯度解调相位角确定方法及装置
CN111650664A (zh) * 2020-06-30 2020-09-11 东南大学 一种航空重力梯度仪实时重力梯度解调方法及装置
CN111830589A (zh) * 2020-06-18 2020-10-27 中国船舶重工集团公司第七0七研究所 旋转加速度计重力梯度敏感器冗余测量信号离线融合方法
CN112327379A (zh) * 2020-09-28 2021-02-05 中国船舶重工集团公司第七0七研究所 一种全张量重力梯度动态测量系统及方法
CN113433596A (zh) * 2021-06-25 2021-09-24 中国船舶重工集团公司第七0七研究所 一种基于空间域的重力梯度动态测量滤波方法
CN113447994A (zh) * 2021-06-25 2021-09-28 中国船舶重工集团公司第七0七研究所 一种基于旋转加速度计原理的重力梯度全张量测量方法
CN113885099A (zh) * 2021-09-28 2022-01-04 中国船舶重工集团公司第七0七研究所 重力梯度仪加速度计标度因数不一致性动态实时估计方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5357802A (en) * 1993-05-07 1994-10-25 Textron, Incorporated Rotating accelerometer gradiometer
CN104898176A (zh) * 2015-06-10 2015-09-09 东南大学 一种旋转加速度计重力梯度仪重力梯度解调方法
US20150301077A1 (en) * 2014-04-16 2015-10-22 Cirrus Logic, Inc. Systems and methods for determining acceleration based on phase demodulation of an electrical signal
CN105044798A (zh) * 2015-06-29 2015-11-11 东南大学 旋转加速度计重力梯度仪加速度计标度因子反馈调整方法
US20160169935A1 (en) * 2013-01-28 2016-06-16 The Regents Of The University Of California Multi-Axis Chip-Scale MEMS Inertial Measurement Unit (IMU) Based on Frequency Modulation
CN106526235A (zh) * 2016-10-28 2017-03-22 西安交通大学 一种低g值电容式MEMS加速度计及其模态局域化测量电路
CN106761711A (zh) * 2016-12-26 2017-05-31 中石化江汉石油工程有限公司 一种井下仪器多功能通讯控制系统
CN108287372A (zh) * 2017-12-21 2018-07-17 中国船舶重工集团公司第七0七研究所 一种基于惯性技术的重力梯度敏感器转速选择方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5357802A (en) * 1993-05-07 1994-10-25 Textron, Incorporated Rotating accelerometer gradiometer
US20160169935A1 (en) * 2013-01-28 2016-06-16 The Regents Of The University Of California Multi-Axis Chip-Scale MEMS Inertial Measurement Unit (IMU) Based on Frequency Modulation
US20150301077A1 (en) * 2014-04-16 2015-10-22 Cirrus Logic, Inc. Systems and methods for determining acceleration based on phase demodulation of an electrical signal
CN104898176A (zh) * 2015-06-10 2015-09-09 东南大学 一种旋转加速度计重力梯度仪重力梯度解调方法
CN105044798A (zh) * 2015-06-29 2015-11-11 东南大学 旋转加速度计重力梯度仪加速度计标度因子反馈调整方法
CN106526235A (zh) * 2016-10-28 2017-03-22 西安交通大学 一种低g值电容式MEMS加速度计及其模态局域化测量电路
CN106761711A (zh) * 2016-12-26 2017-05-31 中石化江汉石油工程有限公司 一种井下仪器多功能通讯控制系统
CN108287372A (zh) * 2017-12-21 2018-07-17 中国船舶重工集团公司第七0七研究所 一种基于惯性技术的重力梯度敏感器转速选择方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨晔 等: "旋转加速度计式重力梯度仪输出解调与滤波", 《中国惯性技术学报》 *
钱学武 等: "一种旋转加速度计重力梯度仪重力梯度解调方法", 《物探与化探》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111830589A (zh) * 2020-06-18 2020-10-27 中国船舶重工集团公司第七0七研究所 旋转加速度计重力梯度敏感器冗余测量信号离线融合方法
CN111830589B (zh) * 2020-06-18 2022-07-26 中国船舶重工集团公司第七0七研究所 旋转加速度计重力梯度敏感器冗余测量信号离线融合方法
CN111624671A (zh) * 2020-06-30 2020-09-04 东南大学 旋转加速度计重力梯度仪重力梯度解调相位角确定方法及装置
CN111650664A (zh) * 2020-06-30 2020-09-11 东南大学 一种航空重力梯度仪实时重力梯度解调方法及装置
CN112327379A (zh) * 2020-09-28 2021-02-05 中国船舶重工集团公司第七0七研究所 一种全张量重力梯度动态测量系统及方法
CN112327379B (zh) * 2020-09-28 2022-04-12 中国船舶重工集团公司第七0七研究所 一种全张量重力梯度动态测量系统及方法
CN113433596A (zh) * 2021-06-25 2021-09-24 中国船舶重工集团公司第七0七研究所 一种基于空间域的重力梯度动态测量滤波方法
CN113447994A (zh) * 2021-06-25 2021-09-28 中国船舶重工集团公司第七0七研究所 一种基于旋转加速度计原理的重力梯度全张量测量方法
CN113433596B (zh) * 2021-06-25 2022-09-16 中国船舶重工集团公司第七0七研究所 一种基于空间域的重力梯度动态测量滤波方法
CN113885099A (zh) * 2021-09-28 2022-01-04 中国船舶重工集团公司第七0七研究所 重力梯度仪加速度计标度因数不一致性动态实时估计方法
CN113885099B (zh) * 2021-09-28 2024-02-27 中国船舶重工集团公司第七0七研究所 重力梯度仪加速度计标度因数不一致性动态实时估计方法

Also Published As

Publication number Publication date
CN109581524B (zh) 2020-06-19

Similar Documents

Publication Publication Date Title
CN109581524A (zh) 一种旋转加速度计式重力梯度敏感器动态测量解调方法
CN104898176B (zh) 一种旋转加速度计重力梯度仪重力梯度解调方法
WO2020140378A1 (zh) 一种旋转加速度计重力梯度仪运动误差事后补偿方法
CN105044798A (zh) 旋转加速度计重力梯度仪加速度计标度因子反馈调整方法
CN109856689B (zh) 一种超导航磁梯度张量数据抑噪处理方法和系统
CN105486331B (zh) 一种具有高精度的光学信号相位解调系统及解调方法
CN112487604B (zh) 海洋重力仪输出数据长时间非线性漂移补偿方法
CN109782023B (zh) 一种通过旋转调制法测量加速度计高阶项数系数的方法
WO2020042696A1 (zh) 一种动基座旋转加速度计重力梯度仪误差补偿装置及方法
CN104570142B (zh) 一种重力梯度仪重力梯度测量信号的解调方法
CN108287372A (zh) 一种基于惯性技术的重力梯度敏感器转速选择方法
CN107576992B (zh) 一种重力梯度仪自标定方法及离心梯度补偿方法
CN110471123A (zh) 一种旋转加速度计重力梯度仪数据诊断及处理方法
CN106705995A (zh) 一种MEMS陀螺仪g值敏感系数的标定方法
CN109085654B (zh) 一种旋转加速度计重力梯度仪数字建模仿真方法
CN109581525B (zh) 旋转加速度计式重力梯度敏感器原始采样频率的选择方法
CN108120439A (zh) 一种三分量感应线圈姿态测量方法及装置
CN104730473B (zh) 绝对磁场测量设备及所适用的绝对磁场测量方法
CN108716925A (zh) 一种九轴传感器的标定方法及装置
CN104569496B (zh) 一种使用引力梯度测试加速度计分辨率的方法
CN109033646B (zh) 基于小量程mems加速度传感器的随钻测斜仪重力工具面角计算方法
CN109212629B (zh) 旋转加速度计重力梯度仪角运动误差补偿装置及方法
CN113447994B (zh) 一种基于旋转加速度计原理的重力梯度全张量测量方法
CN112906482A (zh) 基于随机共振与混沌协同的井下微弱特征信号检测方法
CN111679334A (zh) 一种绝对重力加速度的补偿方法和设备

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