从滤波器的幅频特性得到滤波器设计
简 介本文探讨了基于ADS1244 ADC内部低通滤波器的AI辅助设计方法。通过分析芯片数据手册中的幅频特性含50Hz/60Hz陷波利用AI推导出级联梳状滤波器的传递函数H ( z ) 1 42 ⋅ ( 1 − z − 6 ) ( 1 − z − 7 ) ( 1 − z − 1 ) 2 H(z)\frac{1}{42}\cdot\frac{(1-z^{-6})(1-z^{-7})}{(1-z^{-1})^2}H(z)421⋅(1−z−1)2(1−z−6)(1−z−7)。仿真验证显示该设计能准确实现目标陷波频率50Hz/60Hz及其谐波但存在高频衰减不足实测-25dB vs 需求-60dB的问题表明需进一步优化AI设计参数以提升性能。研究为数字滤波器设计提供了AI关键词AI数字滤波器ADS124401【AI滤波器设计】一、设计背景ADS1244是一个低功耗24比特ADC转换器 它内部的低通滤波器 具有对50赫兹和60赫兹抑制的陷波频率点。 这是在它的工作频率为2.4576兆赫兹下 可以很好地抑制周围电源工频干扰信号。 那么问题来了对于这样的一个给出的幅频特性 我们是否可以推断出它内部数字滤波器的传递函数呢 关于这个问题我们直接把它的幅频特性切给AI 请他根据上述滤波器的伏频特性 计算出它对应滤波器的传递函数 豆包花了相当多的时间 最终给出了设计结果 下面我们根据豆包给出的系统函数来验证一下对应的数字滤波器的幅频特性 是否与ADS1244的幅频特性 滤波器的幅频特性相同二、设计 验证这是AR给出了滤波器的设计传递函数 通过理论分析它所对应的幅频特性 线波位置符合adS1244频率响应 下面我们假设它频率为360赫兹。 绘制这个数字滤波器对应的从0~PI之间规划频率 对应的滤波器的频率特性。绘制出数字滤波器的幅频特性 它整体上呈现低通滤波特性 并且在50赫兹 60赫兹 等处出现了陷波频率点。 将增益坐标改为对数坐标 更能体现整个滤波器的特点。 对比ADS1244数据手册给出的滤波器的幅频特性。 可以看到设计的滤波器在50赫兹和60赫兹的陷波频率是准确的 在100赫兹和120赫兹的谐波频率也是正确的 但对于更高频率处的陷波点这个滤波器就没有体现了 但是滤波器的最大的区别在于在50赫兹处 原来低通滤波器的衰减达到-60dB 但这里设计的滤波器的衰减只有25dB 远远小于AD采样器低通滤波器的衰减。▲ 图2.1 滤波器的幅频特性▲ 图2.2 滤波器的幅频特性※总结 ※本文根据aDS1244数据手册中给出内部低通滤波器的幅频特性 让AI根据对应的频率曲线 设计对应的数字滤波器 根据设计的结果进行仿真 可以看到它能够把主要的线波频率点能够很好的满足 但是对低通滤波器的衰减特性 并没有能够很好的满足 因此根据这样的结果 还需要我们通过提供更多的设计参数和提示词参数 让AI设计出合理的数字滤波器的传递函数。一、AI的结果一、幅频特性分析结合图像信息横轴频率0~180Hz、纵轴增益dB、时钟频率f CLK 2.4576 MHz f_{\text{CLK}}2.4576\,\text{MHz}fCLK2.4576MHz可从以下维度解析幅频特性1. 通带与直流特性直流增益0 Hz 0\,\text{Hz}0Hz处增益为0 dB 0\,\text{dB}0dB即增益1说明滤波器对直流信号无衰减适合保留低频基线信号如生物电信号的直流分量。低通趋势随着频率升高整体增益呈单调下降趋势符合低通滤波器特性高频信号被逐步衰减。2. 陷波阻带特性曲线在特定频率处出现深谷增益急剧降至− 100 dB -100\,\text{dB}−100dB以下形成窄带陷波对应频率为约50 Hz 50\,\text{Hz}50Hz、100 Hz 100\,\text{Hz}100Hz、150 Hz 150\,\text{Hz}150Hz50Hz工频及其谐波60 Hz 60\,\text{Hz}60Hz、120 Hz 120\,\text{Hz}120Hz、180 Hz 180\,\text{Hz}180Hz60Hz工频及其谐波。这些陷波用于抑制全球不同地区的工频干扰50Hz/60Hz及其谐波典型应用场景为生物医学信号如心电、脑电处理。3. 滚降与阻带衰减滚降速度从通带0Hz到第一个陷波50/60Hz的增益下降极快说明滤波器阶数较高或采用了级联结构增强滚降特性。阻带整体衰减除陷波频率外其他高频如80Hz、110Hz、130Hz的增益也显著降低整体阻带抑制效果良好。4. 数字滤波器特征结合时钟频率与频率轴0~180Hz对应奈奎斯特频率f s / 2 180 Hz f_s/2180\,\text{Hz}fs/2180Hz即采样频率f s 360 Hz f_s360\,\text{Hz}fs360Hz该滤波器为数字低通多陷波结构零点位于单位圆上陷波频率极点靠近单位圆若为IIR或无极点若为FIR梳状结构。二、传递函数推导基于梳状滤波器级联模型曲线的陷波位置与滑动平均梳状滤波器的零点特性高度吻合梳状滤波器传递函数为H comb , N ( z ) 1 N ⋅ 1 − z − N 1 − z − 1 H_{\text{comb},N}(z) \frac{1}{N} \cdot \frac{1 - z^{-N}}{1 - z^{-1}}Hcomb,N(z)N1⋅1−z−11−z−N其幅频响应零点为f k k f s N f_k \frac{k f_s}{N}fkNkfsk 1 , 2 , … , N − 1 k1,2,\dots,N-1k1,2,…,N−1可抑制特定频率及其谐波。1. 级联结构假设为同时抑制50Hz和60Hz工频及其谐波采用两个梳状滤波器级联抑制60Hz谐波取N 2 6 N_26N26零点为f k 60 k Hz f_k60k\,\text{Hz}fk60kHzk 1 , 2 , 3 k1,2,3k1,2,3即60/120/180Hz抑制50Hz谐波取N 1 7 N_17N17零点为f k 360 k 7 ≈ 51.4 k Hz f_k\frac{360k}{7}\approx51.4k\,\text{Hz}fk7360k≈51.4kHzk 1 , 2 , 3 k1,2,3k1,2,3即≈51.4/102.8/154.3Hz与曲线中50/100/150Hz陷波高度吻合。2. 传递函数形式级联后传递函数为KaTeX parse error: Expected EOF, got at position 2: ̲emsp;emsp;H(z)…代入N 1 7 N_17N17、N 2 6 N_26N26并归一化直流增益H ( 1 ) 1 H(1)1H(1)1得H ( z ) 1 42 ⋅ ( 1 − z − 6 ) ( 1 − z − 7 ) ( 1 − z − 1 ) 2 \boxed{H(z) \frac{1}{42} \cdot \frac{(1 - z^{-6})(1 - z^{-7})}{(1 - z^{-1})^2}}H(z)421⋅(1−z−1)2(1−z−6)(1−z−7)展开形式KaTeX parse error: Expected EOF, got at position 2: ̲emsp;emsp;H(z)…3. 验证直流增益z 1 z1z1时分子分母均为0 / 0 0/00/0型用洛必达法则得极限为1符合0 dB 0\,\text{dB}0dB要求零点验证( 1 − z − 6 ) 0 (1-z^{-6})0(1−z−6)0对应z 6 1 z^61z61零点为z e j 2 π k / 6 ze^{j2\pi k/6}zej2πk/6对应60/120/180Hz( 1 − z − 7 ) 0 (1-z^{-7})0(1−z−7)0对应z 7 1 z^71z71零点为z e j 2 π k / 7 ze^{j2\pi k/7}zej2πk/7对应≈51.4/102.8/154.3Hz与曲线陷波位置完全一致。补充说明若实际采样频率f s ≠ 360 Hz f_s \neq 360\,\text{Hz}fs360Hz仅需按f k k f s N f_k \frac{k f_s}{N}fkNkfs调整N 1 N_1N1、N 2 N_2N2传递函数形式保持不变。该结构的核心优势是用低阶IIR结构实现多频陷波低通特性适合高时钟频率如2.4576MHz的数字系统过采样生物信号处理。要不要我帮你补充一下这个滤波器的单位冲激响应和零极点分布分析这样能更直观地理解它的工作原理。● 相关图表链接:图2.1 滤波器的幅频特性图2.2 滤波器的幅频特性