频率匹配的原理及其硬件实现

1、频率匹配的原理及其硬件实现作者:北京航空航天大学王建坤夏孙闻 时间:2004-12-14来源: 电子设计应用2004年第11期摘要：本文介绍了一种信号检测的方法一一频率匹配，文中详细地叙述了频率匹配算法的原理，用Matlab从理论上证明了算法的正确，性，在此基础上用可综合的RTLVerilog模型描述了算法，通过对综合后网表进行仿真，验证了频率匹配算法 硬件实现的可行性。最后，根据频率匹配的原理，实现了一种高精度数字频率检测器，并给 出了相关的精度值。关键词：信号检测；频率匹配；Matlab；高精度频率检测器图1实现频率匹配算法的结构框图图2频率匹配算法模块的内部结构曰II言对于未知信号，可以

2、通过傅立叶变换来检测信号中的频谱分布情况。但是，从 傅立叶变换的结果只能看出该信号中包含哪些频率，而对于某频率的幅度在时域的分布情况 却无从知晓，为了弄清楚某一频率的幅度在时域的分布情况，可以通过频率匹配来测定。同 样，如果知道某一信号中包含某一频率，也可以通过频率匹配将信号中该频率的幅度所在的 时域检测出来。所谓频率匹配就是使用一种频率的信号(以下称之为匹配信号)去与未知信号 (以下称之为被匹配信号)进行相关运算。因为相关运算可以用来检测两个信号的相似性，而 做相关运算的目的就是将被匹配信号中与匹配信号频率相同或相近的部分进行放大，这样就 可以弄清楚与匹配信号频率相同或相近的频率在被匹配信号

3、中的分布情兄。匹配信号必须是 一种单一频率的信号，因此最理想的匹配信号当然就是三角函数。本文就以三角函数为例， 叙述频率匹配的原理和方法。频率匹配的实现方法实际应用中的信号通常是离散的，因此，本文只讨论离散信号。可以用频率匹 配来实现离散信号的检测、识别和提取。下面介绍离散信号频率匹配的实现方法。设被匹配离散信号为X(m)，其采样频率为fs，取匹配离散信号为和，离散信 号的频率匹配运算如下：(1)(2)(3)其中N根据实际情况而定，通常取的整数倍。通过匹配运算得到的离散函数M(n)，具有如下性质：若X(n)在区间n1,n2 上的频率为w，则M(n)在区间n1,n2上将出现峰值，即在区间n1,n

4、2内M(n)的值将明显大于其它 地方的值。频率匹配的DSP实现下面简要地介绍频率匹配的DSP实现方法，频率匹配算法的主要部分是计算 M1(n)和M2(n)。频率匹配的DSP实现方法可以分成两种：单采样和块采样。在单采样模式中，每一个采样周期 内只有一个采样值到达，在每一个采样周期内也只产生一个信号输出值。因此，单采样实现 的频率匹配必须实时执行。在块采样模式中，每一个块周期内有多个采样值到达，且每一个 块周期内能产生多个信号输出值。块采样能够根据实际需要，既可以很快也可以很慢地处理 所有输入采样值。在这里将使用块采样模式来实现频率匹配，为了节省存储器和软件开销，大多 数DSP硬件都可实现循环存

5、储器，此处的频率匹配DSP也使用循环缓冲器来实现。对于输入的新的采样值，采用新的采样值进入、 老的采样值溢出的方式，这样可以节省软件开销。频率匹配的DSP实现需要3个循环缓冲器，用来分别存储X(n)、sin(n)和cos(n)。其中 X(n)在每一个块周期内需要接收多个采样值，在新的采样值进入循环缓冲器时，溢 出老的采样值。由于sin(n)和cos(n)是周期函数，因此不需要在每个块周期都向其循环缓冲器输入采样值，只需在 第一个块周期内将一个周期的采样值存到缓冲器中即可。在处理过程中，在每一个块周期内， 不需要再向循环缓冲器中装载sin(n)和cos(n)的值，只需更新循环缓冲器的变址器指针即

