南昌大学数字信号处理实验3、.docVIP

PAGEXXX 数字信号处理实验3 一、实验目的1、理解和掌握基于MATLAB仿真研究离散时间傅里叶变换的时移性质；2、理解和掌握基于MATLAB仿真研究离散时间傅里叶变换的频移性质；3、理解和掌握基于MATLAB仿真研究离散时间傅里叶变换的卷积性质；4、理解和掌握基于MATLAB仿真研究离散时间傅里叶变换的调制性质；5、理解和掌握基于MATLAB仿真研究离散时间傅里叶变换的反转性质。 二、离散时间傅里叶变换 源程序 P3_1 如下: % Program P3_1 % Evaluation of the DTFT clf; % Compute the frequency samples of the DTFT 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,1,1) plot(w/pi,real(h));grid title(Real part of H(e^{j\omega})) xlabel(\omega /\pi); ylabel(Amplitude); subplot(2,1,2) plot(w/pi,imag(h));grid title(Imaginary part of H(e^{j\omega})) xlabel(\omega /\pi); ylabel(Amplitude); pause subplot(2,1,1) plot(w/pi,abs(h));grid title(Magnitude Spectrum |H(e^{j\omega})|) xlabel(\omega /\pi); ylabel(Amplitude); subplot(2,1,2) plot(w/pi,angle(h));grid title(Phase Spectrum arg[H(e^{j\omega})]) xlabel(\omega /\pi); ylabel(Phase in radians); 回答: Q3.1 P3_1的原始序列- H(ejw)=(2+z-1)/(1-0.6z-1) Pause命令作用 – 如果没有带参数则等待用户输入，否则按照参数来等待一定的时间 Q3.2 P3_1图像显示如下: DTFT是w的周期函数 周期 - 2*pi 四个图形表示的对称性: 实部是以2pi为周期的偶函数 虚部是以2pi为周期的奇函数 幅度谱是以2pi为周期的偶函数 相位谱是以2pi为周期的奇函数 Q3.3 修改程序P3_1计算3.3的离散时间傅里叶变换: % Program P3_1 % Evaluation of the DTFT clf; % Compute the frequency samples of the DTFT N = 512; num = [0.7 -0.5 0.3 1]; den = [1 0.3 -0.5 0.7]; [h,w] = freqz(num, den, N); % Plot the DTFT subplot(2,1,1) plot(w/pi,real(h));grid title(Real part of H(e^{j\omega})) xlabel(\omega /\pi); ylabel(Amplitude); subplot(2,1,2) plot(w/pi,imag(h));grid title(Imaginary part of H(e^{j\omega})) xlabel(\omega /\pi); ylabel(Amplitude); pause subplot(2,1,1) plot(w/pi,abs(h));grid title(Magnitude Spectrum |H(e^{j\omega})|) xlabel(\omega /\pi); ylabel(Amplitude); subplot(2,1,2) plot(w/pi,angle(h));grid title(Phase Spectrum arg[H(e^{j\omega})]) xlabel(\omega /\pi); ylabel(Phase in radians); The plots generated by running the modified Program P3_1 are shown below: DTFT是 ?的周期函数. 周期是- 2pi 产生相位谱跳变的原因是 - 计算机中反正切函数规定，在一二象限中的角度为0-pi，三四象限的角度为0-（-pi），若角度从零变到