1.一种基于误差相关熵的三相电力系统的频率估计方法,其步骤如下:A、信号采集
采集三相电力系统的电压信号,并通过克拉克变换得到当前时刻n的复电压信号v(n);
B、复电压信号的估计
由基于宽线性模型的自适应滤波器,输出下一时刻n+1的复电压信号估计值其中,h(n)表示滤波器当前时刻n的标准权系数,其初始值为0.99+0.15i,其中i为虚数单位;g(n)表示滤波器当前时刻n的共轭权系数,其初始值为0;v*(n)是当前时刻n的复电压信号v(n)的共轭值;
C、复电压信号的误差计算
由当前时刻n的复电压信号v(n)与当前时刻n的复电压信号估计值 计算出当前时刻n的复电压信号的误差e(n), 误差e(n)的初值为0.01-0.1;
D、权系数更新
D1、误差相关熵的计算
将前L-1时刻与当前时刻n共L个时刻构成当前时间窗,并将当前时间窗内各个时刻的误差e(n-L+1)、e(n-L+2)、…e(n-1)、e(n),依次按当前时间窗内的时刻顺序重命名为e1、e2、e3…ej…ek…eL,进而构成当前时间窗的误差向量En,En=[e(n-L+1)、e(n-L+2)、…e(n-1)、e(n)]=[e1、e2、e3…ej…ek…eL];其中,ej、ek分别为当前时间窗内的第j个时刻的误差和第k个时刻的误差,ej=e(n-L+j),ek=e(n-L+k),j=1、2……L,k=1、2……L;其中,L为当前时间窗的时刻总数,其取值为15-25;
计算当前时间窗中第j个时刻与第k个时刻的误差差值Δejk,Δejk=ej-ek,进而得到当前时间窗中第j个时刻与第k个时刻的误差相关熵βjk:其中,σ为核宽度,取值为0.8-1.2;
D2、权系数更新
更新计算出下一时刻n+1的标准权系数h(n+1),其中,v*(j)、v*(k)分别为当前时间窗内的第j个和第k个时刻的复电压信号的共轭值,v*j=v*(n-L+j),v*k=v*(n-L+k);
更新计算出下一时刻n+1的共轭权系数g(n+1),其中,μ为步长,其取值为0.01-0.1;v(j)、v(k)分别为当前时间窗内的第j个时刻和第k个时刻的复电压信号,Vj=v(n-L+j),Vk=v(n-L+k);
E、频率估计
由下式得出三相电力系统当前时刻n的频率估计值
其中,Im(*)代表取复数的虚部,ΔT表示采样间隔,其取值为0.0001~0.001s;
F、令n=n+1,重复A~E的操作,即动态实时估计出三相电力系统的频率值。