《数字信号处理》课程设计FIR窗函数滤波器的设计

1、 数字信号处理课程设计说明书课程设计名称：信号分析与处理课程设计课程设计题目：用窗函数法设计fir数字滤波器初始条件：1. matlab6.5以上版本软件；2. 课程设计辅导资料：“matlab语言基础及使用入门”、“信号与系统”、“数字信号处理原理与实现”、“matlab及在电子信息课程中的应用”等；3. 先修课程：信号与系统、数字信号处理、matlab应用实践及信号处理类课程等。要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1. 课程设计时间：1周；2. 课程设计内容：用窗函数法设计fir数字滤波器，具体包括：用窗函数法设计fir数字滤波器基本方法，各种窗

2、函数的应用、比较、选择，线性相位fir滤波器的设计等；3. 本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析，针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结，按要求进行实验演示和答辩等；4. 课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括： 目录； 与设计题目相关的理论分析、归纳和总结； 与设计内容相关的原理分析、建模、推导、可行性分析； 程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结

3、； 课程设计的心得体会（至少500字）； 参考文献（不少于5篇）； 其它必要内容等。时间安排： 1周（第18周） 周一 在图书馆和上网查资料周二 确定大概的设计方案周三 看懂别人写过的精彩例子 周四 自己用matlab编写课程设计 周五 调试程序和向别人请教周六 课程设计答辩附具体设计内容：1. 用矩形窗设计一个fir线性相位数字低通滤波器，已知。求出并画出幅频响应特性曲线。2. 用三角形窗设计一个fir线性相位数字低通滤波器，已知；要求画出滤波器的曲线。3. 用布莱克曼窗设计一个fir线性相位相移的数字带通滤波器设。试求出的表达式，并画出曲线。基于fir窗函数滤波器的设计(1)矩形窗fir滤

4、波器设计clear all;n=21;wc=pi/2;%理想低通滤波器设计(怎样设计理想高通、带通、带阻滤波器？)alpha=(n-1)/2;n=0:(n-1); m=n-alpha+eps;hd=wc/pi*sin(wc*m)./(pi*m);% hd=(sin(pi*m)-sin(wc*m)./(pi*m);b=boxcar(n);%矩形窗string=boxcar ,n=21,num2str(n); h=hd.*(b); %加窗截取%以上过程可直接调用fir1实现%wc=wc/pi;%频率归一化%h=fir1(n-1,wc,boxcar(n);h,m=freqz(h,1,1024,who

5、le); %频率响应 mag=abs(h);db=20*log10(mag+eps)/max(mag);pha=angle(h);figure(1);stem(b);axis(-2 25 0 1);grid;xlabel(n);title(三角窗函数);figure(2);n=0:n-1;stem(n,h,.)axis(0 n-1 -0.1 0.3)hold onn=0:n-1;x=zeros(n);plot(n,x,-)hold offxlabel(n)ylabel(h(n)title(实际低通滤波器的h(n)text(0.3*n,0.27,string)figure(3);plot(m/p

6、i,db)axis(0 1 -100 5)xlabel(w/pi)ylabel(db)title(幅频衰减特性)grid onfigure(4);plot(m,pha)hold onn=0:7;x=zeros(8);plot(n,x,-)hold offaxis(0 3.15 -4 4)xlabel(频率(rad)ylabel(相位(rad)title(相频特性)figure(5);plot(m,mag)axis(0 3.15 0 0.6)xlabel(频率w(rad)ylabel(幅值)title(幅频特性)text(0.9,1.2,string)(2)三角窗函数fir滤波器设计clear

7、all;n=51;wc=pi/2;%理想低通滤波器设计(怎样设计理想高通、带通、带阻滤波器？)alpha=(n-1)/2;n=0:(n-1); m=n-alpha+eps;hd=wc/pi*sin(wc*m)./(pi*m);% hd=(sin(pi*m)-sin(wc*m)./(pi*m);b=bartlett(n);%矩形窗string=boxcar ,n=,num2str(n); h=hd.*(b); %加窗截取%以上过程可直接调用fir1实现%wc=wc/pi;%频率归一化%h=fir1(n-1,wc,boxcar(n);h,m=freqz(h,1,1024,whole); %频率响应

8、 mag=abs(h);db=20*log10(mag+eps)/max(mag);pha=angle(h);figure(1);stem(b);axis(0 52 0 1);grid;xlabel(n);title(三角窗函数);figure(2);n=0:n-1;stem(n,h,.)axis(0 n-1 -0.1 0.3)hold onn=0:n-1;x=zeros(n);plot(n,x,-)hold offxlabel(n)ylabel(h(n)title(实际低通滤波器的h(n)text(0.3*n,0.27,string)figure(3);plot(m/pi,db)axis(0

9、 1 -50 1)xlabel(w/pi)ylabel(db)title(幅频衰减特性)grid onfigure(4);plot(m,pha)hold onn=0:7;x=zeros(8);plot(n,x,-)hold offaxis(0 3.15 -4 4)xlabel(频率(rad)ylabel(相位(rad)title(相频特性)figure(5);plot(m,mag)axis(0 3.15 0 0.6)xlabel(频率w(rad)ylabel(幅值)title(幅频特性)text(0.9,1.2,string)(3)布拉克曼窗函数fir滤波器设计wls = 0.2*pi; wl

10、p = 0.35*pi; whp = 0.65*pi; wc = wlp/pi,whp/pi; b = wlp-wls; n = 51; n=0:n-1; window=blackman(n); h1,w=freqz(window,1); figure(1); stem(window); xlabel(n); title(blackman窗函数); figure(2); plot(w/pi,20*log(abs(h1)/abs(h1(1); grid; xlabel(w/pi); ylabel(幅度(db); title(blackman窗函数的频谱); hn = fir1(n-1,wc,blackman(n); h2,w=freqz(hn,1,512); figure(3); stem(n,hn); xlabel(n); ylabel(h(n); title(blackman窗函数的单位脉冲响应);figure(4);