数字信号处理实验报告 (实验四).docVIP

PAGE PAGE 2 实验四 离散时间信号的DTFT 一、实验目的 １. 运用MATLAB计算离散时间系统的频率响应。 ２. 运用MATLAB验证离散时间傅立叶变换的性质。 二、实验原理 （一）、计算离散时间系统的DTFT 已知一个离散时间系统，可以用MATLAB函数frequz非常方便地在给定的L个离散频率点处进行计算。由于是ω的连续函数，需要尽可能大地选取L的值（因为严格说，在MATLAB中不使用symbolic工具箱是不能分析模拟信号的，但是当采样时间间隔充分小的时候，可产生平滑的图形），以使得命令plot产生的图形和真实离散时间傅立叶变换的图形尽可能一致。在MATLAB中，freqz计算出序列{}和{}的L点离散傅立叶变换，然后对其离散傅立叶变换值相除得到。为了更加方便快速地运算，应将L的值选为2的幂，如256或者512。 例3.1 运用MATLAB画出以下系统的频率响应。 y(n)-0.6y(n-1)=2x(n)+x(n-1) 程序： 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(‘振幅’)； （二）、离散时间傅立叶变换DTFT的性质。 1．时移与频移 设 , 那么 (2.2.6) (2.2.7) 2．时域卷积定理 如果 , 那么 三、实验内容与步骤 １. 已知因果线性时不变离散时间系统 y(n)-0.4y(n-1)+0.75y(n-2)=2.2403x(n)+2.4908x(n-1)+2.2403x(n-2) 运用MATLAB画出该系统的频率响应。 clf; w=-4*pi:8*pi/511:4*pi; num=[2.2403 2.4908 2.2403];den=[1 -0.4 0.75]; 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,2) plot(w/pi,imag(h)); grid title(H(e^{j\omega}的虚部)); xlabel(\omega/ \pi); ylabel(振幅); 运行结果： ２.运行下面程序并显示它，验证离散时间傅立叶变换DTFT的时移性。 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(‘原序列的幅度谱’) subplot(2,2,2) plot(w/pi,abs(h2));grid title(‘时移后序列的幅度谱’) subplot(2,2,3) plot(w/pi,angle (h1));grid title(‘原序列的相位谱’) subplot(2,2,4) plot(w/pi, angle (h2));grid title(‘时移后序列的相位谱’) 运行结果： 3. 运行下面程序并显示它，验证离散时间傅立叶变换DTFT的频移性。 clf; w=-pi:2*pi/255