基于脉冲序列实时生成方波音频的问题排查

核心问题诊断 #

1. 未初始化的音频缓冲区：原回调函数仅在data->play为1时写入1.0f，剩余帧未赋值，缓冲区残留垃圾数据，直接导致爆音或无规律噪音。
2. 线程竞争条件：主线程直接修改data.play，回调线程异步读取，无同步机制，状态值不稳定，无法持续触发方波输出。
3. 方波逻辑错误：原方案试图通过单次信号触发短方波段，而非根据信号标记翻转方波电平，完全不符合"每个信号对应半周期"的需求。
4. Underrun成因：nanosleep精度有限，主线程的信号模拟无法精准匹配PortAudio的回调节奏，导致音频缓冲区填充不及时。

修正方案 #

1. 完善线程同步机制

用互斥锁保护共享状态变量，避免主线程与回调线程的竞争冲突。

2. 重构方波生成逻辑

通过toggle标记记录是否需要翻转电平，每次收到黑盒信号时设置该标记，回调函数中根据当前电平状态填充整个缓冲区，并在标记触发时完成电平翻转。

3. 确保缓冲区完整填充

回调函数必须为每个音频帧赋值，避免垃圾数据，同时明确静音状态（输出0.0f）。

修改后的代码 #

#include <stdio.h>
#include <stdlib.h>
#include <portaudio.h>
#include <time.h>
#include <pthread.h>

void e(PaError err)
{
    if (err != paNoError)
    {
        Pa_Terminate();
        fprintf(stderr, "An error occurred while using the portaudio stream\n");
        fprintf(stderr, "Error number: %d\n", err);
        fprintf(stderr, "Error message: %s\n", Pa_GetErrorText(err));
        exit(1);
    }
}

typedef struct
{
    float currentLevel;  // 当前方波电平：1.0f/-1.0f（0.0f为静音）
    int toggle;          // 电平翻转标记
    pthread_mutex_t mutex; // 线程同步锁
} paData;

PaStreamParameters outputParameters;
static PaStream *stream;
static paData data;

static int paCallback(const void *inputBuffer, void *outputBuffer,
                      unsigned long framesPerBuffer,
                      const PaStreamCallbackTimeInfo *timeInfo,
                      PaStreamCallbackFlags statusFlags,
                      void *userData)
{
    paData *data = (paData *)userData;
    float *out = (float *)outputBuffer;
    (void)inputBuffer;
    (void)timeInfo;
    (void)statusFlags;

    // 加锁读取并更新状态
    pthread_mutex_lock(&data->mutex);
    float level = data->currentLevel;
    if (data->toggle) {
        level = -level;
        data->currentLevel = level;
        data->toggle = 0;
    }
    pthread_mutex_unlock(&data->mutex);

    // 填充整个音频缓冲区
    for (unsigned long i = 0; i < framesPerBuffer; i++)
    {
        *out++ = level;
    }

    return paContinue;
}

int main()
{
    printf("Initializing buzzer\n");

    // 初始化数据结构
    data.currentLevel = 0.0f;
    data.toggle = 0;
    pthread_mutex_init(&data.mutex, NULL);

    e(Pa_Initialize());
    outputParameters.device = Pa_GetDefaultOutputDevice();
    if (outputParameters.device == paNoDevice)
    {
        fprintf(stderr, "Error: No default output device.\n");
        return 1;
    }
    outputParameters.channelCount = 1;
    outputParameters.sampleFormat = paFloat32;
    outputParameters.suggestedLatency = Pa_GetDeviceInfo(outputParameters.device)->defaultLowOutputLatency;
    outputParameters.hostApiSpecificStreamInfo = NULL;
    
    // 增大缓冲区降低underrun概率
    e(Pa_OpenStream(&stream, NULL, &outputParameters, 44100, 512, paClipOff, paCallback, &data));
    e(Pa_StartStream(stream));

    struct timespec interval;
    interval.tv_sec = 0;
    interval.tv_nsec = 2500000; // 2.5ms间隔，对应400次/秒信号，生成200Hz方波

    struct timespec pause;
    pause.tv_sec = 2;
    pause.tv_nsec = 0;

    // 设置初始电平
    pthread_mutex_lock(&data.mutex);
    data.currentLevel = 1.0f;
    pthread_mutex_unlock(&data.mutex);

    // 模拟黑盒信号触发
    for (int i = 0; i < 1000; i++)
    {
        pthread_mutex_lock(&data.mutex);
        data.toggle = 1; // 标记需要翻转电平
        pthread_mutex_unlock(&data.mutex);
        nanosleep(&interval, NULL);
    }

    // 结束后静音
    pthread_mutex_lock(&data.mutex);
    data.currentLevel = 0.0f;
    pthread_mutex_unlock(&data.mutex);

    nanosleep(&pause, NULL);

    e(Pa_StopStream(stream));
    e(Pa_CloseStream(stream));
    Pa_Terminate();
    pthread_mutex_destroy(&data.mutex);

    return 0;
}

关键修改说明 #

• 线程同步：用pthread_mutex_t保护共享状态，避免主线程与回调线程的读写冲突。
• 电平翻转逻辑：每次信号触发时标记翻转，回调中完成电平切换，严格对应"每个信号为方波半周期"的需求。
• 缓冲区处理：确保每个音频帧都被赋值，彻底解决垃圾数据导致的爆音问题。
• 参数调整：将缓冲区帧大小从256改为512，降低音频输出underrun的概率（可根据设备性能微调）。

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