动态演化模型校正方法及系统的制作方法
【技术领域】
[0001 ]本发明涉及在线近红外光谱分析技术领域，更为具体地，涉及一种基于卡尔曼滤 波BP神经网络(KF-BPNN)近红外光谱动态演化模型校正方法及系统。
【背景技术】
[0002] 建立准确性高、自适应性强的校正模型是近红外光谱定量分析方法成功应用的关 键。现代近红外光谱以其分析速度快、重现性好、成本低、不消耗样品、易于实现在线分析等 鲜明的特点正得到越来越多的应用。在近红外光谱分析实际应用中，样品近红外光谱和待 测量存在复杂的间接关系。常通过求解标准曲线，建立测量数据与分析体系中某一待测组 分含量之间的对应关系。BPNN具有出色的非线性关系表达能力，在处理近红外光谱建模问 题上显示出独特的优越性，被广泛用于非线性模型校正过程。
[0003] 但是上诉方法通常假设仪器的环境噪声和部件状态是稳定的，是一种静态建模方 法，其效果仅限于短期相对稳定的模型使用过程。然而光谱获取过程受被测对象变化、仪器 老化、环境扰动等影响具有时变性。当仪器测量过程慢慢发生演变时，基于早期数据的校正 模型预测效果将无法得到保证。传统方法通过利用补充样品数据库定期维护模型，从而保 证其适应性，其特点是工作量大、成本高和不易实现在线分析等。且更新的校正模型中包含 大量光谱仪早期状态信息，不能精确代表仪器最新测量状态。
[0004] 为解决上述问题，本发明需要提供一种新的模型校正方法。
【发明内容】
[0005] 鉴于上述问题，本发明的目的是提供一种基于卡尔曼滤波BP神经网络近红外光谱 动态演化模型校正方法及系统，以保证近红外光谱校正模型具有自适应性，降低重建模型 成本，从而实现在线分析。
[0006] 本发明提供一种基于卡尔曼滤波BP神经网络近红外光谱动态演化模型校正方法， 包括：
[0007] 利用K/S算法从标准样品中选择有代表性的建模样品；
[0008] 采用BPNN法对所述建模样品建立近红外光谱数据与浓度间的非线性关系，形成 BPNN校正模型；
[0009] 定期对待测样品进行化验，获取所述待测样品的样品数据；
[0010]同时，利用所述BPNN校正模型对所述待测样品进行预测，获取所述待测样品的预 测值；
[0011] 获取的所述待测样品的样品数据和所述待测样品的预测值通过采用KF算法修正 所述BPNN校正模型的阈值和权值。
[0012] 本发明还提供一种基于卡尔曼滤波偏最小二乘近红外光谱动态演化模型校正系 统，包括：
[0013] 建模样品选取单元，用于利用K/S算法从标准样品中选择有代表性的建模样品；
[0014] BPNN校正模型建立单元，用于采用BPNN法对所述建模样品建立近红外光谱数据与 浓度间的非线性关系，形成BPNN校正模型；
[0015] 样品数据获取单元，用于定期对待测样品进行化验，获取所述待测样品的样品数 据；
[0016] 预测值获取单元，利用所述BPNN校正模型对所述待测样品进行预测，获取所述待 测样品的预测值；
[0017] BPNN校正模型的阈值和权值修正单元，用于对获取的所述待测样品的样品数据和 所述待测样品的预测值通过采用KF算法修正所述BPNN校正模型的阈值和权值。
[0018] 从上面的技术方案可知，本发明提供的基于卡尔曼滤波BP神经网络近红外光谱动 态演化模型校正方法及系统，利用卡尔曼滤波BP神经网络以保证近红外光谱校正模型具有 自适应性，降低重建模型成本，从而达到实现在线分析的目的。
[0019] 为了实现上述以及相关目的，本发明的一个或多个方面包括后面将详细说明并在 权利要求中特别指出的特征。下面的说明以及附图详细说明了本发明的某些示例性方面。 然而，这些方面指示的仅仅是可使用本发明的原理的各种方式中的一些方式。此外，本发明 旨在包括所有这些方面以及它们的等同物。
【附图说明】
[0020] 通过参考以下结合附图的说明及权利要求书的内容，并且随着对本发明的更全面 理解，本发明的其它目的及结果将更加明白及易于理解。在附图中：
[0021] 图1为根据本发明实施例的基于卡尔曼滤波BP神经网络近红外光谱动态演化模型 校正方法流程示意图；
[0022] 图2为根据本发明实施例的基于卡尔曼滤波BP神经网络近红外光谱动态演化模型 校正系统结构框图。
[0023 ]在所有附图中相同的标号指示相似或相应的特征或功能。
【具体实施方式】
[0024]在下面的描述中，出于说明的目的，为了提供对一个或多个实施例的全面理解，阐 述了许多具体细节。然而，很明显，也可以在没有这些具体细节的情况下实现这些实施例。 [0025]以下将结合附图对本发明的具体实施例进行详细描述。
[0026]为了说明本发明提供的基于卡尔曼滤波BP神经网络近红外光谱动态演化模型校 正方法，图1示出了根据本发明实施例的基于卡尔曼滤波BP神经网络近红外光谱动态演化 模型校正方法流程。
[0027]如图1所示，本发明提供的基于卡尔曼滤波BP神经网络近红外光谱动态演化模型 校正方法包括：
[0028] S110:利用K/S(Kennard-Stone)算法从标准样品中选择有代表性的建模样品；
[0029] S120:采用BPNN法对所述建模样品建立近红外光谱数据与浓度间的非线性，形成 BPNN校正模型；
[0030] S130:定期对待测样品进行化验，获取所述待测样品的样品数据；同时，利用所述 BPNN校正模型对所述待测样品进行预测，获取所述待测样品的预测值；
[0031] S140:获取的所述待测样品的样品数据和所述待测样品的预测值通过采用KF算法 修正所述BPNN校正模型的阈值和权值。
[0032]在上述步骤S120中，（1)设Anxm为η个建模样品在m个波长上的光谱矩阵，CnXi为η个 建模样品1种成分含量构成的浓度矩阵。
[0033] Anxm=[ai，a2，···，Βη；Γ(η=1，2,…，Ν)为输入矢量(光谱数据），或称训练样本。
[0035] 其中，WMI(n)为第η次迭代时输入层与隐层I之间的权值矢量。
[0036] Wjp (η)为第η次迭代时隐层J与输出层之间的权值矢量。
[0037] Yk(n) = [yki(n)，yk2(n)，···，ykP(n)](k=l，2，···，Ν)为第η次迭代时网络的实际输 出;CnXl = [C1，C2，···，CnilYn=!，2,…，Ν)为期望浓度。
[0038] (2)初始化，赋给Wmi (0 )，Wjp (0)各一个较小的随机非零值。
[0039] (3)随机输入样本an，n = 0。
[0040] (4)对输入样本an，前向计算BP网络每层神经元的输入信号和输出信号。
[0041] (5)由期望浓度(神经网络输出）cn和上一步求得的实际输出a/，计算误差E(n)，判 断其是否满足要求;若满足要求转至步骤(8);若不满足要求转至步骤(6)。
[0042] (6)判断n+1是否大于最大迭代次数，若最大迭代次数转至步骤(8)，若不大于最大 迭代次数，对输入样本a n，反向计算每层神经元的局部梯度δ。
[0043] (7)按下式计算权值修正量Δ W，并修正权值;η = n+1，转至步骤(4)
[0045]其中，η为学习速率，η为迭代次数。
[0046] (8)判断是否学完所有的训练样本，若是则结束，若否则转至步骤(3)。
[0047]其中，，设有一N层前向网络，每层神经元数为Sk(k= 1，2，...，Ν)，记输入层为第一 层，输出层为第N层；第k层神经元的连接权值Κ(/ = 1，2，.··，Α ^ = 网络中的所 有权值组成状态向量：
[0051]系统的状态方程和观测方程表示为
[0053] 其中，Vk为观测噪声，其统计特性为：
[0054] E(^) = O,E(^yZ) = Ri
[0056] 则观测方程为：
[0057] Cek = Hkffk+Dk+Vk
[0058] 其中，Cek为标样浓度，ak为第k个样品光谱矢量，Crk为预测浓度。
[0059]在本发明中，需要说明的是，在建立BPNN校正模型后，待测样品就可以使用BPNN校 正模型进行预测，当时此模型使用一段时间后，可以再次对BPNN校正模型进行校正，则校正 的具体方法就是上述步骤S130,校正后的模型在预测时更为精确。
[0060] 在本发明中，由卡尔曼滤波基本公式得BPNN校正模型阈值、权值学习的卡尔曼滤 波算法如下：
[0061] 时间更新(预测）
[0062] Step 1:向前推算权值变量，其中公式为：#A_-
[0063] Step 2:向前推算误差协方差，其中公式为：尽=巧_:1 [0064] 测量更新(修正）
[0065] Step 3:计算卡尔曼增益，其中公式为：《二/〗///' /// + Af1
[0066] Step 4:由期望输出Cek更新估计，其中公式为：虼& [CV-4)
[0067] Step 5:更新误差协方差，其中公式为：巧:
[0068] Step 6:k = k+l转Step 2;
[0069] 其中，斤4为初始估计权值变量;Pk-i为初始估计误差协方差。
[0070] 与上述方法相对应，本发明还提供一种基于卡尔曼滤波BP神经网络近红外光谱动 态演化模型校正系统，图2示出了根据本发明实施例的基于卡尔曼滤波BP神经网络近红外 光谱动态演化模型校正系统逻辑结构。
[0071] 如图2所示，本发明提供的基于卡尔曼滤波BP神经网络近红外光谱动态演化模型 校正系统200,包括:建模样品选取单元210、BPNN校正模型建立单元220、样品数据获取单元 230、预测值获取单元240和BPNN校正模型的阈值和权值修正单元250。
[0072] 具体地，建模样品选取单元210,用于利用K/S (Kennard-Stone)算法从标准样品中 选择有代表性的建模样品；
[0073] BPNN校正模型建立单元220，用于采用BPNN法对所述建模样品建立近红外光谱数 据与浓度间的非线性关系，形成BPNN校正模型；
[0074]样品数据获取单元230,用于定期对待测样品进行化验，获取所述待测样品的样品 数据；
[0075]预测值获取单元240,用于利用所述BPNN校正模型对所述待测样品进行预测，获取 所述待测样品的预测值；
[0076] BPNN校正模型的阈值和权值修正单元250，用于对获取的所述待测样品的样品数 据和所述待测样品的预测值通过采用KF算法修正所述BPNN校正模型的阈值和权值。
[0077] 其中，BPNN校正模型建立单元220在采用BPNN法对所述建模样品建立近红外光谱 数据与浓度间的非线性关系，形成BPNN校正模型的过程中，
[0078] (1)设AnXm为η个建模样品在m个波长上的光谱矩阵，CnXi为η个建模样品1种成分含 量构成的浓度矩阵；
[0081 ]其中，WMi(n)为第η次迭代时输入层与隐层I之间的权值矢量；
[0082] Wjp (η)为第η次迭代时隐层J与输出层之间的权值矢量；
[0083] Yk(n) = [yki(n)，yk2(n)，···，ykP(n)](k=l，2，···，Ν)为第η次迭代时网络的实际输 出；