msp432 adc FFT thd

在MSP432微控制器上实现ADC采样、FFT（快速傅里叶变换）和THD（总谐波失真）计算的完整流程及注意事项如下：

1. ADC配置

采样率选择：根据被测信号最高频率（如1kHz），按奈奎斯特准则（至少2倍），建议选择更高的采样率（如10kHz）以提高精度。
分辨率：MSP432的ADC支持12位分辨率，确保参考电压（如3.3V）稳定，减少量化噪声。
触发方式：使用定时器触发ADC，保证采样间隔均匀，避免时基抖动。
多通道/单通道：根据需求选择，若多通道需配置扫描模式。

2. 数据采集与预处理

缓存设置：开辟数组存储ADC采样值（如1024点），建议使用Ping-Pong缓冲避免数据覆盖。
去直流偏移：对采样数据减去平均值，消除直流分量影响。
窗函数处理：应用汉宁窗（Hanning）、海明窗（Hamming）等，减少频谱泄漏。

3. FFT实现

库选择：推荐使用CMSIS-DSP库（ARM官方优化），MSP432的Cortex-M4F内核支持浮点加速。
点数选择：常用256/512/1024点，点数越多频率分辨率越高（Δf=采样率/点数）。
计算步骤：
1. 将ADC数据转换为浮点数。
2. 应用窗函数。
3. 调用arm_rfft_fast_f32()等函数执行FFT。
4. 计算频域幅值（取复数模值：sqrt(re² + im²)）。

4. THD计算

基波定位：在FFT结果中找到基波频率对应的频点（如50Hz对应索引k = 50 * N / 采样率）。
谐波提取：计算2次、3次等谐波幅值（如100Hz、150Hz）。
有效值计算：
- 基波有效值：V1_rms = (基波幅值) / (窗函数修正因子)
- 谐波总有效值：Vh_rms = sqrt(Σ(谐波幅值²)) / 修正因子
- THD公式：THD (%) = (Vh_rms / V1_rms) * 100%
注：修正因子与窗函数相关（如汉宁窗为0.5）。

5. 优化与调试技巧

实时性优化：使用定点FFT（arm_cfft_q15）或降低点数以节省时间。
噪声抑制：硬件上增加RC滤波，软件上可多次采样取平均。
校准验证：用纯净正弦波信号源测试，验证THD计算结果是否合理。
频谱泄漏处理：若基波不在整周期采样，使用插值算法（如双线性插值）提高频率定位精度。

示例代码片段（基于CMSIS-DSP）

#include "arm_math.h"

#define FFT_SIZE 1024
float32_t adcBuffer[FFT_SIZE];
float32_t fftOutput[FFT_SIZE];

void computeTHD() {
    arm_rfft_fast_instance_f32 fft;
    arm_rfft_fast_init_f32(&fft, FFT_SIZE);

    // 应用汉宁窗
    for (int i=0; i<FFT_SIZE; i++) {
        adcBuffer[i] *= 0.5 * (1 - cosf(2*PI*i/(FFT_SIZE-1)));
    }

    // 执行FFT
    arm_rfft_fast_f32(&fft, adcBuffer, fftOutput, 0);

    // 计算幅值
    float32_t mag[FFT_SIZE/2];
    mag[0] = fftOutput[0]; // DC分量
    for (int i=1; i<FFT_SIZE/2; i++) {
        mag[i] = sqrtf(fftOutput[2*i]*fftOutput[2*i] + fftOutput[2*i+1]*fftOutput[2*i+1]);
    }

    // 定位基波和谐波，计算THD...
}