实验四离散傅里叶变换

《实验四离散傅里叶变换》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验四：离散傅里叶变换实验原理：DFT的快速算法FFT利用了的三个固有特性：(1)对称性(2)周期性(3)可约性。FFT算法基本上可以分为两大类，即按时间抽选法(D1T,Decimation-ln-Time)和按频率抽选法(DIF,Decimation-ln-frequency)。MATLAB中提供了进行快速傅里叶变换的fft函数：X=fft(x),基2时间抽取FFT算法，x是表示离散信号的向量；X是系数向量；X=fft(x,N),补零或截断的N点DFT,当x得长度小于N时，对补零使其长度为N,当x的长度大于N时，对x截断使其长度为N。实验内容：l.

2、k=60;n=[0:1:k/2];xal=2*sin(10*pi*n/k)+cos(18*pi*n/k);subplot(321)stem(n,xal)xlabel('N*);ylabel(*x(n)*);xkl=fft(xal);xkl=abs(xkl)subplot(322)stem(n,xkl)xlabel(*k*);ylabel(，X(k)*);n=[0:l:k*1.5];xal=2*sin(10*pi*n/k)+cos(18*pi*n/k);subplot(323)stem(n,xal)xlabel(W);ylabel(’x(n)’)；xk

3、l=fft(xal);xkl=abs(xkl)subplot(324)stem(n,xkl)xlabel('k*);ylabel(’X(k)*);n=[0:1:k*2];xal=2*sin(10*pi*n/k)+cos(18*pi*n/k);subplot(325)stem(n,xal)xlabel('N*);ylabel(*x(n)*);xkl=fft(xal);xkl=abs(xkl)subplot(326)stem(n,xkl)xlabel(rkf);ylabel(1X(k)’)；N2*Xr642X01020X40SO60708090结果分析：

4、当选取的采样分析长度不为周期信号周期T的整数倍时，就会发生频率泄露，如上阁。反之则不会发生频率泄露2.n=[0:7];L=length(n);x=ones(L,1);N=32;nl=[0:N-1];xk=fft(x,N);subplot(221);n2=(2*pi/N)*[0:(M_l)]，；stem(n2,abs(xk));xlabel(1k1),ylabel(1幅度1);title(=32*);subplot(223);stem(n2,angle(xk));xlabel(1k’)，ylabel(1相位1);N=64;nl=[0:N-l];xk=f

5、ft(x,N);subplot(222);n2=(2*pi/N)*[0:(N-l)]1;stem(n2,abs(xk));xlabel(▼k’)，ylabel(1幅度1);title(*N=641);subplot(224);stem(n2,angle(xk));xlabel(1k’)，ylabel(!相位•)；结果分析：当采样点数越人时，用STEM绘制出的图形包络线越接近于真实的连续幅频和相频特性的连续曲线。