0136-胡国庆-实验4FIR滤波器设计.docxVIP

数字信号处理A 实 验 报 告 实验项目名称: FIR滤波器设计 学 院:______计算机与通信工程____ 专 业:______ _通信工程 _________ 学 号:______201454080136______ 班 级:______ 通信1401________ 报 告 人:________胡国庆 __________ 指 导 老 师:___ 胡双红 _ _______ 实 验 时 间:_______2016-12-5________  实验四 FIR滤波器设计 一、实验目的： 1、掌握不同窗函数的性能指标； 2、根据指标选择窗函数设计FIR滤波器。 3、掌握频率采样法设计FIR滤波器 二、实验内容： 自定滤波器的类型（低通，高通或带通）、设计指标（通带截止频率、通带波纹、阻带截止频率、阻带衰减） 1、根据指标选择合适的窗函数，用窗口设计法设计符合指标的FIR滤波器；并验证其性能是否满足预定指标。 参考书上例题 2、根据指标选择合适的采样点数，用频率采样法设计符合指标的FIR滤波器；并验证其性能是否满足预定指标。 参考书上例题 求出通带内的Rp和阻带内的As，或者用两行两列的子图分别将结果列出来。 3、求输入x(n)=[cos(w1n)+cos(w2n)]u(n)经过系统后的输出y(n)。其中w1位于通带内，w2位于阻带内，要求做一个两行两列的子图，第一个做x的时域波形，第二个做x的幅度谱，第三个做y的时域波形，第四个做y的幅度谱。 三、实验平台： MATLAB集成系统 四、设计流程： 程序清单： 调用函数1 function [Y]=dtft (x,w) N=length(x); Y=0; for n=0:N-1 Y=Y+x(n+1).*(exp(-j*w*n)) end end 调用函数2 function[db,mag,pha,grd,w]=freqz_m(b,a); [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); grd=grpdelay(b,a,w); 主程序 wp=0.3*pi;ws=0.45*pi;deltaw=ws-wp; N0=ceil((6.6*pi)/deltaw); n=0:N0-1; wd=0.54-0.46*cos(2*pi*n/(N0-1)); wd=hamming(N0); tao=(N0-1)/2; wc=(wp+ws)/2; hd=sin(wc*(n-tao))./(pi*(n-tao)); N1=45; n=0:N1-1; tao1=(N1-1)/2; ha=sin(wc*(n-tao1))./(pi*(n-tao1)) ha = Columns 1 through 10 0.0102 -0.0058 -0.0159 -0.0064 0.0125 0.0173 -0.0000 -0.0196 -0.0161 0.0094 Columns 11 through 20 0.0265 0.0111 -0.0225 -0.0327 0.0000 0.0420 0.0375 -0.0244 -0.0796 -0.0406 Columns 21 through 30 0.1125 0.2941 NaN 0.2941 0.1125 -0.0406 -0.0796 -0.0244 0.0375 0.0420 Columns 31 through 40 0.0000 -0.0327 -0.0225 0.0111 0.0265 0.0094 -0.0161 -0.0196 -0.0000 0.0173 Columns 41 through 45 0.0125 -0.0064 -0.0159 -0.0058 0.0102 eps ans = 2.2204e-016 hd=sin(wc*(n-tao1+eps))./(pi*(n-tao1+eps)) hd = Columns 1 through 9 0.0102 -0.0058 -0.0159 -0.0064 0.