dsp實(shí)驗(yàn)1有限長序列、頻譜、dft的性質(zhì)

1、課程名稱：數(shù)字信號(hào)處理實(shí)驗(yàn)名稱：有限長序列、頻譜、DFT的性質(zhì)一、實(shí)驗(yàn)?zāi)康暮鸵笠弧?shí)驗(yàn)?zāi)康暮鸵笤O(shè)計(jì)通過演示實(shí)驗(yàn)，建立對(duì)典型信號(hào)及其頻譜的直觀認(rèn)識(shí)，理解DFT的物理意義、主要性質(zhì)。二、實(shí)驗(yàn)內(nèi)容和步驟二、實(shí)驗(yàn)內(nèi)容和步驟21用MATLAB，計(jì)算得到五種共9個(gè)序列：211實(shí)指數(shù)序列???????otherwiselengthnanxn010)(例如，a=0.5length=10a=0.9length=10a=0.9length=20212復(fù)

2、指數(shù)序列例如，a=0.5b=0.8length=10????????otherwiselengthnjbanxn010)()(213從正弦信號(hào)x(t)=sin(2?ftdelta)抽樣得到的正弦序列x(n)=sin(2?fnTdelta)。如，信號(hào)頻率f=1Hz，初始相位delta=0，抽樣間隔T=0.1秒，序列長length=10。214從余弦信號(hào)x(t)=cos(2?ftdelta)抽樣得到的余弦序列x(n)=cos(2?fnTde

3、lta)。如，信號(hào)頻率f=1Hz，初相位delta=0，抽樣間隔T=0.1秒，序列長length=10。215含兩個(gè)頻率分量的復(fù)合函數(shù)序列x(n)=sin(2?f1nT)deltasin(2?f2nTphi)。如，頻率f1(Hz)頻率f2(Hz)相對(duì)振幅delta初相位phi(度)抽樣間隔T(秒)序列長length130.500.110130.5900.110130.51800.11022用MATLAB，對(duì)上述各個(gè)序列，重復(fù)下列過程。2

4、21畫出一個(gè)序列的實(shí)部、虛部、模、相角；觀察并記錄實(shí)部、虛部、模、相角的特征。222計(jì)算該序列的幅度譜、頻譜實(shí)部、頻譜虛部；觀察和并記錄它們的特征，給予解釋。223觀察同種序列取不同參數(shù)時(shí)的頻譜，發(fā)現(xiàn)它們的差異，給予解釋。subplot(223)stem(nabs(xn)filled)xlabel(n)ylabel(abs(xn))title(xn模)subplot(224)stem(nangle(xn)filled)xlabel(n)

5、ylabel(angle(xn))title(xn相角)figure(2)%畫出頻譜的幅度譜、頻譜實(shí)部、頻譜虛部subplot(311)stem(kabs(xw)filled)xlabel(k)ylabel(abs(xw))title(幅度譜)subplot(312)stem(kreal(xw)filled)xlabel(k)ylabel(real(xw))title(頻譜實(shí)部)subplot(313)stem(kimag(xw)fil

6、led)xlabel(k)ylabel(imag(xw))title(頻譜虛部)13.a=0.9length=20clcclfclear%清除n=0:19%設(shè)置自變量區(qū)間xn=(0.9).^n%計(jì)算相應(yīng)的x(n)k=0:19%設(shè)置DFT的采樣長度xw=dftmtx(20)xn%DFT變換figure(1)%畫出原序列的實(shí)部、虛部、模、相角subplot(221)stem(nreal(xn)filled)xlabel(n)ylabel(r

7、eal(xn))title(xn實(shí)部)subplot(222)stem(nimag(xn)filled)xlabel(n)ylabel(imag(xn))title(xn虛部)subplot(223)stem(nabs(xn)filled)xlabel(n)ylabel(abs(xn))title(xn模)subplot(224)stem(nangle(xn)filled)xlabel(n)ylabel(angle(xn))title(

8、xn相角)figure(2)%畫出頻譜的幅度譜、頻譜實(shí)部、頻譜虛部subplot(311)stem(kabs(xw)filled)xlabel(k)ylabel(abs(xw))title(幅度譜)subplot(312)stem(kreal(xw)filled)xlabel(k)ylabel(real(xw))title(頻譜實(shí)部)subplot(313)stem(kimag(xw)filled)xlabel(k)ylabel(ima

9、g(xw))title(頻譜虛部)2.復(fù)指數(shù)序列a=0.5b=0.8length=10????????otherwiselengthnjbanxn010)()(clcclfclear%清除n=0:9%設(shè)置自變量區(qū)間xn=(0.51j0.8).^n%計(jì)算相應(yīng)的x(n)k=0:9%設(shè)置DFT的采樣長度xw=dftmtx(10)xn%DFT變換figure(1)%畫出原序列的實(shí)部、虛部、模、相角subplot(221)stem(nreal(x