三亩地 三亩地SAN MU DI · CODE DIARY
交流投稿
ARTICLE DETAIL

日记详情

真实记录编程学习的某一天，欢迎挑你感兴趣的翻一翻。

首页/学习日记/日记详情

基于STM32F4的实时FFT与DDS信号处理系统设计

发布时间：2026/8/6 5:23:55 三亩地 编程学习日记
基于STM32F4的实时FFT与DDS信号处理系统设计

1. 项目概述与核心思路

“23年电赛H题方案——（省钱版）”，这个标题一出来，老电赛人估计都会心一笑。每年电赛，H题往往意味着挑战性，而“省钱版”三个字，更是戳中了无数参赛队伍，尤其是学生团队的痛点。官方推荐的方案动辄需要高性能DSP、高速ADC、甚至FPGA，硬件成本直接起飞，让本就紧张的备赛预算雪上加霜。我这个方案，核心思路就是用最“平民”的硬件——STM32F4系列单片机，去实现那些看似需要高端平台才能完成的任务，比如实时FFT频谱分析和DDS信号生成。

为什么是STM32F4？因为它是在性能、外设丰富度和价格之间取得绝佳平衡点的“神U”。以STM32F407或F429为例，它们自带硬件浮点运算单元（FPU），主频能跑到168MHz甚至180MHz，这对于运行FFT算法至关重要。同时，它们通常集成了多个ADC和DAC，以及强大的定时器系统，完全能满足H题中对信号采集、处理和生成的基本要求。用一百多块钱的核心板，去挑战上千元专业板卡才能流畅处理的问题，这就是“省钱版”的精髓。这个方案不仅适用于回顾23年的赛题，其核心的“用通用MCU实现专用功能”的思路，对于准备未来电赛，如26年可能出现的类似题目（从热词看，大家对26年H题已有期待），也具有很强的借鉴意义。它适合那些希望深入理解数字信号处理（DSP）底层原理，而不想被昂贵硬件束缚的嵌入式开发者和学生。

2. 硬件平台选型与最小系统搭建

2.1 MCU选型：为什么是STM32F4？

在“省钱”的前提下，选型的第一原则是：用最少的钱，买最多的算力和最合适的外设。STM32F103（Cortex-M3）虽然便宜，但缺少硬件FPU，做浮点FFT会非常吃力，难以满足实时性要求。STM32H7系列性能爆炸，但价格也上去了，不符合“省钱”主题。因此，STM32F4系列成为了黄金选择。

具体型号上，STM32F407VET6和STM32F429VIT6是两款经典型号。F407主频168MHz，带有FPU，价格亲民，是性价比之王。F429主频180MHz，除了性能稍强，最大的优势是集成了LCD-TFT控制器和Chrom-ART加速器，如果你的题目涉及图形显示（比如绘制频谱图），F429会更有优势，但成本也相应增加二三十元。对于纯信号处理，F407完全足够。购买时，直接选择核心板是最省事省钱的方案，一块F407核心板价格通常在百元以内，集成了基本的电源、晶振和调试接口，我们只需要专注于功能扩展。

2.2 关键外设与电路设计要点

电赛H题通常围绕信号展开，因此ADC和DAC是核心外设。

1. 信号采集（ADC）通道：STM32F407通常有3个ADC，每个ADC最多有16个外部通道。为了获得较高的采样率，我们需要使用ADC的规则组，并配合DMA进行连续搬运，避免CPU频繁中断。这里有个关键技巧：使用定时器触发ADC采样。例如，使用TIM2的更新事件来触发ADC1的采样，这样就能获得一个极其精确、稳定的采样频率（Fs），这是后续进行FFT分析的基础。采样率设置多高？这需要根据题目要求的信号最高频率（Fmax）来定，遵循奈奎斯特采样定理（Fs > 2*Fmax），通常取Fs = (2.5 ~ 5) * Fmax。如果题目信号频率在1kHz以内，设置个10kSPS（每秒采样点数）的采样率就绰绰有余了。

