资源描述:
《数字信号处理实验归纳总结报告三.doc》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、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 xlabel('时间,msec');ylabel('振幅'); title('连续时间信号x_{a}(t)'); axis([01-1.21.2]); subplot(2,1,2) T=0.1;
2、 n=0:T:1; xs=cos(2*pi*f*n); k=0:length(n)-1; stem(k,xs) grid xlabel('时间序列n');ylabel('振幅'); title('离散时间信号x[n])'); axis([0(length(n)-1)-1.21.2]); Q5.8恢复原始显示范围并通过分别改变正弦信号的频率到3Hz和7Hz,重做P5.2相应的等效离散时间信号与原来的离散时间信号有差别吗?若没有,
3、为什么? 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; xlabel('时间,msec');ylabel('振幅'); title('重构连续时间信号y_{a}(t)'); axis([01-1.21
4、.2]); clf; t=0:0.005:10; xa=2*t.*exp(-t); subplot(2,2,1) plot(t,xa); grid xlabel('时间,msec');ylabel('振幅'); title('连续时间信号x_{a}(t)'); subplot(2,2,2); wa=0:10/511:10; ha=freqs(2,[121],wa); plot(wa/(2*pi),abs(ha)); grid
5、 xlabel('时间序列n');ylabel('振幅'); title('
6、x_{a}(j/Omega)
7、'); axis([05/pi02]); subplot(2,2,3) T=1; n=0:T:10;; xs=2*n.*exp(-n); k=0:length(n)-1; stem(k,xs); grid xlabel('时间序列n');ylabel('振幅'); title('离散时间信号x[n]'); subplot(2,2,4); wd=0:pi/225:pi;
8、 hd=freqs(xs,1,wd); plot(wd/(T*pi),T*abs(hd));grid; xlabel('频率,kHz');ylabel('振幅'); title('
9、X(e^{j/Omega})
10、'); axis([01/T02]); Q5.9连续时间函数xa(t)是什么?xa(t)的连续时间傅里叶变换是如何计算的? 答:答:xa=2*t.*exp(-t),将连续时间信号xa(t)转换成为等效的离散时间信号x[n],然后进行傅里叶变换。。 Q
11、5.10运行P5.3,产生并显示离散时间信号及其连续时间等效,以及他们各自的傅里叶变换。有任何明显的混叠影响吗? 答:没有明显的混叠现象。 Q5.11将抽样周期增加到1.5,重复程序P5.3。有任何明显的混叠影响吗? Q5.13在程序P5.4中,通带波纹Rp和最小阻带衰减Rs是多少dB?通带和阻带边界频率是多少Hz? clf; Fp=3500;Fs=4500; Wp=2*pi*Fp;Ws=2*pi*Fs; [N,Wn]=buttord(Wp,Ws,0.5,30,'s'
12、); [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 xlabel('频率,Hz');ylabel('增益,dB'); title('增益响应'); axis([03*Fs-605]); Q5.14运行程序P5.4并增益响应。所设计的滤波器满足给定的指标吗?所设计的滤波器的阶数N和单