[发明专利]基于模糊划分和模型集成的气体传感器阵列浓度检测方法

摘要

一种基于模糊划分和模型集成的气体传感器阵列浓度检测方法，属于气体传感器阵列信号处理技术领域，该方法借助模糊聚类方法对基线漂移数据进行时间段划分，将原始数据集划分成漂移程度不同的多个子数据集；然后建立不同训练数据集的回归模型，得到若干子回归模型；在训练集内获得每个子回归模型的最优权重集合，对聚类中心和最优权重进行拟和，得到最优权重拟和函数；在测试阶段，基于最优权重拟和函数和聚类中心时间计算拟和权重，将这些子回归模型对待测数据的预测结果进行集成，得到最终气体浓度值。此方法能自适应地改变模式识别模型，使其能跟踪漂移的变化，有效减小漂移对浓度检测性能的影响，保证浓度测量的长期准确性。

主权项

1.基于模糊划分和模型集成的气体传感器阵列漂移补偿方法，其步骤如下：步骤1：数据采集：采用样气流量系统测量方法，采集电子鼻检测系统中金属氧化物半导体气体传感器阵列的时间响应信号，提取在标准气体环境下的基线响应数据X₀(t)和不同浓度的单一或混合气体下响应信号的特征数据x(t)＝[x₁(t),x₂(t),…,x_K(t)]，其中，K为气体传感器的个数；步骤2：数据预处理：对步骤1采集的基线响应数据进行预处理，包括异常值剔除和标准化两个步骤：异常值剔除的方法是首先计算样本均值和样本的标准偏差分别如下：然后判断基线响应数据与样本均值之差是否大于标准偏差的3倍，即若大于,则剔除该基线数据样本,并使用样本平均值填充此数据，若小于，则保留；采用阿特曼方法标准化后的数据为其中其中，i＝1,2,3，…，N,j＝1,2,3，…,K；其中，N是基线响应数据的采样个数，即为样本个数；步骤3：在预处理后的基线响应数据中加入对应的时间向量t＝[t₁,t₂,…,t_N]^T，得到增广的基线矩阵数据令则t_i时刻的基线数据可表示为步骤4：采用聚类方法对步骤3得到的增广的基线矩阵数据进行分析，具体步骤如下：步骤①：随机选取C个元素，作为C个类的中心，初始化隶属函数矩阵，给定聚类中心改变量的阈值u_ε；步骤②：确定目标函数；目标函数定义为N为样本的个数，u_ij为第i个数据点属于第j类的隶属度；β_j是第j类的中心；z为隶属度的加权指数，1≤z＜∞；d(x₀(t_i),β_j)为相似度度量方式；步骤③：确定相似度度量方式：设计相似度度量方式为式中ω＝[ω₀,ω₁,ω₂,…,ω_K]是与输入属性相对应的权重向量，ω_k∈[0,1],其值可通过实验进行选择，ω₀表示时间输入属性的向量权重，取值为1；步骤④：最小化目标函数，重新计算C个类各自的中心和隶属度函数；模糊分组的迭代过程中隶属度函数与聚类中心的更新分别为步骤⑤：重复步骤④，直到聚类结果不再变化或相对上次得到的聚类中心和隶属度的改变量小于预先确定的阈值，得到最终确定的聚类中心β_j和隶属度u_ij；步骤⑥：由于增广的基线矩阵数据的第一列元素为响应时间，因此相对应地取每个聚类中心向量β_j的第一个元素按照从小到大排列构成[T₁,T₂,…,T_C]，向量中各元素代表每个时间段划分的聚类中心时刻，体现了在每个时间段内传感器漂移随时间变化的程度是不同的，然后计算与聚类中心时刻相对应的扩展宽度则与聚类中心时刻对应的时间段范围为[T_j‑η·δ_j,T_j+η·δ_j]，η为区间的重叠参数，调整η可以控制时间段的重叠程度；步骤⑦：根据聚类方法得到的时间段范围[T_j‑η·δ_j,T_j+η·δ_j]，将基线响应数据划分为C个不同时段的数据集D₁,D₂,…,D_C，每个数据集的样本数为步骤5：选取D₁,D₂,…,D_L作为训练集，L＜C，D_L+1,D_L+2,…,D_C作为测试集；在每个训练数据集内，采用加权多输出支持向量回归(WM‑SVR)方法得到回归函数在每个训练集内，设这个训练集的样本数为N_Dm，加权多输出支持向量回归方法为，对回归函数目标是寻找W＝[w₁,w₂,…,w_Q]，b＝[b₁,b₂,…,b_Q]^T,求解优化问题:其中，J(W,b)为目标函数；ε为松弛因子；L(υ_i)为ε不敏感损失函数，y_i为单一或混合气体浓度实际值，Q为输出变量个数，即气体种类数，c为惩罚因子，λ_i选取的是步骤4得到的样本点对于当前聚类的隶属度，e_i为目标函数与各分量的误差；步骤6：在每个训练集内按步骤5计算得到的回归函数设为把这些回归函数加权组合得到模型集成的回归预测模型：然后确定每个训练集的权重系数方法为：对每个权重系数在[0,1]区间内进行搜索，搜索间隔为0.02，计算样本浓度实际值y_i与预测值y_model,i之间的均方根误差对权重系数在其搜索范围内的每一种情况得到的均方根误差RMSE最小的那一组为最优权重系数；步骤7：构建权重的拟和函数g_j(t)，j＝1，2，…,L，用以确定以下测试集的回归预测模型的权重系数：本步骤方法具体为：采用步骤6中得到的最优权重系数αⁱ作为机器学习过程的结果，最优权重系数αⁱ对应的各训练集的中心时刻T_i作为数据，即此学习过程的训练样本为基于最小二乘支持向量回归的方法得到权重的拟和函数g_j(t)；最小二乘支持向量回归方法中选取径向基核函数作为核函数；核参数δ、惩罚因子、不敏感参数寻优算法采用基于粒子群智能全局优化算法的参数选择方法；步骤8：测试阶段：根据步骤7得到的权重拟和函数g_j(t)计算拟和权重g₁(T_i),g₂(T_i),…,g_L(T_i)；然后基于各子回归模型与拟和权重进行回归模型的集成i＝L,L+1,…,C,得到每个测试集内单一或混合气体的浓度预测最终值。