時頻分析中的核方法

時頻分析的性質由其分布的核所決定，藉由檢視分布核的限制條件我們能很容易的了解時頻分布的優缺點並讓我們快速的選擇符合需求的時頻分布。因此我們將利用時頻分布和核函數之間的關聯研究相對應的時頻分析性質，並以此作為設計時頻分析的基準。 所有的時頻分布皆可以表示為 C ( t , ω ) = 1 4 π 2 ∫ ∫ ∫ s ∗ ( u − 1 2 τ ) s ( u + 1 2 τ ) ϕ ( θ , τ ) e − j θ t − j τ ω + j θ u d u d τ d θ {\displaystyle C(t,\omega )={\frac {1}{4\pi ^{2}}}\int \int \int s^{*}(u-{\frac {1}{2}}\tau )s(u+{\frac {1}{2}}\tau )\phi (\theta \,,\tau )e^{-j\theta t-j\tau \omega +j\theta u}dud\tau d\theta } ，其中 ϕ ( θ , τ ) {\displaystyle \phi (\theta \,,\tau )} 被稱作時頻分布的核。 下表列出了幾個常見的時頻分布與其對應的核 我們可以將廣義的時頻分布形式做變形以便從不同的角度觀察時頻分布與其核的關係，在時頻分布的性質與核中會看到這些表示式提供對核的刻畫與限制，其中常用的幾種等價表示式為 若令廣義的模稜函數 M ( θ , τ ) {\displaystyle M(\theta \,,\tau )} 為 M ( θ , τ ) = ∫ ∫ s ∗ ( u − 1 2 τ ) s ( u + 1 2 τ ) ϕ ( θ , τ ) e j θ u d u = ϕ ( θ , τ ) A ( θ , τ ) {\displaystyle M(\theta \,,\tau )=\int \int s^{*}(u-{\frac {1}{2}}\tau )s(u+{\frac {1}{2}}\tau )\phi (\theta \,,\tau )e^{j\theta u}du=\phi (\theta \,,\tau )A(\theta \,,\tau )} 其中 A ( θ , τ ) {\displaystyle A(\theta \,,\tau )} 為信號 s ( t ) {\displaystyle s(t)} 本身的模稜函數,則時頻分布可表為 M ( θ , τ ) {\displaystyle M(\theta \,,\tau )} 的特徵函數 C ( t , ω ) = 1 4 π 2 ∫ ∫ M ( θ , τ ) e − j θ t − j τ ω d θ d τ {\displaystyle C(t,\omega )={\frac {1}{4\pi ^{2}}}\int \int M(\theta \,,\tau )e^{-j\theta t-j\tau \omega }d\theta d\tau } 計算廣義的局部自相關函數 R t ( τ ) {\displaystyle R_{t}(\tau )} 為 R t ( τ ) = 1 2 π ∫ ∫ s ∗ ( u − 1 2 τ ) s ( u + 1 2 τ ) ϕ ( θ , τ ) e j θ ( u − t ) d θ d u {\displaystyle R_{t}(\tau )={\frac {1}{2\pi }}\int \int s^{*}(u-{\frac {1}{2}}\tau )s(u+{\frac {1}{2}}\tau )\phi (\theta \,,\tau )e^{j\theta (u-t)}d\theta du} 則時頻分布可表示成 R t ( τ ) {\displaystyle R_{t}(\tau )} 的傅立葉轉換，可類比功率譜密度和自相關函數之間的關係。 C ( t , ω ) = 1 2 π ∫ R t ( τ ) e − j ω τ d τ {\displaystyle C(t,\omega )={\frac {1}{2\pi }}\int R_{t}(\tau )e^{-j\omega \tau }d\tau } 將核函數對參數 θ {\displaystyle \theta } 做傅立葉變換可得 r ( t , τ ) = 1 2 π ∫ ϕ ( θ , τ ) e − j t θ d θ {\displaystyle r(t,\tau )={\frac {1}{2\pi }}\int \phi (\theta ,\tau )e^{-jt\theta }d\theta } C ( t , ω ) = 1 2 π ∫ ∫ r ( t − u , τ ) s ∗ ( u − 1 2 τ ) s ( u + 1 2 τ ) e − j ω τ d τ d u {\displaystyle C(t,\omega )={\frac {1}{2\pi }}\int \int r(t-u,\tau )s^{*}(u-{\frac {1}{2}}\tau )s(u+{\frac {1}{2}}\tau )e^{-j\omega \tau }d\tau du} 對於廣義時頻分布做變數變換 x ′ = u − 1 2 τ ; x = u + 1 2 τ {\displaystyle x'=u-{\frac {1}{2}}\tau ;x=u+{\frac {1}{2}}\tau } 則有 C ( t , ω ) = ∫ ∫ K ( t , ω , x ′ , x ) s ∗ ( x ′ ) s ( x ) d x ′ d x {\displaystyle C(t,\omega )=\int \int K(t,\omega ,x',x)s^{*}(x')s(x)dx'dx} K ( t , ω , x ′ , x ) = 1 2 π r ( t − 1 2 ( x + x ′ ) , x ′ − x ) e − j ω ( x ′ − x ) {\displaystyle K(t,\omega ,x',x)={\frac {1}{2\pi }}r(t-{\frac {1}{2}}(x+x'),x'-x)e^{-j\omega (x'-x)}} 由於任意的二維函數都可以作為核而產生時頻分布，事實上我們要面對的是無限多種的核(與時頻分布)，因此在現實應用與討論中我們常常將注意力放在一部份擁有特定足夠好性質的核函數，以下簡介之 核函數僅依賴的參數 θ τ {\displaystyle \theta \,\tau } 的積 ϕ ( θ τ ) = ϕ ( θ τ ) {\displaystyle \phi (\theta \,\tau )=\phi (\theta \tau )} 核函數可被分解為兩個單變數函數的積 ϕ ( θ τ ) = ϕ 1 ( θ ) ϕ 2 ( τ ) {\displaystyle \phi (\theta \,\tau )=\phi _{1}(\theta )\phi _{2}(\tau )} 核函數不依賴時間參數 t {\displaystyle t} 與頻率參數 ω {\displaystyle \omega } ，函數 K ( t , ω , x ′ , x ) = 1 2 π r ( t − 1 2 ( x + x ′ ) , x ′ − x ) e − j ω ( x ′ − x ) {\displaystyle K(t,\omega ,x',x)={\frac {1}{2\pi }}r(t-{\frac {1}{2}}(x+x'),x'-x)e^{-j\omega (x'-x)}} 對參數 x , x ′ {\displaystyle x,x'} 為雙線性 以下我們討論時頻分布的基本性質與為了達到這些性質核函數必須要滿足的條件 考慮時頻分布函數 C ( t , ω ) {\displaystyle C(t,\omega )} ，我們希望對於頻率參數 ω {\displaystyle \omega } 積分之後能得到在時間 t {\displaystyle t} 的信號能量 ,因此我們必須使下式成立 ∫ C ( t , ω ) d ω = 1 2 π ∫ ∫ ϕ ( θ , 0 ) | s ( u ) | 2 e j θ ( u − t ) = | s ( u ) | 2 {\displaystyle \int C(t,\omega )d\omega ={\frac {1}{2\pi }}\int \int \phi (\theta ,0)|s(u)|^{2}e^{j\theta (u-t)}=|s(u)|^{2}} 因此我們必須要求 ϕ ( θ , 0 ) = 1 {\displaystyle \phi (\theta ,0)=1} 。相同的為了得到 ∫ C ( t , ω ) d t = | s ( ω ) | 2 {\displaystyle \int C(t,\omega )dt=|s(\omega )|^{2}} , 我們必須要求 ϕ ( 0 , τ ) = 1 {\displaystyle \phi (0,\tau )=1} 最後，考慮總能量我們則須要求 ϕ ( 0 , 0 ) = 1 {\displaystyle \phi (0,0)=1} 為了使時頻分布的計算結果為實數，考慮時頻分布的特徵函數計算式，相當於對 M ( θ , τ ) {\displaystyle M(\theta \,,\tau )} 進行兩次傅立葉轉換 C ( t , ω ) = 1 4 π 2 ∫ ∫ M ( θ , τ ) e − j θ t − j τ ω d θ d τ {\displaystyle C(t,\omega )={\frac {1}{4\pi ^{2}}}\int \int M(\theta \,,\tau )e^{-j\theta t-j\tau \omega }d\theta d\tau } 因此廣義模稜函數 M ( θ , τ ) {\displaystyle M(\theta \,,\tau )} 必須滿足共軛對稱性 M ( θ , τ ) = M ∗ ( − θ , − τ ) {\displaystyle M(\theta \,,\tau )=M^{*}(-\theta \,,-\tau )} 在信號本身為實數信號時的 s ( t ) {\displaystyle s(t)} 模稜函數 A ( θ , τ ) {\displaystyle A(\theta \,,\tau )} 直接滿足共軛對稱，因此條件變為 ϕ ( θ , τ ) = ϕ ∗ ( − θ , − τ ) {\displaystyle \phi (\theta \,,\tau )=\phi ^{*}(-\theta \,,-\tau )} 考慮對輸入信號做時頻位移 s ′ = e j ω 0 t s ( t − t 0 ) {\displaystyle s'=e^{j\omega _{0}t}s(t-t_{0})} ，若要對應的時頻分布滿足 C ′ ( t , ω ) = C ( t − t 0 , ω − ω 0 ) {\displaystyle C'(t,\omega )=C(t-t_{0},\omega -\omega _{0})} ，則核函數必須不依賴時間和頻率。 考慮縮放後的信號 s ′ ( t ) = ( a ) s ( a t ) {\displaystyle s'(t)={\sqrt {(}}a)s(at)} , 我們希望對應的時頻分布滿足 C ′ ( t , ω ) = C ( a t , ω a ) {\displaystyle C'(t,\omega )=C(at,{\frac {\omega }{a}})} ，則核函數必須是乘積核 ϕ ( θ τ ) = ϕ ( θ τ ) {\displaystyle \phi (\theta \,\tau )=\phi (\theta \tau )} 我們希望能從時頻分布 C ( t , ω ) {\displaystyle C(t,\omega )} 中重新恢復輸入信號 s ( t ) {\displaystyle s(t)} 由關係式 A ( θ , τ ) = M ( θ , τ ) ϕ ( θ , τ ) {\displaystyle A(\theta \,,\tau )={\frac {M(\theta \,,\tau )}{\phi (\theta \,,\tau )}}} 兩側同時取傅立葉變換可得 s ∗ ( u − 1 2 τ ) s ( u + 1 2 τ ) = 1 2 π ∫ M ( θ , τ ) ϕ ( θ , τ ) e − j θ u d θ {\displaystyle s^{*}(u-{\frac {1}{2}}\tau )s(u+{\frac {1}{2}}\tau )={\frac {1}{2\pi }}\int {\frac {M(\theta \,,\tau )}{\phi (\theta \,,\tau )}}e^{-j\theta u}d\theta } 做變數變換 t = u + 1 2 τ , t ′ = u − 1 2 τ {\displaystyle t=u+{\frac {1}{2}}\tau ,t'=u-{\frac {1}{2}}\tau } ，並取 t ′ = 0 {\displaystyle t'=0} 可得 s ( t ) = 1 2 π s ∗ ( 0 ) ∫ M ( θ , τ ) ϕ ( θ , τ ) e − j θ u d θ = 1 2 π s ∗ ( 0 ) ∫ C ( t ′ , ω ) ϕ ( θ , τ ) e j t ω + j θ ( t ′ − t 2 ) d t ′ d ω d θ {\displaystyle s(t)={\frac {1}{2\pi s^{*}(0)}}\int {\frac {M(\theta \,,\tau )}{\phi (\theta \,,\tau )}}e^{-j\theta u}d\theta ={\frac {1}{2\pi s^{*}(0)}}\int {\frac {C(t',\omega )}{\phi (\theta \,,\tau )}}e^{jt\omega +j\theta (t'-{\frac {t}{2}})}dt'd\omega d\theta } 現在我們考慮對同一個信號使用不同的核函數 ϕ 1 ( θ , τ ) , ϕ 2 ( θ , τ ) {\displaystyle \phi _{1}(\theta \,,\tau ),\phi _{2}(\theta \,,\tau )} ，獲得兩個時頻分布 C 1 , C 2 {\displaystyle C_{1},C_{2}} 這兩個時頻分布的特徵函數為 M i ( θ , τ ) = ϕ i ( θ , τ ) ∫ s ∗ ( u − 1 2 τ ) s ( u + 1 2 τ ) d u , i = 1 , 2 {\displaystyle M_{i}(\theta \,,\tau )=\phi _{i}(\theta \,,\tau )\int s^{*}(u-{\frac {1}{2}}\tau )s(u+{\frac {1}{2}}\tau )du,i=1,2} 因此特徵函數和核函數的關聯為 M 1 ( θ , τ ) = ϕ 1 ( θ , τ ) ϕ 2 ( θ , τ ) M 2 ( θ , τ ) {\displaystyle M_{1}(\theta \,,\tau )={\frac {\phi _{1}(\theta \,,\tau )}{\phi _{2}(\theta \,,\tau )}}M_{2}(\theta \,,\tau )} 因此我們可以將 C 1 {\displaystyle C_{1}} 以 C 2 {\displaystyle C_{2}} 表達為: C 1 ( t , ω ) = 1 4 π 2 ∫ ∫ ∫ ∫ ϕ 1 ( θ , τ ) ϕ 2 ( θ , τ ) C 2 ( t ′ , ω ′ ) e j θ ( t ′ − t ) + j τ ( ω − ω ′ ) {\displaystyle C_{1}(t,\omega )={\frac {1}{4\pi ^{2}}}\int \int \int \int {\frac {\phi _{1}(\theta \,,\tau )}{\phi _{2}(\theta \,,\tau )}}C_{2}(t',\omega ')e^{j\theta (t'-t)+j\tau (\omega -\omega ')}} Leon Cohen, "Generalized phase-space distribution functions," Jour. Math. Phys., vol. 7, pp. 781-786, 1966. Leon Cohen, "Time-frequency analysis," 1995.