數(shù)字信號(hào)處理-時(shí)域離散隨機(jī)信號(hào)處理課件：時(shí)頻分析

時(shí)頻分析5.1引言

5.2短時(shí)傅里葉變換

5.3維格納變換(WD)5.4時(shí)域離散信號(hào)的維格納變換

5.5時(shí)頻分布的統(tǒng)一表示式

5.6時(shí)頻分析在編隊(duì)目標(biāo)架次檢測(cè)中的應(yīng)用

5.1引

言

傳統(tǒng)的信號(hào)分析與處理的數(shù)學(xué)工具是傅里葉變換，它的正變換和逆變換分別用下面兩式表示：

(5.1.1)(5.1.2)式中ω是一個(gè)連續(xù)變量,限制了用計(jì)算機(jī)在頻域進(jìn)行分析與處理，而離散傅里葉變換(DFT)將頻域離散化,使之借助計(jì)算機(jī)可以在時(shí)域也可以在頻域?qū)π盘?hào)進(jìn)行分析與處理。由于傅里葉變換物理概念清晰,同時(shí)也是正交變換,因此長(zhǎng)期以來(lái)科技界及各工程領(lǐng)域廣泛使用傅里葉變換和離散傅里葉變換。

X(ejω)稱(chēng)為信號(hào)x(n)的頻譜，它表示了信號(hào)在頻域的分布規(guī)律。也可以用下面公式表示:（5.1.3）

e(ω)稱(chēng)為信號(hào)x(n)的能量譜，它僅包含信號(hào)的幅度信息。但對(duì)于能量無(wú)限信號(hào)，如周期信號(hào)、平穩(wěn)隨機(jī)信號(hào)等，傅里葉變換不存在，可以用功率譜密度（簡(jiǎn)稱(chēng)功率譜）P(ejω)表示:（5.1.4）

式中rxx(m)是x(n)的自相關(guān)函數(shù)。頻譜、能量譜以及功率譜都是信號(hào)變換到頻域的一種表示方法，對(duì)于頻譜不隨時(shí)間變化的確定性信號(hào)以及平穩(wěn)隨機(jī)信號(hào)都可以用它們進(jìn)行分析和處理。5.2短時(shí)傅里葉變換

5.2.1短時(shí)傅里葉變換的定義及其物理解釋

1.短時(shí)傅里葉變換的定義短時(shí)傅里葉變換的定義有兩種形式,下面分別敘述。

(1)定義一:(5.2.1)

式中w(n)是一個(gè)窗函數(shù)，其作用是取出x(n)在n時(shí)刻附近的一小段信號(hào)進(jìn)行傅里葉變換，當(dāng)n變化時(shí)，窗函數(shù)隨n移動(dòng)，從而得到信號(hào)頻譜隨時(shí)間n變化的規(guī)律。此時(shí)的傅里葉變換是一個(gè)二維域(n,ω)的函數(shù)。窗函數(shù)沿時(shí)間軸移動(dòng)情況如圖5.2.1所示。圖

5.2.1窗函數(shù)的移動(dòng)

令n′=n-m

將n′代入定義一中,再將n′用m代替,可得到第二種定義形式。

(2)定義二：

2.短時(shí)傅里葉變換的物理解釋

對(duì)以上STFT的定義形式，從傅里葉變換和線(xiàn)性濾波兩個(gè)角度，可以有兩種不同的物理解釋。

(1)由傅里葉變換角度解釋。按照(5.2.1)式，STFT可以看作n是參變量，x(m)w(n-m)對(duì)m的傅里葉變換，因此它是（n,ω）的函數(shù)。因?yàn)镾TFT是x(m)w(n-m)的傅里葉變換，可以用x(m)和w(n-m)分別的傅里葉變換的卷積表示。設(shè)：那么

如果再將θ改換成-θ，

得到

(5.2.3)上式是STFT定義的一種頻域表示形式。這里如果x(n)是時(shí)變信號(hào)，式中用了它的傅里葉變換，是不合適的，但可以理解為信號(hào)在時(shí)間窗外變?yōu)?以后，取信號(hào)的傅里葉變換;或者說(shuō)是時(shí)間窗內(nèi)的信號(hào)傅里葉變換的平滑形式。（2）

由線(xiàn)性濾波角度解釋。將定義一重寫(xiě)如下：

上式表明,短時(shí)傅里葉變換可以看成x(n)e-jωn與w(n)的線(xiàn)性卷積，如將w(n)看成一個(gè)低通濾波器的單位脈沖響應(yīng)，短時(shí)傅里葉變換則可用圖5.2.2表示。圖5.2.2表明,首先將信號(hào)x(n)調(diào)制到-ω，然后通過(guò)低通濾波器w(n)，其輸出就是短時(shí)傅里葉變換。實(shí)質(zhì)上是將x(n)在ω附近的頻譜搬移到零頻處，作為短時(shí)傅里葉變換。為使其頻率分辨率高，希望w(n)是一個(gè)低通窄帶濾波器，帶外衰減愈大愈好。

利用定義二可以得到線(xiàn)性濾波的另一種物理解釋?zhuān)?/p>

將定義二重寫(xiě)如下:公式中求和號(hào)部分可看成w(n)ejωn與x(n)的線(xiàn)性卷積，因此上式可以寫(xiě)成式中w(n)是低通濾波器，w(n)ejωn就是以ω為中心的帶通濾波器。按照上式，STFT就是信號(hào)首先通過(guò)帶通濾波器，選出以ω為中心的頻譜，再乘以exp(-jωn)，將選出的頻譜搬移到零頻處。

短時(shí)傅里葉變換如按照定義二的物理解釋?zhuān)瑒t可用圖

5.2.3表示。

圖5.2.2定義一的物理解釋圖

5.2.3定義二的物理解釋

5.2.2短時(shí)傅里葉變換的性質(zhì)短時(shí)傅里葉變換是建立在一般傅里葉變換基礎(chǔ)上的一種變換，因此它具有許多和傅里葉變換相似的性質(zhì)。

1.線(xiàn)性性質(zhì)

設(shè)

z(n)=c·x(n)+d·y(n),c,d

為常數(shù),則

(5.2.4)2.頻移性質(zhì)——調(diào)制特性設(shè) ,則

(5.2.5)3.時(shí)移特性設(shè)x(n)=y(n-n0),則

(5.2.6)證明

以上說(shuō)明STFT具有頻移不變性，但不具有時(shí)移不變性，相差一個(gè)相位因子。

4.共軛對(duì)稱(chēng)性

當(dāng)信號(hào)是實(shí)信號(hào)時(shí),短時(shí)傅里葉變換和一般傅里葉變換一樣具有共軛對(duì)稱(chēng)性,即

(5.2.7)因此,其實(shí)部是偶函數(shù)，虛部是奇函數(shù)。

5.由短時(shí)傅里葉變換恢復(fù)信號(hào)由定義(5.2.1)式得到短時(shí)傅里葉變換的反變換為

設(shè)n=m,則

(5.2.8)

只要w(0)≠0,可以由STFTX(n,ω)準(zhǔn)確地恢復(fù)信號(hào)x(n)。

