實(shí)驗(yàn)五-使用matlab實(shí)現(xiàn)卷積的運(yùn)算

實(shí)驗(yàn)五-使用matlab實(shí)現(xiàn)卷積的運(yùn)算實(shí)驗(yàn)五-使用matlab實(shí)現(xiàn)卷積的運(yùn)算實(shí)驗(yàn)五-使用matlab實(shí)現(xiàn)卷積的運(yùn)算實(shí)驗(yàn)五-使用matlab實(shí)現(xiàn)卷積的運(yùn)算編制僅供參考審核批準(zhǔn)生效日期地址：電話：傳真:郵編:實(shí)驗(yàn)五使用matlab實(shí)現(xiàn)卷積的運(yùn)算一實(shí)驗(yàn)?zāi)康膶W(xué)習(xí)MATLAB語言的編程方法及熟悉MATLAB指令；深刻理解卷積運(yùn)算，利用離散卷積實(shí)現(xiàn)連續(xù)卷積運(yùn)算；二實(shí)驗(yàn)內(nèi)容完成與兩函數(shù)的卷積運(yùn)算其中：在一個(gè)圖形窗口中，畫出、以及卷積結(jié)果。要求每個(gè)坐標(biāo)系有標(biāo)題、坐標(biāo)軸名稱。>>p=;t=0:p:10;f1=exp(-2*t).*u(t);f2=u(t)-u(t-4);f=conv(f1,f2);subplot(1,3,1);plot(t,f1,'r');title('f1(t)=e^-2*t*u(t)');xlabel('t(sec)');ylabel('f1(t)');subplot(1,3,2);plot(t,f2,'g');title('f2(t)=u(t)-u(t-4)');xlabel('t(sec)');ylabel('f2(t)');subplot(1,3,3);plot(f);title('f(t)=f1(t)*f2(t)');xlabel('t(sec)');ylabel('f(t)');若系統(tǒng)模型為：其中求零狀態(tài)響應(yīng)，畫出波形（函數(shù)本身畫出一幅圖，自己再畫出一幅輸入波形圖）。a=[144];b=[13];sys=tf(b,a);td=;t=0:td:10;f=exp(-t).*u(t);y=lsim(sys,f,t);plot(t,y);xlabel('t(sec)');ylabel('y(t)'); a=[144];b=[13];sys=tf(b,a);td=;t=0:td:10;f=exp(-t).*u(t);plot(t,f);xlabel('t(sec)');ylabel('f(t)');001234567891001t(sec)f(t)三實(shí)驗(yàn)原理：離散卷積和：調(diào)用函數(shù)：conv（）為離散卷積和，其中，f1(k),f2(k)為離散序列，K=…-2,-1,0,1,2,…。但是，conv函數(shù)只給出縱軸的序列值的大小，而不能給出卷積的X軸序號。為得到該值，進(jìn)行以下分析：對任意輸入：設(shè)非零區(qū)間n1~n2，長度L1=n2-n1+1；非零區(qū)間m1~m2，長度L2=m2-m1+1。則：非零區(qū)間從n1+m1開始，長度為L=L1+L2-1，所以S（K）的非零區(qū)間為：n1+m1~n1+m1+L-1。連續(xù)卷積和離散卷積的關(guān)系：計(jì)算機(jī)本身不能直接處理連續(xù)信號，只能由離散信號進(jìn)行近似：設(shè)一系統(tǒng)（LTI）輸入為，輸出為，如圖所示。LTILTIt若輸入為f(t):得輸出：當(dāng)時(shí)：所以：如果只求離散點(diǎn)上的f值所以，可以用離散卷積和CONV（）求連續(xù)卷積，只需