2. 信号生成（DAC）通道：STM32F407自带两个12位DAC。我们用其中一个来实现DDS（直接数字频率合成）。DDS的原理后面会细说，这里先讲硬件连接。DAC输出后，通常需要接一个简单的运放电路进行缓冲和滤波。因为DAC输出的是阶梯波，含有高频谐波。一个典型的电路是：DAC输出 → 电压跟随器（提高带载能力）→ 低通滤波器（滤除奈奎斯特频率以上的杂散）。运放可以选择非常便宜的通用型运放，如LM358，完全够用。

3. 时钟与电源：核心板一般已集成8MHz晶振和32.768kHz RTC晶振。我们需要关注的是，在CubeMX配置中，是否正确配置PLL，将系统时钟提升到168MHz，并确保ADC、定时器的时钟源正确。电源部分，核心板的3.3V输出可能驱动能力有限，如果外接模块较多（如屏幕、多个运放），建议从12V或5V输入口单独引一路电源，使用AMS1117-3.3等LDO芯片为数字部分和模拟部分分别供电，并在模拟电源入口加磁珠和电容滤波，减少数字噪声对模拟信号的干扰。

注意：模拟地（AGND）和数字地（DGND）建议在一点共地，通常是电源输入滤波电容的接地端。这是降低系统噪声、提高ADC采样精度的关键一步，很多新手会忽略这一点，导致采样值跳动大。

3. 软件架构与核心算法实现

3.1 基于CMSIS-DSP的FFT实现与优化

FFT是频谱分析的核心。自己写FFT算法不仅容易出错，而且效率低下。STM32的CubeMX软件包中，包含了ARM官方优化的CMSIS-DSP库，这里面就有高度优化的FFT函数，我们必须利用起来。

首先，在CubeMX中安装软件包时，记得勾选“CMSIS DSP”。然后在代码中引入头文件#include “arm_math.h”和#include “arm_const_structs.h”。CMSIS-DSP库提供了多种FFT函数，对于实数序列（我们的ADC采样值就是实数），使用arm_rfft_fast_f32()是最佳选择。它内部会调用复数FFT，但接口更简单高效。

以下是关键步骤：

  1. 数据准备：定义一个浮点数组float32_t adc_buffer[FFT_LENGTH]，用于存放ADC采样值。FFT_LENGTH必须是2的整数次幂，如1024、2048。长度越长，频率分辨率越高（Δf = Fs / N），但计算量也越大。需要根据实时性要求权衡。
  2. 窗函数应用：直接对采样数据进行FFT会发生“频谱泄漏”。为了解决这个问题，需要对数据加窗。常用的有汉宁窗（Hanning）、汉明窗（Hamming）。CMSIS-DSP库也提供了窗函数生成函数，如arm_hanning_f32()。将窗数组与采样数据数组逐点相乘，完成加窗。
  3. 执行FFT：调用arm_rfft_fast_f32(&S, adc_buffer, fft_output, 0)。其中S是初始化好的FFT实例结构体，fft_output是输出数组。注意，这个输出数组的前半部分（N/2+1个点）包含了从0到Fs/2的频率分量对应的复数结果（实部+虚部）。
  4. 计算幅值：FFT输出是复数，我们需要计算每个频率点的幅值来画频谱。调用arm_cmplx_mag_f32(fft_output, magnitude, FFT_LENGTH/2)来计算模值。得到的magnitude数组就是频谱的纵坐标。
  5. 频率刻度：横坐标是频率。第k个点对应的频率是f = k * Fs / FFT_LENGTH。其中k从0到N/2。

优化心得：

  • 将FFT相关的数组（输入、输出、窗、幅值）定义在CCM RAM（如果芯片有）或者使用__attribute__((section(“.dtcm”)))指定到高速RAM中，可以显著提升计算速度，因为FFT是内存访问密集型操作。
  • 如果采样率Fs固定，频率刻度可以预先计算好，存为一个常量数组，避免在每次FFT后都进行浮点乘除运算。
  • 对于固定点数的FFT，arm_rfft_fast_init_f32(&S, FFT_LENGTH)初始化函数只需要执行一次，放在主循环外。