5.2.3短時(shí)傅里葉變換的時(shí)間、頻率分辨率由定義可知，STFT實(shí)際分析的是信號(hào)的局部譜，局部譜的特性決定于該局部?jī)?nèi)的信號(hào)，也決定于窗函數(shù)的形狀和長(zhǎng)度。為了了解窗函數(shù)的影響，假設(shè)窗函數(shù)取兩種極端情況。第一種極端情況是取w(n)=1，

-∞<n<∞,此時(shí)

這種情況下，STFT退化為信號(hào)的傅里葉變換，沒(méi)有任何時(shí)間分辨率，卻有最好的頻域分辨率。第二種極端情況是取w(n)=δ(n)，

此時(shí)

STFT退化為信號(hào)，有理想的時(shí)間分辨率，但不提供任何頻率分辨率。

短時(shí)傅里葉變換由于使用了一個(gè)可移動(dòng)的時(shí)間窗函數(shù)，使其具有一定的時(shí)間分辨率。顯然,短時(shí)傅里葉變換的時(shí)間分辨率取決于窗函數(shù)w(n)的長(zhǎng)度。為了提高信號(hào)的時(shí)間分辨率,希望w(n)的長(zhǎng)度愈短愈好。但頻域分辨率取決于w(n)窗函數(shù)的頻域函數(shù)寬度，也就是低通濾波器w(n)的帶寬或者說(shuō)帶通濾波器w(n)ejωn的帶寬，為了提高頻域分辨率，希望盡量加寬w(n)窗口寬度，這樣必然又會(huì)降低時(shí)域分辨率。因此,

STFT的時(shí)間分辨率和頻率分辨率不能同時(shí)任意提高。這種時(shí)域分辨率和頻域分辨率相互制約的性質(zhì)，也正反映了已為理論所證明了的“不確定原理”：（5.2.9）

式中Δt表示信號(hào)有效持續(xù)時(shí)間，Δf表示信號(hào)的有效帶寬。上面公式說(shuō)明,對(duì)于窗函數(shù)，它的時(shí)間寬度和在頻率域的寬度不能同時(shí)任意小,也就是說(shuō),頻域分辨率和時(shí)域分辨率不能同時(shí)任意小。但可以選擇合適的窗函數(shù)，使Δt和Δf都比較小,其乘積接近于1/(4π)。窗函數(shù)的形式有很多,可以證明從有效時(shí)寬和有效頻寬乘積為最小的意義上講，高斯波形信號(hào)是最好的，但是它在時(shí)間軸和頻率軸上是無(wú)限擴(kuò)張的，因此它并不是一種最好的波形。我們知道，不可能存在既是帶限又是時(shí)限的信號(hào)波形，實(shí)際應(yīng)用中采用放松條件，研究在有限時(shí)寬的情況下，使頻率有效帶寬為最小的波形是什么，或者研究在有限帶寬情況下，使時(shí)寬最小的波形是什么，這部分內(nèi)容可參考文獻(xiàn)［8］、［2］。5.2.4短時(shí)傅里葉變換的實(shí)現(xiàn)

1.短時(shí)傅里葉變換的時(shí)域、頻域采樣

我們已經(jīng)知道,短時(shí)傅里葉變換是低通濾波器w(n)的輸出，假設(shè)其有效帶寬為Δ，對(duì)應(yīng)的模擬濾波器的有效帶寬為B,式中fs是x(n)的采樣頻率，那么對(duì)應(yīng)該模擬濾波器的時(shí)域采樣頻率應(yīng)是帶寬的兩倍(2B)以上,最小采樣頻率為2B即是最小再次采樣率。上式表明該采樣率是信號(hào)采樣率fs的P倍,也就是說(shuō),二次采樣間隔最大為1/P，即窗口每次移動(dòng)的最大間隔是1/P。例如:窗函數(shù)選用哈明窗,長(zhǎng)度為L(zhǎng),帶寬近似為4π/L，設(shè):fs=10kHz,L=100,則

即最大采樣間隔為25，n的取值為0,25,50,75,…。以上計(jì)算的再次采樣率是最小采樣率，采樣間隔是最大采樣間隔。對(duì)于頻率域采樣，假設(shè)在周期2π中采樣M點(diǎn)，為不發(fā)生時(shí)域混疊，

要求M≥L。

2.用FFT計(jì)算STFT假設(shè)在頻率域等間隔采樣M點(diǎn)，

k=0,1,2,3,…,M-1k=0,1,2,3,…,M-1按照定義一，有

（5.2.10）

令m=l+n，

上式變?yōu)?/p>

（5.2.11）

上式的求和區(qū)間是(-∞,∞)，可以按照長(zhǎng)度為M的區(qū)間進(jìn)行劃分。一個(gè)個(gè)區(qū)間計(jì)算后，

再求和，這樣上式變成為

（5.2.12）

令l=m+rM,同時(shí)考慮到 ，得到(5.2.13)

式中

(5.2.14)m=0,1,2,3,…,M-1在(5.2.13)式中，對(duì)任何固定n

值，求和項(xiàng)可以用M點(diǎn)的FFT進(jìn)行計(jì)算，其中信號(hào)用（5.2.14）式計(jì)算。根據(jù)（5.2.13）式和(5.2.14)式,由x(n)計(jì)算STFT的過(guò)程如圖5.2.4所示。圖

5.2.4用FFT計(jì)算短時(shí)傅里葉變換3.用濾波器組法實(shí)現(xiàn)短時(shí)傅里葉變換假設(shè)在頻率域采樣M點(diǎn),采樣點(diǎn)的頻率為

k=0,1,2,3,…,M-1將ωk代入定義(5.2.2)式中，得到

令

(5.2.15)則

(5.2.16)令

(5.2.17)則

(5.2.18)這樣對(duì)應(yīng)M個(gè)采樣點(diǎn)頻率,形成M個(gè)通道。(5.2.15)式即是每個(gè)通道的帶通濾波器的單位取樣響應(yīng),(5.2.17)式即是每個(gè)帶通濾波器的輸出,（5.2.18）式表示每個(gè)通道的STFT輸出。

圖

5.2.5

STFT一個(gè)通道的原理框圖

4.短時(shí)傅里葉變換的綜合

由短時(shí)傅里葉變換恢復(fù)時(shí)域信號(hào)稱(chēng)為綜合。下面先推導(dǎo)短時(shí)傅里葉變換M個(gè)通道的等效傳輸函數(shù)。將M個(gè)帶通濾波器的輸出相加，輸出用y(n)表示,從x(n)到y(tǒng)(n)等效單位取樣響應(yīng)用h(n)表示,(5.2.19)(5.2.20)式中

等效傳輸函數(shù)用H(ejω)表示，即是上式的傅里葉變換。

(5.12.21)式中

(5.12.22)y(n)=x(n)*h(n)(5.2.23)(5.2.24)(5.2.23)式表示等效傳輸函數(shù)是M個(gè)帶通濾波器傳輸函數(shù)的疊加。下面證明公式：

M≥L

(5.2.25)證明

