数字信号处理matlab实验报告材料.docxVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
___计算机__学院_网络工程__专业___3__班 学号_3116005027__ ________铠光________ 教师评定_________________ 实验题目___________数字信号处理实验_____________ 实验三、FIR数字滤波器设计 一、实验目的 1、熟悉FIR滤波器的设计过程。 2、能够根据滤波器设计指标设计不同的滤波器。 3、熟悉FDATool,SPTool的一般操作。 二、实验环境 Matlab信号处理工具箱。 三、实验容 1、用模拟原型法设计法与窗函数法分别设计一个IIR和FIR低通滤波器F1, F2,并满足性能指标:通带归一化频率wp=0.5,阻带归一化频率ws=0.6,阻带衰减不小于40dB,通带波纹不大于3dB。假设一个信号S,其中f1=8Hz,f2=21Hz。信号采样频率为50Hz。试将原信号与通过滤波器的信号进行比较。 要求:画出滤波器的幅频相应图,滤波前后对比图。 wp=0.5*pi;ws=0.6*pi;rs=40;rp=3; T=0.02;Nn=128; %采样间隔 wp=wp/T;ws=ws/T; [N,Wn]=buttord(wp,ws,rp,rs,s); %计算模拟滤波器的最小阶数 [z,p,k]=buttap(N); %设计低通原型数字滤波器 [Bap,Aap]=zp2tf(z,p,k); %零点极点增益形式转换为传递函数形式 [b,a]=lp2lp(Bap,Aap,Wn);%低通滤波器的频率转换 [bz,az]=impinvar(b,a,1/T);%脉冲响应不变法设计数字滤波器传递函数 figure(1) [H,f]=freqz(bz,az,Nn,1/T);%输出幅频响应和相频响应 subplot(2,1,1),plot(f,20*log10(abs(H))) xlabel(频率/Hz);ylabel(振幅/dB);grid on; subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel(频率/Hz);ylabel(相位/^o);grid on; figure(2) f1=8;f2=21;%输入信号含有的频率 N=100;%数据点数 dt=T;n=0:N-1;t=n*dt; %时间序列 x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t);%输入信号 subplot(2,1,1),plot(t,x),title(输入信号) y=filtfilt(bz,az,x);%对信号进行滤波 subplot(2,1,2),plot(t,y),title(输出信号),xlabel(时间/s) wp=0.5*pi;ws=0.6*pi;%滤波器边界频率 wdelta=ws-wp;%过渡带宽 N=ceil(8*pi/wdelta)%计算滤波器所用窗函数最小长度 Nw=N; wc=(wp+ws)/2; %截止频率在同带和阻带边界的中点 n=0:N-1; alpha=(N-1)/2; %求滤波器的相位延迟 m=n-alpha+eps; %eps为MATLAB系统的精度 hd=sin(wc*m)./(pi*m); %求理想滤波器脉冲响应 win=hanning(Nw); %采用汉宁窗 h=hd.*win; %在时间域乘积对应于频率域的卷积 b=h; figure(1) [H,f]=freqz(b,1,512,50);%采用50Hz的采样频率画出该滤波器的幅频和相频响应 subplot(2,1,1),plot(f,20*log10(abs(H))) xlabel(频率/Hz);ylabel(振幅/dB);grid on; f1=8;f2=21; %检测输入信号含有两种频率成分 dt=0.02;t=0:dt:3; %采样间隔和检测信号的时间序列 x=sin(2*pi*f1*t)+cos(2*pi*f2*t); %检测信号 y=filter(b,1,x); %给出滤波器的输出 figure(2) subplot(2,1,1),plot(t,x),title(输入信号) subplot(2,1,2),plot(t,y) hold on;plot([1 1]*(N-1)/2*dt,ylim,r) xlabel(时间/s),title(输出信号); 2、在Workspace生成上面的原始信号S(即8Hz与21Hz组成的信号),利用SPTool导入该信号,调用FDATool设计题1的滤波器F。用滤波器F对信号F做滤波处理。 要求:记录用SPTool及FDATool设计滤波器,滤波,显示滤波效果的全过程。 _计算机__学院 网络工程 专业_3___班 学号 3116005027 铠光

文档评论(0)

linlin921 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档