CN105675018B - 一种姿态检测系统标定方法 - Google Patents
一种姿态检测系统标定方法 Download PDFInfo
- Publication number
- CN105675018B CN105675018B CN201610018024.9A CN201610018024A CN105675018B CN 105675018 B CN105675018 B CN 105675018B CN 201610018024 A CN201610018024 A CN 201610018024A CN 105675018 B CN105675018 B CN 105675018B
- Authority
- CN
- China
- Prior art keywords
- pulse
- axis
- angular speed
- gyroscope
- acceleration transducer
- 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.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000001133 acceleration Effects 0.000 claims abstract description 35
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 claims description 18
- 238000009434 installation Methods 0.000 claims description 4
- 210000001367 artery Anatomy 0.000 claims 1
- 210000003462 vein Anatomy 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Manufacturing & Machinery (AREA)
- Automation & Control Theory (AREA)
- Gyroscopes (AREA)
Abstract
本发明公开了一种姿态检测系统标定方法,属于管道姿态检测技术领域。本发明中,根据姿态检测系统中陀螺仪与加速度传感器误差特性,建立误差模型,在通过位置标定和速率标定来标定上述误差模型中的误差系数,该标校方法简单易用,可操作性强,适于工程应用。
Description
技术领域
本发明涉及管道姿态检测技术领域,尤其涉及一种姿态检测系统标定方法。
背景技术
随着石油天然气等能源的需求量逐渐加大,管道作为最基础、最经济的能源运输方式之一,其安全性越来越受到重视。姿态检测系统随管道检测器在管道内运行,为管道缺陷检测提供沿管道周向旋转的角度和倾角信息。姿态检测系统通常包括陀螺仪和加速度传感器,姿态检测系统的工作性能受许多误差源如安装误差、制造误差等的影响,因此,有必要标定姿态检测系统的各误差参数,并在组成姿态检测系统时对误差做补偿。
发明内容
本发明的目的是提供一种简单易用、操作性强的姿态检测系统标定方法,以修正姿态检测系统的常值误差,提高姿态检测系统的精度。
为达上述目的,本发明提出如下技术方案:
一种姿态检测系统标定方法,包括以下步骤:
步骤一、根据姿态检测系统中陀螺仪与加速度传感器误差特性,建立误差模型:Ngrx=Dfx+Sgrxωx
Ngx=Kgx0+Kgx1ay+Kgx2az+Sgxax (1);
Ngy=Kgy0+Kgy1ax+Kgy2az+Sgyay
Ngz=Kgz0+Kgz1ax+Kgz2ay+Sgzaz
式中,Ngrx为陀螺仪的脉冲输出,脉冲数/s;Ngx、Ngy、Ngz分别为加速度传感器在X轴、Y轴、Z轴的脉冲输出,脉冲数/s;Dfx为陀螺仪的常值漂移,脉冲数/s;Sgrx为陀螺仪的标度因数,(脉冲数/s)/(°/s);Kgx1、Kgx2、Kgy1、Kgy2、Kgz1,Kgz2为加速度传感器的安装误差系数,(脉冲数/S)/(°/s);Kgx0、Kgy0、Kgz0为加速度传感器的偏值,脉冲数/S;Sgx、Sgy、Sgz为加速度传感器的标度因数,(脉冲数/S)/g;
步骤二、标定,包括位置标定和速率标定,位置标定设计6个位置对加速度传感器的误差系数进行标定,6个标定位置分别为令转台的X轴、Y轴、Z轴分别与地球坐标系的天、地重合的六个位置,让姿态检测系统依次处于6个位置,调整好转台后,连续记录一定时间T1内加速度传感器的脉冲输出;
速率标定,设计至少4个角速度对陀螺仪的误差系数进行标定,采集陀螺仪在4个角速度下且在一定时间T2内的脉冲输出,其中两个角速度分别为同一速率下的沿顺时针和逆时针的角速度,另两个角速度分别为另一相同速率下的沿顺时针和逆时针的角速度;
步骤三、解算
首先,利用步骤二的位置标定中获得的数据对加速度传感器误差系数进行解算,再利用下列公式(5)和(6)解算陀螺仪的误差系数:
式中,Ngx(2)、Ngx(3)、Ngx(4)、Ngx(5)为步骤二的位置标定中的除转台的X轴与地球坐标系的地、天重合外的4个位置下加速度计X轴脉冲输出,Ngrx+(1')、Ngrx-(2')、Ngrx+(3')、Ngrx-(4')为陀螺仪在4个角速度下的脉冲输出,且其中Ngrx+(1')、Ngrx-(2')为陀螺仪在同一速率的分别沿顺时针和逆时针的角速度下的的脉冲输出,Ngrx+(3')、Ngrx-(4')为陀螺仪在同一速率的分别沿顺时针和逆时针的角速度下的的脉冲输出。
更进一步地,步骤三中通过下列公式(2)、(3)、(4)对加速度传感器误差系数进行解算:
式中,Ngx(1)、Ngx(2)、Ngx(3)、Ngx(4)、Ngx(5)、Ngx(6)为步骤二的位置标定中的6个位置下加速度计X轴脉冲输出,Ngy(1)、Ngy(2)、Ngy(3)、Ngy(4)、Ngy(5)、Ngy(6)为步骤二的位置标定中的6个位置下加速度计Y轴脉冲输出,Ngz(1)、Ngz(2)、Ngz(3)、Ngz(4)、Ngz(5)、Ngz(6)为步骤二的位置标定中的6个位置下加速度计Z轴脉冲输出,Kgx0、Kgx1、Kgx2,Kgy0、Kgy1、Kgy2,Kgz0、Kgz1、Kgz2分别为加速度计X轴、Y轴、Z轴的误差系数。
更进一步地,步骤二中通过以下四组实验进行速率标定:实验1’、角速度60°/s,顺时针;实验2’、角速度-60°/s,逆时针;实验3’、角速度120°/s,顺时针;实验4’、角速度-120°/s,逆时针;另外,四组实验的T2均为12s。
更进一步地,该方法还包括步骤四、补偿,通过速率标定解算的Sgrx和Dfx,位置标定解算的Sgx,Sgy,Sgz,Kgx0,Kgx1,Kgx2,Kgy0,Kgy1,Kgy2,Kgz0,Kgz1,Kgz2,利用公式(1)补偿陀螺仪与加速度传感器的测量值。
本发明的有益效果为:
本发明中,根据姿态检测系统中陀螺仪与加速度传感器误差特性,建立误差模型,在通过位置标定和速率标定来标定上述误差模型中的误差系数,该标校方法简单易用,可操作性强,适于工程应用。
具体实施方式
下面通过具体实施方式来进一步说明本发明的技术方案。
实施例一:
本实施例提供了一种姿态检测系统标定方法,姿态检测系统包括陀螺仪与三轴加速度传感器,该标定方法具体包括以下步骤:
步骤一、根据姿态检测系统中陀螺仪与加速度传感器误差特性,建立误差模型:Ngrx=Dfx+Sgrxωx
Ngx=Kgx0+Kgx1ay+Kgx2az+Sgxax (1);
Ngy=Kgy0+Kgy1ax+Kgy2az+Sgyay
Ngz=Kgz0+Kgz1ax+Kgz2ay+Sgzaz
式中,Ngrx为陀螺仪的脉冲输出,脉冲数/s;Ngx、Ngy、Ngz分别为加速度传感器在X轴、Y轴、Z轴的脉冲输出,脉冲数/s;Dfx为陀螺仪的常值漂移,脉冲数/s;Sgrx为陀螺仪的标度因数,(脉冲数/s)/(°/s);Kgx1、Kgx2、Kgy1、Kgy2、Kgz1,Kgz2为加速度传感器的安装误差系数,(脉冲数/S)/(°/s);Kgx0、Kgy0、Kgz0为加速度传感器的偏值,脉冲数/S;Sgx、Sgy、Sgz为加速度传感器的标度因数,(脉冲数/S)/g;
步骤二、标定,包括位置标定和速率标定,位置标定设计6个位置对加速度传感器的误差系数进行标定,6个标定位置分别为令转台的X轴、Y轴、Z轴分别与地球坐标系的天、地重合的6个位置,具体地,这6个位置如表1所示,让姿态检测系统依次处于6个位置,调整好转台后,连续记录一定时间T1内加速度传感器的脉冲输出,本实施例中,T1为60s;
表1
位置 | X轴 | Y轴 | Z轴 | ax | ay | az | ωx |
1 | 地 | 北 | 东 | -g | 0 | 0 | -ωieu |
2 | 东 | 天 | 南 | 0 | g | 0 | 0 |
3 | 东 | 北 | 天 | 0 | 0 | g | 0 |
4 | 东 | 地 | 北 | 0 | -g | 0 | 0 |
5 | 东 | 南 | 地 | 0 | 0 | -g | 0 |
6 | 天 | 西 | 南 | g | 0 | 0 | ωieu |
速率标定,设计至少四组实验1’、2’、3’、4’对陀螺仪的误差系数进行标定,四组实验采集陀螺仪在4个角速度下且在一定时间T2内的脉冲输出,其中两个角速度分别为同一速率下的沿顺时针和逆时针的角速度,另两个角速度分别为另一相同速率下的沿顺时针和逆时针的角速度,具体地,该四组实验中的角速度、采样时间和角速度的方向如表2所示:
表2
序号 | 角速度 | T2 | |
1’ | 60°/s | 12s | 顺时针 |
2’ | -60°/s | 12s | 逆时针 |
3’ | 120°/s | 12s | 顺时针 |
4’ | -120°/s | 12s | 逆时针 |
步骤三、解算
首先,利用步骤二的位置标定中获得的数据对加速度传感器误差系数进行解算,具体地,通过下列公式(2)、(3)、(4)对加速度传感器误差系数进行解算:
式中,Ngx(1)、Ngx(2)、Ngx(3)、Ngx(4)、Ngx(5)、Ngx(6)为步骤二的位置标定中的6个位置下加速度计X轴脉冲输出,Ngy(1)、Ngy(2)、Ngy(3)、Ngy(4)、Ngy(5)、Ngy(6)为步骤二的位置标定中的6个位置下加速度计Y轴脉冲输出,Ngz(1)、Ngz(2)、Ngz(3)、Ngz(4)、Ngz(5)、Ngz(6)为步骤二的位置标定中的6个位置下加速度计Z轴脉冲输出,Kgx0、Kgx1、Kgx2,Kgy0、Kgy1、Kgy2,Kgz0、Kgz1、Kgz2分别为加速度计X轴、Y轴、Z轴的误差系数;
其次,再利用下列公式(5)和(6)解算陀螺仪的误差系数:
步骤四、补偿,通过速率标定解算的Sgrx和Dfx,位置标定解算的Sgx,Sgy,Sgz,Kgx0,Kgx1,Kgx2,Kgy0,Kgy1,Kgy2,Kgz0,Kgz1,Kgz2,利用公式(1)补偿陀螺仪与加速度传感器的测量值。
本实施例中,充分利用姿态系统内部传感器信息,依据陀螺误差在短时间内变化不大的特点,首先对加速度传感器角度误差进行修正,在后续检测工作中,以加速度传感器角度为参考,修正陀螺随时间漂移的误差,使信息互用,达到姿态角实时修正的目的,该标校方法简单易用,可操作性强,适于工程应用。
以上结合具体实施例描述了本发明的技术原理。这些描述只是为了解释本发明的原理,而不能以任何方式解释为对本发明保护范围的限制。基于此处的解释,本领域的技术人员不需要付出创造性的劳动即可联想到本发明的其它具体实施方式,这些方式都将落入本发明的保护范围之内。
Claims (4)
1.一种姿态检测系统标定方法,其特征在于,包括以下步骤:
步骤一、根据姿态检测系统中陀螺仪与加速度传感器误差特性,建立误差模型:
式中,Ngrx为陀螺仪的脉冲输出,脉冲数/s;Ngx、Ngy、Ngz分别为加速度传感器在X轴、Y轴、Z轴的脉冲输出,脉冲数/s;Dfx为陀螺仪的常值漂移,脉冲数/s;Sgrx为陀螺仪的标度因数,脉冲数/°;Kgx1、Kgx2、Kgy1、Kgy2、Kgz1,Kgz2为加速度传感器的安装误差系数,脉冲数/°;Kgx0、Kgy0、Kgz0为加速度传感器的偏值,脉冲数/S;Sgx、Sgy、Sgz为加速度传感器的标度因数,脉冲数/S/g;
步骤二、标定,包括位置标定和速率标定,位置标定设计6个位置对加速度传感器的误差系数进行标定,6个标定位置分别为令转台的X轴、Y轴、Z轴分别与地球坐标系的天、地重合的六个位置,让姿态检测系统依次处于6个位置,调整好转台后,连续记录一定时间T1内加速度传感器的脉冲输出;
速率标定,设计至少4个角速度对陀螺仪的误差系数进行标定,采集陀螺仪在4个角速度下且在一定时间T2内的脉冲输出,其中两个角速度分别为同一速率下的沿顺时针和逆时针的角速度,另两个角速度分别为另一相同速率下的沿顺时针和逆时针的角速度;
步骤三、解算
首先,利用步骤二的位置标定中获得的数据对加速度传感器误差系数进行解算,再利用下列公式(5)和(6)解算陀螺仪的误差系数:
式中,Ngx(2)、Ngx(3)、Ngx(4)、Ngx(5)为步骤二的位置标定中的除转台的X轴与地球坐标系的地、天重合外的4个位置下加速度计X轴脉冲输出,Ngrx+(1')、Ngrx-(2')、Ngrx+(3')、Ngrx-(4')为陀螺仪在4个角速度下的脉冲输出,且其中Ngrx+(1')、Ngrx-(2')为陀螺仪在同一速率的分别沿顺时针和逆时针的角速度下的的脉冲输出,Ngrx+(3')、Ngrx-(4')为陀螺仪在同一速率的分别沿顺时针和逆时针的角速度下的的脉冲输出。
2.一种根据权利要求1所述的姿态检测系统标定方法,其特征在于,步骤三中通过下列公式(2)、(3)、(4)对加速度传感器误差系数进行解算:
式中,Ngx(1)、Ngx(2)、Ngx(3)、Ngx(4)、Ngx(5)、Ngx(6)为步骤二的位置标定中的6个位置下加速度计X轴脉冲输出,Ngy(1)、Ngy(2)、Ngy(3)、Ngy(4)、Ngy(5)、Ngy(6)为步骤二的位置标定中的6个位置下加速度计Y轴脉冲输出,Ngz(1)、Ngz(2)、Ngz(3)、Ngz(4)、Ngz(5)、Ngz(6)为步骤二的位置标定中的6个位置下加速度计Z轴脉冲输出,Kgx0、Kgx1、Kgx2,Kgy0、Kgy1、Kgy2,Kgz0、Kgz1、Kgz2分别为加速度计X轴、Y轴、Z轴的误差系数。
3.一种根据权利要求1所述的姿态检测系统标定方法,其特征在于,步骤二中通过以下四组实验进行速率标定:实验1’、角速度60°/s,顺时针;实验2’、角速度-60°/s,逆时针;实验3’、角速度120°/s,顺时针;实验4’、角速度-120°/s,逆时针;另外,四组实验的T2均为12s。
4.一种根据权利要求1至3任一所述的姿态检测系统标定方法,其特征在于,还包括步骤四、补偿,通过速率标定解算的Sgrx和Dfx,位置标定解算的Sgx,Sgy,Sgz,Kgx0,Kgx1,Kgx2,Kgy0,Kgy1,Kgy2,Kgz0,Kgz1,Kgz2,利用公式(1)补偿陀螺仪与加速度传感器的测量值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610018024.9A CN105675018B (zh) | 2016-01-12 | 2016-01-12 | 一种姿态检测系统标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610018024.9A CN105675018B (zh) | 2016-01-12 | 2016-01-12 | 一种姿态检测系统标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105675018A CN105675018A (zh) | 2016-06-15 |
CN105675018B true CN105675018B (zh) | 2018-07-13 |
Family
ID=56300240
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610018024.9A Expired - Fee Related CN105675018B (zh) | 2016-01-12 | 2016-01-12 | 一种姿态检测系统标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105675018B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115752508B (zh) * | 2022-11-14 | 2024-02-06 | 北京自动化控制设备研究所 | 微机电惯性系统高阶耦合误差补偿模型及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1818555A (zh) * | 2006-03-29 | 2006-08-16 | 北京航空航天大学 | 安装误差角与标度因数解耦的微惯性测量单元精确标定方法 |
CN102589573A (zh) * | 2012-02-09 | 2012-07-18 | 黑龙江省博凯科技开发有限公司 | 微型组合导航系统中的传感器野外标定方法 |
CN102853850A (zh) * | 2012-09-11 | 2013-01-02 | 中国兵器工业集团第二一四研究所苏州研发中心 | 基于单轴转台的三轴mems陀螺旋转积分标定方法 |
CN103196462A (zh) * | 2013-02-28 | 2013-07-10 | 南京航空航天大学 | 一种mimu中mems陀螺仪的误差标定补偿方法 |
CN103954299A (zh) * | 2014-04-22 | 2014-07-30 | 北京航天控制仪器研究所 | 一种标定捷联惯性组合陀螺仪组合的方法 |
CN104596543A (zh) * | 2015-01-13 | 2015-05-06 | 北京航天控制仪器研究所 | 一种基准不确定情况下陀螺仪组合误差系数标定方法 |
-
2016
- 2016-01-12 CN CN201610018024.9A patent/CN105675018B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1818555A (zh) * | 2006-03-29 | 2006-08-16 | 北京航空航天大学 | 安装误差角与标度因数解耦的微惯性测量单元精确标定方法 |
CN102589573A (zh) * | 2012-02-09 | 2012-07-18 | 黑龙江省博凯科技开发有限公司 | 微型组合导航系统中的传感器野外标定方法 |
CN102853850A (zh) * | 2012-09-11 | 2013-01-02 | 中国兵器工业集团第二一四研究所苏州研发中心 | 基于单轴转台的三轴mems陀螺旋转积分标定方法 |
CN103196462A (zh) * | 2013-02-28 | 2013-07-10 | 南京航空航天大学 | 一种mimu中mems陀螺仪的误差标定补偿方法 |
CN103954299A (zh) * | 2014-04-22 | 2014-07-30 | 北京航天控制仪器研究所 | 一种标定捷联惯性组合陀螺仪组合的方法 |
CN104596543A (zh) * | 2015-01-13 | 2015-05-06 | 北京航天控制仪器研究所 | 一种基准不确定情况下陀螺仪组合误差系数标定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105675018A (zh) | 2016-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103591949B (zh) | 三轴姿态测量系统非正交性误差的正交补偿方法 | |
CN110160554B (zh) | 一种基于寻优法的单轴旋转捷联惯导系统标定方法 | |
CN103954286B (zh) | 微小卫星磁传感器多误差模型在轨迭代标定方法 | |
CN104344837B (zh) | 一种基于速度观测的冗余惯导系统加速度计系统级标定方法 | |
CN104698485B (zh) | 基于bd、gps及mems的组合导航系统及导航方法 | |
CN107270891B (zh) | 基于抗差估计的惯性地磁匹配定位方法 | |
Li et al. | Analysis and improvement of attitude output accuracy in rotation inertial navigation system | |
CN102788578A (zh) | 一种基于局部重力场逼近的匹配导航方法 | |
CN106017452B (zh) | 双陀螺抗扰动寻北方法 | |
CN102889076A (zh) | 陀螺测斜仪标定方法 | |
CN103954303B (zh) | 一种用于磁力计导航系统航向角动态计算及校正方法 | |
CN109612460B (zh) | 一种基于静止修正的垂线偏差测量方法 | |
CN104406610A (zh) | 一种磁力计实时校正装置及方法 | |
CN104166169A (zh) | 一种捷联于水下潜器的三轴磁场传感器误差补偿方法 | |
CN108507592A (zh) | 一种双轴旋转惯性导航系统转轴非正交角标定方法 | |
CN106840201B (zh) | 一种带双轴转位机构捷联惯导的三位置自对准方法 | |
CN103363989B (zh) | 一种捷联惯导系统内杆臂的估计与误差补偿方法 | |
CN106767925A (zh) | 带双轴转位机构的惯导系统三位置参数辨识对准方法 | |
Bistrov | Performance analysis of alignment process of MEMS IMU | |
Dai et al. | In-field calibration method for DTG IMU including g-sensitivity biases | |
Yang et al. | A robust Mag/INS-based orientation estimation algorithm for measurement while drilling | |
CN110736484B (zh) | 基于陀螺仪及磁传感器融合的背景磁场标定方法 | |
Lin et al. | A high-accuracy method for calibration of nonorthogonal angles in dual-axis rotational inertial navigation system | |
Zhang et al. | Moving state marine SINS initial alignment based on high degree CKF | |
CN113418499B (zh) | 一种旋转飞行器滚转角解算方法及系统 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180713 Termination date: 20190112 |