主权项 |
一种基于改进型目标函数的抗差励磁系统参数辨识方法,其特征在于,具体包括以下步骤:步骤1:利用Matlab/Simulink工具箱搭建待辨识励磁系统的可变参数的传递函数模型,以相量测量单元PMU实测的机端电压和电流相量以及励磁电流作为此模型的输入数据;步骤2:将具有抗差能力的改进型目标函数作为励磁系统参数辨识的目标函数;改进型目标函数的数学形式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>min</mi><mi>J</mi><mrow><mo>(</mo><mi>α</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>ρ</mi><mrow><mo>(</mo><mi>e</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>ρ</mi><mrow><mo>(</mo><msub><mi>E</mi><mi>fdm</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>E</mi><mi>fdc</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>,</mo><mi>α</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000551189720000011.GIF" wi="948" he="305" /></maths>抗差准则函数ρ(·)为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>ρ</mi><mrow><mo>(</mo><mi>e</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mn>2</mn></msup><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><msub><mo>|</mo><mrow><mo>|</mo><mi>e</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>≤</mo><mi>Δσ</mi></mrow></msub><mo>+</mo><msup><mrow><mo>(</mo><mi>Δσ</mi><mo>)</mo></mrow><mn>2</mn></msup><msub><mo>|</mo><mrow><mo>|</mo><mi>e</mi><mrow><mo>(</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>></mo><mi>Δσ</mi></mrow></msub></mrow>]]></math><img file="FDA0000551189720000012.GIF" wi="928" he="142" /></maths>其中,α表示待辨识的参数集合;N为采样点总数,t<sub>k</sub>为采样时刻;E<sub>fdm</sub>(t<sub>k</sub>)为第t<sub>k</sub>时刻的PMU实测励磁电压值,E<sub>fdc</sub>(t<sub>k</sub>,α)为第t<sub>k</sub>时刻的待辨识励磁系统的传递函数模型计算输出的励磁电压值;e(t<sub>k</sub>)为第t<sub>k</sub>时刻的计算输出的励磁电压值与实测励磁电压值的误差;Δσ是用来平衡抗差性和有效性的正实数;步骤3:基于Matlab/Simulink的待辨识励磁系统的传递函数模型,利用遗传算法辨识励磁系统的参数。 |