3.2 定时器触发ADC与DMA双缓冲机制

稳定的数据流是实时处理的前提。我们采用“定时器触发ADC + DMA双缓冲”的经典架构。

  1. CubeMX配置：

    • 配置一个定时器（如TIM2）为PWM输出模式或输出比较模式，但其输出引脚不用接任何东西。我们只利用它的“更新事件”。将定时器的周期值设置为SystemCoreClock / (PSC+1) / Fs - 1，这样它就能以Fs的频率产生更新中断（或触发输出）。
    • 配置ADC1，选择“Regular Conversion”，触发源选择“Timer 2 Trigger Out event”。采样时间根据信号源阻抗设置，一般1.5到239.5个周期，时间越长，采样精度越高，但最高采样率会下降。
    • 启用DMA，模式设为“Circular”（循环模式），数据宽度为半字（对应ADC的12位数据）。内存地址递增。

  2. 双缓冲实现： 单DMA循环缓冲区只能保证数据连续，但处理数据时如果DMA正好写到了处理区域，就会导致数据错乱。双缓冲解决了这个问题。我们定义两个大小均为FFT_LENGTH的数组：BufferA和BufferB。

    • 初始时，DMA指向BufferA。
    • 当DMA搬运完成BufferA（通过DMA半传输完成或传输完成中断判断），我们让一个标志位buf_ready = 1，并记录当前准备好的缓冲区是A。
    • 在主循环或一个低优先级任务中，不断检查buf_ready。一旦为1，就启动FFT处理流程（处理BufferA），同时将DMA的目标地址切换到BufferB。
    • 当BufferB也满了，标志位置位，处理BufferB，DMA切回BufferA，如此循环。 这样，数据采集（DMA）和数据处理（CPU）就实现了流水线并行，互不干扰，极大提高了系统效率和实时性。

3.3 DDS信号生成原理与软件实现

DDS是一种用于生成任意波形的高效数字技术。在STM32上，我们可以用定时器触发DAC，通过查表法来实现。

  1. 原理：DDS的核心是一个相位累加器。我们有一个大的波形查找表（LUT），比如存储了一个正弦波一个周期的1024个采样点。还有一个变量叫“相位累加器”（一个32位或64位的整数），每个时钟周期（定时器触发周期），相位累加器就增加一个固定的值，叫做“频率控制字”（FTW）。

    • 输出频率Fout = (FTW * Fclk) / (2^N)。其中，Fclk是定时器触发DAC的频率（即“更新率”），N是相位累加器的位宽（比如32），2^N就是查找表的长度（如果使用高位索引）。
    • 我们通常用相位累加器的高位（比如高10位）作为索引，去查找表里取出对应的幅度值，送给DAC输出。累加器低位的作用是提高频率分辨率。

  2. STM32实现步骤：

    • 生成查找表：在PC上用Python或MATLAB生成一个正弦波数组（浮点数或整数），然后量化为DAC的输入范围（0-4095对应0-3.3V）。将数组定义为const常量，存储在Flash中。
    • 配置定时器：另一个定时器（如TIM6，基本定时器）配置为以所需的Fclk频率产生更新事件。Fclk决定了DDS输出的最高频率（奈奎斯特极限）和波形阶梯的平滑度，通常至少是目标输出频率的20倍以上。
    • 配置DAC：DAC设置为定时器触发，输出缓冲区禁用（如果驱动外部运放）。
    • 编写中断服务程序：在定时器更新中断中，进行相位累加：phase_accumulator += FTW;。然后，用phase_accumulator的高位（例如>> 22来取高10位）作为索引，从查找表中读取数据，写入DAC的数据保持寄存器。

代码片段示例（概念性）：

