实验指导书 7-8章.doc

实验六 FIR滤波器设计 实验目的 学习FIR数字滤波器窗口函数设计法. 学习FIR数字滤波器频率抽样设计法 进行FIR、IIR滤波器性能比较. 实验内容 1.设计一个线性相位滤波器，其理想特性如下图所示。通带内幅度为1，阻带内幅度为0。数字截止频率为，求。 理想低通特性图 2. 设计如下图所示理想特性FIR线性相位数字滤波器，求。 理想高通特性图 3. 用频率抽样法设计一个FIR数字低通滤波器，截至频率=，取。 三.实验说明 1.设数字截止频率为=,的长度N=21,求解实验内容1的MATLAB实现程序如下： --------------------------------------------------------------------------------------------------------------- clear; %清空所有数 close all; %关闭所有窗口 N=21; %N定义为21 wc=pi/4; n=0:N-1; %定义N的范围且步长为一 r=(N-1)/2; hdn=sin(wc*(n-r))/pi./(n-r); %理想滤波器 if rem(N,2)~=0 % 如果N为奇数，处理不定式 hdn(r+1)=wc/pi; end wn1=boxcar(N); %矩形窗 hn1=hdn.*wn1; %为得到的序列加窗 wn2=hamming(N); %hamming窗 hn2=hdn.*wn2; % 为得到的序列加窗 subplot(2,2,1) %画2*2个图 第一个图为hn2 stem(n,hn1,.) %绘制hn1波形 line([0,20],[0,0]); %画线 title(矩形窗设计的h(n)); %设图的标号名称 xlabel(n);ylabel(h(n)); subplot(2,2,3) stem(n,hn2,.) %绘制单位脉冲响应h(n)波形 line([0,20],[0,0]); title(hamming窗设计的h(n)); xlabel(n);ylabel(h(n)); hn11=fft(hn1,512); %进行离散傅里叶变换 w=2*[0:511]/512; subplot(2,2,2) plot(w,20*log10(abs(hn11))) grid; %添加网格线 axis([0,2,-80,5]); %坐标范围 title(幅度特性); xlabel(w/pi);ylabel(幅度(dB)); hn22=fft(hn2,512); subplot(2,2,4) plot(w,20*log10(abs(hn22))) grid; %添加网格线 axis([0,2,-80,5]); title(幅度特性); xlabel(w/pi);ylabel(幅度(dB)); --------------------------------------------------------------------------------------------------------------- 程序运行结果如下图所示。 用矩形窗和hamming窗设计的FIR低通滤波器