基于FPGA与MATLAB的超声多普勒频移解调应用

本系统基于FPGA技术，围绕快速傅里叶变换（FFT）与直接数字频率合成（DDS）核心功能展开设计，主要应用于超声多普勒频移解调场景。系统通过DDS模块生成特定频率的信号，再经由FFT模块对信号进行频域分析，最终实现对超声多普勒频移信号的处理与解调。整体代码基于Verilog/VHDL硬件描述语言开发，包含测试激励模块、DDS编译器IP核及相关辅助模块，可在Xilinx Vivado开发环境中进行仿真与部署，具备高实时性、高稳定性的特点，满足超声信号处理对精度和速度的双重需求。

一、系统概述

（一）核心流程

DDS IP核生成2MHz和(2MHz+1kHz)的sin频率信号；
乘法IP核实现2MHz和(2MHz+1kHz)的sin频率信号混频处理；
FIR IP实现低通滤波算法（Matlab获取滤波参数）；
FFT IP核实现FFT的算法；
乘法IP核实现FFT处理后数据取模运算；
通过算法实现取模运算后65536数据的峰值搜索；
对峰值结果进行计算，并与1kHz理论值比对。

二、核心模块功能解析

（一）FFT顶层测试模块（tb_fft_top.v）

1. 模块定位

该模块为FFT顶层设计的测试激励模块，主要功能是为FFT顶层模块（fft_top）提供时钟信号，并监测FFT处理后的输出结果，用于验证FFT模块在仿真环境下的功能正确性。

2. 关键参数与信号

参数定义
- SIM：仿真模式控制参数，默认值为1，用于区分仿真环境与实际硬件运行环境，在仿真时可通过该参数启用特定的测试逻辑。
- CLK_PERIOD：时钟周期参数，设置为10ns，对应50MHz时钟频率，与系统主时钟频率保持一致，确保模块间时钟同步。
信号定义
- i_clk_50m：50MHz时钟输出信号，作为FFT顶层模块的主时钟输入，由测试模块内部的时钟生成逻辑驱动。
- o_fft_process_1k：32位宽的FFT处理结果输出信号，用于观测FFT模块对输入信号的频域处理结果，可通过仿真波形查看信号的数值变化，判断FFT运算是否正常。

3. 核心逻辑

时钟生成：通过initial块初始化时钟信号为低电平，随后在always块中以CLK_PERIOD/2（即5ns）为周期翻转时钟信号，持续生成稳定的50MHz时钟，为整个FFT系统提供同步时钟源。
模块例化：例化FFT顶层模块（fft_top），并将测试模块生成的时钟信号i_clk_50m接入，同时获取FFT模块输出的处理结果o_fft_process_1k，实现测试激励与结果监测的闭环。

（二）DDS编译器IP核模块（dds_compiler_2m_2m1k）

1. 模块定位

DDS（直接数字频率合成）编译器IP核是系统的信号生成核心，用于产生高精度、高稳定度的正弦波等周期性信号，为后续FFT模块提供标准的输入信号源，可满足超声多普勒频移解调中对参考信号或激励信号的需求。该IP核由Xilinx Vivado工具生成，包含组件声明、实例化模板及相关配置逻辑。

2. 关键接口信号

信号名称	信号方向	位宽	功能描述
`aclk`	输入	1位	时钟信号，为DDS核提供同步时钟，确保信号生成的时序稳定性
`m_axis_data_tvalid`	输出	1位	数据有效标志信号，高电平时表示`m_axis_data_tdata`和`m_axis_data_tuser`信号中的数据有效
`m_axis_data_tdata`	输出	8位	数据 payload 信号，存储DDS核生成的信号数据（如正弦波采样值），位宽为8位，保证一定的数据精度
`m_axis_data_tuser`	输出	1位	用户自定义信号，可用于标识数据通道或数据帧相关信息，此处为1位宽，满足基础的用户标识需求

3. 组件声明与实例化

组件声明：通过COMPONENT关键字定义DDS核的接口规范，明确各输入输出信号的类型与位宽，为后续模块实例化提供接口模板。
实例化模板：提供了DDS核的实例化代码模板，用户可根据实际需求修改实例名称及信号连接关系，将DDS核集成到系统中，实现信号生成功能的快速调用。

（三）DDS核测试模块（tbdds_compiler_2m_2m1k.vhd）

1. 模块定位

