快速卷積matlab實(shí)現(xiàn)(編輯修改稿)
2025-04-21 02:05 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 x)。 %調(diào)用conv函數(shù)直接計(jì)算線性卷積stem(y)。 %畫出卷積結(jié)果h(n)的序列圖title(39。y(n)39。) 運(yùn)行結(jié)果及分析 程序運(yùn)行后，得到的卷積結(jié)果y(n)。 編程得到的卷積結(jié)果y(n) ，編程得到的卷積結(jié)果序列y(n)與自己先前計(jì)算的結(jié)果相同，說明所設(shè)計(jì)的程序是正確的，得到了正確的結(jié)果。 編程得到的卷積結(jié)果y(n) 題目二用函數(shù)conv和FFT計(jì)算長(zhǎng)為1000序列的卷積，比較其計(jì)算時(shí)間。 設(shè)計(jì)內(nèi)容及原理分析 本題目的在于比較直接卷積和快速卷積的計(jì)算時(shí)間，分析其優(yōu)劣性。：（1）序列補(bǔ)零：將兩序列都補(bǔ)零到L點(diǎn)；（2）計(jì)算X1(k)=FFT[x1(n)]；（3）計(jì)算X2(k)=FFT[x2(n)]；（4）計(jì)算Y(k)=X1(k)X2(k)；（5）計(jì)算y(n)=x1(n)*x2(n)=IFFT[Y(k)]。本題中已經(jīng)給定兩原序列的長(zhǎng)度M=N=1000，則L≥M+N1，而因?yàn)榭焖倬矸e的基礎(chǔ)是FFT，所以要求L滿足L=2n，為減少運(yùn)算量選取L=2048。設(shè)兩原序列分別為： x1=(2*n) ； x2=n^3 計(jì)算時(shí)間的獲取方法為：計(jì)算前先調(diào)用clock函數(shù)讀取瞬時(shí)時(shí)鐘，待計(jì)算結(jié)束后，調(diào)用etime(t1，t2)函數(shù)計(jì)算時(shí)刻t1，t2間所經(jīng)歷的時(shí)間。 程序設(shè)計(jì)框圖 程序設(shè)計(jì)框圖 程序代碼N=1000。L=pow2(nextpow2(1000+10001))。 %計(jì)算L的值n=1:N。x1=*sin(2*n)。 x2=n.^3。t0=clock。 %調(diào)用clock函數(shù)讀取瞬時(shí)時(shí)鐘yc=conv(x1,x2)。 %用函數(shù)conv計(jì)算卷積conv_time=etime(clock,t0) %調(diào)用etime函數(shù)計(jì)算時(shí)間t0=clock。yf=ifft(fft(x1,L).*fft(x2,L))。 %用函數(shù)FFT計(jì)算卷積fft_time=etime(clock,t0)subplot(321),stem(x1,39。.39。)。ylabel(39。x1(n)39。)。subplot(322),stem(x2,39。.39。)。ylabel(39。x2(n)39。)。subplot
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1