基于监控指标切换的多工况过程监控方法和系统的制作方法
【专利摘要】本发明公开了一种基于监控指标切换的多工况过程监控方法和系统，包括：采集不同工况下的正常数据作为训练样本集；基于该训练样本集得到隐马尔科夫模型，并获取隐马尔科夫模型的监控指标所对应的控制限；基于各个工况的训练样本分别建立对应工况的统计模式分析模型，并获取各个统计模式分析模型的监控指标所对应的控制限；基于实时获取的过程数据计算工况向量，进而计算差分工况向量；根据差分工况向量的范数，计算相应的实时监控指标，并将其与对应模型的监控指标所对应的控制限比对，来监控工况的运行状况，该方法实时获取过程数据保证监控的可靠性，且不需要每个工况下的数据服从高斯分布，具有更高的适用性。
【专利说明】基于监控指标切换的多工况过程监控方法和系统
【技术领域】
[0001]本发明涉及过程监控领域，尤其涉及一种基于监控指标切换的多工况过程监控方法和系统。
【背景技术】
[0002]对于过程监控和故障诊断问题，传统的方法大多采用多元统计过程控制技术(Multivariable Statistical Process Control, MSPC),其中以主兀分析(PrincipalComponent Analysis, PCA)和偏最小二乘(Partial Least Squares, PLS)为代表的方法已在工业过程监控中得到了成功的应用。传统的MSPC方法均假设过程运行在单一的操作工况下，但是实际上由于产品改变、产能调整等原因过程常在多个工况中频繁的切换。
[0003]基于主元分析和支持向量数据描述等多工况过程监控方法都假设每个工况的数据服从高斯分布，但这在实际中并不一定成立。并且对多个工况建立统一的模型虽然相比建立多个模型的方法简单，然而缺乏实时对工况的辨识，这会导致无法监控当前设备的工作状况。而基于稀疏表示的多工况过程监控方法虽然不假设数据的高斯性，但该方法没有考虑在化工过程中的动态特性。另外，对单个当前数据判断其所属工况受到噪声的影响可能不准确。
【发明内容】
[0004]本发明所要解决的技术问题之一是需要提供一种基于监控指标切换的多工况过程监控方法，其根据实时获取的过程数据是否属于同一工况来对应地判断该过程是否出现故障。此外，还提供了 一种基于监控指标切换的多工况过程监控系统。
[0005]为了解决上述技术问题，本发明提供了一种基于监控指标切换的多工况过程监控方法，包括:采集步骤，采集不同工况下的正常数据作为训练样本集；第一获取步骤，基于所述训练样本集得到隐马尔科夫模型，并获取所述隐马尔科夫模型的监控指标所对应的控制限；第二获取步骤，基于各个工况的训练样本分别建立对应工况的统计模式分析模型，并获取各个统计模式分析模型的监控指标所对应的控制限；计算步骤，基于实时获取的过程数据计算工况向量，并基于所述工况向量计算差分工况向量；监控步骤，根据所述差分工况向量的范数，选取所述隐马尔科夫模型或各个工况对应的统计模式分析模型相应的监控指标，并对所选取的监控指标进行实时计算，将计算得到的实时监控指标与该模型的监控指标所对应的控制限比对，来监控该工况的运行状况。
[0006]在一个实施例中，所述统计模式分析模型的监控指标进一步包括第一监控指标和第二监控指标，所述第一监控指标对应的控制限为第一控制限，所述第二监控指标对应的控制限为第二控制限；所述监控步骤进一步通过以下步骤判断工况过程是否正常:若所述差分工况向量的范数为零，则计算该差分工况向量对应的工况向量的统计模式分析模型的实时第一监控指标和实时第二监控指标，在所述实时第一监控指标大于第一控制限或所述实时第二监控指标大于第二控制限时，则判断出工况过程发生异常；若所述差分工况向量的范数不为零，则计算所述隐马尔科夫模型的实时NLLP监控指标，在该NLLP监控指标大于所述隐马尔科夫模型的监控指标的控制限时，则判断出工况过程发生异常，其中，NLLP表示所述隐马尔科夫模型的负对数似然概率。
[0007]在一个实施例中，在所述监控步骤中，若差分工况向量的范数IlVIlI= 0，则根据如下表达式计算工况q所对应的统计模式分析模型的实时第一监控指标DJs)和实时第二监控指标Dp (S):
[0008]Dr(S) = W CicjjSW2 =S1CicliS
[0009]Dp(s) = StPw [AwF1Pwts, [0010]其中，r为残差子空间，P为主元子空间，s为实时获取的过程数据的统计模式向量，Cw为工况q下建立的统计模式分析模型中残差空间的投影矩阵，P(Q)为工况q下建立的统计模式分析模型中的负载矩阵，Λω为工况q下建立的统计模式分析模型中主元所对应的协方差矩阵特征值拼成的对角矩阵。
[0011]在一个实施例中，通过以下表达式计算过程数据的统计模式向量s:s =(:01[μ，Σ，Ξ]，其中，μ表示实时获取的W时长内的均值向量，Σ表示二阶矩，Ξ表示高阶矩，col [.]表示将矩阵排列成列向量的形式，所述二阶矩包括方差、协方差以及自相关系数。
[0012]在一个实施例中，在所述监控步骤中，若差分工况向量的范数11 ▽ Il I关0，则根据如下公式计算所述隐马尔科夫模型的实时NLLP监控指标:
[0013]NLLP = -1ogPr (Onew | s*),
[0014]其中，PH.)表示事件发生的概率，Onew为实时获取的最新的过程数据，s*表示Onew所处的工况。
[0015]在一个实施例中，在所述计算步骤中，基于实时获取的过程数据使用Viterbi算法获取所述工况向量I = [ii, i2,…，iw]T，其中，w表示过程数据的获取时长，ij(j = I, *.., w)表示w时长内过程数据所处的工况序号。
[0016]在一个实施例中，利用以下表达式来计算所述差分工况向量:
[0017]V I = [ V I1, V i2,…，V IiwU11K V Ij = l-ψ (iJ+1-1j),函数 Ψ (.)在 O处取值为I，其余点取值为O。
[0018]根据本发明的另一方面，还提供了一种基于监控指标切换的多工况过程监控系统，包括:采集模块，其用于采集不同工况下的正常数据作为训练样本集；第一获取模块，其基于所述训练样本集得到隐马尔科夫模型，并获取所述隐马尔科夫模型的监控指标所对应的控制限；第二获取模块，其基于各个工况的训练样本分别建立对应工况的统计模式分析模型，并获取各个统计模式分析模型的监控指标所对应的控制限；计算模块，其基于实时获取的过程数据计算工况向量，并基于所述工况向量计算差分工况向量；监控模块，其根据所述差分工况向量的范数，选取所述隐马尔科夫模型或各个工况对应的统计模式分析模型相应的监控指标，并对所选取的监控指标进行实时计算，将计算得到的实时监控指标与该模型的监控指标所对应的控制限比对，来监控该工况的运行状况。
[0019]在一个实施例中，所述统计模式分析模型的监控指标进一步包括第一监控指标和第二监控指标，所述第一监控指标对应的控制限为第一控制限，所述第二监控指标对应的控制限为第二控制限；在所述监控模块中进一步通过以下步骤判断工况过程是否正常:若所述差分工况向量的范数为零，则计算该差分工况向量对应的工况向量的统计模式分析模型的实时第一监控指标和实时第二监控指标，在所述实时第一监控指标大于第一控制限或所述实时第二监控指标大于第二控制限时，则判断出工况过程发生异常；若所述差分工况向量的范数不为零，则计算所述隐马尔科夫模型的实时NLLP监控指标，在该NLLP监控指标大于所述隐马尔科夫模型的监控指标的控制限时，则判断出工况过程发生异常，其中，NLLP表示所述隐马尔科夫模型的负对数似然概率。
[0020]在一个实施例中，所述计算模炔基于实时获取的过程数据使用Viterbi算法获取所述工况向量I = Li1, i2,…，iw]T，其中，w表示过程数据的获取时长,ij(j = I,..., w)表示w时长内过程数据所处的工况序号。
[0021]与现有技术相比，本发明的一个或多个实施例可以具有如下优点:
[0022]本发明根据不同工况下的所有正常数据与同一工况的正常数据来分别计算不同的控制限，并且实时获取过程数据计算差分工况向量，最后基于差分工况向量的范数选择切换至合适的监控指标，通过对比该监控指标与其对应的控制限来判断该过程是否正常，该方法实时获取过程数据保证了判断的可靠性，并且不需要每个工况下的数据服从高斯分布，具有更高的适用性。
[0023]本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
【专利附图】
【附图说明】
[0024]附图用来提 供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例共同用于解释本发明，并不构成对本发明的限制。在附图中:
[0025]图1是根据本发明一实施例的基于监控指标切换的多工况过程监控方法的流程图；
[0026]图2是根据本发明一示例的实时过程数据采集的示意图；
[0027]图3是根据本发明一实施例的基于监控指标切换的多工况过程监控系统的框图；
[0028]图4是根据本发明一示例的测试例I中有限混合高斯模型方法的检测结果的曲线图；
[0029]图5是根据本发明一示例的测试例I中SPA方法的检测结果的曲线图；
[0030]图6是根据本发明一示例的测试例I中基于监控指标切换的多工况过程监控方法的检测结果的曲线图；
[0031]图7是根据本发明一示例的测试例2中有限混合高斯模型方法的检测结果的曲线图；
[0032]图8是根据本发明一示例的测试例2中SPA方法的检测结果的曲线图；
[0033]图9是根据本发明一示例的测试例2中基于监控指标切换的多工况过程监控方法的检测结果的曲线图；
[0034]图10是根据本发明的另一示例的连续搅拌加热罐的结构示意图；
[0035]图11是根据本发明的另一示例的有限混合高斯模型方法的检测结果的曲线图；[0036]图12是根据本发明的另一示例的SPA方法的检测结果的曲线图；
[0037]图13是根据本发明的另一示例的基于监控指标切换的多工况过程监控方法的检测结果的曲线图。
【具体实施方式】
[0038]为使本发明的目的、技术方案和优点更加清楚，以下结合附图对本发明作进一步地详细说明。
[0039]第一实施例
[0040]图1是根据本发明一实施例的基于监控指标切换的多工况过程监控方法的流程图。下面结合图1对该方法进行详细说明。
[0041]步骤S110，采集不同工况(即动力设备在一定条件下的工作状态)下的正常数据作为训练样本集。
[0042]具体来说，从化工过程数据库中获取不同工况下的正常数据作为训练样本集:
5?,劣，…，其中，5,= 1，…Μ)是第i个工况的数据样本，f
表示N行m列的实数矩阵，Ni表示第i个工况的样本数目，N表示总的样本数目，m表示传感器的个数。
[0043]步骤S120，基于训练样本集得到隐马尔科夫模型，并获取该隐马尔科夫模型的监控指标所对应的控制限；基于各个工况的训练样本分别建立对应工况的统计模式分析模型(以下简称为SPA模型)，并获取各个的SPA模型的监控指标所对应的不同的控制限。
[0044]详细地说，本实 施例获取隐马尔科夫模型的NLLP监控指标的控制限δ NLLP, NLLP表示隐马尔科夫模型的负对数似然概率。具体地，利用训练样本集B，使用EM算法(即最大期望算法)可训练得到隐马尔科夫模型的参数集λ并获得NLLP监控指标的控制限Smp。通常，利用正常训练样本集计算NLLP监控指标，根据选取的置信水平(比如在实际例子中选取的是98%的置信水平)，即可得到控制限δ mp。
[0045]对于SPA模型，利用各个工况的训练样本Bi来建立相应工况下的SPA模型，SPA模型的监控指标进一步包括第一监控指标W (即残差子空间监控指标)和第二监控指标Dp (BP主元子空间监控指标)，SPA模型的监控指标所对应的控制限进一步包括第一控制限和第二控制限δρ。与Smp的获取相似，根据正常训练样本计算&和Dp监控指标，再根据选取的置信水平(比如在实际例子中选取的是98%的置信水平)，即可得到控制限和δρ。
[0046]此外，SPA方法基于中心极限定理，而中心极限定理的思想为“不论随机变量的分布如何，随机变量的统计量渐进地服从高斯分布”，所以本实施例的方法并不需要每个工况下的数据服从高斯分布。
[0047]步骤S130，基于实时获取的过程数据计算工况向量，并基于该工况向量计算差分工况向量。
[0048]在本实施例中，将实时获取的过程数据形象化在一个可滑动的窗口中，即这里计算的过程数据为窗口内的数据，通过沿采样时间轴滑动窗口向前移动可更新窗口内部的过程数据，实时获取数据的时长则为滑动窗口的长度。具体地，基于实时获取的过程数据使用Viterbi算法获取工况向量Ik = [ilk, i2k,…，iwk]T，其中，w表示过程数据的获取时长,ijk(j=l,-,w)表示第k个滑动窗口(即当前的滑动窗口)内数据所处的工况序号。[0049]图2是根据本发明一示例的实时过程数据采集的示意图。容易理解，图2中的两个滑动窗口并不同时存在，而是滑动窗口 I包含采样时间段[10，50]内的过程数据，滑动窗口 II包含采样时间段[70，110]内的过程数据。
[0050]进一步地，差分工况向量的表达式为▽ Ik = [V ilk, V i2k,…，V i(w_1)k]T，其中，▽ iJk = 1-Ψ (i(j+1)k_ijk)，iJk(j = I, *.., w)表示当前的滑动窗口(即滑动窗口 k)内数据所处的工况序号，函数Ψ(.)在O处取值为I，其余点取值为O。
[0051]步骤S140，基于上述差分工况向量的范数，选取隐马尔科夫模型或各个工况对应的SPA模型相应的监控指标，并对所选取的监控指标进行实时计算，将计算得到的实时监控指标与该模型的监控指标所对应的控制限比对，来监控该工况的运行状况。
[0052]首先，判断差分工况向量的范数I I ▽ Il I是否为零。若差分工况向量的范数
IVIl为零，说明当前所获取的过程数据属于同一工况，即工况并没有发生变化，如图2
中滑动窗口 I中所示。若差分工况向量的范数IlVIlI不为零，则说明当前所获取的过程数据属于不同工况，如图2中滑动窗口 II中所示。
[0053]然后，在差分工况向量的范数IlVIlI为零时，计算该差分工况向量对应的工况的实时监控指标(即SPA模型的实时监控指标DJs)和Dp(s))，在实时第一监控指标DJs)大于其所对应的控制限L时或实时第二监控指标Dp(S)大于其所对应的控制限δρ，则判断出工况过程发生异常。通过如下表达式来计算工况q的SPA模型的监控指标:
[0054]Dr(S)^W dll2 - sTCig、s
[0055]Dp(s) = sTP(q) [A (q)]_1P(q)Ts,`[0056]其中，r为残差子空间，P为主元子空间，s为实时获取的过程数据的统计模式向
量C为工况q下建立的SPA模型中残差空间的投影矩阵，P(q)为工况q下建立的SPA模
型中的负载矩阵，Λ ω为工况q下建立的SPA模型中主元所对应的协方差矩阵特征值拼成的对角矩阵。具体地，通过以下表达式计算过程数据的统计模式向量s:
[0057]s = col [ μ , Σ , Ξ ]
[0058]其中，μ表示实时获取的w时长内的均值向量，Σ表示二阶矩，Ξ表示高阶矩，符号col[.]表示将矩阵排列成列向量的形式，上述二阶矩包括方差、协方差以及自相关系数。
[0059]在差分工况向量的范数IlVIlI不为零时，则计算隐马尔科夫模型的实时NLLP监控指标，在该NLLP监控指标大于隐马尔科夫模型的监控指标所对应的控制限δ 时，则判断出工况过程发生异常。具体通过以下表达式获取实时NLLP监控指标:
[0060]NLLP = -1ogPr (Onew | s*),
[0061]其中，PH.)表示事件发生的概率，Onew为实时获取的最新的过程数据，？表示Onew所处的工况。实际上，由于NLLP监控指标的计算依赖于工况，所以不同的工况下的NLLP指标的控制限也有所区别。
[0062]在一个示例中，也可以通过公式NLLP = -1ogPr (Onew | λ )来计算NLLP监控指标，其中λ为隐马尔科夫模型的完整参数集。
[0063]值得注意的是，由于根据差分工况向量的不同，不同监控指标在当前时刻的适合程度不同，本实施例智能地在NLLP监控指标与SPA模型监控指标间进行指标切换，而不是像现有监控方法中的一个监控指标一直用到底。如上所述，当差分工况向量不为零时，此时新进入滑动窗口的工况数据量不足，即无法“装满” 一个窗口，求出的统计量不准确，导致故障的误报率增高，因此选择选择隐马尔科夫模型的NLLP监控指标。当差分工况向量为零时，说明数据处于同于工况，求出的统计量准确，监控更加全面，所以选择SPA模型中的两个指标。
[0064]最后，输出过程监控结果，以方便对出现故障的工况进行维护。
[0065]示例
[0066]为了进一步理解本发明，下面对两个示例进行说明。
[0067]示例一:数值仿真
[0068]用下面的线性系统产生数据:
X1 0.3723 0.6815「] C1
[0069]X2 = 0.4890 0.2954 1 + e2
—x3」[0.9842 0.1793」匕~」
[0070]其中X = [X1 X2 x3]T表示三个传感器的测量，S= [S1 s2]T表示两个独立的数据源，e= Le1 e2 e3]T表示独立的三个均值为O、标准差为0.1的高斯噪声。下面列出系统处于三个不同工况时，数据源所服 从的概率分布:
[0071]工况l:Sl:N(10, 0.8) ;s2:N(12, 1.3)
[0072]工况2: S1: N (5，0.6) ; S2: N (20，0.7)
[0073]工况3:Sl:N(16, 1.5) ;s2:N(30, 2.5)
[0074]其中Ν(μ，ο2)表示均值为μ，方差为σ 2的高斯分布。对三种不同工况的测试数据的描述如表1所示。其中，偏移故障为传感器上的恒定值偏差，漂移故障为传感器上的幅值从小逐渐变化到大的偏差。
[0075]表1.数值仿真中测试数据的描述
[0076]
测试例序号TM数据描述采样时刻^
?MΠ--
测试例I2正常101-200
2幅值为0.5的偏移故障加在X3201-400
3正常1-100测试例2 2 正常 101-2U0
2漂移故障/⑷= 0.0耿-200)加在X1201-400
[0077]为了验证本示例的有效性以及在故障检测率和误报率方面的优势，选取了有限混合高斯模型(FGMM)方法和统计模式分析(SPA)的方法(没有指标切换)作为对比方法。图4-图9显示了各种方法针对测试例I与测试例2的故障检测结果的示意图。其中，图4、图7中的BIP指标为有限闻斯混合I旲型中的贝叶斯推理概率指标，图5、图8中的Dp指标为SPA模型中主元空间的监控指标,?指标为SPA模型中残差空间的监控指标，图6、图9中的Dp-NLLP指标为Dp指标与NLLP指标切换，Dr-NLLP指标为化指标与NLLP指标切换。
[0078]表2所示为各种方法的故障诊断率和误报率的总结。
[0079]表2.示例一中不同方法故障检测率(FDR)和误报率(FAR)的比较
[0080]
【权利要求】
1.一种基于监控指标切换的多工况过程监控方法，包括: 采集步骤，采集不同工况下的正常数据作为训练样本集； 第一获取步骤，基于所述训练样本集得到隐马尔科夫模型，并获取所述隐马尔科夫模型的监控指标所对应的控制限； 第二获取步骤，基于各个工况的训练样本分别建立对应工况的统计模式分析模型，并获取各个统计模式分析模型的监控指标所对应的控制限； 计算步骤，基于实时获取的过程数据计算工况向量，并基于所述工况向量计算差分工况向量； 监控步骤，根据所述差分工况向量的范数，选取所述隐马尔科夫模型或各个工况对应的统计模式分析模型相应的监控指标，并对所选取的监控指标进行实时计算，将计算得到的实时监控指标与该模型的监控指标所对应的控制限比对，来监控该工况的运行状况。
2.根据权利要求1所述的监控方法，其特征在于， 所述统计模式分析模型的监控指标进一步包括第一监控指标和第二监控指标，所述第一监控指标对应的控制限为第一控制限，所述第二监控指标对应的控制限为第二控制限；所述监控步骤进一步通过以下步骤判断工况过程是否正常: 若所述差分工况向量的范数为零，则计算该差分工况向量对应的工况向量的统计模式分析模型的实时第一监控指标和实时第二监控指标，在所述实时第一监控指标大于第一控制限或所述实时第二监控指标大于第二控制限时，则判断出工况过程发生异常； 若所述差分工况向量的范数不为零，则计算所述隐马尔科夫模型的实时NLLP监控指标，在该NLLP监控指标大于`所述隐马尔科夫模型的监控指标的控制限时，则判断出工况过程发生异常，其中，NLLP表示所述隐马尔科夫模型的负对数似然概率。
3.根据权利要求2所述的监控方法，其特征在于， 在所述监控步骤中，若差分工况向量的范数IlViII = O，则根据如下表达式计算工况q所对应的统计模式分析模型的实时第一监控指标DJs)和实时第二监控指标Dp(S):Dr(^) = II C(?)5ll2=/C(?,5Dp(S) = StPw [Λ wF1Pwts, 其中，r为残差子空间，ρ为主元子空间，s为实时获取的过程数据的统计模式向量，Cw为工况q下建立的统计模式分析模型中残差空间的投影矩阵，Pw为工况q下建立的统计模式分析模型中的负载矩阵，Λω为工况q下建立的统计模式分析模型中主元所对应的协方差矩阵特征值拼成的对角矩阵。
4.根据权利要求3所述的监控方法，其特征在于，通过以下表达式计算过程数据的统计模式向量s:
s = col [ μ , Σ , Ξ ], 其中，μ表示实时获取的w时长内的均值向量，Σ表示二阶矩，Ξ表示高阶矩，col [.]表示将矩阵排列成列向量的形式，所述二阶矩包括方差、协方差以及自相关系数。
5.根据权利要求2所述的监控方法，其特征在于， 在所述监控步骤中，若差分工况向量的范数11 VIl I关0，则根据如下公式计算所述隐马尔科夫模型的实时NLLP监控指标:NLLP = -1ogPr (Onew | s*), 其中，PH.)表示事件发生的概率，Onew为实时获取的最新的过程数据，s*表示Onrat所处的工况。
6.根据权利要求1所述的监控方法，其特征在于，在所述计算步骤中，基于实时获取的过程数据使用Viterbi算法获取所述工况向量I = [i1； i2，…，iw]T，其中，w表示过程数据的获取时长，Ij U = 1，...，《)表示w时长内过程数据所处的工况序号。
7.根据权利要求6所述的监控方法，其特征在于，利用以下表达式来计算所述差分工况向量: VI = [ V I1, V I2,…，V iw_JT,其中，▽ ij = l-ψ (iJ+1-1j)，函数 Ψ (.)在 O 处取值为I，其余点取值为O。
8.一种基于监控指标切换的多工况过程监控系统，包括: 采集模块，其用于采集不同工况下的正常数据作为训练样本集； 第一获取模块，其基于所述训练样本集得到隐马尔科夫模型，并获取所述隐马尔科夫模型的监控指标所对应的控制限； 第二获取模块，其基于各个工况的训练样本分别建立对应工况的统计模式分析模型，并获取各个统计模式分析模型的监控指标所对应的控制限； 计算模块，其基于实时获取的过程数据计算工况向量，并基于所述工况向量计算差分工况向量； 监控模块，其根据所述差分工况向量的范数，选取所述隐马尔科夫模型或各个工况对应的统计模式分析模型相应的监控指标，并对所选取的监控指标进行实时计算，将计算得到的实时监控指标与该模型的监控指`标所对应的控制限比对，来监控该工况的运行状况。
9.根据权利要求8所述的监控系统，其特征在于，所述统计模式分析模型的监控指标进一步包括第一监控指标和第二监控指标，所述第一监控指标对应的控制限为第一控制限，所述第二监控指标对应的控制限为第二控制限； 在所述监控模块中进一步通过以下步骤判断工况过程是否正常: 若所述差分工况向量的范数为零，则计算该差分工况向量对应的工况向量的统计模式分析模型的实时第一监控指标和实时第二监控指标，在所述实时第一监控指标大于第一控制限或所述实时第二监控指标大于第二控制限时，则判断出工况过程发生异常； 若所述差分工况向量的范数不为零，则计算所述隐马尔科夫模型的实时NLLP监控指标，在该NLLP监控指标大于所述隐马尔科夫模型的监控指标的控制限时，则判断出工况过程发生异常，其中，NLLP表示所述隐马尔科夫模型的负对数似然概率。
10.根据权利要求8所述的监控系统，其特征在于， 所述计算模炔基于实时获取的过程数据使用Viterbi算法获取所述工况向量I =[ii, i2,…，iw]T，其中，w表示过程数据的获取时长，ij(j = I,…，w)表示w时长内过程数据所处的工况序号； 在获取所述工况向量后，所述计算模块进一步利用以下表达式计算所述差分工况向量:▽工=[▽ i” ▽ i2，…，V iw_jT，其中，▽ ij = nad.)，函数 Ψ(.)在O处取值为I，其余点取值为O。
【文档编号】G05B17/00GK103631145SQ201310675045
【公开日】2014年3月12日 申请日期:2013年12月11日 优先权日:2013年12月11日
【发明者】周东华, 宁超, 陈茂银 申请人:清华大学