1、.word可編輯.用MATLAB;現信號的時域卷積摘要：利用抽樣方法繪制連續信號和離散信號。編制函數計算連續時間信號，離散時間信號的卷 積。關鍵字：MATLAB離散卷積；連續卷積實驗原理：卷積運算在信號與系統分析中具有非常重要的意義，是信號與系統分析的基本方法之 一。線性是不變系統的零狀態響應就可以通過求解駛入激勵信號與系統單位沖擊響應的卷 積運算獲得。在連續時間系統中，卷積運算成為卷積積分，表示為f t = f1tf2 t =f1f2 t - d-nd在離散時間系統中，卷積運算稱為卷積和，表示為：x n ) u x1n x 2 n)； =x1 m x2 n - mm 二“：；實驗任務及實現方

2、法：編制計算離散時間序列卷積核的M函數dconv()該函數課計算兩離散序列x1(n )和X2 (n )的卷積和x(n )= xi(n3X2(n )。程序計算出卷積和x(n )的同時，還繪制出序列x1m y x2(n押x(nm時域波形圖，并返回x(n )的非零樣值點的對應向量。任務(1)程序：function dconv(x1,x2,y1,y2)%壬意兩序列卷積x1=input('請輸入 x1=');x2=input('請輸入 x2=');y1=input('請輸入 y1=');y2=input('請輸入 y2=');subplot

3、(131),stem(x1,y1),gridon,title('y1' ),axis(min(x1)-1),(max(x1)+1),(min(y1)-1),(max(y1)+1)subplot(132),stem(x2,y2),gridon,title('y2' ),axis(min(x2)-1),(max(x2)+1),(min(y2)-1),(max(y2)+1)%xmin1=min(x1);xmax1=max(x1);xmin2=min(x2);xmax2=max(x2);t=(xmax2+xmax1)-(xmin2+xmin1)+1;xx=0:(t-1)

4、;ft=conv(y1,y2)subplot(133),stem(xx,ft),gridon,title( '*')axis(min(xx)+0.5),(max(xx)+0.5),(min(ft)-0.5),(max(ft)+0.5)專業.專注 調用dconv(原完成下面兩離散序列的卷積和運算，并繪制圖形:hn =xn =un -u n -4任務(2)程序：function f=ddconv(m1,m2)%離散序列1t1=m1-5;t2=m1+5;t=t1:t2;%K 坐標n=length(t);tt=t1:m1;%起始到跳變n1=length(tt);階躍函數序列u(n-m1

5、) '),axis(t1 t2 -0.2f1=zeros(1,n1-1),ones(1,n-n1+1);%subplot(141),stem(t,f1),grid on,title('1.5)%離散序列2p1=m1-5;p2=m1+5;p=p1:p2;%K 坐標z=length(p);pp=p1:m2;炮始到跳變z1=length(pp);離散階躍序列u(n-m2) '),axis(p1 p2 -0.2f2=zeros(1,z1-1),ones(1,z-z1 + 1);%subplot(142),stem(p,f2),grid on,title('1.5)%兩離

6、散序列相減k=min(min(t),min(p):max(max(t),max(p);益同橫坐標s1=zeros(1,length(k);s2=s1;s1(find(k>=min(t)&(k<=max(t)=1)=f1;s2(find(k>=min(p)&(k<=max(p)=1)=f2;f=s1-s2;subplot(131),stem(k,f),gridon,title('x(n)=u(n)-u(n-4)')axis(min(min(t),min(p)-1),(max(max(t),max(p)+1),(min(f)-0.5),(ma

7、x(f)+0.5)subplot(132),stem(k,f),gridon,title('h(n)=u(n)-u(n-4)')axis(min(min(t),min(p)-1),(max(max(t),max(p)+1),(min(f)-0.5),(max(f)+0.5)%xmin1=min(k);xmax1=max(k);xmin2=min(k);xmax2=max(k);t=(xmax2+xmax1)-(xmin2+xmin1)+1;xx=0:(t-1);ft=conv(f,f)subplot(133),stem(xx,ft),gridon,title( 'h(n

8、)*x(n)')axis(min(xx)+0.5),(max(xx)+0.5),(min(ft)-0.5),(max(ft)+0.5)x(nu(n)-u(n-4)h(n)=u(n)-u(n-4)h(n)*x(n)1020 編制計算連續時間信號卷積積分的M函數sconv(，該函數在計算出兩信號f1(t )和f2(t)的卷積積分f(t的數值近似的同時，還繪制出M(t卜f2(t)和f(t)的時域波形圖。實現方法：與第(4)任務用同一個程序。(4)調用函數sconv(原完成下面兩連續時間信號的卷積積分運算，并繪制圖形：111e(t)=u t+- i-u(t -1 ), h(t )= tU(t

9、)-u(t-2)< 2)2任務(3、4)程序：function sconv()%u(t)m1=0;t1=m1-10;t2=m1+10;t=t1:0.01:t2;n=length(t);tt=t1:0.01:m1;n1=length(tt);f1=zeros(1,n1-1),ones(1,n-n1+1);%subplot(131),plot(t,f1),grid on,title('a ?e 0 ?x?D6 a D u(t) '),axis(t1 t2 -0.2 1.5)%u(t+1/2)tx1=t-1/2;fx1=f1;%subplot(241),stairs(tx1,f

10、x1),grid on,title(' u(t+1/2) '),axis(t1 t2 -0.2 1.5)%u(t-1)tx2=t+1;fx2=f1;%subplot(242),stairs(tx2,fx2),grid on,title(' u(t-1) '),axis(t1 t2 -0.2 1.5)k=min(min(tx1),min(tx2):0.01:max(max(tx1),max(tx2); s1=zeros(1,length(k);s2=s1;s1(find(k>=min(tx1)&(k<=max(tx1)=1)=fx1;s2(fi

11、nd(k>=min(tx2)&(k<=max(tx2)=1)=fx2;e=s1-s2;subplot(131),stairs(k,e),gridon,title( 'e(t)=u(t+1/2)-u(t-1)')axis(min(min(tx1),min(tx2)+5),(max(max(tx1),max(tx2)-5),(min(e)-0.5),(max(e)+0.5)tx3=t;fx3=f1;%subplot(243),stairs(tx3,fx3),grid on,title(' u(t) '),axis(t1 t2 -0.2 1.5)t

12、x4=t+2;fx4=f1;%subplot(244),stairs(tx4,fx4),grid on,title(' u(t-2) '),axis(t1 t2 -0.2 1.5)k1=min(min(tx3),min(tx4):0.01:max(max(tx3),max(tx4);s3=zeros(1,length(k1);s4=s1;s3(find(k1>=min(tx3)&(k1<=max(tx3)=1)=fx3;s4(find(k1>=min(tx4)&(k1<=max(tx4)=1)=fx4;h=s3-s4;%subplot(2

13、36),stairs(k1,h),grid on,title(' u(t)-u(t-2)')%axis(min(min(tx3),min(tx4)+5),(max(max(tx3),max(tx4)-5),(min(h)-0.5),(max(h)+0 .5)h(t)=1/2*t*u(t)-u(t-2)kk1=length(k1);ont=0:0.01:kk1;ks=min(min(k1),min(ont):0.01:max(max(k1),max(ont);s5=zeros(1,length(ks);s6=s5;s5(find(ks>=min(ont)&(ks&l

14、t;=max(ont)=1)=ont;s6(find(ks>=min(k1)&(ks<=max(k1)=1)=h;ffu=1/2*(s5.*s6);subplot(132),stairs(ks,ffu),gridon,title( '1/2*t*u(t)-u(t-2)')axis(-5,5,(min(ffu)-0.5),(max(ffu)+0.5)xmin1=min(k);xmax1=max(k);xmin2=min(ks);xmax2=max(ks);ttx=(xmax2+xmax1)-(xmin2+xmin1)+0.01;xx=0:0.01:(ttx-0.01);ft=conv(e,ffu);subplot(133),stairs(xx,ft),gridon,title( 'h(n)*x(n)&