主权项 |
一种基于支持向量机理论的XLPE电缆局部放电信号特征优化方法,其特征是包括以下步骤:步骤1,输入XLPE电缆局部放电信号的特征向量集合;步骤2,将XLPE电缆局部放电信号特征向量集合进行初始化操作设当前特征子集向量:s=[1,2,...,k],特征排序向量:r=[],训练样本矩阵:X<sub>0</sub>=[x<sub>1</sub>,x<sub>2</sub>,...,x<sub>n</sub>]<sup>T</sup>,类标示向量:y=[y<sub>1</sub>,y<sub>2</sub>,...,y<sub>n</sub>]<sup>T</sup>;步骤3,根据特征获取新的训练样本矩阵X=X<sub>0</sub>(:,s);步骤4,训练分类器;步骤5,根据排序准则计算各特征代价分数c,代价函数为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>J</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>α</mi><mi>T</mi></msup><mi>Hα</mi><mo>-</mo><mi>I</mi><msup><mi>α</mi><mi>T</mi></msup><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>α</mi><mi>i</mi></msub><msub><mi>α</mi><mi>j</mi></msub><msub><mi>y</mi><mi>i</mi></msub><msub><mi>y</mi><mi>j</mi></msub><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>α</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000669411200000011.GIF" wi="1366" he="157" /></maths>式中,H是元素为y<sub>i</sub>y<sub>j</sub>K(x<sub>i</sub>,x<sub>j</sub>)的矩阵,K(·)为核函数,采用高斯核函数<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mfrac><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><msup><mi>σ</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000669411200000012.GIF" wi="606" he="209" /></maths>σ为高斯核宽度;各特征的重要性根据其对该代价函数的贡献值即排序代价分数来评判;对于非线性核,为降低消去第h个特征的代价函数计算量,假定α的值保持不变而只重新计算矩阵H,则每个函数的代价分数为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>DJ</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>α</mi><mi>T</mi></msup><mi>Hα</mi><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>α</mi><mi>T</mi></msup><mi>H</mi><mrow><mo>(</mo><mo>-</mo><mi>i</mi><mo>)</mo></mrow><mi>α</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000669411200000013.GIF" wi="1168" he="133" /></maths>步骤6寻找最小代价分数特征:f=argmin(c);步骤7,更新特征排序向量:r=[s(f),r];步骤8,消去得到最小代价分数的特征:s=s(1:f‑1,f+1:length(s));步骤9,返回步骤3,重复循环,直到s为空集;步骤10,输出优化后的XLPE电缆局部放电信号的特征向量集合。 |