主权项 |
1.一种基于线性滤波器的三相四线制谐波电流检测方法,其特征在于,包括如下步骤: (a)利用三相四线制的结构特点,将其各相看成三个独立的单相,不平衡三相负载电流i<sub>a</sub>,i<sub>b</sub>,i<sub>c</sub>,各使用一个二维线性滤波器抽取出基波电流;设电网的基波频率为f,其角频率为ω,滤波器输入的电力信号为i(t),则有: <img file="FDA00002778494200011.GIF" wi="447" he="64" /><img file="FDA00002778494200012.GIF" wi="1888" he="93" />将状态变量选择为x(t)=[x<sub>1</sub>(t) x<sub>2</sub>(t)]<sup>T</sup>:<img file="FDA00002778494200013.GIF" wi="294" he="53" /><img file="FDA00002778494200014.GIF" wi="695" he="62" />则i(t)的估计值<img file="FDA00002778494200015.GIF" wi="93" he="82" />可表示为:<img file="FDA00002778494200016.GIF" wi="1127" he="70" />信号的估计误差等于:<img file="FDA00002778494200017.GIF" wi="1361" he="76" />根据最小方差思想,定义代价函数为:<img file="FDA00002778494200018.GIF" wi="955" he="67" />按照梯度下降法,有: <img file="FDA00002778494200019.GIF" wi="1954" he="111" />矩阵Ν调节x(t)收敛到最优值的速度以及算法的稳定性,令μ<sub>1</sub>,μ<sub>2</sub>>0,选择: <img file="FDA000027784942000110.GIF" wi="1925" he="139" />由此,可以得到周期系数线性动力系统:<img file="FDA000027784942000111.GIF" wi="867" he="60" />其中:<img file="FDA000027784942000112.GIF" wi="434" he="65" />x=[x<sub>1</sub>(t) x<sub>2</sub>(t)]<sup>T</sup><img file="FDA000027784942000113.GIF" wi="975" he="170" /><img file="FDA000027784942000114.GIF" wi="359" he="166" />作旋转变换,令 <img file="FDA00002778494200021.GIF" wi="1866" he="148" />则由式(1)和式(3)知,i(t)的估计值<img file="FDA00002778494200022.GIF" wi="82" he="73" />为:<img file="FDA00002778494200023.GIF" wi="1867" he="66" />将式(9)代入式(8),有: <img file="FDA00002778494200024.GIF" wi="1065" he="160" /><img file="FDA00002778494200025.GIF" wi="1694" he="161" />令μ<sub>1</sub>=μ<sub>2</sub>=μ,η=ω,得到二维线性滤波算法: <img file="FDA00002778494200026.GIF" wi="1865" he="149" />其中:y(t)=[y<sub>1</sub>(t) y<sub>2</sub>(t)]<sup>T</sup>为状态变量;μ,η为滤波器的参数,均大于零的实常数;y<sub>1</sub>(t)为滤波器输出,它表示输入信号中所含的角频率为ω的正弦分量,也即输入信号的基波分量,y<sub>2</sub>(t)为与y<sub>1</sub>(t)正交且超前的同频信号; 设初始条件y<sub>1</sub>(0)=0,y<sub>2</sub>(0)=0,取i(t)为输入,y<sub>1</sub>(t)为输出,则二维线性滤波器的传递函数为 <img file="FDA00002778494200027.GIF" wi="1858" he="128" />令s=jω,有: <img file="FDA00002778494200028.GIF" wi="1831" he="139" />则其幅频特性和相频特性分别为: <img file="FDA00002778494200031.GIF" wi="1831" he="149" /><img file="FDA00002778494200032.GIF" wi="1830" he="131" />当取μ=0.5,η=ω=314rad/s时,有: <img file="FDA00002778494200033.GIF" wi="1831" he="135" /><img file="FDA00002778494200034.GIF" wi="1831" he="119" />(b)不平衡的基波电流i<sub>af</sub>,i<sub>bf</sub>,i<sub>cf</sub>通过瞬时对称分量计算,得到对称的基波正序电流分量i<sub>afp</sub>,i<sub>bfp</sub>,i<sub>cfp</sub>;设不平衡三相负载电流瞬时值为: <img file="FDA00002778494200035.GIF" wi="1757" he="244" />其中:i<sub>a</sub>(t),i<sub>b</sub>(t),i<sub>c</sub>(t)分别为三相电流瞬时值;I<sub>am</sub>,I<sub>bm</sub>,I<sub>cm</sub>分别为三相电流的幅值;<img file="FDA00002778494200036.GIF" wi="191" he="44" />分别为三相电流的初相位;设与三相电流i<sub>a</sub>(t),i<sub>b</sub>(t),i<sub>c</sub>(t)所对应的旋转相量分别为<img file="FDA00002778494200037.GIF" wi="180" he="63" />有:<img file="FDA00002778494200038.GIF" wi="1842" he="257" />由式(20)知,相量<img file="FDA00002778494200039.GIF" wi="170" he="69" />的虚部系数即为三相电流瞬时值,只要能确定出其实部,就可以确定这些旋转相量;设a相电流的瞬时值可表示为:<img file="FDA00002778494200041.GIF" wi="770" he="59" /><img file="2.GIF" wi="403" he="24" /><img file="FDA00002778494200044.GIF" wi="650" he="58" />也即: <img file="FDA00002778494200045.GIF" wi="577" he="61" /><img file="FDA00002778494200046.GIF" wi="1704" he="65" />故有: <img file="FDA00002778494200047.GIF" wi="533" he="57" /><img file="FDA00002778494200048.GIF" wi="1591" he="122" />又因为: <img file="FDA00002778494200049.GIF" wi="390" he="57" /><img file="FDA000027784942000410.GIF" wi="631" he="57" /><img file="3.GIF" wi="403" he="19" /><img file="FDA000027784942000413.GIF" wi="637" he="57" />整理后可得: <img file="FDA000027784942000414.GIF" wi="1817" he="122" />同理可得: <img file="FDA000027784942000415.GIF" wi="1825" he="135" /><img file="FDA000027784942000416.GIF" wi="1816" he="124" />将式(25)~(27)代入式(20)中,构成瞬时旋转相量,则不平衡三相负载电流正序分量值为: <img file="FDA00002778494200051.GIF" wi="1793" he="291" />其中:Im{}为求复数的虚部运算;a=e<sup>j2π/3</sup>,a<sup>2</sup>=e<sup>-j2π/3</sup>,1+a+a<sup>2</sup>=0; (c)与原负载电流运算后,即可分离出谐波电流分量i<sub>ah</sub>,i<sub>bh</sub>,i<sub>ch</sub>。 |