CN109632680A - 一种基于排列熵的水体中磷的检测方法 - Google Patents
一种基于排列熵的水体中磷的检测方法 Download PDFInfo
- Publication number
- CN109632680A CN109632680A CN201811586350.5A CN201811586350A CN109632680A CN 109632680 A CN109632680 A CN 109632680A CN 201811586350 A CN201811586350 A CN 201811586350A CN 109632680 A CN109632680 A CN 109632680A
- Authority
- CN
- China
- Prior art keywords
- arrangement
- subsequence
- entropy
- phosphorus
- absorbance
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 39
- 238000001514 detection method Methods 0.000 title claims abstract description 27
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 title claims abstract description 26
- 229910052698 phosphorus Inorganic materials 0.000 title claims abstract description 26
- 239000011574 phosphorus Substances 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 64
- 238000002835 absorbance Methods 0.000 claims abstract description 32
- 229910000402 monopotassium phosphate Inorganic materials 0.000 claims abstract description 29
- 235000019796 monopotassium phosphate Nutrition 0.000 claims abstract description 29
- PJNZPQUBCPKICU-UHFFFAOYSA-N phosphoric acid;potassium Chemical compound [K].OP(O)(O)=O PJNZPQUBCPKICU-UHFFFAOYSA-N 0.000 claims abstract description 29
- 230000000694 effects Effects 0.000 claims abstract description 21
- 230000003595 spectral effect Effects 0.000 claims abstract description 14
- 238000002211 ultraviolet spectrum Methods 0.000 claims abstract description 12
- 150000001875 compounds Chemical class 0.000 claims abstract description 9
- 238000012360 testing method Methods 0.000 claims abstract description 9
- 238000001228 spectrum Methods 0.000 claims description 9
- 230000007774 longterm Effects 0.000 claims description 8
- 238000005259 measurement Methods 0.000 abstract description 9
- 239000012535 impurity Substances 0.000 abstract description 7
- 239000003153 chemical reaction reagent Substances 0.000 abstract description 6
- 238000012625 in-situ measurement Methods 0.000 abstract description 4
- 239000000243 solution Substances 0.000 description 9
- 238000002474 experimental method Methods 0.000 description 6
- CIWBSHSKHKDKBQ-JLAZNSOCSA-N Ascorbic acid Chemical compound OC[C@H](O)[C@H]1OC(=O)C(O)=C1O CIWBSHSKHKDKBQ-JLAZNSOCSA-N 0.000 description 4
- NBIIXXVUZAFLBC-UHFFFAOYSA-N Phosphoric acid Chemical compound OP(O)(O)=O NBIIXXVUZAFLBC-UHFFFAOYSA-N 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 235000013399 edible fruits Nutrition 0.000 description 4
- 238000002798 spectrophotometry method Methods 0.000 description 4
- 230000001174 ascending effect Effects 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 239000012086 standard solution Substances 0.000 description 3
- ZOKXTWBITQBERF-UHFFFAOYSA-N Molybdenum Chemical compound [Mo] ZOKXTWBITQBERF-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000000862 absorption spectrum Methods 0.000 description 2
- APUPEJJSWDHEBO-UHFFFAOYSA-P ammonium molybdate Chemical compound [NH4+].[NH4+].[O-][Mo]([O-])(=O)=O APUPEJJSWDHEBO-UHFFFAOYSA-P 0.000 description 2
- 229940010552 ammonium molybdate Drugs 0.000 description 2
- 235000018660 ammonium molybdate Nutrition 0.000 description 2
- 239000011609 ammonium molybdate Substances 0.000 description 2
- 229960005070 ascorbic acid Drugs 0.000 description 2
- 235000010323 ascorbic acid Nutrition 0.000 description 2
- 239000011668 ascorbic acid Substances 0.000 description 2
- 239000003638 chemical reducing agent Substances 0.000 description 2
- QXYJCZRRLLQGCR-UHFFFAOYSA-N dioxomolybdenum Chemical compound O=[Mo]=O QXYJCZRRLLQGCR-UHFFFAOYSA-N 0.000 description 2
- 229910052750 molybdenum Inorganic materials 0.000 description 2
- 239000011733 molybdenum Substances 0.000 description 2
- 235000011007 phosphoric acid Nutrition 0.000 description 2
- USHAGKDGDHPEEY-UHFFFAOYSA-L potassium persulfate Chemical compound [K+].[K+].[O-]S(=O)(=O)OOS([O-])(=O)=O USHAGKDGDHPEEY-UHFFFAOYSA-L 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- KQROHCSYOGBQGJ-UHFFFAOYSA-N 5-Hydroxytryptophol Chemical compound C1=C(O)C=C2C(CCO)=CNC2=C1 KQROHCSYOGBQGJ-UHFFFAOYSA-N 0.000 description 1
- 229910019142 PO4 Inorganic materials 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 229910000147 aluminium phosphate Inorganic materials 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 239000008367 deionised water Substances 0.000 description 1
- 229910021641 deionized water Inorganic materials 0.000 description 1
- 239000011964 heteropoly acid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- MEFBJEMVZONFCJ-UHFFFAOYSA-N molybdate Chemical compound [O-][Mo]([O-])(=O)=O MEFBJEMVZONFCJ-UHFFFAOYSA-N 0.000 description 1
- VLAPMBHFAWRUQP-UHFFFAOYSA-L molybdic acid Chemical compound O[Mo](O)(=O)=O VLAPMBHFAWRUQP-UHFFFAOYSA-L 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- AMWVZPDSWLOFKA-UHFFFAOYSA-N phosphanylidynemolybdenum Chemical compound [Mo]#P AMWVZPDSWLOFKA-UHFFFAOYSA-N 0.000 description 1
- NBIIXXVUZAFLBC-UHFFFAOYSA-K phosphate Chemical compound [O-]P([O-])([O-])=O NBIIXXVUZAFLBC-UHFFFAOYSA-K 0.000 description 1
- 239000010452 phosphate Substances 0.000 description 1
- 150000003016 phosphoric acids Chemical class 0.000 description 1
- IIQJBVZYLIIMND-UHFFFAOYSA-J potassium;antimony(3+);2,3-dihydroxybutanedioate Chemical compound [K+].[Sb+3].[O-]C(=O)C(O)C(O)C([O-])=O.[O-]C(=O)C(O)C(O)C([O-])=O IIQJBVZYLIIMND-UHFFFAOYSA-J 0.000 description 1
- 239000000376 reactant Substances 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/33—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using ultraviolet light
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)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Investigating Or Analysing Materials By The Use Of Chemical Reactions (AREA)
- Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
Abstract
本发明提供了一种基于排列熵的水体中磷的检测方法,包括以下步骤:步骤1:取水体样本直接配制浓度为1‑5mg/L的磷酸二氢钾溶液;步骤2:通过紫外光谱测量其吸光度;步骤3:记录下其不同浓度对应的吸光度;步骤4:通过排列熵算法对步骤3中得到的数据进行处理,采用matlab计算的方法不断测试排列熵的值,通过得到的效果图最佳为依据求出最佳嵌入维数和延迟时间。本发明无需添加任何化学试剂,减少二次污染;检测操作简单化、自动化,易于实现原位测量;扩大了测量范围,精度与磷钼蓝方法基本保持一致;排列熵算法减小了光谱检测中相关杂质对水样的影响,并放大了信号的微小变化,更适合检测低浓度的溶液。
Description
技术领域
本发明涉及一种水体中磷的检测方法,特别涉及一种基于排列熵的水体中磷的检测方法。
背景技术
现有水质中检测总磷含量,常用的分光光度法主要为磷钼蓝法,该方法法是常规的化学实验室分析方法,磷钼蓝法的原理是正磷酸先与钼酸盐反应生成淡黄色的磷钼杂多酸,然后再用还原剂在一定条件下将黄色的磷钼杂多酸还原为磷钼蓝络合物。反应式如下:
PO43-+12MoO42-+24H+=PO43-·12MoO3+12H2O
P(Mo3O10)43-+H++还原剂→(MoO2·4MoO3)2·H3PO4+MoO2+H2总磷的测定方法就是用钼酸铵和酒石酸锑钾与磷酸形成三元配合物(俗称磷钼蓝),在波长710nm处测量吸光度,适测范围是0-0.6mg/L,通过绘制吸光度曲线,可以定量分析水样中总磷的含量,需要添加多种化学试剂,检测操作复杂,测量范围小,相关杂质对水样的影响大。
排列熵算法是时间序列复杂性的一种度量方法。算法描述如下:
设一维时间序列:{x(i),i=1,2,…,n},长度为n,对其中任意一个元素x(i)进行相空间重构,根据Takens提出的相空间重构延迟坐标法,即对每个采样点取连续的m个样点,进行坐标延迟,得到m维空间的一个重构向量组:
Xi=[x(i),x(i+l),…,x(i+(m-1)L)]式中:m和L分别为重构参数嵌入维数和延迟时间。再将Xi的m个重构分量进行升序排列,得到:
[x(i+(j1-1)L)≤x(i+(j2-1)L)≤…≤x(i+(jm-1)L)]
这样得到向量Xi的排列方式为{j1,j2,…,jm},它是m!种排列方式中的一种方式,对整个序列中各种排列情况出现的次数进行统计,计算各种排列情况出现的相对频率作为其概率P1,P2,…,Pk,k≤m!,接着由Shannon信息熵的定义计算排列熵:
排列熵H(m)值的大小表示了时间序列{x(i),i=1,2,…,n}的随机程度:H(m)的值越小,说明时间序列越规则;反之,则时间序列越接近随机。H(m)的变化反映并放大了时间序列x(i)的微小细节变化。
发明内容
本发明的目的是解决现有水质中检测总磷含量的方法,存在需要添加多种化学试剂,检测操作复杂,测量范围小,相关杂质对水样的影响大等问题,提供了一种基于排列熵的水体中磷的检测方法,与磷钼蓝分光光度法对比,其优势在于无需添加任何化学试剂,减少二次污染;检测操作简单化、自动化,易于实现原位测量;磷钼蓝检测方法的吸光度与浓度的线性关系为0-0.6mg/L,而基于排列熵算法的检测方法,通过熵值的变化可以得到不同浓度的变化,扩大了测量范围,精度与磷钼蓝方法基本保持一致;排列熵算法减小了光谱检测中相关杂质对水样的影响,并放大了信号的微小变化,更适合检测低浓度的溶液。
本发明提供的技术方案为:
一种基于排列熵的水体中磷的检测方法,包括以下步骤:
步骤1:取水体样本直接配制浓度为1-5mg/L的磷酸二氢钾溶液;
步骤2:通过紫外光谱测量其吸光度;
步骤3:记录下其不同浓度对应的吸光度;
步骤4:通过排列熵算法对步骤3中得到的数据进行处理:基于排列熵算法的具体过程为将长度为T的光谱信号长时间序列分为若干个长度为w的子序列,选取子序列最大重叠序列,即每个子序列向后移动一个数据点得到下一个子序列,这样得到的子序列个数为T-w+1,采用matlab计算的方法不断测试排列熵的值,通过得到的效果图最佳为依据求出最佳嵌入维数和延迟时间。
优选的是,所述紫外光谱为TU1810紫外可见分光光度计。
优选的是,所述步骤1中配制的磷酸二氢钾溶液选取浓度为1mg/L、2mg/L、3mg/L、4mg/L和5mg/L经步骤2、步骤3,将得到的数据经步骤4中对所得数据进行处理,运用matlab仿真进行排列熵计算。
优选的是,步骤4中得出最佳嵌入维数m=5,延迟时间L=1,子序列长度选取为w=100。
排列熵算法是时间序列复杂性的一种度量方法。算法描述如下:
设一维时间序列:{x(i),i=1,2,…,n},长度为n,对其中任意一个元素x(i)进行相空间重构,根据Takens提出的相空间重构延迟坐标法,即对每个采样点取连续的m个样点,进行坐标延迟,得到m维空间的一个重构向量组:
Xi=[x(i),x(i+l),…,x(i+(m-1)L)]式中:m和L分别为重构参数嵌入维数和延迟时间。再将Xi的m个重构分量进行升序排列,得到:
[x(i+(j1-1)L)≤x(i+(j2-1)L)≤…≤x(i+(jm-1)L)]
这样得到向量Xi的排列方式为{j1,j2,…,jm},它是m!种排列方式中的一种方式,对整个序列中各种排列情况出现的次数进行统计,计算各种排列情况出现的相对频率作为其概率P1,P2,…,Pk,k≤m!,接着由Shannon信息熵的定义计算排列熵:
排列熵H(m)值的大小表示了时间序列{x(i),i=1,2,…,n}的随机程度:H(m)的值越小,说明时间序列越规则;反之,则时间序列越接近随机。H(m)的变化反映并放大了时间序列x(i)的微小细节变化。
本发明的有益效果为:
1、无需添加任何化学试剂,减少二次污染。
2、检测操作简单化、自动化,易于实现原位测量。
3、磷钼蓝检测方法的吸光度与浓度的线性关系为0-0.6mg/L,而基于排列熵算法的检测方法,通过熵值的变化可以得到不同浓度的变化,扩大了测量范围,精度与磷钼蓝方法基本保持一致。
4、排列熵算法减小了光谱检测中相关杂质对水样的影响,并放大了信号的微小变化,更适合检测低浓度的溶液。
附图说明
图1为不同浓度的磷酸二氢钾溶液对应的吸光度。
图2为1mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图3为2mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图4为3mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图5为4mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图6为5mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图7为m=4,w=100时1mg/L磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图8为m=8,w=100时1mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图9为m=10,w=100时1mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图10为m=5,w=50时1mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图11为m=5,w=80时1mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
图12为m=5,w=50时1mg/L的磷酸二氢钾溶液运用matlab仿真进行排列熵计算效果图。
具体实施方式
下面对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
实施例1
一种基于排列熵的水体中磷的检测方法,包括以下步骤:
步骤1:取水体样本直接配制浓度为1mg/L的磷酸二氢钾溶液;
步骤2:通过紫外光谱测量其吸光度;
步骤3:记录下其对应的吸光度;
步骤4:通过排列熵算法对步骤3中得到的数据进行处理:基于排列熵算法的具体过程为将长度为T的光谱信号长时间序列分为若干个长度为w的子序列,选取子序列最大重叠序列,即每个子序列向后移动一个数据点得到下一个子序列,这样得到的子序列个数为T-w+1,采用matlab计算的方法不断测试排列熵的值,通过得到的效果图最佳为依据求出最佳嵌入维数和延迟时间,嵌入维数m=5,延迟时间L=1,子序列长度选取为w=100。所述紫外光谱为TU1810紫外可见分光光度计。
实施例2
一种基于排列熵的水体中磷的检测方法,包括以下步骤:
步骤1:取水体样本直接配制浓度为2mg/L的磷酸二氢钾溶液;
步骤2:通过紫外光谱测量其吸光度;
步骤3:记录下其对应的吸光度;
步骤4:通过排列熵算法对步骤3中得到的数据进行处理:基于排列熵算法的具体过程为将长度为T的光谱信号长时间序列分为若干个长度为w的子序列,选取子序列最大重叠序列,即每个子序列向后移动一个数据点得到下一个子序列,这样得到的子序列个数为T-w+1,采用matlab计算的方法不断测试排列熵的值,通过得到的效果图最佳为依据求出最佳嵌入维数和延迟时间,嵌入维数m=5,延迟时间L=1,子序列长度选取为w=100。所述紫外光谱为TU1810紫外可见分光光度计。
实施例3
一种基于排列熵的水体中磷的检测方法,包括以下步骤:
步骤1:取水体样本直接配制浓度为3mg/L的磷酸二氢钾溶液;
步骤2:通过紫外光谱测量其吸光度;
步骤3:记录下其对应的吸光度;
步骤4:通过排列熵算法对步骤3中得到的数据进行处理:基于排列熵算法的具体过程为将长度为T的光谱信号长时间序列分为若干个长度为w的子序列,选取子序列最大重叠序列,即每个子序列向后移动一个数据点得到下一个子序列,这样得到的子序列个数为T-w+1,采用matlab计算的方法不断测试排列熵的值,通过得到的效果图最佳为依据求出最佳嵌入维数和延迟时间,嵌入维数m=5,延迟时间L=1,子序列长度选取为w=100。所述紫外光谱为TU1810紫外可见分光光度计。
实施例4
一种基于排列熵的水体中磷的检测方法,包括以下步骤:
步骤1:取水体样本直接配制浓度为4mg/L的磷酸二氢钾溶液;
步骤2:通过紫外光谱测量其吸光度;
步骤3:记录下其对应的吸光度;
步骤4:通过排列熵算法对步骤3中得到的数据进行处理:基于排列熵算法的具体过程为将长度为T的光谱信号长时间序列分为若干个长度为w的子序列,选取子序列最大重叠序列,即每个子序列向后移动一个数据点得到下一个子序列,这样得到的子序列个数为T-w+1,采用matlab计算的方法不断测试排列熵的值,通过得到的效果图最佳为依据求出最佳嵌入维数和延迟时间,嵌入维数m=5,延迟时间L=1,子序列长度选取为w=100。所述紫外光谱为TU1810紫外可见分光光度计。
实施例5
一种基于排列熵的水体中磷的检测方法,包括以下步骤:
步骤1:取水体样本直接配制浓度为5mg/L的磷酸二氢钾溶液;
步骤2:通过紫外光谱测量其吸光度;
步骤3:记录下其对应的吸光度;
步骤4:通过排列熵算法对步骤3中得到的数据进行处理:基于排列熵算法的具体过程为将长度为T的光谱信号长时间序列分为若干个长度为w的子序列,选取子序列最大重叠序列,即每个子序列向后移动一个数据点得到下一个子序列,这样得到的子序列个数为T-w+1,采用matlab计算的方法不断测试排列熵的值,通过得到的效果图最佳为依据求出最佳嵌入维数和延迟时间,嵌入维数m=5,延迟时间L=1,子序列长度选取为w=100。所述紫外光谱为TU1810紫外可见分光光度计。
对比例
测定总磷常用分析方法:磷钼蓝分光光度法。
实验及实验结果分析
实验一:磷钼蓝分光光度法测水质中总磷含量
1.1实验所用试剂与仪器为:(1)高压灭菌器或一般民用压力锅;(2)50ml比色管;(3)5%过硫酸钾溶液;(4)磷酸盐标准溶液(用于110℃干燥2h的优级纯磷酸二氢钾配制);(5)10%抗坏血酸溶液;(6)钼酸盐溶液;(7)UV紫外光谱;(8)10mm比色皿。
1.2工作曲线的制作:准确吸取0.00、0.125、0.50、1.00、1.50、2.00、2.50mL磷酸盐标准溶液(100mg/L)于50mL比色管中,用水稀释至25mL,加4ml 5%过硫酸钾溶液,摇匀,于高压灭菌器中密塞消解30min,冷却后,定容至标线,按照试验方法操作,分别向消解后的样品加入1ml的抗坏血酸,30秒后加入2ml的钼酸铵溶液,充分搅匀并静置15min后,使用光程为10mm比色皿,在710nm波长处下,以去离子水做参比,测量其吸光度值。
表1为测量水样中总磷紫外吸收光谱曲线和总磷含量的关系表格的数据
实验二:配制不同浓度的磷酸二氢钾溶液通过光谱仪测量其吸光度。
2.1实验所用试剂与仪器为:(1)高压灭菌器或一般民用压力锅;(2)烧杯;(3)磷酸盐标准溶液(用于110℃干燥2h的优级纯磷酸二氢钾配制);(4)紫外光谱(TU1810紫外可见分光光度计);(5)10mm比色皿。
2.2实验步骤:
(1)配制不同浓度的磷酸二氢钾溶液于烧杯中,分为5个样本:1mg/L、2mg/L、3mg/L、4mg/L、5mg/L。
(2)直接用紫外光谱检测5个样品的吸光度,得到190nm-500nm之间对应的数据。
(3)将5个样品浓度得到的数据带入到matlab中已编写好的排列熵算法的程序中。
实验结果分析,直接配制不同浓度的磷酸二氢钾溶液分别为1mg/L、2mg/L、3mg/L、4mg/L、5mg/L,通过光谱仪测量其吸光度,记录下其不同浓度对应的吸光度,如图1所示,得到吸收峰值在360nm处附近,得出磷酸二氢钾溶液在360nm处有相对明显的吸收峰,记录下其不同浓度对应的吸光度,再通过排列熵算法对此数据进行处理,将嵌入维数m=5,子序列长度w=100,子序列向后移动值L=1带入排列熵算法,得到如图2-6所示。
基于排列熵算法的具体过程为:将长度为T的光谱信号长时间序列分为若干个长度为w的子序列,选取子序列最大重叠序列,即每个子序列向后移动一个数据点得到下一个子序列,这样得到的子序列个数为T-w+1,采用matlab计算的方法不断测试排列熵的值,通过得到的效果图最佳为依据求出最佳嵌入维数和延迟时间,如图2-9所示,将分别将m值带入4、5、8、10,由效果可知,嵌入维数m=5为最佳,如图2-6和图10-12所示,子序列长度选取为w=100为最佳,把光谱波长变化每100个nm看作一个序列,嵌入维数为5,这样会得到一个空间重构向量,对这个重构向量的分量进行一个升序排序,它是m!种排列方式中的一种,再对整个序列中各种排列情况出现的次数进行统计,计算出相对频率,进而求得信息熵即排列熵。熵值最大的点代表的是那个点所在序列变化最显著的吸收光谱的峰值,所以吸收光谱峰值变大,熵值也增大。排列熵算法具有计算简单,抗噪能力强;对时间敏感性强,可获得较高的分辨率;算法的输出结果直观,对突变的信息具有较好的识别性。通过排列熵算法的计算,熵值的变化与光谱吸光度变化趋势类似,但可以明显看出熵值的最高峰即是光谱360nm处最大吸光度,且随浓度的增加,吸光度增大,且熵值也不断增大,并放大了信号的微小变化,减小了光谱检测中相关杂质对水样的影响,这些影响都可以通过熵值反映出来,而吸光度变化极小不易被观测到从而影响测量的准确性,同时通过此方法把很小的吸光度转化为熵值来定量分析磷酸二氢钾的浓度,这样更加直观、准确的反应了水样中总磷的含量,同时不通过其它化学试剂的添加就能得到相对准确的总磷含量的浓度值,与磷钼蓝分光光度法对比,精度与磷钼蓝方法基本保持一致,但无需添加任何化学试剂,减少二次污染,检测操作简单化、自动化,易于实现原位测量,磷钼蓝检测方法的吸光度与浓度的线性关系为0-0.6mg/L,而基于排列熵算法的检测方法,通过熵值的变化可以得到不同浓度的变化,扩大了测量范围,精度与磷钼蓝方法基本保持一致,排列熵算法减小了光谱检测中相关杂质对水样的影响,并放大了信号的微小变化,更适合检测低浓度的溶液。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用。它完全可以被适用于各种适合本发明的领域。对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节。
Claims (4)
1.一种基于排列熵的水体中磷的检测方法,其特征在于,包括以下步骤:
步骤1:取水体样本直接配制浓度为1-5mg/L的磷酸二氢钾溶液;
步骤2:通过紫外光谱测量其吸光度;
步骤3:记录下其不同浓度对应的吸光度;
步骤4:通过排列熵算法对步骤3中得到的数据进行处理:基于排列熵算法的具体过程为将长度为T的光谱信号长时间序列分为若干个长度为w的子序列,选取子序列最大重叠序列,即每个子序列向后移动一个数据点得到下一个子序列,这样得到的子序列个数为T-w+1,采用matlab计算的方法不断测试排列熵的值,通过得到的效果图最佳为依据求出最佳嵌入维数和延迟时间。
2.根据权利要求1所述的一种基于排列熵的水体中磷的检测方法,其特征在于,所述紫外光谱为TU1810紫外可见分光光度计。
3.根据权利要求2所述的一种基于排列熵的水体中磷的检测方法,其特征在于,所述步骤1中配制的磷酸二氢钾溶液选取浓度为1mg/L、2mg/L、3mg/L、4mg/L和5mg/L经步骤2、步骤3,将得到的数据经步骤4中对所得数据进行处理,运用matlab仿真进行排列熵计算。
4.根据权利要求1所述的一种基于排列熵的水体中磷的检测方法,其特征在于,步骤4中得出最佳嵌入维数m=5,延迟时间L=1,子序列长度选取为w=100。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811586350.5A CN109632680B (zh) | 2019-01-29 | 2019-01-29 | 一种基于排列熵的水体中磷的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811586350.5A CN109632680B (zh) | 2019-01-29 | 2019-01-29 | 一种基于排列熵的水体中磷的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109632680A true CN109632680A (zh) | 2019-04-16 |
CN109632680B CN109632680B (zh) | 2021-07-20 |
Family
ID=66077024
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811586350.5A Active CN109632680B (zh) | 2019-01-29 | 2019-01-29 | 一种基于排列熵的水体中磷的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109632680B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110104844A (zh) * | 2019-06-11 | 2019-08-09 | 山西师范大学 | 废水中微量钼元素的新型富集工艺 |
CN113380339A (zh) * | 2021-06-17 | 2021-09-10 | 中国环境科学研究院 | 区域空间污染特征的检测方法 |
CN115069398A (zh) * | 2022-07-27 | 2022-09-20 | 江苏思伽循环科技有限公司 | 一种提高废旧电池回收利用率的电池粉碎控制*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100280762A1 (en) * | 2007-02-14 | 2010-11-04 | Chemimage Corporation | System and Method for Analyzing Biological Samples Using Raman Molecular Imaging |
CN107478599A (zh) * | 2017-10-16 | 2017-12-15 | 北京航空航天大学 | 一种饲料原料中非法添加物三聚氰胺甲醛树脂的检测方法 |
CN108195817A (zh) * | 2016-12-08 | 2018-06-22 | 同方威视技术股份有限公司 | 用于去除溶剂干扰的拉曼光谱检测方法 |
CN109765212A (zh) * | 2019-03-11 | 2019-05-17 | 广西科技大学 | 拉曼光谱中不同步褪色荧光的消除方法 |
-
2019
- 2019-01-29 CN CN201811586350.5A patent/CN109632680B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100280762A1 (en) * | 2007-02-14 | 2010-11-04 | Chemimage Corporation | System and Method for Analyzing Biological Samples Using Raman Molecular Imaging |
CN108195817A (zh) * | 2016-12-08 | 2018-06-22 | 同方威视技术股份有限公司 | 用于去除溶剂干扰的拉曼光谱检测方法 |
CN107478599A (zh) * | 2017-10-16 | 2017-12-15 | 北京航空航天大学 | 一种饲料原料中非法添加物三聚氰胺甲醛树脂的检测方法 |
CN109765212A (zh) * | 2019-03-11 | 2019-05-17 | 广西科技大学 | 拉曼光谱中不同步褪色荧光的消除方法 |
Non-Patent Citations (2)
Title |
---|
XIAOLI LI ET AL: "Application of Permutation Entropy in Feature Extraction for Near-Infrared Spectroscopy Noninvasive Blood Glucose Detection", 《JOURNAL OF SPECTROSCOPY》 * |
王光辉 等: "基于高光谱融合神经网络的玉米黄曲霉毒素B1和赤霉烯酮含量预测", 《食品与机械》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110104844A (zh) * | 2019-06-11 | 2019-08-09 | 山西师范大学 | 废水中微量钼元素的新型富集工艺 |
CN110104844B (zh) * | 2019-06-11 | 2021-12-14 | 山西师范大学 | 废水中微量钼元素的富集工艺 |
CN113380339A (zh) * | 2021-06-17 | 2021-09-10 | 中国环境科学研究院 | 区域空间污染特征的检测方法 |
CN113380339B (zh) * | 2021-06-17 | 2022-02-08 | 中国环境科学研究院 | 区域空间污染特征的检测方法 |
CN115069398A (zh) * | 2022-07-27 | 2022-09-20 | 江苏思伽循环科技有限公司 | 一种提高废旧电池回收利用率的电池粉碎控制*** |
CN115069398B (zh) * | 2022-07-27 | 2024-01-12 | 江苏益佳通新能源科技有限公司 | 一种提高废旧电池回收利用率的电池粉碎控制*** |
Also Published As
Publication number | Publication date |
---|---|
CN109632680B (zh) | 2021-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109632680A (zh) | 一种基于排列熵的水体中磷的检测方法 | |
CN104777156A (zh) | 一种基于碳点荧光关-开模式检测肌醇六磷酸的方法 | |
CN110672574B (zh) | 一种用于检测Cu2+的比率荧光传感器及其制备方法和应用 | |
CN103323115A (zh) | 基于波长调制的气体吸收谱线线宽和线型系数的测量方法 | |
CN108613969A (zh) | 一种白酒中甜蜜素的快速检测方法 | |
CN100498287C (zh) | 检测羟基自由基的纳米银分光光度法 | |
CN103868885A (zh) | 基于复合多次谐波的气体浓度在线测量方法 | |
CN107417671A (zh) | 一种含喹啉取代的香豆素衍生物及其制备方法和在比率型pH荧光探针上的应用 | |
CN105115968B (zh) | 纳米金比色生物传感器快速检测热加工食品中丙烯酰胺的方法 | |
CN106083645A (zh) | 一种铁离子荧光探针化合物及其制备方法和应用 | |
CN104614370A (zh) | 一种基于纳米金的快速检测亚硝酸根的方法 | |
CN109781721B (zh) | 一种测定锌电解液中铜钴镍铁含量的测试体系优化方法 | |
Zamora et al. | Trace determination of carbendazim, fuberidazole and thiabendazole in water by application of multivariate calibration to cross-sections of three-dimensional excitation–emission matrix fluorescence | |
CN108152316B (zh) | 一种以马来酸酐含量和磺化率为指标的糖苷磺基琥珀酸酯盐产品的质量检验方法 | |
CN103134759A (zh) | 一种定量检测6-甲基-2-硫代吡啶-N-乙酰-β-D-氨基葡萄糖苷的检测方法 | |
CN103063790A (zh) | 一种用于检测煤气中硫化氢含量的方法 | |
CN102507482A (zh) | 一种定量检测6-甲基-2-硫代吡啶-N-乙酰-β-D-氨基葡萄糖苷(MPT-NAG)的检测方法及试剂 | |
CN101710075A (zh) | 一种测定铝酸钠溶液中微量元素的方法 | |
CN110927091A (zh) | 一种哌嗪水溶液定量检测的方法 | |
CN113340862B (zh) | 荧光分子传感器及制法、水中痕量铀酰离子的检测方法 | |
CN108303388A (zh) | 一种原位定量表征复杂有机物和金属离子络合过程的方法 | |
CN104406947B (zh) | 利用槲皮素作为荧光探针在测定f‑中的应用及其测定方法 | |
CN104569015A (zh) | 测量钢铁中的铅的能量色散x射线荧光光谱法 | |
CN104634773A (zh) | 一种用等离子发射光谱测定三苯基磷中铑含量的方法 | |
CN113758908A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |