数值信号处理教程实验四.docVIP

山东轻工业学院实验报告 课程名称：数字信号处理教程 指导教师：贺爱玲 实验日期：2011.12. 15 院（系）： 理学院 专业班级：信息09-2 实验地点：实验楼北楼 学生姓名： 学号： 200903021076 同组人：无 实验项目名称： 实验4：窗函数法FIR滤波器设计实验 一、实验目的： 利用矩形窗、汉宁窗(Hanning)、海明窗(Hamming)设计线性相位FIR低通滤波器，要求通带截止频率为：w=pi/4 求出分别对应的单位脉冲响应，并进行比较。 画出单位脉冲响应图形和对数幅度响应图形。 二、实验器材： 安装有matlab软件的计算机一台。 三、实验内容 1、画出各种窗函数图像及其单位脉冲响应，图形如下： 由上图可以很形象的比较各窗函数。 其源代码如下： window1=rectwin(n); [h1,w1]=freqz(window1,1); subplot(4,2,1); stem(window1);title(矩形窗);subplot(4,2,2); plot(w1/pi,20*log(abs(h1))/abs(h1(1)));title(矩形窗单位脉冲响应); %三角窗及其频响 window2=triang(n); [h2,w2]=freqz(window2,1); subplot(4,2,3);stem(window2);title(三角窗); subplot(4,2,4); plot(w2/pi,20*log(abs(h2))/abs(h2(1)));title(三角窗单位脉冲响应); %汉宁窗及其频响 window3=hann(n); [h3,w3]=freqz(window3,1); subplot(4,2,5);stem(window3);title(汉宁窗); subplot(4,2,6); plot(w3/pi,20*log(abs(h3))/abs(h3(1)));title(汉宁窗单位脉冲响应); %海明窗频响 window4=hamming(n); [h4,w4]=freqz(window4,1); subplot(4,2,7);stem(window4);title(海明窗); subplot(4,2,8); plot(w4/pi,20*log(abs(h4))/abs(h4(1)));title(海明窗单位脉冲响应); 其单位脉冲响应与对数幅度响应图形如下所示： 实验代码如下所示： function hd=ideal_lp1(wc,N); alpha=(N-1)/2; n=0:1:N-1; m=n-alpha+eps; hd=sin(wc*m)./(pi*m); %滤波器幅值、相位响应函数 function[db,mag,pha,w]=freqz_m2(b,a); %db:相对幅值响应 %mag:绝对幅值响应 %pha:相位响应 %w:采样频率 %b:系统函数H(Z)的分子项(对FIR,b=h) %a:系统函数H(Z)的分母项(对FIR,a=1) [H,w]=freqz(b,a,1000,whole); H=(H(1:1:501)); w=(w(1:1:501)); mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); %主程序 clear all; wp=0.2*pi; ws=0.4*pi; tr_wdith=ws-wp; N=ceil(6.6*pi/tr_wdith)+1; n=0:1:N-1; wc=(ws+wp)/2;%理想低通滤波器的截止频率 hd=ideal_lp1(wc,N);%理想低通滤波器的单位冲激响应 w_ham=(hamming(N));%海明窗 h=hd.*w_ham;%截取得到实际的单位脉冲响应 [db,mag,pha,w]=freqz_m2(h,[1]);%计算实际滤波器的幅度响应 delta_w=2*pi/1000; ap=-(min(db(1:1:wp/delta_w+1)));%实际带通波纹 as=-round(max(db(ws/delta_w+1:1:501)));%实际阻带波纹 subplot(221);stem(n,hd);title(理想单位脉冲响应hd(n)); subplot(222);stem(n,w_ham);title(海明窗w(n)); subplot(223);stem(n,h);title(实际单位脉冲响应hd(n)); subplot(224);stem(w/pi,db);title(幅度响应(dB)); axis([0,1,-10