CN102645186B - 一种基于双线性插值算法的水膜厚度传感器实现方法 - Google Patents

一种基于双线性插值算法的水膜厚度传感器实现方法 Download PDF

Info

Publication number
CN102645186B
CN102645186B CN201210143023.9A CN201210143023A CN102645186B CN 102645186 B CN102645186 B CN 102645186B CN 201210143023 A CN201210143023 A CN 201210143023A CN 102645186 B CN102645186 B CN 102645186B
Authority
CN
China
Prior art keywords
film thickness
water film
value
temperature
current temperatures
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.)
Expired - Fee Related
Application number
CN201210143023.9A
Other languages
English (en)
Other versions
CN102645186A (zh
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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201210143023.9A priority Critical patent/CN102645186B/zh
Publication of CN102645186A publication Critical patent/CN102645186A/zh
Application granted granted Critical
Publication of CN102645186B publication Critical patent/CN102645186B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明公开一种基于双线性插值算法的水膜厚度检测传感器实现方法,该方法包括以下步骤:步骤a、测量不同温度Tk下不同水膜厚度Hi对应的幅度值Aik和相位值Pik,其中k={1,2,…,n};i={1,2,…,m},m代表选取水膜厚度值的个数,水膜厚度Hi值根据水膜厚度测量精度选择;步骤b、测量当前温度T0,对温度分段线性插值,得到当前温度下的幅度值Am0和相位值Pm0,其中T0∈[T1,Tn]。通过本发明的水膜厚度检测传感器的双线性插值算法减少了标定数据的采集量,提高了标定效率,加快了计算速度,改善了测量分辨率,降低了温漂现象,使得水膜厚度检测传感器的性能有了显著提高。

Description

一种基于双线性插值算法的水膜厚度传感器实现方法
技术领域
采用双线性插值算法的水膜厚度检测传感器实现方法属于微波测量技术的参数估计领域,具体用于对不同温度下水膜厚度标定数据的估计。
背景技术
近年来,高速公路建设的规模十分巨大,因此大量公路路面的监测和维护日益重要。公路路面在实际的使用过程中,积水会对整个路面结构产生不利的影响,并极大的增大了交通事故的发生几率。水膜厚度检测传感器能够准确测量路面积水厚度,从而实现公路路面的实时测量和预警,对于高速公路系统的信息化和智能化以及交通事故的预防预警有着十分重要的意义。水膜厚度检测传感器要求具有较高的分辨率,因此需要测量精确的标定数据以提高传感器的性能。然而,实际水膜厚度测试值受到环境温度的显著影响,导致温度漂移现象的发生,误差较大,同时由于实际测量条件限制,从而无法得到大量的不同温度下的水膜厚度标定数据以用于传感器的标定。本发明采用双线性插值的方法分别对温度和水膜厚度数据进行线性插值,插值方法简单,计算速度快,精度高,降低了温度效应对水膜厚度测量的影响,提高了测量分辨率,同时,数据最终通过最小二乘法查找进一步降低了测量误差。
发明内容
技术问题:本发明的目的是提出一种基于双线性插值算法的水膜厚度检测传感器实现方法,应用该方法可以显著减少标定数据的采集量,补偿温度效应误差,提高测量分辨率。
技术方案:为解决上述技术问题,本发明是提出一种基于双线性插值算法的水膜厚度检测传感器实现方法,该方法包括以下步骤:
步骤a、测量不同温度Tk下不同水膜厚度Hi对应的幅度值Aik和相位值Pik,
其中k={1,2,…,n},n代表选取温度值的个数,实际Tk值根据温度测量范围决定;i={1,2,…,m},m代表选取水膜厚度值的个数,水膜厚度Hi值根据水膜厚度测量精度选择;
步骤b、测量当前温度T0,对温度分段线性插值,得到当前温度下的幅度值Am0和相位值Pm0,其中T0∈[T1,Tn];
步骤c、对水膜厚度H分段线性插值,每段插入x个点,x为自然数,得到水膜厚度H下的幅度值A(T0,Hj)和相位值P(T0,Hj); T0为当前温度, Hj为不同水膜的厚度;
步骤d、实时测量得到当前温度T0和水膜厚度H下的幅度值A(T0,H)和相位值P(T0,H),运用最小二乘法查到最合适的二维点,即min={ [A(T0,Hj)- A(T0,H)]2+[P(T0,Hj)- A(T0,H)]2)}对应的j值,得到此时当前水膜厚度H;
水膜厚度H下的幅度值A和相位值P;当前温度T0 ,j代表差值后水膜点个数。
优选的,步骤b中,温度分段线性差值为:当前温度下的幅度值Am0和相位值Pm0分别为:
当T0∈[Ta,Ta+1],a∈[1,n-1]时,
Figure 489238DEST_PATH_IMAGE002
当T0< T1 或 Tn< T0时,设定相应的水膜测试边界条件;
其中,Ama代表温度Ta下水膜厚度Hm时幅度值,Am(a+1) 代表温度Ta+1下水膜厚度Hm时幅度值,Pma代表温度Ta下水膜厚度Hm时相位值,Pm(a+1) 代表温度Ta+1下水膜厚度Hm时相位值。
优选的,步骤c中水膜厚度H分段线性差值为:
H∈[Hb,Hb+1],b∈[1,m-1]
当前温度值T0下的幅度值Aj0和相位值Pj0
Figure 2012101430239100002DEST_PATH_IMAGE003
Figure 556420DEST_PATH_IMAGE004
x为对水膜厚度H分段线性插值,每段插入的点数;m代表选取测量水膜厚度值的个数,y代表每个差值区间计算值,j代表差值后水膜点个数;
当前温度下的幅度值Am0和相位值Pm0;
其中A(b+1)0, Ab0当前温度下的幅度值,P(b+1)0,Pb0均为当前温度下的相位值。
有益效果:本发明的双线性插值算法有效地降低了温度效应对水膜厚度测量值的影响,利用测量少量水膜厚度标定点插值得到全部的不同温度下的水膜厚度标定点,从而极大的减少了标定数据的采集量,提高了标定效率,此算法步骤简单,计算速度快,精度高,降低了温度效应对水膜厚度测量的影响,提高了测量分辨率,同时,数据最终通过最小二乘法查找进一步降低了测量误差,使得水膜厚度检测传感器的性能有了显著改善。
附图说明
图1是双线性插值算法流程图。
具体实施方式
下面结合附图对本发明做进一步说明。
本发明提供的基于双线性插值算法的水膜厚度检测传感器实现方法包括以下步骤:
步骤a、测量不同温度Tk下不同水膜厚度Hi对应的幅度值Aik和相位值Pik,
其中k={1,2,…,n},n代表选取温度值的个数,实际Tk值根据温度测量范围决定;i={1,2,…,m},m代表选取水膜厚度值的个数,水膜厚度Hi值根据水膜厚度测量精度选择;
步骤b、测量当前温度T0,对温度分段线性插值,得到当前温度下的幅度值Am0和相位值Pm0,其中T0∈[T1,Tn];
步骤c、对水膜厚度H分段线性插值,每段插入x个点,x为自然数,得到水膜厚度H下的幅度值A(T0,Hj)和相位值P(T0,Hj); T0为当前温度, Hj为不同水膜的厚度;
步骤d、实时测量得到当前温度T0和水膜厚度H下的幅度值A(T0,H)和相位值P(T0,H),运用最小二乘法查到最合适的二维点,即min={ [A(T0,Hj)- A(T0,H)]2+[P(T0,Hj)- A(T0,H)]2)}对应的j值,得到此时当前水膜厚度H;
水膜厚度H下的幅度值A和相位值P;当前温度T0 ,j代表差值后水膜点个数。
步骤b中,温度分段线性差值为:当前温度下的幅度值Am0和相位值Pm0分别为:
当T0∈[Ta,Ta+1],a∈[1,n-1]时,
Figure 2012101430239100002DEST_PATH_IMAGE005
Figure 620454DEST_PATH_IMAGE002
当T0< T1 或 Tn< T0时,设定相应的水膜测试边界条件;
其中,Ama代表温度Ta下水膜厚度Hm时幅度值,Am(a+1) 代表温度Ta+1下水膜厚度Hm时幅度值,Pma代表温度Ta下水膜厚度Hm时相位值,Pm(a+1) 代表温度Ta+1下水膜厚度Hm时相位值。
步骤c中水膜厚度H分段线性差值为:
H∈[Hb,Hb+1],b∈[1,m-1]
当前温度值T0下的幅度值Aj0和相位值Pj0
Figure 638274DEST_PATH_IMAGE004
x为对水膜厚度H分段线性插值,每段插入的点数;m代表选取测量水膜厚度值的个数,y代表每个差值区间计算值,j代表差值后水膜点个数;
当前温度下的幅度值Am0和相位值Pm0;
其中A(b+1)0, Ab0当前温度下的幅度值,P(b+1)0,Pb0均为当前温度下的相位值。
在具体实现的过程中,对温度T和水膜厚度H进行线性插值,插值过程中均采用分段插值法,分段区间和插值点的选择根据数据点的需要,测量条件的限制和传感器分辨率大小决定。水膜厚度测量区间为0mm-8mm,分辨率为0.1mm,温度选择测量四个温度点为整个温度区间的典型温度点,水膜厚度在大于1mm时每隔0.5mm选一个测量点,在小于1mm时因为测量实际条件限制无法精度测量,所以选择多插值法不测0.5mm数据。
不同温度和水膜厚度下标定数据
Figure 639724DEST_PATH_IMAGE006
表1
当前温度下不同水膜厚度数据
Figure 2012101430239100002DEST_PATH_IMAGE007
表2
当前温度下对水膜厚度线性插值结果
Figure 689589DEST_PATH_IMAGE008
表3
具体步骤如下所示:
a、测量温度T[k]下不同水膜厚度H[i]对应的幅度值A(T[k],H[i])和相位值P(T[k],H[i]),其中k={1,2,3,4},实际测量中选K=1表示T[1]=5°C,同样,T[2]=25°C,T[3]=45°C,T[4]=65°C,温度值选择根据温度测量范围决定;i={0,1,2,3,4,5,6,7,…,15},i代表水膜厚度值i=0代表水膜厚度0mm,i=1表示水膜厚度1mm,后面为每隔0.5mm测一组数据所以i=15代表8mm水膜,实际i值根据水膜测量精度选择;
b、对温度分段线性插值,得到不同温度T下的幅度值A(T,H[i])和相位值P(T,H[i]),其中T∈[T[1],T[4]],i值不变;
当T[1]< T< T[2]< T[3]< T[4]时,
Figure 976213DEST_PATH_IMAGE009
Figure 601492DEST_PATH_IMAGE010
当T[1]< T[2] <T < T[3]< T[4]时,
Figure 948160DEST_PATH_IMAGE011
Figure 485320DEST_PATH_IMAGE012
当T[1] < T[2]< T[3]< T< T[4]时,
Figure 2012101430239100002DEST_PATH_IMAGE013
Figure 77101DEST_PATH_IMAGE014
当T< T[1] 或 T[4]< T时,设定相应的水膜测试边界条件。
c、测量当前温度T0,得到当前温度T0下对应的一组幅度值A(T0,H[i])和相位值P(T0,H[i]);
d、当前温度下,对水膜厚度H分段线性插值得到一定精度下的幅度值A(T0,H[j])和相位值P(T0,H[j]),其中j={0,1,2,3,…,80}共81组数据,j=1代表水膜厚度0.1mm,依次类推得到j=80代表水膜厚度8.0mm,这样就得将水膜厚度分辨率提高到0.1mm,实际插值j根据水膜分辨率设定;
当0mm≤H<1mm时,
Figure 2012101430239100002DEST_PATH_IMAGE015
当1mm≤H<8mm时
Figure 445634DEST_PATH_IMAGE016
    当H=8mm时,
