[发明专利]一种差分滤波器加全周复小波幅值计算方法

摘要

本发明提供了一种分滤波器加全周Morlet复小波幅值计算方法，其包括以下步骤：(1)微机保护装置上电；(2)赋初值给k、频宽调节系数c、在一个周期内对电力元件的周期信号的采样点数N；(3)由Morlet复小波离散表达式和公式，求得比例系数τ；其中，n＝1，2，...，N－1；(4)由公式，求得Morlet复小波的实部；(5)由公式求得Morlet复小波的虚部；(6)以f_s为采样频率，采样获得电力元件的周期信号在一个周期内的离散序列f(n)；其中，f_s＝N×f，f是所述周期信号的基波频率；(7)通过差分滤波器方程g(n)＝f(n)－f(n－N/2)，获得离散序列g(n)；(8)由公式，求出幅值。本发明提供的方法，可以滤除衰减直流分量，且当电网频率发生波动时计算结果误差小。

主权项

1、一种差分滤波器加全周Morlet复小波幅值计算方法，其特征在于，其包括以下步骤：(1)微机保护装置上电；(2)赋初值给k、频宽调节系数c、在一个周期内对电力元件的周期信号的采样点数N；(3)由Morlet复小波离散表达式$\psi \left(n\right)=\tau \times \exp [-{(\frac{6n}{N}-3)}^2/c]\times \exp \left(\frac{j2\mathrm{\pi kn}}{N}\right)和公式$ \tau =N/\underset{n=0}{\overset{N-1}{\Sigma }}|\psi \left(n\right)/\tau |n,求得比例系数\tau ；其中，n＝1，2，...，N-1；(4)由公式$ {\psi }_R\left(n\right)=\tau \times \exp [-{(\frac{6n}{N}-3)}^2/c]\times \cos \left(\frac{2\mathrm{\pi kn}}{N}\right),求得Morlet复小波的实部；(5)由公式$ {\psi }_I\left(n\right)=\tau \times \exp [-{(\frac{6n}{N}-3)}^2/c]\times \sin \left(\frac{2\mathrm{\pi kn}}{N}\right)求得Morlet复小波的虚部；(6)以fs为采样频率，采样获得电力元件的周期信号在一个周期内的离散序列f(n)；其中，fs＝N\times f，f是所述周期信号的基波频率；(7)通过差分滤波器方程g(n)＝f(n)-f(n-N/2)，获得离散序列g(n)；(8)由公式$ F_k\left(n\right)=\sqrt{{(g\left(n\right)*{\psi }_R\left(n\right))}^2+{(g\left(n\right)*{\psi }_I\left(n\right))}^2},求出幅值；当k＝1时，所述幅值是电力元件的周期信号的基波幅值；当k＝i时，所述幅值是电力元件的周期信号的i次谐波幅值，其中，i为整数且i大于1。$$$$$