CN104132967B - 一种基于双阈值随机共振的低浓度气体检测方法 - Google Patents
一种基于双阈值随机共振的低浓度气体检测方法 Download PDFInfo
- Publication number
- CN104132967B CN104132967B CN201410277331.XA CN201410277331A CN104132967B CN 104132967 B CN104132967 B CN 104132967B CN 201410277331 A CN201410277331 A CN 201410277331A CN 104132967 B CN104132967 B CN 104132967B
- Authority
- CN
- China
- Prior art keywords
- dual threshold
- gas
- threshold
- correlation coefficient
- concentration
- 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.)
- Active
Links
- 230000009977 dual effect Effects 0.000 title claims abstract description 92
- 238000001514 detection method Methods 0.000 title claims abstract description 33
- 238000012905 input function Methods 0.000 claims abstract description 16
- 238000012417 linear regression Methods 0.000 claims abstract description 10
- 230000000737 periodic effect Effects 0.000 claims abstract description 8
- 230000004044 response Effects 0.000 claims description 20
- 238000000034 method Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 5
- 229910052756 noble gas Inorganic materials 0.000 claims 1
- 238000006467 substitution reaction Methods 0.000 abstract description 2
- 238000004064 recycling Methods 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 55
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 18
- 229910021529 ammonia Inorganic materials 0.000 description 9
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 6
- 230000009897 systematic effect Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 229910052757 nitrogen Inorganic materials 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004566 IR spectroscopy Methods 0.000 description 1
- 229920001609 Poly(3,4-ethylenedioxythiophene) Polymers 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 150000001336 alkenes Chemical class 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 210000004087 cornea Anatomy 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000023077 detection of light stimulus Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000009774 resonance method Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000002109 single walled nanotube Substances 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于双阈值随机共振的低浓度气体检测方法,包括:测量气敏传感器在不同浓度的参照气体下的电阻值;以测得的电阻值作为信号输入,并对该信号进行归一化处理,再利用傅里叶级数对进行曲线拟合得到周期函数;对周期函数进行离散化处理得到输入函数;将输入函数和外加的高斯白噪声输入双阈值随机共振系统,利用双阈值检测器进行赋值,并以赋值后的信号作为输出信号,并计算双阈值随机共振系统的互相关系数;搜索最优的双阈值使得互相关系数最大,以最大的互相关系数和对应气体浓度建立线性回归方程;选取待测的低浓度气体,得到待测低浓度气体对应的最大互相关系数,代入所述的线性回归方程,即求出待测低浓度气体的浓度。
Description
技术领域
本发明涉及一种在强噪声下的微弱信号检测方法,尤其涉及一种基于双阈值随机共振的低浓度气体检测方法。
背景技术
随着对低浓度气体检测研究的不断深入,呈现出越来越多的检测方法。超声波技术、光干涉原理、红外吸收光谱原理、环形光路的Sagnac效应及腔衰荡测量等技术都可以应用到低浓度的测量。但是由于灵敏度不够,这些测量技术只能在实验室的条件下进行,而且费时、费力、成本高,难以普及。
随机共振理论为强噪声背景下微弱信号检测提供了新思路。随机共振现象在一些领域,比如信号处理、神经信息学、机械误差检测、环境保护等一直备受关注,在信号处理方面的研究,特别将随机共振应用于弱信号的放大识别、传输还原等有其独特的优势。这种现象通常在一定程度上可以描述为提高输出信号的信噪比,是弱信号特征提取的重要组成部分。
如公开号为CN103698371A的发明专利中提供了一种基于自适应随机共振的低浓度气体检测方法,对采集到的传感器响应的电阻信号进行预处理,使之成为小参数周期信号,适合于非线性双稳系统的输入,设置系统参数的初始值和寻优范围,利用加权信噪比作为效果评价指标,自适应的寻找最佳匹配的系统参数,通过比较加权信噪比的最大值来区分不同浓度的确定性混合气体。
随机共振是一种非线性现象,通过加入适当噪声可以增强系统的响应,对系统性能有显著的提高。近些年,单阈值随机共振系统已经得到一定的应用。2000年,Stocks最早提出多阈值随机共振理论,认为多阈值检测器相对于单阈值检测器来说有效的扩展了电平信号的动态范围。多阈值检测器可以由单阈值检测器发展而来,在一个简单的阈值系统中选用多个适当的数值作为阈值,并用互相关系数和信息量作为表征。但是在实际应用中,多阈值随机共振还没有得到很大的发展,特别是在低浓度气体检测方面的研究微乎其微。
在低浓度气体检测方面,由于气体浓度很低,检测信号会淹没在噪声中。采用传感器阵列和光谱法的气体检测技术,存在着灵敏度不够高、噪声大等诸多问题。针对这种不足,本发明提出了一种双阈值随机共振的低浓度气体检测方法。
发明内容
本发明的目的在于克服弱信号在检测过程中淹没在噪声背景中的缺陷,以及阈值随机共振系统评价指标的局限性,因而提出了一种利用双阈值随机共振系统提取最大互相关系数作为评价指标的低浓度气体检测方法。
一种基于双阈值随机共振的低浓度气体检测方法,包括以下步骤:
1)将气敏传感器置于不同浓度的参照气体(浓度已知)所处的环境中,测得气敏传感器在各浓度下的电阻值;
2)以测得的电阻值作为信号输入,并对该信号进行归一化处理,再利用傅里叶级数对归一化处理后的信号进行曲线拟合,得到周期函数f(t);
3)对周期函数f(t)进行离散化处理得到输入函数S(t);
4)将输入函数S(t)和外加的高斯白噪声n(t)输入双阈值随机共振系统,利用双阈值检测器对输入函数S(t)和外加的高斯白噪声n(t)进行赋值,并以赋值后的信号作为输出信号y(t):
其中,θ1和θ2为双阈值随机共振系统的双阈值;
5)利用所述的输入函数S(t)和输出信号y(t),计算双阈值随机共振系统的互相关系数;
6)搜索最优的双阈值θ1和θ2,使得步骤5)中的互相关系数最大,并以最大的互相关系数和对应的参照气体的浓度建立线性回归方程;
7)选取待测的低浓度气体,执行步骤1)~步骤6)中的操作,得到待测低浓度气体对应的最大互相关系数,代入步骤6)中的线性回归方程,即可求出待测低浓度气体的浓度。
在本发明中,参照气体和待测的低浓度气体为同一种气体,参照气体用于数据的采集和线性回归方程的建立,其浓度是已知的,待测的低浓度气体用于验证本发明中技术方案的准确性。
在步骤2)中,利用下式方程对所测的电阻值进行归一化处理,使其值域为[0,1];
f1(t)=(x(t)-xmin)/(xmax-xmin)
其中:x(t)为测量得到的电阻值,xmin为最小电阻值,xmax为最大电阻值。
在这里归一化使数据具有可参照性,而且使数据之间的相对大小不至于差别很大,消除不同数量级带来的影响。
在步骤2)中,进行曲线拟合得到的周期T的函数f(t)为:
其中:a0、an、bn为函数系数,ancos nwt+bnsin nwt为n阶谐波,n=1,2,3,t为气体检测的时间。
曲线拟合不仅可以增大输入的信号量,还可以对原有信号的毛刺进行平滑处理,减少原信号的误差。
在步骤4)中,高斯白噪声n(t)的均值为0,方差为σ,概率密度函数为其中x为连续随机变量。
在步骤5)中,互相关系数的表达式如下:
其中,CCC为互相关系数,S(t)为输入函数,y(t)为输出信号,为S(t)的平均值,为y(t)的平均值。
在步骤6)中,搜索最优双阈值θ1和θ2的步骤如下:
(1)对双阈值随机共振系统的双阈值θ1和θ2进行初始化,并确定阈值的搜索范围和步长;
(2)将输入函数S(t)和外加的高斯白噪声n(t)输入初始化后的双阈值随机共振系统中,计算每一组双阈值所对应的拟合曲线的标准差SD;
(3)在搜索范围内,根据步长进行遍历搜索,得到标准差SD的最小值,以及相对应的双阈值θ1、θ2。
本发明中,线性拟合曲线与最大互相关系数之间的标准差SD为:
其中:N代表所有不同浓度气体的数量,CCCi代表的各浓度对应的最大互相关系数,yi代表的是各浓度经过曲线拟合后的值。
在本发明中,阈值θ1的搜索范围为[-1.0,0.26],系统阈值θ2的搜索范围为[0.89,1.5],搜索步长均为0.02。
优选的,在步骤6)中,建立的线性回归方程为y=0.0028x+0.2051,y为最大互相关系数,x为气体的浓度。
为使测得的电阻值更加准确,在步骤1)中,首先向测试气体的气室中通入惰性气体,使气敏传感器的响应值达到稳定状态,再向气室中通入特定浓度的参照气体,测量3次取平均值,作为气敏传感器在此浓度下的电阻值。
本发明用最大互相关系数表征的双阈值随机共振算法搜索最佳的阈值参数,其特点在于应用多阈值和随机共振相结合,并利用最大互相关系数作为系统表征,可以自动搜索最佳的双阈值参数。设计的双阈值参数调节随机共振方法可以快速、准确的检测未知气体浓度。在实际工程应用中克服了信噪比和信息量作为评价指标的不足和随机共振的系统参数难以确定的限制,有效的实现了低浓度气体的检测。通过最大互相关系数可以标定未知气体浓度,区分不同浓度的确定性混合气体,灵敏而高效。
附图说明
图1为本发明提供的双阈值随机共振实现流程示意图。
图2为传感器气敏测试系统示意图。
图3为6种浓度氨气进过归一化处理,得到传感器响应与时间的关系。
图4为种氨气浓度的传感器响应值和拟合曲线得到的线性关系。
图5为双阈值随机共振的系统模型。
图6和图7分别为不同双阈值条件下气体浓度为3ppb和60ppb时互相关系数与噪声强度的关系图,用于划分随机共振系统双阈值的具体范围,大量缩小搜素阈值的时间。
图8是在阈值θ1=-0.22、θ2=1.08、噪声强度σ=[0,4]条件下得到的不同浓度下的互相关系数与噪声强度关系图。
图9为最大互相关系数与其相对应气体的浓度的线性关系。
具体实施方式
下面用一个实例并结合附图进一步说明本发明的技术方案。
如图1所示,一种基于双阈值随机共振的低浓度气体检测方法,包括以下步骤:
第一步:利用数字万用表和LabView软件记录传感器阻值。
具体为:传感器采用DEP制备PEDOT/PSS-SWCNTs气敏传感器,该传感器的检测范围是3ppb-1ppm。如图2所示,制备的气敏传感器放入气室4中,并且在气室4中充满聚四氟乙烯。在实验过程中,以两个气瓶内的高纯氮气(>99.99%)作为载体和氨气作为目标气体。质量控制器1(VICI,Adjustable span,Model202)的作用是控制所需的浓度,流量计2确定混合气体的流速为0.2L/min,用数字万用表3和LabView软件记录传感器阻值。
所有的气敏测试实验均在室温20℃、空气相对湿度为70%和一个标准大气压的环境下进行。测试前向气室通入一段时间的氮气,使传感器的响应值达到稳定状态,得到传感器的电阻值R0。向气室中通入氨气,通过质量控制器使氨气浓度在3ppb,实验重复3次取平均值,记录此时传感器的电阻值R。然后向气室中通入氮气进行洗气,再向气室中通入6ppb的氨气,重复上面实验。以此类推可以得到6ppb、12ppb、18ppb、40ppb、60ppb氨气的传感器响应值。传感器对目标气体的响应值为:
Res=(R-R0)/R0×100%
六种氨气浓度的传感器响应如图3、图4所示。
图3中横坐标代表时间轴,纵坐标代表的是传感器的电阻响应值,是数据归一化后得来的,原始数据没有经过坐标转换。图4中横坐标代表浓度,纵坐标代表传感器的响应值,根据公式Res计算而来。从图4看到每一浓度的响应值大致趋势呈上升,但是线性关系不是很理想。
第二步:对所有信号进行预处理,利用傅里叶级数进行曲线拟合,得到一个周期函数,离散化处理后得到输入函数。
具体为:在输入信号和高斯白噪声都已知时,可以通过双阈值模型来产生随机共振现象。图5是双阈值随机共振的系统模型,其中输入信号S(t)是一个离散函数,通过对原始数据进行预处理得到。利用下式方程对所测的电阻值进行归一化处理,使其值域为[0,1]。
f1(t)=(x(t)-xmin)/(xmax-xmin)
其中x(t)为测量得到的电阻值,xmin为最小电阻值,xmax为最大电阻值,得到如图3所示的传感器响应值。利用傅里叶级数对f1(t)进行曲线拟合,得到的周期T的函数f(t)为:
其中:a0、an、bn为函数系数,ancosnwt+bnsinnwt为n阶谐波,n=1,2,3,t为气体检测的时间。
然后把上面得到周期函数f(t)进行离散化处理,得到输入信号S(t),方程(输入函数)为:S(t)=(s1,s2,s3,......,sn)。系统添加的噪声信号n(t)是均值为0,方差为σ,x为连续随机变量,概率密度函数为的高斯白噪声。最后把S(t)和n(t)作为系统的输入,利用双阈值交叉模型建立数学模型。
第三步:把输入函数和高斯白噪声代入双阈值交叉模型中,通过双阈值检测器来判断赋予输出函数的值,然后计算系统的互相关系数。
在通过双阈值检测器的过程中可以调节双阈值θ1、θ2来优化系统,使系统性能达到最佳,进而获得最佳的最大互相关系数作为对应浓度的特征参数,找到最佳的线性关系。
具体为:图5中双阈值检测器的函数模型为
式中y(t)为输出信号,式中θ1、θ2是两个阈值,其中θ1<θ2,n(t)是高斯白噪声,和S(t)作为输入信号,通过“双阈值检测器”被赋予一系列的-1、0、1,使其作为输出信号。根据具体情况优化系统的阈值,并判断输入信号是否产生随机共振。
本发明利用互相关系数描述输入信号和输出信号之间的相似度。如果两个信号完全相同,则相关系数为1,如果完全不相关,则相关系数为0。在双阈值随机共振系统中,通过上面叙述可以得到输入信号S(t)和输出信号y(t)。通过var(S(t))和var(y(t))的方程计算出S(t)和y(t)的方差,cov(S(t),y(t))的方程计算S(t)和y(t)的协方差。CCC的方程是互相关系数的表达式,把式var(S(t))、var(y(t))、cov(S(t),y(t))代入CCC的方程中,就可以得到系统的互相关系数。
据上所述,可以计算任意双阈值随机共振系统中互相关系数和噪声强度的关系。选择最佳的双阈值,使系统性能达到最佳状态,才能更全面的研究双阈值随机共振系统检测弱信号的性能。而要找到最佳的双阈值,必须先缩小双阈值的值域。通过仿真实验发现阈值θ1的范围是随着浓度的增加而增大,阈值θ2的范围是随着浓度的增加而减小,所以选择3ppb和60ppb进行分析可以找到双阈值的范围。噪声强度范围设定为[0,4],取输入信号的周期范围为[0,620s]。图6、7给出了3ppb和60ppb时双阈值下的随机共振关系图。
图6给出了氨气浓度为3ppb时,互相关系数随噪声强度的变化规律,并同时绘出了双阈值各不相同时的最大互相关系数的情况。当θ2固定不变,θ1<0.26时可以发生随机共振现象,且随着θ1的减小,最大互相关系数也减小。当θ1固定不变,θ2>0.45的时候可以发生随机共振,它对系统的影响占的比重较大,随着θ2的增加,最大互相关系数相应减小。图7给出了氨气浓度为60ppb时,最大互相关系数随阈值的变化情况。但是只有阈值θ1<0.35和阈值θ2>0.89的时候才能发生随机共振,它的变化趋势和图a是相似的。在对这两幅图进行综合分析,当阈值θ2>1.5时,图6和图7的最大互相关系数已经接近,不利于后续的线性回归。当阈值θ1<-1.0的时候,图6和图7的最大互相关系数都是越来越小,他们之间的差值越来越难区分。综上分析可知:系统阈值θ1的范围为[-1.0,0.26],系统阈值θ2的范围是[0.89,1.5]。
线性拟合曲线与各浓度对应的最大互相关系数之间的标准差定义为:
其中N代表所有浓度的个数,在本示例中N=6,CCCi代表的是各浓度对应的最大互相关系数,yi代表的是各浓度经过曲线拟合后的值。在已知双阈值值域的情况下,可以利用拟合曲线线性度的好坏来评价随机共振系统的优劣。进而调节双阈值的大小,以得到优化结果。当SD达到最小值的时候,双阈值为随机共振系统的最优参数,随机共振效应也是最佳的。
本发明提出的方法可以实现双阈值θ1、θ2的自动调整,并找到最优化的参数。利用标准差SD作为传感器响应特征值的线性度的评价参数,具有快速、精确的优化标准,实现弱信号特征的有效提取。
获取最优化双阈值的步骤如下:
3.1、双阈值系统参数初始化,使噪声强度σ=[0,4],初始化阈值θ1=-1.0、θ2=0.89,确定系统的搜索步长为0.02。
3.2、将预处理后的信号和外加的高斯白噪声作为输入信号输入到初始化以后的双阈值随机共振系统中,计算每一组双阈值所对应的拟合曲线的标准差SD。
3.3、在规定双阈值值域的范围内,进行遍历搜索,寻找标准差SD的最小值,当双阈值都超出搜索范围时,保存SD的最小值和对应的双阈值θ1、θ2。保存最好的线性关系的直线,并记录对应双阈值θ1、θ2的值。当双阈值没有超出搜索范围时,改变双阈值θ1、θ2的值,继续循环进行步骤3.2。
3.4、输出标准差SD的最小值,以及所对应的双阈值θ1、θ2。
在以上优选过程中,得到最优化的双阈值θ1=-0.22、θ2=1.08,标准差SD=0.0126,相关系数R=0.9687。通过双阈值随机共振系统得到的不同浓度下的互相关系数与噪声强度关系图如图8所示。对应的低浓度气体的浓度与最大互相关系数的线性关系如图9所示,方程为:
y=0.0028x+0.2051
式中,y1为最大互相关系数,x为气体的浓度,单位为ppb。在已知双阈值θ1、θ2时,可以求出未知浓度的最大互相关系数,代入上述方程,求出对应的浓度。
第四步:然后利用第三步得到的最大互相关系数,代入步骤3.4得到的回归方程可以求出相关的浓度。
具体为:通过第三步可以知道双阈值随机共振系统的系统参数。双阈值θ1=-0.22、θ2=1.08,噪声强度σ=[0,4],并且知道了对应的浓度与最大互相关系数之间的线性关系。因此,在确定系统参数的情况下,通过步骤三可以求出未知浓度的最大化相关系数,可以利用已知的线性关系求出对应的气体浓度。
在本实例中,经过双阈值随机共振系统后传感器响应的特征值如表1所示,
表1经过双阈值随机共振系统后传感器响应的特征值
在没有进行随机共振的情况下,如图4所示传感器响应值呈现不明显的线性关系。通过实验发现用相关系数进行表征的随机共振系统可以增强检测精度,提高检测限。因此把获取的原始信号进行预处理得到输入信号,并代入双阈值随机共振系统中计算互相关系数,在获取最优化双阈值的情况下可以得到如图8。图8是在阈值θ1=-0.22、θ2=1.08、噪声强度σ=[0,4]条件下的实验结果。通过图8可以看到在每个浓度下都发生了随机共振,并且每种浓度下的随机共振是不同的,可以清晰的看到随着浓度的增加,最大互相关系数也相应增大。表1中,CCC和σ代表的是经过随机共振系统后最大互相关系数以及对应的噪声强度。图9是根据表1的数据得来,星点代表经过随机共振系统后的传感器响应的特征值,直线是通过回归分析和曲线拟合出浓度和最大互相关系数之间的关系。拟合曲线为:y=0.0028x+0.2051。其中,标准差SD=0.0126,相关系数R=0.9687。通过图9经过阈值随机共振系统后传感器的响应值呈现线性关系。拟合曲线的误差在可控范围内,可以认为浓度和最大互相关系数之间存在线性关系。
本发明对双阈值随机共振快速检测低浓度气体进行了研究,利用自制的气敏传感器得到3-60ppb的响应值,采用互相关系数对双阈值随机共振进行表征,并且讨论了双阈值对随机共振系统的影响以及如何选择最佳阈值。结果表明:当通过自动选取的最优双阈值时,最大互相关系数随着气体浓度的增大而增大,通过回归分析可以得到最大互相关系数和气体浓度之间的线性关系。双阈值随机共振在低浓度气体检测方面有比较明显的优势,这为检测低浓度气体提供了很好的方法,同时可以将此方法应用在气体生物传感器和电子鼻系统。
Claims (9)
1.一种基于双阈值随机共振的低浓度气体检测方法,其特征在于,包括以下步骤:
1)将气敏传感器置于不同浓度的参照气体所处的环境中,测得气敏传感器在各浓度下的电阻值;
2)以测得的电阻值作为信号输入,并对该信号进行归一化处理,再利用傅里叶级数对归一化处理后的信号进行曲线拟合,得到周期函数f(t);
3)对周期函数f(t)进行离散化处理得到输入函数S(t);
4)将输入函数S(t)和外加的高斯白噪声n(t)输入双阈值随机共振系统,利用双阈值检测器对输入函数S(t)和外加的高斯白噪声n(t)进行赋值,并以赋值后的信号作为输出信号y(t):
其中,θ1和θ2为双阈值随机共振系统的双阈值;
5)利用所述的输入函数S(t)和输出信号y(t),计算双阈值随机共振系统的互相关系数;
6)搜索最优的双阈值θ1和θ2,使得步骤5)中的互相关系数最大,并以最大的互相关系数和对应参照气体的浓度建立线性回归方程;
搜索最优双阈值θ1和θ2的步骤如下:
(1)对双阈值随机共振系统的双阈值θ1和θ2进行初始化,并确定阈值的搜索范围和步长;
(2)将输入函数S(t)和外加的高斯白噪声n(t)输入初始化后的双阈值随机共振系统中,计算每一组双阈值所对应的拟合曲线的标准差SD;
(3)在搜索范围内,根据步长进行遍历搜索,得到标准差SD的最小值,以及相对应的双阈值θ1、θ2;
7)选取待测的低浓度气体,执行步骤1)~步骤6)中的操作,得到待测低浓度气体对应的最大互相关系数,代入步骤6)中的线性回归方程,即求出待测低浓度气体的浓度。
2.如权利要求1所述的基于双阈值随机共振的低浓度气体检测方法,其特征在于,在步骤2)中,利用下式方程对所测的电阻值进行归一化处理,使其值域为[0,1];
f1(t)=(x(t)-xmin)/(xmax-xmin)
其中:x(t)为测量得到的电阻值,xmin为最小电阻值,xmax为最大电阻值。
3.如权利要求1或2所述的基于双阈值随机共振的低浓度气体检测方法,其特征在于,在步骤2)中,得到的周期T的函数f(t)为:
其中:a0、an、bn为函数系数,ancos nwt+bnsin nwt为n阶谐波,n=1,2,3,t为气体检测的时间。
4.如权利要求1所述的基于双阈值随机共振的低浓度气体检测方法,其特征在于,在步骤4)中,高斯白噪声n(t)的均值为0,方差为σ,概率密度函数为其中x为连续随机变量。
5.如权利要求1所述的基于双阈值随机共振的低浓度气体检测方法,其特征在于,在步骤5)中,互相关系数的表达式如下:
其中,CCC为互相关系数,S(t)为输入函数,y(t)为输出信号,为S(t)的平均值,为y(t)的平均值。
6.如权利要求1所述的基于双阈值随机共振的低浓度气体检测方法,其特征在于,线性拟合曲线与最大互相关系数之间的标准差SD为:
其中:N代表所有不同浓度气体的数量,CCCi代表的各浓度对应的最大互相关系数,yi代表的是各浓度经过曲线拟合后的值。
7.如权利要求6所述的基于双阈值随机共振的低浓度气体检测方法,其特征在于,阈值θ1的搜索范围为[-1.0,0.26],系统阈值θ2的搜索范围为[0.89,1.5],搜索步长均为0.02。
8.如权利要求7所述的基于双阈值随机共振的低浓度气体检测方法,其特征在于,在步骤6)中,建立的线性回归方程为y=0.0028x+0.2051,y为最大互相关系数,x为气体的浓度。
9.如权利要求1所述的基于双阈值随机共振的低浓度气体检测方法,其特征在于,在步骤1)中,首先向测试气体的气室中通入惰性气体,使气敏传感器的响应值达到稳定状态,再向气室中通入特定浓度的参照气体,测量3次取平均值,作为气敏传感器在此浓度下的电阻值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410277331.XA CN104132967B (zh) | 2014-06-20 | 2014-06-20 | 一种基于双阈值随机共振的低浓度气体检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410277331.XA CN104132967B (zh) | 2014-06-20 | 2014-06-20 | 一种基于双阈值随机共振的低浓度气体检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104132967A CN104132967A (zh) | 2014-11-05 |
CN104132967B true CN104132967B (zh) | 2016-12-07 |
Family
ID=51805730
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410277331.XA Active CN104132967B (zh) | 2014-06-20 | 2014-06-20 | 一种基于双阈值随机共振的低浓度气体检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104132967B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6669397B2 (ja) * | 2016-03-31 | 2020-03-18 | キヤノン株式会社 | 信号抽出処理装置および信号抽出処理方法 |
CN106645308B (zh) * | 2016-10-11 | 2020-09-29 | 南京工业大学 | 基于合金二硫化钼钨纳米片的丙酮气体传感器的制作方法 |
EP3698129A1 (en) * | 2017-10-19 | 2020-08-26 | Sanofi | Analyte measurement device |
CN110057973B (zh) * | 2019-04-01 | 2021-09-07 | 浙江工商大学 | 一种实验室用气安全预警系统及方法 |
CN110031155B (zh) * | 2019-04-01 | 2020-12-22 | 浙江工商大学 | 一种实验室气体泄漏检测装置及方法 |
CN110095570B (zh) * | 2019-04-01 | 2021-07-06 | 浙江工商大学 | 一种实验室安全预警系统及方法 |
CN110082480B (zh) * | 2019-04-01 | 2021-11-05 | 浙江工商大学 | 一种实验室安全管理系统及方法 |
CN111220571B (zh) * | 2020-01-20 | 2020-12-22 | 中南大学 | 一种基于幅值离散度的二次谐波信号拟合方法及系统 |
CN114636689B (zh) * | 2022-05-23 | 2022-08-16 | 武汉七斗光电科技有限公司 | 基于libs技术的藏药原矿成分定量检测方法和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0866333B1 (en) * | 1997-03-19 | 2001-11-14 | Ngk Insulators, Ltd. | Method of measuring gas concentration |
CN101672791A (zh) * | 2009-10-13 | 2010-03-17 | 北京化工大学 | 微量氨气气体快速检测管 |
CN102109487A (zh) * | 2009-12-28 | 2011-06-29 | 华瑞科学仪器(上海)有限公司 | 超低浓度气体传感器 |
CN102590288A (zh) * | 2012-01-17 | 2012-07-18 | 浙江工商大学 | 基于电子鼻的食品品质检测系统和检测方法 |
CN103698371A (zh) * | 2013-11-28 | 2014-04-02 | 浙江理工大学 | 一种基于自适应随机共振的低浓度气体检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0618465A (ja) * | 1992-07-02 | 1994-01-25 | Ricoh Co Ltd | 複合センサ |
-
2014
- 2014-06-20 CN CN201410277331.XA patent/CN104132967B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0866333B1 (en) * | 1997-03-19 | 2001-11-14 | Ngk Insulators, Ltd. | Method of measuring gas concentration |
CN101672791A (zh) * | 2009-10-13 | 2010-03-17 | 北京化工大学 | 微量氨气气体快速检测管 |
CN102109487A (zh) * | 2009-12-28 | 2011-06-29 | 华瑞科学仪器(上海)有限公司 | 超低浓度气体传感器 |
CN102590288A (zh) * | 2012-01-17 | 2012-07-18 | 浙江工商大学 | 基于电子鼻的食品品质检测系统和检测方法 |
CN103698371A (zh) * | 2013-11-28 | 2014-04-02 | 浙江理工大学 | 一种基于自适应随机共振的低浓度气体检测方法 |
Non-Patent Citations (1)
Title |
---|
"双阈值系统中随机共振的研究";杨伟等;《电子测量技术》;20061231;第29卷(第6期);第21页第1节 * |
Also Published As
Publication number | Publication date |
---|---|
CN104132967A (zh) | 2014-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104132967B (zh) | 一种基于双阈值随机共振的低浓度气体检测方法 | |
CN102866179B (zh) | 基于人工智能学习机的电子鼻中非目标干扰气味的在线识别和抑制方法 | |
CN103698371B (zh) | 一种基于自适应随机共振的低浓度气体检测方法 | |
CN105334185B (zh) | 基于光谱投影判别的近红外模型维护方法 | |
CN104251846A (zh) | 一种结合判别分析的激光诱导击穿光谱定量分析方法 | |
CN107958267B (zh) | 一种基于光谱线性表示的油品性质预测方法 | |
CN104596957A (zh) | 基于可见光近红外光谱技术的土壤铜含量估算方法 | |
CN103175805A (zh) | 一种近红外光谱测定污水中cod和bod5指标的方法 | |
CN107247034A (zh) | 一种基于吸收光谱技术的双频率波长调制方法 | |
CN102841069B (zh) | 一种由中红外光谱识别原油种类的方法 | |
CN107064177A (zh) | 基于自适应模糊推理的超宽带土壤信号含水量识别方法 | |
CN105300267A (zh) | 基于电磁涡流和数据融合算法的金属板厚度测量方法 | |
CN103196356A (zh) | 基于支持向量机的油管缺陷定量识别方法 | |
CN113095388B (zh) | 基于双层分类算法的固废塑料材质识别方法 | |
CN101650359A (zh) | 一种基于知识融合机器学习算法的煤热值软测量方法 | |
Xu et al. | Temperature and humidity compensation for MOS gas sensor based on random forests | |
CN113158558A (zh) | 一种高速铁路路基连续压实分析方法、装置及分析仪 | |
CN115656239A (zh) | 一种基于机器学习的重金属快速检测方法及系统 | |
Islam et al. | Multi-fault diagnosis of roller bearings using support vector machines with an improved decision strategy | |
CN107220490B (zh) | 水文序列均值跳跃变异的识别与分级方法 | |
Chen et al. | A mixed gas composition identification method based on sample augmentation | |
CN114821322A (zh) | 基于注意力机制的小样本遥感影像分类方法及系统 | |
Kothari et al. | An efficient scheme for water leakage detection using support vector machines (SVM)-Zig | |
Miao et al. | Soil pH value grey relation estimation model based on hyper-spectral | |
Di Lecce et al. | Discriminating gaseous emission patterns in low-cost sensor setups |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20190815 Address after: 222066 West Side Complex Building of Banqiao Industrial Park, Lianyungang Economic Development Zone, Lianyungang City, Jiangsu Province Patentee after: Lianyungang Eurasian Gas Co., Ltd. Address before: Hangzhou City, Zhejiang province 310018 Xiasha Higher Education Park No. 2 Street No. 5 Patentee before: Zhejiang Sci-Tech University |