数字信号处理第二次实验.docxVIP

第三章clf; w=-4*pi:8*pi/511:4*pi; num=[2 1];den=[1 -0.6]; h=freqz(num,den,w); subplot(2,1,1) plot(w/pi,real(h));grid title(‘H(e^{j\omega}的实部’)) xlabel(‘\omega/ \pi’); ylabel(‘振幅’)；subplot(2,1,1) plot(w/pi,imag(h));grid title(‘H(e^{j\omega}的虚部’)) xlabel(‘\omega/ \pi’); ylabel(‘振幅’)； clf; w=-pi:2*pi/255: pi;wo=0.4*pi;D=10; num=[1 2 3 4 5 6 7 8 9]; h1=freqz(num,1,w); h2=freqz([zeros(1,D) num],1,w); subplot(2,2,1)plot(w/pi,abs(h1));gridtitle(原序列的幅度谱)subplot(2,2,2)plot(w/pi,abs(h2));gridtitle(时移后序列的幅度谱)subplot(2,2,3)plot(w/pi,angle (h1));gridtitle(原序列的相位谱)subplot(2,2,4)plot(w/pi, angle (h2));gridtitl