发明名称 一种挖泥船电网谐波监测方法与系统
摘要 本发明实施例公开了一种挖泥船电网谐波监测方法,包括:获取电网中的原始电信号,并提取电网中的原始电信号的各次谐波量测量;根据k时刻第h次谐波的量测量<img file="DDA0000649649500000011.GIF" wi="91" he="99" />计算k时刻第h次谐波幅值的最优估计值<img file="DDA0000649649500000012.GIF" wi="107" he="89" />计算<img file="DDA0000649649500000013.GIF" wi="80" he="90" />的协方差矩阵<img file="DDA0000649649500000014.GIF" wi="113" he="90" />根据最优估计值<img file="DDA0000649649500000015.GIF" wi="74" he="90" />计算k时刻谐波畸变率(THD)<sub>k</sub>;根据每个时刻的谐波畸变率,判断出电网的电信号质量。本发明实施例在常规卡尔曼滤波算法的基础上,引入修正因子λ表示疏浚作业环境,建立模糊自适应寻优推理模型基础上的鲁棒卡尔曼滤波算法,完成电网谐波数据的一步预测与估计,显著提高了电网谐波监测系统的鲁棒性,能够使得谐波畸变率的最优估计值更加接近原始电信号,最后再根据谐波畸变率来判断电网的质量,提高了判断的准确性。
申请公布号 CN104485669B 申请公布日期 2016.10.26
申请号 CN201410855429.9 申请日期 2014.12.31
申请人 广州航海学院 发明人 白明;周晓莹;邹福顺;陈锦华;钟碧良;黄富华;赖伟生;李芳辉;谢程宁;米志红;邱伟成;叶永强;黄康容;郑泳然
分类号 H02J3/01(2006.01)I 主分类号 H02J3/01(2006.01)I
代理机构 广州三环专利代理有限公司 44202 代理人 颜希文
主权项 一种挖泥船电网谐波监测方法,其特征在于,包括:获取电网中的原始电信号,并提取电网中的原始电信号的各次谐波量测量;根据k时刻第h次谐波的量测量<img file="FDA0000933094520000011.GIF" wi="73" he="75" />计算k时刻第h次谐波幅值的最优估计值<img file="FDA0000933094520000012.GIF" wi="100" he="68" /><maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>=</mo><msubsup><mi>A</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mrow><mo>(</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>B</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mrow><mo>(</mo><msubsup><mi>u</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>w</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>=</mo><mo>&lsqb;</mo><msubsup><mi>A</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><msubsup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>h</mi></msubsup><msup><mrow><mo>(</mo><msubsup><mi>A</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup><mo>&rsqb;</mo><mo>+</mo><msub><msup><mi>Q</mi><mi>h</mi></msup><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>G</mi><mi>k</mi><mi>h</mi></msubsup><mo>=</mo><msubsup><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><msubsup><mi>C</mi><mi>k</mi><mi>h</mi></msubsup><mo>/</mo><mo>&lsqb;</mo><msubsup><mi>C</mi><mi>k</mi><mi>h</mi></msubsup><msubsup><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><msup><mrow><mo>(</mo><msubsup><mi>C</mi><mi>k</mi><mi>h</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup><mo>+</mo><msubsup><mi>R</mi><mi>k</mi><mi>h</mi></msubsup><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mover><mi>x</mi><mo>~</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>=</mo><mi>&lambda;</mi><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&lambda;</mi><mo>)</mo></mrow><msubsup><mover><mi>x</mi><mo>&OverBar;</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi></mrow><mi>h</mi></msubsup><mo>=</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>+</mo><msubsup><mi>G</mi><mi>k</mi><mi>h</mi></msubsup><mrow><mo>(</mo><msubsup><mi>z</mi><mi>k</mi><mi>h</mi></msubsup><mo>-</mo><msubsup><mi>D</mi><mi>k</mi><mi>h</mi></msubsup><msubsup><mi>u</mi><mi>k</mi><mi>h</mi></msubsup><mo>-</mo><msubsup><mi>C</mi><mi>k</mi><mi>h</mi></msubsup><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000933094520000013.GIF" wi="814" he="500" /></maths>其中,式中k为时刻序列;h为谐波阶次,h=6·l±1,l为正整数;<img file="FDA0000933094520000014.GIF" wi="101" he="67" />为k时刻谐波幅值的一步预测值;<img file="FDA0000933094520000015.GIF" wi="133" he="68" />为k‑1时刻谐波幅值的最优估计值;<img file="FDA0000933094520000016.GIF" wi="68" he="68" />为k时刻谐波幅值的最优估计值;<img file="FDA0000933094520000017.GIF" wi="80" he="65" />为k‑1时刻的系统传输矩阵;<img file="FDA0000933094520000018.GIF" wi="81" he="65" />为k‑1时刻的控制输入矩阵;<img file="FDA0000933094520000019.GIF" wi="78" he="64" />为k‑1时刻的确定性控制输入;<img file="FDA00009330945200000110.GIF" wi="86" he="70" />为k‑1时刻的系统噪声;<img file="FDA00009330945200000111.GIF" wi="110" he="78" />为<img file="FDA00009330945200000112.GIF" wi="109" he="78" />的协方差矩阵;<img file="FDA00009330945200000113.GIF" wi="165" he="79" />为<img file="FDA00009330945200000114.GIF" wi="137" he="78" />的协方差矩阵;<img file="FDA00009330945200000115.GIF" wi="94" he="79" />为<img file="FDA00009330945200000116.GIF" wi="103" he="86" />的协方差矩阵;<img file="FDA00009330945200000117.GIF" wi="62" he="69" />为k时刻的卡尔曼增益矩阵;<img file="FDA00009330945200000118.GIF" wi="62" he="76" />为k时刻的量测矩阵;<img file="FDA00009330945200000119.GIF" wi="61" he="71" />为量测噪声<img file="FDA00009330945200000120.GIF" wi="53" he="78" />的协方差矩阵;<img file="FDA00009330945200000121.GIF" wi="118" he="84" />为k时刻谐波幅值的鲁棒卡尔曼预测值;λ为修正因子;<img file="FDA00009330945200000122.GIF" wi="93" he="78" />为k‑1时刻谐波幅值的经验预测平均值;<img file="FDA00009330945200000123.GIF" wi="70" he="80" />为k时刻输入输出矩阵;计算<img file="FDA00009330945200000124.GIF" wi="69" he="77" />的协方差矩阵<img file="FDA00009330945200000125.GIF" wi="99" he="71" /><maths num="0002"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi></mrow><mi>h</mi></msubsup><mo>=</mo><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msubsup><mi>G</mi><mi>k</mi><mi>h</mi></msubsup><msubsup><mi>C</mi><mi>k</mi><mi>h</mi></msubsup><mo>)</mo></mrow><msubsup><mi>P</mi><mrow><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>h</mi></msubsup><mo>;</mo></mrow>]]></math><img file="FDA00009330945200000126.GIF" wi="493" he="85" /></maths>I为单位矩阵;根据不同阶次谐波幅值的最优估计值<img file="FDA00009330945200000127.GIF" wi="65" he="68" />计算k时刻的谐波畸变率(THD)<sub>k</sub>:<maths num="0003"><math><![CDATA[<mrow><msub><mrow><mo>(</mo><mi>T</mi><mi>H</mi><mi>D</mi><mo>)</mo></mrow><mi>k</mi></msub><mo>=</mo><mfrac><msqrt><mrow><munder><mo>&Sigma;</mo><mi>h</mi></munder><msup><mrow><mo>(</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi></mrow><mi>h</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>|</mo><mi>k</mi></mrow><mn>1</mn></msubsup></mfrac><mo>&times;</mo><mn>100</mn><mi>%</mi><mo>;</mo></mrow>]]></math><img file="FDA00009330945200000128.GIF" wi="638" he="237" /></maths>其中<img file="FDA00009330945200000129.GIF" wi="74" he="74" />为k时刻基波幅值的最优估计值;根据每个时刻的谐波畸变率,判断出电网的电信号质量。
地址 510725 广东省广州市黄埔区红山三路101号