對(duì)于一定的ω，令ω=ωθ,ωk’=ωθ-ωk帶入上式左邊，得到上式最右邊求和號(hào)是頻域信號(hào) 的IDFT，即w(-n)。

另外,我們知道頻域采樣，使時(shí)域以采樣點(diǎn)數(shù)M為周期進(jìn)行周期性延拓，

因此

n=0,1,2,3,…,M-1這里M≥L

保證不出現(xiàn)時(shí)域混疊,令n=0,r=0,帶入上式，最后得到

(證畢)將上式代入(5.2.23)式,得到

H(ejω)=Mw(0),h(n)=F-1T［H(ejω)］=Mw(0)δ(n)

y(n)=x(n)*h(n)=Mw(0)x(n)(5.2.26)(5.2.26)式說(shuō)明當(dāng)M≥L時(shí),用短時(shí)傅里葉變換可以恢復(fù)原時(shí)域信號(hào),恢復(fù)的可能性與窗函數(shù)的形狀無(wú)關(guān)。但是窗函數(shù)的形狀會(huì)影響其時(shí)間與頻率分辨率。由(5.2.18)式得到

(5.2.27)按照(5.2.19)式，

有

(5.2.28)上式表明，由STFT恢復(fù)原信號(hào)，也就是將每一路的STFT輸出乘以exp(jωkn)，然后進(jìn)行相加得到時(shí)域信號(hào)y(n)。下面將x(n)的短時(shí)傅里葉分析和y(n)信號(hào)的恢復(fù)畫(huà)在一起，如圖5.2.6所示。圖

5.2.6STFT分析與綜合的原理圖

我們知道,STFT的輸出的帶寬決定于低通濾波器的帶寬，其帶寬一般比輸入信號(hào)的帶寬小得多，如果輸入信號(hào)的采樣率滿(mǎn)足采樣定理，那么在STFT的輸出完全可以降低采樣率，對(duì)STFT的輸出進(jìn)行二次采樣，減少數(shù)據(jù)量，即運(yùn)算量。為此,在信號(hào)STFT輸出端加一個(gè)抽取器，信號(hào)綜合時(shí)，再進(jìn)行插值，以恢復(fù)原來(lái)的采樣率。這樣,采用濾波器組實(shí)現(xiàn)STFT分析與綜合的原理框圖如圖5.2.7所示,圖中(a)與(b)分別對(duì)應(yīng)定義一和定義二。二次采樣率用輸入信號(hào)的帶寬與低通濾波器的帶寬之比進(jìn)行計(jì)算，設(shè)輸入信號(hào)帶寬用Δx表示，低通濾波器的帶寬用Δω表示，圖中的L用下式計(jì)算：

(5.2.29)

上面我們研究了STFT的計(jì)算方法，這里主要的問(wèn)題是計(jì)算的有效性和時(shí)間窗的選擇。為了提高STFT計(jì)算的有效性，一些學(xué)者研究了遞歸算法，給出了矩形滑動(dòng)窗遞推算法和矩形擴(kuò)展窗遞推算法（矩形窗的寬度作為遞歸形式進(jìn)行擴(kuò)展）;為了避免矩形窗帶來(lái)的截?cái)嘈?yīng)，又給出哈明（Hamming）窗和漢寧（Hanning）窗的改進(jìn)遞歸算法，這樣又增加了計(jì)算的復(fù)雜性。文獻(xiàn)［3］給出了一種新的遞歸算法，它采用了全極點(diǎn)的滑動(dòng)窗，它不僅克服了采用矩形窗帶來(lái)的截?cái)嘈?yīng)，還比用哈明窗和漢寧窗的計(jì)算量和存儲(chǔ)量少得多。這方面內(nèi)容請(qǐng)參考文獻(xiàn)［2］、［3］。

圖

5.2.7用濾波器組實(shí)現(xiàn)STFT的一個(gè)通道

5.3維格納變換(WD)5.3.1WD的定義

確定性時(shí)間連續(xù)信號(hào)的WD定義如下:(5.3.1)定義表明這種變換是把過(guò)去某一時(shí)間信號(hào)乘上未來(lái)某一時(shí)間信號(hào)，再對(duì)兩個(gè)信號(hào)時(shí)間差τ求傅里葉變換得到。

令

(5.3.2)將rxx(t,τ)稱(chēng)為瞬時(shí)自相關(guān)函數(shù)，那么WD就是信號(hào)瞬時(shí)自相關(guān)函數(shù)的傅里葉變換。x(t)在頻率域的WD分布定義如下:(5.3.3)對(duì)于兩個(gè)連續(xù)時(shí)間信號(hào)x(t)與y(t),互WD定義為

(5.3.4)同樣,它們?cè)陬l率域的互WD定義如下:(5.3.5)式中,X(Ω)是x(t)的傅里葉變換;Y(Ω)是y(t)的傅里葉變換。(5.3.1)、（5.3.4）式是在時(shí)域的定義形式，(5.3.3)、（5.3.5）式是在頻域的定義形式。可以證明,時(shí)域和頻域的定義形式有下面關(guān)系：(5.3.6)(5.3.7)5.3.2WD的性質(zhì)

1.WD的實(shí)數(shù)性和對(duì)稱(chēng)性

(1)WD是t和Ω的實(shí)函數(shù)。

證明

對(duì)(5.3.1）式的兩邊取共軛,得到

令τ′=-τ,則

因此

(5.3.8)(2)如果x(t)是實(shí)信號(hào)，則WD是頻率的偶函數(shù)。

即

(5.3.9)證明

將定義中的Ω換成-Ω,得到

因此

(3)對(duì)于互WD,則具有如下性質(zhì):(5.3.10)

2.邊緣積分特性

(1)在固定時(shí)刻t，WD沿全頻率軸的積分等于在t時(shí)刻信號(hào)的瞬時(shí)功率Px(t)(也稱(chēng)時(shí)間邊緣特性)，即（5.3.11）

證明

由定義(5.3.1)式,得到

令

那么

令t1=t2=t,τ=0,則得到

(2)在固定頻率Ω,WD沿全時(shí)間軸的積分等于該頻率的能量密度Px(Ω)(也稱(chēng)頻率邊緣特性)，即(5.3.12)證明

由定義(5.3.3)式，得到

令

,那么

代入上式，

得到

令Ω1=Ω2=Ω,因此

(3)WD分布在整個(gè)（t，Ω）平面上,對(duì)t，Ω的雙重積分等于信號(hào)的總能量E，即

(5.3.13)利用(5.3.11)式和(5.3.12)式,可以推出

(5.3.14)

用以上性質(zhì)可對(duì)WD進(jìn)行能量化解釋,時(shí)間邊緣特性為信號(hào)的瞬時(shí)功率,頻率邊緣特性為信號(hào)的能譜密度,總能量E將Px(t)和Px(Ω)聯(lián)系起來(lái),因此WD是一種能量化的時(shí)頻表示,但不能把WD解釋為在時(shí)間—頻率上每一點(diǎn)的時(shí)頻能量密度,因?yàn)閃D有時(shí)可能是負(fù)的,且由不確定原理不允許在某個(gè)特定的時(shí)間—頻率處有能量這一概念。