该模块为DDS编译器IP核的专用测试激励模块，基于VHDL语言编写，用于验证DDS核生成信号的正确性、数据有效性及时序稳定性，确保DDS核输出信号符合系统设计要求。

2. 关键参数与信号

时序参数
- CLOCK_PERIOD：时钟周期，设置为100ns，对应10MHz时钟频率，与DDS核的时钟需求匹配，保证DDS核稳定工作。
- T_HOLD：信号保持时间，设置为10ns，确保输入信号在时钟上升沿后保持足够时间，满足时序约束，避免数据采样错误。
- T_STROBE：信号采样时间，设置为CLOCK_PERIOD - 1ns（即99ns），用于在时钟周期内的合适时间点对输出信号进行采样验证。
核心信号
- 时钟与控制信号：aclk为DDS核提供时钟输入；end_of_simulation为仿真结束标志信号，用于控制仿真流程的启停。
- 数据输出信号：m_axis_data_tvalid、m_axis_data_tdata、m_axis_data_tuser分别对应DDS核输出的有效标志、数据及用户信号，用于监测DDS核的输出状态。
- 别名信号：如m_axis_data_tdata_sine（正弦波数据别名）、m_axis_data_tuser_channel_id（通道ID别名）等，用于在仿真波形中更直观地观测信号含义，方便调试分析。

3. 核心测试逻辑

时钟生成：在clock_gen进程中，根据end_of_simulation信号控制时钟生成。仿真未结束时，以CLOCK_PERIOD为周期持续生成时钟信号；仿真结束后，停止时钟生成，避免不必要的资源消耗。
激励生成：stimuli进程用于生成测试激励，在时钟上升沿后延迟T_HOLD时间驱动输入信号，并控制仿真时长（运行足够生成5个输出信号周期的时间），最后设置end_of_simulation为true，结束仿真并输出成功报告。
输出验证：check_outputs进程用于验证DDS核的输出信号。在时钟上升沿后延迟T_STROBE时间采样输出信号，重点检查当m_axis_data_tvalid为高电平时，m_axis_data_tdata和m_axis_data_tuser信号是否为有效电平（非X态），若存在无效数据则输出错误报告，确保DDS核输出数据的有效性。
信号别名映射：通过组合逻辑与时序逻辑，将DDS核输出的m_axis_data_tdata和m_axis_data_tuser信号映射到对应的别名信号，并根据m_axis_data_tvalid信号的有效性更新别名信号，同时对TDM（时分复用）通道数据进行分离与存储，方便在仿真中区分不同通道的输出数据。

三、系统工作流程

（一）整体流程框架

系统工作流程以时钟同步为基础，分为信号生成、信号处理、结果输出三个核心阶段，具体流程如下：

时钟驱动：测试模块（tb_fft_top.v、tbdds_compiler_2m_2m1k.vhd）分别生成50MHz和10MHz时钟信号，为FFT模块和DDS模块提供同步时钟，确保各模块时序一致。
DDS信号生成：DDS编译器IP核在时钟驱动下，生成预设频率的周期性信号（如正弦波），通过m_axis_data_tvalid信号标识数据有效性，将数据通过m_axis_data_tdata输出至后续FFT模块。
FFT信号处理：FFT顶层模块（fft_top）接收DDS核输出的信号，在50MHz时钟同步下对信号进行FFT运算，将时域信号转换为频域信号，完成信号的频域分析与处理。
结果输出与监测：FFT模块处理后的结果通过o_fft_process_1k信号输出，测试模块可监测该信号，验证FFT运算结果的正确性；同时，DDS核的测试模块通过监测输出信号的有效性与数据完整性，确保DDS核工作正常。

（二）关键时序协同

时钟同步：DDS核与FFT模块虽可能使用不同时钟频率（如本系统中DDS为10MHz，FFT为50MHz），但在实际系统集成时需通过时钟管理模块（如PLL）实现时钟同步，避免跨时钟域数据传输错误。
数据握手：DDS核通过m_axis_data_tvalid信号与FFT模块进行数据握手，FFT模块仅在m_axis_data_tvalid为高电平时采样DDS输出数据，确保数据采样的准确性，避免无效数据进入FFT运算流程。

四、系统功能特点与应用场景

（一）功能特点

