實(shí)驗(yàn)一--序列、頻譜、DFT的性質(zhì)教學(xué)教材.doc

1、精品文檔專業(yè)：_信息與通信工程姓名： _ _實(shí)驗(yàn)報(bào)告學(xué)號(hào) _日期： _地點(diǎn) _課程名稱：數(shù)字信號(hào)處理指導(dǎo)老師：成績(jī)： _實(shí)驗(yàn)名稱：FIR 序列、頻譜、 DFT 的性質(zhì) 實(shí)驗(yàn)類型： _演示 _同組學(xué)生姓名： 一、實(shí)驗(yàn)?zāi)康暮鸵笤O(shè)計(jì)通過(guò)演示實(shí)驗(yàn)，建立對(duì)典型信號(hào)及其頻譜的直觀認(rèn)識(shí)， 理解 DFT 的物理意義、主要性質(zhì)。二、實(shí)驗(yàn)內(nèi)容和步驟2-1 用 MATLAB ，計(jì)算得到五種共9 個(gè)序列：n0nlength1例如， a=0.5, length=102-1-1 實(shí)指數(shù)序列 x(n)a0otherwisea=0.9, length=10a=0.9, length=202-1-2 復(fù)指數(shù)序列 x(n)(a

2、jb ) n0 nlength 1 例如， a=0.5, b=0.8, length=100otherwise2-1-3 從正弦信號(hào) x(t)=sin(2 ft+delta)抽樣得到的正弦序列 x(n)=sin(2 fnT+delta)。如，信號(hào)頻率 f=1Hz，初始相位 delta=0，抽樣間隔 T=0.1 秒，序列長(zhǎng)length=10。2-1-4 從余弦信號(hào) x(t)=cos(2 ft + delta)抽樣得到的余弦序列 x(n)=cos(2 fnT + delta)。如，信號(hào)頻率 f=1Hz，初相位 delta=0，抽樣間隔 T=0.1 秒，序列長(zhǎng) length=10。2-1-5 含兩個(gè)

3、頻率分量的復(fù)合函數(shù)序列x(n)=sin(2 f1nT)+deltasin(2 f 2nT+phi)。如，頻率 f1頻率 f2相對(duì)振幅初相位 phi抽樣間隔 T序列長(zhǎng)(Hz)(Hz)delta(度 )(秒)length130.500.110130.5900.110130.51800.1102-2用 MATLAB ，對(duì)上述各個(gè)序列，重復(fù)下列過(guò)程。2-2-1 畫(huà)出一個(gè)序列的實(shí)部、虛部、模、相角；觀察并記錄實(shí)部、虛部、模、相角的精品文檔精品文檔實(shí)驗(yàn)名稱： FIR 序列、頻譜、 DFT 的性質(zhì) 姓名： _特征。2-2-2 計(jì)算該序列的幅度譜、頻譜實(shí)部、頻譜虛部；觀察和并記錄它們的特征，給予解釋。2-2-

