CN106769992A - 基于差分进化的谷氨酰胺太赫兹谱波长选择方法及装置 - Google Patents
基于差分进化的谷氨酰胺太赫兹谱波长选择方法及装置 Download PDFInfo
- Publication number
- CN106769992A CN106769992A CN201611103971.4A CN201611103971A CN106769992A CN 106769992 A CN106769992 A CN 106769992A CN 201611103971 A CN201611103971 A CN 201611103971A CN 106769992 A CN106769992 A CN 106769992A
- Authority
- CN
- China
- Prior art keywords
- glutamine
- population
- terahertz
- initial population
- species
- 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
- ZDXPYRJPNDTMRX-UHFFFAOYSA-N glutamine Natural products OC(=O)C(N)CCC(N)=O ZDXPYRJPNDTMRX-UHFFFAOYSA-N 0.000 title claims abstract description 81
- 238000001228 spectrum Methods 0.000 title claims abstract description 44
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 45
- 238000004445 quantitative analysis Methods 0.000 claims abstract description 16
- 230000035772 mutation Effects 0.000 claims abstract description 14
- ZDXPYRJPNDTMRX-VKHMYHEASA-N L-glutamine Chemical compound OC(=O)[C@@H](N)CCC(N)=O ZDXPYRJPNDTMRX-VKHMYHEASA-N 0.000 claims description 82
- 238000000034 method Methods 0.000 claims description 11
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 claims description 8
- 239000000203 mixture Substances 0.000 claims description 7
- 230000000717 retained effect Effects 0.000 claims description 5
- 150000001408 amides Chemical class 0.000 claims description 4
- 229910021529 ammonia Inorganic materials 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 13
- 238000002474 experimental method Methods 0.000 description 4
- 241001269238 Data Species 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- WHUUTDBJXJRKMK-UHFFFAOYSA-N Glutamic acid Natural products OC(=O)C(N)CCC(O)=O WHUUTDBJXJRKMK-UHFFFAOYSA-N 0.000 description 1
- WHUUTDBJXJRKMK-VKHMYHEASA-N L-glutamic acid Chemical compound OC(=O)[C@@H](N)CCC(O)=O WHUUTDBJXJRKMK-VKHMYHEASA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 235000013922 glutamic acid Nutrition 0.000 description 1
- 239000004220 glutamic acid Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- WJCNZQLZVWNLKY-UHFFFAOYSA-N thiabendazole Chemical compound S1C=NC(C=2NC3=CC=CC=C3N=2)=C1 WJCNZQLZVWNLKY-UHFFFAOYSA-N 0.000 description 1
- 229960004546 thiabendazole Drugs 0.000 description 1
- 235000010296 thiabendazole Nutrition 0.000 description 1
- 239000004308 thiabendazole Substances 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Genetics & Genomics (AREA)
- Toxicology (AREA)
- Physiology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种基于差分进化的谷氨酰胺太赫兹谱波长选择方法及装置。本发明首先对初始种群进行基于差分的变异操作和交叉操作,以得到交叉种群;然后分别利用初始种群和交叉种群从谷氨酰胺样品的太赫兹吸收谱中进行选取,利用所构造的适应度函数分别计算初始种群和交叉种群中每个个体的适应度;并将适应度值较大的个体保留下来,从而得到新一代种群;最后将新一代种群作为新的初始种群进行进化迭代,直至进化代数达到设定阈值,并将最终代种群中适应度值最高的个体作为所选择的谷氨酰胺太赫兹吸收谱波长的最优解。本发明通过对谷氨酰胺样品的太赫兹吸收谱进行逐点深度选择,挑选有用信息,从而显著提高了定量分析的精度,取得了良好的应用效果。
Description
技术领域
本发明涉及一种基于差分进化的谷氨酰胺太赫兹谱波长选择方法及装置,属于太赫兹光谱检测技术领域。
背景技术
在对谷氨酰胺样品进行太赫兹吸收谱定量分析中,通过实验得到的谷氨酰胺样品的原始太赫兹吸收谱通常涵盖一段较宽的频段,包含大量的波长点数据,其中不仅包括信噪比较高的有用数据,也包含信噪比较低的噪声数据以及不属于任一组分特征的冗余数据,若直接将原始吸收谱用于定量分析势必导致较高误差,因此需要进行适当选择。由于吸收谱是由一系列波长点数据组成的,对吸收谱数据的选择实际上就是对波长的选择,因而在光谱学中被定义为波长选择(Wavelength selection)。对于太赫兹光谱定量分析领域而言,波长选择对定量分析的准确度至关重要,若选择不恰当,会导致较大误差。但是目前在太赫兹光谱定量分析中,波长选择常用的做法是人为地依据经验从原始光谱中选取某一波段数据用于定量计算,而对太赫兹光谱波长选择的机理及方法缺乏系统性的深入研究。
中国计量学院的王强教授等人分别利用偏最小二乘法(partial least squares,PLS)、区间偏最小二乘法(interval PLS,iPLS)、向后区间偏最小二乘法(backward iPLS,biPLS)以及移动窗口偏最小二乘法(moving window PLS,mwPLS)对噻苯咪唑位于0.3-1.6THz频段内的太赫兹特征光谱进行了波长选择,并对四种算法的性能进行了细致的比较。桂林电子科技大学的陈涛等人就太赫兹光谱定量分析中的特征谱区筛选进行了相关研究。除上述王强等人提出的波长选择方法外,又采用了联合区间偏最小二乘法(siPLS)并进行了一系列对比。但是基于偏最小二乘的波长选择方法,是通过将原始光谱分割成若干区间加以筛选,难免会将部分无意义数据含入其中,甚至将一些有意义数据错误地抛弃,以至于所选择的波长不合适,从而导致定量分析的误差比较大。
发明内容
本发明的目的是提供一种基于差分进化的谷氨酰胺太赫兹谱波长选择方法,以解决目前谷氨酰胺太赫兹吸收谱波长选择不合适而导致定量分析的误差比较大的问题。同时本发明还提供了一种基于差分进化的谷氨酰胺太赫兹谱波长选择装置。
本发明为解决上述技术问题而提供一种基于差分进化的谷氨酰胺太赫兹谱波长选择方法,该选择方法包括以下步骤:
1)随机生成一个大小为S的初始种群X;
2)对初始种群X进行基于差分的变异操作,以得到变异种群V,并对变异种群V进行交叉操作,以得到交叉种群U;
3)分别利用初始种群X和交叉种群U从谷氨酰胺样品的太赫兹吸收谱中进行选取,以得到初始种群X和交叉种群U中每个个体相对应的经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱;
4)构建适应度函数,利用所构造的适应度函数分别计算初始种群X和交叉种群U中每个个体的适应度;
5)对初始种群X和交叉种群U中相对应个体的适应度值进行比较,将适应度值较大的个体保留下来,从而得到新一代种群Xnext;
6)将新一代种群Xnext作为新的初始种群,重复步骤2)-5),直至进化代数达到设定阈值,并将最终代种群中适应度值最高的个体作为所选择的谷氨酰胺太赫兹吸收谱波长的最优解。
进一步地,所述步骤2)中的变异操作如下::
其中Vi,j表示变异种群V中第i个个体的第j个二进制元素,Xi,a、Xi,b和Xi,c分别表示初始种群X中第i个个体的第a,b和c个二进制元素,并且a,b和c是3个互不相同的正整数。
进一步地,所述步骤2)中的交叉操作为:
其中Ui,j表示交叉种群U中第i个个体的第j个二进制元素,rand为0至1区间内的随机数,CR为预先设定的交叉概率,jrand为1至L区间内的随机整数,Xi,j表示初始种群X中第i个个体的第j个二进制元素。
进一步地,所述步骤1)中的初始种群X由S个长度为L的二进制字符串组成,该二进制字符串与谷氨酰胺样品的太赫兹吸收谱中的L个频率点一一对应。
进一步地,所述步骤3)中的波长选择过程如下:
对于初始种群X中的各个个体,若其某个二进制元素为“1”,则对应谷氨酰胺太赫兹吸收谱的频率点数据被保留,否则该频率点数据则被抛弃,将所有保留下的频率点数据整合在一起,组成经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱;对于交叉种群U中的各个个体,若其某个二进制元素为“1”,则对应谷氨酰胺太赫兹吸收谱的频率点数据被保留,否则该频率点数据则被抛弃,将所有保留下的频率点数据整合在一起,组成经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱。
进一步地,所述步骤4)中构建的适应度函数为:
其中F是适应度值,m是校正集中谷氨酰胺样品的总数量,qe是每个谷氨酰胺样品对应的定量分析误差,n代表校正集中谷氨酰胺样品的某一个,其中ccal和creal分别是谷氨酰胺样品的计算浓度和真实浓度。
进一步地,所述新一代种群Xnext满足下式
其中Xi next为新一代种群Xnext的第i个个体,F(Ui)为交叉种群U中第i个个体的适应度值,F(Xi)为初始种群X中第i个个体的适应度值。
本发明还提供了一种基于差分进化的谷氨酰胺太赫兹谱波长选择装置,该选择装置包括初始种群生成模块、变异交叉操作模块、太赫兹吸收谱重构模块、适应度计算模块、比较模块和选择模块,
所述初始种群生成模块用于随机生成一个大小为S的初始种群X;
所述的变异交叉操作模块用于对初始种群X进行基于差分的变异操作,以得到变异种群V,并对变异种群V进行交叉操作,以得到交叉种群U;
所述的太赫兹吸收谱重构模块用于分别利用初始种群X和交叉种群U从谷氨酰胺样品的太赫兹吸收谱中进行选取,以得到初始种群X和交叉种群U中每个个体相对应的经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱;
所述的适应度计算模块用于构建适应度函数,并利用所构造的适应度函数分别计算初始种群X和交叉种群U中每个个体的适应度;
所述的比较模块用于对初始种群X和交叉种群U中相对应个体的适应度值进行比较,将适应度值较大的个体保留下来,从而得到新一代种群Xnext;
所述的选择模块用于将新一代种群Xnext作为新的初始种群,重复执行变异交叉操作模块、选取模块、适应度计算模块和比较模块,直至进化代数达到设定阈值,并将最终代种群中适应度值最高的个体作为所选择的谷氨酰胺太赫兹吸收谱波长的最优解。
进一步地,所述变异交叉操作模块所采用变异操作如下:
其中Vi,j表示变异种群V中第i个个体的第j个二进制元素,Xi,a、Xi,b和Xi,c分别表示初始种群X中第i个个体的第a,b和c个二进制元素,并且a,b和c是3个互不相同的正整数。
进一步地,所述变异交叉操作模块所采用交叉操作如下:
其中Ui,j表示交叉种群U中第i个个体的第j个二进制元素,rand为0至1区间内的随机数,CR为预先设定的交叉概率,jrand为1至L区间内的随机整数,Xi,j表示初始种群X中第i个个体的第j个二进制元素。
本发明的有益效果是:本发明首先对初始种群X进行基于差分的变异操作和交叉操作,以得到交叉种群U;然后分别利用初始种群X和交叉种群U从谷氨酰胺样品的太赫兹吸收谱中进行选取,利用所构造的适应度函数分别计算初始种群X和交叉种群U中每个个体的适应度;并将适应度值较大的个体保留下来,从而得到新一代种群Xnext;最后将新一代种群Xnext作为新的初始种群进行进化迭代,直至进化代数达到设定阈值,并将最终代种群中适应度值最高的个体作为所选择的谷氨酰胺太赫兹吸收谱波长的最优解。本发明通过对谷氨酰胺样品的太赫兹吸收谱进行逐点深度选择,挑选有用信息,从而显著提高了定量分析的精度,取得了良好的应用效果。
附图说明
图1是基于差分进化的谷氨酰胺太赫兹谱波长选择方法的流程图;
图2是未经波长选择的谷氨酰胺样品的太赫兹吸收谱图;
图3是波长选择后的重构谷氨酰胺太赫兹吸收谱图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步的说明。
本发明基于差分进化的谷氨酰胺太赫兹谱波长选择方法的实施例
本发明首先对初始种群X进行基于差分的变异操作和交叉操作,以得到交叉种群U;然后分别利用初始种群X和交叉种群U从谷氨酰胺样品的太赫兹吸收谱中进行选取,利用所构造的适应度函数分别计算初始种群X和交叉种群U中每个个体的适应度;并将适应度值较大的个体保留下来,从而得到新一代种群Xnext;最后将新一代种群Xnext作为新的初始种群进行进化迭代,直至进化代数达到设定阈值,并将最终代种群中适应度值最高的个体作为所选择的谷氨酰胺太赫兹吸收谱波长的最优解。该方法的流程如图1所示,具体实施过程如下。
1.随机生成一个大小为S的初始种群X。
本实施例中生成的初始种群X由S个长度为L的二进制字符串组成,该二进制字符串与谷氨酰胺样品的太赫兹吸收谱中的L个频率点一一对应。
2.对初始种群x执行基于差分思想的变异操作,从而得到变异种群V。
本实施例中的变异操作后得到的种群V为:
其中Vi,j表示变异种群V中第i个个体的第j个二进制元素,Xi,a、Xi,b和Xi,c分别表示初始种群X中第i个个体的第a,b和c个二进制元素,并且a,b和c是3个互不相同的正整数。
3.对变异种群进行交叉操作,得到交叉种群U。
交叉操作得到的交叉种群U为:
其中Ui,j表示交叉种群U中第i个个体的第j个二进制元素,rand为0至1区间内的随机数,CR为预先设定的交叉概率,jrand为1至L区间内的随机整数,Xi,j表示初始种群X中第i个个体的第j个二进制元素。
4.分别利用初始种群X和交叉种群U中的个体对谷氨酰胺太赫兹吸收谱进行波长选择,得到重构吸收谱。
波长选择通过以下方式进行,对于初始种群X中的各个个体,若其某个二进制元素为“1”,则对应谷氨酰胺太赫兹吸收谱的频率点数据被保留,否则该频率点数据则被抛弃,然后将所有保留下的频率点数据整合在一起,组成经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱;;对于交叉种群U中的各个个体,若其某个二进制元素为“1”,则对应谷氨酰胺太赫兹吸收谱的频率点数据被保留,否则该频率点数据则被抛弃,然后将所有保留下的频率点数据整合在一起,组成经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱。
5.构建适应度函数F。
构建的适应度函数用下列公式表示:
其中F是适应度值,m是校正集中谷氨酰胺样品的总数量(校正集是由若干个成分浓度信息已知的谷氨酰胺样品组成的),qe是每个谷氨酰胺样品对应的定量分析误差,n代表校正集中谷氨酰胺样品的某一个:
其中ccal和creal分别是谷氨酰胺样品的计算浓度和真实浓度;谷氨酰胺样品的计算浓度ccal是利用经过波长选择后的重构太赫兹吸收谱,通过最小二乘线性回归得到,谷氨酰胺样品的真实浓度creal是预先配制的。
6.利用构建的适应度函数F分别对初始种群X和交叉种群U中的个体进行评价,得到它们各自的适应度值。
这里的适应度值是分别利用初始种群X和交叉种群U中的各个个体对谷氨酰胺样品的太赫兹吸收谱进行波长选择,而后通过适应度函数计算种群中各个个体的适应度值。
7.对初始种群X和交叉种群U中相对应个体的适应度值进行比较,将适应度值较大的个体保留下来,从而得到新一代种群Xnext。
本实施例中新一代种群Xnext需满足下式:
其中Xi next为新一代种群Xnext的第i个个体,F(Ui)为交叉种群U中第i个个体的适应度值,F(Xi)为初始种群X中第i个个体的适应度值。
8.将新一代种群Xnext作为新的初始种群,重复步骤2至7,直至进化代数达到上限G,终止进化,并将最终代种群中适应度值最高的个体作为问题的最优解输出。
本发明基于差分进化的谷氨酰胺太赫兹谱波长选择装置的实施例
本实施例中的选择装置包括初始种群生成模块、变异交叉操作模块、太赫兹吸收谱重构模块、适应度计算模块、比较模块和选择模块;初始种群生成模块用于随机生成一个大小为S的初始种群X;变异交叉操作模块用于对初始种群X进行基于差分的变异操作,以得到变异种群V,并对变异种群V进行交叉操作,以得到交叉种群U;太赫兹吸收谱重构模块用于分别利用初始种群X和交叉种群U从谷氨酰胺样品的太赫兹吸收谱中进行选取,以得到初始种群X和交叉种群U中每个个体相对应的经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱;适应度计算模块用于构建适应度函数,并利用所构造的适应度函数分别计算初始种群X和交叉种群U中每个个体的适应度;比较模块用于对初始种群X和交叉种群U中相对应个体的适应度值进行比较,将适应度值较大的个体保留下来,从而得到新一代种群Xnext;选择模块用于将新一代种群Xnext作为新的初始种群,重复执行变异交叉操作模块、选取模块、适应度计算模块和比较模块,直至进化代数达到设定阈值,并将最终代种群中适应度值最高的个体作为所选择的谷氨酰胺太赫兹吸收谱波长的最优解。
这里的波长选择装置可以采用单片机、DSP、PLC或MCU等实现,波长选择装置执行有上述六个模块,这里的模块可以位于RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、移动磁盘、CD-ROM或者本领域已知的任何其他形式的存储介质,可以将该存储介质耦接至波长选择装置,使波长选择装置能够从该存储介质读取信息,或者该存储介质可以是波长选择装置的组成部分。各模块的具体实现手段已在方法的实施例中进行了详细说明,这里不再赘述。
为了验证本发明的效果,下面设计了一系列定量分析的实验。实验选取了9个不同含量的谷氨酰胺样品的太赫兹吸收谱(其中前6个为校正集,后3个为验证集),分别利用不经选择的谷氨酰胺全吸收谱以及经过本发明提出的波长选择方法选择后的谷氨酰胺重构太赫兹吸收谱对谷氨酰胺样品进行定量分析,谷氨酰胺样品含量以及定量分析的误差如表1所示。本实验中,谷氨酰胺样品(具体包括谷氨酸和谷氨酰胺)的原始太赫兹吸收谱范围为0.3-3THz,分辨率约为4.5GHz,共有590个频率点,所以种群中二进制字符串个体的长度为590,种群大小为100,交叉概率为0.9,进化代数上限为30000。
表1
未经波长选择的谷氨酰胺样品的太赫兹吸收谱如图2所示,经本发明波长选择后的重构谷氨酰胺太赫兹吸收谱如图3所示,可见利用本发明提出的波长选择方法,能够有效降低对谷氨酰胺样品太赫兹吸收谱进行定量分析的误差,误差大致在3%以下,取得了优异的效果。
Claims (10)
1.一种基于差分进化的谷氨酰胺太赫兹谱波长选择方法,其特征在于,该选择方法包括以下步骤:
1)随机生成一个大小为S的初始种群X;
2)对初始种群X进行基于差分的变异操作,以得到变异种群V,并对变异种群V进行交叉操作,以得到交叉种群U;
3)分别利用初始种群X和交叉种群U从谷氨酰胺样品的太赫兹吸收谱中进行选取,以得到初始种群X和交叉种群U中每个个体相对应的经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱;
4)构建适应度函数,利用所构造的适应度函数分别计算初始种群X和交叉种群U中每个个体的适应度;
5)对初始种群X和交叉种群U中相对应个体的适应度值进行比较,将适应度值较大的个体保留下来,从而得到新一代种群Xnext;
6)将新一代种群Xnext作为新的初始种群,重复步骤2)-5),直至进化代数达到设定阈值,并将最终代种群中适应度值最高的个体作为所选择的谷氨酰胺太赫兹吸收谱波长的最优解。
2.根据权利要求1所述的基于差分进化的谷氨酰胺太赫兹谱波长选择方法,其特征在于,所述步骤2)中的变异操作如下::
i=1,2,…,S
j=1,2,…,L
其中Vi,j表示变异种群V中第i个个体的第j个二进制元素。
3.根据权利要求1或2所述的基于差分进化的谷氨酰胺太赫兹谱波长选择方法,其特征在于,所述步骤2)中的交叉操作为:
其中Ui,j表示交叉种群U中第i个个体的第j个二进制元素,rand为0至1区间内的随机数,CR为预先设定的交叉概率,jrand为1至L区间内的随机整数,Xi,j表示初始种群X中第i个个体的第j个二进制元素。
4.根据权利要求1所述的基于差分进化的谷氨酰胺太赫兹谱波长选择方法,其特征在于,所述步骤1)中的初始种群X由S个长度为L的二进制字符串组成,该二进制字符串与谷氨酰胺样品的太赫兹吸收谱中的L个频率点一一对应。
5.根据权利要求1所述的基于差分进化的谷氨酰胺太赫兹谱波长选择方法,其特征在于,所述步骤3)中的波长选择过程如下:
对于初始种群X中的各个个体,若其某个二进制元素为“1”,则对应谷氨酰胺太赫兹吸收谱的频率点数据被保留,否则该频率点数据则被抛弃,将所有保留下的频率点数据整合在一起,组成经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱;对于交叉种群U中的各个个体,若其某个二进制元素为“1”,则对应谷氨酰胺太赫兹吸收谱的频率点数据被保留,否则该频率点数据则被抛弃,将所有保留下的频率点数据整合在一起,组成经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱。
6.根据权利要求1所述的基于差分进化的谷氨酰胺太赫兹谱波长选择方法,其特征在于,所述步骤4)中构建的适应度函数为:
其中F是适应度值,m是校正集中谷氨酰胺样品的总数量,qe是每个谷氨酰胺样品对应的定量分析误差,n代表校正集中谷氨酰胺样品的某一个,其中ccal和creal分别是谷氨酰胺样品的计算浓度和真实浓度。
7.根据权利要求1所述的基于差分进化的谷氨酰胺太赫兹谱波长选择方法,其特征在于,所述新一代种群Xnext满足下式
其中Xi next为新一代种群Xnext的第i个个体,F(Ui)为交叉种群U中第i个个体的适应度值,F(Xi)为初始种群X中第i个个体的适应度值。
8.一种基于差分进化的谷氨酰胺太赫兹谱波长选择装置,其特征在于,该选择装置包括初始种群生成模块、变异交叉操作模块、太赫兹吸收谱重构模块、适应度计算模块、比较模块和选择模块,
所述初始种群生成模块用于随机生成一个大小为S的初始种群X;
所述的变异交叉操作模块用于对初始种群X进行基于差分的变异操作,以得到变异种群V,并对变异种群V进行交叉操作,以得到交叉种群U;
所述的太赫兹吸收谱重构模块用于分别利用初始种群X和交叉种群U从谷氨酰胺样品的太赫兹吸收谱中进行选取,以得到初始种群X和交叉种群U中每个个体相对应的经过波长选择的谷氨酰胺样品的重构太赫兹吸收谱;
所述的适应度计算模块用于构建适应度函数,并利用所构造的适应度函数分别计算初始种群X和交叉种群U中每个个体的适应度;
所述的比较模块用于对初始种群X和交叉种群U中相对应个体的适应度值进行比较,将适应度值较大的个体保留下来,从而得到新一代种群Xnext;
所述的选择模块用于将新一代种群Xnext作为新的初始种群,重复执行变异交叉操作模块、选取模块、适应度计算模块和比较模块,直至进化代数达到设定阈值,并将最终代种群中适应度值最高的个体作为所选择的谷氨酰胺太赫兹吸收谱波长的最优解。
9.根据权利要求8所述的基于差分进化的谷氨酰胺太赫兹谱波长选择装置,其特征在于,所述变异交叉操作模块所采用变异操作如下:
i=1,2,…,S
j=1,2,…,L
其中Vi,j表示变异种群V中第i个个体的第j个二进制元素。
10.根据权利要求8所述的基于差分进化的谷氨酰胺太赫兹谱波长选择装置,其特征在于,所述变异交叉操作模块所采用交叉操作如下:
其中Ui,j表示交叉种群U中第i个个体的第j个二进制元素,rand为0至1区间内的随机数,CR为预先设定的交叉概率,jrand为1至L区间内的随机整数,Xi,j表示初始种群X中第i个个体的第j个二进制元素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611103971.4A CN106769992B (zh) | 2016-12-05 | 2016-12-05 | 基于差分进化的谷氨酰胺太赫兹谱波长选择方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611103971.4A CN106769992B (zh) | 2016-12-05 | 2016-12-05 | 基于差分进化的谷氨酰胺太赫兹谱波长选择方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106769992A true CN106769992A (zh) | 2017-05-31 |
CN106769992B CN106769992B (zh) | 2020-03-24 |
Family
ID=58874058
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611103971.4A Active CN106769992B (zh) | 2016-12-05 | 2016-12-05 | 基于差分进化的谷氨酰胺太赫兹谱波长选择方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106769992B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112414996A (zh) * | 2020-07-24 | 2021-02-26 | 北京工商大学 | 基于有限差分和差分进化算法的icp-aes谱线重叠干扰校正方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7434216B1 (en) * | 2002-11-25 | 2008-10-07 | Hewlett-Packard Development Company, L.P. | Update package generator that employs genetic evolution to determine bank order |
US20080307399A1 (en) * | 2007-06-05 | 2008-12-11 | Motorola, Inc. | Gene expression programming based on hidden markov models |
CN102279390A (zh) * | 2011-05-06 | 2011-12-14 | 西南交通大学 | 一种低信噪比雷达辐射源信号脉内调制识别方法 |
CN102325089A (zh) * | 2011-09-19 | 2012-01-18 | 西安电子科技大学 | 基于差分进化和捕食搜索策略的胖树型片上网络映射方法 |
CN104102956A (zh) * | 2014-07-23 | 2014-10-15 | 杭州中威智能科技有限公司 | 一种基于策略自适应差分进化的配电网络扩展规划方法 |
CN104730025A (zh) * | 2015-04-01 | 2015-06-24 | 河南工业大学 | 一种基于太赫兹光谱的混合物定量分析方法 |
CN105136714A (zh) * | 2015-09-06 | 2015-12-09 | 河南工业大学 | 一种基于遗传算法的太赫兹光谱波长选择方法 |
-
2016
- 2016-12-05 CN CN201611103971.4A patent/CN106769992B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7434216B1 (en) * | 2002-11-25 | 2008-10-07 | Hewlett-Packard Development Company, L.P. | Update package generator that employs genetic evolution to determine bank order |
US20080307399A1 (en) * | 2007-06-05 | 2008-12-11 | Motorola, Inc. | Gene expression programming based on hidden markov models |
CN102279390A (zh) * | 2011-05-06 | 2011-12-14 | 西南交通大学 | 一种低信噪比雷达辐射源信号脉内调制识别方法 |
CN102325089A (zh) * | 2011-09-19 | 2012-01-18 | 西安电子科技大学 | 基于差分进化和捕食搜索策略的胖树型片上网络映射方法 |
CN104102956A (zh) * | 2014-07-23 | 2014-10-15 | 杭州中威智能科技有限公司 | 一种基于策略自适应差分进化的配电网络扩展规划方法 |
CN104730025A (zh) * | 2015-04-01 | 2015-06-24 | 河南工业大学 | 一种基于太赫兹光谱的混合物定量分析方法 |
CN105136714A (zh) * | 2015-09-06 | 2015-12-09 | 河南工业大学 | 一种基于遗传算法的太赫兹光谱波长选择方法 |
Non-Patent Citations (4)
Title |
---|
个人图书馆: "《进化算法、遗传算法与粒子算法之间的比较》", 21 November 2016 * |
刘林等: "《现代飞行控制系统的评估与确认方法》", 30 November 2010 * |
孔祥勇等: "无参数变异的二进制差分进化算法", 《东北大学学报( 自然科学版)》 * |
张金玉等: "《红外热波检测及其图像序列处理技术》", 30 June 2015, 国防工业出版社 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112414996A (zh) * | 2020-07-24 | 2021-02-26 | 北京工商大学 | 基于有限差分和差分进化算法的icp-aes谱线重叠干扰校正方法 |
CN112414996B (zh) * | 2020-07-24 | 2022-06-17 | 北京工商大学 | 基于有限差分和差分进化算法的icp-aes谱线重叠干扰校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106769992B (zh) | 2020-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103528990B (zh) | 一种近红外光谱的多模型建模方法 | |
CN105136714B (zh) | 一种基于遗传算法的太赫兹光谱波长选择方法 | |
CN101915744B (zh) | 物质成分含量的近红外光谱无损检测方法及装置 | |
CN104730025B (zh) | 一种基于太赫兹光谱的混合物定量分析方法 | |
CN104020135B (zh) | 基于近红外光谱的校正模型建模方法 | |
CN106124449A (zh) | 一种基于深度学习技术的土壤近红外光谱分析预测方法 | |
CN105225223A (zh) | 基于小波分析和bp神经网络的复合材料损伤检测方法 | |
CN106198482A (zh) | 基于拉曼光谱的检测保健品中是否添加有西药的方法 | |
CN104374738A (zh) | 一种基于近红外提高鉴别结果的定性分析方法 | |
CN110006844A (zh) | 基于函数性主元分析的近红外光谱特征提取方法和系统 | |
CN107169492A (zh) | 基于fcn‑crf主从网络的极化sar目标检测方法 | |
CN102499670A (zh) | 基于稳健估计和固态模函数的心电基线漂移校正方法 | |
CN105784628A (zh) | 一种利用中红外光谱检测土壤有机质化学组成的方法 | |
CN105808689A (zh) | 一种基于人工神经网络的水系实体语义相似性度量方法 | |
CN105758819A (zh) | 一种利用近红外光谱检测土壤的有机组分的方法 | |
CN106769992A (zh) | 基于差分进化的谷氨酰胺太赫兹谱波长选择方法及装置 | |
US20210396729A1 (en) | Small area real-time air pollution assessment system and method | |
CN111428585A (zh) | 一种基于深度学习的超材料太赫兹谱学识别方法 | |
CN111323407A (zh) | 一种快速测定中药多糖含量的拉曼光谱检测方法 | |
CN107796766A (zh) | 一种臭脚盐产地鉴别方法、装置及计算机可读存储介质 | |
CN106706551A (zh) | 基于差分进化的组氨酸太赫兹吸收谱波长选择方法及装置 | |
CN106940290A (zh) | 基于高光谱数据的水库水体氨氮含量估算方法 | |
CN105911003A (zh) | 基于rbm回归的水质toc浓度分析方法 | |
CN108120694A (zh) | 用于晒红烟化学成分分析的多元校正方法及系统 | |
CN116380869A (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 |