cheb1ord 、 chey1 and freqz
切比雪夫Ⅰ型滤波器在通带或阻带上频率响应幅度等波纹波动的滤波器。滤波器通常工作在时域上。
切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小,但是在通频带内存在幅度波动。
语法 cheb1ord
①
[n,Wp] = cheb1ord(Wp,Ws,Rp,Rs)
解释
返回Chebyshev类型I滤波器的最低阶n,该滤波器在通带中损失不超过Rp dB,并且在阻带中具有至少Rs dB的衰减。还返回相应截止频率Wp的标量(或矢量)。
②
[n,Wp] = cheb1ord(Wp,Ws,Rp,Rs,‘s’)
为模拟切比雪夫I型滤波器找到最小阶数n和截止频率Wp。
切比雪夫Ⅰ型滤波器设计
[b,a] = cheby1(n,Rp,Wp)
[b,a] = cheby1(n,Rp,Wp,ftype)
解释:n为阶数,Rp为通带中损失上限,Wp为截止频率—一般通过cheb1ord获得参数
返回的到滤波器系数矩阵[b,a]
[z,p,k]=cheby1(n,Rp,E(f))设计低通、高通、带通或带阻数字Chebyshev I型滤波器,并返回其零点、极点和增益。E(f)表示为输入信号的函数。
数字滤波器的频率响应—自动控制理论的频域响应(没有经过线性化处理)
[h,w] = freqz(b,a,n)
[h,w] = freqz(sos,n)
[h,w] = freqz(d,n)
横坐标:w/pi,频率
纵坐标:20*log(10)*abs(h)
类似于脑电的滤波,选择高通与低通串联或者带通。
例子对于以1000赫兹采样的数据,设计一个低通滤波器,其通带纹波小于3分贝(定义为0至40赫兹),阻带纹波至少60分贝(定义为150赫兹至hmdwdm频率)。
Matlab代码如下