利用Matlab工具生成滤波器

利用Matlab工具生成滤波器

简介

利用Matlab生成IIR型滤波器,IIR型会引入输出项,重新计算出新的输出项;使用matlab的滤波器设计工具,得到系数 a , b a,b a,b,再利用推导出来的公式,带入即可;

Matlab生成系数

以低通滤波器为例:

  • F c F_c Fc截止频率,也就是说从通带开始到截止频率处,dB值会下降3dB
  • F s F_s Fs采样率,采样率要尽可能大一点,过小会造成频谱泄露
  • 指定的阶数这里取2,就是参数 a a a会使用到2个

使用IIR型的巴特沃斯函数来构造的参数,其函数图像如下:

函数图如上,随着频率的上升,dB值呈单调递减的效果,也就是最后经过我们处理后的音频数据也会呈现这种效果;

从上面我们就拿到了 a , b a,b a,b系数,带入下面公式即可

低通滤波器公式

引入IIR型推导出的2阶公式,此公式由前人推导而出,推导相对复杂,后面会有文章简单讲解:
y ( n ) = b 0 ∗ x ( n ) + b 1 ∗ x ( n − 1 ) + b 2 ∗ x ( n − 2 ) − a 1 ∗ y ( n − 1 ) − a 2 ∗ y ( n − 2 ) y(n) = b_0*x(n) + b_1*x(n-1) + b_2*x(n-2) - a_1*y(n-1) - a_2*y(n-2) y(n)=b0∗x(n)+b1∗x(n−1)+b2∗x(n−2)−a1∗y(n−1)−a2∗y(n−2)

最后,将待处理的数据带入公式接口,大致算法如下:

python 复制代码
for index in range(len(data)):
    if index < 2:
        low_signal[index] = data[index]
        continue
    low_signal[index] = b_p[0] * data[index] + b_p[1] * data[index-1] + b_p[2] * data[index-2] - a_p[1] * low_signal[index-1] - a_p[2]*low_signal[index-2]

最后运行的效果如下:

其它的高通滤波器、带通和带阻滤波器也可以按照这种方法做,是不是很简单;还是要理解好原理才重要;

相关推荐
带娃的IT创业者12 分钟前
如何开发一个教育性质的多线程密码猜测演示器
网络·python·算法
Aczone282 小时前
硬件(六)arm指令
开发语言·汇编·arm开发·嵌入式硬件·算法
luckys.one6 小时前
第9篇:Freqtrade量化交易之config.json 基础入门与初始化
javascript·数据库·python·mysql·算法·json·区块链
~|Bernard|7 小时前
在 PyCharm 里怎么“点鼠标”完成指令同样的运行操作
算法·conda
战术摸鱼大师7 小时前
电机控制(四)-级联PID控制器与参数整定(MATLAB&Simulink)
算法·matlab·运动控制·电机控制
Christo38 小时前
TFS-2018《On the convergence of the sparse possibilistic c-means algorithm》
人工智能·算法·机器学习·数据挖掘
好家伙VCC9 小时前
数学建模模型 全网最全 数学建模常见算法汇总 含代码分析讲解
大数据·嵌入式硬件·算法·数学建模
liulilittle10 小时前
IP校验和算法:从网络协议到SIMD深度优化
网络·c++·网络协议·tcp/ip·算法·ip·通信
bkspiderx12 小时前
C++经典的数据结构与算法之经典算法思想:贪心算法(Greedy)
数据结构·c++·算法·贪心算法
中华小当家呐13 小时前
算法之常见八大排序
数据结构·算法·排序算法