CN104330095A - 基于分段拟合的mems陀螺温度补偿方法 - Google Patents
基于分段拟合的mems陀螺温度补偿方法 Download PDFInfo
- Publication number
- CN104330095A CN104330095A CN201410578542.7A CN201410578542A CN104330095A CN 104330095 A CN104330095 A CN 104330095A CN 201410578542 A CN201410578542 A CN 201410578542A CN 104330095 A CN104330095 A CN 104330095A
- Authority
- CN
- China
- Prior art keywords
- temperature
- mems gyro
- zero
- constant multiplier
- compensation
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Gyroscopes (AREA)
Abstract
本发明涉及一种基于分段拟合的MEMS陀螺温度补偿方法,首先,将MEMS陀螺的全温范围划分为若干段温度段;然后,在每个温度段中分别对零偏和标度因数进行拟合分别得到该温度段内的零偏和标度因数的关系式,从而获得全温范围内MEMS陀螺的零偏和标度因数的关系式;最后,计算MEMS陀螺在任意一温度点下的零偏和标度因数,再根据MEMS陀螺在该温度点下的补偿前的输出以及该温度点下的零偏和标度因数计算获得MEMS陀螺在该温度点下的补偿后的输出。本发明对MEMS陀螺全温范围内的温度进行分段拟合补偿,可以不用考虑零偏和标度因数的温度变化曲线不规则以及难以实现建模和多项式拟合的问题,解决了基于建模和多项式拟合的全温零偏和标度因数补偿方法存在的局限性问题。
Description
技术领域
本发明涉及一种MEMS陀螺的零偏和标度因数温度补偿方法。
背景技术
MEMS陀螺的零偏和标度因数是随温度变化的,对MEMS陀螺进行全温下的补偿可以提高MEMS陀螺在全温下的测量精度。
目前,对MEMS陀螺进行温度补偿的方法主要有:对温度点进行性查表、多项式系数拟合、神经网络等方法。由于这些温度补偿方法需要建立MEMS陀螺输出和温度之间的模型,但是MEMS陀螺的零偏和标度因数在全温范围内的变化曲线有时候是很不规则的,难以通过单一模型表达,故而无法通过以上方法来实现全温下的温度补偿。
发明内容
本发明的目的是提供一种能够实现在全温度范围内对MEMS陀螺的零偏和标度因数进行温度补偿的方法。
为达到上述目的,本发明采用的技术方案是:
一种基于分段拟合的MEMS陀螺温度补偿方法,用于对MEMS陀螺的零偏和标度因数进行温度补偿,该方法为:首先,将所述的MEMS陀螺的全温范围划分为若干段温度段,每个温度段的两端点的温度分别为Ti和Tj;然后,在每个温度段中分别对零偏和标度因数进行拟合分别得到该温度段内的零偏和标度因数的关系式,从而获得全温范围内所述的MEMS陀螺的零偏和标度因数的关系式;最后,根据所获得全温范围内所述的MEMS陀螺的零偏和标度因数的关系式计算所述的MEMS陀螺在任意一温度点T下的零偏BiasT和标度因数SFT,再根据所述的MEMS陀螺在该温度点T下的补偿前的输出GyroOutT以及该温度点T下的零偏BiasT和标度因数SFT计算获得所述的MEMS陀螺在该温度点T下的补偿后的输出RateT。
所述的MEMS陀螺的全温范围均分为若干个所述的温度段。
每个所述的温度段内的温度差为10℃。
在每个所述的温度段内,其两端点温度Ti和Tj下的零偏分别为和标度因数分别为和则
该温度段内的零偏
该温度段内的标度因数
其中,Ti~j为该温度段内的任意温度;
零偏的线性拟合斜率
零偏的线性拟合截距或
标度因数的线性拟合斜率
标度因数的线性拟合截距或
所述的MEMS陀螺的补偿后的输出
所述的MEMS陀螺的输出量为角速率。
由于上述技术方案运用,本发明与现有技术相比具有下列优点:本发明对MEMS陀螺全温范围内的温度进行分段拟合补偿,可以不用考虑零偏和标度因数的温度变化曲线不规则以及难以实现建模和多项式拟合的问题,解决了基于建模和多项式拟合的全温零偏和标度因数补偿方法存在的局限性问题。
具体实施方式
下面结合实施例对本发明作进一步描述。
实施例一:一种基于分段拟合的MEMS陀螺温度补偿方法,其用于对MEMS陀螺的零偏和标度因数进行温度补偿。该方法为:
(1)将MEMS陀螺的全温范围平均划分为若干段温度段,每个温度段的两端点的温度分别为Ti和Tj。例如,MEMS陀螺的测试的全温范围设定为-40℃至+90℃,在该全温范围内,以10℃为温差标准进行分段,故而可以划分的温度段分别为:-40℃至-30℃、-30℃至-20℃、-20℃至-10℃、-10℃至0℃、0℃至10℃、10℃至20℃、20℃至30℃、30℃至40℃、40℃至50℃、50℃至60℃、60℃至70℃、70℃至80℃、80℃至90℃。各个温度段的端点温度作为实验温度,分别为:-40℃,-30℃,-20℃,-10℃,0℃,10℃,20℃,30℃,40℃,50℃,60℃,70℃,80℃,90℃。上述每个端点按照MEMS陀螺标度因数测试方法进行实验,计算出每个端点温度下的标度因数和零偏。
(2)在每个温度段中分别对零偏和标度因数进行拟合分别得到该温度段内的零偏和标度因数的关系式,从而获得全温范围内MEMS陀螺的零偏和标度因数的关系式。
在每个温度段内,设其两端点温度Ti和Tj下的零偏分别为和而标度因数分别为和则
该温度段内的零偏
该温度段内的标度因数
其中,Ti~j为该温度段内的任意温度;
零偏的线性拟合斜率
零偏的线性拟合截距或
标度因数的线性拟合斜率
标度因数的线性拟合截距或
不妨以-40℃至-30℃的温度段为例,首先获得端点的-40℃和-30℃的零偏Bias-40和Bias-30,再计算出-40℃至-30℃的温度段的零偏Bias-40~-30和标度因数SF-40~-30的变化关系式:
Bias-40~-30=KB×T+CB;
SF-40~-30=KS×T+CS。
上式中,T为-40℃至-30℃的温度段内的温度,KB为-40℃至-30℃的温度段内的零偏线性拟合斜率,CB为-40℃至-30℃的温度段内的零偏线性拟合截距,KS为-40℃至-30℃的温度段内的标度因数线性拟合斜率,CS为-40℃至-30℃的温度段内的标度因数线性拟合截距。则
CB=Bias-30-KB×(-30℃)或CB=Bias-40-KB×(-40℃);
CS=SF-30-KS×(-30℃)或CS=SF-40-KS×(-40℃)。
得到每个温度段内的零偏和标度因数的关系式并组合即可获得全温范围内MEMS陀螺的零偏和标度因数的关系式。
(3)根据所获得全温范围内MEMS陀螺的零偏和标度因数的关系式计算MEMS陀螺在任意一温度点T下的零偏BiasT=KB×T+CB和标度因数SFT=KS×T+CS,其中KB为该温度点T所在的温度段内的零偏线性拟合斜率,CB为该温度点T所在的温度段内的零偏线性拟合截距,KS为该温度点T所在的温度段内的标度因数线性拟合斜率,CS为该温度点T所在的温度段内的标度因数线性拟合截距;再根据MEMS陀螺在该温度点T下的补偿前的角速率输出GyroOutT以及该温度点T下的零偏BiasT和标度因数SFT计算获得MEMS陀螺在该温度点T下的补偿后的角速率输出RateT。
MEMS陀螺的补偿后的输出
本发明的原理是通过对温度点进行均匀分段,并用线性关系拟合每段温度段中零偏和标度因数的变化关系,实现MEMS陀螺的全温补偿。采用上述温度补偿方法,可以实现MEMS陀螺的全温温度补偿。若减小全温范围内的温差标准,使得温度分段更细致,可以进一步提高全温温度补偿的精度。本发明的方法可以推广到MEMS加速度计的全温温度补偿。
上述实施例只为说明本发明的技术构思及特点,其目的在于让熟悉此项技术的人士能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡根据本发明精神实质所作的等效变化或修饰,都应涵盖在本发明的保护范围之内。
Claims (6)
1.一种基于分段拟合的MEMS陀螺温度补偿方法,用于对MEMS陀螺的零偏和标度因数进行温度补偿,其特征在于:该方法为:首先,将所述的MEMS陀螺的全温范围划分为若干段温度段,每个温度段的两端点的温度分别为Ti和Tj;然后,在每个温度段中分别对零偏和标度因数进行拟合分别得到该温度段内的零偏和标度因数的关系式,从而获得全温范围内所述的MEMS陀螺的零偏和标度因数的关系式;最后,根据所获得全温范围内所述的MEMS陀螺的零偏和标度因数的关系式计算所述的MEMS陀螺在任意一温度点T下的零偏BiasT和标度因数SFT,再根据所述的MEMS陀螺在该温度点T下的补偿前的输出GyroOutT以及该温度点T下的零偏BiasT和标度因数SFT计算获得所述的MEMS陀螺在该温度点T下的补偿后的输出RateT。
2.根据权利要求1所述的基于分段拟合的MEMS陀螺温度补偿方法,其特征在于:所述的MEMS陀螺的全温范围均分为若干个所述的温度段。
3.根据权利要求2所述的基于分段拟合的MEMS陀螺温度补偿方法,其特征在于:每个所述的温度段内的温度差为10℃。
4.根据权利要求1所述的基于分段拟合的MEMS陀螺温度补偿方法,其特征在于:在每个所述的温度段内,其两端点温度Ti和Tj下的零偏分别为和标度因数分别为和则
该温度段内的零偏
该温度段内的标度因数
其中,Ti~j为该温度段内的任意温度;
零偏的线性拟合斜率
零偏的线性拟合截距或
标度因数的线性拟合斜率
标度因数的线性拟合截距或
5.根据权利要求1或4所述的基于分段拟合的MEMS陀螺温度补偿方法,其特征在于:所述的MEMS陀螺的补偿后的输出
6.根据权利要求1所述的基于分段拟合的MEMS陀螺温度补偿方法,其特征在于:所述的MEMS陀螺的输出量为角速率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410578542.7A CN104330095A (zh) | 2014-10-24 | 2014-10-24 | 基于分段拟合的mems陀螺温度补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410578542.7A CN104330095A (zh) | 2014-10-24 | 2014-10-24 | 基于分段拟合的mems陀螺温度补偿方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104330095A true CN104330095A (zh) | 2015-02-04 |
Family
ID=52404866
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410578542.7A Pending CN104330095A (zh) | 2014-10-24 | 2014-10-24 | 基于分段拟合的mems陀螺温度补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104330095A (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105222765A (zh) * | 2015-09-18 | 2016-01-06 | 工业和信息化部电子第五研究所 | Mems陀螺的温度补偿方法及系统 |
CN105277215A (zh) * | 2015-11-19 | 2016-01-27 | 中国兵器工业集团第二一四研究所苏州研发中心 | 基于离心加速度的mems陀螺加速度敏感性补偿方法 |
CN105333888A (zh) * | 2015-11-26 | 2016-02-17 | 湖北三江航天红峰控制有限公司 | 一种利用一次温度实验同时补偿光纤陀螺标度因数及零偏的方法 |
CN105352630A (zh) * | 2015-12-24 | 2016-02-24 | 深圳市博巨兴实业发展有限公司 | 一种用于温度传感器芯片的分段线性校准系统及方法 |
CN105547325A (zh) * | 2015-12-08 | 2016-05-04 | 北京航天时代光电科技有限公司 | 一种基于k均值聚类的光纤惯组温度模型系数确定方法 |
CN105841715A (zh) * | 2016-03-18 | 2016-08-10 | 北京信息科技大学 | 一种高动态双轴角速率陀螺及零偏和标度因数误差补偿 |
CN106123924A (zh) * | 2016-08-23 | 2016-11-16 | 美新微纳传感系统有限公司 | 一种陀螺仪的温度补偿方法 |
CN106546267A (zh) * | 2016-10-27 | 2017-03-29 | 陕西航天时代导航设备有限公司 | 一种速率陀螺的零位温度标定方法 |
CN109163736A (zh) * | 2018-11-12 | 2019-01-08 | 西安微电子技术研究所 | 一种光纤imu全温标定补偿方法 |
CN111238462A (zh) * | 2020-01-19 | 2020-06-05 | 湖北三江航天红峰控制有限公司 | 一种基于深度嵌入聚类的lstm光纤陀螺温补建模方法 |
CN111982101A (zh) * | 2020-08-03 | 2020-11-24 | 中国兵器工业集团第二一四研究所苏州研发中心 | 一种大量程加速度计标度因数非线性补偿方法 |
CN112254742A (zh) * | 2020-10-13 | 2021-01-22 | 天津津航计算技术研究所 | Mems惯性器件温度补偿参数在线拟合方法 |
CN112729333A (zh) * | 2020-11-25 | 2021-04-30 | 北京仿真中心 | 一种适用于硅微陀螺仪的分段多项式数字温度补偿方法 |
CN113418521A (zh) * | 2021-06-10 | 2021-09-21 | 北京航天时代激光导航技术有限责任公司 | 一种提高激光陀螺仪刻度因数长期稳定性的方法 |
CN115047213A (zh) * | 2021-03-09 | 2022-09-13 | 北京大学 | 一种提高mems加速度计长期稳定性的方法 |
CN116429150A (zh) * | 2023-06-15 | 2023-07-14 | 河北美泰电子科技有限公司 | Mems陀螺零偏补偿方法、装置、电子设备及存储介质 |
CN117109566A (zh) * | 2023-08-23 | 2023-11-24 | 长春航盛艾思科电子有限公司 | 一种基于分段多项式拟合的imu温度补偿方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7726188B2 (en) * | 2006-04-24 | 2010-06-01 | Millisensor Systems + Actuators | Scale factor measurement for mems gyroscopes and accelerometers |
CN102230806A (zh) * | 2011-04-26 | 2011-11-02 | 中国科学院软件研究所 | 一种陀螺仪温度漂移补偿方法 |
CN103196462A (zh) * | 2013-02-28 | 2013-07-10 | 南京航空航天大学 | 一种mimu中mems陀螺仪的误差标定补偿方法 |
-
2014
- 2014-10-24 CN CN201410578542.7A patent/CN104330095A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7726188B2 (en) * | 2006-04-24 | 2010-06-01 | Millisensor Systems + Actuators | Scale factor measurement for mems gyroscopes and accelerometers |
CN102230806A (zh) * | 2011-04-26 | 2011-11-02 | 中国科学院软件研究所 | 一种陀螺仪温度漂移补偿方法 |
CN103196462A (zh) * | 2013-02-28 | 2013-07-10 | 南京航空航天大学 | 一种mimu中mems陀螺仪的误差标定补偿方法 |
Non-Patent Citations (5)
Title |
---|
夏敦柱等: "硅微陀螺仪零偏温度性能补控方法设计", 《东南大学学报》 * |
张桂才: "《光纤陀螺原理与技术》", 31 May 2008 * |
杨亮等: "高品质因数微机械陀螺的温度自补偿", 《光学精密工程》 * |
毛奔等: "《微惯性系统及应用》", 31 July 2013 * |
魏泽松等: "温度对IMU微机械陀螺仪零偏影响及标定补偿", 《理论与方法》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105222765B (zh) * | 2015-09-18 | 2018-06-12 | 工业和信息化部电子第五研究所 | Mems陀螺的温度补偿方法及系统 |
CN105222765A (zh) * | 2015-09-18 | 2016-01-06 | 工业和信息化部电子第五研究所 | Mems陀螺的温度补偿方法及系统 |
CN105277215A (zh) * | 2015-11-19 | 2016-01-27 | 中国兵器工业集团第二一四研究所苏州研发中心 | 基于离心加速度的mems陀螺加速度敏感性补偿方法 |
CN105277215B (zh) * | 2015-11-19 | 2017-12-12 | 中国兵器工业集团第二一四研究所苏州研发中心 | 基于离心加速度的mems陀螺加速度敏感性补偿方法 |
CN105333888A (zh) * | 2015-11-26 | 2016-02-17 | 湖北三江航天红峰控制有限公司 | 一种利用一次温度实验同时补偿光纤陀螺标度因数及零偏的方法 |
CN105333888B (zh) * | 2015-11-26 | 2018-02-23 | 湖北三江航天红峰控制有限公司 | 一种利用一次温度实验同时补偿光纤陀螺标度因数及零偏的方法 |
CN105547325A (zh) * | 2015-12-08 | 2016-05-04 | 北京航天时代光电科技有限公司 | 一种基于k均值聚类的光纤惯组温度模型系数确定方法 |
CN105547325B (zh) * | 2015-12-08 | 2018-06-19 | 北京航天时代光电科技有限公司 | 一种基于k均值聚类的光纤惯组温度模型系数确定方法 |
CN105352630A (zh) * | 2015-12-24 | 2016-02-24 | 深圳市博巨兴实业发展有限公司 | 一种用于温度传感器芯片的分段线性校准系统及方法 |
CN105841715A (zh) * | 2016-03-18 | 2016-08-10 | 北京信息科技大学 | 一种高动态双轴角速率陀螺及零偏和标度因数误差补偿 |
CN106123924B (zh) * | 2016-08-23 | 2018-10-30 | 新纳传感系统有限公司 | 一种陀螺仪的温度补偿方法 |
CN106123924A (zh) * | 2016-08-23 | 2016-11-16 | 美新微纳传感系统有限公司 | 一种陀螺仪的温度补偿方法 |
CN106546267A (zh) * | 2016-10-27 | 2017-03-29 | 陕西航天时代导航设备有限公司 | 一种速率陀螺的零位温度标定方法 |
CN109163736A (zh) * | 2018-11-12 | 2019-01-08 | 西安微电子技术研究所 | 一种光纤imu全温标定补偿方法 |
CN109163736B (zh) * | 2018-11-12 | 2020-06-30 | 西安微电子技术研究所 | 一种光纤imu全温标定补偿方法 |
CN111238462A (zh) * | 2020-01-19 | 2020-06-05 | 湖北三江航天红峰控制有限公司 | 一种基于深度嵌入聚类的lstm光纤陀螺温补建模方法 |
CN111982101A (zh) * | 2020-08-03 | 2020-11-24 | 中国兵器工业集团第二一四研究所苏州研发中心 | 一种大量程加速度计标度因数非线性补偿方法 |
CN112254742A (zh) * | 2020-10-13 | 2021-01-22 | 天津津航计算技术研究所 | Mems惯性器件温度补偿参数在线拟合方法 |
CN112254742B (zh) * | 2020-10-13 | 2022-08-09 | 天津津航计算技术研究所 | Mems惯性器件温度补偿参数在线拟合方法 |
CN112729333A (zh) * | 2020-11-25 | 2021-04-30 | 北京仿真中心 | 一种适用于硅微陀螺仪的分段多项式数字温度补偿方法 |
CN115047213A (zh) * | 2021-03-09 | 2022-09-13 | 北京大学 | 一种提高mems加速度计长期稳定性的方法 |
CN113418521A (zh) * | 2021-06-10 | 2021-09-21 | 北京航天时代激光导航技术有限责任公司 | 一种提高激光陀螺仪刻度因数长期稳定性的方法 |
CN116429150A (zh) * | 2023-06-15 | 2023-07-14 | 河北美泰电子科技有限公司 | Mems陀螺零偏补偿方法、装置、电子设备及存储介质 |
CN117109566A (zh) * | 2023-08-23 | 2023-11-24 | 长春航盛艾思科电子有限公司 | 一种基于分段多项式拟合的imu温度补偿方法 |
CN117109566B (zh) * | 2023-08-23 | 2024-01-23 | 长春航盛艾思科电子有限公司 | 一种基于分段多项式拟合的imu温度补偿方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104330095A (zh) | 基于分段拟合的mems陀螺温度补偿方法 | |
CN104765977A (zh) | 考虑局部水头损失的河段综合糙率的计算方法 | |
WO2014181183A9 (en) | System and method for multi-phase fluid measurement | |
Hirschhorn | Overcategories and undercategories of model categories | |
MX2017002640A (es) | Metodo y aparato para determinar densidad diferencial. | |
Weidman et al. | Flows induced by a plate moving normal to stagnation-point flow | |
CN104040378B (zh) | 气象预测装置以及气象预测方法 | |
CN103684350A (zh) | 一种粒子滤波方法 | |
CN106546267A (zh) | 一种速率陀螺的零位温度标定方法 | |
Hameed et al. | Calibration of Manning’s friction factor for rivers in Iraq using hydraulic model (Al-Kufa River as case study) | |
Mehdiyeva et al. | Information-measurement system development for controlling of parameters of measurement | |
CN107025209A (zh) | 一种不等间隔采样的离散数据一阶导数解算方法 | |
Sherstyukov et al. | Reciprocal expansion of modified Bessel function in simple fractions and obtaining general summation relationships containing its zeros | |
Goudar et al. | Variability in propagation path delay for atmospheric remote sensing | |
Ma et al. | Iterative algorithm of conjugate depth for parabolic channels | |
CN105136095B (zh) | 一种树干任意点形状指数的测定方法 | |
Aiguo et al. | Research on application of Rayleigh scattering velocity measurement in hypersonic low density wind tunnel | |
CN104077472B (zh) | 一种利用加速度计组合输出离散度进行精度评估的方法 | |
Lv et al. | An optimization iterative algorithm based on nonnegative constraint with application to Allan variance analysis technique | |
Montecinos | Analytic solutions for the Burgers equation with source terms | |
Zudin | The velocity of gas bubble rise in a tube | |
Nourani et al. | Investigation of the Stage-Discharge Relation and Discharge Coefficient in Sharp-Crested Weirs with Triangular Shape in Plan | |
Kim et al. | Improvement of Uncertainty Assessment of Discharge Estimated by Velocity-Area Method | |
Xin et al. | Supercritical Airfoil Data Corrected for Sidewall Boundary-Layer Effects in Transonic Wind Tunnel | |
Lesovoi et al. | Optimization task of flowmeter designing in view of roughness of pipeline internal surface |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150204 |
|
RJ01 | Rejection of invention patent application after publication |