線性卷積與圓周卷積演示程序的設(shè)計.doc

實驗一 線性卷積與圓周卷積演示程序的設(shè)計 實驗報告姓名 學(xué)號 專業(yè)班級 指導(dǎo)老師 分?jǐn)?shù) 數(shù)字信號處理課程設(shè)計任務(wù)書題目1線性卷積演示程序的設(shè)計（線性移不變離散時間系統(tǒng)的求解）主要內(nèi)容1、動態(tài)演示線性卷積和圓周卷積的完整過程；2、對比分析線性卷積與圓周卷積的結(jié)果。設(shè)計要求1、動態(tài)演示線性卷積和圓周卷積的過程（即翻轉(zhuǎn)、移位、乘積、求和的過程）；2、圓周卷積默認(rèn)使用2序列中的最大長度，且卷積前可設(shè)定用以進(jìn)行混疊分析；3、根據(jù)實驗結(jié)果分析2類卷積的關(guān)系；4、利用FFT實現(xiàn)快速卷積，驗證時域卷積定理，并與直接卷積進(jìn)行效率對比。主要儀器設(shè)備1、計算機(jī)1臺，安裝MATLAB軟件主要參考文獻(xiàn)美維納.K.恩格爾，約翰.G.普羅科斯著，劉樹棠譯.數(shù)字信號處理使用MATLABM.西安：西安交通大學(xué)出版社，2002.飛思科技產(chǎn)品研發(fā)中心編著.MATLAB7輔助信號處理技術(shù)與應(yīng)用M.北京：電子工業(yè)出版社，2005.課程設(shè)計進(jìn)度安排（起止時間、工作內(nèi)容）課程設(shè)計共設(shè)16個設(shè)計題目，每班3至4人為1組，1人1套設(shè)備，每組選作不同的題目，4個班共分4批。完整課程設(shè)計共20學(xué)時，為期1周，具體進(jìn)度如下：5學(xué)時 學(xué)習(xí)題目相關(guān)知識，掌握實現(xiàn)原理；5學(xué)時 用MATLAB語言實現(xiàn)題目要求；5學(xué)時 進(jìn)一步完善功能，現(xiàn)場檢查、答辯；5學(xué)時 完成并提交課程設(shè)計報告。課程設(shè)計開始日期2013.12.30課程設(shè)計完成日期2014.1.5課程設(shè)計實驗室名稱健翔橋校區(qū)計算中心地 點計算中心資料下載地址各班公共郵箱實驗一 線性卷積與圓周卷積演示程序的設(shè)計一、 實驗?zāi)康哪康模?熟練掌握MATLAB工具軟件在工程設(shè)計中的使用； 熟練掌握線性卷積與圓周卷積的關(guān)系及LSI離散時間系統(tǒng)系統(tǒng)響應(yīng)的求解方法。要求： 動態(tài)演示線性卷積的完整過程； 動態(tài)演示圓周卷積的完整過程； 對比分析線性卷積與圓周卷積的結(jié)果。步驟： 可輸入任意2待卷積序列x1(n)、x2(n)，長度不做限定。測試數(shù)據(jù)為：x1(n)=1,1,1,1,0,0,1,1,1,1,0,0，x2(n)=0,1,2,1,0,0,0,1,2,1,0,0； 分別動態(tài)演示兩序列進(jìn)行線性卷積x1(n)x2(n)和圓周卷積x1(n)x2 (n)的過程；要求分別動態(tài)演示翻轉(zhuǎn)、移位、乘積、求和的過程； 圓周卷積默認(rèn)使用2序列中的最大長度，但卷積前可以指定卷積長度N 用以進(jìn)行混疊分析； 根據(jù)實驗結(jié)果分析兩類卷積的關(guān)系。 假定時域序列x1(n)、x2(n)的長度不小于10000，序列內(nèi)容自定義。利用 FFT實現(xiàn)快速卷積，驗證時域卷積定理，并與直接卷積進(jìn)行效率對比。二、實驗原理1、線性卷積：線性時不變系統(tǒng)（Linear Time-Invariant System, or L. T. I系統(tǒng)）輸入、輸出間的關(guān)系為：當(dāng)系統(tǒng)輸入序列為，系統(tǒng)的單位脈沖響應(yīng)為，輸出序列為，則系統(tǒng)輸出為：或 上式稱為離散卷積或線性卷積。圖1.1示出線性時不變系統(tǒng)的輸入、輸出關(guān)系。0L. T. Ih(n) L. T. I 圖1.1 線性時不變系統(tǒng)的輸入、輸出關(guān)系2、圓周卷積D F T設(shè)兩個有限長序列和，均為點長D F T 如果則 N上式稱為圓周卷積。注：為序列的周期化序列；為的主值序列。上機(jī)編程計算時，可表示如下：3、兩個有限長序列的線性卷積序列為點長，序列為點長，為這兩個序列的線性卷積，則為且線性卷積的最大長，也就是說當(dāng)和時。4、圓周卷積與線性卷積的關(guān)系序列為點長，序列為點長，若序列和進(jìn)行N點的圓周卷積，其結(jié)果是否等于該兩序列的線性卷積，完全取決于圓周卷積的長度：當(dāng)時圓周卷積等于線性卷積，即N當(dāng)時，圓周卷積等于兩個序列的線性卷積加上相當(dāng)于下式的時間混疊，即三、實驗步驟已知兩個有限長序列1、實驗前，預(yù)先筆算好這兩個序列的線性卷積及下列幾種情況的圓周卷積 2、編制一個計算圓周卷積的通用程序，計算上述4種情況下兩個序列與的圓周卷積。3、上機(jī)調(diào)試并打印或記錄實驗結(jié)果。4、將實驗結(jié)果與預(yù)先筆算的結(jié)果比較，驗證其正確性。五、實驗報告1、列出計算兩種卷積的公式，列出實驗程序清單（包括必要的程序說明）。2、記錄調(diào)試運行情況及所遇問題的解決方法。3、給出實驗結(jié)果，并對結(jié)果作出分析。驗證圓周卷積兩者之間的關(guān)系實驗結(jié)果（1） 程序clear all;N1=5; N2=4; xn=1,1,1,1,0,0,1,1,1,1,0,0;%生成x(n) hn=0,1,2,1,0,0,0,1,2,1,0,0;%生成h(n) yln=conv(xn,hn);%直接用函數(shù)conv計算線性卷積 ycn=circonv(xn,hn,5);%用函數(shù)circonv計算N1點圓周卷積 ny1=0:1:length(yln)-1; ny2=0:1:length(ycn)-1; subplot(2,1,1);%畫圖stem(ny1,yln);ylabel(線性卷積); subplot(2,1,2);stem(ny2,ycn); ylabel(圓周卷積);題目：已知兩個有限長序列x(n)=(n)+2(n-1)+3(n-2)+4(n-3)+5(n-4)h(n)=(n)+2(n-1)+(n-2)+2(n-3)計算以下兩個序列的線性卷積和圓周卷積（1）x(n)y(n) (2)x(n)y(n) (3)x(n)y(n) (4)x(n)y(n) 調(diào)用函數(shù)circonvfunction yc=circonv(x1,x2,N)%用直接法實現(xiàn)圓周卷積%y=circonv(x1,x2,N)%y:輸出序列%x1,x2:輸入序列%N:圓周卷積的長度if length(x1)N error;endif length(x2)N error;end%以上語句判斷兩個序列的長度是否小于Nx1=x1,zeros(1,N-length(x1);%填充序列x1(n)使其長度為N，序列h(n)的長度為N1,序列x(n)的長度為N2x2=x2,zeros(1,N-length(x2);%填充序列x2(n)使其長度為Nn=0:1:N-1;x2=x2(mod(-n,N)+1);%生成序列x2(-n)N，鏡像，可實現(xiàn)對x(n)以N為周期的周期延拓，加1是因為MATLAB向量下標(biāo)只能從1開始。H=zeros(N,N);%生成N行N列的零矩陣for n=1:1:N H(n,:)=cirshifted(x2,n-1,N);%該矩陣的k行為x2(k-1-n)Nendyc=x1*H;%計算圓周卷積 調(diào)用函數(shù)cirshiftdfunction y=cirshiftd(x,m,N)%直接實現(xiàn)序列x的圓周移位%y=cirshiftd(x,m,N)%x:輸入序列，且它的長度小于N%m:移位位數(shù)%N：圓周卷積的長度%y:輸出的移位序列if length(x)N error(x的長度必須小于N);endx=x,zeros(1,N-length(x);n=0:1:N-1;y=x(mod(n-m,N)+1); 函數(shù)（1）x(n)y(n)clear all;N1=5;N2=4;xn=1 2 3 4 5;%生成x(n)hn=1 2 1 2;%生成h(n)yln=conv(xn,hn);%直接用函數(shù)conv計算線性卷積ycn=circonv(xn,hn,5);%用函數(shù)circonv計算N1點圓周卷積ny1=0:1:length(yln)-1;ny2=0:1:length(ycn)-1;subplot(2,1,1);%畫圖stem(ny1,yln);ylabel(線性卷積);subplot(2,1,2);stem(ny2,ycn);ylabel(圓周卷積); 函數(shù)（2）x(n)y(n)clear all;N1=5;N2=4;xn=1 2 3 4 5;%生成x(n)hn=1 2 1 2;%生成h(n)yln=conv(xn,hn);%直接用函數(shù)conv計算線性卷積ycn=circonv(xn,hn,6);%用函數(shù)circonv計算N1點圓周卷積ny1=0:1:length(yln)-1;ny2=0:1:length(ycn)-1;subplot(2,1,1);stem(ny1,yln);ylabel(線性卷積);subplot(2,1,2);stem(ny2,ycn);ylabel(圓周卷積); 函數(shù)（3）x(n)y(n)clear all;N1=5;N2=4;xn=1 2 3 4 5;%生成x(n)hn=1 2 1 2;%生成h(n)yln=conv(xn,hn);%直接用函數(shù)conv計算線性卷積ycn=circonv(xn,hn,9);%用函數(shù)circonv計算N1點圓周卷積ny1=0:1:length(yln)-1;ny2=0:1:length(ycn)-1;subplot(2,1,1);stem(ny1,yln);ylabel(線性卷積);subplot(2,1,2);stem(ny2,ycn);ylabel(圓周卷積); 函數(shù)（4）x(n)y(n)clear all;N1=5;N2=4;xn=1 2 3 4 5;%生成x(n)hn=1 2 1 2;%生成h(n)yln=conv(xn,hn);%直接用函數(shù)conv計算線性卷積ycn=circonv(xn,hn,10);%用函數(shù)circonv計算N1點圓周卷積ny1=0:1:length(yln)-1;ny2=0:1:length(ycn)-1;subplot(2,1,1);stem(ny1,yln);ylabel(線性卷積);subplot(2,1,2);stem(ny2,ycn);ylabel(圓周卷積);六、思考題：圓周卷積與線性卷積的關(guān)系：若有x1(n)與x2（n）兩個分別為N1與N2的有限長序列，則它們的線性卷積y1（n）為N1+N2-1的有限長序列，而它們的N點圓周卷積y2（n）則有以下兩種情況：1，當(dāng)N N1+N2-1時，y2（n）的前N1+N2-1的點剛好是y1（n）的全部非零序列，而剩下的N-(N1+N2-1)個點上的序列則是補(bǔ)充的零。線性卷積運算步驟：求x1(n)與x2（n） 的線性卷積：對x1(m)或x2（m）先進(jìn)行鏡像移位x1（-m），對移位后的序列再進(jìn)行從左至右的依次平移x(n-m),當(dāng)n=0,1,2.N-1時，分別將x(n-m)與x2（m）相乘，并在m=0,1,2.N-1的區(qū)間求和，便得到y(tǒng)（n）圓周卷積運算步驟：圓周卷積過程中，求和變量為m，n為參變量，先將x2(m)周期化，形成x2(m)N,再反轉(zhuǎn)形成x2(-m)N,取主值序列則得到x2(-m)NRN(m),通常稱之為x2(m)的圓周反轉(zhuǎn)。對x2(m)圓周反轉(zhuǎn)序列圓周右移n,形成x2(n-m)NRN(m),當(dāng)n=0,1,2,N-1時，分別將x1(m)與x2(n-m)NRN(m)相乘，并在m=0到N-1區(qū)間內(nèi)求和，便得到圓周