CN113624710A - 一种液体物质成分变化的太赫兹信息识别方法及其应用 - Google Patents
一种液体物质成分变化的太赫兹信息识别方法及其应用 Download PDFInfo
- Publication number
- CN113624710A CN113624710A CN202111086796.3A CN202111086796A CN113624710A CN 113624710 A CN113624710 A CN 113624710A CN 202111086796 A CN202111086796 A CN 202111086796A CN 113624710 A CN113624710 A CN 113624710A
- Authority
- CN
- China
- Prior art keywords
- terahertz
- entropy
- weighted
- signals
- composite
- 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.)
- Pending
Links
- 239000000126 substance Substances 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000008859 change Effects 0.000 title claims abstract description 24
- 239000007788 liquid Substances 0.000 title claims abstract description 22
- 239000002131 composite material Substances 0.000 claims abstract description 45
- 239000012530 fluid Substances 0.000 claims abstract description 6
- 238000001514 detection method Methods 0.000 claims abstract description 4
- 239000000203 mixture Substances 0.000 claims abstract description 3
- 239000000758 substrate Substances 0.000 claims description 20
- 230000006870 function Effects 0.000 claims description 16
- 239000010453 quartz Substances 0.000 claims description 16
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 16
- 239000000463 material Substances 0.000 claims description 15
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 238000012935 Averaging Methods 0.000 claims description 4
- 239000010432 diamond Substances 0.000 claims description 3
- 229910003460 diamond Inorganic materials 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 239000002210 silicon-based material Substances 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 claims description 2
- 230000001678 irradiating effect Effects 0.000 claims description 2
- 102000004169 proteins and genes Human genes 0.000 claims description 2
- 108090000623 proteins and genes Proteins 0.000 claims description 2
- 230000010365 information processing Effects 0.000 abstract description 2
- 230000008447 perception Effects 0.000 abstract description 2
- 230000000975 bioactive effect Effects 0.000 abstract 1
- 239000003550 marker Substances 0.000 abstract 1
- 238000012544 monitoring process Methods 0.000 abstract 1
- 239000002904 solvent Substances 0.000 abstract 1
- 239000000243 solution Substances 0.000 description 26
- 108010074051 C-Reactive Protein Proteins 0.000 description 10
- 102100032752 C-reactive protein Human genes 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 5
- 230000003993 interaction Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 230000005855 radiation Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 229920002521 macromolecule Polymers 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012353 t test Methods 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 239000012460 protein solution Substances 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3577—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Toxicology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及智能感知信息处理技术领域，尤其涉及一种非接触、无电离破坏的液体物质成分变化的信息识别方法，包括：基于垂直向上反射式太赫兹检测方法获得待识别物质的太赫兹时域信号，进一步获取其复合加权尺度熵，并依据复合加权尺度熵的差异对待识别物质进行区分。本发明提出了一种非接触、无电离破坏的流体物质成分变化的信息获取与识别方法，该方法操作简便、信号重复性好，识别准确率高，适用范围广，可用于识别浓度不同的溶液、生物活性溶剂标志物识别或液体混合物稳定性监测等等。
Description
技术领域
本发明涉及智能感知信息处理技术领域，尤其涉及一种液体物质成分变化的太赫兹信息识别方法及其应用。
背景技术
太赫兹(THz)波是一种介于毫米波和红外之间的电磁波，其频率大致位于0.1～10THz范围内，许多大分子的振转能级落在太赫兹频段，利用太赫兹信号研究与大分子物质的相互作用规律具有更好的一一对应的指认性。THz波辐射能量低，大概在毫电子伏特量级，因而不会对被测物质造成破坏，是一种比较理想的无损检测手段。于是，近年来越来越多的学者使用太赫兹波进行物质检测与识别，特别是在生物、化学、医学等领域。流体或溶液状样品通常无明显特征不易直接通过原始信号区分。
现有技术中的太赫兹照射方式相位稳定性不佳，例如常规采用液体池的方式有多个界面，而且液体池的液面厚度通常较小以便使得太赫兹信号可以有效透过，会出现无法滤除的多波干扰叠加，液体池厚度通常在1mm一下，与太赫兹的波长范围比较接近，多界面后产生许多回波振荡。还有人采用向下照射的方式，此时在更换被测液体时，无法确保液面高度的精确一致，微小的变化导致太赫兹的微扰有可能淹没液体中微量成分的导致的信号变化。
熵是自然界中的一个基本物理量，可以用于表征***的复杂度。对于溶液***或类溶液***，大量分子的相互作用组成了一个复杂的整体体系，太赫兹信号与***发生相互作用时携带了***的信息。样本熵是通过计算条件概率的负对数来确定时间序列的复杂度，可以用于时间信号的特征提取，突破了常规的特征提取的限制，但其计算结果只有一个数值(即建立在单一尺度上)，有时并不能正确的表示长时间信号的变化规律。复合多尺度熵是通过考虑多个尺度的样本熵来量化信号的复杂度，但是默认为不同尺度的信号对样品的反应是一样的，仍然限制了识别的效果。
发明内容
为了解决现有技术和应用存在的问题，本发明提供一种非接触、无电离破坏的流体物质成分变化的太赫兹信息识别方法及其应用。
当物质中的成分发生微量变化，导致整体宏观***发生变化，从而导致观测信号的变化，本发明通过分析，获得微观变化对应的宏观量的变化，从而对微观信息进行有效判断。第一步是高相位精度的信号测试方式。采用高纯度熔融各向同性石英基底作为被测物的承载体；综合各种需求因素选择石英作为首选承载体物质经过大量筛选和对比，需要考虑太赫兹电磁波的透射率、反射率、材料成本、材料的稳定性、耐腐蚀性、硬度、多光谱的透光性、可视性等等，在一些特定的极端条件下，还可以选择金刚石、高阻硅材料作为承载体；被测物依靠自然重力置于基底上表面，并达到一定的厚度和覆盖面积；太赫兹信号从基底下表面垂直上下表面入射，透过下表面到达上表面，在上表面与被承载物发生相互作用，收集被原路反射回的信号用于数据分析。第二步是信号预处理。获得多个待识别物质的太赫兹时域信号；获得所述信号不同频域尺度或不同时域尺度的子波信号。第三是表征信号。获得所述子波信号的复合加权尺度熵；通过所述复合加权尺度熵的差异区分所述多个待识别物质；所述复合加权尺度熵为以太赫兹信号的振幅变化作为权重，对信号进行加权平均后计算样本熵得到。
基于上述原理，本发明提一种出液体物质成分变化的信息识别方法，包括：
将所述待识别物质置于基底上，从所述基底下方，以垂直于所述基底所在平面的方向照射太赫兹波，太赫兹波穿过所述基底，接触待测液体并反射回探测器，即得所述液体的太赫兹信号；
所述基底的材料为各向同性石英、金刚石或高阻硅材料中的一种或多种，优选为各向同性石英。
进一步地，所述基底的厚度大于0.5cm，优选为2～3cm。通常要求规避的石英下表面信号的振荡，厚度越厚，规避越好，但是同时也造成信号衰减增大、信号传输路径过长、波束变形影响提高等负面影响增强。
进一步地，所述获取所述太赫兹时域信号的复合加权尺度熵包括：
对所述太赫兹时域信号进行变分模式分解得到多组模函数；选择一组模函数作为所述太赫兹时域信号的子波信号，获取所述子波信号的复合加权尺度熵作为所述太赫兹时域信号的复合加权尺度熵。
进一步地，所述选择一组模函数作为所述太赫兹时域信号的子波信号包括：
基于模函数的频域和所述太赫兹时域信号的覆盖范围，结合所述待识别物质的物质状态先验信息，选择差异最大的一组模函数；
所述物质状态先验信息包括：预测的浓度范围、溶液的组成成分、环境条件稳定性、设备状态稳定性中的一种或多种。
进一步地，在获取所述子波信号的复合加权尺度熵的过程中，在得到多个加权粗粒化序列后，针对每一个待识别物质，计算将其所有加权粗粒化序列的样本熵，后按照相应的加权尺度之和的平均值作为其的复合加权尺度熵。
进一步地，所述变分模式分解中使用的二次惩罚因子α通过如下方法得到：
(1)将二次惩罚因子α在200～2200范围内寻优，从200开始，依次递增10遍历至2200；
(2)针对每个二次惩罚因子α，实施变分模式分解，结合物质状态先验信息，选择一组模函数作为子波信号，根据设定的尺度s将待识别物质的太赫兹时域信号划分为多个不重叠的窗口，计算每个窗内太赫兹时域信号的平均值作为粗粒化信号，得到s个粗粒化序列；
根据设定的加权尺度因子ss，将每一个粗粒化序列划分为多个不重叠的窗口，获得每个窗内相邻信号点的振幅变化率作为其自身的权重系数ω(由于窗内第一个点无法得到幅值变化故设置其权重为1)，并依据每个粗粒化信号的权重系数进行加权平均得到加权粗粒化信号，针对每一个粗粒化序列，得到ss个加权粗粒化序列；
计算每一个加权粗粒化信号的样本熵后，将ss个加权粗粒化序列的样本熵按照相应的加权尺度之和取平均值即得到复合加权尺度熵；
(3)确定出各尺度下的加权尺度的样本熵值的差异性，获得至少在一个加权尺度下熵值有差异的二次惩罚因子α。
进一步地，步骤(2)中，2≤s≤20，和/或，2≤ss≤20；
进一步地，步骤(2)中，所述计算每一个加权粗粒化信号的样本熵后，复合加权尺度熵可通过如下复合加权尺度熵公式计算得到：
其中，
为所述加权粗粒化序列，m为嵌入维度，r为相似度阈值，ω(k,j)为权重系数，1≤k≤s，1≤kk≤ss。
进一步地，步骤(2)中，所述选择一组模函数作为子波信号为：基于模函数的频域和所述太赫兹时域信号的覆盖范围，结合所述待识别物质的物质状态先验信息，选择差异最大的一组模函数。
进一步地，步骤(3)中，在获得至少在一个加权尺度下熵值有差异的二次惩罚因子α后，比较具有差异性的样本熵值的加权尺度个数，个数最多的对应的二次惩罚因子α即为最优的二次惩罚因子α。
进一步地，步骤(3)中，通过二值化图像获得确定出各加权尺度下的样本熵值的差异性。
进一步地，所述依据所述复合加权尺度熵的差异对所述待识别物质进行识别为，根据t检验，确定不同待识别物质在时间尺度上的统计学差异，以区分不同待识别物质。
本发明进一步提供所述信息识别方法在鉴别流体或溶液样品的应用。
进一步地，所述溶液样品为不同浓度的同种蛋白样品。
进一步地，所述蛋白样品为C-反应蛋白。
本发明提供一种向上垂直入射获取高度一致性太赫兹原始信号的方法，使得精细的分析高复杂度流体或溶液***成为可能。采用向上垂直入射才能够保证在基底材料上界面与被测目标发生相互作用时，时域太赫兹信号的相位一致性非常的高，而其他的太赫兹照射方式通常无法保证如此高的相位稳定性。
此外，本发明不仅在计算复合尺度熵时，考虑对信号进行粗粒化处理(即对划分的每个区间内的点数取平均值从而得到相应的序列)，而且在粗粒化的基础上进一步添加了加权平均处理过程，提出一种复合加权尺度熵算法，根据不同尺度以及加权尺度下的样本熵对性质相似的物质的太赫兹信号进行区分，有效提升了识别的可靠性。在获得高度稳定的原始信号后，通过变分模式分解的方法，将原始信号变分为不同时域分量或频域分量的子波，根据不同的被测目标，高效的相互作用区间通常是不同的，因此选区合适分量的子波可以更好的识别出痕量物质对***的细微影响。分解后的子波受参数α取值的影响，经过优化提取出特定子波后，计算子波信号的复合加权尺度熵，用于对性质相似的物质的太赫兹信号进行区分。
本发明具备如下有益效果：
本发明使用基于各向同性石英基底的垂直向上反射式太赫兹测量方法，获得高质量原始信号，针对不同物质的太赫兹时域信号的差异，在原有时间序列的复合多尺度熵的基础上，通过模式分析获得优化的子波信号，并计算其复合加权尺度熵，通过计算复合加权尺度熵进行不同物质的区分，施加的权重系数结合太赫兹信号的幅值和时延特性的变化并增加了相应的权重尺度，从而可以获得更多的样本熵，提升微小差异的识别效果，这在物质识别领域有重要意义，例如实现不同浓度蛋白溶液的区分。
附图说明
图1为本发明所使用的基于各向同性石英基底的垂直向上反射式太赫兹测量装置和方法的示意图。
图2为本发明实施例1提供的不同浓度C-反应蛋白校准品溶液的太赫兹时域信号；其中a为浓度为19.7mg/L 52.2mg/L的C-反应蛋白校准品溶液的太赫兹时域信号差异，b为浓度为19.7mg/L与143mg/L的C-反应蛋白校准品溶液的太赫兹时域信号差异，c为浓度为19.7mg/L与281mg/L的C-反应蛋白校准品溶液。
图3为本发明实施例1提供的通过幅值变化对粗粒化序列进行加权的示意图。
图4为本发明实施例1提供的浓度为19.7mg/L 52.2mg/L的C-反应蛋白校准品溶液的复合尺度熵和复合加权尺度熵的计算结果；其中a为α＝139时第一模态子波信号的复合加权尺度熵p值分布二值图像；b为α＝139时第一模态子波信号在第2尺度下的复合多尺度熵曲线。
图5为本发明实施例1提供的浓度为19.7mg/L 143mg/L的C-反应蛋白校准品溶液的复合加权尺度熵的计算结果；其中a为α＝9时第一模态子波信号的复合加权尺度熵p值分布二值图像；b为α＝9时第一模态子波信号在第3尺度下的复合多尺度熵曲线。
图6为本发明实施例1提供的浓度为19.7mg/L 281mg/L的C-反应蛋白校准品溶液的复合加权尺度熵的计算结果；其中a为α＝169时第一模态子波信号的复合加权尺度熵p值分布二值图像；b为α＝169时第一模态子波信号在第1尺度下的复合多尺度熵曲线。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
本实施例提供一种基于太赫兹信号复合加权尺度熵的不同浓度的C-反应蛋白校准品溶液的识别方法，具体包括如下流程：
1、选取第一组的两种浓度(19.7和52.2mg/L)的C-反应蛋白校准品溶液(购于德赛诊断***有限公司)进行测试。
2、对于19.7mg/L溶液，使用一次性注射器抽取1ml，轻轻滴放至各向同性石英片上，并采集太赫兹时域信号。采集之后，用酒精和清水清洗各向同性石英片。然后再使用不同的注射器依次提取1ml该溶液采集相应的信号。总共测试3次。同理对于52.2mg/L的溶液采取同样的方式进行测试。
该步骤具体可参照图1所示流程，具体可以将待测液体滴放在各向同性石英基底上，各向同性石英对于太赫兹波具有较高的透过率。太赫兹波垂直向上穿过各向同性石英，接触待测液体并反射回探测器，获得所述液体的太赫兹信号。采用向上垂直入射才能够保证在基底材料上界面与被测目标发生相互作用时，时域太赫兹信号的相位一致性非常的高，而其他的太赫兹照射方式通常无法保证如此高的相位稳定性。各向同性石英厚度为2cm，可以有效去除各向同性石英下表面的反射信号对于样品反射信号的干扰。
3、依次选取第二组(52.2与281mg/L)、第三组(143与281mg/L)的溶液并按照步骤(2)进行测试，图2表示每种浓度溶液的平均信号。
4、计算以上信号的复合加权尺度熵的数值。
(1)利用变分模式分解将样品的太赫兹时域信号进行分解，取第一模态函数作为太赫兹信号频域或时域的子波信号，将二次惩罚因子α在200～2200范围内寻优，从200开始，依次递增10遍历至2200，故每个样品信号共分解201个子波信号。
(2)根据设定的尺度s将上述分解子波信号{xi}划分为若干个不重叠的窗口(窗口的个数由尺度数s确定)，计算窗内点数的平均值，各窗内的平均值所构成的序列(粗粒化序列)为：
(2)将上述序列再次划分为若干不重叠的窗口(其数量依赖于设置的加权尺度)，取窗口内相邻两点间的幅值变化率作为权重系数ω(由于窗内第一个点无法得到幅值变化故设置其权重为1)，得到带权重的序列：
图3表示对该算法权重尺度序列说明(尺度和加权尺度均取2)。
(3)根据下面的公式对上述序列计算样本熵，得到复合加权尺度熵。
其中，m为嵌入维度，r为相似度阈值，一般情况下，m＝2，r＝0.15σ，σ为该时间信号的标准差。
最终得到的结果是在不同尺度和不同加权尺度下的熵值。
5、根据t检验，确定出每两组浓度溶液在相应尺度下样本熵值的统计学差异性，得到加权尺度下熵值差异性较大的二次惩罚因子α，并得到具有差异性的熵值，从而达到识别的目的。
步骤(2)、(3)中每种浓度溶液测试3次并采集相应的信号，所需测试量为1ml。
步骤(4)中复合加权尺度熵的参数设置：尺度s与加权尺度ss均设置为20，嵌入维度为2，阈值r＝0.15σ，σ为时间信号的标准差。
步骤(5)中，对于复合加权尺度熵，本实施例采用二值化图像来表示识别结果(以p＝0.05作为阈值)。对于一维曲线结果的表示，本实施例选取差异性熵值最多的某一尺度，得到该尺度下的不同加权尺度的熵值分布曲线。
从图4～6的计算结果来看，复合加权尺度熵具有更多统计学差异性的熵值。由此说明该方法可以对太赫兹信号相似的物质进行有效区分，该方法在相关领域进行相似的物质识别具有重要的价值。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种液体物质成分变化的太赫兹信息识别方法，其特征在于，包括：
基于垂直向上反射式的太赫兹检测方法获得待识别物质的太赫兹时域信号；获取所述太赫兹时域信号的复合加权尺度熵，并依据所述复合加权尺度熵的差异对所述待识别物质进行识别。
2.根据权利要求2所述的太赫兹信息识别方法，其特征在于，所述垂直向上反射式的太赫兹检测方法为：
将所述待识别物质置于基底上，从所述基底下方，以垂直于所述基底所在平面的方向照射太赫兹波，太赫兹波穿过所述基底，接触待测液体并反射回探测器，即得所述液体的太赫兹信号；
所述基底的材料为各向同性石英、金刚石或高阻硅材料中的一种或多种，优选为各向同性石英。
3.根据权利要求1或2所述的太赫兹信息识别方法，其特征在于，所述基底的厚度大于0.5cm，优选为2～3cm。
4.根据权利要求1-3任一项所述的太赫兹信息识别方法，其特征在于，所述获取所述太赫兹时域信号的复合加权尺度熵包括：
对所述太赫兹时域信号进行变分模式分解得到多组模函数；选择一组模函数作为所述太赫兹时域信号的子波信号，获取所述子波信号的复合加权尺度熵作为所述太赫兹时域信号的复合加权尺度熵。
5.根据权利要求4所述的太赫兹信息识别方法，其特征在于，所述选择一组模函数作为所述太赫兹时域信号的子波信号包括：
基于模函数的频域和所述太赫兹时域信号的覆盖范围，结合所述待识别物质的物质状态先验信息，选择差异最大的一组模函数；
所述物质状态先验信息包括：预测的浓度范围、溶液的组成成分、环境条件稳定性、设备状态稳定性中的一种或多种。
6.根据权利要求4或5任一项所述的太赫兹信息识别方法，在获取所述子波信号的复合加权尺度熵的过程中，在得到多个加权粗粒化序列后，针对每一个待识别物质，计算将其所有加权粗粒化序列的样本熵，后按照相应的加权尺度之和的平均值作为其的复合加权尺度熵。
7.根据权利要求4-6任一项所述的太赫兹信息识别方法，其特征在于，所述变分模式分解中使用的二次惩罚因子α通过如下方法得到：
(1)将二次惩罚因子α在200～2200范围内寻优，从200开始，依次递增10遍历至2200；
(2)针对每个二次惩罚因子α，实施变分模式分解，结合物质状态先验信息，选择一组模函数作为子波信号，根据设定的尺度s将待识别物质的太赫兹时域信号划分为多个不重叠的窗口，计算每个窗内太赫兹时域信号的平均值作为粗粒化信号，得到s个粗粒化序列；
根据设定的加权尺度因子ss，将每一个粗粒化序列划分为多个不重叠的窗口，获得每个窗内相邻信号点的振幅变化率作为其自身的权重系数ω(由于窗内第一个点无法得到幅值变化故设置其权重为1)，并依据每个粗粒化信号的权重系数进行加权平均得到加权粗粒化信号，针对每一个粗粒化序列，得到ss个加权粗粒化序列；
计算每一个加权粗粒化信号的样本熵后，将ss个加权粗粒化序列的样本熵按照相应的加权尺度之和取平均值即得到复合加权尺度熵；
(3)确定出各尺度下的加权尺度的样本熵值的差异性，获得至少在一个加权尺度下熵值有差异的二次惩罚因子α。
8.根据权利要求7所述的太赫兹信息识别方法，其特征在于，
步骤(2)中，2≤s≤20，和/或，2≤ss≤20；和/或，
步骤(3)中，在获得至少在一个加权尺度下熵值有差异的二次惩罚因子α后，比较具有差异性的样本熵值的加权尺度个数，个数最多的对应的二次惩罚因子α即为最优的二次惩罚因子α。
9.权利要求1-8任一项所述太赫兹信息识别方法在鉴别流体或溶液样品的应用。
10.根据权利要求9所述的应用，其特征在于，所述溶液样品为不同浓度的同种蛋白样品。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111086796.3A CN113624710A (zh) | 2021-09-16 | 2021-09-16 | 一种液体物质成分变化的太赫兹信息识别方法及其应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111086796.3A CN113624710A (zh) | 2021-09-16 | 2021-09-16 | 一种液体物质成分变化的太赫兹信息识别方法及其应用 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113624710A true CN113624710A (zh) | 2021-11-09 |
Family
ID=78390216
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111086796.3A Pending CN113624710A (zh) | 2021-09-16 | 2021-09-16 | 一种液体物质成分变化的太赫兹信息识别方法及其应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113624710A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116908148A (zh) * | 2023-09-12 | 2023-10-20 | 中国科学技术大学 | 一种基于复杂度分析的光谱检测控制***及方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070073115A1 (en) * | 2005-09-23 | 2007-03-29 | Samsung Electronics Co., Ltd. | Apparatus and method for measuring concentration of blood component using terahertz |
| JP2009019925A (ja) * | 2007-07-10 | 2009-01-29 | Iwate Prefectural Univ | 分光測定試料、分光測定基板、及び、分光測定方法 |
| CN104215776A (zh) * | 2014-09-15 | 2014-12-17 | 深圳大学 | 一种基于太赫兹时域光谱的无标记血凝素蛋白检测方法 |
| US20150193376A1 (en) * | 2014-01-06 | 2015-07-09 | National Central University | Method of multi-scales intrinsic entropy analysis |
| CN107468250A (zh) * | 2017-08-21 | 2017-12-15 | 北京理工大学 | 基于多尺度熵的生物组织太赫兹成像方法、***和设备 |
| CN113155773A (zh) * | 2021-04-07 | 2021-07-23 | 中国科学院重庆绿色智能技术研究院 | 一种利用太赫兹光谱技术检验液体中标志物蛋白的*** |
-
2021
- 2021-09-16 CN CN202111086796.3A patent/CN113624710A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070073115A1 (en) * | 2005-09-23 | 2007-03-29 | Samsung Electronics Co., Ltd. | Apparatus and method for measuring concentration of blood component using terahertz |
| JP2009019925A (ja) * | 2007-07-10 | 2009-01-29 | Iwate Prefectural Univ | 分光測定試料、分光測定基板、及び、分光測定方法 |
| US20150193376A1 (en) * | 2014-01-06 | 2015-07-09 | National Central University | Method of multi-scales intrinsic entropy analysis |
| CN104215776A (zh) * | 2014-09-15 | 2014-12-17 | 深圳大学 | 一种基于太赫兹时域光谱的无标记血凝素蛋白检测方法 |
| CN107468250A (zh) * | 2017-08-21 | 2017-12-15 | 北京理工大学 | 基于多尺度熵的生物组织太赫兹成像方法、***和设备 |
| CN113155773A (zh) * | 2021-04-07 | 2021-07-23 | 中国科学院重庆绿色智能技术研究院 | 一种利用太赫兹光谱技术检验液体中标志物蛋白的*** |
Non-Patent Citations (2)
| Title |
|---|
| HAISHUN LIU等: "Diagnosis of hepatocellular carcinoma based on a terahertz signal and VMD-CWSE", 《BIOMEDICAL OPTICS EXPRESS》, vol. 11, no. 9, pages 1 - 10 * |
| 徐德刚等: "《光学太赫兹辐射源及其生物医学应用》", 30 April 2021, 华东理工大学出版社, pages: 6 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116908148A (zh) * | 2023-09-12 | 2023-10-20 | 中国科学技术大学 | 一种基于复杂度分析的光谱检测控制***及方法 |
| CN116908148B (zh) * | 2023-09-12 | 2024-01-09 | 中国科学技术大学 | 一种基于复杂度分析的光谱检测控制***及方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6421121B1 (en) | Method and apparatus for rapid particle identification utilizing scattered light histograms | |
| US7072038B2 (en) | Particle ID with narrow angle detectors | |
| US10379027B2 (en) | Method for identifying blood particles using a photodetector | |
| CN109581546B (zh) | 基于微波链路衰减及偏振信息的降雨类型识别方法 | |
| EP2373957B1 (en) | Methods for determining a liquid level in a container using imaging | |
| Li et al. | Determination of non-spherical particle size distribution from chord length measurements. Part 2: Experimental validation | |
| Suzuki et al. | The evaluation of macroscopic and microscopic textures of sand grains using elliptic Fourier and principal component analysis: implications for the discrimination of sedimentary environments | |
| US20170322137A1 (en) | Method and system for characterizing particles using a flow cytometer | |
| CN101981446A (zh) | 用于使用支持向量机分析流式细胞术数据的方法和*** | |
| Richardson et al. | Bragg scatter detection by the WSR-88D. Part I: Algorithm development | |
| CN113624710A (zh) | 一种液体物质成分变化的太赫兹信息识别方法及其应用 | |
| CN116399850B (zh) | 一种用于光信号处理的光谱探测识别***及其探测方法 | |
| CN111366600B (zh) | 一种检测水果甜度的方法、装置、***及存储介质 | |
| CN114578040B (zh) | 样本检测方法、样本分析仪及其控制装置 | |
| CN110567900A (zh) | 试样反应中抗原过量的判断方法、装置及光学检测*** | |
| EP3244191A1 (en) | Method and system for characterizing particles using a flow cytometer | |
| Wang et al. | Nondestructive detecting cracks of preserved eggshell using polarization technology and cluster analysis | |
| CN114088817A (zh) | 基于深层特征的深度学***板陶瓷膜超声缺陷检测方法 | |
| CN113406644A (zh) | 天气雷达数据质量控制方法、装置及设备 | |
| Maffezzoli et al. | Detection of ice core particles via deep neural networks | |
| CN113624759A (zh) | 一种基于机器学习的苹果无损检测方法 | |
| RU2498275C2 (ru) | Дистанционный способ классификации нефтяных загрязнений на поверхности воды | |
| CN1469111A (zh) | 采用散射光直方图快速辨别颗粒的方法及设备 | |
| CN114720541B (zh) | 一种提高红花椒麻感强度分级准确性的方法 | |
| CN110276282B (zh) | 危险液体识别方法、装置及计算机可读存储介质 |
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 |