数字信号处理-第三次作业.docVIP

数字信号处理作业三 时域离散系统和信号的频域表示 1-14 (a) X()=1+2 H()=2- (b) Y() =(1+2)*( 2-)=2+3-2 y(n)=2δ(n)+3δ(n-1)-2δ(n-2) (c) n=0:20 x=delta(n)+2*delta(n-1); h=2*delta(n)-delta(n-1); y=conv(x,h); m=length(y); figure; stem([0:m-1],y); xlabel(n);ylabel(y(n));title(y(n)=x(n)*h(n)); grid on; 1-15 (a) H()=++ n=-20:20; h=0.5*delta(n)+delta(n-1)+0.5*delta(n-2); H=fft(h); subplot(211),plot(n,abs(H)), title(·幅频特性);xlabel(freq(rad/s));grid on; subplot(212),plot(n,angle(H)), title(相频特性);xlabel(freq(rad/s));grid on; (b) n=-20:20; h=0.5*delta(n)+delta(n-1)+0.5*delta(n-2); x1=5*cos(n*pi*0.25); y1=conv(x1,h); figure,stem([0:length(y1)-1],y1);grid on; title(输入x(n)=5(cos)的稳态响应);xlabel(n);ylabel(y1( n)); (c) n=-20:20; h=0.5*delta(n)+delta(n-1)+0.5*delta(n-2); x2=5*cos(n*pi*0.75); y2=conv(x2,h); %m=length(y2); figure,stem([0:length(y2)-1],y2);grid on; title(输入x(n)=5(cos)的稳态响应);xlabel(n);ylabel(y1( n)); 1-20 (a) n=-20:20; h=(-0.25)*delta(n+1)+0.5*delta(n)-0.25*delta(n-1); H=fft(h); zplane(H); title(零极点图);grid on; 系统不稳定 (b) h(n)= -δ(n+1)+δ(n)-δ(n-1) 当n=-1时，h(n)= -, 即n0时，h(n)不恒为0，所以不是因果系统。 (c) H()=-+- (d) n=-20:20; h=(-0.25)*delta(n+1)+0.5*delta(n)-0.25*delta(n-1); H=fft(h); subplot(211),plot(n,abs(H)),title(幅频特性| H()|);grid on; subplot(212),plot(n,angle(H)),title(相频特性arg H());grid on; (e) 低通滤波器 1-21 | H()|=/ arg H()=arctan(4/3*sinω) 1-24 (a) H()= h(n)= (b) 冲击响应为无限长的序列,所以h(n)表示IIR滤波器。 1-25 (a) (t)=(t) X()== (b) 因为Ts=0.1T＜0.5T，所以能从X(n)恢复到xa(t)。 (c) fsmin=2×4Ω0=8Ω0 1-26 (a) H()= H()=3 y(n)=15cos() (b) H()= H()=H()=3 y(n)=15cos() 1-27 (a)系统为IIR系统。系统的冲激响应是一个无限长序列，所以是IIR系统。 (b) ω=rad/s (c) 最大幅度：50 (d) 模拟频率：500/s Ω0 10Ω0 -10Ω0 2Ω0 -Ω0 -2Ω0