3.WD的運(yùn)算性質(zhì)

(1)時(shí)移與頻移的不變性。將x(t)時(shí)移τ，相應(yīng)的其WD分布也時(shí)移τ，用公式表示如下:

如果x(t)=y(t-τ)，那么（5.3.15）

x(t)用調(diào)制，其WD分布也頻移Ω0，用公式表示如下:

如果 ,那么（5.3.16）

(2)兩信號(hào)的時(shí)域卷積等于兩信號(hào)分別的WD在時(shí)間軸上的卷積，即:

如果y(t)=x(t)*h(t),則

（5.3.17）

(3)如果兩信號(hào)的相乘，它們的傅里葉變換服從卷積關(guān)系，則和它們對(duì)應(yīng)的WD在頻率軸上也服從卷積關(guān)系，用公式表示如下:

如果y(t)=x(t)h(t),則

如果x(t)表示信號(hào)，h(t)表示窗函數(shù)，此性質(zhì)表明信號(hào)加窗處理時(shí)，只影響頻率分辨率，不影響時(shí)間分辨率。（5.3.18）

(4)兩信號(hào)相加,設(shè)z(t)=x(t)+y(t),則

式中第三項(xiàng)稱(chēng)為交叉(干擾)項(xiàng)，

其性質(zhì)意義在后面介紹。

（5.3.19）

4.WD的時(shí)限性和帶限性——區(qū)域性信號(hào)的維格納分布的時(shí)寬與頻寬，與信號(hào)本身的時(shí)寬與頻寬相同，即:若x(t)限制在t1≤t≤t2中，則它的WD分布也限制在同一時(shí)間域中；若x(t)的傅里葉變換X(jΩ)限制在Ω1≤Ω≤Ω2中，則它的WD分布也限制在同一頻率域中，用公式表示如下：如果

則

（5.3.20）

如果

則

（5.3.21）

利用該性質(zhì)，又可以得到下面結(jié)論：（1）

因果信號(hào)x(t)的WD也是因果的，即:若

t≥0t<0則

（5.3.22）

（2）解析信號(hào)z(t)的傅里葉變換限制在頻率的正半軸，解析信號(hào)z(t)的WD也限制在Ω≥0的上半平面，

即

（5.3.23）

5.可逆性

由定義得到

令

則

代入上式，

得到

令t2=0,再將t1用t代替，得到

上式說(shuō)明,信號(hào)可以由其WD分布進(jìn)行重建，僅僅缺少初相位信息。

(5.3.24)

5.3.3常用信號(hào)的WD舉例 例

5.3.1求其WD。

解

下面確定對(duì)τ的積分限：

因此

|t|<T

|t|>T

(5.3.25)

上式表明WD在時(shí)間軸上限制在-T～T之間，在頻率軸上是sinx/x形式，最大值在（T,0）處。

波形圖如圖5.3.1所示。

圖

5.3.1例

5.3.1圖

例5.3.2 ,求其WD。解

該例題的信號(hào)是一個(gè)復(fù)正弦信號(hào)，可以看作平穩(wěn)隨機(jī)信號(hào)，其WD分布與時(shí)間無(wú)關(guān)，對(duì)任意時(shí)間都是一個(gè)在Ω=Ω0處的δ函數(shù)，

如圖

5.3.2所示。

(5.3.26)圖

5.3.2例

5.3.2圖

例5.3.3

已知x(t)=Acos(Ω0t),求其WD。

解按照例5.3.2和(5.3.19)式，推導(dǎo)如下:上面結(jié)果中的第一項(xiàng)表示信號(hào)，第二項(xiàng)即是交叉干擾項(xiàng)。交叉干擾項(xiàng)將在后面介紹。x(t)的WD如圖

5.3.3所示。

圖

5.3.3例

5.3.3圖

例5.3.4

已知 ,求其WD。解

(5.3.28)

x(t)是一個(gè)線(xiàn)性調(diào)頻信號(hào)，其WD清楚地表示出功率譜隨時(shí)間線(xiàn)性變化的性質(zhì)。WD如圖

5.3.4所示。

圖

5.3.4例

5.3.4圖

例5.3.5

已知 ,求互WD。

解

(5.3.29)例5.3.6已知高斯信號(hào),求其WD。解

上式表明，

高斯信號(hào)的WD在時(shí)間上和頻率上有相同的波形。

(5.3.30)5.3.4關(guān)于二次時(shí)頻分布中的交叉(干擾)項(xiàng)前面曾說(shuō)過(guò)WD是在時(shí)間、頻率二維域中的能量分布函數(shù)，必然是信號(hào)的二次型，或者說(shuō)是一種雙線(xiàn)性變換。這里雙線(xiàn)性變換指的是x(t+τ/2)x*(t-τ/2)的傅里葉變換。這樣就使能量分布函數(shù)不服從線(xiàn)性疊加原理。我們知道,短時(shí)傅里葉變換是一種線(xiàn)性變換，它可以表示信號(hào)頻譜隨時(shí)間變化的規(guī)律。如果取它的模的平方，則可以粗略表示信號(hào)在時(shí)間、頻率二維域中的能量分布。一般把短時(shí)傅里葉變換模的平方稱(chēng)為譜圖。譜圖是一種能量分布函數(shù)，也不服從線(xiàn)性疊加原理，兩個(gè)信號(hào)之和的譜圖并不等于它們分別的譜圖的和，還存在第三項(xiàng)即交叉項(xiàng)。下面推導(dǎo)說(shuō)明。假設(shè)信號(hào)x(t)譜圖用SPECx(t,Ω)表示,（5.3.31）

設(shè)

式中,a,b為實(shí)常數(shù)。

式中

(5.3.32)(5.3.32)式表明,兩信號(hào)相加的譜圖并不等于兩個(gè)信號(hào)的譜圖之和，其譜圖包括三部分，第一部分和第二部分分別是兩個(gè)信號(hào)的譜圖，第三部分即是交叉項(xiàng)。交叉項(xiàng)的最大幅度是兩個(gè)信號(hào)譜圖幅度乘積的兩倍，而且交叉項(xiàng)幅度受到一個(gè)余弦函數(shù)的調(diào)制，其相位變化服從兩個(gè)譜圖相位差的變化規(guī)律。如果兩個(gè)信號(hào)自身的譜圖沒(méi)有重疊部分，則第三部分為零，即沒(méi)有交叉項(xiàng)，

這是譜圖的優(yōu)點(diǎn)。

事實(shí)上,任何二次時(shí)頻分布都不服從線(xiàn)性疊加原理，而服從二次疊加原理。下面介紹二次疊加原理。信號(hào)x(t)的二次時(shí)頻分布用Tx(t,Ω)表示，設(shè)

（5.3.33）

則

（5.3.34）

式中：和分別稱(chēng)為x1(t)和x2(t)的自時(shí)譜;和 分別稱(chēng)為x1(t)對(duì)x2(t)和x2(t)對(duì)x1(t)的互時(shí)譜。這種互時(shí)譜形成了二次時(shí)頻分布的交叉項(xiàng)。

下面再分析WD中的交叉項(xiàng)，

