1.一种基于非线性动态成分分析的动态过程监测方法，其特征在于，包括以下步骤：n×m

步骤(1)采集生产过程正常运行状态下的样本，组成训练数据矩阵X∈R ，并计算矩阵X中各列向量的均值μ1，μ2，…，μm以及标准差δ1，δ2，…，δm，对应组成均值向量μ＝[μ1，μ2，…，T

μm] 与标准差向量δ＝[δ1，δ2，…，δm]，其中，n为训练样本数，m为过程测量变量数，R为实数n×m

集，R 表示n×m维的实数矩阵，上标号T表示矩阵或向量的转置；

步骤(2)根据如下所示公式对矩阵X实施标准化处理得到n×m T

上式①中，U∈R 是由n个相同的均值向量μ组成的矩阵，即U＝[μ，μ，…，μ] ，对角矩阵Φ中对角线上的元素由标准差向量δ组成；

步骤(3)记 设置自相关阶数为D后，依据如下所示公式计算核矩阵Kn×n

∈R ：

其中，xi为第i个标准化后的样本数据，上标号T表示矩阵或者向量的转置，Ki，j表示核矩阵K中的第i行第j列元素，c为核函数参数，下标号i＝1，2，…，n，j＝1，2，…，n；

步骤(4)根据公式 中心化处理核矩阵K从而得到中心化后的核矩阵 其中，矩阵 中所有元素都等于1；

步骤(5)利用非线性动态成分分析算法求解得到 个非线性动态成分向量 以及相应的特征向量 具体的实施过程如下所示：T D×1

步骤(5.1)初始化τ＝1、矩阵 以及向量β＝[1，1，…，1] ∈R ，并定义操作表示将矩阵 中第a行至第b行、第g列至第h列的元素作为一个矩阵的操作；

步骤(5.2)设置 后，根据公式 分别构造矩T

阵k1，k2，…，kD，并设置K2＝[k1，k2，…，kD]，其中d＝1，2，…，D；

步骤(5.3)根据公式β＝β/||β||单位化处理向量β后，设置Kφ＝β1k1+β2k2+…+βDkD；

步骤(5.4)并求解广义特征值问题： 中最大特征值所对应的特征向量pτ，需保证向量pτ的长度满足条件：步骤(5.4)根据公式 更新向量β，其中 表示Kronecker内积，其具体的计算方式如下所示：步骤(5.5)判断向量β是否收敛；若否，则返回步骤(5.3)；若是，则得到第τ个特征向量pτ，并执行步骤(5.6)步骤(5.6)根据公式 计算第τ个非线性动态成分向量tτ，并根据公式 更T T

新核矩阵 其中矩阵θτ＝In‑tτtτ/(tτtτ)，In为n×n维的单位矩阵；

步骤(5.7)判断是否满足条件： 若是，则置τ＝τ+1后返回步骤(5.2)；若否，则输出矩阵 个非线性动态成分向量 以及相应的特征向量步骤(6)根据公式 计算残差中各样本的监测统计量Q，其中矩阵 diag{}表示将矩阵对角线上的元素转变成向量的操作；

步骤(7)依据如下所示公式分别计算监测统计量的控制上限ψlim与Qlim：上两式中， 表示置信度为α、自由度分别为 与 的F分布所对应的值， 表示自由度为ξ、置信度为α为卡方分布所对应的值，u和v分别为步骤(6)中监测统计量Q的均值和方差；

1×m

步骤(8)收集新采样时刻的数据样本x∈R ，并根据公式 对x实施标准化处理得到

1×n

步骤(9)根据如下所示公式计算核向量κ∈R ：其中，κi表示核向量κ中的第i个元素，且i＝1，2，…，n；

步骤(10)根据公式 对κ实施中心化处理得到中心化后的核向量

1×n

其中向量Ix∈R 中所有元素都等于1；

步骤(11)根据步骤(5)中的输出结果计算对应于当前样本数据的非线性动态成分具体的实施过程如下所示：

步骤(11.1)初始化τ＝1与初始化步骤(11.2)若τ＝1，则依据公式 计算s中第一个非线性动态成分；若τ＞1，则依据公式 更矩阵 后，在根据公式 计算第τ个非线性动态成分；

步骤(11.3)判断是否满足条件： 若是，则置τ＝τ+1后返回步骤(11.2)；若否则将得到的 组成向量

步骤(12)根据如下所示公式计算监测统计量ψ与Q上式⑦中，矩阵 矩阵

步骤(13)判断是否满足条件ψ≤ψlim且Q≤Qlim；若是，则当前样本数据采集自正常工况，返回步骤(8)继续监测下一时刻的样本数据；若否，则当前监测样本数据采集自故障工况。