;  
Figure 261406DEST_PATH_IMAGE018
e、硬件测量得到当前温度T0和水膜厚度H0下的幅度值A(T0,H0)和相位值P(T0,H0),运用查找最小二乘法在幅度值A(T0,H[j])和相位值P(T0,H[j])查找到最合适的二维点,即min={ [A(T0,H[j])- A(T0,H0)]2+[P(T0,H[j])- A(T0,H0)]2)}对应的j值,得到此时当前水膜厚度H0 =j×0.1mm。
最小二乘法:
a、设初始最小值min={ [A(T0,H[j])- A(T0,H0)]2+[P(T0,H[j])- A(T0,H0)]2)},j=0;
b、顺序在表3中查找,当min>{ [A(T0,H[j])- A(T0,H0)]2+[P(T0,H[j])- A(T0,H0)]2)},j∈{1,2,3,…,80}时,替换最小值min={ [A(T0,H[j])- A(T0,H0)]2+[P(T0,H[j])- A(T0,H0)]2)},此时最小值代表的j点即为和当前水膜厚度最接近的点,即为标定点;
以上所述仅为本发明的较佳实施方式,本发明的保护范围并不以上述实施方式为限,但凡本领域普通技术人员根据本发明所揭示内容所作的等效修饰或变化,皆应纳入权利要求书中记载的保护范围内。

