在混频周期信号中提取特定频率的波的数字滤波方法
【技术领域】
[0001]本发明涉及一种数字滤波方法，特别涉及一种混频信号中提取特定频率波的幅值及相位的数字滤波算法。
【背景技术】
[0002]在电力设备检测系统中，消弧线圈及小电流接地选线多采用注入信号法，因其高效越来越被广泛应用。所谓注入信号法，实际上就是在电网上注入一个几十毫安不同于电网工频的特定频率正弦信号电流，然后在有关节点安装电流或电压传感器，通过采集注入的特定频率信号，从而计算出各电力系统参数或作出故障诊断。所以采用高效精确的方法将注入信号在工频电网信号中提取出来尤其重要。由于注入信号微弱，而工频信号较强，在工程上不外乎采用高选择性(高Q值)及高放大倍数的硬件和软件滤波法，其基本原理就是利用电路的频率特性实现对信号中频率成分的选择。频率滤波时，是把信号看成是由不同频率正弦波叠加而成的模拟信号，通过选择不同的频率成分来实现信号滤波。此方式硬件投入较大、电路设计复杂、可靠性和稳定性较差，同时对不同频率的干扰信号还存在阻抗匹配的问题；此外，过高的Q值和过高的放大倍数容易引起自激振荡，振荡周期长，采样值稳定所需要的时间也就长，影响采样数据的实时性。
【发明内容】
[0003]为了解决现有技术的不足，本发明提供了一种通过对采集的数据进行滤波处理，从而消除或削弱干扰和噪声的影响，提高测量的可靠性和精度，实现在有噪声的周期信号中提取特定频率的波的数字滤波方法。
[0004]本发明所采用的技术方案如下:一种在混频周期信号中提取特定频率的波的数字滤波方法，包括以下步骤:
[0005](I)通过SPI通讯从AD采样芯片中读取需要进行数字滤波的数据；
[0006](2)根据采样频率和所求的信号分辨率，以及各混频信号的频率，在步骤⑴中读取的需要进行数字滤波的数据中确定采样周期和采样点个数；
[0007](3)利用数字锁相放大滤波算法，将步骤(2)中提取数据与设计的锁相滤波矩阵相乘；
[0008](4)将步骤(3)中的数据进行周期累加后得出常数，然后求常数的均方根和商的反正切，得到最终需要的滤波并输出。
[0009]其中，SPI为串行外设接口 ；AD为电路专业术语，AD采样芯片是模拟信号转换为数字信号的芯片。
[0010]作为优选，所述数字锁相放大滤波算法包括以下步骤:
[0011](I)根据采样频率和所求的信号分辨率，以及各混频信号的频率选择N个采样点，得到一个NX I阶矩阵F ;
[0012](2)根据最终所求的信号频率，设计与所要求信号频率相同的正弦NX I阶矩阵和余弦NX I阶矩阵f2，其中匕和f 2作为传递函数矩阵；
[0013](3)将步骤⑴中矩阵F分别与步骤⑵中的传递函数矩阵f#Pf2点乘后得到相同个数的N个点，分别存储为NX I阶矩阵FjP F 2;
[0014](4)再分别将步骤(3)中的NX I阶矩阵FJPFj^N个元素累加后求平均值，因采样时间为各混频周期的公倍数，可得到常数QJP Q 2。
[0015]本发明提供的技术方案带来的有益效果是:克服了上百倍的噪声干扰，获得真实的幅值，且电路简单，便于操作，可靠性高，并且不存在阻抗匹配、特性波动、非一致性等问题；只要适当改变数字滤波程序有关参数，就能方便的改变滤波特性，因此数字滤波使用时方便灵活。通过数字锁相放大滤波算法，可以有效地消除其他周期信号对所需频率信号的干扰，并能够根据所求信号频率调整滤波参数，从而提高了数字滤波的方便性，灵活性以及可靠性。
[0016]由此可见，本发明与现有技术相比，具有突出的实质性特点和显著的进步，其实施的有益效果也是显而易见的。
【具体实施方式】
[0017]为了使本发明的目的、技术方案及优点更加清楚明白，以下结合具体实施例，对本发明进行进一步详细说明。此处所描述的具体实施例仅用以解释本发明，但并不用于限定本发明。
[0018]实施例一
[0019]—种在混频周期信号中提取特定频率的波的数字滤波方法，包括以下步骤:
[0020](I)混频周期信号 Y = B+A*sin(2 Ji cot+O+AfsinQ Ji c^t+Cl)+...+Ak*sin(2 3i ?kt+Ck)(其中，B为信号零点，即信号以该点上下对称.Mk1……Ak为振幅；C、
C1......Ck为初始相位；t为时间，ω为频率，Y为混频信号)先经过AD采样芯片采样后，再通过SPI将采样数据传输到CPU芯片。
[0021](2)根据采样频率和所求的信号分辨率，以及各混频信号的频率选择N个采样点，得到一个NX I阶矩阵F ;如何确定点数N，需要根据所求的信号分辨率η来确定。其中，所谓信号分辨率就是能够用该算法分析得到所需频率信号幅值和相位的每周期最小采样点数。根据所需要提取的信号分辨率并结合AD芯片的采样频率，将每周期最小采样点数定为η完全能够满足该算法的分辨率要求。
[0022](3)米样时间由混频彳目号的周期决定，取各混频频率和所需频率彳目号频率的差及和，再算出各差值及和值频率对应周期的最小公倍数作为最终采样时间。
[0023](4)若要在混频中提取出频率为ω的信号，则将Y分别与频率为ω的正弦参考信号yi= sin (2 Ji cot)和余弦参考信号y 2= cos (2 π ω t)的NX I阶矩阵f f2，作为传递函数矩阵，该步骤中参考信号的频率ω由CPU时钟频率和采样频率控制决定；将矩阵F分别与传递函数矩阵匕和f 2点乘后得到相同个数的N个点，分别存储为NX I阶矩阵F
F20
[0024](5)再分别将步骤⑷中的NX I阶矩阵FJPFJ^N个元素累加后求平均值，因采样时间为各混频周期的公倍数，可得到常数QJP Q 2。
[0025](6)最后求QjP Q 2的均方根和商的反正切，作为最终需要的信号的幅值和初始相位。
[0026]实施例二:在50HZ和83.3HZ的混频信号中提取出83.3HZ信号的幅值和相位。
[0027](I)混频周期信号 Y = B+A*sin(2 Ji *83.St+O+AfsinQ Ji *5(^+0 (其中，B 为信号零点，即信号以该点上下对称；A和A1为振幅；C和C i为初始相位；t为时间；Y为混频信号)先经过AD采样芯片采样后，再通过SPI将采样数据传输到CPU芯片。
[0028](2)根据采样频率和要求的信号分辨率，以及各混频信号的频率选择N个采样点，得到一个NX I阶矩阵F ;如何确定点数N，需要根据要求的信号分辨率η来确定。所谓信号分辨率就是能够用该算法分析得到所需频率信号幅值和相位的每周期最小采样点数；在本实施例中根据83.3ΗΖ信号分辨率并结合AD芯片的采样频率，将每周期最小采样点数定为η = 128完全能够满足该算法的分辨率要求。
[0029](3)米样时间由混频彳目号的周期决定，取各混频频率和所需频率彳目号频率的差及和，再算出各差值及和值频率对应周期的最小公倍数作为最终采样时间；本实施例中采样时间确定为83.3ΗΖ、33.3ΗΖ、133.3ΗΖ频率对应周期的最小公倍数60ms ；60ms是83.3HZ信号周期的5倍，因此可得出采样点N = n*5 = 640。
[0030](4)设计与要提取83.3HZ信号频率相同的正弦信号yi= sin (2 Ji *83.3t)和余弦信号y2= COS (2 31 *83.3t)的NX I阶矩阵f JP f 2，作为传递函数矩阵；将矩阵F分别与传递函数矩阵匕和f 2点乘后得到相同个数的N个点，分别存储为NX I阶矩阵F F 2；
[0031](5)分别将NXl阶矩阵FjPFJ^N个元素累加，由于采样周期为各混频83.3HZ、33.3HZU33.