#define LUT_SIZE 1024 #define PHASE_BITS 32 const uint16_t sine_lut[LUT_SIZE] = { ... }; // 正弦波查找表 volatile uint32_t phase_accumulator = 0; uint32_t frequency_tuning_word = 0; // 需要根据目标频率计算 void TIM6_DAC_IRQHandler(void) { if(TIM6->SR & TIM_SR_UIF) { TIM6->SR &= ~TIM_SR_UIF; // 清除中断标志 phase_accumulator += frequency_tuning_word; uint16_t index = (phase_accumulator >> 22) & 0x3FF; // 取高10位作为索引 DAC->DHR12R1 = sine_lut[index]; // 写入DAC通道1 } }

通过改变frequency_tuning_word，就可以实时、平滑地改变输出信号的频率，这是模拟PLL电路难以做到的。

4. 系统集成与调试实战

4.1 多任务调度与实时性保障

一个完整的电赛系统，通常需要同时完成信号采集、FFT计算、结果显示、按键响应、DDS频率控制等多个任务。如何让它们在小小的单片机上和谐共处？这里不建议上重量级的RTOS，对于F4来说，用一个简单的“前后台+状态机”模型或者轻量级调度器足以应对，而且更节省资源，确定性更强。

我常用的方法是基于SysTick中断的“时间片轮询”调度。将SysTick中断设置为1ms一次。在中断服务程序中，不进行复杂操作，只对一系列标志位进行累加计数。例如：

volatile uint32_t g_ticks_10ms = 0; volatile uint32_t g_ticks_100ms = 0; void SysTick_Handler(void) { static uint8_t cnt_10 = 0， cnt_100 = 0; if(++cnt_10 >= 10) { g_ticks_10ms++; cnt_10 = 0; } if(++cnt_100 >= 100) { g_ticks_100ms++; cnt_100 = 0; } }

在主循环中，通过检查这些全局的标志位来执行不同周期的任务：

