设计题目与题目要求.doc

_一.二. 设计题目与题目要求设计一个模拟巴特沃斯滤波器，技术指标：通带频率2000-3000HZ,通带波纹1dB;两侧过渡带宽300HZ，阻带最小衰减50dB，画出其幅度和相位谱。三. 设计原理有分析可知，题目要求的是巴特沃斯带通滤波器，须先由低通变换为带通。巴特沃斯低通滤波器原理：1 设计原理a) 滤波器分类滤波器按所处理的信号分为模拟滤波器和数字滤波器。按所通过信号的频段分为低通、高通、带通和带阻滤波器四种。低通滤波器：它允许信号中的低频或直流分量通过，抑制高频分量或干扰和噪声。高通滤波器：它允许信号中的高频分量通过，抑制低频或直流分量。带通滤波器：它允许一定频段的信号通过，抑制低于或高于该频段的信号、干扰和噪声。带阻滤波器：它抑制一定频段内的信号，允许该频段以外的信号通过。b) 模拟滤波器的设计指标模拟滤波器的设计指标参数有, p, s, s。其中p和s分别称为通带边界频率和阻带边界频率，p称为通带最大衰减（即通带0, p中允许A（）的最大值）s称为通带最小衰减（即阻带s上允许A（）的最小值），p和s的单位dB。p =-10lg|Ha(jp)| s =-10lg|Ha(js| c称为3dB截止频率|Ha(jc)=|2 /2=0.707-20lg|Ha(jc)|=3dB滤波器的技术指标给定后，需要设计一个传输函数Ha（s），希望其幅度平方函数满足给定的指标p和s，一半滤波器的单位冲击响应为实数，有：Ha（j）=Ha（j）Ha*（j）=Ha（s）Ha（-s）s=jc) 巴特沃斯型滤波器i. 巴特罗夫带通滤波器的设计原理P=p（s 0）Bws（6.2.53）在p平面与s平面虚轴上的频率关系为=p（0 ）Bw（6.2.54）式中，Bw=pu-pl，表示带通滤波器的通带宽度，pl和pu分别为带通滤波器的通带下截止频率和通带上截止频率；0称为带通滤波器的中心频率。根据式（6.2.54）的映射关系，频率=0映射为频率=，频率=p映射为频率pu和-pl，频率=p映射为频率pu和-pl，也就是说，将低通滤波器G（p）的通带-p, p映射为带通滤波器的通带-pu,- pl和pu，pl。同样道理，频率=p映射为频率su和-，频率=-s映射为频率-su和sl。所以将式（6.2.53）带入式（6.2.48），就将Q(P)转换为带通滤波器的系统函数，即HBP（s）=Q(p) p=p（s +0 ）Bw 可以证明plpu=slsu=0 所以，带通滤波器的通带（阻带）边界频率关于中心频率0几何对称。如果原指标给定的边界频率不满足式（6.2.56），将要改变其中一个边界频率，以便满足式（6.2.56），丹药保证翻遍后的指标高于原始指标。具体方法是，如果plpuslsu，则减小pl（或增大su）使式（6.2.56）得到满足，具体计算公式为pl=slsupu或sl=plpusu减小pl使铜带宽度大于原指标要求的通带宽度，增大sl或减小pl都使左边的过渡带宽度小于原指标要求的过渡带宽度；反之，如果plpuslsu，则减小su（或增大pu）使式（6.2.56）得到满足，而且在关于中心频率0几何对称的两个正频率点上，带通滤波器的幅值相等。综上所述，低通到带通的边界频率及幅频响应特性的映射关系如图6.2.14所示，低通原型的每一个边界频率都映射为带通滤波器的两个相应的边界频率，途中标出了设计时有用的频率对应关系。ii. 巴特沃斯带通滤波器的设计步骤A、确定模拟带通滤波器的技术指标即：带通上限频率u，带通下限频率l下阻带上限频率s1 ,上阻带下限频率s2通带中心频率0 =lu,通带宽度B=u-l与以上边界频率对应的归一化边界频率如下:sl=slB, s2=s2B, l=lBu=uB, =luB、确定归一化低通技术要求： p=1，s=(s2 0 ) s2,s=(sl 0 ) s1s与-s的绝对值可能不相等,一般取绝对值小的s,这样保证在较大的s处更能满足要求。通带最大衰减仍为p,阻带最小衰减亦为s。C、设计归一化低通G（p）。D、将G(p)转换成带通H（s）。d)用MATLAB工具箱函数设计巴特沃斯滤波器MATLAB信号处理工具箱函数buttap，buttord和butter是巴特沃斯滤波器设计函数，其5中调用格式如下。1） Z,P,K=buttap(N)该格式用于计算N阶巴特沃斯归一化（3dB截止频率c=1）模拟抵用原型滤波器喜用函数的零点、极点和增益因子。返回长度为N的列向量Z和P，分别给出N个零点和极点的位置，表示滤波器增益。得到的系统函数为如下形式：（）（）（）（）（）（）（）式中，（）和分别为向量和的第个元素。如果要从计算得到的零极点得到喜用函数的根子和分母多项式系数向量和，可以调用函数2)N,wc=buttord(wp,ws,Rp,As,s)该格式用于计算巴特沃斯数字滤波器的阶数和截止频率。调用参数wp和ws分别为数字滤波器的通带边界频率和阻带边界频率的归一化值。要求0wp1,0ws1，表示数字频率（对应模拟频率，表示采样频率）。和分别为通带最大衰减和阻带最小衰减。当wswp时，为高通滤波器，当wp和ws为二院矢量时，为带通或带阻滤波器，这是也是二院向量。和作为函数的调用参数。3)BB,AB=butter(N,wc,ftype,s)计算巴特沃斯数字滤波器喜用函数分子和分母多项式的系数向量和。调用参数和分别为巴特沃斯模拟滤波器的阶数和截止频率（实际角频率）。由系数向量和A写出模拟滤波器的系统函数为Ha（s）=B(s) A(s)=B(1) B(2)z-1+B(N)z-（N-1）+B（N+1）z-nA(1) A(2)z-1+A(N)z-（N-1）+A（N+1）z-n由于高通滤波器和低通滤波器都只有一个3dB截止频率wc，因此仅由调用参数wc不能区别要设计的是高通还是低通滤波器。当然仅由二维向量wc也不能区分带通和带阻。所以用参数ftype来区分。Ftype=high时，设计3dB截止频率为wc的带阻滤波器，此时wc为二元向量wcl,wcu,wcl/wcu分别为带阻滤波器的通带3dB下截止频率和上截止频率，缺省ftype时设计带通滤波器，通带为频率区间wclwwcu.因当注意，设计的带通和带阻滤波器系统函数时2N阶。这是因为带通滤波器相当于N阶低通滤波器与N阶高通滤波器3.设计内容a）用MATLAB编程实现b） 设计结果分析从仿真过程和结果可以得出：IIR模拟滤波器的阶数为：n=17截止频率wc= 1.0e+004 *1.2426 1.9062系统函数H（z）的分子系数为：BB = 1.0e+064 * Columns 1 through 7 0 0 0 0 0 0 0 Columns 8 through 14 0 0 0 0 0 0 0 Columns 15 through 21 0 0 0 9.3880 0 0 0 Columns 22 through 28 0 0 0 0 0 0 0 Columns 29 through 35 0 0 0 0 0 0 0系统函数H（z）的分母系数为：AB= 1.0e+142 * Columns 1 through 7 0 0 0 0 0 0 0 Columns 8 through 14 0 0 0 0 0 0 0 Columns 15 through 21 0 0 0 0 0 0 0 Columns 22 through 28 0 0 0 0