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