while(1) { // 任务1：高频任务，如检查DMA缓冲是否就绪，就绪则处理FFT（非阻塞式） if(adc_buf_ready_flag) { process_fft(); adc_buf_ready_flag = 0; } // 任务2：10ms任务，如扫描按键、更新部分显示 static uint32_t last_10ms = 0; if(g_ticks_10ms != last_10ms) { last_10ms = g_ticks_10ms; key_scan(); update_some_display(); } // 任务3：100ms任务，如刷新整个频谱图界面、串口发送数据 static uint32_t last_100ms = 0; if(g_ticks_100ms != last_100ms) { last_100ms = g_ticks_100ms; refresh_spectrum_display(); uart_send_data(); } // 空闲时可以进入低功耗模式 __WFI(); }

这种架构清晰、响应及时，且能保证FFT这种计算密集型任务得到及时处理。关键是要确保process_fft()函数的执行时间远小于ADC缓冲区填满的时间（即FFT_LENGTH / Fs），否则系统会崩溃。

4.2 频谱显示与交互设计

显示部分，如果为了极致省钱，可以用串口打印到电脑，用PC软件（如SerialPlot、MATLAB）绘图。但电赛通常要求独立运行，所以一块廉价的SPI接口OLED屏（如0.96寸 SSD1306）或TFT LCD屏（如ILI9341）是更好的选择。

在屏幕上绘制频谱图，本质上是将计算得到的magnitude数组进行可视化。有几个要点：

  1. 坐标映射：需要将频率数组（0 ~ Fs/2）映射到屏幕的X轴（0 ~ 屏幕宽度-1），将幅值数组映射到屏幕的Y轴。幅值通常取对数（20*log10）转换成分贝（dB）显示，这样能更好地观察动态范围。
  2. 绘图优化：避免在循环中频繁调用单点绘制函数，那样太慢。可以先将所有点计算好，然后使用画线函数（ILI9341_DrawLine）或一次性填充缓冲区（对于OLED）的方式进行绘制。对于FFT_LENGTH=1024，屏幕宽度可能只有128或320，需要进行降采样显示，比如每8个点取一个最大值显示。
  3. 交互控制：通过旋转编码器或按键来调整DDS的输出频率、幅度，或者调整FFT分析的参数（如采样率、窗函数类型）。编码器的读数处理建议使用定时器编码器接口模式，非常稳定。按键则采用状态机进行消抖和长短按识别。

4.3 系统校准与性能测试

系统搭建好后，必须进行校准，否则测量结果毫无意义。

  1. ADC直流偏移校准：短接ADC输入引脚到地（GND），启动ADC连续采样一段时间（比如采集1024个点），计算其平均值，这个值就是“零点偏移”。在后续所有采样值中，都需要减去这个偏移量。
  2. ADC幅度校准：输入一个已知幅度和频率的纯净正弦波（可以用信号发生器或另一个已校准的DDS通道产生），测量其通过FFT计算得到的主频幅值。与理论值对比，可以得到一个幅度校正系数。注意，这个系数可能与频率有关，所以最好在多个频点进行校准，存入一个校正表。
  3. DAC输出校准：用万用表测量DAC在不同输出码值下的实际电压，与理论值对比。由于运放偏移、电阻精度等问题，可能存在增益和偏移误差。可以在软件中建立一个简单的线性校正公式：实际输出码值 = a * 期望码值 + b，通过测量两个点解出a和b。
  4. 性能测试：
    • 最高实时采样率：逐步提高ADC的采样定时器频率，直到DMA搬运开始出现丢数（通过检查DMA缓冲区数据是否连续可判断）。这个频率就是系统能稳定工作的最高采样率。
    • FFT计算耗时：在process_fft()函数开始和结束处读取SysTick的计数器值，计算差值。这个时间必须小于缓冲区更新周期。
    • 信噪比与动态范围测试：输入一个单频正弦波，观察频谱图上主频的幅值与其他频点（噪声基底）的幅值之差，大致就是系统的动态范围。好的设计应该在70dB以上。

5. 常见问题排查与深度优化技巧

5.1 信号采集与FFT的典型问题

问题1：频谱图上有明显的杂散峰（毛刺），不是单一的信号频率线。

  • 原因A：频谱泄漏。这是最常见的原因。没有加窗，或者加窗函数选择不当。解决方案：务必对采样数据加窗（汉宁窗适用于大多数情况）。确保窗函数数组与采样数据正确逐点相乘。
  • 原因B：ADC采样时钟或触发不稳定。检查定时器配置是否正确，是否被其他高优先级中断频繁打断。确保定时器时钟源稳定（通常用APB总线时钟）。可以在定时器触发引脚上用示波器观察，看脉冲间隔是否均匀。
  • 原因C：电源噪声或接地不良。模拟部分电源纹波过大，数字地噪声串扰到了模拟信号。检查电源滤波电容是否足够（建议在模拟电源入口加10uF钽电容并联0.1uF陶瓷电容），确保模拟地和数字地单点连接。
  • 原因D：信号本身不纯净。使用信号发生器输入一个纯净正弦波进行测试，如果问题消失，说明问题在外部信号源或前端调理电路。

问题2：测量得到的频率值总是有固定偏差。

  • 原因：采样率（Fs）不准确。Fs是由定时器频率决定的。计算公式Fs = Timer_Clock / (PSC+1) / (ARR+1)。请仔细检查：
    1. 定时器的时钟源频率是多少？是否经过了倍频？在CubeMX的Clock Configuration界面确认。
    2. PSC和ARR的值计算是否正确？ARR = (Timer_Clock / (PSC+1) / Fs) - 1。
  • 验证方法：用定时器触发一个GPIO翻转，用示波器测量该GPIO的频率，它应该是Fs的一半（因为每次触发翻转一次）。用这个实测频率作为Fs代入FFT频率计算公式。

问题3：幅值测量不准，尤其是低频和高频部分。

  • 原因A：ADC采样时间不足。对于高阻抗信号源，ADC采样电容充电需要时间。增加ADC的采样周期（SMPR寄存器设置）。
  • 原因B：未进行系统幅频响应校准。ADC和前端运放电路对不同频率的信号的增益可能不同。需要对多个标准频率点进行测量，建立幅值-频率校正曲线或查找表。
  • 原因C：FFT点数不足，频率分辨率低。如果信号频率正好落在两个FFT频点之间，其能量会泄漏到相邻频点，导致主频幅值测量偏低。增加FFT点数（N）可以提高分辨率，但会增加计算量。也可以采用插值算法（如Rife算法）来估计两个频点之间的精确频率和幅值。

5.2 DDS输出信号的优化

问题1：DDS输出正弦波有台阶感，不够平滑。

  • 原因：DAC更新率（Fclk）不够高，或者未使用滤波器。根据奈奎斯特采样定理，要无失真地重建一个频率为Fout的信号，DAC更新率至少需要2*Fout。实际上，为了获得平滑波形，通常要求Fclk > 20 * Fout。提高定时器触发DAC的频率。同时，必须在DAC输出后接入一个截止频率略高于最大Fout的低通滤波器，以滤除高频量化噪声（台阶带来的谐波）。

问题2：改变频率控制字（FTW）时，输出信号有相位跳变。

  • 原因：直接给phase_accumulator赋值新的FTW。这会导致相位累加器的值不连续，从而在输出波形上产生一个相位跳变。正确的做法是只改变frequency_tuning_word变量，而保持phase_accumulator自然累加。这样，频率切换是平滑的，相位是连续的。如果需要相位重置，可以在特定时刻将phase_accumulator清零。

问题3：输出信号幅度随频率变化。

  • 原因：这是DDS的固有缺陷——sinc函数效应。由于DAC是零阶保持输出，其频谱会乘以一个sinc函数。频率越高，衰减越大。可以在软件中进行反sinc补偿：在生成查找表时，对每个点的幅值预先除以sinc(π * f / Fclk)。更简单的方法是，在输出幅值前乘以一个与当前输出频率相关的补偿系数。

5.3 资源与实时性的极限压榨

当FFT点数很大（如4096），或者采样率很高时，系统可能会卡顿。以下是一些压榨性能的技巧：

  1. 使用整数FFT：如果对精度要求不是极高，可以使用CMSIS-DSP库中的Q15或Q31格式的定点FFT（如arm_rfft_q15）。定点运算速度远快于浮点，尤其在没有FPU的芯片上。但需要处理好数据定标（Q格式），防止溢出。
  2. 启用CPU缓存和ART加速器：对于STM32F429等带有ART加速器的芯片，务必在系统初始化时启用它。将程序代码特别是FFT库代码放在Flash中，ART加速器能大幅提升取指速度。同时，启用指令缓存和数据缓存。
  3. 优化内存访问：这是最关键的一点。确保FFT的输入/输出数组、旋转因子表都放在DTCM或CCM RAM（核心耦合内存）中。这类内存与CPU同速，没有总线竞争。可以通过修改链接脚本（.ld文件）或使用__attribute__指定段来实现。
  4. 降低显示刷新率：频谱图不需要每秒更新60次。根据人眼视觉暂留，10-20Hz的刷新率已经足够流畅。将频谱显示任务周期设为50ms或100ms，能节省大量CPU时间用于信号处理。
  5. 分时复用ADC：如果题目需要多通道采集，但不需要严格同步，可以配置ADC在扫描模式下工作，用DMA搬运多通道数据。然后在对各通道数据进行FFT时，错开时间片进行处理，避免集中计算导致CPU过载。

最后，分享一个调试利器：SEGGER的SystemView。它是一个基于J-Link的实时系统可视化工具。即使不用RTOS，你也可以插入一些自定义事件，来观察ADC中断、DMA传输、FFT计算、显示刷新等任务的执行时序和耗时。它能帮你一眼看出系统瓶颈在哪里，是优化实时系统不可或缺的“上帝视角”。虽然J-Link比ST-Link贵，但对于深入调试复杂嵌入式系统来说，这笔投资是值得的。

编程学习 学习日记 实战经验
← 返回列表