將(5.3.33)式帶入信號(hào)的雙線(xiàn)性變換中，

得到

(5.3.35)式中，第一、二項(xiàng)形成自時(shí)譜，第三、四項(xiàng)則形成互時(shí)譜，即交叉項(xiàng)。對(duì)照(5.3.34)式，WD服從二次疊加原理。對(duì)于有N個(gè)分量的信號(hào)，

二次疊加原理用下式表示：

設(shè)

,則

（5.3.36）

k≠l

例5.3.7已知 ,式中c1,c2為實(shí)數(shù)，求其WD。解

(5.3.37)

該例表明，兩個(gè)復(fù)單頻信號(hào)頻譜不重疊（Ω1≠Ω2）時(shí)，它的時(shí)頻分布除了自時(shí)譜以外，仍有交叉項(xiàng);交叉項(xiàng)處在兩個(gè)信號(hào)頻率連線(xiàn)的中點(diǎn)，幅度受到兩信號(hào)頻率差值的余弦波調(diào)制，不能保證時(shí)頻分布是非負(fù)函數(shù)，其最大幅度是兩信號(hào)幅度乘積的4π倍。

例5.3.3是一個(gè)高頻實(shí)信號(hào)，由于分成了兩個(gè)復(fù)指數(shù)信號(hào)，即存在負(fù)頻率分量，在零頻率處形成交叉干擾項(xiàng)。因解析信號(hào)沒(méi)有負(fù)頻率分量，先將實(shí)信號(hào)轉(zhuǎn)變成解析信號(hào)，再進(jìn)行WD分析，可消除這種頻譜正負(fù)部分之間的交叉干擾項(xiàng)。下面對(duì)解析信號(hào)作簡(jiǎn)單介紹。假設(shè)x(t)是實(shí)的連續(xù)時(shí)間信號(hào)，用z(t)表示對(duì)應(yīng)的解析信號(hào)。

z(t)定義為

（5.3.38）

式中是x(t)的Hilbert變換，或者說(shuō)是x(t)通過(guò)一個(gè)Hilbert變換器形成的。Hilbert變換器的傳輸函數(shù)為(5.3.39)或者

H(jΩ)=-jsgn(Ω)(5.3.40)式中

(5.3.41)Hilbert變換器的單位沖激響應(yīng)為

(5.3.42)x(t)通過(guò)Hilbert變換器后，其輸出為

(5.3.43)這樣解析信號(hào)z(t)和x(t)之間的關(guān)系為

(5.3.44)將上式進(jìn)行傅里葉變換，

得到

Z(jΩ)=X(jΩ)+jH(jΩ)X(jΩ)=X(jΩ)［1+jH(jΩ)］

(5.3.45)將（5.3.39）式帶入上式，

得到

(5.3.46)上式表明,解析信號(hào)的頻譜只分布在正頻率范圍，是由實(shí)信號(hào)頻譜的正的部分乘以2構(gòu)成的;負(fù)頻率部分為0。

用上式可以求解析信號(hào)。例如ejΩt，其頻譜是在Ω處的δ函數(shù)，如果Ω是負(fù)的，那么也就沒(méi)有正頻率存在，因此它的解析信號(hào)為

Ω<0Ω>0例

5.3.8

求

f0≠0的解析信號(hào)。

解觀(guān)察（5.3.39）式，它的Hilbert變換當(dāng)Ω>0時(shí)，需將信號(hào)的相位變化-90°，幅度不變，因此 ,那么，x(t)的解析信號(hào)為

或者由

取其正頻率部分的兩倍，得到的解析信號(hào)和上式一樣。

5.4時(shí)域離散信號(hào)的維格納變換

5.4.1時(shí)域離散信號(hào)的WD定義

1.Classen提出的定義

按照連續(xù)時(shí)間信號(hào)維格納變換定義（5.3.1）式，可以引申出Classen提出的時(shí)域離散信號(hào)的WD定義:（5.4.1）

令k′=2k，

得到

（5.4.2）

注意式中的ω是數(shù)字頻率，但頻域的重復(fù)周期不是2π，而是π。因此下面公式成立:

（5.4.3）

對(duì)應(yīng)的頻域定義為

（5.4.4）

Classen提出的定義應(yīng)用較廣，保持了連續(xù)時(shí)間信號(hào)WD定義中有關(guān)時(shí)間上的一些特性，但頻域上的一些特性被破壞了，

例如,時(shí)間域的邊緣特性為

（5.4.5）

它不同于連續(xù)時(shí)間信號(hào)WD在時(shí)間上的積分等于信號(hào)在某一頻率的能量密度((5.3.12)式)。

2.Peyrin提出的定義將（5.4.1）式中的系數(shù)2去掉，

寫(xiě)成下式:(5.4.6)令:n+m′/2=k,n-m′/2=n′-k,則k=(n′+m′)/2,這樣(5.4.6)式變成

令:n+m′/2=k,n-m′/2=n′-k,則k=(n′+m′)/2,這樣(5.4.6)式變成

(5.4.7)上式就是Piyrin提出的定義。它與Classen提出的定義相比，下面關(guān)系式成立:（5.4.8）

式中,上標(biāo)c代表Classen;上標(biāo)p代表Peyrin。該式表明等式左面的點(diǎn)只是等式右邊時(shí)間變量為偶數(shù)點(diǎn)的結(jié)果。因此,Peyrin提出的定義中包含有更多的信息，其中時(shí)間變量為奇數(shù)點(diǎn)的信息是Classen定義中缺少的。Peyrin提出的定義還有其它性質(zhì)和優(yōu)點(diǎn)，請(qǐng)參考文獻(xiàn)［2］。

5.4.2利用FFT計(jì)算維格納分布維格納分布的計(jì)算量很大，目前各種快速算法還不能從根本上解決實(shí)時(shí)處理的問(wèn)題。下面介紹用FFT計(jì)算維格納分布的方法。如果用離散哈特萊變換（DHT）計(jì)算，計(jì)算的復(fù)雜性可由三倍復(fù)FFT減少到三倍實(shí)FFT的計(jì)算量，這部分內(nèi)容可參考文獻(xiàn)［2］。

離散時(shí)間信號(hào)的WD為

(5.4.9)

為用FFT進(jìn)行計(jì)算，需對(duì)信號(hào)進(jìn)行加窗處理，并且將頻率域離散化。假設(shè)窗函數(shù)用w(l)表示，它的時(shí)寬為2L-1。且當(dāng)|l|≥L時(shí)，

w(l)=0。加窗后的WD為

(5.4.10)

WDx(n,ω)的頻域周期是π，若在一個(gè)周期內(nèi)采樣N點(diǎn)(N=2L-1)，采樣間隔為Δω=π/N，為便于計(jì)算，再在尾部加個(gè)

0，使N=2L，令

(5.4.11)由于FFT的計(jì)算域是正的（l=0,1,…,N-1)，重新排列序列如下:l=0,1,…,L-1l=L,…,2L-1(5.4.12)最后得到

(5.4.13)上式就是用FFT計(jì)算WD的公式。

