发明名称 一种快速识别发动机燃料种类的中红外光谱方法
摘要 本发明公开了一种快速识别发动机燃料种类的中红外光谱方法,该方法包括收集具有一定数量的样品作为数据库样品;测定数据库样品中红外吸收光谱,建立燃料红外光谱数据库。对于待测发动机燃料样品的识别,首先测定其红外吸收光谱,然后由计算机自动根据红外光谱数据,选择合适区间的光谱数据,采取移动窗口法,计算待测样品与数据库中各个样品的相关系数,选择相关系数最大的样品类别作为未知样品的类别。本发明方法可在数分钟内快速判断发动机种类,确定其牌号,提高燃料打假能力,避免用户误用燃料或者使用掺假燃料带来的损失。
申请公布号 CN102221534B 申请公布日期 2014.06.18
申请号 CN201110139729.3 申请日期 2011.05.27
申请人 中国人民解放军总后勤部油料研究所 发明人 田高友;褚小立;熊春华;鲁长波;周友杰;任连岭;安高军
分类号 G01N21/35(2014.01)I 主分类号 G01N21/35(2014.01)I
代理机构 北京知本村知识产权代理事务所 11039 代理人 刘江良;周自清
主权项 1.一种快速识别发动机燃料种类的中红外光谱方法,包括如下步骤:第一步:收集一定数量的样品作为数据库样品;第二步:测定数据库样品的中红外吸收光谱;第三步:将每个样品的种类信息和与其对应的红外光谱数据组合起来,建立燃料红外光谱数据库;其特征在于还包括:第四步:对于待测发动机燃料样品的识别,首先测定其红外吸收光谱,并选择合适区间的光谱数据,采取移动窗口法,计算待测样品红外光谱与数据库中各个样品红外光谱的相关系数,选择最大相关系数所对应样品的类别作为待测样品的类别,所述合适区间为719cm<sup>-1</sup>~1298cm<sup>-1</sup>和1491cm<sup>-1</sup>~2763cm<sup>-1</sup>,移动窗口法计算待测样品与数据库中样品的相关系数的过程如下:两个样品的光谱数据为X和Y,给定一移动窗口wn和波长i,则选定i-wn~i+wn区间的吸光度数据计算X和Y的相关系数CR<sub>i</sub><maths num="0001"><![CDATA[<math><mrow><msub><mi>CR</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>-</mo><mi>wn</mi></mrow><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mi>wn</mi></mrow></munderover><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>-</mo><mover><msub><mi>x</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>-</mo><mover><msub><mi>y</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><mrow><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>-</mo><mi>wn</mi></mrow><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mi>wn</mi></mrow></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>-</mo><mover><msub><mi>x</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></msqrt><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>-</mo><mi>wn</mi></mrow><mrow><mi>j</mi><mo>=</mo><mi>i</mi><mo>+</mo><mi>wn</mi></mrow></munderover><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>-</mo><mover><msub><mi>y</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></mfrac></mrow></math>]]></maths>其中,<img file="FDA0000443922520000012.GIF" wi="59" he="83" />和<img file="FDA0000443922520000013.GIF" wi="67" he="83" />分别为X和Y的光谱在波长i-wn~i+wn的吸光度平均值;x<sub>j</sub>和y<sub>j</sub>为X和Y的光谱在波长i-wn~i+wn的吸光度,wn为奇数:29;移动波长点i,每移动一个波长点后均选择该波长点±wn区间的吸光度计算相关系数CR<sub>i</sub>,直到计算完所选区间的所有波长点的CR<sub>i</sub>后,计算所有CR<sub>i</sub>的平均值记为两个样品X和Y的相关系数CR,即<img file="FDA0000443922520000014.GIF" wi="302" he="196" />n为波长数目;其中,第三步所述数据库构建过程如下:将样品类型数据矩阵Y,并列在中红外光谱数据矩阵X中,组成新的数据库矩阵Z,即:Y=[1  y<sub>1</sub>];<maths num="0002"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>w</mi><mn>1</mn></msub></mtd><mtd><msub><mi>A</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>w</mi><mn>2</mn></msub></mtd><mtd><msub><mi>A</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>w</mi><mi>n</mi></msub></mtd><mtd><msub><mi>A</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced><mo>;</mo><mi>Z</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>X</mi></mtd></mtr><mtr><mtd><mi>Y</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>w</mi><mn>1</mn></msub></mtd><mtd><msub><mi>A</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>w</mi><mn>2</mn></msub></mtd><mtd><msub><mi>A</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>w</mi><mi>n</mi></msub></mtd><mtd><msub><mi>A</mi><mi>n</mi></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><msub><mi>y</mi><mn>1</mn></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,y<sub>1</sub>为类型代码,汽油类型代码为1,柴油类型代码为2,3号喷气燃料代码为3,90号、93号和97号的车用汽油类型代码依次为10、11、12,90号、93号和97号的乙醇汽油类型代码依次为13、14、15,0号、-10号和-35号不同牌号的轻柴油类型代码依次为20、21、22,-10号和-35号的军用柴油类型代码依次为23、24;w<sub>i</sub>为第i个波长,(i=1,2,......,n)A<sub>i</sub>为第i个吸光度(i=1,2,......,n),n为波长数目;
地址 102300 北京市门头沟区三家店东街1号