實驗2離散序列的基本運(yùn)算

1、實驗卡和實驗報告 信息科學(xué)與工程學(xué)院 課程編號 實驗項目序號本科學(xué)生實驗卡和實驗報告信息科學(xué)與工程學(xué)院通信工程專業(yè) 2013級 1301 班課程名稱：數(shù)字信號處理實驗項目：離散序列的基本運(yùn)算20152016學(xué)年 第 二 學(xué)期學(xué)號： 201308030104_ 姓名：_王少丹_ 專業(yè)年級班級： _通信1301_四合院_ 實驗室 組別_ 實驗日期 _2016 年_ 4_ 月_28 日課程名稱數(shù)字信號處理實驗課時4實驗項目名稱和編號離散序列的基本運(yùn)算同組者 姓 名實驗?zāi)康膶W(xué)會運(yùn)用MATLAB進(jìn)行離散序列的運(yùn)算，并掌握程序編寫的方法實驗環(huán)境MATLAB實驗內(nèi)容和原理實驗步驟方 法關(guān)鍵代碼1、 序列的移

2、位X1(n)=x(n-m)例：x0=u(n)（顯示-10<n<10）X1=u(n+6) （顯示-10<n<10）X2=u(n-4) （顯示-10<n<10）例：已知x(n)=2sin(2*pi*n/10)，求x(n-2)和x(n+2)在-210區(qū)間的波形2、 序列的相加情況1 兩序列有相同的維數(shù)例：求x(n)=(n-2)+ (n-4) (0<n<10) 情況2 兩序列具有不同的維數(shù) 例：x1(n)=u(n+2) (-4<n<6) X2(n)=u(n-4) (-5<n<8) 求x(n)=x1(n)+x2(n)=u(n-4)3

3、、 序列相乘例：4、 序列反折例：x(n)=e(-0.3n) (-4<n<4)求x(-n)5、 序列的尺度變換例：x(n)=sin(2*pi*n),求x(2n)和x(n/2)測試記錄分 析結(jié) 論  實驗數(shù)據(jù)與理論相符。（后附代碼）小 結(jié)對于離散序列的基本運(yùn)算，移位、相加、相乘、反折還有尺度變換，都有了更深的理解，相加相乘等運(yùn)算是兩個向量之間的運(yùn)算，因此參加運(yùn)算的兩個序列向量必須具有相同的維數(shù)，否則應(yīng)進(jìn)行相應(yīng)的處理。尤其是如何用matlab實現(xiàn)有了一些基本的了解。以下由實驗教師填寫記 事評 議  成績評定  平時成績_ 實驗報告成績_ 綜合成績 _ 指導(dǎo)

4、教師簽名：代碼：%序列的移位n1=-10;n2=10;k0=0;k1=-6;k2=4;n=n1:n2;%生成離散信號的時間序列x0=n>=k0;x1=(n-k1)>=0;x2=(n-k2)>=0;figuresubplot(3,1,1),stem(n,x0,'filled','k');axis(n1,n2,1.1*min(x0),1.1*max(x0);ylabel('u(n)');subplot(3,1,2),stem(n,x1,'filled','k')axis(n1,n2,1.1*min(x

5、1),1.1*max(x1);ylabel('u(n+6)');subplot(3,1,3),stem(n,x2,'filled','k');axis(n1,n2,1.1*min(x2),1.1*max(x2);ylabel('u(n-4)'); %x(n)=2sin(2*pi*n)/10,Çóx(n-2)x(n+2) -210n=-2:10;n0=2;n1=-2;x=2*sin(2*pi*n/10);x1=2*sin(2*pi*(n-n0)/10);x2=2*sin(2*pi*(n-n1)/10);figur

6、esubplot(3,1,1),stem(n,x,'filled','k');ylabel('x(n)');subplot(3,1,2),stem(n,x1,'filled','k');ylabel('x(n-2)');subplot(3,1,3),stem(n,x2,'filled','k');ylabel('x(n+2)'); %序列的相加（相同維度）clear all;n1=0;n2=10;n01=2;n01=4;n=n1:n2;x1=(n-n0

7、1)=0;x2=(n-n01)=0;x3=x1+x2;figuresubplot(3,1,1);stem(n,x1,'filled','k');axis(n1,n2,1.1*min(x1),1.1*max(x1);ylabel('¦Ä(n-2)');subplot(3,1,2);stem(n,x2,'filled','k');axis(n1,n2,1.1*min(x2),1.1*max(x2);ylabel('¦Ä(n-4)');subplot(3,1,3);

8、stem(n,x3,'filled','k');axis(n1,n2,1.1*min(x3),1.1*max(x3);ylabel('¦Ä(n-2)+¦Ä(n-4)'); %序列相乘clear all;n1=-4:10;n2=-2:6;n01=-1;n=min(n1,n2):max(n1,n2);x1=3*exp(-0.25*n);x2=(n2-n01)>=0;N=length(n);y1=zeros(1,N);y2=zeros(1,N);y1(find(n>=min(n1)&(n<

9、;=max(n1)=x1;y2(find(n>=min(n2)&(n<=max(n2)=x2;x=y1.*y2;figuresubplot(3,1,1),stem(n1,x1,'filled','k');ylabel('x1(n)');axis(min(n),max(n),1.1*min(x1),1.1*max(x1);subplot(3,1,2),stem(n2,x2,'filled','k');ylabel('x2(n)');axis(min(n),max(n),1.1*mi

10、n(x2),1.1*max(x2);subplot(3,1,3),stem(n,x,'filled','k');ylabel('x(n)');axis(min(n),max(n),1.1*min(x),1.1*max(x) %序列相加（不同維度）clear all;n1=-4:10;n01=-2;n2=-5:8;n02=4;n=min(n1,n2):max(n1,n2);x1=(n1-n01)>=0;x2=(n2-n02)>=0;N=length(n);y1=zeros(1,N);y2=zeros(1,N);y1(find(n>

11、=min(n1)&(n<=max(n1)=x1;y2(find(n>=min(n2)&(n<=max(n2)=x2;x=y1+y2;figuresubplot(3,1,1),stem(n1,x1,'filled','k');ylabel('x1(n)');axis(min(n),max(n),1.1*min(x),1.1*max(x);subplot(3,1,2),stem(n2,x2,'filled','k');ylabel('x2(n)');axis(min(n

12、),max(n),1.1*min(x),1.1*max(x);subplot(3,1,3),stem(n,x,'filled','k');ylabel('x(n)');axis(min(n),max(n),1.1*min(x),1.1*max(x); %序列反折clear all;n=-4:4;x=exp(-0.3*n);x1=fliplr(x);n1=-fliplr(n);figuresubplot(1,2,1),stem(n,x,'filled','k');title('x(n)');subplot(1,2,2),stem(n1,x1,'filled','k');title('x(-n)'); %序列的尺度變換clear all;n=0:40;tn=n./20;x=sin(2*pi*tn);x1=sin(2*pi*tn*2);x2=sin(2*pi*tn/2);figuresubplot(3,1,1),stem(tn,x,'filled','k');ylabel('x(n)');axis(0,2,1.1*min(x),1.1*max(x);subplot(3,1,2),stem