数字信号课程设计第六组源程序加图形.docVIP

下载本文档

1
0
约1.16千字
约 4页
2017-06-08 发布于河南
举报
版权申诉

数字信号课程设计第六组源程序加图形.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

数字信号课程设计第六组源程序加图形

%题目六、假设有方波信号与正弦信号相乘并混有白噪声；设计滤波器对其进行滤波， %提出信号，对其进行频谱分析，观察其频谱，看一看这个信号的频谱有什么特点。 Fs 50000; t 0:1/Fs:0.2; y square 60*pi*t .*sin 1000*pi*t +randn 1,10001 ; yn2 square 60*pi*t .*sin 1000*pi*t ; %设计椭圆低通滤波器，根据要求，截取频率小于1000Hz的信号 %通带内最大衰减纹波Rp＝0.1db，阻带内最小衰减Rs＝60dB %通带为0～1000Hz，阻带为1200Hz至fs/2 Rp 0.1; Rs 60; Wp 1000; Ws 1200; %确定滤波器的阶数 [n,Wn] ellipord Wp/ Fs/2 ,Ws/ Fs/2 ,Rp,Rs,z ; %设计滤波器 [b1,a1] ellip n,Rp,Rs,Wn %计算设计低通滤波器的幅度频率响应 [h1,w1] freqz b1,a1,1024,Fs ; %用设计的低通滤波器对信号y进行滤波得到信号yn1 yn1 filter b1,a1,y ; figure; subplot 4,1,1 ; plot w1 1:200 ,abs h1 1:200 ; grid; xlabel 频率/Hz ; ylabel 幅值 ; title 低通滤波器 ; subplot 4,1,2 ; plot t 1:1000 ,yn2 1:1000 ; xlabel 频率/Hz ; ylabel 幅值 ; title 不加噪声的信号 ; subplot 4,1,3 ; plot t 1:1000 ,y 1:1000 ; xlabel 频率/Hz ; ylabel 幅值 ; title 加噪声的信号 ; subplot 4,1,4 ; plot t 1:1000 ,yn1 1:1000 ; xlabel 频率/Hz ; ylabel 幅值 ; title 提取信号 ; %对滤波后的信号进行频谱分析， %验证对信号的衰减是否达到要求 figure; subplot 2,1,1 Xk1 fft y 400:1000 ; M abs Xk1 plot M ; title 原信号频谱 ; xlabel 频率/Hz ; ylabel 幅度 ; subplot 2,1,2 Xk2 fft yn1 400:1000 ; N abs Xk2 plot N ; title 滤波后信号频谱 ; xlabel 频率/Hz ; ylabel 幅度 ;