切比雪夫1型滤波器讲解.docx
《切比雪夫1型滤波器讲解.docx》由会员分享，可在线阅读，更多相关《切比雪夫1型滤波器讲解.docx（14页珍藏版）》请在冰豆网上搜索。
切比雪夫1型滤波器讲解
目录
1课题描述……………………………………………………………………1
2设计原理……………………………………………………………………1
2.1滤波器的分类…………………………………………………………1
2.2模拟滤波器的设计指标………………………………………………1
2.3切比雪夫1型滤波器…………………………………………………2
2.3.1切比雪夫1型滤波器的设计原理……………………………3
2.3.2切比雪夫1型滤波器的设计步骤……………………………3
3脉冲响应不变法
3.1脉冲响应不变法原理……………………………………6
4设计内容………………………………………………………………………6
4.1设计步骤6
4.2用MATLAB编程实现…………………………………………………6
4.3设计结果分析…………………………………………………………10
5总结……………………………………………………………………………10
6参考文献……………………………………………………………………10
1课题描述
数字滤波器是数字信号处理的重要工具之一，它通过数值运算处理改变输入信号所含频率成分的相对比例或者滤出某些频率成分的数字器件或程序，二数字滤波器处理精度较高，体积小，稳定，重量轻，灵活，不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。
故本书课题使用MATLAB信号处理箱和运算用切比雪夫法设计数字低通滤波器。
利用脉冲响应不变法设计切比雪夫Ⅰ数字低通滤波器，通带截止频率100hz，阻带截止频率150Hz，采样频率1000hz，通带最大衰减为0.5dB，阻带最小衰减为10dB，画出幅频、相频响应曲线，并设计信号验证滤波器设计的正确性。
2设计原理
2.1.滤波器的分类
（1）从功能上分；低、带、高、带阻。
（2）从实现方法上分:
FIR、IIR
（3）从设计方法上来分：
Chebyshev（切比雪夫）,Butterworth（巴特沃斯）
（4）从处理信号分：
经典滤波器、现代滤波器
2.2模拟滤波器的设计指标
设ha（jΩ）是一个模拟滤波器的频率响应，则基于平方幅度响应J（Ω）=Ha（jΩ）的低通滤波器技术指标为：
0≤∣Ha（jΩ）∣≤1/A^2,Ωs≤∣Ω∣
其中ε为通带波动系数，pΩ和sΩ是通带和阻带边缘频率。
A为阻带衰减系数
从图知必须满足
其中参数ε和A是数字滤波器指标
2.3切比雪夫1型滤波器
2.3.1切比雪夫滤波器介绍
在巴特沃兹滤波器中，幅度响应在通带和阻带内都是单调的。
因此，若滤波器的技术要求是用最大通带和阻带的逼近误差来给出的话，那么，在靠近通带低频端和阻带截止频率以上的部分都会超出技术指标。
一种比较有效的途径是使逼近误差均匀地分布于通带或阻带内，或同时在通带和阻带内都均匀分布，这样往往可以降低所要求的滤波器阶次。
通过选择一种具有等波纹特性而不是单调特性的逼近方法可以实现这一点。
切比雪夫型滤波器就具有这种性质：
其频率响应的幅度既可以在通带中是等波纹的，而在阻带中是单调的（称为I型切比雪夫滤波器），也可以在通带中是单调的，而在阻带中是等波纹的（称为II型切比雪夫滤波器）。
I型切比雪夫滤波器的幅度平方函数是
∙其中∣∈∣＜1。
而
是滤波器在截止频率
的放大率（注意:
常用的以幅度下降3分贝的频率点作为截止频率的定义不适用于切比雪夫滤波器!
）
∙
是
阶切比雪夫多项式：
从定义切比雪夫多项式可以直接得出由Tn（x）和Tn-1（x）求Tn+1（x）的递推公式。
将三角恒等式得Tn+1（x）=2（x）T（x）-Tn-1（x）
当01时，cos−是虚数，所以Tn（x）像双曲余弦一样单调地增加;∣Hn（w）∣对于0≤w/w0≤1呈现出在1和1/（2ε^2）之间的波动；而对于w/w0〉1单调地减小。
需要用三个参量来确定该滤波器：
ε，w0和N。
在典型的设计中，用容许的通带波纹来确定ε，而用希望的通带截止频率来确定w0。
然后选择合适的阶次n，以便阻带的技术要求得到满足。
2.3.2切比雪夫1型滤波器的设计原理
切比雪夫滤波器（又译车比雪夫滤波器）是在通带或阻带上频率响应幅度等波纹波动的滤波器。
在通带波动的为“I型切比雪夫滤波器”，在阻带波动的为“II型切比雪夫滤波器”。
切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。
切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。
切比雪夫I型为例介绍其设计原理幅度平方函数用︱H（j
）︱表示
ε为小于1的正数，表示通带内幅度波动的程度，ε越大，波动幅度也越大。
Ωp称为通带截止频率。
令λ=Ω/Ωp，称为对Ωp的归一化频率。
定义允许的通带内最大衰减αp用下式表示
定义允许的通带内最大衰减pα用下式表示αp=10lgmax︱Ha（jΩ）︱^2/min︱Ha（jΩ）︱^2︱Ω︱≤Ωp
式中max︱Ha（jΩ）︱^2=1min︱Ha（jΩ）︱^2=1/（1+ε^2）
因此αp=10lg（1+ε^2）ε^2=10^（0.1αp）-1
这样，可以根据通带内最大衰减αp，可求出参数ε。
阶数N影响过渡带的宽度，同时也影响通带内波动的疏密，因为N等于通带内的最大值和最小值的总个数。
设阻带的起点频率为Ωs，则有
2.3.3切比雪夫低通滤波器的设计步骤
（1）确定低通滤波器的技术指标：
边带频率
p，它们满足
（2）求滤波器阶数N和参数
先求出1/k1，再求出阶数N，取符合条件的N的最小整数。
（3）求归一化系统函数Ga（p）
（4）将去归一化，得到实际的Ha（s）
函数说明：
用MATLAB设计切比雪夫1型低通滤波器
（1）[z,p,k]=cheb1ap（N,Rs）
该格式用于计算N阶切比雪夫1型归一化模拟滤波器的零极点和增益因子，返回长度为N的列向量Z和P，分别给出N个零点和极点的位置。
Rs是阻带最小衰减。
（2）[N,wso]=cheb1ord（wp,ws,Rp,As）
该格式用于计算切比雪夫1型数字滤波器的阶数N和阻带截止频率wso，调用参数分别为数字滤波器的通带频率和阻带边界频率的归一化值。
（3）[N,wso]=cheb1ord（wp,ws,Rp,As）
该格式用于计算切比雪夫1型模拟滤波器的阶数N和阻带频率wso。
wp.ws是实际模拟角频率。
（4）[B,A]=cheby1（N,Rs,wso,’ftype’）
该格式用于计算N阶切比雪夫1型数字滤波器系统函数的分子和分母的多项式系数向量B和A。
调用参数N和wso分别为切比雪夫1型数字滤波器的阶数和阻带截止频率的归一化值。
（5）[B,A]=cheby1（N,Rp,wso,’ftype’,’s’）
该格式用于计算N阶切比雪夫1型模拟滤波器系统函数的分子和分母多项式系数向量B和A。
调用参数N和wso分别为切比雪夫1型模拟滤波器的阶数和阻带截止频率的归一化值。
Impinvar
功能：
用“脉冲响应不变法”将模拟滤波器离散化。
[BZ,AZ]=impinvar（B,A,Fs）
Plot（X,Y）功能：
绘制线性二维图形
gridon
功能：
对当前坐标系添加主要网格线。
3脉冲响应不变法脉冲响应不变法原理
脉冲响应不变法是实现模拟滤波器数字化的一种直观而常用的方法。
它特别适合于对滤波器的时域特性有一定要求的场合。
脉冲响应不变法的核心原理是通过对连续函数ha（t）等间隔采样得到离散序列ha（nT）。
令h（n）=ha（nT），T为采样间隔。
它是一种时域上的转换方法。
一个模拟滤波器的传递函数可以用有理分式表达式表示为：
（3-1）
通过反拉普拉斯变换我们就可以得到它的冲激相应：
（3-2）
脉冲响应不变法就是要保证脉冲响应不变，即：
（3-3）
对上式的冲激相应序列
作
变换，就可以得到数字滤波器的传递函数：
（3-4）
一般来说，
的分母多项式阶次总是大于分子多项式的阶次。
假定
的没有多重极点，则式
（1）就可分解为：
（3-5）
式中
，
均为复数，
是
的极点。
其拉氏变换为
一般来说，
的分母多项式阶次总是大于分子多项式的阶次。
数字滤波器的传递函数H（z）经过合并简化，成为一般形式的有理分式传递函数
（3-6）
在讨论采样序列z变换与模拟信号拉氏变换之间关系的有关章节中，我们已经知道
（3-7）
按照z=
的关系，每一个s平面上宽度为
/T的水平条带将重迭映射到z平面上。
因此脉冲响应不变法将s平面映射到z平面，不是一个简单的一一对应的关系。
对于高采样频率（T小）的情况，数字滤波器在频域可能有极高的增益。
为此我们采用
（3-8）
（3-9）
在脉冲响应不变法设计中，模拟频率与数字频率之间的转换关系是线性的（
）。
同时，它可以保持脉冲响应不变，
=
。
因此，这一方法往往用于低通时域数字滤波器设计及相应的模拟系统数字仿真设计。
4实验内容
4.1设计步骤
（1）按一定规则将数字滤波器的技术指标转换为模拟低通滤波器的技术指标；
（2）根据转换后的技术指标使用滤波器阶数函数，确定滤波器的最小阶数N和截止频率Wc；
（3）利用最小阶数N产生模拟切比雪夫低通滤波原型。
（4）利用冲激响应不变法或双线性不变法把模拟切比雪夫滤波器转换成数字切比雪夫低通滤波器。
4.2试验程序
低通滤波器程序
fp=100;fs=150;
Rp=0.5;Rs=10;
Fs=1000;Wp=2*pi*fp;Ws=2*pi*fs;
[N,wc]=cheb1ord（Wp,Ws,Rp,Rs,'s'）;%计算阶数N和截止频率wc
N
wc
[z,p,k]=cheb1ap（N,Rp）;%计算系统函数零极点和增益
[Bap,Aap]=zp2tf（z,p,k）;%将系统函数的零极点转化为系统函数一般形式的系数
[b,a]=lp2lp（Bap,Aap,wc）;%将模拟低通滤波器原型，转换为模拟低通滤波器
[bz,az]=impinvar（b,a,Fs）;%采用脉冲响应不变法，将模拟低通滤波器，转换为数字低通滤波器
figure
（1）
[H,W]=freqz（bz,az,Fs）;%计算频率H和频率响应W
subplot（2,1,1）,plot（W/pi,20*log10（abs（H）））
xlabel（'幅频/Hz'）;ylabel（'幅度/dB'）;gridon;
subplot（2,1,2）,plot（W/pi,180/pi*unwrap（angle（H）））
xlabel（'相频/Hz'）;ylabel（'相位/^o'）;gridon;
验证滤波器的正确性程序
figure
（2）;
f1=50;
n=150;m=0:
（n-1）;t=m/Fs;%采样点数
subplot（2,3,1）
x=sin（2*pi*f1*t）;
plot（t,x）;title（'f1输入信号'）;
f2=200;
x=sin（2*pi*f2*t）;
subplot（2,3,2）
plot（t,x）;title（'f2输入信号'）;
subplot（2,3,3）
x=sin（2*pi*f1*t）+sin（2*pi*f2*t）;
plot（t,x）;title（'叠加后输入信号'）;
y=filtfilt（bz,az,x）;
subplot（2,3,4）
plot（t,y）;title（'滤波之后的信号'）
subplot（2,3,5）
plot（t,sin（2*pi*f1*t））;title（'希望输出的信号'）;
4.3实验结果分析
求得阶数N=3，通带边界频率wc=628.3125，频率在0-100hz（0-0.2
）完全通过，在100-150hz（0.2
-0.3
）通过幅度递减，150hz（0.3
-∞）后被截止。
验证信号f1=50hz,f2=200hz，观察验证后的波形，可以看出f1通过，而f2被截止，数字低通滤波器的阻带截止频率为150hz，所以设计的数字低通滤波器符合技术指标
5实验总结与体会
通过对数字信号处理的课程设计，熟悉了MATLAB的运行环境，了解更多有关于Matlab软件的知识，初步掌握了MATLAB语言在数字信号处理中一些基本库函数的调用和编写基本程序等应用；通过对数字低通滤波器的设计让我熟悉了滤波器设计的一般原理，对滤波器有了一个感性的认识；学会了数字低通滤波器设计的一般步骤；加深了对脉冲响应不变法的理解和认识。
在设计的过程中，我也认识到了自己所学知识的不足。
以前上课都是学一些基本的东西，自以为知识掌握的很熟练，实际只是了解了最表面的东西，好多理论只是不能很好的应用于实践，现在运用学到得的东西做出有实际应用价值的东西,对所学知识点进一步的理解，并进行系统化。
这也让我再次认识到知识是无尽的，只有不断的充实自己、完善自己的知识理论体系，才能够更好的胜任自己以后的工作。
设计过程中知识的不足也让我更加坚定了终身学习的决心。
6参考文献
数字信号处理（第三版）高西全丁玉美
MATLAB基础及应用（第三版）于润伟朱晓慧
MATLAB辅助现代工程数字信号处理（第二版）李益华
升级会员 