[发明专利]基于平均影响值数据变换的特征加权K近邻故障诊断方法

权利要求书

1.基于MIV数据变换的特征加权电网信息故障诊断方法，其特征在于该方法具体包括以下各步骤：

步骤(1)依据智能电网信息物理系统信号的前向传播和误差的反向传播来构建BP神经网络：

步骤(1-1)网络的初始化：

假设输入层的节点个数为n，隐含层的节点个数为l，输出层的节点个数为m；输入层到隐含层的权重ω_ij，隐含层到输出层的权重为ω_jk'，输入层到隐含层的偏置为a_j，隐含层到输出层的偏置为b_k'；以上权重及偏置参数的初始化是随机的，在接下来的多次迭代中进行一次次的更新修正；学习速率为η，取[0.01,0.1]之间，激励函数g(x)的输出结果是[-1,1]或者[0,1]，其中激励函数g(x)将其取作Sigmoid函数；形式为：

步骤(1-2)正向传递：有了输入和初始参数以及激励函数，一个节点的输出表示如下，对应输入乘以连接隐含层的对应权值加上偏置，最后经过激励函数得到隐含层节点的一个输出，假设系统输入为X＝[x₁,x₂,....x_n]，计算表达式如下；

隐含层节点的输入：

隐含层节点的输出：

输出节点的输出：

其中， i＝1,2,…,n；j＝1,2,…,l；k'＝1,2,…,m；

步骤(1-3)反向传递：误差计算，定义损失函数如下：

以上，Y为期望输出，O为实际输出，权值更新即是隐含层和输出层权值更新；其表达式如下

隐含层到输出层的权值更新计算：

ω_jk'＝ω_jk'+ηH_je_k' (8)

输出层到隐含层的权值更新计算：

步骤(1-4)判断算法迭代是否结束，一种是迭代次数达到上限，另一种是判断相邻的两次误差之间的差别是否小于指定的值；

步骤(2)将待输入的智能电网信息物理系统样本数据通过公式(13)进行标准化处理

步骤(3)求取数据集的各输入变量对输出的MIV，MIV是用于确定输入神经元对输出神经元影响大小的一个指标，其符号代表相关的方向，绝对值大小代表影响的相对重要性；具体计算过程如下：

IV_i＝O_i,+-O_i,- (17)

综上，取0.1≤δ≤0.3；表示所有样本的第i个特征分量分别增减δ后的矩阵；表示第k'个样本数据的第i个特征分量分别增减δ后的向量；x_i(k')表示第k'样本数据的第i个特征分量的值；网络训练终止后，将训练样本集和分别作为仿真样本利用已建成的网络进行仿真，得到两个仿真结果O_i,+和O_i,-，IV_i即为样本数据中第i个变量的平均影响值；同理可求得样本数据中其它各变量的平均影响值；

IV＝[IV₁,IV₂,…,IV_n]^T (18)

IV中的各个参数即为样本数据中各个特征变量对输出的影响值；最后，多次仿真取m组IV值，求取其中对应特征变量IV值的平均值，记作MIV；

MIV绝对值的大小即为各自变量对网络输出影响的相对重要性，从而实现样本数据特征的加权；

步骤(4)K近邻算法设计

a.训练样本中找k个近邻

计算测试样本中每一个样本x(k')，k'＝1,2,…,N，与训练样本之间的欧氏距离，并找到距离它最近的k'个样本；

d_i,j＝||x(i)-x(j)||₂,i≠j (20)

上式中，x(i)和x(j)分别表示两个样本，i,j＝1,2,…,N；

b.统计计算k个近邻中属于各类的概率

以上k表示k个近邻，k_α表示k个近邻中属于α类的样本个数，α＝1，2，…，c，其中c为样本中类别个数；

c.故障诊断

对于一个新的样本x(k)，根据式(20)计算其k个近邻；根据式(21)统计计算p(α)；将p(α)从大到小排序，p_max(α)所属的类就是测试样本x(k)所属的故障类。