CN108414007A - 一种土壤温湿度传感器温度相关非线性补偿算法 - Google Patents

一种土壤温湿度传感器温度相关非线性补偿算法 Download PDF

Info

Publication number
CN108414007A
CN108414007A CN201810107576.6A CN201810107576A CN108414007A CN 108414007 A CN108414007 A CN 108414007A CN 201810107576 A CN201810107576 A CN 201810107576A CN 108414007 A CN108414007 A CN 108414007A
Authority
CN
China
Prior art keywords
temperature
nonlinear
equation
soil
relevant
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
Application number
CN201810107576.6A
Other languages
English (en)
Other versions
CN108414007B (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.)
Zhejiang University City College ZUCC
Original Assignee
Zhejiang University City College ZUCC
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 Zhejiang University City College ZUCC filed Critical Zhejiang University City College ZUCC
Priority to CN201810107576.6A priority Critical patent/CN108414007B/zh
Publication of CN108414007A publication Critical patent/CN108414007A/zh
Application granted granted Critical
Publication of CN108414007B publication Critical patent/CN108414007B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00—Measuring or testing not otherwise provided for
    • G01D21/02—Measuring two or more variables by means not covered by a single other subclass

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

本发明涉及一种土壤温湿度传感器温度相关非线性补偿算法，包括以下步骤：步骤一、非线性建模和系数估计是根据实验数据得到的，实验数据包括土壤湿度、土壤温度和真实的土壤湿度，得到表1；步骤二、根据实证Topps方程，体积含水率θ可以用水的介电常数εb作非线性描述：步骤三、用表1中每一行测量湿度θi,meas求方程，计算得到εb(i)；步骤四、用表1中每一行真实湿度θi,nl求方程，计算得到εb(i,T)。本发明的有益效果是：本发明提出了一个三阶的温度相关的非线性系统模型，并且提供了相应的非线性参数估计算法；本发明的土壤温湿度传感器温度相关非线性补偿算法精确度高，普适性强，参数少。

Description

