CN106769996A - 一种小麦粉中山梨酸钾检测方法 - Google Patents
一种小麦粉中山梨酸钾检测方法 Download PDFInfo
- Publication number
- CN106769996A CN106769996A CN201710087815.1A CN201710087815A CN106769996A CN 106769996 A CN106769996 A CN 106769996A CN 201710087815 A CN201710087815 A CN 201710087815A CN 106769996 A CN106769996 A CN 106769996A
- Authority
- CN
- China
- Prior art keywords
- sample
- potassium sorbate
- concentration
- matrix
- absorption peak
- 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
- CHHHXKFHOYLYRE-UHFFFAOYSA-M 2,4-Hexadienoic acid, potassium salt (1:1), (2E,4E)- Chemical compound [K+].CC=CC=CC([O-])=O CHHHXKFHOYLYRE-UHFFFAOYSA-M 0.000 title claims abstract description 109
- 239000004302 potassium sorbate Substances 0.000 title claims abstract description 109
- 235000010241 potassium sorbate Nutrition 0.000 title claims abstract description 109
- 229940069338 potassium sorbate Drugs 0.000 title claims abstract description 109
- 241000209140 Triticum Species 0.000 title claims abstract description 39
- 235000021307 Triticum Nutrition 0.000 title claims abstract description 39
- 238000001514 detection method Methods 0.000 title claims abstract description 31
- 235000013312 flour Nutrition 0.000 title claims abstract description 22
- 238000010521 absorption reaction Methods 0.000 claims abstract description 107
- 239000000203 mixture Substances 0.000 claims abstract description 41
- 238000000034 method Methods 0.000 claims abstract description 19
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 14
- 238000001328 terahertz time-domain spectroscopy Methods 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 69
- 238000001228 spectrum Methods 0.000 claims description 25
- 238000012417 linear regression Methods 0.000 claims description 17
- 238000004611 spectroscopical analysis Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 5
- 235000013305 food Nutrition 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 4
- 238000003908 quality control method Methods 0.000 abstract description 2
- 238000010238 partial least squares regression Methods 0.000 description 17
- 239000000843 powder Substances 0.000 description 14
- 239000004698 Polyethylene Substances 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 8
- -1 polyethylene Polymers 0.000 description 8
- 229920000573 polyethylene Polymers 0.000 description 8
- 238000003384 imaging method Methods 0.000 description 7
- 238000004445 quantitative analysis Methods 0.000 description 7
- 230000003595 spectral effect Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 239000000654 additive Substances 0.000 description 4
- 230000000996 additive effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 108090000623 proteins and genes Proteins 0.000 description 4
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 229910052700 potassium Inorganic materials 0.000 description 3
- 239000011591 potassium Substances 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 238000004451 qualitative analysis Methods 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 230000002421 anti-septic effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 235000013339 cereals Nutrition 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 150000001875 compounds Chemical group 0.000 description 2
- BEFDCLMNVWHSGT-UHFFFAOYSA-N ethenylcyclopentane Chemical compound C=CC1CCCC1 BEFDCLMNVWHSGT-UHFFFAOYSA-N 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000004334 sorbic acid Substances 0.000 description 2
- 235000010199 sorbic acid Nutrition 0.000 description 2
- 229940075582 sorbic acid Drugs 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- WCUXLLCKKVVCTQ-UHFFFAOYSA-M Potassium chloride Chemical compound [Cl-].[K+] WCUXLLCKKVVCTQ-UHFFFAOYSA-M 0.000 description 1
- 238000005411 Van der Waals force Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005251 capillar electrophoresis Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 235000021472 generally recognized as safe Nutrition 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- VIKNJXKGJWUCNN-XGXHKTLJSA-N norethisterone Chemical compound O=C1CC[C@@H]2[C@H]3CC[C@](C)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1 VIKNJXKGJWUCNN-XGXHKTLJSA-N 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 235000011164 potassium chloride Nutrition 0.000 description 1
- 239000001103 potassium chloride Substances 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000004809 thin layer chromatography Methods 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/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
- 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/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (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
技术领域
本发明涉及农产品/食品品质控制与质量安全技术领域，更具体的涉及一种小麦粉中山梨酸钾检测方法。
背景技术
由于小麦是我国主要的粮食作物之一，作为一种公认安全、高效防腐的食品添加剂，山梨酸及钾盐在食品/农产品行业的应用越来越广泛。但如果山梨酸钾含量超标，造成人体中含量过多，在一定程度上对肾、肝脏的健康产生不利影响。
目前，在农产品/食品行业常用的山梨酸钾检测方法主要有：薄层层析、高效毛细管电泳法、高效益相色谱法等。然而，这些技术在实际应用中存在前期处理复杂，操作专业性强，检测时间长问题。
综上所述，现有技术中的山梨酸钾检测方法，存在检测过程复杂困难和耗时长的问题。
发明内容
本发明实施例提供一种小麦粉中山梨酸钾检测方法，用以解决现有技术中存在检测过程复杂困难耗和时长的问题。
本发明实施例提供一种小麦粉中山梨酸钾检测方法，包括：
获取山梨酸钾与小麦粉混合物样品；
采用太赫兹时域光谱系统，确定所述样品的吸收光谱；
根据所述吸收光谱，确定所述样品的特征吸收峰的吸收系数；
根据所述样品的特征吸收峰的吸收系数，通过所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系，确定所述样品中山梨酸钾的浓度值。
较佳地，所述采用太赫兹时域光谱系统，确定所述样品的吸收光谱图；包括：
采用太赫兹时域光谱系统，确定所述样品的时域光谱；
根据所述样品的时域光谱，确定所述样品的频域光谱；
根据所述样品的频域光谱，确定所述样品的吸收光谱。
较佳地，采用简单一元线性回归法，确定所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系；所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系，由下式确定：
Y＝a0+a1X+ε
其中，Y为样品中山梨酸钾的浓度矩阵；X为样品的特征吸收峰的吸收系数矩阵；a0和a1为回归系数，ε为残差。
较佳地，采用偏最小二乘回归法，确定所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系；所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系，由下式确定：
Y＝TBQ
其中，T为样品的特征吸收峰的吸收系数矩阵的得分矩阵；B为样品的特征吸收峰的吸收系数矩阵的得分矩阵和样品中山梨酸钾的浓度矩阵的得分矩阵作线性回归处理的回归因子；Q为样品中山梨酸钾的浓度矩阵的主成分矩阵。
本发明实施例中，提供一种小麦粉中山梨酸钾检测方法，该方法采用太赫兹时域光谱系统完成小麦粉中山梨酸钾的检测，由于太赫兹光谱具有反映化合物结构的指纹特征的特点，许多生物分子之间弱的相互作用(氢键、范德华力等)、生物大分子的骨架振动等低频振动吸收正好处于THz频带范围，并且THz光谱技术对探测物质结构存在的微小差异和变化非常灵敏，使得其具有太赫兹指纹谱，太赫兹成像技术不但能够分辨出物体的形貌，还能够鉴别物体的组成成分，基于THz辐射的以上特点，利用其对山梨酸钾及其混合物进行无损检测，结合定性定量分析算法，可以将山梨酸钾从混合物样品中分离出来。该方法解决了山梨酸钾的光谱表征与种类区分问题，无需添加其他化学试剂，其操作过程简单，数据处理快速，检测准确性高，效果良好；并且THz电磁辐射的光子能量很低，只有毫电子伏特，不会因为电离而破坏生物分子，太赫兹环境下对人体不会带来辐射上的危害，因而是一种安全有效的无损检测方法。
附图说明
图1为本发明实施例提供的一种小麦粉中山梨酸钾检测方法流程示意图；
图2为本发明实施例提供的山梨酸钾折射谱图；
图3为本发明实施例提供的山梨酸钾吸收谱图；
图4为本发明实施例提供的山梨酸钾与聚乙烯混合物的简单一元线性回归实际浓度与预测浓度相关图；
图5为本发明实施例提供的山梨酸钾与小麦粉混合物的简单一元线性回归实际浓度与预测浓度相关图；
图6为本发明实施例提供的山梨酸钾和聚乙烯混合物与山梨酸钾和小麦粉混合物的偏最小二乘回归实际浓度与预测浓度的相关图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1示例性的示出了本发明实施例提供的一种小麦粉中山梨酸钾检测方法流程示意图。如图1所示，本发明实施例提供的一种小麦粉中山梨酸钾检测方法，包括以下步骤：
步骤101，获取山梨酸钾与小麦粉混合物样品。
需要说明的是，选取小麦粉和山梨酸钾，然后将其研磨，实验前需要将样品粉末进行压片。
步骤102，采用太赫兹时域光谱系统，确定样品的吸收光谱。
较佳地，采用太赫兹时域光谱系统，确定样品的吸收光谱图；具体包括：采用太赫兹时域光谱系统，确定样品的时域光谱(时间与幅值光谱图)；根据样品的时域光谱，确定样品的频域光谱(太赫兹THz与幅值光谱图)；根据样品的频域光谱，确定样品的吸收光谱(太赫兹THz与吸收度光谱图)。
需要说明的是，太赫兹波由于其能量低、穿透性强及指纹谱特性等独特优势，使得太赫兹成像技术是毫米波和X射线成像技术的重要补充。根据样品测量方式的不同，太赫兹成像可以分为透射成像和反射成像两种。在太赫兹反射成像中，不仅能够获得样品表面反射的THz脉冲，并且能够获得穿透样品内部反射的THz脉冲，因此，可以获得样品内部的有用的物理和化学信息。
需要说明的是，本发明采用一种利用太赫兹光谱技术表征了添加剂山梨酸钾，获得其宽频谱段内的特征吸收峰，利用在宽频段内的太赫兹指纹峰特征来区别不同种类的添加剂。接着对山梨酸钾和聚乙烯混合物，以及山梨酸钾和小麦粉混合物的THz光谱进行定性与定量分析，利用简单一元线性发SLR和偏最小二乘回归法PLS建立了吸收系数和浓度之间的回归模型，最后对混合物中山梨酸钾含量进行计算。由此可知，利用THz技术进行山梨酸钾定性与定量分析的可行性，具有广阔的应用前景。
需要说明的是，样品通过太赫兹时域光谱系统后得到样品的时域谱、频域谱，通过简单的计算，还可以得到样品的折射率谱。
步骤103，根据吸收光谱，确定样品的特征吸收峰的吸收系数。
需要说明的是，由于吸收光谱图的横坐标为频率，纵坐标为吸收度(即吸收系数)，分析可知，特征吸收峰的频率与添加剂的种类有关(0.98THz处是山梨酸钾添加剂)，特征吸收峰的吸收系数与添加剂的浓度有关。
步骤104，根据样品的特征吸收峰的吸收系数，通过样品的特征吸收峰的吸收系数与样品中山梨酸钾的浓度之间的关系，确定所述样品中山梨酸钾的浓度值。
需要说明的是，样品的特征吸收峰的吸收系数与样品中山梨酸钾的浓度之间的关系，具体包括定性关系、采用简单一元线性回归法的定量关系和采用偏最小二乘回归法的定量关系。
1)定性分析：
通过折射谱图和吸收谱图均可获知，山梨酸钾的浓度随着山梨酸钾特征吸收峰的吸收系数的增加而增大。
2)采用简单一元线性回归法的定量分析：
采用简单一元线性回归法，确定样品的特征吸收峰的吸收系数与样品中山梨酸钾的浓度之间的关系；样品的特征吸收峰的吸收系数与样品中山梨酸钾的浓度之间的关系，由下式确定：
Y＝a0+a1X+ε
其中，Y为样品中山梨酸钾的浓度矩阵；X为样品的特征吸收峰的吸收系数矩阵；a0和a1为回归系数，ε为残差。
需要说明的是，残差ε服从N(0,σ2)分布，使用最小二乘法使残差的平方最小，需要满足n为样品数量。
需要说明的是，简单一元线性回归法是直接找出输出变量(一维特征吸收峰的吸收系数值)与输入变量(浓度值)的线性相关性。通过简单一元线性回归法获得样品的特征吸收峰的吸收系数与样品中山梨酸钾的浓度之间的关系，其操作简单，速度快，适用于对浓度检测要求不高的场合。
3)采用偏最小二乘回归法的定量分析：
采用偏最小二乘回归法，确定样品的特征吸收峰的吸收系数与样品中山梨酸钾的浓度之间的关系；样品的特征吸收峰的吸收系数与样品中山梨酸钾的浓度之间的关系，由下式确定：
Y＝TBQ
其中，T为样品的特征吸收峰的吸收系数矩阵的得分矩阵；B为样品的特征吸收峰的吸收系数矩阵的得分矩阵和样品中山梨酸钾的浓度矩阵的得分矩阵作线性回归处理的回归因子；Q为样品中山梨酸钾的浓度矩阵的主成分矩阵。
需要说明的是，PLS采用逐步提取的方法来选取主成分，此次类推，直到满足精度要求为止。此方法提取的主成分，不仅要求其方差尽可能的大，还要求与浓度的相关性尽可能的大，主成分也是按照对原光谱信息概括能力从高到低依次排列的。PLS主要用于矩阵X和Y之间的线性分解，建模步骤如下：
a、对矩阵X和矩阵Y进行分解，其公式如下：
其中，X为样品的特征吸收峰的吸收系数矩阵；Y为样品中山梨酸钾的浓度矩阵；T为样品的特征吸收峰的吸收系数矩阵的得分矩阵；U为样品中山梨酸钾的浓度矩阵的得分矩阵，P为样品的特征吸收峰的吸收系数矩阵的主成分矩阵；Q为样品中山梨酸钾的浓度矩阵的主成分矩阵，E为样品的特征吸收峰的吸收系数矩阵的残差矩阵；F为样品中山梨酸钾的浓度矩阵的残差矩阵；tk为样品的特征吸收峰的吸收系数矩阵的第k个主因子的得分；pk为样品的特征吸收峰的吸收系数矩阵的第k个主因子的载荷；uk为样品中山梨酸钾的浓度矩阵的第k个主因子的得分；qk为样品中山梨酸钾的浓度矩阵的第k个主因子的载荷；f为主因子数。
b、对得分矩阵T和U作线性回归处理，其公式如下：
U＝TB
其中，回归因子B可以由上式求得：
B＝(TT)-1TTU
c、根据样品的特征吸收峰的吸收系数矩阵的主成分矩阵求出样品的特征吸收峰的吸收系数矩阵的得分矩阵，然后得出预测样品浓度：
Y＝TBQ
需要说明的是，PLS的本质可归结为分别建立样品的特征吸收峰的吸收系数矩阵对其得分矩阵和样品中山梨酸钾的浓度矩阵对其得分矩阵的一元线性回归模型。通过对样品的特征吸收峰的吸收系数矩阵和样品中山梨酸钾的浓度矩阵的主成分分解，利用分解后的得分矩阵来建立线性回归方程，确定系数矩阵。
需要说明的是，PLS模型在建立过程中，在保证样品的特征吸收峰的吸收系数矩阵的得分矩阵和样品中山梨酸钾的浓度矩阵有最大相关性的同时，使样品的特征吸收峰的吸收系数矩阵的得分矩阵的方差尽可能大，这样使尽可能多的有用信息被提取出来，使其与样品中山梨酸钾的浓度矩阵产生最大线性关系。
需要说明的是，PLS回归模型预测的混合物中山梨酸钾浓度的散点对于零误差线更加密集，表明该回归模型对混合物定量分析结果要优于SLR回归模型分析的结果。其中，在SLR和PLS模型中，山梨酸钾与聚乙烯混合物的散点相对于零误差线变化不明显，而山梨酸钾与小麦混合物的散点变化较明显。因此，利用THz波技术并结合化学计量学方法进行储粮中防腐剂含量检测是可行的，为防腐剂含量的检测提供了一种新的快速无损检测方法。
本发明提供的一种小麦粉中山梨酸钾检测方法的分析实例，如下所示：
1)样品准备
选取小麦粉和山梨酸钾，然后将其研磨，混合物样品还需按照一定比例进行混合均匀后保存并标识，实验前需要将所有样品粉末进行压片。
山梨酸钾含量分别为2％、3％、7％、10％、30％、40％、50％等7种浓度的山梨酸钾与聚乙烯混合物样品和纯聚乙烯(山梨酸钾含量为0)样品，以及山梨酸钾含量分别为2％、3％、7％、10％、30％、40％、50％等7种浓度的山梨酸钾与小麦粉末混合物样品和纯小麦粉末(山梨酸钾含量为0)、纯山梨酸钾(含量为100％)样品。
2)样品的太赫兹光谱测量
采用太赫兹时域光谱系统对山梨酸钾、不同浓度混合物样品进行测量，获得它们的THz时域光谱，即获得参考信号和不同浓度山梨酸钾与小麦粉混合物经过太赫兹时域光谱系统后的时域光谱，从时域光谱分析可知，样品的THz光谱存在明显的区别(振幅和时间延迟)，说明它们的折射率和吸收系数不同。
对参考信号和不同浓度山梨酸钾与小麦粉混合物经过太赫兹时域光谱系统后的时域光谱进行频域处理，获得参考信号和不同浓度山梨酸钾与小麦粉混合物经过太赫兹时域光谱系统后的频域光谱；从频域光谱分析可知，与参考信号(未放置样品)相比，不同小麦混合物信号的相位和幅度均存在一定程度的延迟和衰减；不同小麦混合物信号的延迟时间还存在差异。表明混合物样品的折射率不同，且大于空气在THz波段的折射率；不同小麦混合物的频谱曲线表明其吸收系数和折射率有所不同，体现出不同样品的物理和化学特性；样品对高频部分的太赫兹吸收较多，导致频谱变窄。
3)定性分析
图2为本发明实施例提供的山梨酸钾折射谱图；由图2可以看出，山梨酸钾在0.2～2THz的范围内折射率波动幅度在1.96～2.15之间，折射率平均为2左右。波动比较小，在0.98THz处有一个小的幅度跳动，这和吸收谱中吸收峰所处的位置一致。
图3为本发明实施例提供的山梨酸钾吸收谱图；由图3可以看出，在0.2～1THz频段内时，吸收系数谱存在着明显的线性关系。吸收系数随着频率增加而增大。在0.98THz处存在明显的强吸收峰，说明样品在该频率的对太赫兹波吸收特别显著。当频率在1.6THz以后，随着频率增加，吸收谱曲线呈缓慢上升趋势。
通过分析不同浓度山梨酸钾与小麦粉混合物吸收谱可知：吸收系数随着频率的增加而逐渐增大，在0.98THz处存在明显的吸收峰，处当频率在1.6THz以后，吸收系数谱的质量开始明显下降，在1.8THZ之后，频谱出现巨大波动，山梨酸钾样品的色散现象反常，不同浓度的混合物的吸收峰强度值存在一定的影响，从1％到50％增加时，吸收谱增多，且强度明显变化。
4)定量分析
为了检测山梨酸钾和小麦粉混合物的吸收系数和混合物中山梨酸钾的质量分数的关系，本发明采用混合物特征吸收峰作为样本数据，结合简单一元线性回归(SimpleLinearRegression,SLR)方法对样品进行定量分析，可以得到混合物中山梨酸钾的含量即浓度。
图4为本发明实施例提供的山梨酸钾与聚乙烯混合物的简单一元线性回归实际浓度与预测浓度相关图；图5为本发明实施例提供的山梨酸钾与小麦粉混合物简单一元线性回归的实际浓度与预测浓度相关图；由图4和图5可知，山梨酸钾浓度的预测整体效果较好，但也有预测结果存在较大误差，主要由于实验过程中混合物样品数量不大，并且样品制备过程中难以保证混合物混合均匀。
上述采用混合物峰值与浓度存在线性关系进行拟合预测，为了与上述方法进行对照，提高模型的预测性能，采用偏最小二乘回归法(PLS)对山梨酸钾混合物有效光谱频段(0.2-1.6THz)信息来建立回归模型，将预测集中的样品代入PLS回归模型中对混合物中山梨酸钾进行计算。
图6为本发明实施例提供的山梨酸钾和聚乙烯混合物与山梨酸钾和小麦粉混合物的偏最小二乘回归实际浓度与预测浓度的相关图；从图6可以看出，直线表示零误差，即实际浓度与预测浓度一致。
以上公开的仅为本发明的几个具体实施例，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
Claims (4)
1.一种小麦粉中山梨酸钾检测方法，其特征在于，包括：
获取山梨酸钾与小麦粉混合物样品；
采用太赫兹时域光谱系统，确定所述样品的吸收光谱；
根据所述吸收光谱，确定所述样品的特征吸收峰的吸收系数；
根据所述样品的特征吸收峰的吸收系数，通过所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系，确定所述样品中山梨酸钾的浓度值。
2.如权利要求1所述的小麦粉中山梨酸钾检测方法，其特征在于，所述采用太赫兹时域光谱系统，确定所述样品的吸收光谱图；包括：
采用太赫兹时域光谱系统，确定所述样品的时域光谱；
根据所述样品的时域光谱，确定所述样品的频域光谱；
根据所述样品的频域光谱，确定所述样品的吸收光谱。
3.如权利要求1所述的小麦粉中山梨酸钾检测方法，其特征在于，采用简单一元线性回归法，确定所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系；所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系，由下式确定：
Y＝a0+a1X+ε (1)
其中，Y为样品中山梨酸钾的浓度矩阵；X为样品的特征吸收峰的吸收系数矩阵；a0和a1为回归系数，ε为残差。
4.如权利要求1所述的小麦粉中山梨酸钾检测方法，其特征在于，采用偏最小二乘回归法，确定所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系；所述样品的特征吸收峰的吸收系数与所述样品中山梨酸钾的浓度之间的关系，由下式确定：
Y＝TBQ (2)
其中，T为样品的特征吸收峰的吸收系数矩阵的得分矩阵；B为样品的特征吸收峰的吸收系数矩阵的得分矩阵和样品中山梨酸钾的浓度矩阵的得分矩阵作线性回归处理的回归因子；Q为样品中山梨酸钾的浓度矩阵的主成分矩阵。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111271400.2A CN114002180A (zh) | 2017-02-18 | 2017-02-18 | 一种获得小麦粉中山梨酸钾检测曲线的方法和应用 |
| CN201710087815.1A CN106769996A (zh) | 2017-02-18 | 2017-02-18 | 一种小麦粉中山梨酸钾检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710087815.1A CN106769996A (zh) | 2017-02-18 | 2017-02-18 | 一种小麦粉中山梨酸钾检测方法 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111271400.2A Division CN114002180A (zh) | 2017-02-18 | 2017-02-18 | 一种获得小麦粉中山梨酸钾检测曲线的方法和应用 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106769996A true CN106769996A (zh) | 2017-05-31 |
Family
ID=58958100
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710087815.1A Pending CN106769996A (zh) | 2017-02-18 | 2017-02-18 | 一种小麦粉中山梨酸钾检测方法 |
| CN202111271400.2A Pending CN114002180A (zh) | 2017-02-18 | 2017-02-18 | 一种获得小麦粉中山梨酸钾检测曲线的方法和应用 |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111271400.2A Pending CN114002180A (zh) | 2017-02-18 | 2017-02-18 | 一种获得小麦粉中山梨酸钾检测曲线的方法和应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN106769996A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110553983A (zh) * | 2019-09-05 | 2019-12-10 | 东北农业大学 | 一种牛奶中山梨酸钾的快速检测方法 |
| CN118518620A (zh) * | 2024-07-25 | 2024-08-20 | 中国工程物理研究院激光聚变研究中心 | 检测小麦粉中甜蜜素含量的方法、设备、存储介质及产品 |
-
2017
- 2017-02-18 CN CN201710087815.1A patent/CN106769996A/zh active Pending
- 2017-02-18 CN CN202111271400.2A patent/CN114002180A/zh active Pending
Non-Patent Citations (6)
| Title |
|---|
| 史永刚等: "《化学计量学方法及MATLAB实现》", 30 April 2010, 中国石化出版社 * |
| 汪景荣等: "偏最小二乘法和THzTDS在正品大黄鉴别中的应用", 《光谱学与光谱分析》 * |
| 王彩云等: "近红外光谱法快速检测牛奶中氯霉素残留", 《食品科学》 * |
| 罗洪群等: "《市场调查与预测 第2版》", 31 January 2016, 清华大学出版社 * |
| 葛宏义: "储粮品质的THz波检测理论与分析方法研究", 《中国科学院大学博士学位论文》 * |
| 陆婉珍等: "《现代近红外光谱分析技术》", 30 April 2000, 中国石化出版社 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110553983A (zh) * | 2019-09-05 | 2019-12-10 | 东北农业大学 | 一种牛奶中山梨酸钾的快速检测方法 |
| CN110553983B (zh) * | 2019-09-05 | 2022-03-04 | 东北农业大学 | 一种牛奶中山梨酸钾的快速检测方法 |
| CN118518620A (zh) * | 2024-07-25 | 2024-08-20 | 中国工程物理研究院激光聚变研究中心 | 检测小麦粉中甜蜜素含量的方法、设备、存储介质及产品 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114002180A (zh) | 2022-02-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103335976B (zh) | 利用太赫兹时域光谱检测硫酸盐、硝酸盐溶液浓度的方法 | |
| CN102192890A (zh) | 利用近红外光谱分析技术快速测定木材化学成分 | |
| CN104181122A (zh) | 利用太赫兹时域光谱技术检测粮食中农药含量的方法 | |
| CN103969214A (zh) | 利用太赫兹频段红外光谱技术检测粮食中农药含量的方法 | |
| CN106525761A (zh) | 基于太赫兹光谱扫描的亚硝酸盐检测方法 | |
| CN103543132B (zh) | 一种基于小波变换的煤质特性测量方法 | |
| Liu et al. | Prediction of soil organic carbon with different parent materials development using visible-near infrared spectroscopy | |
| CN103969212A (zh) | 利用太赫兹频段ftir技术定量检测粮食中农药残留的方法 | |
| CN104215599A (zh) | 一种利用THz-TDS技术定量检测粮食中农药残留的方法 | |
| CN104406923A (zh) | 一种利用THz-TDS技术定量检测粮食中氨基酸含量的方法 | |
| CN108020526A (zh) | 一种丁羟推进剂药浆组分近红外检测方法 | |
| CN106769996A (zh) | 一种小麦粉中山梨酸钾检测方法 | |
| CN104297202B (zh) | 利用THz‑TDS频域谱定量检测粮食中农药残留的方法 | |
| CN102042967A (zh) | 一种基于近红外光谱技术的葡萄糖水溶液快速识别方法 | |
| CN104020124B (zh) | 基于吸收率择优的分光波长筛选方法 | |
| Jiang et al. | Determination of potassium sorbate and sorbic acid in agricultural products using THz time-domain spectroscopy | |
| CN106706554A (zh) | 一种快速无损测定玉米单穗籽粒直链淀粉含量的方法 | |
| CN101303295A (zh) | 红外光谱定量分析混合溶液组成的快速分析方法 | |
| CN102262055B (zh) | 一种测定聚丙烯酰胺类物质中丙烯酰胺单体残留量的方法 | |
| JP2019158560A5 (zh) | ||
| CN107328733A (zh) | 一种快速检测鱼糜中添加的淀粉含量的方法 | |
| CN106872398A (zh) | 一种hmx炸药水分含量快速测量方法 | |
| CN109580443A (zh) | 采用太赫兹技术检测物质中金属颗粒含量的方法 | |
| CN115436315A (zh) | 一种基于近红外光谱的水泥外加剂检测方法 | |
| CN110455742B (zh) | 利用太赫兹时域光谱技术鉴别黄曲霉毒素b1和b2的方法 |
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 |