数字信号处理matlab实验第三章浅析.docxVIP

第三章 离散时间信号的频域分析 一、实验目的 通过matlab仿真进一步研究离散时间系统的时域特性。 二、实验要求 学习并调试本章所给的例子。 回答书后给出的问题。 三、实验过程 Q.2 运行程序P3.1,求离散时间傅里叶变换的实部、虚部以及幅度和相位谱。离散时间傅里叶变换是w的周期函数吗?若是周期是多少？描述这四个图形表示的对称性。 clf; w = -4*pi:8*pi/511:4*pi; num = [2 1];den = [1 -0.6]; h = freqz(num, den, w); % Plot the DTFT subplot(2,2,1) plot(w/pi,real(h));grid title(H(e^{j\omega})的实部) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,2) plot(w/pi,imag(h));grid title(H(e^{j\omega})的虚部) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,3) plot(w/pi,abs(h));grid title(|H(e^{j\omega})|幅度谱) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,4) plot(w/pi,angle(h));grid title([H(e^{j\omega})]相位谱) xlabel(\omega /\pi); ylabel(以弧度为单位的相位); 由图得出是周期函数，周期为2，实部是和幅度谱是实偶对称，虚部和相位谱是实奇对称 Q.4修改程序P3.1，计算如下有限长序列的离散时间傅里叶变换： g[n]=[1 3 5 7 9 11 13 15 17] 并重做习题Q.2。讨论你的结果。你能解释相位谱中的跳变吗？ clf; w = -4*pi:8*pi/511:4*pi; num = [1 3 5 7 9 11 13 15 17];den = [1]; h = freqz(num, den, w); % Plot the DTFT subplot(2,2,1) plot(w/pi,real(h));grid title(H(e^{j\omega})的实部) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,2) plot(w/pi,imag(h));grid title(H(e^{j\omega})的虚部) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,3) plot(w/pi,abs(h));grid title(|H(e^{j\omega})|幅度谱) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,4) plot(w/pi,angle(h));grid title([H(e^{j\omega})]相位谱) xlabel(\omega /\pi); ylabel(以弧度为单位的相位); 1.由图得出是周期函数，周期为2，实部是和幅度谱是实偶对称，虚部和相位谱是实奇对称。 2.要计算一个系统相频特性，就要用到反正切函数，计算机中反正切函数规定，在一、二象限中的角度为0～pi，三四???限的角度为0～-pi。若一个角度从0变到2pi，但实际得到的结果是0～pi，再由-pi～0，在w=pi处发生跳变。 Q.6通过加入合适的注释语句和程序语句，修改程序P3.2，对程序生成的图形的两个轴加标记。哪个参数控制时移量？ 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));grid title(原序列的幅度谱) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,2) plot(w/pi,abs(h2));grid title(时移后序列的幅度谱) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,3) plot(w/pi,angle(h1));grid title(原序列的相位谱) xlabel(\omega /\pi); ylabel(以弧度为单位的相位); subplot(2,2,4) plot(w/pi,angle(h2));grid title(时移后序列的相位谱) xlabel(\omeg