高精度信号生成：DDS编译器IP核支持生成高精度、低相位噪声的周期性信号，可通过配置灵活调整输出信号的频率、相位与幅度，满足不同场景下的信号源需求。
高效频域处理：FFT模块采用硬件实现，具备高速并行运算能力，可快速将时域信号转换为频域信号，实时性强，适用于超声多普勒频移等对实时信号处理要求较高的场景。
完善的测试验证：各核心模块均配备专用测试模块，可在仿真环境中验证模块功能的正确性，包括时钟稳定性、数据有效性、运算结果准确性等，降低硬件调试难度。
高可扩展性：系统采用模块化设计，DDS核与FFT模块可根据实际需求独立配置（如调整DDS输出信号频率、FFT运算点数等），同时支持与其他信号处理模块（如FIR滤波器）集成，扩展系统功能。

（二）应用场景

本系统主要应用于超声多普勒频移解调领域，具体应用如下：

超声信号激励：DDS核生成特定频率的超声激励信号，驱动超声探头发射超声波。
频移信号采集与处理：超声探头接收反射回波信号（含多普勒频移信息），经前端调理后输入FFT模块，FFT模块对信号进行频域分析，提取频移信息。
频移解调：通过FFT处理后的频域结果，计算多普勒频移量，进而实现血流速度等生理参数的测量，为医学超声诊断等应用提供数据支持。

此外，系统还可扩展应用于其他需要高精度信号生成与频域分析的场景，如通信系统中的信号调制解调、雷达信号处理等领域。

五、仿真与验证要点

（一）仿真环境要求

开发工具：推荐使用Xilinx Vivado（版本需支持DDS Compiler 6.0 IP核及FFT相关IP核），仿真工具可选用ModelSim或QuestaSim。
库文件依赖：仿真时需加载Xilinx Vivado的VHDL/Verilog仿真库，确保DDS IP核及其他硬件原语的正确编译。
时序约束：需在仿真配置中添加时钟周期、信号保持时间等时序约束，与实际硬件时序要求一致，避免仿真与硬件运行结果差异。

（二）验证核心要点

DDS核验证
- 数据有效性：监测m_axis_data_tvalid信号，确保其在数据输出期间保持高电平，无异常跳变。
- 信号频率准确性：通过仿真波形测量m_axis_data_tdata信号的周期，计算输出信号频率，验证是否与预设频率一致。
- TDM通道分离：若DDS核配置为多通道输出，需验证各TDM通道数据分离的正确性，确保不同通道数据无串扰。
FFT模块验证
- 运算结果正确性：向FFT模块输入已知频率的信号（如DDS核生成的正弦波），对比FFT输出结果与理论频域结果，验证FFT运算精度。
- 时序稳定性：监测FFT模块输入输出信号的时序关系，确保无数据延迟、丢失或错拍现象。
系统集成验证
- 跨模块数据传输：验证DDS核与FFT模块之间的数据传输正确性，确保m_axis_data_tvalid信号有效时，FFT模块能正确接收并处理数据。
- 长时间运行稳定性：进行长时间仿真（如生成1000个信号周期），验证系统无死锁、数据溢出等异常情况，确保系统稳定运行。

六、总结与扩展建议

（一）系统总结

本系统基于FPGA实现了DDS信号生成与FFT频域处理的核心功能，通过模块化设计与完善的测试验证，确保了信号生成的高精度与信号处理的高效性，可满足超声多普勒频移解调等场景的需求。系统各模块功能明确、时序协同清晰，为后续硬件部署与功能扩展奠定了基础。

（二）扩展建议

功能扩展
- 集成FIR滤波器：在DDS核与FFT模块之间添加FIR滤波器，对DDS生成的信号进行滤波去噪，进一步提升FFT处理的精度。
- 多通道扩展：配置DDS核为多通道输出模式，结合FFT模块的多通道处理能力，支持多路超声信号的同时处理，提升系统吞吐量。
性能优化
- 时钟频率提升：通过FPGA的PLL模块提升系统时钟频率（如将FFT模块时钟提升至100MHz），进一步提高FFT运算速度，满足更高实时性需求。
- 数据位宽调整：根据实际信号精度需求，调整DDS核输出数据位宽（如从8位扩展至16位）与FFT模块运算位宽，平衡数据精度与资源消耗。
硬件部署
- 时序约束优化：在硬件部署时，通过Vivado的时序约束工具添加更详细的时序约束（如建立时间、保持时间约束），确保硬件运行时无时序违规。
- 资源占用优化：对FFT模块与DDS核进行资源占用分析，通过调整IP核配置（如FFT运算点数、DDS核的相位累加器位宽），在满足性能需求的前提下减少FPGA资源消耗。

基于FPGA与MATLAB的超声多普勒频移解调应用