5.5時(shí)頻分布的統(tǒng)一表示式5.5.1模糊函數(shù)及其和WD之間的關(guān)系模糊函數(shù)也是一種常用的時(shí)頻表示，它廣泛應(yīng)用于雷達(dá)、水聲等領(lǐng)域。本節(jié)主要介紹它的定義及其和WD之間的關(guān)系。

WD分布是對(duì)信號(hào)的雙線(xiàn)性變換x(t+τ/2)x*(t-τ/2)關(guān)于τ作傅里葉變換，如果對(duì)該雙線(xiàn)性變換關(guān)于時(shí)間t作傅里葉變換，則得到模糊函數(shù)的定義，公式為（5.5.1）

另外，對(duì)應(yīng)WD的頻率域定義（5.3.3）式，模糊函數(shù)在頻率域的定義是

（5.5.2）

而且

（5.5.3）

式中,X(f)是x(t)的傅里葉變換;t為時(shí)間，τ為時(shí)延;f為頻率;θ為頻偏。為分析模糊函數(shù)和WD之間的關(guān)系，定義:（5.5.4）

（5.5.5）

rx(t,τ)稱(chēng)為瞬時(shí)自相關(guān)函數(shù)，Rx(f,θ)稱(chēng)為瞬間頻自相關(guān)函數(shù)。這樣,WD是瞬時(shí)自相關(guān)函數(shù)關(guān)于τ的傅里葉變換，模糊函數(shù)是瞬時(shí)自相關(guān)函數(shù)關(guān)于時(shí)間t的傅里葉變換。按照定義，模糊函數(shù)和WD還可以表示成以下各式：

（5.5.6）

（5.5.7）

（5.5.8）

（5.5.9）

按照定義，可以證明模糊函數(shù)具有以下性質(zhì):

（1）

時(shí)移性。令x(t)=y(t-t0),則

（5.5.10）

(2)頻移性。令

,則

（5.5.11）

(3)濾波。令

（5.5.12）

則

（5.5.13）

(4)調(diào)制。令y(t)=x(t)m(t),則

(5.5.14)下面推導(dǎo)信號(hào)的WD和模糊函數(shù)之間的關(guān)系。

將(5.5.6)式重寫(xiě)如下:(5.5.15)對(duì)上式進(jìn)行傅里葉反變換，得到瞬時(shí)自相關(guān)函數(shù)和WD之間的關(guān)系:(5.5.16)再對(duì)上式t作傅里葉變換，得到模糊函數(shù)和WD之間的關(guān)系:（5.5.17）

按照(5.5.8)式，還可以得到下式:（5.5.18）

(5.5.17)式和(5.5.18)式表明了WD和模糊函數(shù)之間存在如公式那樣的二維傅氏變換關(guān)系。另外,還可以推導(dǎo)出瞬時(shí)自相關(guān)和瞬間頻自相關(guān)之間的關(guān)系，如下式:(5.5.19)(5.5.20)

總結(jié)以上，(5.5.6)式～(5.5.9)式和(5.5.17)式～(5.5.20)式，將rx(t,τ)、WDx(t,f)、Rx(θ,τ)以及Ax(θ,τ)聯(lián)系起來(lái)，如圖5.5.1所示。圖中 表示將τ映射為f的傅氏變換，反過(guò)來(lái),則表示由f映射為τ的傅氏反變換，其它類(lèi)似。圖5.5.1rx(t,τ)、WDx(t,f)、Rx(θ,f)、Ax(θ,τ)之間的關(guān)系

這四個(gè)函數(shù)有四個(gè)變量，即時(shí)間變量t、時(shí)間延遲τ、頻率f、頻偏θ，共形成了四個(gè)域，即：（1）時(shí)頻域(t,f)，對(duì)應(yīng)WDx(t,f);

（2）瞬時(shí)相關(guān)域（t,τ），對(duì)應(yīng)rx(t,τ);

（3）譜相關(guān)域（θ,f），對(duì)應(yīng)Rx(f,θ);

（4）模糊域（θ,τ），對(duì)應(yīng)Ax(τ,θ)。WDx(t,f)和Ax(τ,θ)是信號(hào)的兩個(gè)不同的時(shí)頻表示，按照上述分析，還應(yīng)注意它們下面的兩個(gè)不同點(diǎn)：（1）WD是能量化的時(shí)頻表示，存在時(shí)間邊緣特性Px(t)和頻率邊緣特性Px(Ω)，公式重寫(xiě)如下:(5.5.21)(5.5.22)Px(t)稱(chēng)為信號(hào)的瞬時(shí)功率，Px(f)稱(chēng)為信號(hào)的能譜密度，信號(hào)的總能量為

(5.5.23)模糊函數(shù)是相關(guān)化的時(shí)頻表示，將模糊函數(shù)的定義重寫(xiě)如下:(5.5.24)(5.5.25)在(5.5.24)式中,令θ=0,得到時(shí)間相關(guān)化邊緣特性rx(τ),(5.5.26)在(5.5.25)式中,令τ=0,得到頻率相關(guān)化邊緣特性Rx(θ),rx(τ)稱(chēng)為瞬時(shí)相關(guān),Rx(θ)稱(chēng)為譜相關(guān)，因此模糊函數(shù)將瞬時(shí)相關(guān)和譜相關(guān)兩個(gè)概念綜合在一起。

（2）WD滿(mǎn)足時(shí)頻移不變性質(zhì)，即滿(mǎn)足（5.3.15）式和(5.3.16)式,或者用下式表示:如果 ,則

(5.5.28)而模糊函數(shù)滿(mǎn)足相關(guān)化移不變性質(zhì),用公式表示如下:

如果 ,則

(5.5.29)相關(guān)化移不變性質(zhì)來(lái)源于瞬時(shí)相關(guān)和譜相關(guān)移位性質(zhì)得到的術(shù)語(yǔ),用公式表示如下:如果 ,則

(5.5.30)總結(jié)起來(lái),WD和模糊函數(shù)是信號(hào)的兩種不同的時(shí)頻表示,性質(zhì)不同,相互關(guān)系是(5.5.17)、(5.5.18)式表示的二維傅氏變換對(duì)。這兩種時(shí)頻表示的一些對(duì)偶關(guān)系如表5.5.1所示。

表

5.5.1WD和模糊函數(shù)的一些對(duì)偶關(guān)系

5.5.2Cohen類(lèi)時(shí)頻分布

維格納分布和模糊函數(shù)是實(shí)際中比較常用的兩種時(shí)頻分析,此外還有許多種時(shí)頻分析,盡管它們的形式和性質(zhì)不盡相同,但彼此常有一定的聯(lián)系和共同點(diǎn),因此可用一個(gè)統(tǒng)一的表達(dá)式表示出來(lái),這就是L.Cohen提出的廣義雙線(xiàn)性時(shí)頻表示,公式為

(5.5.31)式中φ(θ,τ)表示核函數(shù)。采用不同的核函數(shù)可以得到不同的時(shí)頻分布，時(shí)頻分布的各種性質(zhì)要求,則反映在對(duì)核函數(shù)的約束條件上。把滿(mǎn)足（5.5.31）式的時(shí)頻分布，統(tǒng)稱(chēng)為Cohen類(lèi)時(shí)頻分布。

