数字信号处理第三版西安科大出版高西全丁玉美课后第3和4章课题.ppt

第3章　离散傅里叶变换(DFT) 及其快速算法(FFT) 3.3.3 用DFT/FFT进行频谱分析　　对序列进行N点的DFT/FFT就是对序列频域的N点离散采样， 采样点的频率为ωk=2πk/N, k=0, 1, 2, …, N－1。 　　对信号进行频谱分析要关心三个问题： 频谱分辨率、 频谱分析范围和分析误差。 　　DFT的分辨率指的是频域采样间隔2π/N， 用DFT/FFT进行频谱分析时， 在相邻采点之间的频谱是不知道的， 因此频率分辨率是一个重要指标， 希望分辨率高， 即2π/N要小， DFT的变换区间N要大。 图3.4.2 　　显然， 由于频域采样点数MN， 不满足频域采样定理， 所以， 不能由X(k)恢复x(n)，即丢失了x(n)的频谱信息。 　［例3.4.8］ 已知序列　　　　　x(n)={1, 2, 2, 1}, h(n)={3, 2, －1, 1} （1）计算5点循环卷积y5(n)=x(n) L h(n); 　　（2）用计算循环卷积的方法计算线性卷积y(n)=x(n)*h(n)。  　　解：(1)这里是2个短序列的循环卷积计算， 可以用矩阵相乘的方法(即用教材第82页式(3.2.7)）计算， 也可以用类似于线性卷积的列表法。 因为要求5点循环卷积， 因此每个序列尾部加一个零值点， 按照教材式(3.2.7)写出 （2） 由DFT的共轭对称性可知， 如果  　　　　　x(n)=xep(n)+xop(n)且　　　　　X(k)=Re［X(k)］+j Im［X(k)］则　　Re［X(k)］=DFT［xep(n)］, j Im［X(k)］=DFT［xop(n)］所以， 当x(n)=x(N－n)时， 等价于上式中xop(n)=0, x(n)中只有xep(n)成分， 所以X(k)只有实部， 即X(k)为实函数。 又由（1）证明结果知道， 实序列的DFT必然为共轭对称函数， 即X(k)=X*(N－k)=X(N－k)， 所以X(k)实偶对称。 　　23*． 已知序列x(n)={1, 2, 3, 3, 2, 1}。 　　（1） 求出x(n)的傅里叶变换X(ejω)， 画出幅频特性和相频特性曲线(提示： 用1024点FFT近似X(ejω))； 　　（2） 计算x(n)的N（N≥6）点离散傅里叶变换X(k)， 画出幅频特性和相频特性曲线； 　　（3） 将X(ejω)和X(k)的幅频特性和相频特性曲线分别画在同一幅图中， 验证X(k)是X(ejω)的等间隔采样， 采样间隔为2π/N； 　　（4） 计算X(k)的N点IDFT， 验证DFT和IDFT的惟一性。 　　程序ex325.m如下： 　　%程序ex325.m 　　hn=［1 1 1 1］； xn=［0 1 2 3］； 　　%用DFT计算4点循环卷积yc4n： 　　H4k=fft(hn， 4)； %计算h(n)的4点DFT　　X4k=fft(xn， 4)； %计算x(n)的4点DFT　　Yc4k=H4k.*X4k； yc4n=ifft(Yc4k， 4)； 　　%用DFT计算8点循环卷积yc8n： 　　H8k=fft(hn， 8)； 　 %计算h(n)的8点DFT　　X8k=fft(xn， 8)； 　%计算x(n)的8点DFT　　Yc8k=H8k.*X8k； 　yc8n=ifft(Yc8k， 8)； 　　yn=conv(hn， xn)；　 %时域计算线性卷积yn： 题26*解图 　　x1(n)的频谱如题27*解图(a)和(b)所示， x2(n)的频谱如题27*解图(c)和(d)所示。 用32点DFT对x3(n)的谱分析结果见题27*解图(e)、 (f)和(g)， 用64点DFT对x3(n)的谱分析结果见题27*解图(h)、 (i)和(j)。 比较可知， 仅用32点分析结果就可以了。 　　请注意， x3(n)的相频特性曲线的幅度很小， 这是计算误差引起的。 实质上， x3(n)是一个实偶对称序列， 所以其理论频谱应当是一个实偶函数， 其相位应当是零。  x3anp(m)=x3n(m+10)；x3anp(N3a+1-m)=x3n(11－m)；　　end　　x3bnp=zeros(1， N3b)； 　　　　%构造x3n的周期延拓序列， 周期为N3b　　for m=1： 10， 　　 x3bnp(m)=x3n(m+10)； 　　　　x3bnp(N3b+1－m)=x