欢迎来到课桌文档! | 帮助中心 课桌文档-建筑工程资料库
课桌文档
全部分类
  • 党建之窗>
  • 感悟体会>
  • 百家争鸣>
  • 教育整顿>
  • 文笔提升>
  • 热门分类>
  • 计划总结>
  • 致辞演讲>
  • 在线阅读>
  • ImageVerifierCode 换一换
    首页 课桌文档 > 资源分类 > DOC文档下载  

    数字信号处理实验5-FIR数字滤波器设计与软件实现.doc

    • 资源ID:6051       资源大小:399KB        全文页数:5页
    • 资源格式: DOC        下载积分:10金币
    快捷下载 游客一键下载
    会员登录下载
    三方登录下载: 微信开放平台登录 QQ登录  
    下载资源需要10金币
    邮箱/手机:
    温馨提示:
    用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)
    支付方式: 支付宝    微信支付   
    验证码:   换一换

    加入VIP免费专享
     
    账号:
    密码:
    验证码:   换一换
      忘记密码?
        
    友情提示
    2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,就可以正常下载了。
    3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
    4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰。
    5、试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。

    数字信号处理实验5-FIR数字滤波器设计与软件实现.doc

    . . 实验五:FIR数字滤波器设计与软件实现一、实验目的:(1)掌握用窗函数法设计FIR数字滤波器的原理和方法。(2)掌握用等波纹最正确逼近法设计FIR数字滤波器的原理和方法。(3)掌握FIR滤波器的快速卷积实现原理。(4)学会调用MATLAB函数设计与实现FIR滤波器。二、实验容与步骤:(1)认真复习第七章中用窗函数法和等波纹最正确逼近法设计FIR数字滤波器的原理;(2)调用信号产生函数xtg产生具有加性噪声的信号xt,并自动显示xt与其频谱,如图1所示;图1 具有加性噪声的信号x(t)与其频谱如图(3)请设计低通滤波器,从高频噪声中提取xt中的单频调幅信号,要求信号幅频失真小于0.1dB,将噪声频谱衰减60dB。先观察xt的频谱,确定滤波器指标参数。(4)根据滤波器指标选择适宜的窗函数,计算窗函数的长度N,调用MATLAB函数fir1设计一个FIR低通滤波器。并编写程序,调用MATLAB快速卷积函数fftfilt实现对xt的滤波。绘图显示滤波器的频响特性曲线、滤波器输出信号的幅频特性图和时域波形图。(4)重复(3),滤波器指标不变,但改用等波纹最正确逼近法,调用MATLAB函数remezord和remez设计FIR数字滤波器。并比较两种设计方法设计的滤波器阶数。友情提示:MATLAB函数fir1和fftfilt的功能与其调用格式请查阅本课本;采样频率Fs=1000Hz,采样周期T=1/Fs;根据图10.6.1(b)和实验要求,可选择滤波器指标参数:通带截止频率fp=120Hz,阻带截至频率fs=150Hz,换算成数字频率,通带截止频率,通带最大衰为0.1dB,阻带截至频率,阻带最小衰为60dB。实验程序框图如图2所示。Fs=1000,T=1/Fsxt=xtg产生信号xt, 并显示xt与其频谱用窗函数法或等波纹最正确逼近法设计FIR滤波器hn对信号xt滤波:yt=fftfilt(hn,xt)1、计算并绘图显示滤波器损耗函数2、绘图显示滤波器输出信号ytEnd图2 实验程序框图3、 实验程序:1、 信号产生函数xtg程序清单:%xt=xtg(N) 产生一个长度为N,有加性高频噪声的单频调幅信号xt,采样频率Fs=1000Hz%载波频率fc=Fs/10=100Hz,调制正弦波频率f0=fc/10=10Hz.function xt=xtgN=1000;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,频率为f0ct=cos(2*pi*fc*t); %产生载波正弦波信号ct,频率为fcxt=mt.*ct; %相乘产生单频调制信号xtnt=2*rand(1,N)-1; %产生随机噪声nt%=设计高通滤波器hn,用于滤除噪声nt中的低频成分,生成高通噪声=fp=120; fs=150;Rp=0.2;As=60; % 滤波器指标fb=fp,fs;m=0,1; % 计算remezord函数所需参数f,m,devdev=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函数进行设计,用于滤除噪声nt中的低频成分yt=filter(hn,1,10*nt); %滤除随机噪声中低频成分,生成高通噪声yt%=xt=xt+yt; %噪声加信号fst=fft(xt,N);k=0:N-1;f=k/Tp;subplot(3,1,1);plot(t,xt);grid;xlabel('t/s');ylabel('x(t)');axis(0,Tp/5,min(xt),max(xt);title('(1) 信号加噪声波形')subplot(3,1,2);plot(f,abs(fst)/max(abs(fst);grid;title('(2) 信号加噪声的频谱')axis(0,Fs/2,0,1.2);xlabel('f/Hz');ylabel('幅度')2、 主程序/实验程序清单:clear all;clear all;%=调用xtg产生信号xt, xt长度N=1000,并显示xt与其频谱,=N=1000;xt=xtg;fp=120; fs=150;Rp=0.2;As=60;Fs=1000; % 输入给定指标% (1) 用窗函数法设计滤波器wc=(fp+fs)/Fs; %理想低通滤波器截止频率(关于pi归一化)B=2*pi*(fs-fp)/Fs; %过渡带宽度指标Nb=ceil(11*pi/B); %blackman窗的长度Nhn=fir1(Nb-1,wc,blackman(Nb);Hw=abs(fft(hn,1024); % 求设计的滤波器频率特性ywt=fftfilt(hn,xt,N); %调用函数fftfilt对xt滤波%以下为用窗函数法设计法的绘图部分(滤波器损耗函数,滤波器输出信号波形)f=0:1023*Fs/1024;figure(2)subplot(2,1,1)plot(f,20*log10(Hw/max(Hw);grid;title('(3) 低通滤波器幅频特性')axis(0,Fs/2,-120,20);xlabel('f/HZ');ylabel('幅度')t=0:N-1/Fs;Tp=N/Fs;subplot(2,1,2)plot(t,ywt);grid;axis(0,Tp/2,-1,1);xlabel('t/s');ylabel('y_w(t)');title('(4) 滤波噪声后的信号波形')% (2) 用等波纹最正确逼近法设计滤波器fb=fp,fs;m=1,0; % 确定remezord函数所需参数f,m,devdev=(10(Rp/20)-1)/(10(Rp/20)+1),10(-As/20);Ne,fo,mo,W=remezord(fb,m,dev,Fs); % 确定remez函数所需参数hn=remez(Ne,fo,mo,W); % 调用remez函数进行设计Hw=abs(fft(hn,1024); % 求设计的滤波器频率特性yet=fftfilt(hn,xt,N); % 调用函数fftfilt对xt滤波%以下为用等纹波设计法的绘图部分(滤波器损耗函数,滤波器输出信号波形)f=0:1023*Fs/1024;figure(3)subplot(2,1,1)plot(f,20*log10(Hw/max(Hw);grid;title('(5) 低通滤波器幅频特性')axis(0,Fs/2,-80,10);xlabel('f/HZ');ylabel('幅度')%t=0:N-1/Fs;Tp=N/Fs;subplot(2,1,2)plot(t,yet);grid;axis(0,Tp/2,-1,1);xlabel('t/s');ylabel('y_e(t)');title('(6) 滤波噪声后的信号波形')四、思考题:(1)如果给定通带截止频率和阻带截止频率以与阻带最小衰减,如何用窗函数法设计线性相位低通滤波器?请写出设计步骤。答:用窗函数法设计线性相位低通滤波器的设计步骤教材中有详细的介绍。(2) 如果要求用窗函数法设计带通滤波器,且给定通带上、下截止频率为和,阻带上、下截止频率为和,试求理想带通滤波器的截止频率。答:希望逼近的理想带通滤波器的截止频率分别为:。(3)解释为什么对同样的技术指标,用等波纹最正确逼近法设计的滤波器阶数低?用窗函数法设计的滤波器,如果在阻带截止频率附近刚好满足,则离开阻带截止频率越远,阻带衰减富裕量越大,即存在资源浪费; 几种常用的典型窗函数的通带最大衰减和阻带最小衰减固定,且差别较大,又不能分别控制。所以设计的滤波器的通带最大衰减和阻带最小衰减通常都存在较大富裕。如本实验所选的blackman窗函数,其阻带最小衰减为74dB,而指标仅为60dB。用等波纹最正确逼近法设计的滤波器,其通带和阻带均为等波纹特性,且通带最大衰减和阻带最小衰减可以分别控制,所以其指标均匀分布,没有资源浪费,所以期阶数低得多。5、 实验波形图:5 / 5

    注意事项

    本文(数字信号处理实验5-FIR数字滤波器设计与软件实现.doc)为本站会员(夺命阿水)主动上传,课桌文档仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知课桌文档(点击联系客服),我们立即给予删除!

    温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




    备案号:宁ICP备20000045号-1

    经营许可证:宁B2-20210002

    宁公网安备 64010402000986号

    课桌文档
    收起
    展开