哈工程數(shù)字信號處理實驗報告1

數(shù)字信號處理實驗實驗一：基本信號班級：姓名：學(xué)號：指導(dǎo)教師：2012年10月實驗一：基本信號一：實驗原理：本節(jié)專注于用MATLAB產(chǎn)生一些基本離散信號的問題。主要是有那個MATLAB內(nèi)部向量程序來產(chǎn)生信號。用MATLAB的stem指令會出離散時間信號。依據(jù)MATLAB的編址約定，標(biāo)號n=0必須對應(yīng)nn（1）；必須給指定向量的第一個參數(shù)以得到正確的n軸。二：實驗內(nèi)容：沖擊信號產(chǎn)生并繪出下面的序列。在每種情況下，水平n軸應(yīng)該只在指定的區(qū)間上展開并應(yīng)該相應(yīng)標(biāo)注。使用stem指令使每個序列顯示成離散時間信號。x[n]=0.9δ[n-5] 1<=n<=20x[n]=0.8δ[n] -15<=n<=15x[n]=1.5δ[n-333] 300<=n<=350x[n]=4.5δ[n+7] -10<=n<=0L=20;nn=1:(L);imp=zeros(L,1);imp(5)=0.9;stem(nn,imp))L=31;nn=-15:(L-16);imp=zeros(L,1);imp(16)=0.8;stem(nn,imp))L=51;nn=300:350;imp=[zeros(L,1)]';imp(34)=1.5stem(nn,imp)L=11;nn=-10:(L-11);imp=zeros(L,1);imp(4)=4.5;stem(nn,imp)實驗分析：所得4個圖形均符合題目要求3、指數(shù)信號衰減的指數(shù)信號是數(shù)字信號是數(shù)字信號處理的基本信號。因為它是線性常系數(shù)差分方程的解。使用函數(shù)在區(qū)間n=0,1,2，。。。，20上繪出信號x[n]=(0.9)?。在許多推導(dǎo)中，指數(shù)信號序列a?u[n]須在有限區(qū)間上求和。使用（a）中的函數(shù)產(chǎn)生一個指數(shù)信號然后對其求和并比較結(jié)果。指數(shù)序列在信號處理中常常出現(xiàn)的一個原因是，時移并不改變其信號特征。證明一有限長指數(shù)信號滿足移位關(guān)系：y[n]=ay[n-1], 1<=n<=L-1比較向量y（2：L）和a*y(1:L-1)。在MATLAB中移位有限長度信號時，因為不能自動補零，所以必須留意信號終點。D.產(chǎn)生指數(shù)信號另外的方法是使用差分方程給出的遞歸表示式。當(dāng)輸入x[n]是一個沖擊信號時，信號y[n]=a?u[n]是下面差分方程的解y[n]-ay[n-1]=x[n],初始條件y[-1]=0,由于假定差分方程以因果方式遞歸（即n增大），n=-1的初始條件是必需的。用MATLAB的filter函數(shù)可實現(xiàn)差分方程。首先創(chuàng)建并保存如下M文件：functiony=genexp(b,n0,L)if(L<=0)errorendnn=n0+[1:L]'-1;y=b.^nn;stem(nn,y)enda．y=genexp(0.9,0,21)nn=0:20stem(nn,y)b．y=genexp(0.9,0,21)x=sum(y)y=1.00000.90000.81000.72900.65610.59050.53140.47830.43050.38740.34870.31380.28240.25420.22880.20590.18530.16680.15010.13510.1216x=8.9058實驗分析：用題目所給公式所求結(jié)果為8.9058，結(jié)果一致。c．y1=[genexp(0.9,2,20)]'nn=2:21stem(nn,y1)y1=Columns1through160.81000.72900.65610.59050.53140.47830.43050.38740.34870.31380.28240.25420.22880.20590.18530.1668Columns17through200.15010.13510.12160.1094x2=[genexp(0.9,1,20)]';y2=0.9*x2nn=[1:20];stem(nn,y2)y2=Columns1through160.81000.72900.65610.59050.53140.47830.43050.38740.34870.31380.28240.25420.22880.20590.18530.1668Columns17through200.15010.13510.12160.1094實驗分析：由上圖比較可知，證明成立d．a(chǎn)=[1,-0.9];b=1;xin=zeros(21,1);xin(1)=1;nn=0:20;yout=[filter(b,a,xin)]'stem(nn,yout)yout=Columns1through121.00000.90000.81000.72900.65610.59050.53140.47830.43050.38740.348