一种基于最小二乘支持向量机的水松纸透气度检测方法
【技术领域】
[0001] 本发明涉及一种基于最小二乘支持向量机的水松纸透气度检测方法，属于信息处 理技术领域。
【背景技术】
[0002] 当前，随着国家经济和社会的发展，人们对健康问题的日益关注，自我保护意识不 断加强。我们知道吸烟对健康有很大的危害性，为了减轻吸烟对人体健康的危害，烟草行业 投入了大量的人力物力进行研究。
[0003] 打孔水松纸透气度的大小是有效控制卷烟中焦油含量的一个重要指标，而焦油量 是香烟对人体健康造成伤害的一个重要的指标，我国是一个烟草大国，拥有世界上最多的 烟民，研究高效可靠的水松纸透气度检测方法对国民经济和人民生活水平的提高都具有特 别重要的意义。
[0004] 现有的水松纸检测方法有如下几种：
[0005] ( -）离线检测
[0006] 本特生法检测原理：在一种规定的压差下，在一定的时间内测量通过试样一定的 表面积的空气量。通常试样的测定表面积为2cm2，透气度最终精确到1CU。
[0007] 肖伯尔法检测原理：在规定的压差下，测试一定时间透过试样的空气体积。测试的 面积为 2. 0±0. 2cm2。
[0008] 通过对本特生法和肖伯尔法这两种传统方法的分析和调研可知，它们在实际应用 中已经取得了很好的效果，但是仍然存在一些不足。首先，传统的方法需要使空气压缩机产 生的气体通过水松纸，然后通过检测水松纸两端的压差来计算透气度，故所需要较多的硬 件设备，导致造价比较昂贵；其次，传统方法的测量速度比较慢，最快也是要到4s，不利于 提高生产率；此外，采用这种方法设计的透气度检测是只能用离线测量。
[0009] (二）在线检测
[0010] 国内使用的水松纸透气度在线检测仪大多来自德国和瑞典，德国BASTAN公司生 产的在线透气度检测仪，其原理是使用普通光源光穿过孔后被聚焦到光电二极管阵列，使 光信号转变为电信号，该电信号即为与透气度相对应的一个测量值，但光源的不稳定性及 光的不均匀性，会使测量产生误差。瑞典L&W公司的2000型在线透气度检测仪，其原理则是 使移动的纸张紧贴到测量辐板表面，在辐板上一定面积的测量槽内施加恒定吸力，然后测 量一定体积的空气穿透纸张所需的时间，再将时间转换成与透气度指标相对应的参数。由 于纸张移动速度极快。这种测量方法会产生较大的测量误差。
【发明内容】
[0011] 本发明提供了一种基于最小二乘支持向量机的水松纸透气度检测方法，以用于现 有的水松纸透气度检测设备造价昂贵，测量速度慢及在线检测测量误差较大的问题。
[0012] 本发明基于最小二乘支持向量机的水松纸透气度检测方法是这样实现的：所述基 于最小二乘支持向量机的水松纸透气度检测方法的具体步骤如下：
[0013] S1、获取已知的打孔水松纸透气度样本数据；
[0014] S2、对采集的数据进行归一化处理，得到样本集；
[0015] S3、将样本集的数据分为两部分，一部分作为LS-SVM(最小二乘支持向量机）训练 样本，一部分作为测试样本，并通过PSO (粒子群优化算法）算法对LS-SVM参数进行优化， 即粒子群算法对LS-SVM参数进行优化；将训练样本用于支持向量回归机的训练，得到最小 二乘支持向量机数学模型；测试样本用于测试模型的精度；
[0016] S4、获得LS-SVM最优参数；
[0017] S5、利用训练样本集，对优化参数的LS-SVM进行训练，得到最优化的LS-SVM数学 模型，即打孔水松纸透气度检测模型；
[0018] S6、采用建立的打孔水松纸透气度检测模型，对测试样本进行检测。
[0019] 所述Sl步骤中，所述已知的打孔水松纸透气度样本数据包括水松纸的打孔面积 及水松纸透气度。
[0020] 所述S2步骤中，所述归一化处理的公式为：
[0021]
[0022] 其中，Sk为获取的第k个打孔水松纸的打孔面积，s_为所述Sk中的最小值，s_ 为所述s k中的最大值，η为所述打孔水松纸透气度样本的数目，g (s k)为％归一化处理后得 到的数据。
[0023] 所述S3步骤中，所述的将样本集的数据分为两部分包括：随机选择样本集中2/3 的数据为训练样本，构成训练样本集；其余1/3的数据为测试样本，构成测试样本集。
[0024] 所述S3步骤中，所述的通过PSO算法对LS-SVM参数进行优化，其具体步骤如下：
[0025] S3. 1、初始化粒子群，即随机产生一组初始化值，包括微粒的速度和位置；设定粒 子群参数，在空间Rn中随机产生η个粒子X i，x2, ...，Xn，组成初始化种群X (t);随机产生各 粒子的初始速度V1, v2, ...，Vn，组成速度矩阵V (t);每个粒子的个体最优值P1的初始值为 Xi的初始值，所述Pi= (Pil, Pi2);
[0026] S3. 2、计算所述粒子群中所有粒子的适应度；
[0027] 其中，所述适应度值公式如下：
[0028]
[0029] 其中，y]为训练样本集中的第j个样本的透气度真实值，^为训练样本集中的第j 个样本的透气度预测值，m为训练样本集中样本向量的个数；
[0030] S3. 3、对于每个粒子，将当前各微粒的适应度值f (X1)与自身的最优值适应度值 f (P1)进行比较，若f (X1) <f (P1)，则用适应度值f (X1)更新f (P1);
[0031] S3. 4、对于每个粒子，将每个粒子的最优适应度值f (X1)与所有粒子的最优适应度 值f(Pg)进行比较，若f(Xl) <f(Pg)，则用适应度值f(Xl)更新f(Pg)，同时保存粒子的当前 位置和速度；
[0032] S3. 5、更新粒子位置、速度，产生新种群X(t+1)，其中，速度调整方式为：
[0033] 当 V1 > v nax时，视为 v i = v nax;当 v A -v nax时，视为 v i = -v nax;
[0034] 其中，所述的更新粒子的速度和位置公式如下：
[0035]
[0036]
[0037] 其中，ω为惯性权因子，(^和c 2为正的加速常数，r挪r 2为在0到1之间均匀分 布的随机数，另外，通过设置微粒的速度区间[-v_，v_]和位置范围[1_，1_]，则对微粒 的移动进行限制；
[0038] 对种群中每一个粒子i而言，它的位置表示为X1= (X u，xl2, ...，xld)，它的飞行速 度表示为Vi= (V vi2, · · ·，vid)，粒子自身经历过的最优位置记为Pi= (P Pi2, · · ·，Pid)， 种群中所有粒子目前为止所经历过的最优位置记为Pg= (p gl，pg2,...，Pgd)。
[0039] 对于公式中i，j的解释：其中j = 1，2, . . .，d，d为解空间的维数，（由于本方法中 只需要优化LS-SVM中的两个参数，故d = 2) ;i = 1，2, ...，η，η为粒子种群规模，即随机产 生的粒子数，（本发明中可取η = 20)。
[0040] 其中，更新惯性因子ω公式如下：
[0041 ]
[0042]其中，
[0043] Tniax为最大迭代次数，ω e [ ω _，ω _]，将ω初始化为〇· 9,并使其随迭代次数增 加而线性递减至0. 4 ;
[0044] S3. 6、检查停止条件，如果达到最大迭代次数Tniax则停止迭代；否则，返回步骤 S3. 2 ；
[0045] S3. 7、输出最优参数。
[0046] 本发明的有益效果是：
[0047] 基于最小二乘支持向量机水松纸透气度检测方法对硬件需求低，结构简单，具有 良好的可靠性。
[0048] 基于最小二乘支持向量机的水松纸透气度检测方法能够对其进行快速，准确的拟 合预测。
[0049] 基于最小二乘支持向量机的水松纸透气度检测方法，提高了模型精度和泛化能 力，克服了现有检测方法精度低的缺点，对实现水松纸透气度在线检测和智能控制具有重 要意义。
【附图说明】
[0050] 图1为本发明中的流程图；
[0051] 图2为本发明采用PSO-LS-SVM模型的测试集预测结果图；
[0052] 图3为本发明采用PSO-LS-SVM模型的测试集误差结果图。
【具体实施方式】
[0053] 实施例1 :如图1-3所示，一种基于最小二乘支持向量机的水松纸透气度检测方 法，所述基于最小二乘支持向量机的水松纸透气度检测方法的具体步骤如下：
[0054] S1、获取已知的打孔水松纸透气度样本数据；
[0055] S2、对采集的数据进行归一化处理，得到样本集；
[0056] S3、将样本集的数据分为两部分，一部分作为LS-SVM训练样本，一部分作为测试 样本，并通过PSO算法对LS-SVM参数进行优化，即粒子群算法对LS-SVM参数进行优化；将 训练样本用于支持向量机的训练，得到最小二乘支持向量机数学模型；测试样本用于测试 模型的精度；
[0057] S4、获得LS-SVM最优参数；
[0058] S5、利用训练样本集，对优化参数的LS-SVM进行训练，得到最优化的LS-SVM数学 模型，即打孔水松纸透气度检测模型；
[0059] S6、采用建立的打孔水松纸透气度检测模型，对测试样本进行检测。
[0060] 所述Sl步骤中，采集到某烟厂的108个打孔水松纸样本数据，所述已知的打孔水 松纸透气度样本数据包括水松纸的打孔面积及水松纸透气度。
[0061] 所述S2步骤中，所述归一化处理的公式为：
[0062]
[0063] 其中，Sk为获取的第k个打孔水松纸的打孔面积，s_为所述Sk中的最小值，s_ 为所述s k中的最大值，η为所述打孔水松纸透气度样本的数目，g (s k)为％归一化处理后得 到的数据。
[0064] 所述S3步骤中，所述的将样本集的数据分为两部分包括：随机选择样本集中2/3 的数据（即72个样本）为训练样本，构成训练样本集；其余1/3的数据（即36个样本）为 测试样本，构成测试样本集。
[0065] 利用训练样本集，建立最小二乘支持向量机数学模型，具体步骤如下：
[0066] 给定一个有N个训练样本的集合{xk, yk}，k = 1，2, . . .，n，其中训练样本η维向量 Xke RnJke
[0067] 首先用一个非线性映射〇把原空间样本从1^映射到特征空间Ρ(Χ,)，这样就把低 维空间的非线性逼近问题转化为高维空间的线性化逼近问题，在这个高维特征空间中构造 最优决策函数：
[0068]
[0069] 依据结构风险最小化原则，寻找ω，b就是最小化：
[0070]
[0071] 其中，I I ω I I2控制模型的复杂度；γ是正规化参数，控制对超出误差样本的惩罚 程度；ω为权矢量；b为偏差量，1