4、3 觀察同種序列取不同參數(shù)時(shí)的頻譜，發(fā)現(xiàn)它們的差異，給予解釋。三、主要儀器設(shè)備MATLAB編程。四、操作方法和實(shí)驗(yàn)步驟（參見(jiàn)“二、實(shí)驗(yàn)內(nèi)容和步驟”）五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理列出 MATLAB程序清單，加注釋。2-1-1a (a=0.5, length=10)程序n=0:9;xn=(0.5).n).*(0=n&n=9);xw=dftmtx(10)*xn;%用 DFT 求頻譜f=n/10.*(0=n&n=5)+(10-n)/10.*(6=n&n=9);%求出對(duì)應(yīng)頻率figure(1);%畫(huà)出序列的實(shí)部、虛部、模、相角subplot(2,2,1);stem(n,real(xn);xlabel(n);y

5、label(real(xn);subplot(2,2,2);stem(n,imag(xn);xlabel(n);ylabel(imag(xn);subplot(2,2,3);stem(n,abs(xn);xlabel(n);ylabel(abs(xn);subplot(2,2,4);stem(n,angle(xn);xlabel(n);ylabel(angle(xn);figure(2);%畫(huà)出序列的幅度譜、頻譜實(shí)部、頻譜虛部subplot(3,1,1);stem(f,abs(xw);xlabel(f/Hz);ylabel(abs(xw);實(shí)驗(yàn)名稱： FIR 序列、頻譜、 DFT 的性質(zhì) 姓名

6、： _ 3精品文檔精品文檔subplot(3,1,2);stem(f,real(xw);xlabel(f/Hz);ylabel(real(xw);subplot(3,1,3);stem(f,imag(xw);xlabel(f/Hz);ylabel(imag(xw);2-1-1b（a=0.9, length=10）程序n=0:9;xn=(0.9).n).*(0=n&n=9);xw=dftmtx(10)*xn;%用 DFT 求頻譜f=n/10.*(0=n&n=5)+(10-n)/10.*(6=n&n=9);%求出對(duì)應(yīng)頻率figure(1);%畫(huà)出序列的實(shí)部、虛部、模、相角subplot(2,2,1

7、);stem(n,real(xn);xlabel(n);ylabel(real(xn);subplot(2,2,2);stem(n,imag(xn);xlabel(n);ylabel(imag(xn);subplot(2,2,3);stem(n,abs(xn);xlabel(n);ylabel(abs(xn);subplot(2,2,4);stem(n,angle(xn);xlabel(n);ylabel(angle(xn);figure(2); %畫(huà)出序列的幅度譜、頻譜實(shí)部、頻譜虛部 subplot(3,1,1);stem(f,abs(xw); xlabel(f/Hz);ylabel(abs

8、(xw); subplot(3,1,2);stem(f,real(xw); xlabel(f/Hz);ylabel(real(xw); subplot(3,1,3);stem(f,imag(xw); xlabel(f/Hz);ylabel(imag(xw);2-1-1cn=0:19;xn=(0.9).n).*(0=n&n=19);實(shí)驗(yàn)名稱： FIR 序列、頻譜、 DFT 的性質(zhì) 姓名： _邵振江 _學(xué)號(hào) _3080102350_P. 4精品文檔精品文檔xw=dftmtx(20)*xn;%用 DFT 求頻譜f=n/20.*(0=n&n=10)+(20-n)/20.*(11=n&n=0&n=9);

9、f=n/10.*(0=n&n=5)+(10-n)/10.*(6=n&n=0&n=9);xw=dftmtx(10)*xn;%用 DFT 求頻譜f=n.*(0=n&n=5)+(10-n).*(6=n&n=0&n=9);xw=dftmtx(10)*xn;%用 DFT 求頻譜f=n.*(0=n&n=5)+(10-n).*(6=n&n=0&n=9);xw=dftmtx(10)*xn;%用 DFT 求頻譜f=n.*(0=n&n=5)+(10-n).*(6=n&n=0&n=9);xw=dftmtx(10)*xn;%用 DFT 求頻譜f=n.*(0=n&n=5)+(10-n).*(6=n&n=0&n=9);x

10、w=dftmtx(10)*xn;%用 DFT 求頻譜f=n.*(0=n&n=5)+(10-n).*(6=n&n=9);%求出對(duì)應(yīng)頻率figure(1);%畫(huà)出序列的實(shí)部、虛部、模、相角subplot(2,2,1);stem(n,real(xn);xlabel(n);ylabel(real(xn);subplot(2,2,2);stem(n,imag(xn);xlabel(n);ylabel(imag(xn);subplot(2,2,3);stem(n,abs(xn);xlabel(n);ylabel(abs(xn);subplot(2,2,4);stem(n,angle(xn);xlabel(

11、n);ylabel(angle(xn);w_begin=0;w_step=pi/1600;w_end=2*pi;figure(2);%畫(huà)出序列的幅度譜、頻譜實(shí)部、頻譜虛部subplot(3,1,1);stem(f,abs(xw);xlabel(f/Hz);ylabel(abs(xw);subplot(3,1,2);stem(f,real(xw);xlabel(f/Hz);ylabel(real(xw);subplot(3,1,3);stem(f,imag(xw);xlabel(f/Hz);ylabel(imag(xw);實(shí)驗(yàn)名稱： FIR 序列、頻譜、 DFT 的性質(zhì) 姓名： _ _P. 9精

12、品文檔精品文檔六、實(shí)驗(yàn)結(jié)果與分析觀察實(shí)驗(yàn)結(jié)果（數(shù)據(jù)及圖形）的特征，做必要的記錄，做出解釋。包括：6-1各種序列的圖形（時(shí)域）和頻譜（頻域）各有何特征，給予解釋。6-2DFT 物理意義。 X(0) 、X(1) 和 X(N 1)的物理意義。6-3DFT 的主要性質(zhì)。實(shí)驗(yàn)結(jié)果：2-1-1ax(n)0.5n0 n 90otherwise11)2wx0.5(1)sbnnaxx(0.5g0l0aaem0123456789rik-0.5)20-1510510w00x(1nnlaer1101234567890)0.5k)1nnx)x(0.5e0wslbgxn(ag0a-0.5ami-10510-1510012

13、345678900knn序列的實(shí)部、虛部、模、相角序列的 DFT 結(jié)果)wx(sba)w(xlare21000.050.10.150.20.250.30.350.40.450.5f/Hz21000.050.10.150.20.250.30.350.40.450.5f/Hz)wx(gami10-100.050.10.150.20.250.30.350.40.450.5f/Hz序列的頻譜實(shí)驗(yàn)名稱： FIR 序列、頻譜、 DFT 的性質(zhì) 姓名： _ _P. 102-1-1b0.9n0 n 9x(n)otherwise0精品文檔精品文檔11)10wx)0.5(5)snbnxax(0.5g0laa0em

14、r0123456789i-0.5k100-1)05100510wx(nnl5aer1101234567890)0.5k)n1nx)x(e(0.50wslbgxan(0ag-0.5ami-10-1012345678905100510knn序列的實(shí)部、虛部、模、相角序列的 DFT 結(jié)果)10wx(5sba000.050.10.150.20.250.30.350.40.450.5f/Hz)10wx(l5aer000.050.10.150.20.250.30.350.40.450.5f/Hz)1wx(g0ami-100.050.10.150.20.250.30.350.40.450.5f/Hz序列的頻

15、譜n0 n1 92-1-1cx(n)0. 90otherwise)1011wx(5s0.5b)annx0x(l0.5g002468101214161820aaemkri-0.510)0-1wx(0510152005101520l5anner110246810121416182000.5k)5n)nxx(w(0.5e0slxbg(0ngaaa-0.5mi-50-1024681012141618200510152005101520knn實(shí)驗(yàn)名稱： FIR 序列、頻譜、 DFT 的性質(zhì) 姓名： _ _P. 11精品文檔精品文檔)10wx(5sba00.050.10.150.20.250.30.350.40.450.50f/Hz)10wx(5laer000.050.10.150.20.250.30.350.40.450.5f/Hz)5wx(g0ami-50.050.10.150.20.250.30.350.40.450.50f/Hz序列的頻譜觀察以上三個(gè)序列， 發(fā)現(xiàn)它們都為正的實(shí)序列，所以序列的虛部和相角都為零。觀察它