遗传程序设计用于非分散红外吸收光谱的CO浓度测量-天津大学研究生e.doc

最小二乘拟合用于非分散红外吸收光谱的浓度测量 佟颖非分散红外吸收光谱法(NDIR)是CO浓度测量的主要手段之一。设计了NDIR的CO气体测量系统及连续测量实验装置，对10种不同浓度的标准CO气体进行测量，分析获取的0个实测数值。将用于CO光强采样数据的建模。实测数据建模结果表明，可以用于的气体浓度建模，拟合误差较小，性能较好。引言 CH4等气体的检测[5-8]。但是，目前CO浓度测量的主要方法仍然是非分散红外吸收光谱法(NDIR)[9-12]。文献[11]通过气体相关滤波(GFC)的NDIR方法设计了公路隧道的CO浓度测量系统，该方法选取的红外吸收中心波长为4.65μm，采用开放式的测量空间，是本文前期的研究基础。 最小二乘法（又称最小平方法）是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。最小二乘拟合是一种数学上的近似和优化，利用已知的数据得出一条直线或者曲线，使之在坐标系上与已知数据之间的距离的平方和最小。 (1) 其中为出射光强，为入射光强，β为气体吸收系数，L为光程，C为气体浓度。 1.2 测量系统 主要包括红外光源及光束准直、开放式测量空间、红外反射镜、GFC轮以及红外探测器等组成部分。 2 实验过程 2.1 实验装置 CO浓度连续测量的实验装置如图1所示。测量系统被固定在 3 m 长的气密舱内，A为气密舱的进气口，B为出气口，CO浓度的测量数据通过RS485总线传送至上位机。实验在天津市南开区华苑产业园区进行，标准气源为天津东祥特种气体有限公司不同浓度的CO、N2混合气体。实验开始前，首先在B处用气泵将气密舱抽成真空状态；然后从A处充入固定浓度的混合气体；充气的时间为15分钟，B处的泄压阀将气密舱内的气压保持在单位大气压；充气完成后，上位机读取并保存不同CO浓度的光强采样数值。 图 1 实验过程示意图 2.2 实验数据获取 上位机读取的原始数据是未经过GFC轮参考气室的测量光强采样数值 Dm 和经过参考气室的参考光强采样数值 Dr 。原始数据需要经过处理才能用于CO浓度的数据建模。设定光强吸收系数γ的计算式如下： (2) 上位机以1 s的周期连续接收CO测量系统传送的光强采样数值。考虑到实验装置中的CO浓度在一段时间内基本保持恒定，因此可以扩大数据的抽取间隔。在采集的过程中，实际的测量时间为90分钟。建模数据从测量数据中随机抽取。 3 数据处理与分析 此次实验的数据处理平台为Matlab，采用最小二乘线性拟合进行建模，建模所用数据为光强数据γ及一氧化碳浓度值C。 通过在Matlab环境下进行拟合（代码见附录），得到结果如下： p = 1.0e+004 * 0.8914 -1.3782 0.5331 样本数据的分布及曲线如图2所示。其中，横轴为光强吸收系数γ，纵轴为CO浓度，单位是ppm 。从图中可以看出，样本数据与曲线的一致性好。 图2 实验数据与曲线拟合结果 经模型计算，样本的拟合结果如表1所示。 表1.最小二乘的CO浓度实测数据与拟合值 True value (ppm) Fitting value (ppm) Relative error (%) 0.774285057 0 4 Inf 0.769831547 0 4 Inf 0.767627223 0 4 Inf 0.709470795 51 40 21.6 0.713334863 51 36 29.4 0.708896525 51 41 19.6 0.666973533 100 104 -4.0 0.669301005 100 100 0.0 0.666999386 100 104 -4.0 0.654649103 125 129 -3.2 0.653899507 125 130 -4.0 0.653594268 125 131 -4.8 0.645106514 154 150 2.6 0.645392808 154 150 2.6 0.645385149 154 150 2.6 0.631485348 175 183 -4.6 0.631164013 175 184 -5.1 0.630842927 175 184 -5.1 0.621791392 200 208 -4.0 0.621122829 200 210 -5.0 0.622749845 200 205 -2.5 0.607198871 254 249 2.0 0.608832808 254 244 3.9