數(shù)字信號處理：dsp12 有限離散傅氏變換(III)

1、數(shù)字信號處理7.有限離散傅氏變換(III) 8/11/20221馬盡文7.4應用快速傅氏變換進行頻譜分析 頻譜分析：計算出信號的頻譜，并由此計算出振幅譜、 能譜（或功率譜）、相位譜。A. 頻譜分析的步驟 設連續(xù)實信號為 ，其中 ，以間隔 抽樣得到離散信號 ， 可能是無限長離散信號 也可能為有限長度離散信號進行其頻譜分析的步驟如下：8/11/20222馬盡文7.4應用快速傅氏變換進行頻譜分析一、數(shù)據(jù)準備令 ， 為正整數(shù)。對于有限長度離散信號，取 滿足并且將 改造為：對于無限長離散信號或 非常大的離散信號，我們只能截取 項進行分析，截取信號記為：8/11/20223馬盡文7.4應用快速傅氏變換進行

2、頻譜分析二、用FFT計算頻譜當 時， 表示對應于頻率 的頻譜值。頻譜 是由實部 和虛部 組成的復數(shù)，即三、由頻譜求振幅譜、相位譜，功率譜 其中8/11/20224馬盡文7.4應用快速傅氏變換進行頻譜分析四、對振幅譜或功率譜進行平滑處理 由于在信號中往往含有干擾成分，使振幅譜或功率譜極 不平滑，因此需做平滑處理。 平滑因子 通常取3個點或5個點，如 或8/11/20225馬盡文五、求振幅譜或功率譜的最大值頻率點 和中間值頻率點 7.4應用快速傅氏變換進行頻譜分析 現(xiàn)以振幅譜 來說明如何求 和 ，令 其中 使 這時我們稱 為最大頻率點 令 其中 使8/11/20226馬盡文7.4應用快速傅氏變換進

3、行頻譜分析這時我們稱 為中間值頻率點，其直觀意義是：在頻率范圍 之內，以 為中間分界點，在 的兩邊振幅譜 所占的比重是一樣的。嚴格講， 往往不存在，實際計算時，令B. 頻譜分析中參數(shù)的選取8/11/20227馬盡文7.4應用快速傅氏變換進行頻譜分析對于連續(xù)信號 進行頻譜分析，實際上是對長度為 的離散信號 進行頻譜分析。容易看出，頻譜分析中參數(shù) 和 的選取很重要。 對于連續(xù)信號 ，我們用兩個參數(shù)來刻畫其頻率特點 的截頻或 最高可達的頻率（單位為Hz）； 的頻率分辨間隔（單位為Hz）。對 的振幅譜 取離散值觀察時，離散值的間隔不能大于 ，故當兩個頻率之差大于 時，對 所包含的這兩個頻率成分可分辨開

4、。當 的振幅譜 曲線擺動比較大時，8/11/20228馬盡文7.4應用快速傅氏變換進行頻譜分析 就要取得小些，當 曲線較平滑時， 就可取得大 些。 與 的選取原則： 即 在 部分的能量充分地小。8/11/20229馬盡文7.4應用快速傅氏變換進行頻譜分析由有限離散傅氏變換知，有限離散頻譜的頻率間隔為因此，還要求 滿足從而由于 表示進行頻譜分析的信號記錄長度， 稱為最小記錄長度。8/11/202210馬盡文7.4應用快速傅氏變換進行頻譜分析 ， 的選取規(guī)則：例 已知某信號的截頻 Hz，頻率分辨間隔 Hz 現(xiàn)在要對信號作頻譜分析，問： （1）要求最小記錄長度 等于多少？ （2）抽樣間隔 應滿足什么

