matlab sprintf_如何在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))

503c97644cabe00fdb187acb6a1de710.png

三.如何使用带通滤波器对一组数据进行滤波

y = filter(b,a,x) 

利用分子和分母系数 b 和 a 定义的有理传递函数对输入数据 x 进行滤波。

设计一个matlab带通滤波器代码_百度知道​zhidao.baidu.com
74bd057aed5e4fe288c5f906dae9f7e5.png

四.如何合理确定滤波器的参数?

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

Matlab fir滤波(conv)​blog.csdn.net
8315466819a1555dfef0fe6dbc836dc4.png
matlab中fdatool使用说明 - 图文 - 百度文库​wenku.baidu.com
  • 滤波的阶数

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

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

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值