如何用ESP32配合MAX98357A生成无失真正弦波？

我正尝试通过I2S驱动MAX98357A生成正弦波，目前已能输出正弦波，但波形存在失真（如图所示）。

以下是我的实现代码：

#include <Arduino.h>
#include <I2S.h>

#define I2S_SAMPLE_RATE     (44100)
#define I2S_SAMPLE_BITS     16
#define PIN_LRC             8
#define PIN_BCLK            7
#define PIN_DOUT            25

int count = 0;

void setup() 
{
  I2S.setSckPin(PIN_BCLK);
  I2S.setFsPin(PIN_LRC);
  I2S.setDataPin(PIN_DOUT);
  if (!I2S.begin(I2S_PHILIPS_MODE, I2S_SAMPLE_RATE, I2S_SAMPLE_BITS))
  {
    Serial.println("Failed to initialize I2S!");
    while (1); // do nothing
  }
}

int16_t GenerateSineWave()
{
    double rad = 2 * M_PI * 1000 * count++ / I2S_SAMPLE_RATE;
    int16_t sineVal = 32767 * sin(rad);
    return sineVal;
}

void loop()
{
  I2S.write(GenerateSineWave());
  I2S.write(GenerateSineWave());
}

请问有什么方法可以生成更完美的正弦波？

优化方案 #

• 修复相位累加器溢出问题：当前count为int类型，持续累加会溢出，导致相位计算混乱，进而引发波形失真。改用无符号类型并添加模运算限制范围：

  uint32_t count = 0;
  // 1000Hz正弦波的周期采样数，避免浮点运算
  const uint32_t phaseCycle = I2S_SAMPLE_RATE / 1000;
  
  int16_t GenerateSineWave()
  {
      double rad = 2 * M_PI * count++ / phaseCycle;
      // 重置计数器，防止溢出
      if (count >= phaseCycle) count = 0;
      int16_t sineVal = static_cast<int16_t>(32767 * 0.8 * sin(rad));
      return sineVal;
  }
  

• 预生成正弦波样本，批量输出：实时调用sin函数和单次写入I2S会引入时序波动，预生成一个周期的正弦波样本，循环批量输出能大幅提升稳定性：

  const int SINE_CYCLE_SAMPLES = 44; // 44100/1000≈44.1，取整为44
  int16_t sineBuffer[SINE_CYCLE_SAMPLES * 2]; // 左右声道各一份
  
  void setup() 
  {
    I2S.setSckPin(PIN_BCLK);
    I2S.setFsPin(PIN_LRC);
    I2S.setDataPin(PIN_DOUT);
    if (!I2S.begin(I2S_PHILIPS_MODE, I2S_SAMPLE_RATE, I2S_SAMPLE_BITS))
    {
      Serial.println("Failed to initialize I2S!");
      while (1);
    }
  
    // 预生成正弦波样本
    for (int i = 0; i < SINE_CYCLE_SAMPLES; i++) {
        double rad = 2 * M_PI * i / SINE_CYCLE_SAMPLES;
        int16_t val = static_cast<int16_t>(32767 * 0.8 * sin(rad));
        sineBuffer[i*2] = val;    // 左声道数据
        sineBuffer[i*2+1] = val;  // 右声道数据
    }
  }
  
  void loop()
  {
    I2S.write(sineBuffer, sizeof(sineBuffer));
  }
  

• 使用查找表替代实时浮点运算：sin函数的实时浮点计算存在精度误差和延迟，预计算正弦值存入数组，直接读取可消除这类失真：

  // 预计算256点正弦查找表（可根据需求调整点数）
  const uint16_t SINE_LUT_SIZE = 256;
  int16_t sineLUT[SINE_LUT_SIZE];
  
  void setup()
  {
    // ... 原有初始化代码
    // 生成查找表
    for (int i = 0; i < SINE_LUT_SIZE; i++) {
        double rad = 2 * M_PI * i / SINE_LUT_SIZE;
        sineLUT[i] = static_cast<int16_t>(32767 * 0.8 * sin(rad));
    }
  }
  
  int16_t GenerateSineWave()
  {
      static uint8_t lutIndex = 0;
      return sineLUT[lutIndex++];
  }
  

• 调整输出幅度避免削波：当前代码使用满幅32767，可能触发芯片的削波失真，将幅度乘以0.8左右，留足余量可改善波形。
• 验证I2S时序匹配：确认MAX98357A的BCLK、LRC信号时序与Arduino输出一致，部分平台的I2S采样率存在微小误差，可尝试切换为48000Hz等标准采样率测试。

内容的提问来源于stack exchange，提问作者Arun K