5、條件？ （3）抽樣點數(shù) 應滿足什么條件？8/11/202211馬盡文7.4應用快速傅氏變換進行頻譜分析 解：（1） （2） （3） 若要求 為 形式 ，則可取 。 下頁為一個信號作頻譜分析的圖例。8/11/202212馬盡文7.4應用快速傅氏變換進行頻譜分析8/11/202213馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)A. 有限離散哈特利變換設 為實數(shù)，令其中cas為“cos and sin”三個字的縮寫。一、函數(shù)由于則8/11/202214馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)同樣，由有二、有限離散哈特利變換（FDHT）離散信號 ，令則 的離散哈特利變換定義為顯

6、然周期為 。8/11/202215馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)下面求哈特利逆變換由前面問題中cas的含義，知由本章第一節(jié)所證明的等式8/11/202216馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)有則有該式稱為有限離散逆哈特利變換從而有下面有限離散哈特利和逆哈特利變換公式其中，8/11/202217馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)B. 有限離散余弦變換一、有限離散 型余弦變換離散信號：其有限離散 型余弦變換和有限離散 型逆余弦變換為其中8/11/202218馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)二、有限離散 型余弦變換

7、離散信號則其有限離散 型余弦變換和有限離散 型逆余弦變換為其中8/11/202219馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)三、有限離散逆 型余弦變換公式的證明問題：已知（5-1）式，證明（5-2）式。首先給出幾個初等公式8/11/202220馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)取上式的實部和虛部，得將（5-1）式代入（5-2）式右邊，得8/11/202221馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)其中8/11/202222馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)由上式知當 時， 8/11/202223馬盡文7.5有限離散哈特利變換、

8、余弦變換和廣義中值函數(shù)當 時，由 ，故 與 同偶或同奇當 ， 為奇數(shù)時， 也為奇數(shù)，因此由（5-4）式知當 ， 為偶數(shù)時， ，這時有 8/11/202224馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)因此同樣有8/11/202225馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)因此由（5-8）、（5-9）和（5-4）知由（5-5）、（5-6）、（5-7）、（5-10）和（5-4）知有限離散逆 型變換公式成立。8/11/202226馬盡文C. 廣義中值函數(shù)7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)一、廣義中值函數(shù)若對函數(shù) ，存在函數(shù) ，使或者則稱 為廣義中值函數(shù)，其中 為

9、整數(shù)， 為實數(shù)。容易驗證， ， ， 和 都是廣義中值函數(shù)，而且它們的 還是相同的，即 。8/11/202227馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)性質1 為廣義中值函數(shù)， 為正整數(shù)， 為一組離散數(shù)據(jù)，則 其中 ；或者其中 ，上兩式中要求 。8/11/202228馬盡文證明：由中值函數(shù)的定義知7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)因此由上式可知性質1兩式成立。8/11/202229馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)性質2 為廣義中值函數(shù)， 為正偶數(shù)， 為一組離散數(shù)據(jù)，則 其中 ；或者其中 ，上兩式中要求 。8/11/202230馬盡文7.5有限離散

10、哈特利變換、余弦變換和廣義中值函數(shù)證明：根據(jù)自然分解式： 對上式的第二項，利用性質1（用 代替性質1中的 ）就得到性質2中的兩式。 性質2中的兩式把長度為 的變換轉化為兩個長度為 的同一變換，體現(xiàn)了快速二分法的思想。二、一種余弦變換的快速算法設 為 個實數(shù)，當討論快速算法時，取 ， 為正整數(shù)。8/11/202231馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)考慮一種余弦變換 ：令 ， ， ，由于是廣義中值函數(shù)，由性質2中第二式可得其中 。8/11/202232馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)由上式可得其中 8/11/202233馬盡文7.5有限離散哈特利變換、余弦變換和廣義中值函數(shù)由上兩組等式可知， 點的余弦變換 轉換成兩個 點余弦變換 ，這就構成了一種快速算法?，F(xiàn)在計算 點余弦變換的乘法次數(shù) 和加法次數(shù)由（5-11）式知 易知 ， ，由上式可計算 ： 8/11/202234馬盡文7.5有限離散哈特利變換、余弦變換和廣