CN102654515B - Calibration algorithm for z sensitive shaft of three-shaft acceleration transducer - Google Patents

Calibration algorithm for z sensitive shaft of three-shaft acceleration transducer Download PDF

Info

Publication number
CN102654515B
CN102654515B CN201110051829.0A CN201110051829A CN102654515B CN 102654515 B CN102654515 B CN 102654515B CN 201110051829 A CN201110051829 A CN 201110051829A CN 102654515 B CN102654515 B CN 102654515B
Authority
CN
China
Prior art keywords
formula
sensitive
acceleration sensor
axis acceleration
shaft
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
Application number
CN201110051829.0A
Other languages
Chinese (zh)
Other versions
CN102654515A (en
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.)
Meixin Semiconductor Wuxi Co Ltd
Original Assignee
Memsic Transducer Systems Co Ltd
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 Memsic Transducer Systems Co Ltd filed Critical Memsic Transducer Systems Co Ltd
Priority to CN201110051829.0A priority Critical patent/CN102654515B/en
Publication of CN102654515A publication Critical patent/CN102654515A/en
Application granted granted Critical
Publication of CN102654515B publication Critical patent/CN102654515B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Force Measurement Appropriate To Specific Purposes (AREA)

Abstract

The invention relates to a calibration method for a z sensitive shaft of a three-shaft acceleration transducer, which comprises the following steps of: setting the zero points of three sensitive shafts x, y and z of the three-shaft acceleration transducer as bx, by and bz, and setting sensitive degrees of the three sensitive shafts as sx, sy and sz, wherein the unit of the sensitive degree is the gravitational acceleration, and the zero point and the sensitive degree of the sensitive shaft z is bz and sz. Because the accuracy degrees of bz and sz are poor, bz and sz are regarded as unknown numbers. The formula is set as a*z<2>+b*z+c=g(x,y) (1). The calibration method also comprises the following steps of: adopting a fixed sampling frequency to sample three groups of initial data (namely x1, y1, z1, x2, y2, z2, x2, x3, y3 and z3) of the three-shaft acceleration transducer under different static conditions, substituting three groups of initial data into the formula (1) to obtain parameters a, b and c, and then calculating the actual zero point bz and the actual sensitive degree sz of the sensitive shaft Z of the three-shaft acceleration transducer.

Description

