发明名称 一种基于神经网络的煤质特性在线测量方法
摘要 一种基于神经网络的煤质特性在线测量方法,用在激光诱导等离子光谱的煤质检测装置上。先用激光诱导等离子体系统得到定标样品的光谱,用PCA或PLS-DA提取光谱强度矩阵的主成分,在主成分得分图上把定标样品分为不同类别。随后对各类别的定标样品分别建立神经网络模型并进行训练,其中,光谱强度数据作为神经网络的输入,各元素浓度作为输出。对于未知的待测样品,先用PCA或PLS-DA判断其所属类别,然后把待测样品的光谱强度数据输入该类别的训练好的神经网络模型中,即可得到各元素浓度。该方法可实现煤质的多元素测量并减小由于煤的种类不同带来的误差,同时有较强的抗干扰能力和适应性,提高了激光诱导等离子光谱测量的准确度。
申请公布号 CN102004088A 申请公布日期 2011.04.06
申请号 CN201010537027.6 申请日期 2010.11.09
申请人 清华大学 发明人 李政;王哲;侯宗余
分类号 G01N21/63(2006.01)I;G01J3/28(2006.01)I;G06N3/00(2006.01)I 主分类号 G01N21/63(2006.01)I
代理机构 北京鸿元知识产权代理有限公司 11327 代理人 邸更岩
主权项 1.一种基于神经网络的煤质特性在线测量方法,其特征在于该方法包括如下步骤:1)首先使用各元素质量浓度已知的一组煤炭样品作为定标样品,利用安装在输煤皮带上方的激光诱导等离子光谱系统对定标样品进行检测,得到该组定标样品的光谱谱线,即得到了每种定标样品中各种元素的激光诱导等离子特征光谱谱线强度,形成谱线强度矩阵E<sub>0</sub>,E<sub>0</sub>矩阵的结构如下,<maths num="0001"><![CDATA[<math><mrow><msub><mi>E</mi><mn>0</mn></msub><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>1</mn></msub><mn>1</mn></msubsup></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>2</mn></msub><mn>1</mn></msubsup></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>3</mn></msub><mn>1</mn></msubsup></mtd><mtd><mi>L</mi></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mi>m</mi></msub><mn>1</mn></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>1</mn></msub><mn>2</mn></msubsup></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>2</mn></msub><mn>2</mn></msubsup></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>3</mn></msub><mn>2</mn></msubsup></mtd><mtd><mi>L</mi></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mi>m</mi></msub><mn>2</mn></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>1</mn></msub><mn>3</mn></msubsup></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>2</mn></msub><mn>3</mn></msubsup></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>3</mn></msub><mn>3</mn></msubsup></mtd><mtd><mi>L</mi></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mi>m</mi></msub><mn>3</mn></msubsup></mtd></mtr><mtr><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd><mtd><mi>M</mi></mtd><mtd><mi>O</mi></mtd><mtd><mi>M</mi></mtd></mtr><mtr><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>1</mn></msub><mi>n</mi></msubsup></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>2</mn></msub><mi>n</mi></msubsup></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mn>3</mn></msub><mi>n</mi></msubsup></mtd><mtd><mi>L</mi></mtd><mtd><msubsup><mi>I</mi><msub><mi>&lambda;</mi><mi>m</mi></msub><mi>n</mi></msubsup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FDA0000031337570000012.GIF" wi="49" he="67" />表示第i种样品在波长λ<sub>j</sub>处对应的谱线强度,i=1,2,L,n;j=1,2,L,m;2)利用主成分分析或者偏最小二乘判别分析方法对定标样品进行分类,两种方法都需要先提取主成分;a.主成分分析是从矩阵E<sub>0</sub>中提取主成分,首先得到矩阵E<sub>0</sub>的协方差矩阵,求取协方差矩阵的特征值,特征值从大到小依次为A<sub>1</sub>,A<sub>2</sub>,L,A<sub>h</sub>;最大特征值λ<sub>1</sub>对应的特征向量为第一主轴a<sub>1</sub>,第二特征值λ<sub>2</sub>对应的特征向量为第二主轴a<sub>2</sub>,由此求得第一、第二主成分t<sub>1</sub>、t<sub>2</sub>,t<sub>1</sub>=E<sub>0</sub>a<sub>1</sub>                                    (2)t<sub>2</sub>=E<sub>0</sub>a<sub>2</sub>                                    (3)依次类推可以求得第h个主成分t<sub>h</sub>;b.偏最小二乘判别分析法则以矩阵E<sub>0</sub>为自变量矩阵,各定标样品的挥发份含量组成因变量矩阵F<sub>0</sub>,矩阵F<sub>0</sub>的结构如下,F<sub>0</sub>=(V<sup>1</sup>  V<sup>2</sup>  V<sup>3</sup>  L  V<sup>n</sup>)<sup>T</sup>                    (4)其中,V<sup>i</sup>表示第i种定标样品的挥发份含量,i=1,2,L,n;从矩阵E<sub>0</sub>中提取主成分,使得主成分既能代表E<sub>0</sub>的变异信息,又能使E<sub>0</sub>与F<sub>0</sub>的相关程度达到最大;w<sub>1</sub>是对应于矩阵E<sub>0</sub><sup>T</sup>F<sub>0</sub>F<sub>0</sub><sup>T</sup>E<sub>0</sub>最大特征值的特征向量,得到E<sub>0</sub>的第一主成分为,t<sub>1</sub>=Xw<sub>1</sub>(5)实施E<sub>0</sub>在t<sub>1</sub>上的回归以及F<sub>0</sub>在t<sub>1</sub>上的回归,即E<sub>0</sub>=t<sub>1</sub>p′<sub>1</sub>+E<sub>1</sub>             (6)F<sub>0</sub>=t<sub>1</sub>r<sub>1</sub>+F<sub>1</sub>               (7)p′<sub>1</sub>、r<sub>1</sub>为相应的回归系数;E<sub>1</sub>、F<sub>1</sub>为残差矩阵,w<sub>2</sub>是对应于矩阵<img file="FDA0000031337570000021.GIF" wi="199" he="58" />最大特征值的特征向量,第二主成分为,t<sub>2</sub>=E<sub>1</sub>w<sub>2</sub>                  (8)依次类推可以求得第h个成分t<sub>h</sub>,最后可用交叉有效性确定偏最小二乘回归中成分t<sub>h</sub>的提取个数,停止迭代;利用方法a或方法b分别提取第一、第二主成分,两种方法分别做主成分得分图,在图上以每个定标样品对应的主成分的值作为该样本的坐标点,绘出所有坐标点的位置,观察坐标点在二维空间的分布特征和相似性结构;如果二维图的分类效果不明显,则加入第三主成分,观察坐标点在三维空间的分布特征和相似性结构;通过主成分得分图,定标样品被明显分类;3)经过主成分分析或偏最小二乘判别分析法对定标样品分类后,对各类定标样品分别建立神经网络模型,神经网络的输入层为特征谱线强度,输入层节点数目为特征谱线的数目;输出层为各元素浓度,输出层节点数目为需要测定的元素种类数;4)由于定标样品的各元素浓度已知,故用定标样品的特征谱线强度和元素浓度对神经网络模型进行训练,训练结束后则调整隐含层的层数或各层的节点数目,然后重新训练;5)重复步骤4),直至神经网络模型的训练误差达到最小,即得到最优的神经网络模型;6)对于各元素浓度未知的待测样品,首先通过安装在输送皮带上方的激光诱导等离子光谱系统得到该待测样品的特征谱线强度矩阵;然后按照步骤2)作出该待测样品在主成分得分图上的坐标点,根据该坐标点在得分图上的位置确定该待测样品所属的类别,最后将特征谱线强度数据输入到所属类别的已训练好的最优的神经网络模型中,即求得待测样品中各元素的浓度。
地址 100084 北京市海淀区北京市100084信箱82分箱清华大学专利办公室