一种土壤温湿度传感器温度相关非线性补偿算法
技术领域
本发明涉及土壤温湿度传感器，更具体说，它涉及一种土壤温湿度传感器温度相关非线性补偿算法。
背景技术
土壤湿度是农业产生中最重要的参数之一。采用传感器来测量土壤水分，可以提高水资源管理的效率。然而，由于湿度传感器的基本原理，传感器的响应值随湿度的变化呈非线性变化。并且，随着土壤真实湿度的增加和环境温度的变化，传感器测量的非线性也发生剧烈的变化。因此，土壤湿度传感器的温度相关的非线性特性难以标定。
在传统的湿度传感器中，校准通常是通过大量严格的测量来实现的。在每个测量中，有两个输入：真实湿度和温度，和一个输出：感测湿度。一个查找表(LUT)可以是通过改变真实湿度和温度的值来测量得到。在校准过程中，根据感测湿度和环境温度这两个坐标，可以在LUT中找到相应的真实湿度。基于LUT的标定方法简单有效，但其精度与复杂度之间存在严重的矛盾关系。如需达到较高的校准精度，则需要一个非常复杂和庞大的查找表，在内存有限的微控制器中难以实现。
为了克服这些缺点，有研究者建立了描述湿度传感器非线性的数学模型。例如，一个建立了真实体积含水率与传感器测得的探针电压之间非线性关系的三阶非线性数学表达式。然而，现有的非线性模型与环境温度无关，这对测量结果有很大影响。
发明内容
本发明的目的是克服现有技术的不足，提供一种土壤温湿度传感器温度相关非线性补偿算法。
这种土壤温湿度传感器温度相关非线性补偿算法，包括如下步骤：
步骤一、非线性建模和系数估计是根据实验数据得到的，实验数据包括土壤湿度、土壤温度和真实的土壤湿度，得到表1；
步骤二、根据实证Topps方程，体积含水率θ可以用水的介电常数εb作非线性描述如下：
水的介电常数随温度的升高而降低，因此温度相关的非线性可以由以下方程得出：
εb(i，T)＝εb(i)·C(i，T) (3)
C(i，T)＝1-a1(i)×(T-25)-a2(i)×(T-25)2-a3(i)×(T-25)3 (4)
其中θi，n1是修正后的温度相关的非线性模型。εb(i，T)是水的温度补偿介电常数，而C(i，T)是温度相关非线性因子。εb(i，T)和C(i，T)都是摄氏温度T的不同数值。am(m＝1，2，3)是估计的非线性系数；
步骤三、用表1中每一行测量湿度θi，meas求方程(1)，计算得到εb(i)；
步骤四、用表1中每一行真实湿度θi，n1求方程(2)，计算得到εb(i，T)；
步骤五、代入εb(i)和εb(i，T)到方程(3)，求得C(i，T)；
步骤六、应用最小二乘法算法求解方程(4)，得到估计的非线性系数am(m＝1，2，3)。
作为优选：所述步骤四中的最小二乘法算法可以用如下形式表示：
对于第i个测试，方程(4)可以用矩阵形式重写：
Ci＝ones(8，1)-T·ai (5)
其中矩阵如下：
(T)表示矩阵转置；
通过Ci和T，采用最小二乘法估计得到非线性系数ai，目标函数是
方程(6)的最小二乘解是
其中估计的非线性系数矩阵。
本发明的有益效果是:提出了一个3阶的温度相关的非线性系统模型，并且提供了相应的非线性参数估计算法。实验结果表明，本方法显著减少土壤湿度传感器测量误差。从而证明，土壤湿度传感器的非线性可以通过少量系数的数学模型得到很好的补偿，而没有必要采用大规模的查找表(LUT)。本发明的土壤温湿度传感器温度相关非线性补偿算法精确度高，普适性强，参数少。
附图说明
图1是系统架构比较图；
图2是温度非线性补偿方法的误差性能验证图。
具体实施方式
下面结合实施例对本发明做进一步描述。下述实施例的说明只是用于帮助理解本发明。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
本文提出了基于温度相关的非线性补偿模型，开发了一个集成的土壤湿度和温度传感器。非线性补偿模块包含一个温度相关的非线性模型，并给出了相应的系数估计算法。温度相关的非线性模型是通过对土壤湿度传感器测量结果的分析而得到的。采用最小二乘法来估计非线性系统模型的参数。采用本非线性补偿系统的集成式土壤湿度和温度传感器，其精度改进得到了实验结果的验证。
系统架构
常规温湿度传感器如图1(a)所示。该传感器集成了两个传感器元件，分别用于测量湿度和温度。在模数转换器(ADC)进行量化后，通过基于LUT的线性化模块对感测信号进行校准。校准内存应该足够大，以存储一个巨大的查找表，以实现高精度的湿度传感。但是，当传感器的非线性特性发生变化时，需要重写LUT。然而，大多数传感器不会被重新校准，从而随着使用时间的增长，它们的敏感度越低，可用生命周期变短。
如图1(b)所示，本文提出的传感器采用温度补偿非线性模型代替LUT。随温度变化的非线性建模需要三组数据，包括测量的湿度、温度和实际湿度。测量的湿度和温度由内部传感器提供。实际湿度数据通过RS485接口从外部传输。当非线性模型及其阶数确定后，最小二乘法算法用于估计模型系数。
非线性建模与系数估计
非线性建模和系数估计是根据实验数据得到的，实验数据包括土壤湿度、土壤温度和真实的土壤湿度，如表1所示。
根据实证Topps方程，体积含水率θ可以用水的介电常数εb作非线性描述如下：
表1湿度传感的实验数据
水的介电常数随温度的升高而降低，因此温度相关的非线性可以由以下方程得出：
εb(i，T)＝εb(i)·C(i，T) (3)
C(i，T)＝1-a1(i)×(T-25)-a2(i)×(T-25)2-a3(i)×(T-25)3 (4)
其中θi，n1是修正后的温度相关的非线性模型。εb(i，T)是水的温度补偿介电常数，而C(i，T)是温度相关非线性因子。εb(i，T)和C(i，T)都是摄氏温度T的不同数值。am(m＝1，2，3)是估计的非线性系数。
表1中数据用来进行系数估计，计算过程执行如下：
步骤1，用表1中每一行测量湿度θi,meas求方程(1)，计算得到εb(i)。
步骤2，用表1中每一行真实湿度θi,nl求方程(2)，计算得到εb(i,T)。
步骤3，代入εb(i)和εb(i,T)到方程(3)，求到C(i,T)。
步骤4，应用最小二乘法算法求解方程(4)，得到估计的非线性系数am(m＝1，2，3)。
具体而言，步骤4中的最小二乘法算法可以用如下形式表示：
对于第i个测试，方程(4)可以用矩阵形式重写：
Ci＝ones(8，1)-T·ai (5)
其中矩阵如下：
(T)表示矩阵转置。
通过Ci和T，采用最小二乘法估计得到非线性系数ai，目标函数是
方程(6)的最小二乘解是
其中估计的非线性系数矩阵。
利用表1的数据，估计非线性系数列于表2。
表2估计的非线性系数
实验结果与讨论
表2中估计的非线性系数被应用于验证与温度相关的非线性补偿的性能。代入am(m＝1，2，3)和εb(i)进入方程(2-4)，得到土壤真实湿度的估计对于每一个试验i，无非线性补偿和有非线性补偿的湿度误差率分别定义如下：
如图2所示，通过比较湿度误差比，验证了温度非线性补偿方法的性能。由此可见，这种非线性补偿可以显著降低土壤湿度传感器的误差。

