CN106679848A - 一种利用普通温度传感器快速测温的方法 - Google Patents
一种利用普通温度传感器快速测温的方法 Download PDFInfo
- Publication number
- CN106679848A CN106679848A CN201710188091.XA CN201710188091A CN106679848A CN 106679848 A CN106679848 A CN 106679848A CN 201710188091 A CN201710188091 A CN 201710188091A CN 106679848 A CN106679848 A CN 106679848A
- Authority
- CN
- China
- Prior art keywords
- temperature
- coefficient
- time
- polynomial
- sensor
- 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
- G01K13/00—Thermometers specially adapted for specific purposes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
Abstract
本发明涉及工业测温技术领域,具体涉及一种利用普通温度传感器快速测温的方法,包括以下步骤:S1、建立温度‑时间响应模型,得出温度传感器输出温度与时间的关系曲线,用多项式拟合关系曲线,得到由多项式的系数构成的系数矩阵;S2、将温度传感器放入被测环境,测量得到多个时刻下温度传感器的实测温度,多个时刻均远小于热平衡时间;S3、根据系数矩阵、多项式,以及步骤S2中的时刻、实测温度通过计算方法得出最终温度。本发明在正式测温时仅测量远小于热平衡时间的多个时刻的实测温度,再结合系数矩阵和多项式即可计算得出待测温度,整个测量和计算过程所需的时间远小于热平衡时间,实现了利用普通温度传感器的快速测温。
Description
技术领域
本发明属于工业测温技术领域,具体涉及一种利用普通温度传感器快速测温的方法。
背景技术
在工业应用中,大量使用温度传感器测量制造工艺过程中的实时温度,用温度数据来监控制造过程执行情况或作为改造工艺参数的重要依据。在温度测量过程中,温度传感器的热平衡过程总需要一定的时间,我们在此将传感器开始工作到热平衡状态所消耗的时间称为热平衡时间。严格来说,从测温开始必须等待大于热平衡时间之后得到的温度才是被测件的准确温度。因此在实际生产过程中,使用普通温度传感器的企业往往需要在测温环节等待较长时间,或者抛弃普通传感器转而寻求其他成本高昂的测温方式。
发明内容
本发明的目的在于提供一种利用普通温度传感器快速测温的方法,该方法可以很好地解决采用普通温度传感器测温耗时较长的问题。
为达到上述要求,本发明采取的技术方案是:提供一种利用普通温度传感器快速测温的方法,包括以下步骤:
S1、建立温度-时间响应模型,通过实验方法得到温度传感器在不同条件下输出温度与时间的关系曲线,采用多项式拟合所述关系曲线,得到由多项式系数构成的系数矩阵;
S2、将温度传感器放入被测环境,测量得到多个时刻下温度传感器的实测温度,所述多个时刻均远小于热平衡时间;
S3、根据所述系数矩阵、多项式,以及步骤S2中的多个时刻和实测温度计算得出待测温度。
与现有技术相比,本发明具有以下优点:在测量之前先计算出指定温度传感器的系数矩阵,在正式测温时仅测量远小于热平衡时间的多个时刻的输出温度,再结合系数矩阵和多项式即可计算得出待测温度,整个测量和计算过程所需的时间远小于热平衡时间,实现了利用普通温度传感器的快速测温。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,在这些附图中使用相同的参考标号来表示相同或相似的部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1为本发明的流程示意图;
图2为温度随时间变化的二维曲线。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,以下结合附图及具体实施例,对本申请作进一步地详细说明。为简单起见,以下描述中省略了本领域技术人员公知的某些技术特征。
如图1所示,本实施例提供一种利用普通温度传感器快速测温的方法,包括以下步骤:
S1、建立温度-时间响应模型,通过实验方法得到温度传感器在不同条件下输出温度与时间的关系曲线,采用多项式拟合关系曲线,得到由多项式的系数构成的系数矩阵;
S11、建立温度-时间响应模型,通过实验方法得到温度传感器输出温度与时间的关系曲线,采用多项式拟合关系曲线;
根据数学理论可知,任意的一条连续二维曲线,可以用多项式f(x)=a0+a1x+a2x2+...+anxn的形式进行拟合逼近。将温度传感器的输出温度值随时间变化的关系可以画出一条曲线,曲线的一般形式如图2所示,表征了温度传感器在某待测温度条件下的测温过程中,输出温度随时间变化的关系,其中T0表示温度传感器已达到热平衡状态,传感器的输出温度基本保持稳定。需要特别注意的是,同一型号传感器在不同待测温度条件下的曲线形状一般是不同的。根据二维曲线的多项式表示方法,可以将温度传感器在待测温度Tm下的输出温度表示为:
C(t)|Tm=am0+am1t+am2t2+...+amntn;
其中,为温度传感器的输出温度,t为测量时间,am0,am1,...amn为Tm温度下多项式的系数;
S12、在需要测量的温度范围内选择温度点[T1,T2,...,Tm],将外部环境温度依次固定在已选择的温度点,并在每个温度下记录温度传感器在多个时刻[tm1,tm2,...,tmn]下对应的输出温度[cm1,cm2,...,cmn],tmn应大于温度计热平衡所需的时间;
S13、根据多项式及步骤S12中记录的多个时刻和输出温度,计算得出在不同待测温度条件下多项式的系数,并由系数构成系数矩阵;
步骤S12、S13具体操作如下:
控制外部环境温度固定为T1,放入指定型号传感器并开始计时,得到时间[t11,t12,...,t1n]下对应的输出温度[c11,c12,...,c1n],采用最小二乘法拟合算法,可以计算得到在条件T1下,多项式的系数[a10,a11,...,a1n];
将外部环境温度固定为T2,放入指定型号传感器并开始计时,得到时间[t21,t22,...,t2n]下对应的输出温度[c21,c22,...,c2n],采用最小二乘法拟合算法,可以计算得到在T2条件下,多项式的系数[a20,a21,...,a2n];
以此类推,得到Tm条件下多项式的系数[am0,am1,...amn],将以上系数构成系数矩阵H如下:
在实际应用中,一般取n=5即可得到较为准确的曲线逼近效果,也可根据实际情况调整n的取值。系数矩阵H完整表示了该型号传感器在不同待测温度条件下,输出温度信号随时间的变化特征。
S2、将温度传感器放入被测环境,测量得到多个时刻Ttime=[t1,t2,...,tn]下温度传感器的实测温度C=[c1,c2,...,cn],多个时刻均远小于热平衡时间T0;此处的时刻数量与步骤S12中的时刻数量相同;
S3、根据系数矩阵、多项式,以及步骤S2中的多个时刻和实测温度计算得出待测温度;
S31、根据系数矩阵和多项式,计算得出步骤S2中每个时刻下的计算温度;具体如下:
将步骤S2中的t1和系数矩阵的每一行代入多项式,得到在t1时刻下的计算温度C1=[c11,c12,...,c1m];
将步骤S2中的t2和系数矩阵的每一行代入多项式,得到在t2时刻下的计算温度C2=[c21,c22,...,c2m];
同理,得到tn时刻下的计算温度Cn=[cn1,cn2,...,cnm];
S32、计算同一时刻下计算温度与步骤S2中实测温度的差值,并将步骤S2中所有时刻下计算温度与实测温度的差值构成偏离矩阵;具体如下:
在t1时刻下,实测温度c1与计算温度C1=[c11,c12,...,c1m]的差值为C1_delta=[c1-c11,c1-c12,...,c1-c1m];
在t2时刻下,实测温度c2与计算温度C2=[c21,c22,...,c2m]的差值为C2_delta=[c2-c21,c2-c22,...,c2-c2m];
同理,在tn时刻下实测温度cn与计算温度Cn=[cn1,cn2,...,cnm]的差值为Cm_delta=[cn-cn1,cn-cn2,...,cn-cnm];
将上述差值组成偏离矩阵Cdelta=[C1_delta T,C2_delta T,...,Cm_delta T];
S33、对偏离矩阵中的每一行所有元素进行平方和运算,选取平方和最小的那一行系数,平方和最小的一行下标记为b,对应系数矩阵中的系数为[ab0,ab1,...,abn],1<b<m;
S34、多项式中t取远大于热平衡时间T0的t0,多项式的系数取步骤S33得到的[ab0,ab1,...,abn],根据下式:
Cmeas=ab0+ab1t0+ab2t0 2+...+abnt0 n,
计算即可得出待测温度。
以上实施例仅表示本发明的几种实施方式,其描述较为具体和详细,但并不能理解为对本发明范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明保护范围。因此本发明的保护范围应该以权利要求为准。
Claims (3)
1.一种利用普通温度传感器快速测温的方法,其特征在于,包括以下步骤:
S1、建立温度-时间响应模型,通过实验方法得到温度传感器在不同条件下输出温度与时间的关系曲线,采用多项式拟合所述关系曲线,得到由多项式的系数构成的系数矩阵;
S2、将温度传感器放入被测环境,测量得到多个时刻下温度传感器的实测温度,所述多个时刻均远小于热平衡时间;
S3、根据所述系数矩阵、多项式,以及步骤S2中的多个时刻和实测温度计算得出待测温度。
2.根据权利要求1所述的利用普通温度传感器快速测温的方法,其特征在于,所述步骤S1具体包括:
S11、建立温度-时间响应模型,通过实验方法得到温度传感器输出温度与时间的关系曲线,用多项式拟合所述关系曲线:
C(t)|Tm=am0+am1t+am2t2+...+amntn;
其中,为温度传感器的输出温度,t为测量时间,am0,am1,...amn为Tm温度下多项式的系数;
S12、在需要测量的温度范围内选择温度点[T1,T2,...,Tm],将外部环境温度依次固定在已选择的温度点,并在每个温度下记录温度传感器在多个时刻[tm1,tm2,...,tmn]下对应的输出温度[cm1,cm2,...,cmn],tmn应大于温度计热平衡所需的时间;
S13、根据多项式及步骤S12中记录的多个时刻和输出温度,计算得出在不同待测温度条件下多项式的系数,并由所述系数构成系数矩阵;
其中,H为系数矩阵,a10,a11,...a1n为T1温度下多项式的系数,a20,a21,...a2n为T2温度下多项式的系数,am0,am1,...amn为Tm温度下多项式的系数。
3.根据权利要求2所述的利用普通温度传感器快速测温的方法,其特征在于,所述步骤S3具体包括:
S31、根据所述系数矩阵和多项式,计算得出步骤S2中每个时刻下的计算温度;
S32、计算同一时刻下计算温度与步骤S2中实测温度的差值,并将步骤S2中所有时刻下计算温度与实测温度的差值构成偏离矩阵;
S33、对偏离矩阵中的每一行所有元素进行平方和运算,选取平方和最小的系数[ab0,ab1,...,abn],1<b<m;
S34、根据所述多项式、系数[ab0,ab1,...,abn]和远大于热平衡时间的时间值计算得出待测温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710188091.XA CN106679848B (zh) | 2017-03-27 | 2017-03-27 | 一种利用普通温度传感器快速测温的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710188091.XA CN106679848B (zh) | 2017-03-27 | 2017-03-27 | 一种利用普通温度传感器快速测温的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106679848A true CN106679848A (zh) | 2017-05-17 |
CN106679848B CN106679848B (zh) | 2018-12-25 |
Family
ID=58826359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710188091.XA Expired - Fee Related CN106679848B (zh) | 2017-03-27 | 2017-03-27 | 一种利用普通温度传感器快速测温的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106679848B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109405989A (zh) * | 2018-10-16 | 2019-03-01 | 天津百腾信科技有限公司 | 一种高温升温环境下石墨坩埚温场的延展测量方法 |
CN111096649A (zh) * | 2018-10-26 | 2020-05-05 | 佛山市顺德区美的电热电器制造有限公司 | 运行控制方法、装置、烹饪器具和计算机可读存储介质 |
CN112179037A (zh) * | 2019-07-05 | 2021-01-05 | 青岛海尔电冰箱有限公司 | 冰箱间室温度的计算方法和装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4541734A (en) * | 1982-06-24 | 1985-09-17 | Terumo Kabushiki Kaisha | Electronic clinical thermometer, and method of measuring body temperature |
EP0338076A1 (en) * | 1986-12-24 | 1989-10-25 | Terumo Kabushiki Kaisha | Thermometric method for living organism and apparatus therefor |
CN1075205A (zh) * | 1992-12-11 | 1993-08-11 | 中国康复研究中心 | 快速皮肤温度计 |
CN102252770A (zh) * | 2011-04-29 | 2011-11-23 | 中冶赛迪工程技术股份有限公司 | 一种测温补偿方法及修正型高精度温度计 |
CN104535222A (zh) * | 2015-01-22 | 2015-04-22 | 哈尔滨工业大学 | 一种基于三价镨离子发光特性的高灵敏度温度测量方法 |
CN105823575A (zh) * | 2016-05-05 | 2016-08-03 | 北京航空航天大学 | 像素级热色液晶的标定方法和温度测量方法 |
-
2017
- 2017-03-27 CN CN201710188091.XA patent/CN106679848B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4541734A (en) * | 1982-06-24 | 1985-09-17 | Terumo Kabushiki Kaisha | Electronic clinical thermometer, and method of measuring body temperature |
EP0338076A1 (en) * | 1986-12-24 | 1989-10-25 | Terumo Kabushiki Kaisha | Thermometric method for living organism and apparatus therefor |
CN1075205A (zh) * | 1992-12-11 | 1993-08-11 | 中国康复研究中心 | 快速皮肤温度计 |
CN102252770A (zh) * | 2011-04-29 | 2011-11-23 | 中冶赛迪工程技术股份有限公司 | 一种测温补偿方法及修正型高精度温度计 |
CN104535222A (zh) * | 2015-01-22 | 2015-04-22 | 哈尔滨工业大学 | 一种基于三价镨离子发光特性的高灵敏度温度测量方法 |
CN105823575A (zh) * | 2016-05-05 | 2016-08-03 | 北京航空航天大学 | 像素级热色液晶的标定方法和温度测量方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109405989A (zh) * | 2018-10-16 | 2019-03-01 | 天津百腾信科技有限公司 | 一种高温升温环境下石墨坩埚温场的延展测量方法 |
CN111096649A (zh) * | 2018-10-26 | 2020-05-05 | 佛山市顺德区美的电热电器制造有限公司 | 运行控制方法、装置、烹饪器具和计算机可读存储介质 |
CN111096649B (zh) * | 2018-10-26 | 2021-11-19 | 佛山市顺德区美的电热电器制造有限公司 | 运行控制方法、装置、烹饪器具和计算机可读存储介质 |
CN112179037A (zh) * | 2019-07-05 | 2021-01-05 | 青岛海尔电冰箱有限公司 | 冰箱间室温度的计算方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106679848B (zh) | 2018-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Evaluation of different calibration equations for NTC thermistor applied to high-precision temperature measurement | |
Chen | Evaluation of resistance–temperature calibration equations for NTC thermistors | |
CN106679848A (zh) | 一种利用普通温度传感器快速测温的方法 | |
CN109000879A (zh) | 一种风洞天平温度漂移修正方法 | |
Martín et al. | Methodology for the calculation of response factors through experimental tests and validation with simulation | |
CN102620833B (zh) | 红外测温方法和红外测温系统 | |
CN103761370B (zh) | 一种板带热轧过程表面换热系数的预测方法 | |
CN102521439B (zh) | 结合有限元法和反传热法计算淬火介质换热系数的方法 | |
WO2020078213A1 (zh) | 红外光学材料均匀性测试的温度影响评估和控制方法 | |
CN113483900B (zh) | 基于黑体点在线校准的红外辐射铝合金板温度场测量方法 | |
CN104215660B (zh) | 一种可同时测固体材料导热系数及热扩散率的方法及系统 | |
CN110096755B (zh) | 固体蓄热炉内高温加热元件在线温度软测量方法及系统 | |
Cai et al. | Heat transfer prediction with unknown thermal boundary conditions using physics-informed neural networks | |
CN107630723A (zh) | 汽轮机转子热应力实时监测系统 | |
CN112965364A (zh) | 一种高超声速飞行器气动热地面模拟系统建模及iPID方法 | |
CN102879421A (zh) | 导热系数的测量装置及其测量方法 | |
CN104535257B (zh) | 一种硅压阻温度补偿评估方法 | |
CN110006554A (zh) | 一种温度计校准装置和方法 | |
CN112883624B (zh) | 一种拉制msf预制棒时对流换热模型的构建方法及装置 | |
CN104913862B (zh) | 基于圆箔热流计的对流热流测量方法 | |
CN103995017B (zh) | 一种测定周期性换热系数的实验方法 | |
CN114964562A (zh) | 一种薄膜热电阻和同轴热电偶的在线标定方法 | |
RU2566641C2 (ru) | Способ учета тепловой энергии, отдаваемой отопительным прибором | |
CN104897306A (zh) | 基于圆箔热流计的瞬态辐射热流测量方法 | |
François et al. | Building thermal bridge heat losses quantification by infrared thermography. Steady-state evaluation and uncertainty calculation |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181225 |