CN110068876B - 基于载体自振动航空重力梯度仪运动误差补偿方法 - Google Patents
基于载体自振动航空重力梯度仪运动误差补偿方法 Download PDFInfo
- Publication number
- CN110068876B CN110068876B CN201910464038.7A CN201910464038A CN110068876B CN 110068876 B CN110068876 B CN 110068876B CN 201910464038 A CN201910464038 A CN 201910464038A CN 110068876 B CN110068876 B CN 110068876B
- Authority
- CN
- China
- Prior art keywords
- accelerometer
- vertical
- signals
- frequency
- signal
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Navigation (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种基于载体自振动航空重力梯度仪运动误差补偿方法,其技术特点在于:包括以下步骤:步骤1、对加速度计组合信号在飞机螺旋桨旋转频率上进行正余弦解调;步骤2、对垂向加速度计输出信号在飞机螺旋桨旋转频率上进行正余弦解调;步骤3、求得两组加速度计敏感轴垂向安装误差实时信号;步骤4、对垂向加速度计输出信号在二倍频频点上进行正余弦解调;步骤5、得到重力梯度信号中加速度计敏感轴俯仰角安装误差予以补偿。本发明合理补偿航空重力梯度仪的垂向加速度误差,提高了系统的测量精度。
Description
技术领域
本发明属于旋转加速度计式航空重力梯度仪技术领域,涉及旋转加速度计式航空重力梯度仪垂向运动误差补偿方法,尤其是一种基于载体自振动航空重力梯度仪运动误差补偿方法。
背景技术
目前,旋转加速度计式重力梯度仪是对地球表面微小重力梯度变化进行连续测量的仪器。如图2所示,作为核心敏感器的重力梯度测量组件基于加速度计位置差分测量原理,通过机械旋转的方式将重力梯度张量分量调制到系统旋转频率的二倍频处,加速度计输出与重力梯度张量分量之间的关系可表示为:
(a1+a3)-(a2+a4)=2R(Γxx-Γyy)sin2ωt+4RΓxycos2ωt (1)
式中a1、a2、a3、a4是四个加速度计的输出,R是加速度计检测质心到旋转中心的距离,Γxx、Γyy、Γxy是重力梯度张量分量,ω是旋转装置的旋转角速度。将系统输出的加速度计组合信号进行2ω频率上的同步解调得到最终的重力梯度张量信号Γxx-Γyy和Γxy。同时定义Γuv=Γxx-Γyy。进行航空测量时,由惯性稳定平台承载重力梯度敏感器,隔离载体角运动对重力梯度测量的影响,同时为敏感器提供基准坐标系。
利用飞机等机动平台开展重力梯度勘查具有效率高、成本低、环境适应性强的优点,可在人员难以到达的复杂地理环境中部署使用。但在机载测量条件下,飞机在飞行过程中垂向振动的二倍频分量与四只加速度计的敏感轴垂向安装误差角相耦合,进入加速度计组合信号的二倍频分量上,对重力梯度测量信号造成干扰。经测量,机载垂向运动加速度的二倍频分量约为2×10-3g量级,在每只加速度计敏感轴垂向安装误差精度优于5″的条件下,对重力梯度信号的干扰为5000E。因此,机载条件下对载体垂向加速度误差的补偿是影响仪器测量性能的重要环节,需通过特定的方法予以补偿。
发明内容
本发明的目的在于克服现有技术的不足,提出一种基于载体自振动航空重力梯度仪运动误差补偿方法,能够合理补偿旋转加速度计式重力梯度仪在机载测量时,输出的重力梯度信号中存在由载体垂向加速度与加速度计敏感轴垂向安装误差相耦合引起的测量误差。
本发明解决其现实问题是采取以下技术方案实现的:
一种旋转加速度计式航空重力梯度仪垂向线运动测量误差的补偿方法,包括以下步骤:
步骤3、将由从步骤1的加速度计组合信号得到的和信号与从步骤2的垂向加速度计输出信号得到的和信号进行相除处理,得到两组加速度计敏感轴垂向安装误差实时信号Δβs和Δβc,将两组信号Δβs和Δβc求平均得到更加真实的Δβ实时信号。其中Δβs是由飞机螺旋桨旋转频率正弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβc是由飞机螺旋桨旋转频率余弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβ是两个安装误差实时信号的平均值;
本发明的优点和有益效果:
1、本发明能够依据固定翼飞机飞行测量时的运动特点和重力梯度测量的误差机理,利用本是误差源的飞机螺旋桨振动作为加速度计敏感轴垂向安装误差的观测激励项,对其进行实时监测,并依据该项安装误差实时监测值对重力梯度测量信号中的垂向加速度误差进行补偿。
2、本发明针对动态飞行测量的特点和误差原理,提出利用本是误差源的飞机自振动作为观测激励项,合理补偿航空重力梯度仪的垂向加速度误差,提高了系统的测量精度;
3、本发明在不依靠外界信息、不增加系统的复杂程度的基础上,仅使用系统稳定平台的惯性导航数据对重力梯度信号中的垂向加速度误差进行补偿,提高了航空重力梯度仪系统的可靠性。
附图说明
图1为本发明的机载试验中飞机平飞过程中的垂向振动谱示意图;
图2为本发明的旋转加速度计式重力梯度敏感器测量原理示意图。
具体实施方式
以下结合附图对本发明实施例作进一步详述:
一种基于载体自振动航空重力梯度仪运动误差补偿方法,可以用于补偿机载测量时由固定翼飞机垂向运动引起的重力梯度测量误差。
在本实施例中,机载条件下重力梯度仪由加速度计敏感轴垂向安装误差带来的重力梯度敏感器输出信号误差为:
ΔAout(Δβi)=a0z[(Δβ1+Δβ3)-(Δβ2+Δβ4)]
+[Γxz(Δβ2-Δβ4)-Γyz(Δβ1-Δβ3)]sinωt
+[Γxz(Δβ1-Δβ3)+Γyz(Δβ2-Δβ4)]cosωt (2)
式中,ΔAout(Δβi)是由加速度计敏感轴垂向安装误差引起的重力梯度敏感器输出加速度计组合信号的误差,βi(i=1,2,3,4)是第i个加速度计敏感轴的垂向角度安装误差,Γxz、Γyz是对应方向的重力梯度张量,a0z是载体的垂向加速度。
从中可以看出,误差第二项、第三项其载波频率为一倍频,可由同步解调技术除去,误差第一项表示载体垂向加速度的二倍频分量与敏感轴垂向安装误差角相耦合,进入重力梯度测量通道。因此,重力梯度测量信号中由载体垂向加速度引起的误差可表示为:
式中:为垂向加速度的正弦二倍频分量,为垂向加速度的余弦二倍频分量。经测量,机载垂向运动加速度的二倍频分量约为2×10-3g量级,在安装精度|Δβi|<5″的条件下对解调后重力梯度信号造成的影响量级为5000E,在系统动态使用时需对该项误差进行补偿。
考虑在测量过程中重力梯度仪处于高动态的飞行条件下,重力梯度敏感器中加速度计敏感轴垂向安装误差(Δβ=(Δβ1+Δβ3)-(Δβ2+Δβ4))在机载高机动过程中变化0.1″,对重力梯度信号的影响为20E量级,不可忽略,因此需对加速度计敏感轴垂向安装误差进行实时监测并采取后处理方式进行补偿。
图1为机载试验中固定翼飞机平飞运动中的垂向振动谱,可以观察到在113.2Hz处有一个0.01g的分量,为固定翼飞机螺旋桨的工作频率。可以利用此频点上的振动激励对加速度计敏感轴垂向安装误差进行实时观测并在后处理的过程中予以补偿。
一种旋转加速度计式航空重力梯度仪垂向线运动测量误差的补偿方法,如图1所示,包括以下步骤:
步骤3、将由从步骤1的加速度计组合信号得到的和信号与从步骤2的垂向加速度计输出信号得到的和信号进行相除处理,得到两组加速度计敏感轴垂向安装误差实时信号Δβs和Δβc,将两组信号Δβs和Δβc求平均得到更加真实的Δβ实时信号。其中Δβs是由飞机螺旋桨旋转频率正弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβc是由飞机螺旋桨旋转频率余弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβ是两个安装误差实时信号的平均值;
需要强调的是,本发明所述实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (1)
1.一种基于载体自振动航空重力梯度仪运动误差补偿方法,其特征在于:包括以下步骤:
步骤3、将由从步骤1的加速度计组合信号得到的和信号与从步骤2的垂向加速度计输出信号得到的和信号进行相除处理,得到两组加速度计敏感轴垂向安装误差实时信号Δβs和Δβc,将两组信号Δβs和Δβc求平均得到更加真实的Δβ实时信号;其中Δβs是由飞机螺旋桨旋转频率正弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβc是由飞机螺旋桨旋转频率余弦分量求得的加速度计敏感轴垂向安装误差实时信号,Δβ是两个安装误差实时信号的平均值;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910464038.7A CN110068876B (zh) | 2019-05-30 | 2019-05-30 | 基于载体自振动航空重力梯度仪运动误差补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910464038.7A CN110068876B (zh) | 2019-05-30 | 2019-05-30 | 基于载体自振动航空重力梯度仪运动误差补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110068876A CN110068876A (zh) | 2019-07-30 |
CN110068876B true CN110068876B (zh) | 2021-01-26 |
Family
ID=67372135
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910464038.7A Active CN110068876B (zh) | 2019-05-30 | 2019-05-30 | 基于载体自振动航空重力梯度仪运动误差补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110068876B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111650664B (zh) * | 2020-06-30 | 2022-08-26 | 东南大学 | 一种航空重力梯度仪实时重力梯度解调方法及装置 |
CN113433600B (zh) * | 2021-06-23 | 2022-04-12 | 中国船舶重工集团公司第七0七研究所 | 一种重力仪安装误差角标定方法 |
CN113447994B (zh) * | 2021-06-25 | 2022-04-12 | 中国船舶重工集团公司第七0七研究所 | 一种基于旋转加速度计原理的重力梯度全张量测量方法 |
CN113885098B (zh) * | 2021-09-07 | 2023-04-28 | 中国船舶重工集团公司第七0七研究所 | 一种重力敏感器低频频率响应误差在线建模及补偿方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2612677A1 (en) * | 2006-11-23 | 2008-05-23 | Technological Resources Pty. Limited | Compensation for unwanted accelerations in a gravity gradiometer |
DE102008054739A1 (de) * | 2008-12-16 | 2010-06-17 | Robert Bosch Gmbh | Navigationsverfahren |
CN103743395A (zh) * | 2014-01-17 | 2014-04-23 | 哈尔滨工程大学 | 一种惯性重力匹配组合导航系统中时间延迟的补偿方法 |
CN105204083A (zh) * | 2015-09-14 | 2015-12-30 | 武汉大学 | 基于绝对重力仪的垂直重力梯度提取方法及系统 |
CN108931824A (zh) * | 2018-04-27 | 2018-12-04 | 东南大学 | 一种动基座旋转加速度计重力梯度仪误差增益系数标定方法 |
CN109470241A (zh) * | 2018-11-23 | 2019-03-15 | 中国船舶重工集团公司第七0七研究所 | 一种具备重力扰动自主补偿功能的惯性导航系统及方法 |
CN109709628A (zh) * | 2019-02-15 | 2019-05-03 | 东南大学 | 一种旋转加速度计重力梯度仪标定方法 |
CN109766812A (zh) * | 2018-12-31 | 2019-05-17 | 东南大学 | 一种旋转加速度计重力梯度仪运动误差事后补偿方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AUPS114702A0 (en) * | 2002-03-18 | 2002-04-18 | Bhp Billiton Innovation Pty Ltd | Enhancement of sensors for airborne operation |
-
2019
- 2019-05-30 CN CN201910464038.7A patent/CN110068876B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2612677A1 (en) * | 2006-11-23 | 2008-05-23 | Technological Resources Pty. Limited | Compensation for unwanted accelerations in a gravity gradiometer |
DE102008054739A1 (de) * | 2008-12-16 | 2010-06-17 | Robert Bosch Gmbh | Navigationsverfahren |
CN103743395A (zh) * | 2014-01-17 | 2014-04-23 | 哈尔滨工程大学 | 一种惯性重力匹配组合导航系统中时间延迟的补偿方法 |
CN105204083A (zh) * | 2015-09-14 | 2015-12-30 | 武汉大学 | 基于绝对重力仪的垂直重力梯度提取方法及系统 |
CN108931824A (zh) * | 2018-04-27 | 2018-12-04 | 东南大学 | 一种动基座旋转加速度计重力梯度仪误差增益系数标定方法 |
CN109470241A (zh) * | 2018-11-23 | 2019-03-15 | 中国船舶重工集团公司第七0七研究所 | 一种具备重力扰动自主补偿功能的惯性导航系统及方法 |
CN109766812A (zh) * | 2018-12-31 | 2019-05-17 | 东南大学 | 一种旋转加速度计重力梯度仪运动误差事后补偿方法 |
CN109709628A (zh) * | 2019-02-15 | 2019-05-03 | 东南大学 | 一种旋转加速度计重力梯度仪标定方法 |
Non-Patent Citations (2)
Title |
---|
基于旋转加速度计原理的重力梯度测量技术研究与试验;杨晔 等;《导航定位与授时》;20170731;第4卷(第4期);第20-28页 * |
旋转加速度计重力梯度仪重力梯度信号仿真;蔡体菁 等;《物探与化探》;20151231;第39卷;第76-79页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110068876A (zh) | 2019-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110068876B (zh) | 基于载体自振动航空重力梯度仪运动误差补偿方法 | |
CN110207697B (zh) | 基于角加速度计/陀螺/加速度计的惯性导航解算方法 | |
Crawford et al. | A sensitive fast-response probe to measure turbulence and heat flux from any airplane | |
CN102023051B (zh) | 一种测量星上有效载荷三轴角位移高频微振动的方法 | |
JP5474938B2 (ja) | 加速度計および誤差補償方法 | |
CN108759845A (zh) | 一种基于低成本多传感器组合导航的优化方法 | |
JPH06288771A (ja) | 飛行機の慣性航行計器の校正方法 | |
CN201561759U (zh) | 惯性姿态方位测量装置 | |
US10654564B2 (en) | Aircraft assembly including deflection sensor | |
CN102620892A (zh) | 一种转动部件的动平衡测试方法 | |
CN103323625A (zh) | 一种mems-imu中加速度计动态环境下的误差标定补偿方法 | |
CN102322873A (zh) | 一种分布式pos地面演示验证系统 | |
CN109766812A (zh) | 一种旋转加速度计重力梯度仪运动误差事后补偿方法 | |
JP6629691B2 (ja) | センサパッケージおよび自動運転車両 | |
CN111650664B (zh) | 一种航空重力梯度仪实时重力梯度解调方法及装置 | |
Rodi et al. | Correction of static pressure on a research aircraft in accelerated flight using differential pressure measurements | |
CN102072729A (zh) | 用于测量旋转飞行器姿态的测量装置及其测量方法 | |
CN106885587A (zh) | 旋翼扰动下惯性/gps组合导航外杆臂效应误差补偿方法 | |
KR101658473B1 (ko) | Mems자이로스코프의 가속도 민감도 보정 방법 | |
US6810739B1 (en) | Accelerometer augmented leveling device | |
CN109931959B (zh) | 硅微陀螺仪正交误差校正方法 | |
US7461548B2 (en) | Method and a device for measuring the speed of an aircraft, in particular a rotorcraft at low speed | |
Pérez-D'Arpino et al. | Development of a low cost inertial measurement unit for uav applications with kalman filter based attitude determination | |
CN114370930A (zh) | 利用mhd微角振动传感器测量卫星结构微角振动的方法 | |
EA037452B1 (ru) | Гравиметрический узел |
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 |