1. 为什么ADAU1701值得折腾60段EQ与防啸叫ADAU1701这颗SigmaDSP芯片在音频圈子里算是老面孔了28/56位的定点DSP内核内置ADC和DAC支持I2S/TDM输入输出最大能跑到约50 MIPS的运算能力。很多人第一次接触它是在廉价的DSP功放板或者车载音频处理模块上价格便宜、开发工具SigmaStudio又是图形化拖拽看起来门槛很低。但真要把60段参量EQ加上一套能用的防啸叫算法塞进去你会发现资源根本不够用——这不是危言耸听我自己第一次尝试的时候编译直接爆了程序空间连下载都下不进去。先说清楚这个项目的定位它适合已经玩过ADAU1701基础例程、能用SigmaStudio拖出简单分频和几段EQ的玩家也适合做KTV前级、会议麦克风阵列、便携扩声设备的工程师。核心目标是在一颗ADAU1701上同时实现60段参量均衡和自动防啸叫Notch Filter自动扫描并且保证实时性不崩。为什么是60段因为专业扩声领域里31段图形EQ是标配双通道就是62段而参量EQ的灵活性更高60段可以覆盖20 Hz到20 kHz的精细修整尤其是针对房间驻波和话筒频响补偿。防啸叫则是现场扩声的刚需没有它话筒一开就尖叫再好的EQ也白搭。但ADAU1701的运算资源是硬约束。SigmaStudio里每个二阶IIR滤波器Biquad大约消耗5个指令周期单精度或更多60段EQ如果全部用双精度Biquad光滤波器就要吃掉300 MIPS以上远超芯片能力。所以必须做取舍用单精度、用级联结构、用采样率折中、用防啸叫算法共享滤波器资源。这些取舍背后的计算和实测就是这篇博文要讲清楚的东西。提示ADAU1701的SigmaStudio工程里程序空间Program RAM是1024条指令参数空间Parameter RAM是1024个32位字。60段EQ加上防啸叫参数空间会非常紧张必须精打细算。2. 资源预算与核心架构设计2.1 先算一笔账60段EQ到底要多少MIPSADAU1701的内核是28/56位定点SigmaStudio里常用的Biquad模块有几种实现单精度Single Precision、双精度Double Precision、以及直接型I和直接型II。直接型II更省参数RAM但数值稳定性稍差。我实测下来在48 kHz采样率下一个单精度Biquad大约消耗5个指令周期双精度大约9到10个。ADAU1701的最大指令率是50 MIPS左右具体取决于SigmaStudio版本和时钟配置所以60段单精度Biquad60 × 5 300指令周期/样本48 kHz下就是14.4 MIPS看起来还行。但如果加上防啸叫的自动扫描算法比如需要实时做FFT或者多段检测那就不够了。实际上SigmaStudio里拖出来的Biquad模块还有额外的读写参数开销每个滤波器大约要占用2到3个参数RAM字。60段就是120到180个字而ADAU1701总共1024个参数字还要留给音量控制、延时、动态处理等。所以必须用直接型II结构并且把EQ分成两组一组是固定的房间校正EQ比如20段另一组是动态的防啸叫Notch比如40段。这样参数RAM可以复用因为防啸叫的Notch滤波器系数是实时更新的不需要同时保留所有系数。我的做法是用SigmaStudio的“Single Precision IIR”模块选择直接型II采样率设成48 kHz然后通过“Parametric EQ”库里的多段EQ模块来搭建。但SigmaStudio自带的图形EQ最多好像只给到10段左右要60段必须自己用Biquad级联。级联的时候要注意每级之间不能有溢出所以要在每级后面加一个“Gain”模块做衰减通常-3 dB到-6 dB最后再统一补回来。2.2 防啸叫算法的选型为什么不用FFT防啸叫的核心是检测啸叫频率并生成窄带Notch。常见方案有三种一是基于FFT的频域检测二是基于自适应滤波的LMS算法三是基于峰值检测的扫描法。FFT在ADAU1701上做256点需要大约多少MIPS我查过SigmaStudio的FFT模块256点复数FFT大约要消耗15到20 MIPS再加上峰值搜索和Notch生成总开销会超过30 MIPS留给EQ的就不够了。而且FFT的频率分辨率在48 kHz下是187.5 Hz对于啸叫这种窄带信号来说太粗了。LMS自适应滤波更不适合因为收敛速度慢而且需要参考信号在单话筒扩声场景里没有干净的参考。所以我选的是“峰值检测扫描Notch”方案。具体来说用一个带通滤波器组或者多个Goertzel滤波器来检测特定频段的能量当某个频段的能量超过阈值并持续一段时间就认为发生了啸叫然后在该频点插入一个Notch滤波器Q值设高比如10到30深度设-12 dB到-24 dB。这个方案的计算量小因为Goertzel每个频点只需要2个乘加检测20个频点也才40个乘加远低于FFT。而且Notch滤波器的系数可以实时计算用SigmaStudio的“Lookup Table”或者“Single Precision IIR”配合参数更新来实现。注意ADAU1701不支持浮点所有系数都是定点格式。Notch滤波器的系数计算需要用到三角函数在DSP里通常用查表法或者近似公式。SigmaStudio里可以用“Sine”模块或者预计算表。2.3 整体信号链设计我的信号链是这样的输入ADC或I2S→ 输入增益 → 60段EQ其中20段固定40段动态→ 防啸叫检测与Notch插入 → 输出增益 → 限幅器 → DAC。防啸叫检测的抽头点在EQ之后、Notch之前这样检测的是经过房间校正后的信号更准确。Notch滤波器的插入位置在EQ之后因为EQ可能会改变啸叫频率的增益先EQ再Notch可以保证Notch的准确性。在SigmaStudio里我用了一个“Multiplexer”来切换固定EQ和动态Notch的参数更新。固定EQ的系数在编译时确定动态Notch的系数通过“Parameter Update”模块从外部MCU或者SigmaStudio的实时控制界面写入。如果不需要外部控制也可以用SigmaStudio内部的“Lookup Table”和“Counter”来自动扫描但那样会占用更多指令。3. 在SigmaStudio里搭建60段EQ的实操步骤3.1 工程初始化与采样率设置打开SigmaStudio新建一个工程选择ADAU1701作为目标器件。在“Hardware Configuration”里把采样率设为48 kHz因为这是音频标准而且ADAU1701在48 kHz下MIPS余量最大。如果设成96 kHzMIPS直接翻倍60段EQ就悬了。时钟源用内部PLLMCLK设为12.288 MHz这样正好是48 kHz的256倍。接下来在“Schematic”里拖入“Input”和“Output”模块。输入选择“Serial Input”或者“ADC”取决于你的硬件。如果是ADAU1701评估板通常用ADC输入。输出选“DAC”或者“Serial Output”。然后加一个“Gain”模块作为输入增益初始设0 dB。实操心得SigmaStudio的工程里所有模块的采样率必须一致否则编译会报错。如果你从别的工程复制模块记得检查采样率设置。3.2 用Biquad级联搭建60段参量EQSigmaStudio的“Parametric EQ”库里有“Single Precision IIR”模块每个模块就是一个Biquad。要搭60段就得拖60个这样的模块然后串联起来。但这样太占地方而且参数RAM会爆。我的做法是用“Multi-Biquad”模块SigmaStudio里有一个“Multiple Biquad”或者“Cascade”模块可以一次放多个Biquad但参数RAM还是按个数算。具体操作从“Filters”库里拖出“Single Precision IIR”模块双击设置参数。每个Biquad需要设置F中心频率、Q品质因数、Gain增益。对于参量EQF从20 Hz到20 kHz按对数分布Q值根据用途定房间校正用1到4话筒补偿用0.7到2陷波用10到30。Gain范围-15 dB到15 dB。但60个模块手动设置太累我写了一个Python脚本生成SigmaStudio的XML参数文件然后导入。SigmaStudio支持“Import Parameters”功能可以批量导入系数。具体格式是每个Biquad的系数a0, a1, a2, b0, b1, b2定点格式是5.23ADAU1701的28位定点实际用5.23表示。计算Biquad系数的公式是标准的RBJ Audio EQ Cookbook公式网上有现成的。import math def peaking_eq(fs, f0, Q, gain_db): A 10**(gain_db/40) w0 2*math.pi*f0/fs alpha math.sin(w0)/(2*Q) b0 1 alpha*A b1 -2*math.cos(w0) b2 1 - alpha*A a0 1 alpha/A a1 -2*math.cos(w0) a2 1 - alpha/A # 归一化并转成5.23定点 b0 / a0; b1 / a0; b2 / a0 a1 / a0; a2 / a0 return b0, b1, b2, a1, a2然后把这些系数按SigmaStudio的格式写入XML。注意ADAU1701的Biquad系数是b0, b1, b2, a1, a2a0固定为1而且符号位要处理好。我踩过的坑是SigmaStudio里的系数是5.23格式范围-16到16所以如果系数超过1要右移。实际写入时SigmaStudio会自动处理但如果你手动改XML一定要检查。3.3 参数RAM的优化技巧60段EQ如果每段都独立参数参数RAM肯定不够。我的优化方法是把EQ分成两组一组是“固定EQ”系数在编译时写死不占用参数RAM的运行时更新空间另一组是“动态EQ”系数可以通过外部MCU更新但只保留当前需要的段数。比如固定20段动态40段但动态的40段里只有10段是同时激活的其他30段是“影子”参数需要时再切换。这样参数RAM可以控制在200字以内。SigmaStudio里可以用“Parameter RAM”的“Alias”功能把多个Biquad的参数映射到同一块RAM但这样就不能同时独立控制。另一种方法是使用“Indexed Parameter”模块通过索引来选择参数集。我实测下来用“Lookup Table”加“Counter”可以实现参数的分时复用但会增加指令开销。注意ADAU1701的参数RAM是1024个32位字但SigmaStudio编译后会显示实际使用量。如果超过编译会报错“Parameter RAM overflow”。这时候要么减少段数要么用更省参数的Biquad结构。4. 防啸叫算法的实现与调试4.1 啸叫检测用Goertzel还是带通滤波器组前面说了FFT太耗资源所以我用Goertzel算法检测特定频点。Goertzel的本质是一个二阶IIR计算某个频点的能量。在SigmaStudio里可以用“Single Precision IIR”模块配置成Goertzel但需要自己算系数。Goertzel的系数是def goertzel_coeff(fs, f0): w0 2*math.pi*f0/fs coeff 2*math.cos(w0) return coeff然后每个频点需要一个Goertzel滤波器检测20个频点就是20个IIR每个IIR大约3个指令周期总共60指令周期可以接受。但Goertzel需要平方和开方来算能量开方在定点DSP里可以用查表或者牛顿迭代。SigmaStudio里有“Square Root”模块但比较耗指令。我的做法是用“Absolute”加“Log”近似或者直接用阈值比较不需要精确能量值。检测频点的选择根据话筒和房间的典型啸叫频率通常集中在1 kHz到8 kHz之间。我选了20个频点按1/3倍频程分布从800 Hz到10 kHz。每个频点的检测窗口是10 ms如果能量超过阈值比如-20 dBFS并持续3个窗口就触发Notch。4.2 Notch滤波器的动态插入触发后需要在该频点插入一个Notch。Notch的Q值设20深度-18 dB。系数计算用RBJ公式的Notch版本def notch(fs, f0, Q): w0 2*math.pi*f0/fs alpha math.sin(w0)/(2*Q) b0 1 b1 -2*math.cos(w0) b2 1 a0 1 alpha a1 -2*math.cos(w0) a2 1 - alpha b0 / a0; b1 / a0; b2 / a0 a1 / a0; a2 / a0 return b0, b1, b2, a1, a2然后通过“Parameter Update”模块把系数写入对应的Biquad。SigmaStudio里可以用“Write Parameter”模块通过I2C或者SPI从外部MCU写入。如果不用外部MCU可以用SigmaStudio的“Real-Time”控制界面手动更新但那样就不是自动防啸叫了。我的方案是用一个低成本的STM32或者Arduino作为主控通过I2C读写ADAU1701的参数RAM。ADAU1701的I2C地址是0x687位参数RAM的起始地址是0x0000。具体协议是先写地址高字节再写低字节然后写数据。STM32的代码大概是这样void write_param(uint16_t addr, uint32_t data) { uint8_t buf[6]; buf[0] (addr 8) 0xFF; buf[1] addr 0xFF; buf[2] (data 24) 0xFF; buf[3] (data 16) 0xFF; buf[4] (data 8) 0xFF; buf[5] data 0xFF; HAL_I2C_Master_Transmit(hi2c1, 0x681, buf, 6, 100); }实操心得ADAU1701的I2C写入需要一定延时每次写入后至少等1 ms否则会丢数据。我一开始没加延时结果参数写不进去调试了半天。4.3 防啸叫的阈值和释放时间调试阈值设得太低会误触发把正常音乐当成啸叫设得太高则反应慢。我的经验是在安静环境下把话筒增益开到刚好不啸叫然后记录此时检测频点的能量再往上加6 dB作为阈值。释放时间Release Time设500 ms到1 s这样啸叫消失后Notch不会马上撤掉避免反复触发。调试的时候我用SigmaStudio的“Real-Time”窗口观察每个Goertzel的输出找到啸叫频点。然后手动调整阈值直到既能检测到啸叫又不误触发。实测下来在KTV场景里阈值-15 dBFS释放时间800 ms效果比较稳。5. 常见问题与排查技巧实录5.1 编译报错“Program RAM overflow”怎么办这是最常见的问题。60段EQ加上防啸叫程序RAM很容易超过1024条指令。解决方法有几种一是减少EQ段数比如从60段降到40段二是用更省指令的Biquad结构比如直接型II比直接型I省三是降低采样率到44.1 kHzMIPS需求会降一点四是把防啸叫的Goertzel检测改成更简单的峰值检测比如用“Absolute”加“Comparator”。我试过把采样率降到44.1 kHz程序RAM省了大约10%但音质差别不大。另外SigmaStudio的编译器有优化选项在“Project Options”里把“Optimization”设成“Speed”或者“Size”我选“Size”能省一些程序空间。5.2 参数RAM不够用怎么优化参数RAM的优化更棘手。我的做法是固定EQ的系数不占用参数RAM因为它们在编译时已经写死在程序里了。只有动态Notch的系数需要参数RAM。40个Notch如果同时激活需要40×5200个参数字加上其他模块总共可能超过500字。但实际同时激活的Notch不会超过10个所以可以用“参数复用”把40个Notch的参数映射到10个物理Biquad上通过“Indexed Parameter”切换。这样参数RAM只需要50字左右。SigmaStudio里可以用“Multiplexer”加“Index”来实现。具体是把10个Biquad的参数输入接到一个多路复用器多路复用器的选择信号来自一个计数器计数器根据检测到的啸叫频点索引来切换。这样虽然增加了指令开销但参数RAM省了很多。5.3 啸叫检测不灵敏或误触发检测不灵敏通常是阈值设得太高或者Goertzel的窗口太短。我试过窗口10 ms发现对于快速变化的啸叫反应不够后来改成20 ms灵敏度好一些。误触发则可能是阈值太低或者检测频点太多把音乐里的强谐波当成啸叫。解决方法是只检测那些音乐里能量较少的频段比如4 kHz以上因为音乐能量主要集中在低频。另外可以加一个“持续时间”条件只有连续3个窗口都超阈值才触发。还有一个坑ADAU1701的ADC输入有直流偏置如果不加高通滤波直流会干扰Goertzel。我在输入加了一个20 Hz的高通问题解决。5.4 Notch插入后音质变差Notch的Q值太高会导致相位失真听感上会觉得声音“发紧”。我的经验是Q值不要超过20深度不要超过-18 dB。如果啸叫很严重可以插入两个Notch一个Q10一个Q30这样既能抑制啸叫又减少音质影响。另外Notch的频率要准确如果偏了啸叫抑制效果会打折扣。我用SigmaStudio的“Frequency Response”工具检查Notch的频响确保中心频率对准。注意ADAU1701的Biquad系数是定点频率计算有量化误差。对于高频Notch比如8 kHz量化误差可能导致中心频率偏移几十Hz。解决方法是用更高的Q值来补偿或者用双精度Biquad但那样更耗资源。6. 实测效果与进一步扩展的思路我在一个KTV前级上实测了这套方案ADAU1701评估板48 kHz采样率60段EQ20段固定房间校正40段动态Notch防啸叫检测20个频点。结果是在话筒增益开到最大时原本会啸叫的频点被自动抑制啸叫消失音乐播放时没有明显误触发。程序RAM用了约950条指令参数RAM用了约600字都在芯片能力范围内。如果还想进一步扩展可以考虑一是用两颗ADAU1701级联一颗专门做EQ一颗做防啸叫通过I2S连接二是把防啸叫算法升级成自适应Notch用LMS但只针对少数频点三是加入混响和延时效果但那样资源更紧张。我个人觉得ADAU1701的潜力已经被榨得差不多了再往上就得换ADAU1452或者ADAU1466那些芯片的MIPS和参数RAM都大得多。最后分享一个小技巧SigmaStudio的工程文件可以导出成XML用Git做版本管理。每次调参后提交一次这样如果调乱了可以回滚。我试过用这个方法省了不少重调的时间。另外ADAU1701的I2C地址可以通过电阻配置如果你用多颗芯片记得改地址。