嵌入式项目无相位跳变可重复正弦波采样生成方案咨询

你的方案核心问题在于：

1. 单周期采样数为采样率/频率，多数情况下非整数，截断后导致周期不完整，重复播放时出现相位跳变；
2. 单周期采样数随频率大幅波动（比如200Hz对应40个采样、2000Hz对应4个采样），播放缓冲块大小变化引发感知上的音调/音量异常。

以下提供两种无需在频率不变时重复生成采样的解决思路，优先推荐第二种，更适配嵌入式场景。

解决方案一：预生成整数倍周期的采样块 #

针对每个频率，生成包含整数个完整正弦周期的采样块，确保重复播放时相位连续；同时通过最大公约数计算，让块长度尽可能规律，减少感知差异。

原理 #

对目标频率f，计算采样率SAMPLE_RATE与f的最大公约数g：

• 块内包含的完整周期数为f/g，总采样数为SAMPLE_RATE/g（必为整数）；
• 该块播放完毕后相位恰好回到初始值，重复播放无跳变。

Rust代码示例 #

use std::f32::consts::TAU;
use num_integer::gcd; // 需要在Cargo.toml添加num_integer依赖

const SAMPLE_RATE: u32 = 8000;

struct ToneGenerator {
    phase: f32,
}

impl ToneGenerator {
    fn generate_integer_cycle_block(&mut self, frequency: u32, volume_db: f32) -> Vec<f32> {
        let g = gcd(SAMPLE_RATE, frequency);
        let total_samples = SAMPLE_RATE / g;
        let delta = frequency as f32 * TAU / SAMPLE_RATE as f32;
        let amplitude = libm::powf(10.0, volume_db / 20.).min(1.0);
        
        let mut block = Vec::with_capacity(total_samples as usize);
        self.phase = 0.0; // 重置相位，保证块从0相位开始
        
        for _ in 0..total_samples {
            block.push(self.phase.sin() * amplitude);
            self.phase += delta;
            if self.phase >= TAU {
                self.phase -= TAU;
            }
        }
        
        block
    }
}

// 使用示例
fn main() {
    let mut generator = ToneGenerator { phase: 0.0 };
    // 495Hz的采样块，重复播放无相位跳变
    let block_495 = generator.generate_integer_cycle_block(495, 0.0);
    // 切换频率时仅重新生成一次
    let block_496 = generator.generate_integer_cycle_block(496, 0.0);
}

优缺点 #

• 优点：重复播放无相位跳变，仅在频率切换时生成采样块；
• 缺点：块长度仍随频率变化，可能残留感知差异，需额外计算最大公约数。

解决方案二：正弦查找表+相位累加器（推荐） #

放弃预生成单周期采样块，改用固定长度的正弦查找表+实时相位累加的方式生成采样。切换频率时仅更新相位增量，频率不变时直接查表，彻底解决两个核心问题。

原理 #

1. 程序启动时预生成一个覆盖0~TAU完整周期的正弦查找表（建议取2的幂次长度，如4096，方便取模运算）；
2. 播放时维护浮点型相位累加器，每次采样累加delta = 频率 * TAU / 采样率；
3. 将相位转换为查找表索引，通过线性插值获取高精度采样值；
4. 频率切换时仅更新delta，无需重新生成任何采样数据。

Rust代码示例 #

use std::f32::consts::TAU;

const SAMPLE_RATE: u32 = 8000;
const SIN_TABLE_SIZE: usize = 4096; // 2的幂次，位运算替代取模提升效率
const BUFFER_SIZE: usize = 100; // 固定播放缓冲大小

// 预生成正弦查找表，仅初始化一次
fn build_sin_table() -> [f32; SIN_TABLE_SIZE] {
    let mut table = [0.0; SIN_TABLE_SIZE];
    for i in 0..SIN_TABLE_SIZE {
        let phase = i as f32 * TAU / SIN_TABLE_SIZE as f32;
        table[i] = phase.sin();
    }
    table
}

struct TonePlayer {
    sin_table: [f32; SIN_TABLE_SIZE],
    phase: f32,
    delta: f32,
    amplitude: f32,
}

impl TonePlayer {
    fn new(sin_table: [f32; SIN_TABLE_SIZE]) -> Self {
        Self {
            sin_table,
            phase: 0.0,
            delta: 0.0,
            amplitude: 1.0,
        }
    }

    // 切换频率/音量时调用，仅更新参数
    fn set_tone(&mut self, frequency: u32, volume_db: f32) {
        self.delta = frequency as f32 * TAU / SAMPLE_RATE as f32;
        self.amplitude = libm::powf(10.0, volume_db / 20.).min(1.0);
    }

    // 填充固定大小的播放缓冲，频率不变时重复调用
    fn fill_buffer(&mut self, buffer: &mut [f32; BUFFER_SIZE]) {
        for value in buffer.iter_mut() {
            // 计算查找表索引（含小数部分）
            let table_pos = self.phase / TAU * SIN_TABLE_SIZE as f32;
            let idx = table_pos.floor() as usize;
            let frac = table_pos - idx as f32;

            // 线性插值提升精度，性能不足时可省略
            let val1 = self.sin_table[idx];
            let val2 = self.sin_table[(idx + 1) & (SIN_TABLE_SIZE - 1)]; // 位运算替代取模
            *value = (val1 * (1.0 - frac) + val2 * frac) * self.amplitude;

            // 更新相位，保持连续
            self.phase += self.delta;
            if self.phase >= TAU {
                self.phase -= TAU;
            }
        }
    }
}

// 使用示例
fn main() {
    let sin_table = build_sin_table();
    let mut player = TonePlayer::new(sin_table);
    
    // 设置初始频率
    player.set_tone(495, 0.0);
    let mut buffer = [0.0; BUFFER_SIZE];
    // 填充缓冲，重复调用即可播放
    player.fill_buffer(&mut buffer);
    
    // 切换频率时仅更新参数
    player.set_tone(496, 0.0);
    player.fill_buffer(&mut buffer);
}

优缺点 #

• 优点：
  • 彻底解决相位跳变问题，相位累加连续无截断；
  • 缓冲大小固定，频率切换时感知一致；
  • 仅需初始化一次正弦表，切换频率时仅更新参数，生成耗时极低；
  • 线性插值保证采样精度，省略插值也能满足多数音频场景需求。
• 缺点：需要少量实时计算，但现代MCU完全可胜任。

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