最小平方頻譜分析法

最小平方頻譜分析法是一種利用最小平方法尋找適配於資料點之最佳正弦曲線，以估算頻譜的方法。其數學原理與科學界中最常用的傅立葉分析相似。一般而言，傅立葉分析會將間隔較長之訊號的長週期雜訊放大，而最小平方頻譜分析法則解決了這個問題。最小平方頻譜分析法也稱為凡尼切克法以及等人也曾開發出與之關係密切的其它方法。歷史背景.長久以來，科學界對於傅立葉分析、週期圖法與正弦曲線之最小平方擬合間的密切關聯早有所知。然而，大多數以上述理論為基礎開發的方法僅適用於取樣間距相等的訊號。1963年，荷蘭數學和電腦科學研究學會的弗裡克·巴寧以避免過適、或使用正交匹配追蹤法。新不倫瑞克大學的加拿大大地測量學家於1969年也提出了匹配追蹤法，他稱之為「連續頻譜分析」，並稱其所得的結果為「最小平方週期圖」，且同時適用於取樣間距相同或不同的訊號。他更於1971年將該方法一般化，使其能夠分析出訊號中單一平均值以外的各種系統成分，包含可預測但大小未知的線性、二次方或指數趨勢。雪梨大學的尼可拉斯·隆布於1976年簡化了凡尼切克法，並指出該方法與週期圖法有相當密切的關連。後來，傑佛瑞·史卡構也定義了非等距取樣訊號的週期圖，並對其做了分析；他並發現，只要將該週期圖稍作修改，便能夠得到與隆布的最小平方擬合法完全相同的結果。史卡構曾表示他的論文「並未提出一個新的測量方法，而是對現有最常用的方法——週期圖法——在取樣間距不相等時的可靠性以及效率進行研究」。他在論文中同時引用了最小平方弦波擬合以及週期圖法分析，並指出該兩種方法是完全相等的。媒體對相關研究的發展則有如下的摘要：1989年，皇后大學的麥可·科恩伯格開發了所謂的「快速正交搜尋法」；該方法藉由將目標轉化為線性規劃問題，並使各項係數的L1-距離最小化，以更快速地得到解答。方法.凡尼切克法.凡尼切克法乃是利用標準的線性回歸或最小平方法，將一個離散訊號近似為許多不同頻率的弦波之加權總和。這些弦波的頻率則是以類似巴寧所提出的方法來決定，並另外藉由減少最小平方擬合後的殘差以使每次決定的新頻率得到最佳化。假設formula_1為一個由許多不同頻率的基底弦波加權總和所得的訊號向量，其基底函式以矩陣formula_2表示，而每個基底所對應的權重則以formula_3向量表示，也就是：formula_4會希望權重向量formula_3使得上式在近似formula_1時的平方誤差和為最小。利用標準的線性回歸法即可求得formula_3的閉合形式解為：formula_8其中，矩陣formula_2可由任意彼此獨立的基底函式所組成；對頻譜分析而言，這些基底函式一般為某特定頻率範圍內均勻分佈的正弦或餘弦波。如果在某過窄的頻寬內選擇了過多的頻率，則這些所選的函式可能不夠獨立、矩陣formula_2本身較為病態，而使得最後求得的頻譜不具意義。當矩陣formula_2的所有基底函式皆彼此正交，而這情況同時也是離散傅立葉變換的特例之一。formula_18     隆布所提出上述的簡化式可用於大多數的情況，因為一般而言不同弦波彼此之間的相關係數較小，惟同一頻率下成對的正弦和餘弦波彼此間並非獨立而不適用。本質上，這和傳統週期圖法的公式無異，但可用於非等距取樣的訊號。其中，formula_3向量可為訊號的潛在頻譜做有效的估計，但由於基底之間的相關性被忽略，formula_21便不再能夠良好地近似原始的訊號，這個方法也因此不再是真正的「最小平方法」，惟學界仍繼續如是稱之。隆布—史卡構法.在標準的週期圖法公式中，一般是直接計算訊號與正弦或餘弦波之間的內積。史卡構則對其稍作修改：首先，先找到一時間延遲項formula_22使得成對的正弦與餘弦波在取樣時間為formula_23時將會彼此正交；同時，他也對可能具有不同能量的兩個基底之間做了標準化，以取得在某頻率下更加準確的能量估計值。這些修改使得史卡構的週期圖法和隆布的最小平方法即完全相同。時間延遲項formula_22的定義如下：formula_25頻率formula_26的週期圖則可估算如下：formula_27史卡構並表示，上式所得的週期圖與等距取樣情況下的結果具有相同的統計分佈。在任一單獨的頻率formula_26下，這個方法求得的頻譜能量與最小平方法使用該頻率之弦波擬合後所得的能量相同。擬合後函式的形式如下：formula_29一般化.標準的隆布—史卡構週期圖法適用於平均值為零的模型，因此一般而言在計算週期圖之前會先減去訊號的平均值並假設其為零。然而，當模型的平均值為非零時，這項假設並不準確。一般化的隆布—史卡構週期圖法便將此假設移除，並在計算頻譜時一併求得其平均值。如此一來，擬合後函式的形式如下：formula_30此一般化的隆布—史卡構週期圖法也被稱為「浮動平均週期圖法」。快速正交搜尋法.科恩伯格於1989年提出瞭如下的方法：先從一的函式集合中選出一稀疏子集等其它的問題。基底追蹤法.陳與多諾霍提出的法用於擬合，惟此法定義能夠使各項係數的L1-距離得到最小化者即為最佳解。如此一來，頻譜分析便可藉此轉化為已知的線性規劃問題，並可利用現有的快速演演算法得到解答。卡方法.大衛·帕爾默提出的卡方法則能夠找到任意數量之諧波的最佳擬合函式，提高了尋找非弦波諧和函式的自由度。這是一項基於快速傅立葉變換、對取樣間距任意且標準誤差不均的訊號進行分析的快速演演算法。實作此方法的原始碼也已被公開。由於離散訊號往往不具有相等的取樣間距，此方法先將訊號「柵格化」，即稀疏地在數個取樣點間填入一串時間序列。所有時間軸上互相交會的柵格點之統計權重都將被設為零，亦即取樣點之間的為無窮大。應用.最小平方頻譜分析法最主要的用處是對未完整記錄的訊號進行頻譜分析，且不須對其竄改或摻入不存在的資料。利用最小平方頻譜分析法所得的頻譜中之各強度量值即代表了某個頻率或週期對於整個時間序列的變異數之貢獻。一般而言，經由上述方法所定義的頻譜強度可用於推算輸出結果的顯著水準狀態。此外，凡尼切克法求得的頻譜強度也可由分貝值表示。值得注意的是，凡尼切克頻譜在統計上具有Β分佈。凡尼切克的最小平方頻譜之反運算可由下述方法實現：若將正運算的過程表示為一矩陣的乘法，那麼求取它的反矩陣或偽反矩陣便可進行反運算；當所選的基底弦波在取樣點時皆彼此獨立且其總數和取樣點數量相等時，反運算所得的結果便會與原始的訊號相同。相對地，週期圖法並不具有已知的反運算方法。實作.最小平方頻譜分析法可用一頁以內的MATLAB程式碼實作完成。簡言之：「欲求得最小平方頻譜，我們必須計算formula_31個頻譜值……可以藉由formula_31次的最小平方法近似，每次求得一種不同頻率。」也就是說，對一所欲擬合之頻率集合中的每個頻率formula_33，將其對應的正弦於訊號的取樣時間點formula_23取值，計算訊號與弦波之取樣值的內積並給予適當的標準化；若是使用隆布—史卡構的週期圖法，則在取內積之前先計算每個頻率formula_33的時間延遲項formula_38使得該頻率下的正弦與餘弦成分能夠彼此正交；最後，從振幅的內積值便可求得該兩項弦波成分的頻譜能量。若將這項程式套用於取樣間距相等、且擬合弦波頻率選為取樣時間全長倒數之整數倍的訊號時，即相當於離散傅立葉變換。上述即凡尼切克的原始方法，它在處理每個弦波成分時都將其視為彼此獨立，即使在取樣點時這些弦波有可能並不彼此正交。相對地，透過矩陣方程式求解、並將原始訊號的總變異數分割至指定的弦波頻率之間的方法，也能夠在一次考慮所有弦波成分的情況下求得頻譜值。這種矩陣求解的最小平方法可直接透過MATLAB內建的反斜線運運算元得到。麥克·克雷默，並指出：相對地，隆布提出的方法便如同標準的週期圖法，可使用任意多數或任意高密度的頻率成分做為擬合用的基底；也就是說，這個方法可在頻域下以任意的倍率進行過取樣。在傅立葉分析中，用於擬合訊號的弦波皆為彼此正交，所以前述的獨立投影內積法和一次計算的矩陣求解最小平方法並沒有差異；也就是說，使用最小平方法將變異數分割至頻率不同且彼此正交的弦波時即不須計算反矩陣。由於此類分析可用快速傅立葉變換的方法實作，當訊號完整且取樣間距相等時，便不會使用相對耗時的最小平方頻譜分析法。