6、可，而对于 X(n)循环缓冲器，每一个块周期都要向其中输入新的采样值，并且需要更新循环缓 冲器的变址器指针。为了理解上述实现思想，举一个简单的例子：设采样的块长度为4,做相关运 算的点数为5,保存X(n)值的循环缓冲器长度为8;设三角离散函数的周期为5(在实际应用中，通常取做相关运算的点数=三 角离散函数的周期)，考虑到循环缓冲器的长度必须是2的幂次，所以保存sin(n)和cos(n)的循环缓冲器的长度也取8。在第一个 块周期内，第一次输出集合M1(0),M1(1),M1(2),M1(3) 按下式计算：(4)从式(4)可以看出，在第一个块周期内，需要从缓冲器得到8个采样值X(0)X(7)。同理

7、可以推出在第二个块周期内，需要从缓冲器得到 X(4)X(11)，在第三个块周期内，需要从缓冲器得到X(8广X(15)。因此在第二个块周期后， 应按顺序向循环缓冲器中装入4个新的采样值。一个块周期开始时，循环缓冲器的变址器指 针总是指向在该块周期内使用的第一个采样值，因此X(n)的循环缓冲器变址指针应该按照式(5)来计算。变址器指针=前一个变指针+块长度mod缓冲器长度(5)由于sin(n)是周期为5的函数，因此只需在第一个块周期内向sin(n)的循环 缓冲器装载5个采样值。sin(n)的循环缓冲器变址指针应该按照式(6)来计算。变址器指针=前一个变指针+块长度mod三角离散函数的周期(6)在计

8、算M2(n)时，X(n)循环缓冲器实现方法如前面所述，cos(n)循环缓冲器的 实现方法与sin(n)的循环缓冲器相同，其循环缓冲器变址指针计算也遵循式。使用Verilog HDL实现频率匹配算法频率匹配算法并不是很复杂，因此可以使用Verilog HDL实现该算法，然后将 生成的网表文件下载到FPGA中。使用VerilogHDL实现频率匹配算法的结构框图如图1所示。图1包含了两个模块，其中数据转换缓存器的作用是实现数据格式转换和数据 缓存，数据格式转换是将从外部接收到数据的格式转换成频率匹配算法模块可以处理的数据 格式，数据缓存就是将接收到的大量数据存到寄存器中，然后在每个时钟周期向频率匹配

9、算 法模块输送一个数据。频率匹配算法模块是整个系统的核心，在实现该模块时，使用了两级 流水线，大大加快了数据的处理速度。频率匹配算法模块的内部结构如图2所示，其中reg1和reg2用来保存sin(n)和cos(n)的值。频率匹配算法模块的可综合RTL Verilog模型描述的代码保存在freqmatch.v 文件中。另外，在设计的过程中，还编写了顶层测试模块文件t_freqmatch.v， 在该文件中，首先将Matlab生成的仿真数据(保存在data.txt中)读入到缓存器中，然后在每个时钟周期 内向 freqmatch模块发送一个数据。data.txt文件的数据包含三个部分，前一部分为随机噪

10、声， 中间是频率为500Hz的正弦信号，后一部分是频率为1000Hz的正弦信号。在freqmatch.v文件中，将匹配信号的频率设为500Hz。 通过仿真，freqmatch模块能够有效地将中间部分频率为500Hz的正弦信号提取出来。结语在设计的过程中，作者为该数字频率检测器编写了可综合的RTL Verilog模块，并通过了仿真和综合，其规模可以控制在合理的范围之内。另外，可以 通过增加频率匹配中作相关运算的点数来提高频率检测器的精度。并且在设计中使用了一种 改进的频率匹配算法，可以使硬件规模不会随着作相关运算的点数的增加而增大。由于频率 匹配检测器与频率滤波器很相似，因此在这里借用频率滤波器的参数来描述频率检测器，表 1和表2给出了根据频率匹配原理已经实现的几种频率检测器的参数(表中的通带指以该匹配 频率进行检测时所能检测出的频率段)。最后，需要指出的是本文所讨论的频率检测器与频率滤波器很相似，但它们有 不同之处。使用理想频率检测器对信号进行处理时，能够将某一频率信号提取出来，但抹掉 了原始信号中该频率信号振幅上的差异，即在输出信号中，该频率信号振幅是一个恒定值。 而使用理想滤波