数字信号处理课程设计-窗函数法线性相位FIR数字滤波器设计.doc

汕 头 大 学 工 学 院 三级项目报告 课 程 名 称 ： 数字信号处理 三级项目题目：窗函数法线性相位FIR数字滤波器设计 指 导 教 师 ： **** 系 别： 电子工程系 专 业： 姓名 学号 完成时间： 2010年 12月12 日 成绩： 评阅人：**** 内容与要求 1、线性相位FIR数字滤波器应满足的条件。 2、频率选择性线性相位FIR数字滤波器的设计算法。 报告正文 一、线性相移FIR 数字滤波器的条件： 第一类线性相位： 由偶对称h(n)=h(N-1-n),可得 n=(N–1)/2为h(n)的偶对称中心 幅度函数： 1）h(n)N为奇数 2）h(n)N为偶数 第二类线性相位 由奇对称h(n)=-h(N-1-n),可得 n =(N–1)/2为h(n)的奇对称中心 幅度函数： 3）h(n)N为奇数 4）h(n)N为偶数 总结四种线性相位FIR的特点： 当h(n)为实数且偶对称时，FIR为恒相时延，相位曲线是一条过原点、以-(N-1)/2为斜率的直线。信号通过这类滤波器后，各种频率分量的时延都是(N-1)/2。当N为奇数时，时延(N-1)/2是整数，是采样间隔的整数倍，采样点时延后仍是采样点。但当N为偶数时，时延(N-1)/2不是整数，采样点时延后不在采样点位置上。同时，N为偶数时，π点是幅度的零点，不能做高通、带阻滤波器。一般情况下，第一类FIR特别适合做各种滤波器。 当h(n)为实数且奇对称时，FIR仅是恒群时延。相位曲线是一条截距为π／2，以-(N-1)/2为斜率的直线。信号通过该滤波器产生的时延也是(N-1)/2个采样周期，但另外对所有频率分量均有一个附加的90度的相移。 程序代码： 图一： h=[-4 1 -1 -2 5 6 5 -2 -1 1 -4] M=length(h); n=0:M-1; [Hr,w,a,L]=Hr_type1(h); subplot(2,2,1); stem(n,h); xlabel(n);ylabel(h(n));title(脉冲响应) subplot(2,2,3); stem(0:L,a); xlabel(n);ylabel(a(n));title(a(n)系数) subplot(2,2,2); plot(w/pi,Hr); xlabel(频率 单位pi);ylabel(Hr);title(1型幅度响应) subplot(2,2,4); pzplotz(h,1);%画极零图 图二： h=[-4 1 -1 -2 5 6 6 5 -2 -1 1 -4] M=length(h); n=0:M-1; [Hr,w,b,L]=Hr_type2(h); subplot(2,2,1); stem(n,h); xlabel(n);ylabel(h(n));title(脉冲响应) subplot(2,2,3); stem(1:L,b); xlabel(n);ylabel(b(n));title(b(n)系数) subplot(2,2,2); plot(w/pi,Hr); xlabel(频率 单位pi);ylabel(H);title(2型幅度响应) axis([0 2 -20 40]); subplot(2,2,4); pzplotz(h,1);%画极零图 图三： h=[-4 1 -1 -2 5 0 -5 2 1 -1 4] M=length(h); n=0:M-1; [Hr,w,c,L]=Hr_type3(h); subplot(2,2,1); stem(n,h); xlabel(n);ylabel(h(n));title(脉冲响应) subplot(2,2,3); stem(0:L,c); xlabel(n);ylabel(c(n));title(c(n)系数) subplot(2,2,2); plot(w/pi,Hr); xlabel(频率 单位pi);ylabel(H);title(3型幅度响应) subplot(2,2,4); pzplotz(h,1);%画极零图 图四： h=[-4 1 -1 -2 5 6 -6 -5 2 1 -1 4] M=length(h); n=0:M-1; [Hr,w,d,L]=Hr_type4(h); subplot(2,2,1); stem(n,h); xlabel(n);ylabel(h(n));title(脉冲响应) subplot(2,2,3); stem(1:L,d