有限長(zhǎng)序列、頻譜、DTFT的性質(zhì)

1、實(shí)驗(yàn)報(bào)告課程名稱： 數(shù)字信號(hào)處理 指導(dǎo)老師： 劉英 成績(jī)：_實(shí)驗(yàn)名稱： 有限長(zhǎng)序列、頻譜、DTFT的性質(zhì) 實(shí)驗(yàn)類型： 演示 同組學(xué)生姓名： 一、實(shí)驗(yàn)?zāi)康暮鸵笸ㄟ^演示實(shí)驗(yàn)，建立對(duì)典型信號(hào)及其頻譜的直觀認(rèn)識(shí)，理解DTFT的物理意義、主要性質(zhì)。二、實(shí)驗(yàn)內(nèi)容和步驟2-1用MATLAB，計(jì)算得到五種共9個(gè)序列： 2-1-1單位階躍序列 x(n)=2u(n-5)2-1-2實(shí)指數(shù)序列例如，a=0.5, length=10 a=0.9, length=10 a=0.9, length=202-1-3復(fù)指數(shù)序列例如，a=0.4, b=0.6, length=202-1-4從正弦信號(hào)x(t)=sin(2pft+

2、delta)抽樣得到的正弦序列x(n)=sin(2pfnT+delta)。如，信號(hào)頻率f=1Hz，初始相位delta=0，抽樣間隔T=0.1秒，序列長(zhǎng)length=10。2-1-5含兩個(gè)頻率分量的復(fù)合函數(shù)序列x(n)=sin(2pf1nT)+delta×sin(2pf2nT+phi)。如，頻率f1(Hz)頻率f2(Hz)相對(duì)振幅delta初相位phi (度)抽樣間隔T(秒)序列長(zhǎng)length130.500.110130.5900.110130.51800.1102-2 用MATLAB，對(duì)上述各個(gè)序列，重復(fù)下列過程。2-2-1畫出一個(gè)序列的實(shí)部、虛部、模、相角；觀察并記錄實(shí)部、虛部、模

3、、相角的特征。2-2-2 計(jì)算該序列的幅度譜、頻譜實(shí)部、頻譜虛部；觀察和并記錄它們的特征，給予解釋。2-2-3 觀察同種序列取不同參數(shù)時(shí)的頻譜，發(fā)現(xiàn)它們的差異，給予解釋。實(shí)驗(yàn)名稱：_姓名：_學(xué)號(hào)：_ 2.3 用MATLAB，對(duì)下面系統(tǒng)進(jìn)行分析因果系統(tǒng), 求（1） 求出, 畫出零極點(diǎn)示意圖（2） 畫出和三、主要儀器設(shè)備MATLAB編程。四、操作方法和實(shí)驗(yàn)步驟（參見“二、實(shí)驗(yàn)內(nèi)容和步驟”）五、實(shí)驗(yàn)數(shù)據(jù)記錄和處理列出MATLAB程序清單，加注釋。2-12-1-1N=15;x=ones(1,N);%生成行矩陣xn=5:N-1+5;%橫坐標(biāo)figure(1);xlabel('x(n)')

4、;ylabel('n');stem(xn,2*x);%縱坐標(biāo)變?yōu)?倍axis on;axis(-1 14 0 4);set(gca,'YTick',0:1:4);set(gca,'XTick',-1:1:14);title('單位階躍序列');2-1-2N=10;%可修改參數(shù)a=0.5;%可修改參數(shù)xn=0:N-1;%橫坐標(biāo)figure(1);xlabel('x(n)');ylabel('n');stem(xn,a.(xn);%畫指數(shù)函數(shù) axis on;axis(-3 N+3 0 a0);set(

5、gca,'XTick',-1:1:N+2);title('a=0.5 length=10');2-1-3N=20;%可修改參數(shù)a=0.4;%可修改參數(shù)b=0.6;%可修改參數(shù)xn=0:N-1;%橫坐標(biāo)figure(1);xlabel('x(n)');ylabel('n');stem(xn,(a+j*b).(xn);%畫復(fù)指數(shù)函數(shù) axis on;axis(-3 N+3 -0.5 (a+j*b)0);set(gca,'XTick',-1:1:N+2);title('a=0.4 b=0.6 length=20&

6、#39;);2-1-4f=1;%可修改參數(shù)delta=0;%可修改參數(shù)T=0.1;%可修改參數(shù)N=10;%可修改參數(shù)xn=0:T:T*(N-1);%橫坐標(biāo)figure(1);stem(xn,sin(2*pi*f*xn+delta);%畫采樣后的結(jié)果hold on;%停留t=-0.5:0.001:1.5;%畫原始函數(shù)plot(t,sin(2*pi*f*t);axis on;axis(-0.5 1.5 -1.2 1.2);grid on;title('Sample');2-1-5f1=1;%可修改參數(shù)f2=3;%可修改參數(shù)delta=0.5;%可修改參數(shù)phi=0;%可修改參數(shù)T=

7、0.1;%可修改參數(shù)N=10;%可修改參數(shù)xn=0:T:T*(N-1);%橫坐標(biāo)figure(1);stem(xn,sin(2*pi*f1*xn)+delta*sin(2*pi*f2*xn+phi);%畫采樣后的結(jié)果hold on;%停留t=-0.5:0.001:1.5;%畫原始函數(shù)plot(t,sin(2*pi*f1*t)+delta*sin(2*pi*f2*t+phi);axis on;axis(-0.5 1.5 -1.2 1.2);grid on;title('2-1-5(1)');2-2-1(5)f1=1;%可修改參數(shù)f2=3;%可修改參數(shù)delta=0.5;%可修改參

8、數(shù)phi=180;%可修改參數(shù)T=0.1;%可修改參數(shù)N=10;%可修改參數(shù)xn=0:T:T*(N-1);%橫坐標(biāo)figure(1);yn=sin(2*pi*f1*xn)+delta*sin(2*pi*f2*xn+phi);yn0=real(yn);%實(shí)部yn1=imag(yn);%虛部yn2=abs(yn);%模yn3=angle(yn);%角subplot(2,2,1);stem(xn,yn0,'filled');title('real');subplot(2,2,2);stem(xn,yn1,'filled');title('imag');subplot(2,2,3);stem(xn,yn2,'filled');title('abs');subplot(2,2,4);stem(xn,yn3,'filled');title('angle');六、實(shí)驗(yàn)結(jié)果與分析觀察實(shí)驗(yàn)結(jié)果（數(shù)據(jù)及圖形）的特征，做必要的記錄，做出解釋。包括：6-1 各種序列的圖形（時(shí)