重新分佈法

重新分佈法等別稱。以時頻譜或短時距傅立葉變換而言，重新分佈法可藉由估算區域性的瞬時頻率以及群延遲，使模糊的時頻資料點重新定位並清晰化。當訊號可藉由分析視窗進行時域和頻域的分離時，這項時頻座標的重新對映是相當精準的。概述.許多訊號的能量分佈都具有隨著時間和頻率變化的特性。舉例而言，任何有起始或者結束的音訊訊號之能量分佈都將隨著時間改變，而大多數聲音的能量不論在時間上還是頻率上都有相當大的變化；時頻分析便經常用於尋找此類訊號的特徵，這種分析法將原為一維的時域訊號轉換為二維的時域–頻域函式。而經分析轉換後的結果便稱為，它可以用於描述一個訊號隨時間的頻譜能量分佈變化，就像樂譜描述樂曲中高低音隨時間的變化一樣。在音訊訊號分析中，時頻譜是最常用的法；這可能是因為該方法長久以來為學界所熟知，且能夠避免產生不易解讀的「交叉項」。不過，計算時頻譜所須使用的視窗函式使得分析結果必須在時域或頻域的解析度中做一取捨，也就是時頻譜在時域及頻域的解析度必定至少有其中一者不理想。這項時頻表示的重新分佈法便是透過先把資料對映至較接近原始訊號之真實支撐區域的座標，以將時頻譜或其它時頻表示下較模糊的資料點重新定位並聚焦。時頻譜.時頻譜為最為人所熟知的時頻表示法之一，其定義為短時距傅立葉變換值的平方：formula_1雖然短時距的相位頻譜確實含有與訊號相關的重要時態資訊，但這些資訊通常難以解讀，故一般而言在短時距頻譜分析中只考慮其短時距強度頻譜。以時頻表示法而言，時頻譜具有相對較差的解析度。它在時域與頻域的解析度取決於分析時所選擇的視窗函式，且其中一者愈清晰將導致另一者愈模糊。相對於時頻譜，維格納–魏爾分佈則是另一種解析度較高的時頻表示法；本質上，維格納–魏爾分佈的做法與短時距傅立葉變換類似，惟其選擇的視窗函式與訊號之間有完好的對應。維格納–魏爾分佈在時域和頻域中皆可得到高解析度，但它同時也具有高度的非線性和非區域性性等特質，因此對雜訊相當敏感且容易產生覆蓋訊號本身的交叉項，使得含有多成分的訊號分析後的資訊難以解讀。利用雙線性轉換得到時頻表示的科恩系列分佈函式則可視為「平滑化」的維格納–魏爾分佈。它藉由一平滑核降低分佈對於雜訊的敏感度，並抑制交叉項的產生，但此舉將使時域和頻域的分佈模糊化；這不僅會降低頻譜的解析度，同時也可能使原始維格納準機率分佈中實際上沒有能量的區域得到非零的能量值。時頻譜便屬於科恩系列分佈之一；它也可視為一種平滑化的維格納–魏爾分佈，惟其使用之平滑核即相等於其分析視窗函式的維格納–魏爾分佈。重新分佈法便是先將維格納–魏爾分佈平滑化，再將它於原始訊號成分之真實支撐區域附近的座標點重新聚焦，進而達到銳化的目的。實驗結果顯示，此方法可以減少任一科恩系列分佈在時域與頻域的模糊現象。對時頻譜而言，重新分佈法便是利用訊號的短時距相位頻譜來修正其時頻座標點，再將其重新映射回較靠近所分析訊號之真實支撐區域的位置。方法.重新分佈法首先由小寺邦彥的時頻解析度。標準的滑動視窗法乃是將一時域訊號formula_2分解為一組基底訊號formula_3的加權總和，其中代表每項基底的係數則為formula_4。基底訊號的定義如下：formula_5其中，formula_6為一實值之低通核函式，如短時距傅立葉變換中的視窗函式。此般分解下的各項係數則可定義如下：formula_7其中，formula_8為訊號formula_2在時間位移為formula_10、使用formula_6做為視窗函式下所得到的傅立葉變換；而formula_12與formula_13則分別代表formula_8的強度以及相位值。formula_2可透過滑動視窗的係數以下式重建：formula_16當一個訊號的強度頻譜為formula_17，且其時間上較相位上的變化來得慢時，上述重建式中積分值最主要的貢獻滿足下列穩相條件的formula_18及其鄰近區域：formula_19亦即滿足下列定義的formula_20及其鄰近區域：formula_21此現象在光學等領域中稱為：對週期或準週期函式而言，其傅立葉相位頻譜在其振盪頻率附近相對於時間的變化較慢，反之在外圍區域則較快。以時域上能量較為集中的脈衝訊號為例，其相位頻譜在脈衝時間點附近的變化較慢，而較遠的區域則變化較快。在訊號重建時，合成波形之正值與負值於相位快速變化的頻率區域會因破壞性干涉而互相抵消。只有相位變化緩慢即發生在相位隨時間與頻率變化最緩慢的地方。由此知，可以直接從特徵。修正滑動視窗法便是將計算formula_24的時頻點從原始的formula_18重新移至區域性極大貢獻點formula_26。由於本質上這個區域性極大點和質量分佈中心點有相似之處，故也常類比成「重心」。該類比也有助於做出下述的理解：一個訊號的頻譜能量分佈對重心的貢獻只在「有能量」的區域才具有意義，因此對時頻譜上量值為零的座標點使用重新分佈法是沒有意義的。快速時頻重新分佈法.在數位訊號處理中，最常見的做法是先在時域與頻域分別對訊號進行取樣，並使用離散傅立葉變換計算時域訊號formula_27的頻域變換formula_28。因離散的時頻資料點無法直接計算偏微分，小寺邦彥等人所提出的重新分佈法並不能直接套用於離散形式的短時距傅立葉變換；甚至有人指出，這項難處就是使得重新分佈法無法得到更廣泛應用的主因。然而，偏微分的值可以利用有限差分法來估計。舉例而言，在相位頻譜上兩個鄰近的時間點：formula_32在頻譜上，某些區域的相位變化主要源自於附近某單一弦波成分的旋轉；這些區域內的相位值為一線性函式，因此對於足夠微小的formula_31及formula_34而言，上式的有限差分可為相位值的偏微分做一良好的近似。道格拉斯·尼爾森所計算之偏微分近似值和前述的有限差分法無異。弗朗索瓦·俄歇則指出，小寺邦彥等人所提出的時頻譜重新分佈法可透過下述的一般化，延用於任一科恩系列的時頻表示法：formula_35其中，formula_36為formula_2的維格納–魏爾分佈，而formula_38則是定義該分佈的核函式。他們更提出了一項不須直接計算相位之偏微分，便可快速且準確地得到重新分佈後之時頻譜的方法。以時頻譜而言，其重新分佈法的計算如下：formula_39其中，formula_40是訊號以視窗函式formula_6所得到的短時距傅立葉變換、formula_42是訊號以時間加權之視窗函式formula_43所得到的短時距傅立葉變換、formula_44則是訊號以對時間取導數之視窗函式formula_45所得到的短時距傅立葉變換；故知，利用輔助視窗函式formula_46以及formula_47，可得到含原始在內共三種傅立葉變換的結果。接著，將這三種變換分別在formula_18處取值並做代數上的組合，便可在任一時頻座標點formula_18計算重新分佈法的結果。由於上述的方法只是在短時距頻譜上的某單一時頻點取值，而不須做任何導數的計算，因此可視作離散短時距傅立葉變換之重新分佈的快速演演算法。上述方法有一限制，即formula_50之值必須為非零；不過，這項限制一般而言皆會被滿足，因為只在頻譜有能量的區域才有值可以「重新分佈」，而方法本身在頻譜值為零處也沒有意義。可分離性.短時距傅立葉變換經常用於估算一個多成分訊號中個別單一成分的振幅和相位，而時頻重新分佈則能將變換後的時頻點集中至能量分佈的區域性重心，使頻譜更加清晰。對單一成分的訊號而言，其瞬時頻率可由該訊號經任意短時距傅立葉變換後得到的相位之偏微分估算。若訊號可拆解為多個成分：formula_51且定義各成分的瞬時頻率為其相位對時間的偏微分，即：formula_52則可以藉由設計適當的帶通濾波器使其通帶只包含某一特定訊號成分，再從該濾波器經變換後之頻譜的相位值，計算該訊號成分的瞬時頻率。尼爾森稱這項頻域中的條件為「可分離性」，且所有透過此方法分析的訊號皆須滿足該條件。若該條件並未滿足，則因不能透過短時距傅立葉變換估算個別成分的引數，而無法得到欲求的多成分訊號拆解；在這樣的情況下，必須採用其它的分析視窗函式，以滿足可分離性的條件。若某特定的短時距頻譜分析視窗可用於拆解訊號中不同頻率的成分，則每個濾波器經短時距傅立葉變換後的脈衝響應便是經濾波後的結果，故formula_53之相位值對時間的偏微分值，即相等於頻率formula_54下主要成分之相位值對時間的偏微分值。因此，若某瞬時頻率為formula_55的訊號成分formula_56 在formula_54附近為主要成分時，則該成分的瞬時頻率可透過原始訊號經短時距傅立葉變換後位於formula_54的相位值求得，即：formula_59頻域的可分離性要求每個帶通濾波器經短時距傅立葉變換後，最多隻可透過一項訊號成分，而在時域同樣也有類似的限制：訊號中兩個不同時間的成分彼此間分隔須夠大，使得任一時域視窗下的分段，都不得包含超過一個事件。此即時域下的可分離性條件，且等義於任兩事件間的間隔須大於短時距傅立葉變換濾波器之脈衝響應長度，也就是formula_6中非零取樣點所涵蓋的範圍。一般而言，任一多成分訊號有無限多種相同合理的拆解方式；不過，可分離性必須是在所欲求得的拆解下進行考慮。舉例而言，在分析一個語音訊號時，若使用時間範圍較聲門脈衝間隔長的分析視窗，則足以將不同諧波分離；但個別的聲門脈衝則會因單一視窗包含了多個脈衝而被模糊，也就是所選的分析視窗使得訊號在時域中不具可分離性。反之，若使用時間範圍遠小於聲門脈衝間隔的分析視窗，則能夠成功分離每個脈衝，因為每個視窗中最多隻包含一個脈衝；但不同諧波便會彼此模糊交疊，因為視窗在頻域中的範圍比單一諧波還要寬，也就是所選的分析視窗使得訊號在頻域中不具可分離性。