| 如何在matlab用带通滤波器进行滤波？ | 您所在的位置：网站首页 › matlab中fs › 如何在matlab用带通滤波器进行滤波？ |
如何在matlab用带通滤波器进行滤波？
|
以带通滤波器——巴特沃斯滤波器为例 matlab给出的代码如下：[n,Wn] = buttord(Wp,Ws,Rp,Rs) 是不是显得很不人性化，即便是看了help中的英文解释，仍旧是对输入输出参数一知半解~那么我就用简单的语言来描述一下~ 首先要先理解一个带通滤波器的幅频特性图中，分为通带、阻带、过渡带三大部分。 通带截止频率和阻带截止频率对这三部分进行划分，二者之间的部分即为过渡带。 一.如何设计带通滤波器?1.Wp,Ws分别是通带频率和阻带截止频率 这两个频率是指同一侧的通带和阻带频率，比如设计一个40Hz~100Hz的带通滤波器 Wp1=40Hz,Wp2=100Hz,代表Wp1~Wp2之间是通带 Ws1=20Hz,Ws2=120Hz,代表0~20Hz和120Hz~正无穷为阻带 过渡带宽度=阻带边缘频率-通带边缘频率（高频的一侧）=Ws2-Wp2=20Hz 2.Rp,Rs是描述通带、阻带幅值的两个重要参量 Rp是描述通带波纹（起伏程度）的一个参量，通带纹波是指在滤波器的频响中通带的最大幅值和最小幅值之间的差值，正常的纹波一般小于1db。 通带波纹当然越小越好，这样通带内频率的幅度都基本稳定在单倍幅度上，因此Rp是允许的通带波纹的最大值。 Rs是描述阻带衰减的一个参量 阻带衰减越大越好，衰减越大代表对不想要的信号频率成分的滤除效果越好，因此Rs是允许的需要达到的阻带衰减的最小值。 3.n,Wn是设计得到滤波器的阶数和截止频率 通道滤波器的设计代码例程如下： Wp = [40 100]/500; Ws = [20 120]/500; Rp = 3; Rs = 40; [n,Wn] = buttord(Wp,Ws,Rp,Rs)上述在Wp,Ws的赋值过程中进行了归一化，因此run之后得到的Wn的值也是归一化过的 n = 12 Wn = 0.0799 0.2002二.如何绘制带通滤波器的幅频响应?可以先调用butter()函数返回巴特沃斯滤波器传递函数的零极点系数和增益。 [z,p,k] = butter(n,Wn); sos = zp2sos(z,p,k); freqz(sos,512,1000) title(sprintf('n = %d Butterworth Lowpass Filter',n)) 三.如何使用带通滤波器对一组数据进行滤波y = filter(b,a,x) 利用分子和分母系数 b 和 a 定义的有理传递函数对输入数据 x 进行滤波。 四.如何合理确定滤波器的参数？我们已知所需要的通带截止频率，阻带截止频率，可以通过一个便捷的方式由matlab直接确定滤波器的参数。输入fdatool命令,具体方法可参见以下博客。 滤波的阶数滤波器的阶数，就是滤波几次的意思，一阶，就是对谐波过滤一次，滤波器的阶数越高，滤波效果越好，但是，成本也会成倍的增加，因此，选择合适的阶数是非常重要的。 （如果对你帮助就点个赞呀，关注我会有更多干货分享哦~） |
【本文地址】