數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告-離散時(shí)間序列卷積和MATLAB實(shí)現(xiàn)

數(shù)字信號(hào)處理實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)一：離散時(shí)間序列卷積和MATLAB實(shí)現(xiàn)一、實(shí)驗(yàn)?zāi)康耐ㄟ^本實(shí)驗(yàn)學(xué)會(huì)用MATLAB對(duì)信號(hào)與系統(tǒng)進(jìn)行分析，實(shí)現(xiàn)離散序列卷積和的計(jì)算。二、實(shí)驗(yàn)儀器：（1）微機(jī)（2）MATLAB編程工具三：實(shí)驗(yàn)內(nèi)容題一：已知序列f1(k)= f2(k)= 調(diào)用conv()函數(shù)求上述兩序列的卷積和題二：編寫計(jì)算兩離散序列卷積和f(k)=f1(k)*f2(k)的實(shí)用函數(shù)dconv().要求該程序在計(jì)算出卷積和f(k)的同時(shí)，還繪出序列f1(k),f2(k)和f(k)的時(shí)域波形圖。 function[f,k]=dconv(f1,f2,k1,k2) %f1(k),f2(k)及f(k)的對(duì)應(yīng)序號(hào)向量分別為k1，k2和k。題三：試用MATLAB計(jì)算如下所示序列f1(k)與f2(k)的卷積和f(k)，繪出它們的時(shí)域波形，并說明序列f1(k)與f2(k)的時(shí)域?qū)挾扰c序列f(k)的時(shí)域?qū)挾鹊年P(guān)系。 f1(k)= f2(k)=提示：可用上述dconv()的函數(shù)來解決。題四：已知某LTI離散系統(tǒng)，其單位響應(yīng)h(k)=e(k)-e(k-4)，求該系統(tǒng)在激勵(lì)為f(k)=e(k)-e(k-3)時(shí)的零狀態(tài)響應(yīng)，并繪出其時(shí)域波形圖。 四、實(shí)驗(yàn)原理、方法和手段1、離散時(shí)間序列f1(k)和f2(k)的卷積和定義： f(k)=f1(k)*f2(k)= 2、在離散信號(hào)與系統(tǒng)分析中有兩個(gè)與卷積和相關(guān)的重要結(jié)論： a、f(k)==f(k)*δ(k)即離散序列可分解為一系列幅度由f(k)決定的單位序列δ(k)及其平移序列之積。 b、對(duì)線性時(shí)不變系統(tǒng)，設(shè)其輸入序列為f(k)，單位響應(yīng)為h(k)，其零狀態(tài)響應(yīng)為y(k)，則有：y(k)=五、matlab程序?qū)崿F(xiàn)卷積和計(jì)算：1、題一(conv.m)clc;k1=3;k2=3;k=k1+k2-1;f1=[1,1,1];f2=[0,1,2,3];f=conv(f1,f2);nf=0:k;stem(nf,f,'*r');xlabel('n');ylabel('f(n)');gridon;2、題二(dconv.m)function[f,k]=dconv(f1,f2,k1,k2);f=conv(f1,f2);k(1)=k1(1)+k2(1);k(length(f))=k1(length(f1))+k2(length(f2));k=k(1):k(length(f));subplot(1,3,1);stem(k1,f1,'*r');xlabel('k1');ylabel('f1(k1)');gridon;subplot(1,3,2);stem(k2,f2,'*r');xlabel('k2');ylabel('f2(k2)');gridon;subplot(1,3,3);stem(k,f,'*r');xlabel('k');ylabel('f(k)');gridon;3、題三(san.m)clc;f1=[1,2,1];f2=[1,1,1,1,1];k1=-1:1;k2=-2:2;dconv(f1,f2,k1,k2);4、題四(si.m)clc;f1=[1,1,1,1];f2=[1,1,1];k1=0:3;k2=0:2;dconv(f1,f2,k1,k2);六、實(shí)驗(yàn)心得：(1)通過這次實(shí)驗(yàn)，讓我知道序列的卷積和可以用matlab編程實(shí)現(xiàn)，并且實(shí)現(xiàn)的方法比較簡單，而且matlab本身有豐富的函數(shù)庫，畫圖、求卷積和時(shí)只要調(diào)用相應(yīng)的函數(shù)即可。(2)這次實(shí)驗(yàn)的難點(diǎn)主要就在寫dconv()函數(shù)和各序列長度的計(jì)算，如果不注意其中的計(jì)算關(guān)系則很容易出錯(cuò)，而dconv()只要寫錯(cuò)，后面的題目的圖形就不能準(zhǔn)確的繪制出來。對(duì)于后面的3、4兩個(gè)小題只要把相應(yīng)的序列寫出來，寫入用matlab表示的矩陣中，在調(diào)用dconv()畫出圖形即可。(3)由于matlab這門課的學(xué)習(xí)距現(xiàn)在隔