发明名称 一种基于BUD谱峭度的暂态电能质量扰动分类识别方法
摘要 本发明为一种基于BUD谱峭度的暂态电能质量扰动分类识别方法。该方法该算法采用BUD谱峭度方法计算暂态脉冲和暂态振荡两类扰动信号的谱峭度,选取谱峭度的最大值、最小值和均值作为特征量,输入PSO优化参数的SVM进行训练测试。通过PSCAD/EMTDC获得仿真数据,用该方法进行分析。结果表明,基于BUD谱峭度方法能够有效提取扰动特征量,且具有良好的抗噪性能;利用SVM分类器在小样本和叠加有其他扰动时,有效识别两类扰动,识别率较高。
申请公布号 CN102750543B 申请公布日期 2015.10.14
申请号 CN201210174343.0 申请日期 2012.05.31
申请人 西南交通大学 发明人 刘志刚;陈刚;张巧革
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 成都信博专利代理有限责任公司 51200 代理人 张澎
主权项 一种基于BUD谱峭度的暂态电能质量扰动分类识别方法,在电力系统中对暂态扰动进行识别,其主要步骤是:A、提取扰动特征信号输入含有扰动信号的输入电压信号u(t),t为时间,将u(t)进行小波变换,提取扰动特征信号u<sub>r</sub>(t);B、计算扰动特征信号谱峭度由A获得的扰动特征信号u<sub>r</sub>(t)求算其BUD,即Butterworth分布Cu<sub>r</sub>(t,f),再根据Cu<sub>r</sub>(t,f)的2n阶瞬时谱矩得出2阶和4阶瞬时谱距,最后根据谱峭度的定义,求得u<sub>r</sub>(n)的谱峭度<img file="FDA0000691695030000011.GIF" wi="184" he="94" />C、提取扰动特征向量选取谱峭度<img file="FDA0000691695030000012.GIF" wi="162" he="94" />的最大值、最小值和均值组成特征向量T,表示为T=(K<sub>max</sub>,K<sub>min</sub>,K<sub>avg</sub>)   (10)其中K<sub>max</sub>为谱峭度的最大值,K<sub>min</sub>为谱峭度的均值,K<sub>avg</sub>为谱峭度的均值;D、PSO优化参数的SVM分类识别将特征向量和特征向量对应的类别输入PSO优化参数的SVM进行训练,然后对测试样本进行分类,得出并向后续处理设备输出所求分类结果;A步计算扰动特征信号谱峭度的具体作法为:A1、扰动特征信号u<sub>r</sub>(n),计算其Butterworth分布,结果为Cu<sub>r</sub>(t,f),<maths num="0001" id="cmaths0001"><img file="FDA0000691695030000013.GIF" wi="1513" he="181" /></maths>A2、根据Butterworth分布Cu<sub>r</sub>(t,f)求取的2阶瞬时谱距<img file="FDA0000691695030000014.GIF" wi="126" he="87" />和4阶瞬时谱距<img file="FDA0000691695030000015.GIF" wi="167" he="91" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mover><mi>S</mi><mo>^</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mrow><mo>&lt;</mo><msup><mrow><mo>|</mo><msub><mi>Cu</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>></mo></mrow><mi>k</mi></msub></mtd></mtr><mtr><mtd><msub><mover><mi>S</mi><mo>^</mo></mover><mn>4</mn></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><msub><mrow><mo>&lt;</mo><msup><mrow><mo>|</mo><msub><mi>Cu</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>f</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>4</mn></msup><mo>></mo></mrow><mi>k</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000691695030000016.GIF" wi="804" he="241" /></maths>其中,&lt;·&gt;<sub>k</sub>表示的是k阶时平均;A3、根据谱峭度定义,求得基于BUD的谱峭度,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>K</mi><mo>^</mo></mover><msub><mi>u</mi><mi>r</mi></msub></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mover><mi>S</mi><mo>^</mo></mover><msub><mrow><mn>4</mn><mi>u</mi></mrow><mi>r</mi></msub></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><mrow><msubsup><mover><mi>S</mi><mo>^</mo></mover><msub><mrow><mn>2</mn><mi>u</mi></mrow><mi>r</mi></msub><mn>2</mn></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mn>2</mn><mo>,</mo><mi>f</mi><mo>&NotEqual;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA0000691695030000021.GIF" wi="831" he="177" /></maths>
地址 610031 四川省成都市二环路北一段111号西南交通大学科技处