[发明专利]一种物质气味嗅频提取方法

摘要

本发明公开了一种物质气味嗅频提取方法，给出了物质气味嗅频的定义；提出叠加降维算法(Superposition Mapping Analysis，SMA)对物质气味样品进行判别训练，结合物质气味化学配比表，以后序遍历规则的平衡二叉树为存储结构，建立物质气味样品属性库，对未知物质气味进行识别；构建物质气味成分比例及浓度模型对物质气味实现浓度实时提取。本发明提供了一种利用仿生嗅觉系统提取物质气味嗅频的方法，具有快速检测、高效准确、性能稳定等优点，为实现物质气味远距离网络传输提供前提基础。

主权项

一种物质气味嗅频提取方法，其特征在于，物质气味嗅频是指用于表征物质气味的特征信息，包含物质气味的种类名称记为Rn，成分记为Lc，比例浓度记为Pi，对于嗅频的提取按照以下步骤进行：步骤1：选取物质气味样品备用，样品备份为若干份，固体每份为M克、液体每份为M毫升、气体为M立方毫米，固体物质研磨成粉末或切割成体积小于1立方毫米的块状，气体和液体物质密度均匀；样品置于恒温恒湿的实验箱中静置70分钟；步骤2：应用仿生嗅觉系统对已知物质气味样品进行采集检测，测量时间t为100秒，采集速率为7.749毫升/分钟，保存所测数据Fs(S1，S2，…，SN)，并对Fs保存至计算机；步骤3：对步骤2中采集所得到的测量数据信息采用叠加映射降维算法SMA进行判断训练，并将所建立的样本序号OSMA与气味名称成分属性L(Oname，Cn1，Cn2，…，Cnk)，其中Oname是物质气味种类名称，Cn1是物质气味第1种成分名称，Cn2是物质气味第2种成分名称，Cnk是物质气味第k种成分名称；以R(OSMA，Oname，Cn1，Cn2，…，Cnk)为结点的后序遍历规则的平衡二叉树形式保存至气味成分信息库；所述步骤3中使用的叠加映射降维方法，具体步骤如下：步骤1)：将采集并测量得到的一个气味样本数据信息矩阵Fs向量化，选取其中达到稳态响应值的测量值，组成训练样本矩阵PT∈Rr×n，其中r代表训练样本矩阵的行数，由选取的样本类，样本个数及稳态时间值个数决定，j代表训练样本类别个数，j∈[1,c]，i代表第j类训练样本个数，i∈[1,ni]；则训练样本的均值为μ=Σtji/Σni---(3)]]>步骤2)：由步骤1)设获取训练样本去均值后形成的样本矩阵为T＝PT‑μ，T∈Rr×n则T的协方差矩阵为Q＝T×TT，Q∈Rr×r                               (4)其中TT是T的转置矩阵，并计算Q矩阵的特征值及特征向量，按特征值累计贡献率不小于99.5％，从大到小顺序选取前a个特征值所对应的特征向量组成第一降维特征系数矩阵Pc，Pc∈Ra×n，a≤n‑1，并将训练样本矩阵T投影到Pc中，得到第一降维识别矩阵：Pf＝T×PcT，Pf∈Rr×a                              (5)其中PcT是Pc的转置矩阵；步骤3)：将第一降维识别矩阵Pf作为第二降维的输入矩阵，即Pf∈Rr×a，其中k代表新训练样本矩阵的类别数k∈[1,c′]，m代表每个训练样本的个数m∈[1,a]，则的样本矩阵均值矩阵为训练的总样本Pf的均值矩阵μc∈R1×a，υ∈R1×a，并计算Pf类内散布矩阵Sω和类间散布矩阵Sb，即Sω=Σm=1aΣk=1c′(Vkm-μc)T×(Vkm-μc),Sω∈Ra×a---(6)]]>Sb=Σm=1a(μc-υ)T×(μc-υ),Sb∈Ra×a---(7)]]>则根据Fisher准则函数JF(ω)=ωTSbωωTSωω---(8)]]>由式(8)可知，当选取的矢量ω使JF(ω)取最大值时具有最优分析，其物理意义为以ω为投影方向，投影后的样本空间具有最大的类间离散度和最小的类内离散度；则可对式(8)应用拉格朗日乘法，设其存在特征根λ，λ是最优解即最佳投影矩阵，则有JF(ω)＝ωTSbω‑λ(ωTSωω‑1)                          (9)则对式(9)等号两边同时对ω求导，可得dJF(ω)dω=2Sbω-2λSωω=0---(10)]]>则有Sbω＝λSωω                                    (11)继而得到则对于λ的求解即可转化为求解特征矩阵的特征向量，由于Sω∈Ra×a且Sb∈Ra×a，所以求得Sω‑1Sb的特征向量Lc，Lc∈Ra×a，用来构建第二降维特征系数矩阵，最终获得叠加映射特征系数矩阵C＝Lc×Pc，C∈Ra×n                                (13)则其训练样本的叠加映射算法识别特征矩阵为：Cf＝T×CT Cf∈Rr×a，                              (14)再依据欧式距离判别式来判别n维空间的中样本点间的相似度；步骤4：输入未知物质气味样品，应用仿生嗅觉系统采集气味样品，并采用SMA算法对未知物质气味进行判别分析获得样本序号OSMA′，并在已建立的后序遍历规则的平衡二叉树内进行查找，若找到与之匹配的OSMA，则输出R(OSMA，Oname，Cn1，Cn2，…，Cnk)，若找不到与之匹配的OSMA，则返回无此信息；同时对测量值进行成分比例计算P(Tn1，Tn2，…，Tnk)，Tn1是第1种成分所占比例，Tn2是第2种成分所占比例，Tnk是第k种成分所占比例，其计算公式为Tni=A(i)Σi=1kA(i)---(1)]]>A(i)是k种成分对应传感器测量平均值，是第1种成分至第k种成分的测量平均值加权和；则物质气味浓度为Pi=Σi=1kA(i)·TniΣi=1kA(i)×V′V---(2)]]>其中V′是采集预压缩后的气味体积，V是采集的原始气味体积；步骤5：依据步骤4中所得R(OSMA，Oname，Cn1，Cn2，…，Cnk)可以求知物质气味的名称Rn及成分Lc，所得P(Tn1，Tn2，…，Tnk)可知气味的成分比例及浓度Pi，则可得到物质气味嗅频。