动态演化模型校正方法及系统的制作方法
【技术领域】
[0001] 本发明涉及在线近红外光谱分析技术领域，更为具体地，涉及一种基于KF-PLS(卡 尔曼滤波偏最小二乘)近红外光谱动态演化模型校正方法及系统。
【背景技术】
[0002] 建立准确性高、自适应性强的校正模型是近红外光谱定量分析方法成功应用的关 键。现代近红外光谱以其分析速度快、重现性好、成本低、不消耗样品、易于实现在线分析等 鲜明的特点正得到越来越多的应用。在近红外光谱分析实际应用中，样品近红外光谱和待 测量存在复杂的间接关系。常通过求解标准曲线，建立测量数据与分析体系中某一待测组 分含量之间的对应关系。PLS在处理近红外光谱建模问题上显示出独特的优越性，被广泛用 于模型校正过程。
[0003] 但是，上述方法通常假设仪器的环境噪声和部件状态是稳定的，是一种静态建模 方法。其效果仅限于短期相对稳定的模型使用过程。然而光谱获取过程受被测对象变化、仪 器老化、环境扰动等影响具有时变性。当仪器测量过程慢慢发生演变时，基于早期数据的校 正模型预测效果将无法得到保证。传统方法通过利用补充样品数据库定期维护模型，从而 保证其适应性，其特点是工作量大、成本高和不易实现在线分析等。且更新的校正模型中包 含大量光谱仪早期状态信息，不能精确代表仪器最新测量状态。
[0004] 为解决上述问题，本发明需要提供一种新的模型校正方法。
【发明内容】
[0005] 鉴于上述问题，本发明的目的是提供一种基于卡尔曼滤波偏最小二乘法近红外光 谱动态演化模型校正方法及系统，以保证近红外光谱校正模型具有自适应性，降低重建模 型成本，从而实现在线分析。
[0006] 本发明提供一种基于卡尔曼滤波偏最小二乘近红外光谱动态演化模型校正方法， 包括：
[0007] 利用K/S算法从标准样品中选择有代表性的建模样品；
[0008] 采用PLS法对所述建模样品建立近红外光谱数据与浓度间的线性关系；
[0009] 所述线性关系利用去一交互检验法确定PLS因变量数，形成初始PLS校正模型，然 后通过所述标准样品确定PLS校正模型；
[0010]定期对待测样品进行化验，获取所述待测样品的样品数据；
[0011] 同时，利用所述PLS校正模型对所述待测样品进行预测，获取所述待测样品的预测 值；
[0012] 获取的所述待测样品的样品数据和所述待测样品的预测值通过采用KF算法修正 所述PLS校正模型的主因子系数。
[0013] 本发明还提供一种基于卡尔曼滤波偏最小二乘近红外光谱动态演化模型校正系 统，包括：
[0014] 建模样品选取单元，用于利用K/S算法从标准样品中选择有代表性的建模样品；
[0015] 线性关系建立单元，用于采用PLS法对所述建模样品建立近红外光谱数据与浓度 间的线性关系；
[0016] PLS校正模型建立单元，用于利用去一交互检验法确定PLS因变量数，形成初始PLS 校正模型，然后通过所述标准样品确定PLS校正模型；
[0017] 样品数据获取单元，用于定期对待测样品进行化验，获取所述待测样品的样品数 据；
[0018] 预测值获取单元，用于利用所述PLS校正模型对所述待测样品进行预测，获取所述 待测样品的预测值；
[0019] PLS校正模型的主因子系数修正单元，用于对获取的所述待测样品的样品数据和 所述待测样品的预测值通过采用KF算法修正所述PLS校正模型的主因子系数。
[0020] 从上面的技术方案可知，本发明提供的基于卡尔曼滤波偏最小二乘近红外光谱动 态演化模型校正方法及系统，利用卡尔曼滤波偏最小二乘算法以保证近红外光谱校正模型 具有自适应性，降低重建模型成本，从而达到实现在线分析的目的。
[0021] 为了实现上述以及相关目的，本发明的一个或多个方面包括后面将详细说明并在 权利要求中特别指出的特征。下面的说明以及附图详细说明了本发明的某些示例性方面。 然而，这些方面指示的仅仅是可使用本发明的原理的各种方式中的一些方式。此外，本发明 旨在包括所有这些方面以及它们的等同物。
【附图说明】
[0022] 通过参考以下结合附图的说明及权利要求书的内容，并且随着对本发明的更全面 理解，本发明的其它目的及结果将更加明白及易于理解。在附图中：
[0023] 图1为根据本发明实施例的基于卡尔曼滤波偏最小二乘近红外光谱动态演化模型 校正方法流程示意图；
[0024] 图2为根据本发明实施例的基于卡尔曼滤波偏最小二乘近红外光谱动态演化模型 校正系统结构框图。
[0025 ]在所有附图中相同的标号指示相似或相应的特征或功能。
【具体实施方式】
[0026] 在下面的描述中，出于说明的目的，为了提供对一个或多个实施例的全面理解，阐 述了许多具体细节。然而，很明显，也可以在没有这些具体细节的情况下实现这些实施例。
[0027] 以下将结合附图对本发明的具体实施例进行详细描述。
[0028] 为了说明本发明提供的基于卡尔曼滤波偏最小二乘近红外光谱动态演化模型校 正方法，图1示出了根据本发明实施例的基于卡尔曼滤波偏最小二乘近红外光谱动态演化 模型校正方法流程。
[0029] 如图1所示，本发明提供的基于卡尔曼滤波偏最小二乘近红外光谱动态演化模型 校正方法包括：
[0030] S110:利用K/S(Kennard-Stone)算法从标准样品中选择有代表性的建模样品；
[0031] S120:采用PLS法对所述建模样品建立近红外光谱数据与浓度间的线性关系；
[0032] S130:线性关系利用去一交互检验法确定PLS因变量数，形成初始PLS校正模型，然 后通过所述标准样品确定PLS校正模型；
[0033] S140:定期对待测样品进行化验，获取所述待测样品的样品数据；同时，利用所述 PLS校正模型对所述待测样品进行预测，获取所述待测样品的预测值；
[0034] S150:获取的所述待测样品的样品数据和所述待测样品的预测值通过采用KF算法 修正所述PLS校正模型的主因子系数。
[0035] 在上述步骤S120中，设AnXm为η个样品在m个波长上的光谱参数矩阵，Cnxi为η个样 品1种成分含量构成的浓度矩阵。PLS法不直接建立每种成份与光谱参数向量的关系方程， 而是考虑AnXm与CnXl的外部关系和联系二者的内部关系，将AnXm和CnXl分解为如下形式：
[0036] AnXm一TnXhPhXm+EnXm [0037 ] CnXl = UnXhQhXl+FnXl
[0038] 其中，h为样品的抽象组分数，TnXh为光谱参数特征因子阵，UnXh为浓度特征因子 阵，Phx m为光谱参数载荷阵，Qhxi为浓度载荷阵，EnxdPFnxi分别为光谱参数阵和浓度矩阵的 残差矩阵。
[0039] 然后建立TnXt^UnXh的关系矩阵B(浓度与光谱参数间的内部关系）：
[0040] UnXh = TnXhBhXh
[0041 ]此时，浓度与光谱参数间的外部关系为：
[0042] CnXl = TnXhBhXhQhXl+FnXl
[0043] 其中要求I |Fnxi| I达到最小。
[0044] 在上述步骤S130中，在建立PLS初始校正模型中，Tnxh与UnXh的关系矩阵Bhxh，设B矩 阵元素分别为b(i，j)(i = l，2,…，h;j = l，2,…，h)。为使回归关系系数的计算转化为滤波 递推估计形式，将模型中的所有系数值组成状态向量：
[0045] ff=[b(l,l)---b(l,h)---b(h,l)---b(h,h)]T
[0046] 则系统的状态方程和观测方程可表示为：
[wk=w,,
[0047] ^ Iq-c,+^=^,4) + ?
[0048] 其中，Cek为标样浓度，Ak为第k个样品光谱矢量，Crk为预测浓度;Vk为观测噪声，即 为随机白噪声，其统计特性为：
[0049] 印;）= 〇，￡(W ) = /?,' _〇]令盖U=~；啊，為)-;一
[0051 ] 则观测方程为：
[0052] Cek = Hkffk+Dk+Vk
[0053] 在本发明中，需要说明的是，在建立PLS校正模型后，待测样品就可以使用PLS校正 模型进行预测，当时此模型使用一段时间后，可以再次对PLS校正模型进行校正，则校正的 具体方法就是上述步骤S140至S150,校正后的模型在预测时更为精确。
[0054]在本发明中，由卡尔曼滤波基本公式得PLS校正模型主因子系数学习的卡尔曼滤 波算法如下：
[0055] 时间更新(预测）
[0056] Step 1:向前推算权值变量，其中公式为：七
[0057] Step 2:向前推算误差协方差，其中公式为：4二6>:1 [0058]测量更新(修正）
[0059] Step 3:计算卡尔曼增益，其中公式为：6 //，（/?.//，+& )-〖
[0060] Step 4:由期望输出Yek更新估计，其中公式为：#,-=#;+心[(^--/|@--,4)
[0061] Step 5:更新误差协方差，其中公式为心//;)6
[0062] Step 6:k = k+l转Step 2
[0063] 其中，^^为初始估计权值变量;Pk-i为初始估计误差协方差。
[0064] 与上述方法相对应，本发明还提供一种基于卡尔曼滤波偏最小二乘近红外光谱动 态演化模型校正系统，图2示出了根据本发明实施例的基于卡尔曼滤波偏最小二乘近红外 光谱动态演化模型校正系统逻辑结构。
[0065] 如图2所示，本发明提供的基于卡尔曼滤波偏最小二乘近红外光谱动态演化模型 校正系统200,包括:建模样品选取单元210、线性关系建立单元220、PLS校正模型建立单元 230、样品数据获取单元240、预测值获取单元250和PLS校正模型的主因子系数修正单元 260 〇
[0066] 具体地，建模样品选取单元210，用于利用1(/5(1(6111^1(1-51:〇116)算法从标准样品中 选择有代表性的建模样品；
[0067]线性关系建立单元220,用于采用PLS法对所述建模样品建立近红外光谱数据与浓 度间的线性关系；
[0068] PLS校正模型建立单元230,用于利用去一交互检验法确定PLS因变量数，形成初始 PLS校正模型，然后通过所述标准样品确定PLS校正模型；
[0069]样品数据获取单元240,用于定期对待测样品进行化验，获取所述待测样品的样品 数