Claims (3)

1. 一种基于双线性插值算法的水膜厚度检测传感器实现方法,其特征在于:该方法包括以下步骤:
步骤a、测量不同温度Tk下不同水膜厚度Hi对应的幅度值Aik和相位值Pik,
其中k={1,2,…,n},n代表选取温度值的个数,实际Tk值根据温度测量范围决定;i={1,2,…,m},m代表选取水膜厚度值的个数,水膜厚度Hi值根据水膜厚度测量精度选择;
步骤b、测量当前温度T0,对温度分段线性插值,得到当前温度下的幅度值Am0和相位值Pm0,其中T0∈[T1,Tn];
步骤c、对水膜厚度H分段线性插值,每段插入x个点,x为自然数,得到水膜厚度H下的幅度值A(T0,Hj)和相位值P(T0,Hj); T0为当前温度, Hj为不同水膜的厚度;
步骤d、实时测量得到当前温度T0和水膜厚度H下的幅度值A(T0,H)和相位值P(T0,H),运用最小二乘法查到最合适的二维点,即min={ [A(T0,Hj)- A(T0,H)]2+[P(T0,Hj)- A(T0,H)]2)}对应的j值,得到此时当前水膜厚度H;
水膜厚度H下的幅度值A和相位值P;当前温度T0 ,j代表差值后水膜点个数。
2. 根据权利要求1所述的基于双线性插值算法的水膜厚度检测传感器实现方法,步骤b中,温度分段线性差值为:当前温度下的幅度值Am0和相位值Pm0分别为:
当T0∈[Ta,Ta+1],a∈[1,n-1]时,
Figure 233467DEST_PATH_IMAGE002
当T0< T1 或 Tn< T0时,设定相应的水膜测试边界条件;
其中,Ama代表温度Ta下水膜厚度Hm时幅度值,Am(a+1) 代表温度Ta+1下水膜厚度Hm时幅度值,Pma代表温度Ta下水膜厚度Hm时相位值,Pm(a+1) 代表温度Ta+1下水膜厚度Hm时相位值。
3. 根据权利要求1所述的基于双线性插值算法的水膜厚度检测传感器实现方法,步骤c中水膜厚度H分段线性差值为:
H∈[Hb,Hb+1],b∈[1,m-1]
当前温度值T0下的幅度值Aj0和相位值Pj0
Figure 46231DEST_PATH_IMAGE006
Figure 922920DEST_PATH_IMAGE008
x为对水膜厚度H分段线性插值,每段插入的点数;m代表选取测量水膜厚度值的个数,y代表每个差值区间计算值,j代表差值后水膜点个数;
当前温度下的幅度值Am0和相位值Pm0;
其中A(b+1)0, Ab0当前温度下的幅度值,P(b+1)0,Pb0均为当前温度下的相位值。
CN201210143023.9A 2012-05-10 2012-05-10 一种基于双线性插值算法的水膜厚度传感器实现方法 Expired - Fee Related CN102645186B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210143023.9A CN102645186B (zh) 2012-05-10 2012-05-10 一种基于双线性插值算法的水膜厚度传感器实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210143023.9A CN102645186B (zh) 2012-05-10 2012-05-10 一种基于双线性插值算法的水膜厚度传感器实现方法