The calibration steps of the z sensitive axes of 3-axis acceleration sensor
Technical field
The present invention relates to a kind of calibration algorithm of z sensitive axes of 3-axis acceleration sensor.
Background technology
The popularity rate of current 3-axis acceleration sensor in consumer electronics is more and more higher, and thereupon its competition is also more and more fierce.Thus, the cost pressure of 3-axis acceleration sensor is increasing.And then the 3-axis acceleration sensor that also just result in some low costs more and more occurs.
In some low cost 3-axis acceleration sensors, acceleration sensitive axle performance parameter (zero point and sensitivity) on x-axis and y-axis direction generally can be guaranteed, but z sensitive axes is due to singularity such as techniques, its performance can be more far short of what is expected than x sensitive axes and y sensitive axes.That is, there is very large error and instability in zero point of its z sensitive axes and sensitivity.This can bring very large error to use, sometimes even cannot meet actual user demand.
Therefore, when using this type of acceleration transducer, due to the z sensitive axes that it is poor, bring a lot of errors and uncertainty can to follow-up calculating and process.
Summary of the invention
Technical matters to be solved by this invention is, the z sensitive axes poor-performing of existing low cost 3-axis acceleration sensor, exists the defect such as influence error and instability.The present invention aims to provide a kind of z sensitive axes calibration algorithm of 3-axis acceleration sensor, its can be larger to 3-axis acceleration sensor error in use z sensitive axes calibrate, obtain its real zero point and sensitivity, thus when not increasing processing step and cost, improve the performance of 3-axis acceleration sensor.
In order to solve the problems of the technologies described above, technical scheme proposed by the invention is:
The calibration algorithm of the z sensitive axes of 3-axis acceleration sensor, it includes following steps:
The x of setting 3-axis acceleration sensor, y, z tri-zero point of sensitive axes be followed successively by bx, by, bz, sensitivity is followed successively by sx, sy, sz, its unit is an acceleration of gravity, wherein z-axis zero point bz and sensitivity sz due to its poor accuracy, be regarded as unknown number at this;
Setting formula:
a*z 2+b*z+c=g(x,y) (1)
Wherein:
a = 1 sz 2
b = - 2 * bz sz 2 - - - ( 2 )
c = bz 2 sz 2
With 3 group raw data: the xs of fixed sampling frequency sampling 3-axis acceleration sensor under different attitude stationary state 1, y 1, z 1, x 2, y 2, z 2, x 3, y 3, z 3; Then can be obtained by formula (1):
z 1 2 z 1 1 z 2 2 z 2 1 z 3 2 z 3 1 a b c = g ( x 1 , y 1 ) g ( x 2 , y 2 ) g ( x 3 , y 3 ) - - - ( 3 )
Calculate the value of parameter a, b, c according to formula (3), then can calculate 3-axis acceleration sensor z sensitive axes bz at real zero point and sensitivity sz by formula (2).
Further, in different embodiments, wherein three axle raw data x, y, z meet following formula:
( x - bx sx ) 2 + ( y - by sy ) 2 + ( z - bz sz ) 2 = 1 - - - ( 1 - 1 )
Further, in different embodiments, wherein following formula is set by formula (1-1):
g ( x , y ) = 1 - ( x - bx sx ) 2 + ( y - by sy ) 2 - - - ( 1 - 2 )
Formula (1-2) is substituted into formula (1-1) can obtain:
( z - bz sz ) 2 = g ( x , y ) - - - ( 1 - 3 )
And formula (1) is namely converted by formula (1-3).
Further, in different embodiments, wherein formula (3) can change into further:
a b c = z 1 2 z 1 1 z 2 2 z 2 1 z 3 2 z 3 1 - 1 g ( x 1 , y 1 ) g ( x 2 , y 2 ) g ( x 3 , y 3 ) - - - ( 3 - 1 )
Calculate the value of parameter a, b, c according to formula (3-1), then can calculate 3-axis acceleration sensor z sensitive axes bz at real zero point and sensitivity sz by formula (2).
Compared with prior art, the invention has the beneficial effects as follows: the calibration algorithm that use the present invention relates to can very easily in the use procedure of 3-axis acceleration sensor, and the z sensitive axes larger to its error is calibrated, and obtains real zero point and sensitivity.Thus, the performance of 3-axis acceleration sensor is improved when not increasing processing step and cost.
Embodiment
An embodiment of the invention provide the calibration algorithm that a kind of acceleration transducer to having poor z sensitive axes carries out z sensitive axes zero point and sensitivity, and it includes following steps:
According to the databook of 3-axis acceleration sensor, can x be obtained, y, z tri-zero point of sensitive axes be followed successively by bx, by, bz, sensitivity is followed successively by sx, sy, sz, and its unit is an acceleration of gravity.Wherein zero point of z-axis and sensitivity are due to its poor accuracy, and be regarded as unknown number at this, the object of this algorithm is exactly the actual value calculating this Two Variables.Because it is subject to the effect of acceleration of gravity under 3-axis acceleration sensor static position, therefore, its three axle raw data x, y, z meet following formula:
( x - bx sx ) 2 + ( y - by sy ) 2 + ( z - bz sz ) 2 = 1 - - - ( 1 - 1 )
3-axis acceleration sensor three axle raw data x, the unit of y, z is an acceleration of gravity.If:
g ( x , y ) = 1 - ( x - bx sx ) 2 + ( y - by sy ) 2 - - - ( 1 - 2 )
Formula (1-2) can be obtained on behalf of formula (1-1):
( z - bz sz ) 2 = g ( x , y ) - - - ( 1 - 3 )
Formula (1-3) is converted to following form:
a*z 2+b*z+c=g(x,y) (1)
Wherein:
a = 1 sz 2
b = - 2 * bz sz 2 - - - ( 2 )
c = bz 2 sz 2
With 3 group raw data: the xs of fixed sampling frequency sampling 3-axis acceleration sensor under different attitude stationary state 1, y 1, z 1, x 2, y 2, z 2, x 3, y 3, z 3; Then can be obtained by formula (2):
z 1 2 z 1 1 z 2 2 z 2 1 z 3 2 z 3 1 a b c = g ( x 1 , y 1 ) g ( x 2 , y 2 ) g ( x 3 , y 3 ) - - - ( 3 )
Formula (3) can change into further:
a b c = z 1 2 z 1 1 z 2 2 z 2 1 z 3 2 z 3 1 - 1 g ( x 1 , y 1 ) g ( x 2 , y 2 ) g ( x 3 , y 3 ) - - - ( 3 - 1 )
Calculate the value of parameter a, b, c according to formula (3-1), then can calculate 3-axis acceleration sensor z sensitive axes bz at real zero point and sensitivity sz by formula (2).
The calibration algorithm that the present invention relates to can easily in the use procedure of 3-axis acceleration sensor, and the z sensitive axes larger to its error is calibrated, and obtains real zero point and sensitivity.Thus, the performance of 3-axis acceleration sensor is improved when not increasing processing step and cost.
The foregoing is only better embodiment of the present invention; protection scope of the present invention is not limited with above-mentioned embodiment; in every case those of ordinary skill in the art modify or change according to the equivalence that disclosure of the present invention is done, and all should include in the protection domain recorded in claims.

