




下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、实验四 用窗函数法设计FIR滤波器(一)实验目的1. 掌握窗函数法设计FIR滤波器的原理和方法,观察用几种常用窗函数设计的FIR数字滤波器技术指标;2. 掌握FIR滤波器的线性相位特性;3. 了解各种窗函数对滤波特性的影响。(二)实验原理如果所希望的滤波器的理想频率响应函数为Hd(ej),则其对应的单位脉冲响应为,用窗函数wN(n)将hd(n)截断,并进行加权处理,得到实际滤波器的单位脉冲响应h(n)=hd(n)wN(n),其频率响应函数为。如果要求线性相位特性,则h(n)还必须满足。可根据具体情况选择h(n)的长度及对称性。(三)实验内容1、生成四种窗函数:矩形窗、三角窗、汉宁窗、海明窗,并
2、观察其频率响应。实验代码以及运行结果%矩形窗及其频响 n=15;window1=rectwin(n);h1,w1=freqz(window1,1);subplot(2,1,1);stem(window1);title('矩形窗');subplot(2,1,2);plot(w1/pi,20*log(abs(h1)/abs(h1(1);title('矩形窗频响');%三角窗及其频响n=15;window2=triang(n);h2,w2=freqz(window2,1);subplot(2,1,1);stem(window2);title('三角窗
3、');subplot(2,1,2);plot(w2/pi,20*log(abs(h2)/abs(h2(1);title('三角窗频响');%汉宁窗及其频响n=15; window3=hann(n);window3=hann(n);h3,w3=freqz(window3,1);subplot(2,1,1);stem(window3);title('汉宁窗');subplot(2,1,2);plot(w3/pi,20*log(abs(h3)/abs(h3(1);title('汉宁窗频响');%海明窗频响n=15;window4=ha
4、mming(n);h4,w4=freqz(window4,1);subplot(2,1,1);stem(window4);title('海明窗');subplot(2,1,2);plot(w4/pi,20*log(abs(h4)/abs(h4(1);title('海明窗频响');运行结果: 2、根据下列技术指标,设计一个FIR数字低通滤波器:wp=0.2,ws=0.4,ap=0.25dB, as=50dB,选择一个适当的窗函数,确定单位冲激响应,绘出所设计的滤波器的幅度响应。提示:根据窗函数最小阻带衰减的特性表,可采用海明窗可提供大于50dB的衰减,其过渡带为6
5、.6/N,因此具有较小的阶次。实验代码%用窗函数法设计FIR滤波器 clear all;N=33;Wc=0.3*pi;%理想低通滤波器设计(怎样设计理想高通、带通、带阻滤波器?)alpha=(N-1)/2;n=0:(N-1); m=n-alpha+eps;hd=sin(Wc*m)./(pi*m);B=hamming(N);%海明窗string='Hamming ','N=',num2str(N); h=hd.*(B)' %加窗截取%以上过程可直接调用FIR1实现%wc=Wc/pi;%频率归一化%h=fir1(N-1,wc,hamming(N);H,m=f
6、reqz(h,1,1024,'whole'); %频率响应mag=abs(H);db=20*log10(mag+eps)/max(mag);pha=angle(H);subplot(2,2,1)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,strin
7、g)subplot(2,2,2)plot(m/pi,db)axis(0 1 -100 0)xlabel('w/pi')ylabel('dB')title('副频衰减特性')grid onsubplot(2,2,3)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('相频特性')subplot(2,2,4)plot(m,mag)axis(0 3.15 0 1.5)xlabel('频率W(rad)')ylabel('幅值')title('幅频特性')text(0.9,1.2,string)运行结果(四)实验心得本次实验是使用窗函数来设计FIR滤波器,在MATLAB中只需要将给定的参数输入到函数中即可马上得到处结果。通过本次实验,我对使用MATLAB快速设计滤波器的流程更为熟练,同时,也由衷地感叹MATLAB的功能强大之处,它让我们在设计时能节
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 科技园区场地租赁分成及人才引进合同
- 草场租赁与草原生态补偿及资源保护合同
- 清算还款协议书范本
- 建筑工程测量员专业服务协议
- 出租车乘客安全保障合同协议书
- 花艺沙龙培训
- 2024年“工会杯”职工技能竞赛化学检验员赛项理论考试题库(浓缩500题)
- 高三化学一轮复习 训练题-物质结构与性质
- 工厂成本方面培训
- 植物护理幼儿园
- 执业药师 中药一笔记
- 新科hg5300功放说明书
- 2023-2024学年湖南省常德市小学语文六年级期末评估试卷附参考答案和详细解析
- 华北理工环境监测教案02水和废水监测-4金属化合物-铝、汞、镉、铅、铜、锌、铬、砷等测定
- 气污染源自动监控设施台账记录模版校准记录
- JJF 1169-2007汽车制动操纵力计校准规范
- 新高考高中物理竞赛专题1力学50题竞赛真题强化训练原卷版
- 2023年成都市成华发展集团有限责任公司招聘笔试模拟试题及答案解析
- 晒纹资料大全
- 山东省泰安市2021-2022学年高一下学期期末考试数学试题(Word版含答案)
- ××团支部换届选举选票
评论
0/150
提交评论