離散傅里葉變換(DFT)快速傅里葉變換(FFT)-實(shí)驗(yàn)四

1、精選優(yōu)質(zhì)文檔-傾情為你奉上課程名稱(chēng)：數(shù)字信號(hào)處理 實(shí)驗(yàn)成績(jī)： 指導(dǎo)教師： 實(shí) 驗(yàn) 報(bào) 告院系： 信息工程學(xué)院 班級(jí)： 學(xué)號(hào)： 姓名： 日期： 2011. 10.30 實(shí)驗(yàn)四：DFS、DFT與FFT1、已知某周期序列的主值序列為x(n)=0,1,2,3,2,1,0，編程顯示2個(gè)周期的序列波形。要求： 用傅里葉級(jí)數(shù)求信號(hào)的幅度譜和相位譜，并畫(huà)出圖形 求傅里葉級(jí)數(shù)逆變換的圖形，并與原序列進(jìn)行比較程序清單： N=7; xn=0,1,2,3,2,1,0; xn=xn,xn; n=0:2*N-1; k=0:2*N-1; Xk=xn*exp(-1i*2*pi/N).(n'*k); x=(Xk*exp

2、(1i*2*pi/N).(n'*k)/N; subplot(2,2,1);stem(n,xn); title('x(n)');axis(-1,2*N,1.1*min(xn),1.1*max(xn); subplot(2,2,2);stem(n,abs(x); title('IDFS|X(k)|');axis(-1,2*N,1.1*min(x),1.1*max(x); subplot(2,2,3),stem(k,abs(Xk); title('|X(k)|');axis(-1,2*N,1.1*min(abs(Xk),1.1*max(abs(

3、Xk); subplot(2,2,4),stem(k,angle(Xk); title('arg|X(k)|');axis(-1,2*N,1.1*min(angle(Xk),1.1*max(angle(Xk);程序運(yùn)行結(jié)果如下圖：2、已知有限長(zhǎng)序列x(n)=1,0.5,0,0.5,1,1,0.5,0，要求： 求該序列的DFT、IDFT的圖形；程序清單： xn=1,0.5,0,0.5,1,1,0.5,0; N=length(xn); n=0:N-1; k=0:N-1; Xk=xn*exp(-1i*2*pi/N).(n'*k); x=(Xk*exp(1i*2*pi/N).(

4、n'*k)/N; subplot(2,2,1);stem(n,xn); title('x(n)');axis(-1,N,1.1*min(xn),1.1*max(xn); subplot(2,2,2);stem(n,abs(x); title('IDFT|X(k)|');axis(-1,N,1.1*min(x),1.1*max(x); subplot(2,2,3),stem(k,abs(Xk); title('|X(k)|');axis(-1,N,1.1*min(abs(Xk),1.1*max(abs(Xk); subplot(2,2,4)

5、,stem(k,angle(Xk); title('arg|X(k)|');axis(-1,N,1.1*min(angle(Xk),1.1*max(angle(Xk);程序運(yùn)行結(jié)果如下圖： 用FFT算法求該序列的DFT、IDFT的圖形；程序清單： xn=1,0.5,0,0.5,1,1,0.5,0; N=length(xn); subplot(2,2,1);stem(n,xn); title('x(n)'); k=0:N-1; Xk=fft(xn,N); subplot(2,1,2);stem(k,abs(Xk); title('Xk=DFT(xn)

6、9;); xn1=ifft(Xk,N); subplot(2,2,2);stem(n,xn1);title('x(n)=IDFT(Xk)程序運(yùn)行結(jié)果如下圖： 假定采用頻率Fs=20Hz，序列長(zhǎng)度N分別取8、32和64，用FFT計(jì)算其幅度譜和相位譜。程序清單： Ts=0.05;C=8,32,64; for r=0:2; N=C(r+1); n=0:N-1; xn=1,0.5,0,0.5,1,1,0.5,0,zeros(1,N-8); D=2*pi/(N*Ts); k=floor(-(N-1)/2:(N-1)/2); X=fftshift(fft(xn,N); subplot(3,2,2*

7、r+1);plot(k*D,abs(X); subplot(3,2,2*r+2);stairs(k*D,angle(X); end程序運(yùn)行結(jié)果如下圖：3、 已知一個(gè)無(wú)限長(zhǎng)序列x(n)=0.5n（n0），采樣周期Ts=0.2s，要求序列長(zhǎng)度分別取8、32和64，用FFT求其頻譜。程序清單： Ts=0.2;C=8,32,64; for r=0:2; N=C(r+1); n=0:N-1; xn=0.5.n; D=2*pi/(N*Ts); k=floor(-(N-1)/2:(N-1)/2); X=fftshift(fft(xn,N); subplot(3,2,2*r+1);plot(k*D,abs(X

8、); axis(-80,80,0,3); subplot(3,2,2*r+2);stairs(k*D,angle(X); axis(-80,80,-1,1); end 程序運(yùn)行結(jié)果如下圖： 4. 思考題： 離散傅里葉級(jí)數(shù)與連續(xù)性周期信號(hào)的傅里葉級(jí)數(shù)有何不同？周期序列的頻譜有何特點(diǎn)？答:傅里葉級(jí)數(shù)僅有N個(gè)獨(dú)立的諧波分量，展成離散傅立葉級(jí)數(shù)時(shí)，只能取k=0N-1的N個(gè)獨(dú)立的諧波分量，而連續(xù)性周期信號(hào)的傅里葉級(jí)數(shù)由無(wú)窮多個(gè)與基波頻率成整數(shù)倍的諧波分量疊加而成。周期序列的頻譜特點(diǎn)：以N為周期的周期序列的頻譜也是一個(gè)以N為周期的周期序列。 DFS、DFT、FFT有何聯(lián)系？答:DFS是求周期序列（周期為N)的離散傅里葉級(jí)數(shù)，也是一個(gè)以N為周期的周期序列，DFT是求有限長(zhǎng)序列的離散傅立葉變換X(k)，而X(k)是的主值序列，也就是說(shuō)對(duì)DFS的變換結(jié)果只取0N-1的序列值就得到DFT的變換。FFT只是DFT的一種快速高效的算法，它適用序列長(zhǎng)度很大時(shí)，使用其可以大大縮短運(yùn)算時(shí)間。5. 實(shí)驗(yàn)總結(jié)： 答：通過(guò)本次實(shí)驗(yàn)，我知道了對(duì)有限長(zhǎng)的序列求DFT時(shí)一定注意序列長(zhǎng)度的取值范圍；矩陣的求冪運(yùn)算與數(shù)組的求冪運(yùn)算的表達(dá)格式是不同的，前者為:an,后者為：a.n,另外參與矩陣冪運(yùn)算的量必須有一個(gè)是標(biāo)量，底數(shù)矩陣必須是方陣，才能反復(fù)相乘；進(jìn)行頻