WD時(shí)頻分布是Cohen類(lèi)時(shí)頻分布中最基本的時(shí)頻分布。下面證明當(dāng)核函數(shù)φ(θ,τ)=1時(shí)，

Cohen類(lèi)時(shí)頻分布將轉(zhuǎn)換成WD。

將φ(θ,τ)=1帶入（5.5.31）式,得到

下面介紹時(shí)頻分布的性質(zhì)對(duì)核函數(shù)要求的約束條件：（1）如要求Cx(t,Ω)=C*x(t,Ω)，即要求時(shí)頻分布是實(shí)的，充要條件是要求核函數(shù)滿(mǎn)足下式:

（2）

如要求時(shí)頻分布具有時(shí)間邊緣特性和頻率邊緣特性，

即若要求：

則要求：

φ(θ,0)=1若要求

則要求：

（3）

如要求:

則要求：

（4）

如要求時(shí)移不變和頻移不變，

即若要求：

則要求φ(θ,τ)獨(dú)立于t和f。

（5）

如能抑制Cx(t,f)的交叉項(xiàng)，則要求φ(θ,τ)為低通濾波。

（6）如Cx(t,f)具有時(shí)間支持特性，即當(dāng)|t|>ta時(shí)，x（t）=0，Cx(t,f)=0,則要求:

（7）如Cx(t,f)具有頻率支持特性，即當(dāng)|f|>fc時(shí),X(f)=0,Cx(t,f)=0,則要求:|θ|<2|fc|

下面介紹Cohen類(lèi)時(shí)頻分布的四種表示形式。

(1)二維濾波的維格納分布。如果對(duì)核函數(shù)φ(θ,τ)作二維傅里葉變換,得到

(5.5.32)Φ(t,f)稱(chēng)為時(shí)頻域核函數(shù),φ(θ,τ)稱(chēng)為模糊域核函數(shù)。由上式得到

(5.5.33)又由WD的定義得到

(5.5.34)將(5.5.33)式和(5.5.34)式帶入(5.5.31)式,得到

(5.5.35)上式就是維格納分布與時(shí)頻域核函數(shù)Φ(t,f)的二維褶積,也稱(chēng)為廣義維格納分布,因此Cx(t,f)類(lèi)時(shí)頻分布可以理解為二維濾波的維格納分布。(2)廣義模糊函數(shù)M(θ,τ)的M(-θ,τ)的二維傅里葉變換。

(5.5.36)

式中

(5.5.37)將上式帶入(5.5.31)式,得到(5.5.36)式。M（θ,τ）稱(chēng)為廣義模糊函數(shù),它是由模糊域核函數(shù)φ(θ,τ)對(duì)Ax(θ,τ)加權(quán)所得的模糊函數(shù)。因此,(5.5.36)式可理解為廣義模糊函數(shù)M(θ，τ)的M(-θ,τ)的二維傅里葉變換。

(3)廣義自相關(guān)函數(shù)對(duì)時(shí)延τ的一維傅里葉變換。

(5.5.38)式中

(5.5.39)rx′(t,τ)稱(chēng)為廣義自相關(guān)函數(shù),而

(5.5.40)ψ(t,τ)稱(chēng)為瞬時(shí)相關(guān)域的核函數(shù)。

由(5.4.40)式得到

(5.5.41)(5.5.39)式和(5.5.41)式表明廣義自相關(guān)函數(shù)為信號(hào)的瞬時(shí)自相關(guān)函數(shù)rx(t,τ)與瞬時(shí)相關(guān)域的核函數(shù)ψ(t,τ)對(duì)t的褶積,因此(5.5.38)式可以理解為廣義自相關(guān)函數(shù)rx′(t,τ)對(duì)時(shí)延τ的一維傅氏變換。

將（5.5.39）式帶入(5.5.38)式,得到

(5.5.42)(4)廣義譜相關(guān)函數(shù)的對(duì)θ的一維傅里葉變換。

式中

(5.5.43)(5.5.44)Rx′(θ,f)稱(chēng)為廣義譜相關(guān)函數(shù),而

(5.5.45)Ψ(θ,f)稱(chēng)為譜相關(guān)域的核函數(shù),由上式得到

(5.5.46)譜相關(guān)函數(shù)Rx(θ,ξ)為

(5.5.47)因此

(5.5.48)

上面四種等價(jià)的表示形式可歸納為對(duì)維格納分布、模糊函數(shù)，瞬時(shí)相關(guān)函數(shù)和瞬間頻相關(guān)函數(shù)的四種表示形式。它們是

(5.5.49)(5.5.50)(5.5.51)(5.5.52)四種核函數(shù)之間的關(guān)系為

(5.5.53)(5.5.54)(5.5.55)四種核函數(shù)之間的關(guān)系如圖5.5.3所示。

圖

5.5.2Cx(t,f)在(t,τ)、

(θ,f)、

(θ,τ)及(t,f)域之間的關(guān)系

圖

5.5.3四種核函數(shù)之間的關(guān)系

WD是能量化的時(shí)頻表示，而模糊函數(shù)是相關(guān)化的時(shí)頻表示，相應(yīng)的前者是時(shí)頻域（t，f），后者是模糊域（θ,τ）。相應(yīng)的可以將Cohen類(lèi)時(shí)頻分布分成兩類(lèi)時(shí)頻表示，即能量化時(shí)頻表示和相關(guān)化時(shí)頻表示。能量化時(shí)頻表示將瞬時(shí)功率和譜能量密度兩種概念綜合在一起，能量化的解釋用邊緣特性(5.5.21)式和(5.5.22)式表示，它的原始公式是(5.5.31)式，(5.5.49)～(5.5.52)式是它的四種等價(jià)形式。能量化時(shí)頻表示簡(jiǎn)稱(chēng)CE類(lèi)(下標(biāo)E表示能量化)。相關(guān)化時(shí)頻表示將瞬時(shí)相關(guān)和相關(guān)綜合在一起，并用相關(guān)化邊緣特性描述(即(5.5.26)、(5.5.27)式)，相關(guān)化時(shí)頻表示簡(jiǎn)稱(chēng)Cc類(lèi)(下標(biāo)c表示相關(guān)化)。CE類(lèi)時(shí)頻表示Cx(t,f)中任何一種等價(jià)形式的對(duì)偶相關(guān)化時(shí)頻表示Pdual,x(θ,τ)均有下面公式：

(5.5.56)CE類(lèi)和Cc類(lèi)的傅氏變換關(guān)系用圖

5.5.4表示。

圖

5.5.4CE類(lèi)和Cc類(lèi)的傅氏變換關(guān)系

5.5.3廣義雙線(xiàn)性時(shí)頻分布舉例

1.指數(shù)分布（ED）這種廣義雙線(xiàn)性時(shí)頻分布的核函數(shù)是指數(shù)形式的，因此稱(chēng)為指數(shù)分布。又因?yàn)槭怯蒀hoi

Williams提出的，也稱(chēng)Choi

