matlab中buttord用法_如何在matlab用带通滤波器进行滤波？

最新推荐文章于 2024-09-27 10:58:25 发布

weixin_39600837

最新推荐文章于 2024-09-27 10:58:25 发布

阅读量7.4k

点赞数 6

文章标签： matlab中buttord用法

本文链接：https://blog.csdn.net/weixin_39600837/article/details/112207746

版权

以带通滤波器——巴特沃斯滤波器为例

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带通滤波器代码_百度知道zhidao.baidu.com

我们已知所需要的通带截止频率，阻带截止频率，可以通过一个便捷的方式由matlab直接确定滤波器的参数。输入fdatool命令,具体方法可参见以下博客。

Matlab fir滤波(conv)blog.csdn.net

matlab中fdatool使用说明 - 图文 - 百度文库wenku.baidu.com

滤波器的阶数，就是滤波几次的意思，一阶，就是对谐波过滤一次，滤波器的阶数越高，滤波效果越好，但是，成本也会成倍的增加，因此，选择合适的阶数是非常重要的。

（如果对你帮助就点个赞呀，关注我会有更多干货分享哦~）

关注