[0005] 本发明所要解决的主要技术问题是：实施数据驱动的风力发电机故障检测时，能实时提取出有利于分离故障的单变量特征，从而可以将监测多个测量数据转变成单变量的异常变化监测思路，以实现实时检测风力发电机是否出现故障的目的。具体来讲，本发明方法通过在线即时特征提取，通过最优化的思想即时提取最有利于分离故障的单变量特征，并通过上下限范围判断来完成风力发电机的故障检测任务。
[0006] 本发明方法解决上述问题所采用的技术方案为：一种基于单变量特征提取策略的风力发电机故障检测方法，包括以下所示步骤：
[0007] 步骤(1)：在风力发电机正常运行状态下，利用数据采集系统每间隔1分钟采集一次样本数据，具体包括11个数据，并将每次采集的这11个数据组成一个列向量；其中，列向量中的11个数据依次为：风速，旋翼转速，发电机转速，机械转矩，发电功率，叶片螺距角，叶片方位角，叶片根部力矩，顶部水平轴加速度，顶部纵轴加速度，偏航误差。
[0008] 步骤(2)：按照采集的时间先后顺序，将风速在5m/s至20m/s之间的N个列向量x1，11×N
x2，...，xN组成矩阵X＝[x1，x2，...，xN]，并对X∈R 中各个行向量实施标准化处理从而得
11×N 11×1
到新矩阵 其中，R 表示11×N维的实数矩阵，R 表示11×1维的实数向量，标
准化处理的方式具体如步骤(2.1)至步骤(2.2)所示。
[0009] 步骤(2.1)：设zj∈R1×N表示矩阵X中的第j行的行向量；其中，j∈{1，2，...，11}。
[0010] 步骤(2.2)：计算行向量zj中所有元素的平均值μj和标准差δj后，根据计算得到新矩阵 中第j行的行向量
[0011] 步骤(3)：根据公式 计算基矩阵C；其中，上标号T表示矩阵或向量的转置。
[0012] 步骤(4)：确定单变量特征的上限ζmax和下限ζmin，具体的实施过程如步骤(4.1)至步骤(4.5)所示。
[0013] 步骤(4.1)：设 表示新矩阵 中第i列的列向量；其中，i∈{1，2，...，N}。
[0014] 步骤(4.2)：根据公式 计算矩阵Gi后，再求解Gi最大特征值所对应的特11×1
征向量pi∈R 。
[0015] 步骤(4.3)：根据公式wi＝Cpi计算投影变换向量wi∈R11×1后，再根据 计算单变量特征yi。
[0016] 步骤(4.4)：分别设置i等于1，2，...，N，并执行步骤(4.1)至步骤(4.3)从而对应得到单变量特征y1，y2，...，yN，并计算y1，y2，...，yN的平均值μ0和标准差δ0。
[0017] 步骤(4.5)：根据ζmax＝μ0+3×δ0和ζmin＝μ0‑3×δ0分别确定单变量特征的上限ζmax和下限ζmin。
[0018] 步骤(5)：利用风力发电机的数据采集系统采集最新采样时刻的样本数据，具体包11×1
括步骤(1)中所列的11个数据，并将这11个数据组成列向量xnew∈R 。
[0019] 步骤(6)：判断xnew中的第一个数据(即：风速)是否小于等于20且大于等于5；若是，则执行步骤(7)；若否，则风力发电机未产生电能，处于待机状态，并返回步骤(4)继续实施对最新采样时刻的风力发电机故障检测。
[0020] 步骤(7)：根据公式 对xnew中各行的元素实施标准化处理，从而得到新数据向量 其中， 和 分别表示xnew和 中第j行的元素。
[0021] 步骤(8)：根据公式 计算矩阵Gnew后，再求解Gnew最大特征值所对应11×1
的特征向量pnew∈R 。
[0022] 步骤(9)：根据公式wnew＝Cpnew计算投影变换向量wnew∈R11×1后，再计算单变量特征[0023] 步骤(10)：判断是否满足条件：ζmin≤ynew≤ζmax；若是，则风力发电机正常运行，并返回步骤(5)继续利用最新采样时刻的样本数据实施对风力发电机的故障检测；若否，则执行步骤(11)。
[0024] 步骤(11)：返回步骤(5)继续利用最新采样时刻的样本数据实施故障检测；若连续6个采样时刻的样本数据所对应的单变量特征皆不满足步骤(10)中的判断条件，则触发故障警报；反之，则风力发电机正常运行，并返回步骤5)继续利用最新采样时刻的样本数据实施对风力发电机的故障检测。
[0025] 值得指出的是，上述步骤(4.2)与步骤(4.3)以及步骤(8)与步骤(9)中求矩阵及其最大特征来源于如下所示算法的推理过程，该算法旨在搜寻投影变换向量，从而使新数据与训练数据进行最大化区分，该算法的的目标函数及其约束条件如下所示：
[0026]
[0027] 上式中，正常工况数据经投影变换向量w变换后，方差或长度为1。而目标函数旨在使在线采样的数据向量经同样的投影变换后，距离远点越远越好。换句话讲，就是尽可能的将在线数据与正常工况数据拉开，这样可以提取到最适合监测故障的特征成分。
[0028] 上式①的优化求解可通过经典的拉格朗日乘子法完成，需要先引入中间量p＝T 1/2(XX) w进行过渡，并且注意 其中tr( )表示求括号内矩阵的迹。
[0029] 通过以上所述实施步骤，本发明方法的优势介绍如下。
[0030] 本发明方法针对不同时刻的样本数据实时计算出对应的投影变换向量，从而可以提取出能最大区分出与正常数据间差异的单变量特征。因此，在风力发电机正常工况数据充裕的情况下，本发明方法能及时的根据采样数据实时检测出风力发电机的故障。在即将陈述的具体实施案例中，通过实际应用验证了本发明方法的可行性与优势。