Williams分布（CWD）。這種分布對(duì)含有多個(gè)頻率分量的信號(hào)，可以有效地抑制交叉項(xiàng)，并保持時(shí)頻分布的一些期望特性，是一種常用的廣義雙線(xiàn)性時(shí)頻分布，下面作簡(jiǎn)單介紹。

指數(shù)分布在（θ,τ）域的核函數(shù)為

(5.5.57)在(t,τ)域的核函數(shù)為

(5.5.58)式中σ>0，是尺度因子。將上式帶入（5.5.38）、(5.5.39)式中,得到它的時(shí)頻分布Ex(t,f)為

(5.5.59)下面說(shuō)明按照核函數(shù)的約束條件,這種時(shí)頻分布具有的性質(zhì)。(1)按照(5.5.57)式,核函數(shù)滿(mǎn)足下式:因此該時(shí)頻分布是實(shí)值的。

（2）

按照(5.5.57)式,核函數(shù)滿(mǎn)足下列式:因此指數(shù)分布具有邊緣特性和能量特性，

即

（3）信號(hào)的時(shí)移或頻移在指數(shù)分布中產(chǎn)生相應(yīng)的時(shí)移或頻移，這是因?yàn)樗械腃E類(lèi)廣義時(shí)頻分布均有這種性質(zhì)。

（4）指數(shù)分布對(duì)時(shí)間t的一階矩和對(duì)頻率的一階矩分別等于群延遲T（Ω）和瞬時(shí)頻率Ω（t）特性。這是因?yàn)楹撕瘮?shù)滿(mǎn)足

（5）指數(shù)分布可以有效地抑制多頻率成分信號(hào)的交叉項(xiàng)。假設(shè)信號(hào)是一個(gè)多頻率信號(hào)，用下式表示:(5.5.60)可以推出它的廣義模糊函數(shù)用下式表示:(5.5.61)上式中第一項(xiàng)是廣義自模糊函數(shù)項(xiàng)，第二項(xiàng)是廣義互模糊函數(shù)項(xiàng)，

它們分別用下式表示：l≠m(5.5.62)(5.5.63)廣義自模糊函數(shù)主要集中在原點(diǎn)附近,而廣義互模糊函數(shù)或者說(shuō)交叉項(xiàng)則離開(kāi)原點(diǎn)一段距離。這樣,如果要突出廣義自模糊函數(shù)項(xiàng)，并且抑制廣義互模糊函數(shù)項(xiàng)，則要求核函數(shù)在（θ,τ）平面的原點(diǎn)附近有較大的加權(quán)值，對(duì)交叉項(xiàng)的加權(quán)值盡可能小，從而達(dá)到抑制交叉項(xiàng)的目的。指數(shù)分布的核函數(shù)則具備這種性質(zhì)，因此它可以有效地抑制交叉項(xiàng)。假設(shè)兩頻率信號(hào)如下式所示：

(5.5.64)該信號(hào)的指數(shù)分布為

(5.5.64)式中加權(quán)值

(5.5.65)上面兩式表明交叉項(xiàng)的幅度由WEIGHT控制,當(dāng)Ω=(Ω1+Ω2)/2時(shí),交叉項(xiàng)的幅度與成反比,另外,它以距離［Ω-(Ω1+Ω2)/2］2式指數(shù)地衰減。

2.廣義指數(shù)分布（GED）與巴特沃斯分布（BUD）

1)廣義指數(shù)分布廣義指數(shù)分布的核函數(shù)為

（5.5.66）

式中,N、M為正的冪次;θ1和τ1分別是正的頻率和時(shí)間尺度常數(shù)。下面說(shuō)明該核函數(shù)具有低通特性。

先分析下面低通濾波函數(shù)θM(x)的特性:（5.5.67）

該函數(shù)的濾波特性如圖5.5.5所示。此圖表明,當(dāng)x1一定時(shí),隨M的增大,通帶變得愈平坦,過(guò)渡帶愈窄。當(dāng)M→∞時(shí),變成在x=x1

處截止的理想低通濾波器。

圖

5.5.5θM(x)的濾波特性

令 ,則

(5.5.68)這種情況下,廣義指數(shù)核函數(shù)變成一維低通濾波函數(shù)θM(x),當(dāng)N=M=1時(shí), ,x=τθ,上式變?yōu)?5.5.69)即θ1(x)是指數(shù)分布的核函數(shù)。此時(shí),核函數(shù)具有最差的低通特性。當(dāng)M加大時(shí)，濾波性能變好。因此說(shuō),廣義指數(shù)分布是可以改善時(shí)頻分布性能的一種分布，通過(guò)改善,使其自項(xiàng)失真小，

交叉項(xiàng)衰減大。

2）巴特沃斯分布巴特沃斯分布的核函數(shù)定義為

N,M,θ1,τ1>0(5.5.70)上式表明,該核函數(shù)具有可變的、平坦的通帶及窄的過(guò)渡帶的低通特性。當(dāng) ,

上式變成

(5.5.71)上式即是濾波器中的最大平坦平方幅度的巴特沃斯濾波器的幅度平方函數(shù),當(dāng)M加大時(shí)，它趨近于理想低通濾波器。因此，這種分布也可以改善指數(shù)分布的時(shí)頻分布的性能。除了上面介紹的廣義雙線(xiàn)性時(shí)頻分布以外，常用的還有平滑偽維格納分布（SPWD）、錐形（Cone）核分布（CKD）、降低干擾項(xiàng)分布（RID）、貝塞爾（Bessel）分布（BD）等等。

5.6時(shí)頻分析在編隊(duì)目標(biāo)架次檢測(cè)中的應(yīng)用

5.6.1編隊(duì)飛機(jī)目標(biāo)多普勒特性分析本節(jié)我們分析編隊(duì)飛機(jī)目標(biāo)的多普勒頻率特性，從而探討在多普勒域?qū)崿F(xiàn)多目標(biāo)分辨的可能性。眾所周知，雷達(dá)目標(biāo)回波的多普勒頻率fd取決于雷達(dá)工作波長(zhǎng)λ、目標(biāo)運(yùn)動(dòng)速度v以及目標(biāo)飛行方向與雷達(dá)視線(xiàn)的夾角φ，

即

(5.6.1)

圖5.6.1顯示了兩架飛機(jī)編隊(duì)情況。假設(shè)編隊(duì)飛行時(shí)，兩架飛機(jī)的速度相同，則它們回波的多普勒頻率的差別就取決于它們各自與雷達(dá)視線(xiàn)夾角的不同，

即

(5.6.2)φ1與φ2的差別是由兩飛機(jī)的間距d和目標(biāo)距雷達(dá)的距離R所決定的。一般有R>>d,

因此，φ1與φ2的差別是很小的，進(jìn)一步化簡(jiǎn)，

有

(5.6.3)圖5.6.1兩架飛機(jī)編隊(duì)情況式中，利用了當(dāng)R>>d

時(shí)，sin(φ1-φ2)≈φ1-φ2≈(d/R)sinφ,其中φ=(φ1+φ2)/2是平均的雷達(dá)視線(xiàn)與飛行方向的夾角。由(5.6.3)式可以看出，