dsp第十一講ppt課件-資料下載頁
2025-05-05 12:10本頁面
【正文】 [1， 2， 3， 4， 4， 3， 2， 1] [1， 2， 3， 4， 4， 3， 2， 1] [3， 3， 3， 4， 4， 3 直接做圓周卷積 利用圓周卷積與線性卷積的關系 圖 用 DFT計算線性卷積框圖 補L － N個零點L點D F T補L － M個零點L點D F TL 點I D F Ty ( n )h ( n )x ( n ) MN。 若仍選取 L≥N＋ M－ 1， 以 L為循環(huán)卷積區(qū)間 ， 并用上述快速卷積法計算線性卷積 ， 則要求對短序列補很多零點 ， 而且長序列必須全部輸入后才能進行快速計算 。 因此要求存儲容量大 ， 運算時間長 ， 并使處理延時很大 ， 不能實現(xiàn)實時處理 。 通常采用分段卷積 。 若 將 x(n)均勻分段 ， 每段長度取 M， 則 于是 ， h(n)與 x(n)的線性卷積可表示為 圖 重疊相加法卷積示意圖 M0NM Mx1( n )x0( n ) x2( n )N ＋ M － 1N ＋ M － 1y0( n )y1( n )N ＋ M － 1y2( n )2 M M3 M ＋ N － 10N － 1y ( n ) ＝ y0( n ) ＋ y1( n ) ＋ y2( n ) ＋ …nnnnnnh ( n ) 用 DFT計算分段卷積 yk (n)的方法： （ 1） i=0； L=N＋ M－ 1； 計算并保存 H(k)=DFT［ h(n)] L。 （ 2） 讀入 x i(n)=x(n)R M(n－ kM)， 構造變換區(qū)間［ 0， L1］上的序列 ，實際中就是將 x i (n)的 M個值存放在長度為 M的數(shù)組中 , 并計算 （ 3） （ 4） ， n = 0,1,2,…, L－ 1； （ 5） 計算： （ 6） i =i＋ 1，返回 (2)。 應當說明，一般 x(n)是因果序列，假設初始條件 y1(n)=0。 MATLAB信號處理工具箱中提供了一個函數(shù) fftfilt，該函數(shù)用重疊相加法實現(xiàn)線性卷積的計算。調(diào)用格式為： y=fftfilt(h, x， M)。式中 , h是系統(tǒng)單位脈沖響應向量； x是輸入序列向量； y是系統(tǒng)的輸出序列向量（ h與x的卷積結果）； M是由用戶選擇的輸入序列 x的分段長度，缺省 M時，默認輸入序列 x的分段長度 M=512 【 例 】 假設 h(n)=R5(n)， x(n)=［ cos(πn/10) +cos(2πn/5)］ u(n)，用重疊相加法計算 y(n)=h(n)*x(n)，并畫出 h(n)、 x(n)和 y(n) 解 h(n)的長度為 N=5, 對 x(n)進行分段，每段長度為M=10。計算 h(n)和 x(n)的線性卷積的 MATLAB程序如下： %例 重疊相加法的 MATLAB實現(xiàn)程序： Lx=41。 N=5。 M=10。 %Lx為信號序列 x(n) hn=ones(1, N)。 hn1=［ hn zeros(1, LxN)］ 。 %產(chǎn)生 h(n)，其后補零是為了繪圖好看 n=0:L1。 xn=cos(pi*n/10)+cos(2*pi*n/5)。 %產(chǎn)生 x(n)的 Lx yn=fftfilt(hn, xn, M)。 %調(diào)用 fftfilt用重疊相加法計算卷積 %以下為繪圖部分 , 省略 運行程序畫出 h(n)、 x(n)和 y(n)的波形如圖 示。請讀者從理論上證明 y(n)的穩(wěn)態(tài)波形是單一頻率的 運行繪圖程序 本例題的 xk(n), yk(n)和 y(n)=y0(n)+y1(n)+y2(n)+y3(n), 如圖 。 圖 例 圖 重疊相加法時域波形
教學課件相關推薦
鄂ICP備17016276號-1