Claims (4)

1. the calibration steps of the z sensitive axes of 3-axis acceleration sensor, is characterized in that: it includes following steps:
The zero point of x, y, z three sensitive axes of setting 3-axis acceleration sensor is followed successively by bx, by, bz, sensitivity is followed successively by sx, sy, sz, its unit is an acceleration of gravity, wherein z-axis zero point bz and sensitivity sz due to its poor accuracy, be regarded as unknown number at this;
Setting formula:
a*z 2+b*z+c=g(x,y) (1)
Wherein:
a = 1 sz 2 b = - 2 * bz sz 2 c = bz 2 sz 2 - - - ( 2 )
With 3 group raw data: the xs of fixed sampling frequency sampling 3-axis acceleration sensor under different attitude stationary state 1, y 1, z 1, x 2, y 2, z 2, x 3, y 3, z 3; Then can be obtained by formula (1):
z 1 2 z 1 1 z 2 2 z 2 1 z 3 2 z 3 1 a b c = g ( x 1 , y 1 ) g ( x 2 , y 2 ) g ( x 3 , y 3 ) - - - ( 3 )
Calculate the value of parameter a, b, c according to formula (3), then can calculate 3-axis acceleration sensor z sensitive axes bz at real zero point and sensitivity sz by formula (2).
2. the calibration steps of the z sensitive axes of 3-axis acceleration sensor as claimed in claim 1, it is characterized in that: wherein three axle raw data x, y, z meet following formula:
( x - bx sx ) 2 + ( y - by sy ) 2 + ( z - bz sz ) 2 = 1 - - - ( 1 - 1 ) .
3. the calibration steps of the z sensitive axes of 3-axis acceleration sensor as claimed in claim 2, is characterized in that: wherein set following formula by formula (1-1):
g ( x , y ) = 1 - ( x - bx sx ) 2 + ( y - by sy ) 2 - - - ( 1 - 2 )
Formula (1-2) is substituted into formula (1-1) can obtain:
( z - bz sz ) 2 = g ( x , y ) - - - ( 1 - 3 )
And formula (1) is namely converted by formula (1-3).
4. the calibration steps of the z sensitive axes of 3-axis acceleration sensor as claimed in claim 1, is characterized in that: wherein formula (3) can change into further:
a b c = z 1 2 z 1 1 z 2 2 z 2 1 z 3 2 z 3 1 - 1 g ( x 1 , y 1 ) g ( x 2 , y 2 ) g ( x 3 , y 3 ) - - - ( 3 - 1 )
Calculate the value of parameter a, b, c according to formula (3-1), then can calculate 3-axis acceleration sensor z sensitive axes bz at real zero point and sensitivity sz by formula (2).
CN201110051829.0A 2011-03-04 2011-03-04 Calibration algorithm for z sensitive shaft of three-shaft acceleration transducer Active CN102654515B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110051829.0A CN102654515B (en) 2011-03-04 2011-03-04 Calibration algorithm for z sensitive shaft of three-shaft acceleration transducer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110051829.0A CN102654515B (en) 2011-03-04 2011-03-04 Calibration algorithm for z sensitive shaft of three-shaft acceleration transducer

Publications (2)

Publication Number Publication Date
CN102654515A CN102654515A (en) 2012-09-05
CN102654515B true CN102654515B (en) 2015-04-15

Family

ID=46730170

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110051829.0A Active CN102654515B (en) 2011-03-04 2011-03-04 Calibration algorithm for z sensitive shaft of three-shaft acceleration transducer

Country Status (1)

