数字信号处理实验报告三.doc

下载文档

5
0
约1.43千字
约 8页
2017-02-24 发布于江西
举报
版权申诉
保障服务

数字信号处理实验报告三.doc

1、本文档共8页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

Q5.5 通过将正弦信号频率变为3Hz和7Hz，重做习题P5.1。相应的等效离散时间信号与之前的离散时间信号有差别嘛？若没有，为什么？　　　clf;　　　 t=0:0.0005:1;　　　 f=13;　　　 xa=cos(2*pi*f*t);　　　 subplot(2,1,1)　　　 plot(t,xa);　 grid　　 axis([0 1 -1.2 1.2]);　 subplot(2,1,2)　 T=0.1;　　 n=0:T:1;　 xs=cos(2*pi*f*n);　　 k=0:length(n)-1;　　 stem(k,xs)　　 grid　　　 axis([0 (length(n)-1) -1.2 1.2]);　　　　　　　　　　　　　　　　　　　　　　　　clf;　　 T=0.1;f=13;　　　 n=(0:T:1);　 xs=cos(2*pi*f*n);　 t=linspace(-0.5,1.5,500);　　 ya=sinc((1/T)*t(:,ones(size(n)))-(1/T)*n(:,ones(size(t))))*xs;　 plot(n,xs,o,t,ya);grid;　　　　　　　　axis([0 1 -1.2 1.2]);　　　　　　　　　　　　　　　　clf;　　　 t=0:0.005:10;　 xa=2*t.*exp(-t);　　 subplot(2,2,1)　 plot(t,xa);　 grid　　　　　　　subplot(2,2,2);　　 wa=0:10/511:10;　　　 ha=freqs(2,[1 2 1],wa);　　 plot(wa/(2*pi),abs(ha));　 grid　　 title(|x_{a}(j/Omega)|);　 axis([0 5/pi 0 2]);　　 subplot(2,2,3)　　 T=1;　　　 n=0:T:10;;　　 xs=2*n.*exp(-n);　 k=0:length(n)-1;　　 stem(k,xs);　 grid　　　　　　　　　subplot(2,2,4);　 wd=0:pi/225:pi;　　　 hd=freqs(xs,1,wd);　　　 plot(wd/(T*pi),T*abs(hd));grid;　　　　　　title(|X(e^{j/Omega})|);　　 axis([0 1/T 0 2]);　　　　 Q5.9 连续时间函数xa（t）是什么？xa（t）的连续时间傅里叶变换是如何计算的？　　　xa=2*t.*exp(-t)，将连续时间信号xa(t)转换成为等效的离散时间信号x[n]，然后进行傅里叶变换。。　　　　　　　　　　　　clf;　　　 Fp=3500;Fs=4500;　　　 Wp=2*pi*Fp;Ws=2*pi*Fs;　　　 [N,Wn]=buttord(Wp,Ws,0.5,30,s);　　 [b,a]=butter(N,Wn,s);　　　 wa=0:(3*Ws)/511:3*Ws;　　　 h=freqs(b,a,wa);　　 plot(wa/(2*pi),20*log10(abs(h)));grid　　 axis([0 3*Fs -60 5]); 　　　　　　　　　　　　　　　　　　　