窗函数设计FIR滤波器技巧.doc

(式3.2.1) 由于是矩形频率响应特性，故一定是无限长序列，且是非因果的，而FIR滤波器的必然是有限长的，所以要用有限长的来逼近无限长的，最有效的方法是截断或者说用一个有限长度的窗口函数序列来截取，即 (式3.2.2) 因而窗函数序列的形状及长度的选择就是关键。 我们以一个截止频率为的线性相位的理想矩形幅度特性的低通滤波器为例来讨论。设低通特性的群延时为，即 (式3.2.3) 这表明，在通带≤范围内，的幅度是均匀的，其值为1，相位是。利用(1)式可得 (式3.2.4) 是中心点在的偶对称无限长非因果序列，要得到有限长的，一种最简单的方法就是取矩形窗，即 但是按照线形相位滤波器的约束，必须是偶对称的，对称中心应为长度的一半(N-1)/2，因而必须=(N-1)/2，所以有 (式3.2.5) 将(式3.2.4)代入(式3.25)，可得 (式3.2.6) 此时，一定满足这一线性相位的条件。 下面求的傅里叶变换，也就是找出待求FIR滤波器的频率特性，以便能看出加窗处理后究竟对频率响应有何影响。 按照复卷积公式，在时域是相乘、频域上是周期性卷积关系，即 (式3.2.7) 因而逼近的好坏，完全取决于窗函数的频率特性。 窗函数的频率特性为 (式3.2.8) 对矩形窗，则有 (式3.2.9) 也可表示成幅度函数与相位函数 (式3.2.10) 其中 (式3.2.11) 就是频域抽样内插函数，其幅度函数在之内为一个主瓣，两侧形成许多衰减振荡的旁瓣，如果将理想频率响应也写成 (式3.2.12) 则其幅度函数为 (式3.2.13) 波器； 调用格式：b = fir1(N,Wn); b = fir1(N,Wn,‘high’); b = fir1(N,Wn, ‘stop’); 参数说明：N：阶次，滤波器长度为N＋1； Wn：通带截止频率，其值在0～1之间，1对应Fs/2； b:滤波器系数。 在上述所有格式中，若不指定窗函数的类型，fir1自动选择Hamming窗。 2）fir2 ： 用来设计具有任意幅度响应的FIR滤波器。 调用格式：b = fir2(N, F, M); 参数说明：F是频率向量，其值在0～1之间； M是和F相对应的所希望的幅频相应。 如同fir1, 缺省时自动选用Hamming窗。 3）为了观测到设计出来的滤波器的特性，用freqz得到频率响应。 其中在画频率响应的时候我们分为幅度和相位画出。 又因为我们要观测的是衰减的大小程度，以dB为单位，所以我们在画幅度的时候纵坐标应该转换成dB。 观测[H,w]=freqz(b,a,N) b和分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H则包含了离散系统频响在 0pi范围内N个频率等分点的值（其中N为正整数），w则包含了范围内N个频率等分点。调用默认的N时，其值是512。可以先调用freqz()函数计算系统的频率响应，然后利用abs()和angle()函数及plot()函数，绘制出系统的频响曲线。 4.设计内容 4.1用MATLAB设计程序如下 clear;fs=500;t=(1:250)/fs; x=10*cos(2*pi*20*t)+cos(2*pi*80*t)+10*cos(2*pi*120*t); L=length(x);N=2^(nextpow2(L));Hw=fft(x,N); figure(2);subplot(2,1,1);plot(t,x); grid on;title(滤波前信号x);xlabel(时间/s);% 原始信号 subplot(2,1,2);plot((0:N-1)*fs/L,abs(Hw));% 查看信号频谱 grid on;title(滤波前信号频谱图);xlabel(频率/Hz);ylabel(振幅|H(e^jw)|); Ap=;As=20