通常我們分析的信號是一個純實數(shù)信號,沒有虛部。經(jīng)典的例子當然是用正弦波。當我們用傅立葉變換分析,通常是FFT算法,分析信號時,絕大多數(shù)人都知道我們將得到的頻率結果位于0Hz到奈奎斯特頻率(采樣頻率的一半)之間。同時也明白相應的FFT輸出結果是復數(shù)形式,因此可以用幅值和相位形式或實部和虛部形式表示。
出于解釋目的,假設我們使用的正弦波的頻率為128Hz,幅值為5,初相位為60度。信號的開始部分如圖1所示?,F(xiàn)在假設我們以1024個樣本點每秒的速率對這個信號采集了4s的數(shù)據(jù)。顯然,這些樣本點稍微人為地作了些選擇,為了避免截取成非周期信號,并且要確保信號的頻率正好與FFT的某一條譜線重合。
對這個信號進行常規(guī)的FFT計算,得到幅值和相位形式的結果如圖2所示。
圖2 FFT的半譜
注意到計算得到的幅值是2.5,剛好是原信號幅值的一半。這從數(shù)學上是正常的,將在后面進行說明。然而,一些商業(yè)軟件給出的結果幅值是5,這是因為這些軟件自動對FFT結果進行加倍。我們也注意到相位是-30度,但我們開始使用的正弦波的初相位是60度。但這-30度卻是正確的,這是因為FFT是用余弦定義實部,正弦定義虛部。因此,F(xiàn)FT的0度,用作余弦波在正峰開始時的參考。因此,初相位為60度的正弦波與初相位為-30度的余弦波是相同的。
從數(shù)學上,有
sin(60)=sin(90-30) =sin(90)cos(-30)+cos(90)sin(-30)=cos(-30)
展示的頻率范圍0到奈奎斯特頻率經(jīng)常被稱作半范圍的傅立葉變換。這是一個全范圍的變換,我們可以選擇的范圍從0至采樣頻率,或者從負奈奎斯特頻率到正奈奎斯特頻率(-采樣頻率/2至+采樣頻率/2)。后一種形式實際上更易于解釋,同時包含了負頻率的概念,這似乎初看起來更直觀。實際上的確非常直觀,讓我們回顧一下基本知識。
按常規(guī),正旋轉是按逆時針方向。如果我們觀察一個向量的端點以角頻率ω旋轉,那么相應的軌跡為正弦波Asinωt=Asin2**π ft 。 波形的過去和未來為圖所示。
圖3 逆時針方向旋轉的向量
假設我們現(xiàn)在觀察相同的設置,但現(xiàn)在這個向量按順時針方向旋轉。
圖4 順時針方向旋轉的向量
兩個波形實際上是相同的。因此,如果我們測量這個正弦波,我們不知道波形何時開始,我們也不能確定它是否是順時針方向還是逆時針方向。實際上,正頻率相當于逆時針方向旋轉,負頻率相當于順時針方向旋轉。由測量到的實數(shù)信號,我們不同區(qū)分出正負頻率。
圖5 全范圍的FFT
假設現(xiàn)在我們要對相同的正弦波進行一個從-采樣頻率/2至+采樣頻率/2有全范圍的FFT計算。給出結果如圖5所示。這有兩個分量,二者的幅值為2.5,但一個是負頻率-128,相位為+30度,另一個為+128Hz,相位為-30度。這二者組合得到了開始時的正弦波。對于一個純實數(shù)信號而言,負頻率的幅值是正頻率幅值的鏡像,而相位鏡像成負值。因為我們能從正頻率部分完全預測出負頻率部分的幅值和相位,因而,沒有必要去計算負頻率,因為它們并不包含更多的信息。
從數(shù)學公式上用幅值 A ,相位 ? ,描述一個信號
Ae ^i?^ =A(cos ? +isin?),而?=2πf+θ
因此,我們有正頻率成分
Acos (2πf-30)+i Asin (2πf-30)
和負頻率成分
Acos (-2πf+30)+i Asin (-2πf+30)
因此,整個信號有
Acos (2πf-30)+i Asin (2πf-30)+ Acos
(-2πf+30)+i Asin (-2πf+30)
我們記得
cos (-?)= cos (?)
sin(-?)=-sin(?)
因此,我們可以將負頻率成分寫成
Acos (2πf-30)-i Asin (2πf-30)
將信號的整個表達示進行化簡,得到
2Acos (2πf-30)=2Asin (2πf+30)
這表明,信號的相位與正頻率成分的相位相同,幅值2倍于正頻率成分的幅值。
-
奈奎斯特
+關注
關注
0文章
9瀏覽量
9271 -
傅立葉變換
+關注
關注
3文章
105瀏覽量
32482 -
FFT算法
+關注
關注
0文章
31瀏覽量
13410
發(fā)布評論請先 登錄
相關推薦
評論