《4.離散信號(hào)的DTFT和DFT - 數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告》由會(huì)員分享，可在線閱讀，更多相關(guān)《4.離散信號(hào)的DTFT和DFT - 數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告（5頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
1、
計(jì)算機(jī)與信息工程學(xué)院驗(yàn)證性實(shí)驗(yàn)報(bào)告
專業(yè)：通信工程 年級(jí)/班級(jí)：2011級(jí) 2013—2014學(xué)年第一學(xué)期
課程名稱
數(shù)字信號(hào)處理
指導(dǎo)教師
段新濤
本組成員
學(xué)號(hào)姓名
實(shí)驗(yàn)地點(diǎn)
計(jì)科樓111
實(shí)驗(yàn)時(shí)間
周五7-8節(jié)
項(xiàng)目名稱
離散信號(hào)的DTFT和DFT
實(shí)驗(yàn)類型
驗(yàn)證性
一、實(shí)驗(yàn)?zāi)康模?
加深對(duì)離散信號(hào)的DTFT和DFT的及其相互關(guān)系的理解。
二、實(shí)驗(yàn)原理及方法
在各種信號(hào)序列中，有限長(zhǎng)序列信號(hào)處理占有很重要地位，對(duì)有限長(zhǎng)序列，我們可以使用離散Fouier變換(DFT)。這一變換不但可以很好的反映序列的頻譜特性，而且
2、易于用快速算法在計(jì)算機(jī)上實(shí)現(xiàn)。
DTFT和DFT的主要區(qū)別就是DFT在時(shí)域和頻域都是離散的，它帶來(lái)的最大好處就是適合于數(shù)值計(jì)算，適合于計(jì)算機(jī)處理，DTFT和DFT有許多相似的性質(zhì)。
利用MATLAB工程計(jì)算語(yǔ)言按要求編寫程序算法，實(shí)現(xiàn)對(duì)有限長(zhǎng)序列的離散時(shí)間傅立葉變換（DTFT）和離散傅立葉變換（DFT）的求解。
序列x[n] 的DTFT定義：
N點(diǎn)序列x[n] 的DFT定義：
在MATLAB中，對(duì)形式為
的DTDFT可以用函數(shù)H=Freqz（num，den，w）計(jì)算；可以用函數(shù)U=fft（u，N）和u=ifft（U，N）計(jì)算N點(diǎn)序列的DFT正、反變換。
三、實(shí)驗(yàn)內(nèi)
3、容：
1、已知序列分別計(jì)算16點(diǎn)序列的16點(diǎn)和32點(diǎn)DFT，繪出幅度譜圖形，并繪出該序列的DTFT圖形。
（1）x(n) 的 16 點(diǎn)和 32 點(diǎn) DTFT，繪出 X(e jw) 幅度譜圖形;
（2）x(n) 的 16 點(diǎn)和 32 點(diǎn) DFT，繪出 X (k) 幅度譜圖形;
2、已知序列: x(n)={1，2，3，4，5，6，5，4，3，2，1}
（1）計(jì)算x(n) 的 DFT 為 X (k) ，繪出它的幅度和相位圖;
（2）計(jì)算x(n) 的 DTFT 為 ，繪出它的幅度和相位圖;
（3）利用 hold 函數(shù)，比較并驗(yàn)證 X (k) 是 的采樣。
參考流程圖：
開(kāi)始
4、
輸入序列x(n)
計(jì)算x(n)的DFT并計(jì)算其幅度和相位
計(jì)算x(n)的DFT并計(jì)算其幅度和相位
繪出圖形
開(kāi)始
四、實(shí)驗(yàn)步驟：
（1）16點(diǎn)序列X(n)的16點(diǎn)及32點(diǎn)DFT:
clc;clear all
N=16;
n=1:16;
x= sin(5*pi*n/16);
X1=fft(x,16);
X11=abs(X1);
subplot(2,1,1);
stem(X11);
xlabel('頻率');
ylabel('幅度');
title('16點(diǎn)序列x（n）的16點(diǎn)DFT');
X2=fft(x,32);
X22=abs(X2
5、);
subplot(2,1,2);
stem(X22);
xlabel('頻率');
ylabel('幅度');
title('16點(diǎn)序列x（n）的32點(diǎn)DFT');
（2）序列x（n）的DTFT:
clc;clear all
N=16;
n=1:.01:N;
x= sin(5*pi*n/16);
X1=fft(x);
X11=abs(X1);
plot(x)
xlabel('頻率');
ylabel('幅度');
title('序列x（n）的DTFT');
（3）：
clc;clear all
x=[1 2 3 4 5 6 5 4 3 2
6、 1];
nx=0:length(x)-1;
K=128;dw=2*pi/K;
k=floor((-K/2+0.5):(K/2-0.5));
X=x*exp(-j*dw*nx'*k);
subplot 311;plot(k*dw,abs(X))
hold on
xlabel('\omega');ylabel('幅度響應(yīng)');
title('11點(diǎn)序列的DTFT和FFT');
grid
Xd=fft([1 2 3 4 5 6 5 4 3 2 1]);
plot([0:length(x)-1]*2*pi/11,abs(Xd),'r.')
Xd1=fftshift(Xd);
7、
subplot 312;plot(k*dw,abs(X));
hold on
xlabel('\omega');ylabel('幅度響應(yīng)');
title('FFT移位后');
plot([-5:5]*2*pi/11,abs(Xd1),'r.');
grid
subplot 313;plot(k*dw,angle(X));
hold on
xlabel('\omega');ylabel('相位響應(yīng)');
title('FFT移位后');
grid
五、DTFT和DFT之間的區(qū)別和關(guān)系:
1、DTFT是離散時(shí)間傅里葉變換，DFT是離散傅里葉變換。
2、DTFT變
8、換后的圖形中的頻率是一般連續(xù)的（cos(wn)等這樣的特殊函數(shù)除外，其變換后是沖擊串)，而DFT是DTFT的等間隔抽樣，是離散的點(diǎn)，其函數(shù)表示為X(k）,而DTFT的函數(shù)表示為（DFT是DTFT的等間隔抽樣，DTFT變化后的頻率響應(yīng)一般是連續(xù)的，DFT變換后的頻率響應(yīng)是離散的）。
3、DTFT是以為周期的。而DFT的序列X(k）是有限長(zhǎng)的。
4、DTFT是以復(fù)指數(shù)序列{}的加權(quán)和來(lái)表示的，而DFT是等間隔抽樣，抽樣間隔為(N為離散序列的長(zhǎng)度)。
5、DTFT和DFT都能表征原序列的信息。由于現(xiàn)在計(jì)算主要使用計(jì)算機(jī)，必需要是離散的值才能參與運(yùn)算，因此在工程中DFT應(yīng)用比較廣泛，F(xiàn)FT是DFT的快速算法。
教師簽名：
年 月 日
4.離散信號(hào)的DTFT和DFT - 數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告
