CN106706153A - 一种低温热敏电阻温度特性数据的处理方法 - Google Patents
一种低温热敏电阻温度特性数据的处理方法 Download PDFInfo
- Publication number
- CN106706153A CN106706153A CN201710160919.0A CN201710160919A CN106706153A CN 106706153 A CN106706153 A CN 106706153A CN 201710160919 A CN201710160919 A CN 201710160919A CN 106706153 A CN106706153 A CN 106706153A
- Authority
- CN
- China
- Prior art keywords
- temperature
- resistance
- thermistor
- polynomial
- fitting
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
- G01K7/16—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using resistive elements
- G01K7/22—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using resistive elements the element being a non-linear resistance, e.g. thermistor
Landscapes
- Physics & Mathematics (AREA)
- Nonlinear Science (AREA)
- General Physics & Mathematics (AREA)
- Thermistors And Varistors (AREA)
- Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
Abstract
本发明涉及一种低温热敏电阻温度特性数据的处理方法,步骤一:根据已获得的温度传感器分度表中电阻R及温度数据T,幂函数对已获得的电阻和温度进行变量代换,将变量代换后电阻和温度记为R’和T’;步骤二:取变量代换后的数据R’和T’,利用最小二乘法拟合其函数关系,得到拟合的多项式R*(T′);步骤三:根据步骤一二的结果,求出拟合相关系数r;步骤四,判断拟合相关系数r是否满足热敏电阻使用时的精度要求,若满足,则步骤二中的多项式即可表示该电阻的温度特性曲线;否则,调整或者重新拟合电阻变换值多项式R*(T′),从步骤三开始执行。
Description
技术领域
本发明属于温度测量领域,尤其涉及一种低温热敏电阻温度特性数据的处理方法。
背景技术
在工业生产和自动化控制系统中,温度参数的测量是普遍而重要的,常用的温度传感器按照传感器材料及电子元件特性可分为热电偶、热敏电阻等,通过材料或元件感受温度并转换为可用输出信号。温度传感器的输出一般为模拟量的电信号,需要通过变送装置将其显示为温度数值。
热敏电阻的重要性质之一就是电阻温度系数,电阻温度系数的定义是温度变化导致的电阻随之变化的相对速率,电阻温度系数越大,热敏电阻对温度变化的反应越灵敏。按照电阻温度系数的不同性质,热敏电阻温度传感器还分为正温度系数和负温度系数两种类型,正温度系数是指在电阻随温度升高而升高,负温度系数是指电阻随温度升高而降低。电阻温度系数TCR的定义式为:
式中,T为温度;R(T)为热敏电阻阻值,与温度有关。
目前常用的温度传感器是通过将传感器分度表输入数据采集系统,利用插值的方法求出待测温度。分度表是指传感器温度与相关电信号之间的对应关系表格,也可以利用分度表来求出热敏电阻的温度系数。
发明内容
本发明的目的提供一种能够提高温度显示精度的低温热敏电阻温度特性数据处理方法。
本发明的上述目的主要是通过如下技术方案予以实现的:一种低温热敏电阻温度特性数据的处理方法,包括以下步骤:
步骤一:根据已获得的温度传感器分度表中电阻R及温度数据T,幂函数对已获得的电阻和温度进行变量代换,将变量代换后电阻和温度记为R’和T’;
步骤二:取变量代换后的数据R’和T’,利用最小二乘法拟合其函数关系,得到拟合的多项式R*(T′);
步骤三:根据步骤一二的结果,求出拟合相关系数r,其公式为
式中,N为电阻温度特性曲线中二维数据组的个数,R′为电阻变换值的平均值;
步骤四,判断拟合相关系数r是否满足热敏电阻使用时的精度要求,若满足,则步骤二中的多项式即可表示该电阻的温度特性曲线;否则,调整或者重新拟合电阻变换值多项式R*(T′),从步骤三开始执行。
进一步的,步骤一中变量代换的函数形式如下:
式中,幂级n为一正实数。
进一步的,步骤二中选用的拟合多项式形式为
式中,ai为拟合系数,m为多项式次数。
进一步的,步骤四中通过调整多项式的次数或者幂函数的幂级,以适应曲线拟合的不同精度要求。
进一步的,r值越接近于1,则温度特性曲线的显函数表示方法精度越高。
进一步的,电阻温度特性曲线的函数表达
或者写成
式中,ai、bi为拟合系数,m为多项式次数。
进一步的,将步骤四中确定的多项式进行导数计算,确定电阻温度系数TCR。
进一步的,利用步骤四中确定的多项式表示热敏电阻阻值与其温度的一一对应关系,将热敏电阻加工制成温度传感器后,利用多项式函数关系计算得到温度数值。
本发明与现有技术相比具有如下有益效果:
1、基于最小二乘法的多项式拟合,能够利用简单函数关系表示热敏电阻阻值与温度的对应关系,将温度传感器的输出电信号迅速变送显示为温度数值,有利于工业计算机及其他测量控制设备的读取与处理。
2、取得了低温热敏电阻温度特性曲线的显函数表达式,可以通过微分计算的方法得到有效测温区间内任一温度点时热敏电阻的温度系数解析表达式,与常用的查表、插值方法相比,解析表达式计算得出的温度系数精度更高。
3、将分度表所描述的温度特性曲线利用显函数表示,可以提高传感器的输出精度,真实地反映温度变化,简化设计计算,便于数据的后续处理。
附图说明
图1为基于最小二乘法的低温热敏电阻温度特性数据处理方法流程图;
图2为本实施例中提供的热敏电阻温度特性曲线;
图3为本实施例中经过变量代换的电阻温度特性曲线;
图4为本实施例中拟合得到的电阻温度特性曲线。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
以下结合附图和具体实施例对本发明提出的一种低温热敏电阻温度特性数据的处理方法作进一步详细说明。根据下面说明和权利要求书,本发明的优点和特征将更清楚。需说明的是,附图均采用非常简化的形式且均使用非精准的比率,仅用于方便、明晰地辅助说明本发明实施例的目的。
本发明的核心思想在于,本发明提供的低温热敏电阻温度特性数据处理方法,在进行拟合之前,通过变量代换提高了拟合精度,该方法操作简便,保证了测温元件电阻值与对应温度的线性度,进而保障了测温精度,并能够获得可靠的电阻温度系数。
本发明实施例提供的一种低温热敏电阻温度特性数据的处理方法实现如图1所示,步骤如下:
(1)将某薄膜型热敏电阻温度特性数据以分度表的形式录入计算机,如下表所示;
电阻R/Ω | 温度T/K |
2150 | 5.0353 |
1500 | 7.2465 |
1000 | 12.0248 |
465 | 35.2828 |
292 | 65.4305 |
216 | 95.4198 |
143 | 155.3876 |
107 | 215.9580 |
86 | 275.3754 |
74 | 325.0130 |
(2)将上述数据在二维坐标轴中绘制出来,如图2所示;
(3)进行变量代换,令再次将上述数据作图如图3所示;
(4)使用计算机软件,利用最小二乘法原理,进行4次多项式拟合,经过计算,得到的拟合相关系数r为0.999984,满足精度要求,因此可以根据变量代换关系的反函数,求出电阻温度特性曲线的函数关系为:
该拟合函数与原数据的关系见图4,从图中可见,利用拟合函数,可以避免在分度表上插值读数所引起的误差,可提高热敏电阻处于某一阻值时所显示温度的精度,同时,电阻温度曲线还可以表示为:
此时可以方便地利用导数计算,求出该电阻的温度系数TCR。
以上的实施例说明仅为本发明的较佳实施例说明,本领域技术人员可依据本发明的上述实施例说明而做出其它种种等效的替换及修改。然而这些依据本发明实施例所作的种种等效替换及修改,属于本发明的发明精神及由权利要求所界定的保护范围内。
Claims (8)
1.一种低温热敏电阻温度特性数据的处理方法,其特征在于包括以下步骤:
步骤一:根据已获得的温度传感器分度表中电阻R及温度数据T,幂函数对已获得的电阻和温度进行变量代换,将变量代换后电阻和温度记为R’和T’;
步骤二:取变量代换后的数据R’和T’,利用最小二乘法拟合其函数关系,得到拟合的多项式R*(T′);
步骤三:根据步骤一二的结果,求出拟合相关系数r,其公式为
式中,N为电阻温度特性曲线中二维数据组的个数,为电阻变换值的平均值;
步骤四,判断拟合相关系数r是否满足热敏电阻使用时的精度要求,若满足,则步骤二中的多项式即可表示该电阻的温度特性曲线;否则,调整或者重新拟合电阻变换值多项式R*(T′),从步骤三开始执行。
2.根据权利要求1所述的方法,其特征在于:步骤一中变量代换的函数形式如下:
式中,幂级n为一正实数。
3.根据权利要求1所述的方法,其特征在于:步骤二中选用的拟合多项式形式为
式中,ai为拟合系数,m为多项式次数。
4.根据权利要求1或2所述的方法,其特征在于:步骤四中通过调整多项式的次数或者幂函数的幂级,以适应曲线拟合的不同精度要求。
5.根据权利要求1所述的方法,其特征在于:r值越接近于1,则温度特性曲线的显函数表示方法精度越高。
6.根据权利要求1所述的方法,其特征在于:电阻温度特性曲线的函数表达
或者写成
式中,ai、bi为拟合系数,m为多项式次数。
7.根据权利要求1所述的方法,其特征在于:将步骤四中确定的多项式进行导数计算,确定电阻温度系数TCR。
8.根据权利要求1所述的方法,其特征在于:利用步骤四中确定的多项式表示热敏电阻阻值与其温度的一一对应关系,将热敏电阻加工制成温度传感器后,利用多项式函数关系计算得到温度数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710160919.0A CN106706153B (zh) | 2017-03-17 | 2017-03-17 | 一种低温热敏电阻温度特性数据的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710160919.0A CN106706153B (zh) | 2017-03-17 | 2017-03-17 | 一种低温热敏电阻温度特性数据的处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106706153A true CN106706153A (zh) | 2017-05-24 |
CN106706153B CN106706153B (zh) | 2019-09-06 |
Family
ID=58887253
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710160919.0A Active CN106706153B (zh) | 2017-03-17 | 2017-03-17 | 一种低温热敏电阻温度特性数据的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106706153B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107462342A (zh) * | 2017-07-31 | 2017-12-12 | 上海辛格林纳新时达电机有限公司 | 热电偶的温度获取方法、电子设备和计算机可读存储介质 |
CN107505843A (zh) * | 2017-09-15 | 2017-12-22 | 中国科学院长春光学精密机械与物理研究所 | 一种空间光学有效载荷的主动热控优化方法 |
CN108225594A (zh) * | 2017-12-25 | 2018-06-29 | 青岛海尔科技有限公司 | 获取ntc电阻温度的方法、装置、设备及可读存储介质 |
CN110823405A (zh) * | 2019-10-17 | 2020-02-21 | 五邑大学 | 一种温度传感器的温度标定方法 |
CN111143911A (zh) * | 2019-12-12 | 2020-05-12 | 上海卫星工程研究所 | 卫星温度遥测的数据处理方法及系统 |
CN111552275A (zh) * | 2020-05-20 | 2020-08-18 | 深圳开立生物医疗科技股份有限公司 | 一种温控校准装置、温控设备及其方法 |
CN111623898A (zh) * | 2020-06-24 | 2020-09-04 | 广东纽恩泰新能源科技发展有限公司 | Ntc热敏电阻温度与电阻检测方法以及温度传感器检测方法 |
CN112268633A (zh) * | 2020-10-22 | 2021-01-26 | 合肥中科离子医学技术装备有限公司 | 用于30至73k内工业铂电阻的温度与电阻值转换方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002257641A (ja) * | 2001-03-02 | 2002-09-11 | Ricoh Co Ltd | 温度履歴表示材料 |
CN2816948Y (zh) * | 2005-09-05 | 2006-09-13 | 张卫社 | 基于串联热敏特性的模拟量线型感温探测线缆 |
CN101922981A (zh) * | 2010-09-10 | 2010-12-22 | 中冶南方(武汉)自动化有限公司 | 一种基于dsp的热敏电阻测温方法 |
CN104048777A (zh) * | 2014-06-24 | 2014-09-17 | 中国东方电气集团有限公司 | 用于ntc型热敏电阻的温度与电阻值转换方法 |
CN104697659A (zh) * | 2015-03-04 | 2015-06-10 | 徐州工程学院 | 一种采用ntc热敏电阻的温度测量装置与方法 |
-
2017
- 2017-03-17 CN CN201710160919.0A patent/CN106706153B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002257641A (ja) * | 2001-03-02 | 2002-09-11 | Ricoh Co Ltd | 温度履歴表示材料 |
CN2816948Y (zh) * | 2005-09-05 | 2006-09-13 | 张卫社 | 基于串联热敏特性的模拟量线型感温探测线缆 |
CN101922981A (zh) * | 2010-09-10 | 2010-12-22 | 中冶南方(武汉)自动化有限公司 | 一种基于dsp的热敏电阻测温方法 |
CN104048777A (zh) * | 2014-06-24 | 2014-09-17 | 中国东方电气集团有限公司 | 用于ntc型热敏电阻的温度与电阻值转换方法 |
CN104697659A (zh) * | 2015-03-04 | 2015-06-10 | 徐州工程学院 | 一种采用ntc热敏电阻的温度测量装置与方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107462342A (zh) * | 2017-07-31 | 2017-12-12 | 上海辛格林纳新时达电机有限公司 | 热电偶的温度获取方法、电子设备和计算机可读存储介质 |
CN107505843A (zh) * | 2017-09-15 | 2017-12-22 | 中国科学院长春光学精密机械与物理研究所 | 一种空间光学有效载荷的主动热控优化方法 |
CN107505843B (zh) * | 2017-09-15 | 2020-05-15 | 中国科学院长春光学精密机械与物理研究所 | 一种空间光学有效载荷的主动热控优化方法 |
CN108225594A (zh) * | 2017-12-25 | 2018-06-29 | 青岛海尔科技有限公司 | 获取ntc电阻温度的方法、装置、设备及可读存储介质 |
CN110823405A (zh) * | 2019-10-17 | 2020-02-21 | 五邑大学 | 一种温度传感器的温度标定方法 |
CN110823405B (zh) * | 2019-10-17 | 2020-12-01 | 五邑大学 | 一种温度传感器的温度标定方法 |
CN111143911A (zh) * | 2019-12-12 | 2020-05-12 | 上海卫星工程研究所 | 卫星温度遥测的数据处理方法及系统 |
CN111552275A (zh) * | 2020-05-20 | 2020-08-18 | 深圳开立生物医疗科技股份有限公司 | 一种温控校准装置、温控设备及其方法 |
CN111623898A (zh) * | 2020-06-24 | 2020-09-04 | 广东纽恩泰新能源科技发展有限公司 | Ntc热敏电阻温度与电阻检测方法以及温度传感器检测方法 |
CN112268633A (zh) * | 2020-10-22 | 2021-01-26 | 合肥中科离子医学技术装备有限公司 | 用于30至73k内工业铂电阻的温度与电阻值转换方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106706153B (zh) | 2019-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106706153A (zh) | 一种低温热敏电阻温度特性数据的处理方法 | |
CN103946672B (zh) | 物理量传感器或其装置的输出值修正方法及装置 | |
CN100533087C (zh) | 高精度压力传感器的误差补偿方法 | |
CN103575414B (zh) | 一种热电偶温度测量误差补偿方法 | |
CN101706346B (zh) | 智能力传感器测量非线性温度漂移的温度补偿方法 | |
CN104048777B (zh) | 用于ntc型热敏电阻的温度与电阻值转换方法 | |
CN103961066B (zh) | 一种测温方法及测温装置 | |
CN104283479B (zh) | 基于三维电机参数表的永磁同步电机负载转矩监测系统 | |
CN206601203U (zh) | 温度检测电路 | |
CN104697659A (zh) | 一种采用ntc热敏电阻的温度测量装置与方法 | |
CN106092358A (zh) | 一种电阻式传感器信号的处理方法和装置 | |
CN103063233B (zh) | 一种采用多传感器降低测量误差的方法 | |
CN104062025A (zh) | 一种节能自动ad温度采集监控系统 | |
CN205102947U (zh) | 一种三线式铂热电阻pt100测温电路 | |
CN112946554B (zh) | 一种电流传感器的标定方法 | |
CN116296047B (zh) | 一种单晶硅压力变送器的温度补偿改进方法 | |
CN203133019U (zh) | 电化学气体传感器的干扰补偿装置 | |
CN108875218B (zh) | 一种基于sinc函数的变压器涡流场有限元求解方法 | |
CN1308225A (zh) | 基于新测温原理的快速预测型温度计 | |
CN103925868A (zh) | 一种没有长导线电阻影响的电阻应变仪 | |
CN204649327U (zh) | 一种热电阻信号转换电路 | |
Dong et al. | Curve fitting and piecewise linear method for z-type temperature sensor | |
CN105092914A (zh) | 数字式电量变送器及仪表的温漂补偿算法 | |
CN105354385A (zh) | 一种半导体工艺仿真中的网格排序实现方法 | |
CN104197975B (zh) | 一种基于测量值微分约束的传感器测量精度提高方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220803 Address after: 100074, Lane 1, Yungang Tian Town, Beijing, Fengtai District Patentee after: BEIJING INSTITUTE OF AEROSPACE TESTING TECHNOLOGY Address before: 100074, No. 1, Zhong Li, Tian City, Beijing, Fengtai District Patentee before: BEIJING INSTITUTE OF AEROSPACE TESTING TECHNOLOGY Patentee before: BEIJING AEROSPACE RATE MECHANICAL & ELECTRICAL ENGINEERING Co.,Ltd. |
|
TR01 | Transfer of patent right |