FIR数字滤波器的设计及matlab实现.ppt

（3）b=fir1(n,wn,window); 它与(1)的区别在于可以自由地选择不同的窗来完成线性相位FIR滤波器的设计。 w = boxcar(M) 数组w中返回M点矩形窗函数w = triang(M) 数组w中返回M点Bartlett(三角)窗函数w = hanning(M) 数组w中返回M点汉宁窗函数w = hamming(M) 数组w中返回M点哈明窗函数w = blackman(M) 数组w中返回M点布莱克曼窗函数w = kaiser(M,beta) 数组w中返回beta值M点凯泽窗函数 注意：窗的长度应为M=n+1，它应该与h(n)的点数相同。 （4）b=fir1(n,wn,’ftype’,window); 它可以用任意的窗函数设计不同的线性相位的高通或带阻FIR滤波器 。 例4：根据下列技术指标，设计一个数字FIR低通滤波器： wp =0.2π，ws=0.3π,Rp=0.25dB，As=50dB。 （1）若选择哈明窗设计：  wp=0.2*pi;ws=0.3*pi;  tr_width=ws-wp;  M=ceil(6.6*pi/tr_width) n=[0:1:M-1];  wc=(ws+wp)/2;  h = fir1(M,wc/pi); [H,W] = freqz(h,1); plot(W/pi, 20*log10(abs(H))); （2）若选择布莱克曼窗设计：  wp=0.2*pi;ws=0.3*pi;  tr_width=ws-wp;  M=ceil(6.6*pi/tr_width) n=[0:1:M-1];  wc=(ws+wp)/2;  h = fir1(M,wc/pi, blackman(M+1)); [H,W] = freqz(h,1); plot(W/pi, 20*log10(abs(H))); 2019 ppt资料 * 欢迎批评指导！！ 快乐工作，快乐生活！ 感谢您的聆听！ ！ 加矩形窗处理后，对理想频率响应产生了两点影响： 1）使理想频率特性不连续点ω=ωc 处，形成了 一 个过渡带, 过渡带的宽度等于矩形窗的频率 响应WR(ω)的主瓣宽度△ω=4π/N ； 2) 在截止频率ωc的两边ω=ωc ±2π/N 处(即过 渡带的两边)，H(ω)出现最大的肩峰值，肩峰 的两侧形成起伏振荡，其振荡幅度取决于旁瓣 的相对幅度，而振荡的快慢，则取决于WR(ω) 波动的快慢。 若增加截取长度N，则在主瓣附近的窗的频率响应为 该函数的性质：随着x加大（即N加大），函数曲线波动的 频率加快，主瓣幅度加高，旁瓣幅度也同 样加高，主瓣与旁瓣的相对比例保持不变。 这个相对比例是由sinx/x决定的, 也就是说是由矩形窗 函数的形状决定的。 因而，当长度N增加时，只会减小过渡带宽(4π/N )，而不会改变肩峰的相对值。 在矩形窗情况下，最大相对肩峰值为8.95%，N增加时，4π/N减小，起伏振荡变密，但最大肩峰则总是8.95%， 这就是吉布斯(Gibbs)效应。 由于窗谱肩峰的存在，影响到H(ω)通带的平坦和阻带的衰减，使阻带最小衰减只有21dB左右，因此在实际中，矩形窗很少采用。 为了消除吉布斯效应，取得较好频率特性，一般采用其他类型的窗函数 ，对 进行加窗处理。 2、常用的窗函数1. 三角形窗(Bartlett Window) (6-2-11a) (6-2-11b) 其频率响应为 主瓣宽度为( ). 2. 汉宁(Hanning)窗，又称升余弦窗 (6-2-12a) (6-2-12b) 其频率响应 和幅度响应 分别为 是三项矩形窗的幅度响应 的移位加权和，使旁瓣相互抵消，能量更集中在主瓣，但主瓣宽度比矩形窗的主瓣加宽了一倍，为( ). 3. 汉明(Hamming)窗，又称改进 的升余弦窗