Claims (2)

1.一种土壤温湿度传感器温度相关非线性补偿算法，其特征在于，包括如下步骤：
步骤一、非线性建模和系数估计是根据实验数据得到的，实验数据包括土壤湿度、土壤温度和真实的土壤湿度，得到表1；
步骤二、根据实证Topps方程，体积含水率θ可以用水的介电常数εb作非线性描述如下：
水的介电常数随温度的升高而降低，因此温度相关的非线性可以由以下方程得出：
εb(i，T)＝εb(i)·C(i，T) (3)
C(i，T)＝1-a1(i)×(T-25)-a2(i)×(T-25)2-a3(i)×(T-25)3 (4)
其中θi，n1是修正后的温度相关的非线性模型；εb(i，T)是水的温度补偿介电常数，而C(i，T)是温度相关非线性因子；εb(i，T)和C(i，T)都是摄氏温度T的不同数值；am(m＝1，2，3)是估计的非线性系数；
步骤三、用表1中每一行测量湿度θi，meas求方程(1)，计算得到εb(i)；
步骤四、用表1中每一行真实湿度θi，n1求方程(2)，计算得到εb(i，T)；
步骤五、代入εb(i)和εb(i，T)到方程(3)，求得C(i，T)；
步骤六、应用最小二乘法算法求解方程(4)，得到估计的非线性系数am(m＝1，2，3)。
2.根据权利要求1所述的土壤温湿度传感器温度相关非线性补偿算法，其特征在于，所述步骤四中的最小二乘法算法可以用如下形式表示：
对于第i个测试，方程(4)可以用矩阵形式重写：
Ci＝ones(8，1)-T·ai (5)
其中矩阵如下：
Ci＝[C(i，T＝5)，C(i，T＝10)，…，C(i，T＝40)]T
ones(8，1)＝[1，1，1，1，1，1，1，1]T
T＝[(T-25)|T＝5，(T-25)2|T＝5，(T-25)3|T＝5；
(T-25)|T＝10，(T-25)2|T＝10，(T-25)3|T＝10；
(T-25)|T＝40，(T-25)2|T＝40，(T-25)3|T＝40]
ai＝[a1(i)，a2(i)，a3(i)]T
T表示矩阵转置；
通过Ci和T，采用最小二乘法估计得到非线性系数ai，目标函数是
方程(6)的最小二乘解是
其中估计的非线性系数矩阵。
CN201810107576.6A 2018-02-02 2018-02-02 一种土壤温湿度传感器温度相关非线性补偿算法 Active CN108414007B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810107576.6A CN108414007B (zh) 2018-02-02 2018-02-02 一种土壤温湿度传感器温度相关非线性补偿算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810107576.6A CN108414007B (zh) 2018-02-02 2018-02-02 一种土壤温湿度传感器温度相关非线性补偿算法

Publications (2)

Publication Number Publication Date
CN108414007A true CN108414007A (zh) 2018-08-17
CN108414007B CN108414007B (zh) 2020-04-28

Family

ID=63127644

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810107576.6A Active CN108414007B (zh) 2018-02-02 2018-02-02 一种土壤温湿度传感器温度相关非线性补偿算法

Country Status (1)

