資源描述：
《用MATLAB實(shí)現(xiàn)信號(hào)的時(shí)域卷積》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、.word可編輯.用MATLAB實(shí)現(xiàn)信號(hào)的時(shí)域卷積摘要：利用抽樣方法繪制連續(xù)信號(hào)和離散信號(hào)。編制函數(shù)計(jì)算連續(xù)時(shí)間信號(hào)，離散時(shí)間信號(hào)的卷積。關(guān)鍵字：MATLAB；離散卷積；連續(xù)卷積實(shí)驗(yàn)原理：卷積運(yùn)算在信號(hào)與系統(tǒng)分析中具有非常重要的意義，是信號(hào)與系統(tǒng)分析的基本方法之一。線性是不變系統(tǒng)的零狀態(tài)響應(yīng)就可以通過求解駛?cè)爰?lì)信號(hào)與系統(tǒng)單位沖擊響應(yīng)的卷積運(yùn)算獲得。在連續(xù)時(shí)間系統(tǒng)中，卷積運(yùn)算成為卷積積分，表示為在離散時(shí)間系統(tǒng)中，卷積運(yùn)算稱為卷積和，表示為：實(shí)驗(yàn)任務(wù)及實(shí)現(xiàn)方法：⑴編制計(jì)算離散時(shí)間序列卷積核的M函數(shù)
2、,該函數(shù)課計(jì)算兩離散序列和的卷積和。程序計(jì)算出卷積和的同時(shí)，還繪制出序列、和的時(shí)域波形圖，并返回的非零樣值點(diǎn)的對(duì)應(yīng)向量。任務(wù)（1）程序：functiondconv(x1,x2,y1,y2)%任意兩序列卷積x1=input('請(qǐng)輸入x1=');x2=input('請(qǐng)輸入x2=');y1=input('請(qǐng)輸入y1=');y2=input('請(qǐng)輸入y2=');subplot(131),stem(x1,y1),gridon,title('y1'),axis([(min(x1)-1),(max(x1)+1
3、),(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
4、+xmin1)+1;.專業(yè).專注..word可編輯.xx=0:(t-1);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)])⑵調(diào)用來完成下面兩離散序列的卷積和運(yùn)算，并繪制圖形：任務(wù)（2）程序：functionf=ddconv(m1,m2)%-----------------------------%離散序列1t
5、1=m1-5;t2=m1+5;t=t1:t2;%橫坐標(biāo)n=length(t);tt=t1:m1;%起始到跳變n1=length(tt);f1=[zeros(1,n1-1),ones(1,n-n1+1)];%subplot(141),stem(t,f1),gridon,title('階躍函數(shù)序列u(n-m1)'),axis([t1t2-0.21.5])%-------------------------------%離散序列2.專業(yè).專注..word可編輯.p1=m1-5;p2=m1+5;p=p1
6、:p2;%橫坐標(biāo)z=length(p);pp=p1:m2;%起始到跳變z1=length(pp);f2=[zeros(1,z1-1),ones(1,z-z1+1)];%subplot(142),stem(p,f2),gridon,title('離散階躍序列u(n-m2)'),axis([p1p2-0.21.5])%------------------------------------%兩離散序列相減k=min(min(t),min(p)):max(max(t),max(p));%共同橫坐標(biāo)s1
7、=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),(max(f)+0.5)])su
8、bplot(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:
