数字信号处理实验指导书（审）(6)

实验五 离散傅立叶变换DFT

4． 调用离散傅立叶正变换函数运行下面MATLAB程序，验证DFT的共轭对

称性。

x=[1 2 4 2 6 32 6 4 2 zeros(1,247)];

x1=[x(1) x(256:-1:2)];

xep=0.5*(x+x1); xop=0.5*(x-x1); XF=dft(x,256); XEPF=dft(xep,256); XOPF=dft(xop,256); clf k=0:255; subplot(2,2,1);

plot(k/128,real(XF));grid; ylabel(‘振幅’)；

title(‘原序列DFT的实部’)；

subplot(2,2,2);

plot(k/128,imag(XF));grid; ylabel(‘振幅’)；

title(‘原序列DFT的虚部’)； subplot(2,2,3);

plot(k/128,real(XEPF));grid; ylabel(‘振幅’)；

title(‘其共轭对称序列DFT’)；

subplot(2,2,4);

plot(k/128,imag(XOPF));grid; ylabel(‘振幅’)；

title(‘其共轭反对称序列DFT’)；

四、实验仪器设备

计算机，MATLAB软件

五、实验注意事项

课前预先阅读并理解实验程序；

22

实验五 离散傅立叶变换DFT

六、思考题

１.讨论实验程序1中的.^代表什么含义？

２.讨论实验程序2中圆周移位与线性移位的关系。指令mod的含义是什么？ 中极点与系统稳定性的关系？根据程序运行结果判断该系统的稳定性。 3 .计算实验程序3的理论计算结果，对比两种程序运行结果，结果是否一致？ 4 .讨论实验程序4说明实序列的共轭对称部分的离散傅立叶变换与原序列的傅立叶变换实部的关系，共轭反对称部分的离散傅立叶变换与原序列的傅立叶变换虚部的关系。

23

实验六 快速傅立叶变换FFT及其应用

实验六 快速傅立叶变换FFT及其应用

一、实验目的

１. 利用MATLAB的快速傅立叶变换来计算信号的离散傅立叶变换。 ２. 利用MATLAB程序，理解进一步离散傅立叶变换的物理意义。 3. 利用MATLAB程序，理解快速卷积算法。

二、实验原理

在MATLAB中，使用函数fft可以很容易地计算有限长序列x(n)的离散傅立叶变换X[k]。此函数有两种形式，fft(x)计算序列x(n) 的离散傅立叶变换X(k)，这里X(k)的长度与x(n)的长度相等。fft(x，L)计算序列x(n) 的L点离散傅立叶变换，其中L≥N。若L>N，在计算离散傅立叶变换之前，对x(n)尾部的L-N个值进行补零。同样，离散傅立叶变换序列X(k)的离散傅立叶逆变换x(n)用函数ifft计算，它也有两种形式。 （一）、基本序列的离散傅立叶变换计算

N点离散傅立叶变换的一种物理解释就是，X[k]是x(n)以N为周期的周期延拓序列的

~~离散傅立叶级数系数X(k)的主值区间序列，即X(k)?X(k)RN(k)。例如序列

cos(n)RN(n),当N=16时，cos(n)RN(n)正好是cos(n)的一个周期，所以

888cos(????8n)RN(n)的周期延拓序列就是这种单一频率的正弦序列。而当N=8时，

cos(n)RN(n)正好是cos(n)的半个周期，cos(n)RN(n)的周期延拓就不再是单一频888???率的正弦序列，而是含有丰富的谐波成分，其离散傅立叶级数的系数与N=16时的差别很大，因此对信号进行谱分析时，一定要截取整个周期，否则得到错误的频谱。

（二）、验证N点DFT的物理意义

假如x(n)非周期、有限长，则傅立叶变换存在，那么对X(e?k=2πk/N, k=0,1,…,,N-1取样,则可得X(k)。

24

j?)在N个等间隔频率

实验六 快速傅立叶变换FFT及其应用

X(k)?X(?)??2?k/N???x(n)en????j2?kn/N 0?k?N?1

序列x(n)的N点DFT的物理意义是对X(ω)在[0，2π]上进行N点的等间隔采样。

（三）、利用FFT计算序列的线性卷积

直接计算线性卷积计算量大，并且计算机无法判断y(n)的长度，需要计算多少的y(n)值，若输入为无限长，就更无法计算，其运算量随长度成级数增长。由于可以利用FFT对DFT进行有效的计算，我们希望能够利用DFT来计算线性卷积。

设 x(n) 和 h(n) 是长度分别为M和N的有限长序列， 令 L=M+N-1，定义两个长度L的有限长序列： x'(n)??

h'(n)???h(n),?0,0?n?N?1N?n?L?1?x(n),?0,0?n?M?1M?n?L?1 (3.4.8)

(3.4.9)

通过对x(n) 和 h(n)补充零样本值得到上面两个序列。那么：

yl(n)?x(n)?h(n)?yc(n)?x'(n) h'(n) (3.4.10) 上面的过程如下图所示：

计算线性卷积也可以直接调用函数con来计算，因为MATLAB中的计时比较粗糙，所以只有M和N较大的时候，才能比较两种方法的执行时间快慢。

25

实验六 快速傅立叶变换FFT及其应用

三、实验内容与步骤

１. 对复正弦序列x(n)?ej?8nTLAB程序求当N=16和N=8时的离散傅RN(n)，利用MA

立叶变换，并显示其图形。 程序：

N=16；N1=8； n=0:N-1;k=0: N1-1; x=exp(j*pi*n/8); X1=fft(x,N); X2=fft(x,N1); subplot(2,1,1); stem(n,abs(X1)); axis([0,20,0,20]); ylabel(‘∣X1(k)∣’); title(‘16点的DFT’)；

subplt(2,1, 2);

stem(n, abs(X2)); axis([0, 20, 0, 20]); ylabel(‘∣X2(k)∣’); title(‘8点的DFT’)；

２.已知x(n)?R4(n),X(?)?N=16时的DFT。

1?e?j4??j?1?e, 绘制相应的幅频和相频曲线，并计算N=8和

程序：

N1=8;N2=16;

n=0:N1-1;k1=0:N1-1;k2=0:N2-1; w=2*pi*(0:2047)/2048;

26

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库数字信号处理实验指导书（审）(6)在线全文阅读。