数字滤波器的设计(下).pptVIP

求取截断得h(n) 设h(n)满足要求。滤波器输出为： 本例没有给出T，T包含在ωc=0.5π＝ΩcT 中。 如果要求截止频率为1000Hz，则Ωc=1000×2π T＝ 0.5π/(1000×2π)=0.00025(秒) 采样频率=4000Hz如果要求截止频率为2000Hz，则Ωc=2000×2π T＝ 0.5π/(2000×2π)=0.000125(秒) 采样频率=8000Hz 矩形窗计算函数 汉宁窗计算函数 海明窗计算函数 布莱克曼窗计算函数 基于窗函数法FIR滤波器的matlab设计 以下函数的频率单位均为归一化频率、加hamming窗。 fir1(n,wc)——设计截止频率为wc的低通、带通滤波器，滤波器长度长度为n+1。设计带通滤波器时wc=[wc1,wc2]。 fir1(n,wc,’high’)——设计截止频率为wc的高通滤波器。如果n为偶数，滤波器长度为n+1；如果n为奇数，滤波器长度为n+2。 fir1(n,wc,’stop’)——设计截止频率为wc=[wc1,wc2]的带阻滤波器。如果n为偶数，滤波器长度为n+1；如果n为奇数，滤波器长度为n+2。 函数fir1(…，window)——可指定窗函数。 close all;clear all; n=33;wc=[0.4,0.6]; b = fir1(n,wc,‘stop’); %带阻滤波器。boxcar(n+2) [h,w]=freqz(b); figure;plot(w/pi,abs(h));grid on; n=33;wc=[0.4,0.6]; b = fir1(n,wc); %带通滤波器 [h,w]=freqz(b); figure;plot(w/pi,abs(h));grid on; n=33;wc=0.4; b = fir1(n,wc); %低通滤波器 [h,w]=freqz(b); figure;plot(w/pi,abs(h));grid on; n=33;wc=0.4; b = fir1(n,wc,high); %高通滤波器 [h,w]=freqz(b); figure;plot(w/pi,abs(h));grid on; 因为序列的离散傅立叶变换是该序列的离散频率响应值，因此将给定的理想数字滤波器的频率响应进行采样，对采样序列进行离散傅立叶反变换，即可得到所需的FIR滤波器冲击响应序列。 3.8.3 用频率抽样设计法设计FIR滤波器 序列h(n)的傅立叶变换既是数字滤波器频谱Hd(ejω)的采样值。 例：设计一线性相位FIR低通滤波器，给定h(n)长度N=15，幅频响应的采样值为： 若希望有一个序列的傅立叶变换的频谱等于Hd(ejω)，则： 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 π k H(k) 2π 为了保证线性相位，当N为奇数时，应有 h(n)=h(N－1－n) 即： h(0)= h(14)，h(1)= h(13)…h(6)= h(8)，h(7)= h(7) 可以看出序列h(n)的相频特性是线性的。 剩下的问题就是如何确定h(n)。根据幅频特性的采样值已知的条件，列写方程： 8个方程解8个未知数： 基于频率采样法fir滤波器的matlab设计 fir2(n,w,m) fir2(n,w,m,window) n——滤波器阶数。设计低通、带通滤波器时，滤波器长度为n+1。设计高通、带阻滤波器时，如果n为偶数，滤波器长度为n+1；如果n为奇数，滤波器长度为n+2。 w——频率采样点的归一化频率值矢量，[0，… ,1] m——对应频率采样点的幅频特性值矢量 window——窗函数，默认为hamming窗。 close all;clear all; n=16; w=[0, 0.5,0.5,1];f=[1,1,0,0]; %w=[0, 0.5,0.5,1];f=[0,0,1,1]; %w=[0,0.3,0.3,0.5,0.5,1];f=[1,1,0,0,1,1]; %w=[0,0.3,0.3,0.5,0.5,1];f=[0,0,1,1,0,0]; b=fir2(n,w,f,boxcar(n+1)); figure;stem([0:length(b)-1],b); [h,w]=freqz(b,1,n); figure;plot(w/pi,abs(h)); hold on;stem(w/pi,abs(h),r);grid on; 例：用频率采样法生成各类滤波器的演示 1、FIR滤波器的直接实现 由系统函数的一般形式 得到: 则: FIR滤波器的实现方法 2、FIR滤波器的快速卷积实现 FFT FFT × IFFT x(n) h(n) y(n) X(k) H(k) Y(k)= X(k)H(k) 直接计