- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验5__FIR数字滤波器设计与软件实现
实验五:FIR数字滤波器设计与软件实现
10.5.1 实验指导
1.实验目的
(1)掌握用窗函数法设计FIR数字滤波器的原理和方法。
(2)掌握用等波纹最佳逼近法设计FIR数字滤波器的原理和方法。
(3)掌握FIR滤波器的快速卷积实现原理。
(4)学会调用MATLAB函数设计与实现FIR滤波器。
2. 实验内容及步骤
(1)认真复习第七章中用窗函数法和等波纹最佳逼近法设计FIR数字滤波器的原理;
(2)调用信号产生函数xtg产生具有加性噪声的信号xt,并自动显示xt及其频谱,如图10.5.1所示;
?
图10.5.1 具有加性噪声的信号x(t)及其频谱如图
(3)请设计低通滤波器,从高频噪声中提取xt中的单频调幅信号,要求信号幅频失真小于0.1dB,将噪声频谱衰减60dB。先观察xt的频谱,确定滤波器指标参数。
(4)根据滤波器指标选择合适的窗函数,计算窗函数的长度N,调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序,调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。
(4)重复(3),滤波器指标不变,但改用等波纹最佳逼近法,调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。
提示:1MATLAB函数fir1和fftfilt的功能及其调用格式请查阅本书第7章和第?章;
2采样频率Fs=1000Hz,采样周期T=1/Fs;
3根据图10.6.1(b)和实验要求,可选择滤波器指标参数:通带截止频率fp=120Hz,阻带截至频率fs=150Hz,换算成数字频率,通带截止频率,通带最大衰为0.1dB,阻带截至频率 ,阻带最小衰为60dB。]
4实验程序框图如图10.5.2所示,供读者参考。
?
?
Fs=1000,T=1/Fs
xt=xtg
产生信号xt, 并显示xt及其频谱
用窗函数法或等波纹最佳逼近法
设计FIR滤波器hn
对信号xt滤波:yt=fftfilt(hn,xt)
1、计算并绘图显示滤波器损耗函数
2、绘图显示滤波器输出信号yt
End?
图10.5.2? 实验程序框图
?
4.思考题
(1)如果给定通带截止频率和阻带截止频率以及阻带最小衰减,如何用窗函数法设计线性相位低通滤波器?请写出设计步骤.
(2)如果要求用窗函数法设计带通滤波器,且给定通带上、下截止频率为 和 ,阻带上、下截止频率为 和,试求理想带通滤波器的截止频率 。
(3)解释为什么对同样的技术指标,用等波纹最佳逼近法设计的滤波器阶数低?
5.实验报告要求
(1)对两种设计FIR滤波器的方法(窗函数法和等波纹最佳逼近法)进行分析比较,简述其优缺点。
(2)附程序清单、打印实验内容要求绘图显示的曲线图。
(3)分析总结实验结果。
(4)简要回答思考题。
?
?
=========================================================================================
?
?
?
信号产生函数xtg程序清单
function xt=xtg(N)
%实验五信号x(t)产生,并显示信号的幅频特性曲线
%xt=xtg(N) 产生一个长度为N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz
%载波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz.
Fs=1000;T=1/Fs;Tp=N*T;
t=0:T:(N-1)*T;
fc=Fs/10;f0=fc/10; %载波频率fc=Fs/10,单频调制信号频率为f0=Fc/10;
mt=cos(2*pi*f0*t);? %产生单频正弦波调制信号mt,频率为f0
ct=cos(2*pi*fc*t);? %产生载波正弦波信号ct,频率为fc
xt=mt.*ct;????????? %相乘产生单频调制信号xt
nt=2*rand(1,N)-1;?? %产生随机噪声nt
%=======设计高通滤波器hn,用于滤除噪声nt中的低频成分,生成高通噪声=======
fp=150; fs=200;Rp=0.1;As=70;???? % 滤波器指标
fb=[fp,fs];m=[0,1]; ???????????? ??% 计算remezord函数所需参数f,m,dev
dev=[10^(-As/20),(10^(Rp/20)-1)/(10^(Rp/20)+1)];
[n,fo,mo,W]=remezord(fb,m,dev,Fs);?? % 确定remez函数所需参数
hn=remez(n,fo,mo,W);?? ? % 调用remez函数进行设
您可能关注的文档
最近下载
- 胆总管结石小讲课.pptx VIP
- 基于stm32嵌入式技术的智能药盒设计.docx VIP
- 广西公务员考试《行测》通关模拟试题及答案解析【2022】:70 - 行测模拟题.docx VIP
- 学堂课程在线“一带一路”沿线国家跨文化商务交际(大连外国语大学)期末测试答案.docx VIP
- 教职工篮球组队训练方案.docx
- 第29课对手是一面镜子ppt.ppt
- 桥头搭板施工技术交底.docx VIP
- JTG 3432—2024 公路工程集料试验规程.docx
- 05s502-砖砌矩形水表井.pdf
- 09作文-湖南省长沙市五年(2018-2022)小升初语文卷真题分题型分层汇编(共28题).doc VIP
文档评论(0)