CN111273058A - 一种加速度计标定方法 - Google Patents
一种加速度计标定方法 Download PDFInfo
- Publication number
- CN111273058A CN111273058A CN202010264315.2A CN202010264315A CN111273058A CN 111273058 A CN111273058 A CN 111273058A CN 202010264315 A CN202010264315 A CN 202010264315A CN 111273058 A CN111273058 A CN 111273058A
- Authority
- CN
- China
- Prior art keywords
- axis
- accelerometer
- vertical plane
- measured
- axes
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明实施例公开了一种加速度计标定方法,获取加速度计在三轴旋转台上的静态数学模型,并将加速度计被测轴上的重力加速度分量作为参数的输入;利用在静态数学模型位于三轴旋转台不同垂直平面的位置时加速度计被测轴的参数输出结果,由谐波系数法解算出静态数学模型的模型参数,获得标定的静态数学模型,解决加速度计在没有测量加速度计与转台安装偏差的光学器件以及没有北向测量的条件下加速度计标定问题,无需测量加速度计与转台的安装偏差,无需测量指北方向,使得加速度计标定模型的获取更加简单快捷。
Description
技术领域
本发明实施例涉及加速度计标定技术领域,具体涉及一种加速度计标定方法。
背景技术
所谓传感器的标定,是指通过试验监理传感器输出与输入之间的关系并确定不同使用条件下的误差的过程,一般来说,对传感器进行标定时,必须以国家和地方计量部分的有关检定规程为依据,选择正确的标定条件和适当的仪器设备,按照一定的程序进行。
传感器标定的基本方法是将已知的被测量作为待标定传感器的输入,同时用输出量测量环节将待标定传感器的输出信号测量并显示出来,对所获得的传感器输入量和输出量进行处理和比较,从而获得一系列表征两者对应关系的标定曲线,进而获得传感器性能指标的实测结果。
而通常在实验中,没有测量加速度计与转台安装偏差的光学器件以及没有北向测量条件下的条件下,则无法对加速度计进行标定。
发明内容
为此,本发明实施例提供一种加速度计标定方法,解决在没有测量加速度计与转台安装偏差的光学器件以及没有北向测量条件下的加速度计标定的问题。
为了实现上述目的,本发明的实施方式提供如下技术方案:
一种加速度计标定方法,包括如下步骤:
S100、获取加速度计在三轴旋转台上的静态数学模型,并将加速度计被测轴上的重力加速度分量作为参数的输入;
S200、利用在静态数学模型位于三轴旋转台不同垂直平面的位置时加速度计被测轴的参数输出结果,由谐波系数法解算出静态数学模型的模型参数,获得标定的静态数学模型。
作为本发明的一种优选方案,所述的加速度计被测轴具体包括X、Y和Z轴,并分别计算X、Y和Z轴分别位于垂直于三轴旋转平台的垂直面内的模型参数。
作为本发明的一种优选方案,所述X、Y和Z轴分别位于垂直于三轴旋转平台的垂直面的具体摆放情况有六种:
X轴位于垂直面内:
X、Y轴位于垂直面内和X、Z轴位于垂直面内;
Y轴位于垂直面内:
Y、X轴位于垂直面内和Y、Z轴位于垂直面内;
Z轴位于垂直面内:
Z、X轴位于垂直面内和Z、Y位于垂直面内。
作为本发明的一种优选方案,以X轴为例:
式中:Ax为实际输出的加速度值;ax、ay和az分别为X、Y和Z轴方向的加速度值分量;K0x为零偏;K1x为标度因数;K2x为二次非线性系数;K3x为三次的非线性系数;Kxy、Kxz分别为x轴与y轴和x轴与z轴间的交叉耦合系数;Ky、Kz分别为x轴相对x轴绕y轴和z轴的失准角;
当X、Y轴位于平面内,并加入加速度计在三轴平台上的安装误差角,加速度计被测轴上的重力加速度分量作为参数的输入的加速度为:
式中:θ为输入轴与水平面的夹角;θ0为初始安装误差角;则Ax为:
作为本发明的一种优选方案,利用谐波系数法计算谐波系数与静态数学模型的方程系数之间的关系:
作为本发明的一种优选方案,在计算X、Y和Z轴分别位于垂直于三轴旋转平台的垂直面内的模型参数,垂直面的圆周上均匀分布至少8个测试点,并通过转动至少3圈获取24个位置的加速度计被测轴的参数输出结果。
作为本发明的一种优选方案,在垂直面的圆周上均匀分布的测试点与三轴平台的水平轴的夹角为:
并同时提取加速度计被测轴在两个旋转的垂直平面上的平均值,记录对应角度。
作为本发明的一种优选方案,在获得静态数学模型位于三轴旋转台不同垂直平面的位置时加速度计被测轴的参数输出结果后,利用最小二乘法对输出结果数据进行处理。
作为本发明的一种优选方案,在以X轴为例的计算过程中,由于z轴朝上,所以在输出数据中应该减去z轴上重力加速度的输出,并将采集的数据中除x,y,z轴的输出以外的量除去。
作为本发明的一种优选方案,将加速度计安装在三轴转台上,并使得加速度计各个被测轴与三轴转台的坐标轴平行,并采集至少2分钟的垂直面的圆周上均匀分布的测试点的位置数据。
本发明的实施方式具有如下优点:
本发明由谐波系数法解算出静态数学模型的模型参数,获得标定的静态数学模型,解决加速度计在没有测量加速度计与转台安装偏差的光学器件以及没有北向测量的条件下加速度计标定问题,无需测量加速度计与转台的安装偏差,无需测量指北方向,使得加速度计标定模型的获取更加简单快捷。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引伸获得其它的实施附图。
本说明书所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明可实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容得能涵盖的范围内。
图1为本发明实施方式中一种加速度计标定方法流程图;
图2为本发明实施方式中x,y轴位于垂直平面时的加速度输入示意图;
图3为本发明实施方式中加速度计X、Y和Z轴分别位于垂直于三轴旋转平台的垂直面的具体摆放情况六种位置示意图。
具体实施方式
以下由特定的具体实施例说明本发明的实施方式,熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本发明的其他优点及功效,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供了一种加速度计标定方法,包括加速度计建模与标定,具体包括
S100、获取加速度计在三轴旋转台上的静态数学模型,并将加速度计被测轴上的重力加速度分量作为参数的输入;
S200、利用在静态数学模型位于三轴旋转台不同垂直平面的位置时加速度计被测轴的参数输出结果,由谐波系数法解算出静态数学模型的模型参数,获得标定的静态数学模型。
实验设备为三轴速率转台,没有测量加速度计与转台安装偏差的光学器件;没有北向,所以无法计算出地球自转对速率实验的影响。
于是,根据GJB719-89《单轴摆式伺服线加速度计通用规范》的加速度计测量模型,结合现有实验条件,忽略不可测因素以及该型号不需考虑的因素。
进一步由GJB719-89《单轴摆式伺服线加速度计通用规范》获取加速度计的静态数学模型。
以X轴为例:
式中:Ax为实际输出的加速度值;ax、ay和az分别为X、Y和Z轴方向的加速度值分量;K0x为零偏;K1x为标度因数;K2x为二次非线性系数;K3x为三次的非线性系数;Kxy、Kxz分别为x轴与y轴和x轴与z轴间的交叉耦合系数;Ky、Kz分别为x轴相对x轴绕y轴和z轴的失准角;
为了标定上述模型系数,垂直平面内不同位置时,重力加速度在被测轴上的加速度分量作为参考的输入,采用多位置的方法标定模型中参数。
所述的加速度计被测轴具体包括X、Y和Z轴,并分别计算X、Y和Z轴分别位于垂直于三轴旋转平台的垂直面内的模型参数。
如图2和图3所示,所述X、Y和Z轴分别位于垂直于三轴旋转平台的垂直面的具体摆放情况有六种:
X轴位于垂直面内:
X、Y轴位于垂直面内和X、Z轴位于垂直面内;
Y轴位于垂直面内:
Y、X轴位于垂直面内和Y、Z轴位于垂直面内;
Z轴位于垂直面内:
Z、X轴位于垂直面内和Z、Y位于垂直面内。
被测轴x必须位于垂直平面内,而三轴加速度计使得x轴在垂直平面内的摆放情况有两类,即x、y轴位于垂直平面内和x、z轴位于垂直平面内。当x、y轴位于垂直平面内时,z轴位于水平面内且加速度值为0,此时与z轴相关的误差模型参数无法标定出。
所以,加速度计的标定方案中必须包括上述两类摆放情况,分别旋转至某几个位置,进行位置实验。
为了便于数据分析,对上式进行归一化处理,即加速度单位采用g。
当X、Y轴位于平面内,并加入加速度计在三轴平台上的安装误差角,加速度计被测轴上的重力加速度分量作为参数的输入的加速度为:
式中:θ为输入轴与水平面的夹角;θ0为初始安装误差角;若忽略二阶及其以上小量,则Ax为:
谐波系数法求解模型参数
为了利用不同垂直平面的位置时加速度计被测轴的输出,求解模型参数,可采用谐波系数法,简化求解上述方程。
首先将上式改写成:
Ax=B0+S1 sinθ+C1 cosθ+S2 sin2θ+C2 cos2θ+S3 sin3θ
比较各次谐波系数与模型方程系数之间的关系为:
谐波系数一共有六个谐波系数,所以x,y轴在垂直平面内时和x,z轴在垂直平面内时分别至少有六个位置的测量才能标定出系数。由于标定实验一般都采用对称测量,以消除一些干扰因素的影响,所以确定每个垂直平面使用八位置法标定参数。
为了标定出加速度计静态模型,必须在以上方程中选择忽略一项。K1x不是小量,不能忽略;实验的安装偏差可能会较大,无法达到角秒级,安装偏差不能忽略。Ky是小量,选择忽略Ky,于是可以解出其他参数。当x,z轴位于垂直平面内时,同理可以解出Kxz,则模型中所有参数即可确定。
图3中1和6,2和3,4和5分别是相同的旋转平面,所以在实际实验中,为了降低实验成本,一共需要旋转3圈即可,每圈8个实验位置,一共24位置。
所以,三轴加速度计标定方法为采用24位置法,即可得出加速度计的模型系数。
虽然在同一旋转平面上,但是不同的轴,对应不同的角度输入。数据分析时,需要同时提取在旋转平面上的两个轴平均值,并记录对应角度。
加速度计标定步骤:
①安装:将加速度计安装在三轴转台上。使得加速度计各个坐标轴平行,X轴——内环轴,Y轴——中环轴,Z轴——外环轴;
②调平:X、Y轴在水平面上,Z轴垂直指天;
③上电、稳定:连接好数据线,板卡,电池,将产品稳定三分钟后,开始实验;
④采集此位置数据2分钟,x,y,z轴的输出,此处是测试的x,y轴。加速度计模型参数系数解算。
由于无离心机设备条件,加速度计标定过程中,忽略动态模型的影响,只标定静态参数,采用最小二乘法对输出数据进行处理。
最小二乘法是求得模型参数的一类方法,它是以模型与样本数据的误差的平方和最小为标准。
解算结果是所有测量值的共同作用,并且保证了结果的一致性。
根据模型方程可以解算出各参数,得到加速度计三轴测量模型。
本发明还提供了一种标定结果分析测试方法:
将产品安装在三轴转台上,z轴朝上放置在调平的平台上,连接好采集PC。打开直流稳压电源,待产品稳定后,开始记录产品的输出数据,记录时间为2小时。
数据预处理:由于z轴朝上,所以在输出数据中应该减去z轴上重力加速度的输出,并将采集的数据中除x,y,z轴的输出以外的量除去。零偏稳定性的计算,实为计算输出数据的标准偏差。
虽然,上文中已经用一般性说明及具体实施例对本发明作了详尽的描述,但在本发明基础上,可以对之作一些修改或改进,这对本领域技术人员而言是显而易见的。因此,在不偏离本发明精神的基础上所做的这些修改或改进,均属于本发明要求保护的范围。
Claims (10)
1.一种加速度计标定方法,其特征在于,包括如下步骤:
S100、获取加速度计在三轴旋转台上的静态数学模型,并将加速度计被测轴上的重力加速度分量作为参数的输入;
S200、利用在静态数学模型位于三轴旋转台不同垂直平面的位置时加速度计被测轴的参数输出结果,由谐波系数法解算出静态数学模型的模型参数,获得标定的静态数学模型。
2.根据权利要求1所述的一种加速度计标定方法,其特征在于,所述的加速度计被测轴具体包括X、Y和Z轴,并分别计算X、Y和Z轴分别位于垂直于三轴旋转平台的垂直面内的模型参数。
3.根据权利要求2所述的一种加速度计标定方法,其特征在于,所述X、Y和Z轴分别位于垂直于三轴旋转平台的垂直面的具体摆放情况有六种:
X轴位于垂直面内:
X、Y轴位于垂直面内和X、Z轴位于垂直面内;
Y轴位于垂直面内:
Y、X轴位于垂直面内和Y、Z轴位于垂直面内;
Z轴位于垂直面内:
Z、X轴位于垂直面内和Z、Y位于垂直面内。
6.根据权利要求2或5所述的一种加速度计标定方法,其特征在于,在计算X、Y和Z轴分别位于垂直于三轴旋转平台的垂直面内的模型参数,垂直面的圆周上均匀分布至少8个测试点,并通过转动至少3圈获取24个位置的加速度计被测轴的参数输出结果。
8.根据权利要求2所述的一种加速度计标定方法,其特征在于,在获得静态数学模型位于三轴旋转台不同垂直平面的位置时加速度计被测轴的参数输出结果后,利用最小二乘法对输出结果数据进行处理。
9.根据权利要求6所述的一种加速度计标定方法,其特征在于,在以X轴为例的计算过程中,由于z轴朝上,所以在输出数据中应该减去z轴上重力加速度的输出,并将采集的数据中除x,y,z轴的输出以外的量除去。
10.根据权利要求6所述的一种加速度计标定方法,其特征在于,将加速度计安装在三轴转台上,并使得加速度计各个被测轴与三轴转台的坐标轴平行,并采集至少2分钟的垂直面的圆周上均匀分布的测试点的位置数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010264315.2A CN111273058A (zh) | 2020-04-07 | 2020-04-07 | 一种加速度计标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010264315.2A CN111273058A (zh) | 2020-04-07 | 2020-04-07 | 一种加速度计标定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111273058A true CN111273058A (zh) | 2020-06-12 |
Family
ID=70998235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010264315.2A Pending CN111273058A (zh) | 2020-04-07 | 2020-04-07 | 一种加速度计标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111273058A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113376404A (zh) * | 2021-05-19 | 2021-09-10 | 中寰卫星导航通信有限公司 | 一种标定验证系统、设备、方法及存储介质 |
CN114324978A (zh) * | 2021-12-17 | 2022-04-12 | 兰州空间技术物理研究所 | 一种加速度计捕获范围的地面静态标定方法 |
CN114397479A (zh) * | 2021-12-07 | 2022-04-26 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种石英振梁加速度计to部件的两自由度测试装置及方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW528854B (en) * | 2001-09-04 | 2003-04-21 | American Gnc Corp | Digital signal processing method and system thereof for precision orientation measurements |
US6640609B1 (en) * | 2001-07-27 | 2003-11-04 | Trimble Navigation Limited | Spectral method for calibrating accelerometers |
CN1818555A (zh) * | 2006-03-29 | 2006-08-16 | 北京航空航天大学 | 安装误差角与标度因数解耦的微惯性测量单元精确标定方法 |
CN112166685B (zh) * | 2008-09-12 | 2012-10-31 | 中国电子科技集团公司第二十六研究所 | 加速度计组合及其控制方法 |
CN102889076A (zh) * | 2012-09-14 | 2013-01-23 | 西安思坦仪器股份有限公司 | 陀螺测斜仪标定方法 |
CN104133081A (zh) * | 2014-07-31 | 2014-11-05 | 工业和信息化部电子第五研究所 | 加速度计性能参数的检测方法和系统 |
CN106813680A (zh) * | 2016-12-28 | 2017-06-09 | 兰州空间技术物理研究所 | 一种高精度、高分辨率石英扰性传感器的静态标定方法 |
CN109813336A (zh) * | 2017-11-22 | 2019-05-28 | 广东虚拟现实科技有限公司 | 惯性测量单元标定方法 |
CN110006450A (zh) * | 2019-04-15 | 2019-07-12 | 哈尔滨工业大学 | 一种激光捷联惯导系统在卧式三轴转台上的标定方法 |
-
2020
- 2020-04-07 CN CN202010264315.2A patent/CN111273058A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6640609B1 (en) * | 2001-07-27 | 2003-11-04 | Trimble Navigation Limited | Spectral method for calibrating accelerometers |
TW528854B (en) * | 2001-09-04 | 2003-04-21 | American Gnc Corp | Digital signal processing method and system thereof for precision orientation measurements |
CN1818555A (zh) * | 2006-03-29 | 2006-08-16 | 北京航空航天大学 | 安装误差角与标度因数解耦的微惯性测量单元精确标定方法 |
CN112166685B (zh) * | 2008-09-12 | 2012-10-31 | 中国电子科技集团公司第二十六研究所 | 加速度计组合及其控制方法 |
CN102889076A (zh) * | 2012-09-14 | 2013-01-23 | 西安思坦仪器股份有限公司 | 陀螺测斜仪标定方法 |
CN104133081A (zh) * | 2014-07-31 | 2014-11-05 | 工业和信息化部电子第五研究所 | 加速度计性能参数的检测方法和系统 |
CN106813680A (zh) * | 2016-12-28 | 2017-06-09 | 兰州空间技术物理研究所 | 一种高精度、高分辨率石英扰性传感器的静态标定方法 |
CN109813336A (zh) * | 2017-11-22 | 2019-05-28 | 广东虚拟现实科技有限公司 | 惯性测量单元标定方法 |
CN110006450A (zh) * | 2019-04-15 | 2019-07-12 | 哈尔滨工业大学 | 一种激光捷联惯导系统在卧式三轴转台上的标定方法 |
Non-Patent Citations (1)
Title |
---|
毛奔,张晓宇: "《微惯性系统及应用》", 31 July 2013, 哈尔滨工程大学出版社 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113376404A (zh) * | 2021-05-19 | 2021-09-10 | 中寰卫星导航通信有限公司 | 一种标定验证系统、设备、方法及存储介质 |
CN114397479A (zh) * | 2021-12-07 | 2022-04-26 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种石英振梁加速度计to部件的两自由度测试装置及方法 |
CN114397479B (zh) * | 2021-12-07 | 2024-05-10 | 华中光电技术研究所(中国船舶重工集团公司第七一七研究所) | 一种石英振梁加速度计to部件的两自由度测试装置及方法 |
CN114324978A (zh) * | 2021-12-17 | 2022-04-12 | 兰州空间技术物理研究所 | 一种加速度计捕获范围的地面静态标定方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110006450B (zh) | 一种激光捷联惯导系统在卧式三轴转台上的标定方法 | |
CN111273058A (zh) | 一种加速度计标定方法 | |
CN106052595B (zh) | 基于激光陀螺捷联惯导的三轴转台轴线垂直度检测方法 | |
CN106017507B (zh) | 一种用于中低精度的光纤惯组快速标定方法 | |
CN110108300B (zh) | 一种基于卧式三轴转台的imu正六面体标定方法 | |
CN105043412B (zh) | 一种惯性测量单元误差补偿方法 | |
CN101290326B (zh) | 石英挠性加速度计测量组件的参数辨识标定方法 | |
CN106969783B (zh) | 一种基于光纤陀螺惯性导航的单轴旋转快速标定技术 | |
CN107121707A (zh) | 一种三轴磁传感器测量基准与结构基准的误差校正方法 | |
CN104316037B (zh) | 一种电子罗盘的校正方法及装置 | |
WO2020164206A1 (zh) | 一种旋转加速度计重力梯度仪标定方法 | |
CN110006454B (zh) | 一种imu标定三轴转台垂直度和初始姿态的方法 | |
CN112698258B (zh) | 一种三轴磁强计的一体化误差校正方法 | |
CN107490809A (zh) | 一种重力敏感器的惯组级标定方法 | |
CN101629830A (zh) | 一种三轴一体高精度光纤陀螺的标定方法及装置 | |
US6758080B1 (en) | Calibration of sensors | |
CN108088431B (zh) | 一种自校正电子罗盘及其校正方法 | |
CN110631605B (zh) | 一种陀螺阵列标定方法及系统 | |
CN109059917B (zh) | 一种动态水平仪及其动态调整测量方法 | |
CN103823083B (zh) | 提高加速度计校准精度的方法与系统 | |
CN108716925A (zh) | 一种九轴传感器的标定方法及装置 | |
CN108036756B (zh) | 一种利用加速度计进行双轴旋转惯性测量装置相邻轴线垂直度检查的方法 | |
CN108168516B (zh) | 基于光纤陀螺测量待测台面与基准水平面之间倾斜夹角的方法 | |
CN112902954A (zh) | 塔架晃动传感器和塔架晃动角度测量方法 | |
CN106370206B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200612 |