Publications (2)

Publication Number Publication Date
CN102645186A CN102645186A (zh) 2012-08-22
CN102645186B true CN102645186B (zh) 2014-05-07

Family

ID=46658161

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210143023.9A Expired - Fee Related CN102645186B (zh) 2012-05-10 2012-05-10 一种基于双线性插值算法的水膜厚度传感器实现方法

Country Status (1)

Country Link
CN (1) CN102645186B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103674285B (zh) * 2013-11-30 2016-04-13 江苏物联网研究发展中心 具有自校正功能的红外热电堆温度传感器的温度计算模块及温度计算方法
CN104714967B (zh) * 2013-12-14 2017-10-24 中国航空工业集团公司第六三一研究所 一种确定汽车风扇控制扭矩的方法
CN110244565B (zh) * 2019-06-21 2020-08-04 清华大学 一种scr系统分区控制方法和装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0432360A3 (en) * 1989-12-09 1991-10-23 Tzn Forschungs- Und Entwicklungszentrum Unterluess Gmbh Procedure for determining the mean thickness of a water film on the street surface
CN202041495U (zh) * 2010-12-17 2011-11-16 交通运输部公路科学研究院 一种被动式路面天气传感器
CN102283709A (zh) * 2011-05-26 2011-12-21 福建师范大学 一种水膜层厚度的精确控制调节装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0432360A3 (en) * 1989-12-09 1991-10-23 Tzn Forschungs- Und Entwicklungszentrum Unterluess Gmbh Procedure for determining the mean thickness of a water film on the street surface
US5105157A (en) * 1989-12-09 1992-04-14 Tzn Forschungs- Und Entwicklungszentrum Unterluss Gmbh Apparatus for determining the average water film thickness of road surfaces
CN202041495U (zh) * 2010-12-17 2011-11-16 交通运输部公路科学研究院 一种被动式路面天气传感器
CN102283709A (zh) * 2011-05-26 2011-12-21 福建师范大学 一种水膜层厚度的精确控制调节装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
季天剑.道路表面水膜厚度预测模型.《交通运输工程学报》.2004,第4卷(第3期),第1-3页.
道路表面水膜厚度预测模型;季天剑;《交通运输工程学报》;20040930;第4卷(第3期);第1-3页 *

