发明名称 一种基于改进粒子群优化算法的电力系统故障诊断方法
摘要 本发明公开了一种基于改进粒子群优化算法的电力系统故障诊断方法,所述粒子群优化算法初始化为一群随机粒子,然后通过迭代找到最优解,在每一次迭代中,粒子通过跟踪两个“极值”来更新自己:第一个极值就是粒子本身所找到的最优解,这个极值称为个体极值pbest;另一个极值是整个种群目前找到的最优解,这个极值是全局极值gbest。改进的粒子群优化算法在迭代过程引入进化速度因子H,寻求极小值时,令H为上一次迭代的全局最优值/当前迭代的全局最优值,计算粒子的进化速度,当H低于Δ时,在速度更新方程中加入负扰动项以改变粒子运动方向。与遗传算法方法相比,粒子群优化算法具有稳定性高、收敛特性好、运行速度快的突出优点。
申请公布号 CN104820142A 申请公布日期 2015.08.05
申请号 CN201510126409.2 申请日期 2015.03.20
申请人 国家电网公司;国网河北省电力公司;国网河北省电力公司鹿泉市供电分公司 发明人 薛飞;王珂;石立志
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 石家庄元汇专利代理事务所(特殊普通合伙) 13115 代理人 刘闻铎
主权项 一种基于改进粒子群优化算法的电力系统故障诊断方法,其特征在于,所述粒子群优化算法的电力系统故障诊断步骤如下:①初始化粒子,粒子个数,pbest,gbest,迭代次数,w,C<sub>1</sub>,C<sub>2</sub>;②若未达到最大迭代次数,则根据(2)式计算每个粒子的f(S)适应值,更新pbest和gbest;若达到最大迭代次数,则跳转到步骤7执行;x<sub>i+1</sub>=x<sub>i</sub>+V<sub>i+1</sub>     (2)其中,下标i表示迭代次数,x<sub>1</sub>表示第i次迭代时的粒子空间位置,V<sub>1</sub>表示第i次迭代时的粒子速度。③为防止陷入局部最优,若某个粒子的位置等于最佳位置,则对该粒子重新初始化赋值;④根据式(4)计算H,若H小于Δ时在速度更新方程中加入负扰动项以改变粒子运动方向,即用式(4)和(5)产生新的速度,并更新粒子群的位置;若H大于或等于Δ,则用式(3)和(4)产生新的速度,并更新粒子群的位置;H=f(t‑1)/f(t)        (3)式中:H为速度因子。V<sub>i+1</sub>=W*V<sub>i</sub>+C<sub>1</sub>*rand()*(pbest‑x<sub>i</sub>)+C<sub>2</sub>*rand()*(gbest‑x<sub>i</sub>)‑C<sub>d</sub>*(pbest+gbest)/2      (4)式中cd为扰动因子,可在0.30~0.65之间取值;(pbest+gbest)为群体和个体最优解的平均值;惯性权重w采取线性递减策略。<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>S</mi><mo>)</mo></mrow><mo>=</mo><mi>U</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mi>r</mi></msub></munderover><mo>|</mo><msub><mi>r</mi><mi>k</mi></msub><mo>-</mo><msup><msub><mi>r</mi><mi>k</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>S</mi><mo>)</mo></mrow><mo>|</mo><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mi>c</mi></msub></munderover><mo>|</mo><msub><mi>c</mi><mi>j</mi></msub><mo>-</mo><msup><msub><mi>c</mi><mi>j</mi></msub><mo>*</mo></msup><mrow><mo>(</mo><mi>S</mi><mo>,</mo><mi>R</mi><mo>)</mo></mrow><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000685751200000011.GIF" wi="1159" he="217" /></maths>⑤由于粒子的位置只有0、1两种状态,因此对函数值以0.5为分界点进行离散化,x(t)&lt;0.5取值为0,否则取为1;⑥迭代次数加1,返回步骤2;⑦输出结果。
地址 100000 北京市西城区西长安街86号