Country Link
CN (1) CN102654515B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103823084A (en) * 2014-03-21 2014-05-28 苏州纳芯微电子有限公司 Method for calibrating three-axis acceleration sensor
CN105759078B (en) * 2016-04-12 2019-01-04 北京荣之联科技股份有限公司 The three axis method for self-calibrating and device of vehicle-mounted acceleration sensor
CN106199758A (en) * 2016-06-30 2016-12-07 联想(北京)有限公司 Measurement data calibration steps and electronic equipment
CN106491138B (en) * 2016-10-26 2019-04-09 歌尔科技有限公司 A kind of motion state detection method and device
CN109521941A (en) * 2018-10-25 2019-03-26 西安闻泰电子科技有限公司 Data sampling device and method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101246024A (en) * 2008-03-26 2008-08-20 北京航空航天大学 Method for external field fast calibrating miniature multi-sensor combined navigation system
CN101685102A (en) * 2009-04-17 2010-03-31 幻音科技(深圳)有限公司 Device and method for regulating precision of triaxial accelerometer
CN101887068A (en) * 2010-06-01 2010-11-17 中北大学 Calibration compensation method for triaxial vector sensor and biaxial vector sensor

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7644602B2 (en) * 2004-03-31 2010-01-12 National Institute Of Advanced Industrial Science And Technology Method of measuring transverse sensitivity of sensor for detecting acceleration and acceleration measuring method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101246024A (en) * 2008-03-26 2008-08-20 北京航空航天大学 Method for external field fast calibrating miniature multi-sensor combined navigation system
CN101685102A (en) * 2009-04-17 2010-03-31 幻音科技(深圳)有限公司 Device and method for regulating precision of triaxial accelerometer
CN101887068A (en) * 2010-06-01 2010-11-17 中北大学 Calibration compensation method for triaxial vector sensor and biaxial vector sensor

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于有限单元法三轴体微加速度传感器固有频率的灵敏度分析;吴振亭等;《功能材料与器件学报》;20100831;第16卷(第4期);第399-402页 *

Also Published As

Publication number Publication date
CN102654515A (en) 2012-09-05

Similar Documents

Publication Publication Date Title
US10605619B2 (en) Precision calibration method of attitude measuring system
CN104884902B (en) For three axle magnetometer and the method and apparatus of the data fusion of three axis accelerometer
CN102654515B (en) Calibration algorithm for z sensitive shaft of three-shaft acceleration transducer
US20110301897A1 (en) Calibration of a triaxial magnetic field sensor
CN105387859B (en) MEMS sensor combination temp drift error compensation method
CN102168991B (en) Calibration and compensation method for mounting errors between triaxial vector sensor and mounting carrier
CN108195404B (en) Calibration method for zero drift of gyroscope
CN109550219B (en) Method and system for determining motion information and mobile device
CN103808349B (en) The error calibration method of vector sensor and device
CN103196462A (en) Compensation method for error calibration of MEMS gyroscope in MIMU
CN108534744A (en) A kind of attitude angle acquisition methods, device and handle
CN103823084A (en) Method for calibrating three-axis acceleration sensor
KR20140025319A (en) Apparatuses and methods for dynamic tracking and compensation of magnetic near field
CN111189474A (en) Autonomous calibration method of MARG sensor based on MEMS
CN116147624B (en) Ship motion attitude calculation method based on low-cost MEMS navigation attitude reference system
CN106705995A (en) Calibration method of MEMS gyroscope g value sensitive coefficient
CN103913179A (en) Ship-based single-axial rotation modulation strapdown inertial navigation attitude angle correction method
CN111982155B (en) Calibration method and device of magnetic sensor, electronic equipment and computer storage medium
CN105277210A (en) Calibration method for multi-axis integrated gyro installed in any way
JP2012037405A (en) Sensor device, electronic apparatus, and offset correction method of angular velocity sensor
CN102042833B (en) Motion tracing method and system
CN110631605B (en) Gyro array calibration method and system
CN102322844B (en) Pointing correction method of telescope based on planar interpolating model
CN110657821B (en) Calibration method of track detection inertia measurement device
CN107389092A (en) A kind of Gyro Calibration method based on Magnetic Sensor auxiliary

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160923

Address after: 214000 Wuxi national high tech Industrial Development Zone, Hui Hui Road, No. 2, No.

Patentee after: Meixin Semi-conductor (Wuxi) Co., Ltd.

Address before: 214101 Dacheng Road, Xishan Economic Development Zone (East District), Nanan, Jiangsu, Wuxi, Sheng Ludong

Patentee before: Memsic Transducer Systems Co., Ltd.