CN108732380A - 一种三角函数温度补偿方法 - Google Patents
一种三角函数温度补偿方法 Download PDFInfo
- Publication number
- CN108732380A CN108732380A CN201810560713.1A CN201810560713A CN108732380A CN 108732380 A CN108732380 A CN 108732380A CN 201810560713 A CN201810560713 A CN 201810560713A CN 108732380 A CN108732380 A CN 108732380A
- Authority
- CN
- China
- Prior art keywords
- trigonometric function
- temperature
- fitting
- data
- coefficient
- 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
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P15/00—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration
- G01P15/02—Measuring acceleration; Measuring deceleration; Measuring shock, i.e. sudden change of acceleration by making use of inertia forces using solid seismic masses
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
- G06F7/38—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation
- G06F7/48—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices
- G06F7/544—Methods or arrangements for performing computations using exclusively denominational number representation, e.g. using binary, ternary, decimal representation using non-contact-making devices, e.g. tube, solid state device; using unspecified devices for evaluating functions by calculation
- G06F7/548—Trigonometric functions; Co-ordinate transformations
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
Abstract
本发明公开了一种三角函数温度补偿方法,所述方法包括:在温箱中对需要温度补偿的***进行温循环实验;将实验数据进行预处理,N组待拟合数据;用三角函数对N组待拟合数据点进行拟合,得到三角函数表达式:利用三角函数表达式对原始数据进行补偿得到补偿结果;实现了能够达到高阶数的多项式拟合效果,也能弥补双指数函数对复杂曲线拟合效果较差问题,同时计算量小的技术效果。
Description
技术领域
本发明涉及温度补偿方法领域,具体地,涉及一种三角函数温度补偿方法。
背景技术
随着MEMS技术的发展,惯性传感器在过去的几年中成为应用最广泛的微机电***器件之一,其中微加速度计已经成为被广泛用于测量物体加速度的惯性装置。但是环境温度的变化是MEMS加速度计产生误差的最重要原因之一,而且其导致MEMS加速度计的性能大幅度的降低。
现有的温度补偿软件方法大多是采用多项式函数:
如y=a0+a1×x+a2×x2+a3×x3......,由于多项式函数计算简单,可以采用线性最小二乘法拟合因此被广泛应用。然而多项式函数也存在一定的局限性,由于温度曲线具有非线性的特点,当面对较为复杂的温度曲线时,需要采用较高的阶次。因此参数较多,计算量增加。神经网络和卡尔曼滤波等算法也可以有效地对温度曲线进行补偿,不过算法比较复杂,仅限于理论研究,很难集成到传感器芯片中。
发明内容
本发明提供了一种三角函数温度补偿方法,解决了现有的温度补偿方法计算量大,计算复杂的技术问题,实现了能够达到高阶数的多项式拟合效果,也能弥补双指数函数对复杂曲线拟合效果较差问题,同时计算量小的技术效果。
为实现上述发明目的,本申请提供了一种三角函数温度补偿方法,所述方法包括:
在温箱中对需要温度补偿的***进行温循环实验;
将实验数据进行预处理,N组待拟合数据;
用三角函数对N组待拟合数据点进行拟合,得到三角函数表达式:
y=a0+a1×sin(x)+a2×cos(x);其中,N为大于等于2的正整数,a0为线性系数;a1为正弦系数;a2为余弦系数;x为函数自变量即对应的温度值;
利用三角函数表达式对原始数据进行补偿得到补偿结果。
进一步的,本方法采用三角函数最小二乘法对N组待拟合数据点进行拟合,具体包括:
对原始数据进行预处理,获得散点图(Xi Yi),Xi为温度电压值,Yi为输出电压值;
确定拟合函数表达式:y=a0+a1×sin(x)+a2×cos(x);
采用线性最小二乘法计算系数;
1.现在需要构造一组系数,使其满足与原始数据点的偏差程度最小,即
2.对函数进行替换,如下
3.的I为系数a0,a1,...,an的多元函数,所以上述问题可以转换为求解多元函数I极值的问题。根据多元函数求解极值的必要条件,对每一个系数求偏导数,使其等于0可得
4.化简后可得法方程组
将数据(Xi Yi)带入法方程组,求解法方程组,得到系数和函数表达式。
进一步的,对数据进行预处理包括:利用小波分析进行降噪处理、平滑处理、等间隔采样剔除奇异值。
本申请提供的一个或多个技术方案,至少具有如下技术效果或优点:
本申请中的三角函数的泰勒展开式具有无穷阶数,在应对非线性的温度曲线时,能够达到高阶数的多项式拟合效果,也能弥补双指数函数对复杂曲线拟合效果较差问题。
本方法中的表达式系数比较少,只有a0,a1,a2三个系数,采用最小二乘法进行计算,计算量较小。
如果需要设计出易集成到MEMS加速度计中的补偿电路,可以使用CORDIC算法或查找表等方法设计补偿电路模块,仅需要加减法迭代就可以计算三角函数,硬件资源少。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定;
图1是本申请中三角函数温度补偿方法框架示意图;
图2是本申请中三角函数最小二乘法流程示意图;
图3是本申请中三角函数拟合效果示意图;
图4是本申请中温度补偿效果示意图。
具体实施方式
本发明提供了一种三角函数温度补偿方法,解决了现有的温度补偿方法计算量大,计算复杂的技术问题,实现了能够达到高阶数的多项式拟合效果,也能弥补双指数函数对复杂曲线拟合效果较差问题,同时计算量小的技术效果。
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在相互不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述范围内的其他方式来实施,因此,本发明的保护范围并不受下面公开的具体实施例的限制。
一种三角函数温度补偿方法,流程框图如图1所示,具体流程如下:
在温箱中对需要温度补偿的***进行温循环实验;
将实验数据进行预处理,包括利用小波分析、平滑处理、等间隔采样;
用三角函数对N组数据点进行拟合,得到三角函数表达式;
利用三角函数表达式对原始数据进行补偿得到补偿结果。
三角函数拟合算法流程图如图2所示。
算法包括三角函数模型y=a0+a1×sin(x)+a2×cos(x)。
拟合系数为a0,a1,a2。
补偿算法主要采用最小二乘法,首先得到预处理后的N组待拟合数据(Xi,Yi)(i=1,2,3,4……N);其次确定拟合函数表达式即y=a0+a1×sin(x)+a2×cos(x);接下来计算得到关于系数的线性方程组,即法方程组;最后将N组数据(Xi,Yi)(i=1,2,3,4……N)带入法方程组解出系数得到函数表达式。
具体的实施方式可以采用一种补偿***实现:
将含有温度传感器的MEMS传感器进行温循环实验,得到在一定温度范围内的温度-输出数据。
将数据采用MATLAB进行处理,并进行等间隔得到采样得到需要拟合的数据点。
确定三角函数并采用最小二乘法进行拟合,得到拟合系数。
设计一款集成到MEMS传感器中的补偿电路,三角函数可以采用CORDIC算法或查找表等方法实现,补偿电路中采用rom存储系数,直接对输出进行补偿即可。
本发明采用利用三角函数的函数特点对MEMS加速度计的温度曲线进行补偿,由于最终补偿结果(如下图3,4所示),拟合精度达到RSME=4.7024E-05、Rsquare=0.9947,补偿后温度曲线峰峰值减少一个数量级以上。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (3)
1.一种三角函数温度补偿方法,其特征在于,所述方法包括:
在温箱中对需要温度补偿的***进行温循环实验;
将实验数据进行预处理,N组待拟合数据;
用三角函数对N组待拟合数据点进行拟合,得到三角函数表达式:
y=a0+a1×sin(x)+a2×cos(x);其中,N为大于等于2的正整数,a0为线性系数;a1为正弦系数;a2为余弦系数;x为函数自变量即对应的温度值;
利用三角函数表达式对原始数据进行补偿得到补偿结果。
2.根据权利要求1所述的三角函数温度补偿方法,其特征在于,本方法采用三角函数最小二乘法对N组待拟合数据点进行拟合,具体包括:
对原始数据进行预处理,获得散点图(Xi Yi),Xi为温度电压值,Yi为输出电压值;
确定拟合函数表达式:y=a0+a1×sin(x)+a2×cos(x);
计算得到关于系数的线性方程组;
将数据(Xi Yi)带入线性方程组,求解线性方程组,得到系数和函数表达式。
3.根据权利要求1或2所述的三角函数温度补偿方法,其特征在于,对数据进行预处理包括:利用小波分析、平滑处理、等间隔采样。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810560713.1A CN108732380A (zh) | 2018-05-25 | 2018-05-25 | 一种三角函数温度补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810560713.1A CN108732380A (zh) | 2018-05-25 | 2018-05-25 | 一种三角函数温度补偿方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108732380A true CN108732380A (zh) | 2018-11-02 |
Family
ID=63931545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810560713.1A Pending CN108732380A (zh) | 2018-05-25 | 2018-05-25 | 一种三角函数温度补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108732380A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111271810A (zh) * | 2020-01-08 | 2020-06-12 | 宁波奥克斯电气股份有限公司 | 一种空调控制方法、控制***、空调器及计算机存储介质 |
CN113063964A (zh) * | 2021-03-23 | 2021-07-02 | 西安微电子技术研究所 | 一种温度补偿型石英挠性加速度计伺服电路及石英挠性加速度计 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101135723A (zh) * | 2007-07-26 | 2008-03-05 | 北京航空航天大学 | 通用电力互感器校验装置 |
CN101142049A (zh) * | 2005-02-15 | 2008-03-12 | 电子科学工业公司 | 修正激光处理***中的***误差的方法 |
CN101718584A (zh) * | 2008-10-10 | 2010-06-02 | 电子科技大学 | 一种基于多项式拟合的全光采样的线性度和转换效率计算方法 |
CN101980453A (zh) * | 2010-10-16 | 2011-02-23 | 中国科学院上海微***与信息技术研究所 | 多普勒频率估计与补偿方法及*** |
CN102519471A (zh) * | 2011-12-08 | 2012-06-27 | 北京控制工程研究所 | 基于三角函数拟合的成像型地球敏感器地球扁率补偿方法 |
CN103453898A (zh) * | 2013-09-09 | 2013-12-18 | 无锡慧联信息科技有限公司 | 在恶劣环境下提高光纤陀螺寻北精度的方法 |
CN104467816A (zh) * | 2014-12-29 | 2015-03-25 | 电子科技大学 | 一种晶体振荡器的温度补偿*** |
CN104535091A (zh) * | 2014-12-29 | 2015-04-22 | 华中科技大学 | 光纤传感单元及其在同时测量折射率温度的应用 |
CN105868475A (zh) * | 2016-03-30 | 2016-08-17 | 中国科学院国家天文台 | 一种望远镜反射面的准实时辅助控制方法及*** |
CN106406830A (zh) * | 2015-07-29 | 2017-02-15 | 腾讯科技(深圳)有限公司 | 一种准周期信号的预测方法和装置 |
CN106603068A (zh) * | 2016-12-09 | 2017-04-26 | 电子科技大学 | 一种晶体振荡器的温度补偿*** |
-
2018
- 2018-05-25 CN CN201810560713.1A patent/CN108732380A/zh active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101142049A (zh) * | 2005-02-15 | 2008-03-12 | 电子科学工业公司 | 修正激光处理***中的***误差的方法 |
CN101135723A (zh) * | 2007-07-26 | 2008-03-05 | 北京航空航天大学 | 通用电力互感器校验装置 |
CN101718584A (zh) * | 2008-10-10 | 2010-06-02 | 电子科技大学 | 一种基于多项式拟合的全光采样的线性度和转换效率计算方法 |
CN101980453A (zh) * | 2010-10-16 | 2011-02-23 | 中国科学院上海微***与信息技术研究所 | 多普勒频率估计与补偿方法及*** |
CN102519471A (zh) * | 2011-12-08 | 2012-06-27 | 北京控制工程研究所 | 基于三角函数拟合的成像型地球敏感器地球扁率补偿方法 |
CN103453898A (zh) * | 2013-09-09 | 2013-12-18 | 无锡慧联信息科技有限公司 | 在恶劣环境下提高光纤陀螺寻北精度的方法 |
CN104467816A (zh) * | 2014-12-29 | 2015-03-25 | 电子科技大学 | 一种晶体振荡器的温度补偿*** |
CN104535091A (zh) * | 2014-12-29 | 2015-04-22 | 华中科技大学 | 光纤传感单元及其在同时测量折射率温度的应用 |
CN106406830A (zh) * | 2015-07-29 | 2017-02-15 | 腾讯科技(深圳)有限公司 | 一种准周期信号的预测方法和装置 |
CN105868475A (zh) * | 2016-03-30 | 2016-08-17 | 中国科学院国家天文台 | 一种望远镜反射面的准实时辅助控制方法及*** |
CN106603068A (zh) * | 2016-12-09 | 2017-04-26 | 电子科技大学 | 一种晶体振荡器的温度补偿*** |
Non-Patent Citations (3)
Title |
---|
冯林 等: ""标准击实试验的曲线拟合方法"", 《路基工程》 * |
刘一兵 等: ""一种基于传感器温度补偿的双指数函数拟合算法"", 《电子应用技术》 * |
朱冬冬 等: ""基于硬件算法库的三角函数二次拟合查表法"", 《中国仪器仪表》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111271810A (zh) * | 2020-01-08 | 2020-06-12 | 宁波奥克斯电气股份有限公司 | 一种空调控制方法、控制***、空调器及计算机存储介质 |
CN113063964A (zh) * | 2021-03-23 | 2021-07-02 | 西安微电子技术研究所 | 一种温度补偿型石英挠性加速度计伺服电路及石英挠性加速度计 |
CN113063964B (zh) * | 2021-03-23 | 2023-07-14 | 西安微电子技术研究所 | 一种温度补偿型石英挠性加速度计伺服电路及石英挠性加速度计 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108732380A (zh) | 一种三角函数温度补偿方法 | |
CN108051018A (zh) | 基于fpga的分布式光纤光栅传感信号寻峰控制*** | |
Barczy et al. | Explicit formulas for Laplace transforms of certain functionals of some time inhomogeneous diffusions | |
CN107645493B (zh) | 一种ip组相似度计算方法 | |
CN109001489A (zh) | 一种双指数函数温度补偿方法 | |
CN114595946A (zh) | 一种海域理论最低潮面计算方法、***、设备及介质 | |
CN105357335B (zh) | 一种dns权威日志信息挖掘处理方法 | |
Abodayeh et al. | Generalized Ω-distance mappings and some fixed point theorems | |
CN115060306B (zh) | 模拟量校准方法、装置及模拟量采集设备 | |
Zhang et al. | Spectral analysis based on fast Fourier transformation (FFT) of surveillance data: the case of scarlet fever in China | |
Tan et al. | Gaussian recursive filter for nonlinear systems with finite-step correlated noises and packet dropout compensations | |
Zheng et al. | An Enhanced Differential Private Protection Method Based on Adaptive Iterative Wiener Filtering in Discrete Time Series | |
Jayakumar | Generalized Simpson-Newton's Method for Solving Nonlinear Equations with Cubic Convergence | |
CN110596425B (zh) | 一种无人机mems加速度传感器噪声消除方法 | |
Xu et al. | Fpga implementation of a best-precision fixed-point digital pid controller | |
Li et al. | A best approximation for the solution of one‐dimensional variable‐coefficient Burgers' equation | |
Agarwal et al. | Boundary Value Problems | |
CN114254487A (zh) | 一种光纤陀螺温度漂移建模及实时补偿方法 | |
CN106777834B (zh) | 一种高精度的传递函数采样方法 | |
Hao et al. | Adaptive weighted data fusion of muti-sensor based on fuzzy preference relations | |
CN108875142A (zh) | 一种无限阶数易实现温度补偿*** | |
Li et al. | Uniform convergence analysis for singularly perturbed elliptic problems with parabolic layers | |
Xu et al. | Frequency detection of self-adaption control based on chaotic theory | |
Geng et al. | Analytical approximation to solutions of singularly perturbed boundary value problems. | |
CN113687345B (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 |
Application publication date: 20181102 |
|
RJ01 | Rejection of invention patent application after publication |