Also Published As

Publication number Publication date
CN102645186A (zh) 2012-08-22

Similar Documents

Publication Publication Date Title
CN108254032B (zh) 河流超声波时差法流量计算方法
CN102879639A (zh) 一种电力系统中频率的实时测量方法
CN109543356A (zh) 考虑空间非平稳性的海洋内部温盐结构遥感反演方法
CN108333617B (zh) 一种海水中放射性物质检测的快速寻峰方法
CN102645186B (zh) 一种基于双线性插值算法的水膜厚度传感器实现方法
CN103235883A (zh) 一种基于非接触式河流表面测速的瞬时流量回归计算方法
CN107179568A (zh) 一种降水预报精度评定方法
CN103954842A (zh) 一种确定大型接地系统接地阻抗值的方法
CN109271466A (zh) 一种基于层次聚类与k均值算法的气象数据分析方法
CN102426352B (zh) 风廓线雷达风计算方法
CN105021663A (zh) 一种盐度测量方法
CN104392113B (zh) 一种近海海面冷空气大风风速的估算方法
CN103559411A (zh) 一种基于数字高程模型的洪水演进数值模拟计算修正方法
CN105068075B (zh) 一种近地面大风的计算方法
CN106228021A (zh) 农田风蚀量预测模型及风蚀量预测方法
CN103033274B (zh) 一种日照温度概率密度的测定方法
CN104422391A (zh) 传感式物体体积识别技术
CN109946765B (zh) 风电场的流场的预测方法和系统
CN107202661B (zh) 一种基于伪最大公约数的拉索振动频阶识别方法
CN105389466A (zh) 一种校正尺度效应的中低分辨率遥感产品真值获取方法
CN101556255A (zh) 一种围护结构热阻现场检测数据分析方法
CN108919313A (zh) 利用最优数值导数的gnss多普勒观测值生成方法
CN104032630A (zh) 一种基于角度传递的轨道沉降连续测量方法
CN107014354B (zh) 一种倾角传感器数据采集方法
CN104713631B (zh) 一种油井内平均声速的检测方法

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
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: 20140507

Termination date: 20180510