Country Link
CN (1) CN108414007B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109405884A (zh) * 2018-12-03 2019-03-01 无锡华润矽科微电子有限公司 基于温湿度传感器的实现湿度校准功能的系统及方法
CN110568153A (zh) * 2019-08-21 2019-12-13 浙江大学城市学院 一种基于自适应阶数调节非线性模型的温湿度非线性补偿方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072922A (zh) * 2009-11-25 2011-05-25 东北林业大学 基于粒子群优化神经网络模型的木材含水率检测方法
CN103792340A (zh) * 2014-03-07 2014-05-14 国家海洋局第一海洋研究所 一种基于物联网的盐渍土监测预警系统及方法
CN105158106A (zh) * 2015-05-07 2015-12-16 中国科学院寒区旱区环境与工程研究所 一种寒区微型蒸渗仪测量系统及其对环境温度的补偿方法
CN107179392A (zh) * 2017-06-15 2017-09-19 江南大学 土壤探测器中基于最小二乘法的参数标定方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072922A (zh) * 2009-11-25 2011-05-25 东北林业大学 基于粒子群优化神经网络模型的木材含水率检测方法
CN103792340A (zh) * 2014-03-07 2014-05-14 国家海洋局第一海洋研究所 一种基于物联网的盐渍土监测预警系统及方法
CN105158106A (zh) * 2015-05-07 2015-12-16 中国科学院寒区旱区环境与工程研究所 一种寒区微型蒸渗仪测量系统及其对环境温度的补偿方法
CN107179392A (zh) * 2017-06-15 2017-09-19 江南大学 土壤探测器中基于最小二乘法的参数标定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈韦名等: "《一种湿度传感器温度补偿的非线性校正方法》", 《传感技术学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109405884A (zh) * 2018-12-03 2019-03-01 无锡华润矽科微电子有限公司 基于温湿度传感器的实现湿度校准功能的系统及方法
CN109405884B (zh) * 2018-12-03 2024-02-27 华润微集成电路(无锡)有限公司 基于温湿度传感器的实现湿度校准功能的系统及方法
CN110568153A (zh) * 2019-08-21 2019-12-13 浙江大学城市学院 一种基于自适应阶数调节非线性模型的温湿度非线性补偿方法

Also Published As

Publication number Publication date
CN108414007B (zh) 2020-04-28

Similar Documents

Publication Publication Date Title
CN103698056B (zh) 标定温度误差修正装置及修正方法
CN103162901A (zh) 一种压力传感器的多温度点非线性校准方法
CN102374902B (zh) 一种提高辐射温度计测温准确度的量子论修正方法
CN102096057B (zh) 一种电容测量电路的校准方法及装置
CN103575414A (zh) 一种热电偶温度测量误差补偿方法
CN106248296A (zh) 压力变送器的多变元、变温标定方法
CN206057424U (zh) 一种电流测量装置
CN109443601A (zh) 基于热电阻的高精度多点温度测量系统及温度校准方法
CN103983361A (zh) 在线式网络测温热像仪标定方法
CN107063509A (zh) 一种基于神经网络的热敏温度计校准方法
CN108414007A (zh) 一种土壤温湿度传感器温度相关非线性补偿算法
CN110514788A (zh) 一种气体传感器的标定方法、装置、设备和存储介质
CN104374514A (zh) 一种硅压阻式压力变送器的温度补偿方法
CN104390932A (zh) 基于红外差谱技术的木材含水率检测方法
CN105092097A (zh) 一种光纤光栅温度传感器的标定方法
CN102735152B (zh) 一种吸波涂层测厚仪的校准和测量方法
CN115265724A (zh) 一种液位计现场校准装置及获得校准结果的方法
CN109668813A (zh) 智能非饱和土土水特性曲线与渗透系数检测系统及方法
CN109253781A (zh) 节流式压差流量传感器的校准方法及校准系统
CN201218838Y (zh) 一种地基微波辐射计自然变温智能定标装置
CN116296047A (zh) 一种单晶硅压力变送器的温度补偿改进方法
CN115639335A (zh) 水质监测数据校准方法、系统及计算机可读存储介质
CN110568153A (zh) 一种基于自适应阶数调节非线性模型的温湿度非线性补偿方法
CN105389466A (zh) 一种校正尺度效应的中低分辨率遥感产品真值获取方法
Tong et al. Improving soil heat flux accuracy with the Philip correction technique

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