发明名称 一种对电力系统电压电流信号进行分离的方法及装置
摘要 一种对电力系统电压电流信号进行分离的方法及装置,属于信号分离技术领域,装置由下位机和上位机组成,下位机包括传感器、信号调理模块、A/D转换模块、双单片机系统单元、存储模块、键盘模块、通信模块,其中传感器、信号调理模块、A/D转换模块、双单片机系统单元依次连接,存储模块、键盘操作单元模块、通信模块分别与双单片机系统单元相连,上位机与通信模块相连。方法由软件部分完成,软件部分嵌入双单片机系统单元,分离矩阵与源信号相乘得到最终分离信号表达式。本发明降低了存储空间并提高了记录的精度,可以在任意小的时间里以任意精度逼近实际波形,可以很好的分析信号中谐波的组成,为谐波分析提供了数据基础。
申请公布号 CN101334428A 申请公布日期 2008.12.31
申请号 CN200810011419.1 申请日期 2008.05.16
申请人 东北大学;鞍山华夏电力电子设备有限公司 发明人 孙秋野;张化光;王占山;杨东升;王云爽;季策;李腾;李治平;刘国威
分类号 G01R19/00(2006.01);G01R5/12(2006.01);G06F17/00(2006.01) 主分类号 G01R19/00(2006.01)
代理机构 沈阳东大专利代理有限公司 代理人 朱光林
主权项 1、一种对电力系统电压电流信号进行分离的方法,其特征在于依赖装置,包括以下步骤:步骤1、开始;步骤2、采集源信号利用装置采集电力系统的电压、电流;步骤3、确定分离网络结构采集到信号后,根据要分离的混叠信号的个数来选择网络结构进行分离;步骤4、设定非线性混叠模型矩阵分离出信号作为源信号的估计信号存在n个来自信号源的统计独立的信号,s<sub>i</sub>(t),i=1~n以及相同数目的混合信号x<sub>i</sub>(t),信号的混合是线性瞬时的,混合信号为<maths num="0001"><![CDATA[<math><mrow><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mi>ij</mi></msub><msub><mi>s</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>式中a<sub>ij</sub>为混合系数;估计信号y(t)=Wx(t)=WAs(t)=Gs(t)全局矩阵G=PD其中:P为nxn维置换矩阵,D为nxn维对角矩阵,则y(t)=PDs(t),首先是建立一个以W为变元的目标函数F(W),若某个<img file="A2008100114190002C2.GIF" wi="48" he="71" />能使F(W)达到极大或者极小值,该<img file="A2008100114190002C3.GIF" wi="46" he="72" />即为所需的解;步骤5、计算源信号与其估计信号的误差把源信号s与其估计信号y的误差e=s-y作为噪声信号;步骤6、建立恢复信号的信噪比函数建立恢复信号s的信噪比函数:<maths num="0002"><![CDATA[<math><mrow><mi>SNR</mi><mo>=</mo><mn>10</mn><mi>log</mi><mfrac><mrow><mi>s</mi><mo>*</mo><msup><mi>s</mi><mi>T</mi></msup></mrow><mrow><mi>e</mi><mo>*</mo><msup><mi>e</mi><mi>T</mi></msup></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>s</mi><mo>*</mo><msup><mi>s</mi><mi>T</mi></msup></mrow><mrow><mrow><mo>(</mo><mi>s</mi><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mo>*</mo><msup><mrow><mo>(</mo><mi>s</mi><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mi>T</mi></msup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>步骤7、建立最大信噪比目标函数用估计信号的滑动平均y代替源信号,得到最大信噪比目标函数如下:<maths num="0003"><![CDATA[<math><mrow><mi>F</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mi>SNR</mi><mo>=</mo><mn>10</mn><mi>log</mi><mfrac><mrow><mi>y</mi><mo>*</mo><msup><mi>y</mi><mi>T</mi></msup></mrow><mrow><mrow><mo>(</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mo>*</mo><msup><mrow><mo>(</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mi>T</mi></msup></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>式中y=Wx,y=Wx,x为混合信号x经滑动平均处理后的信号;步骤8、优化目标函数最大信噪比目标函数如下形式:<maths num="0004"><![CDATA[<math><mrow><mi>F</mi><mrow><mo>(</mo><mi>W</mi><mo>,</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mn>10</mn><mi>log</mi><mfrac><mrow><mi>y</mi><mo>*</mo><msup><mi>y</mi><mi>T</mi></msup></mrow><mrow><mrow><mo>(</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mo>*</mo><msup><mrow><mo>(</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mi>T</mi></msup></mrow></mfrac></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mo>=</mo><mn>10</mn><mi>log</mi><mfrac><mrow><mi>Wx</mi><mo>*</mo><msup><mi>x</mi><mi>T</mi></msup><msup><mi>W</mi><mi>T</mi></msup></mrow><mrow><mrow><mo>(</mo><mi>s</mi><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mo>*</mo><msup><mrow><mo>(</mo><mi>s</mi><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mi>T</mi></msup></mrow></mfrac></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mo>=</mo><mn>10</mn><mi>log</mi><mfrac><mrow><mi>WC</mi><msup><mi>W</mi><mi>T</mi></msup></mrow><mrow><mi>W</mi><mover><mi>C</mi><mo>&OverBar;</mo></mover><msup><mi>W</mi><mi>T</mi></msup></mrow></mfrac></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mo>=</mo><mn>10</mn><mi>log</mi><mfrac><mi>V</mi><mi>U</mi></mfrac></mrow></math>]]></maths>步骤9、计算分离矩阵梯度分离矩阵W的函数F的梯度为<maths num="0008"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>F</mi></mrow><mrow><mo>&PartialD;</mo><mi>W</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mn>2</mn><mi>W</mi></mrow><mi>V</mi></mfrac><mi>C</mi><mo>-</mo><mfrac><mrow><mn>2</mn><mi>W</mi></mrow><mi>U</mi></mfrac><mover><mi>C</mi><mo>&OverBar;</mo></mover><mo>,</mo></mrow></math>]]></maths>求得分离矩阵<img file="A2008100114190003C6.GIF" wi="69" he="76" />步骤10、分离矩阵与源信号相乘得到最终分离信号表达式分离矩阵与源信号相乘得到最终分离信号表达式;步骤11、将电压电流信号的基波及其谐波的传递到上位机并显示。
地址 110004辽宁省沈阳市和平区文化路3号巷11号
您可能感兴趣的专利