正弦曲線擬合若干問題探討

1、，第 29 卷 第 14 期計算機工程與設(shè)計Computer Engineering and Design2008 年 7 月July 2008Vol . 2 9No . 1 4正弦曲線擬合若干問題探討齊國清 1， 呂健 2(1. 大連海事大學(xué) 信息工程學(xué)院，遼寧 大連 116024；2. 大連海事大學(xué) 航海學(xué)院，遼寧 大連 116024)摘 要：研究了測量噪聲較小情況下正弦曲線的最小二乘多項式擬合誤差與擬合階數(shù)的關(guān)系，分別采用均方誤差和誤差平 方和分析了測量噪聲以及測量數(shù)據(jù)有效位數(shù)對擬合誤差的影響，對多周期正弦曲線擬合以及正弦曲線的外推存在的問題進 行了探討，指出了正弦曲線的最小二乘多項式擬

2、合方法的局限性。最后，提出了一種基于傅利葉變換的頻率已知正弦曲線 擬合方法, 仿真結(jié)果表明其性能優(yōu)于最小二乘多項式擬合方法。關(guān)鍵詞：曲線擬合; 最小二乘擬合; 正弦曲線; 誤差分析; 傅利葉變換中圖法分類號：TP311.11文獻標(biāo)識碼：A文章編號：1000-7024 (2008) 14-3677-04Investigation of sine wave fitting algorithms. .QI Guo-qing1,LU Jian2(1. College of Information Engineering, Dalian Maritime University, Dalian 11602

3、4, China;2. College of Navigation, Dalian Maritime University, Dalian 116024, China)Abstract：The relation between the order of the polynomials and fitting error of LS sine-wave fitting algorithm is discussed. The effects of measurement noise and the significant digits of the measured data on fitting

4、 error are also studied in terms of mean square error and sum of squared error. Problems encountered in fitting multi-cycle sine wave and extrapolation of sine-wave are investigated. The limitations of applying LS fitting to sine-wave fitting are pointed out. Finally, a simple fitting algorithm base

5、d on Fourier transform for sine wave withknown frequency is proposed. Simulation results show that the performance of the proposed method is better than LS polynomial fitting.Key words：curve fitting; LS fitting; sine wave; error analysis;Fourier transform多于所需要的最小數(shù)目。前者只是一般建立和求解方程的過0引言程，得到解析表達式 =的目的僅僅

6、是為了得到不在采樣根據(jù)有限的離散測量點進行曲線擬合是工程實踐中經(jīng)常遇到的問題。設(shè)變量 與自變量 ( 通常為時間) 之間的關(guān)系點上的 所對應(yīng)的 值，沒有濾除測量噪聲的作用，也不涉及最小二乘法或其它優(yōu)化準(zhǔn)則；只有后者(當(dāng)測量數(shù)據(jù)多于擬合所 選定的曲線的階數(shù)所需要的最小數(shù)目時) 才涉及最小二乘曲 線擬合問題。而對于第 3 種情況，即采樣數(shù)據(jù)的平滑處理，為 了有效地濾除測量噪聲一般采樣數(shù)據(jù)個數(shù)遠遠超過擬合曲線 的階數(shù)，通常采用最小二乘法(或其它優(yōu)化準(zhǔn)則)。可以用 =來表示，實際當(dāng)中一般無法直接得到該解析表達式，通常只能通過測量獲得自變量 離散采樣點 對應(yīng)的函數(shù)值 ，而且測量值一般都不同程度的帶有測量噪

7、聲。曲線擬合的目的就是根據(jù)有限的測量值( , )得到解析表達式=(或者是表達式中的參數(shù))，或者得到最小二乘法的思想是尋找合適的函數(shù)之間的誤差平方和使其與測量值的多項式近似表達式。根據(jù)獲得解析表達式的目的不同曲線擬合通常分以下幾種情況：插值：為了根據(jù)離散的測量數(shù)據(jù)獲得位于 離散采樣點 之間的任意點的 值；或者由于某種原因致使測量過程中 的 取值不準(zhǔn)確(不均勻)，根據(jù)擬合曲線經(jīng)過插值獲得所需要的 點對應(yīng)的 值；外推：利用測量的某一范圍的 對應(yīng)的 值， 根據(jù)擬合曲線獲得測量范圍之外的 對應(yīng)的 值；濾波 (平 滑)：通過曲線擬合減小測量噪聲的影響，獲得離散采樣點 對 應(yīng)的 的更精確的值。上述第 1 和

8、第 2 種情況，擬合過程中所 利用的測量數(shù)據(jù)個數(shù)可以正好滿足擬合所選定的曲線的階數(shù) 所需要的最小數(shù)目，如利用 3 個測量值擬合二次曲線，也可以(1)2 =2= 0達到最小1。當(dāng)信號模型為正弦函數(shù)時(考慮到初始相位可以任意，在信處理中通常余弦和正弦統(tǒng)稱為正弦)，即= sin2 是和 的高度非線性函數(shù)，無法求出使 2 最小的閉+合解，通常只能通過搜索或迭代最小化法1，如三參數(shù)法2-3 和四參數(shù)3-5。搜索迭代法雖然精度較高，但當(dāng)未知參數(shù)較多時，速度較慢。當(dāng)?shù)臏?zhǔn)確表示式未知或直接擬合困難時，也可以采用多項式擬合方法。即在一定范圍之內(nèi)用一個 階多項收稿日期：2007-08-07E-mail：qgq作者

9、簡介：齊國清 (1960)， 男， 遼寧凌海人，博士，教授，研究方向為雷達、通信及圖像信號處理； 呂健 (1957)，男，北京人，副教授，研究方向為電子海圖及船舶導(dǎo)航雷達。：，-式來近似，并將帶入式(1)求出使 2 達到最=1.510.50-0.5-1-1.51.510.50-0.5-1-1.5= 0小的多項式系數(shù)= 0。此時問題的求解過程簡化為求解線性方程組，復(fù)雜程度遠遠低于直接用正弦函數(shù)擬合。最小二乘法多項式擬合是一種常用的方法，當(dāng)測量誤差可以用高 斯白噪聲來表示時，最小二乘準(zhǔn)則與最大似然準(zhǔn)則等價7。文 獻6對最小二乘法多項式擬合應(yīng)用于正弦曲線的擬和進行了 研究。本文對最小二乘法正弦曲線擬

10、合中的一些問題進行研 究和分析。02460246(a) 3 次擬合(b) 5 次擬合1.510.50-0.5-1-1.51.510.50-0.5-1-1.51無測量噪聲情況下擬合多項式階數(shù)對誤差的影響02460246文獻 6 中給出了正弦曲線 19 階多項式擬合系數(shù)，其中9 階擬合系數(shù)是錯誤的，對應(yīng)的擬合曲線和擬合誤差都是錯誤 的，并由此得出了(無測量噪聲條件下，正弦曲線)擬合多項式 階數(shù)越高擬合效果不一定越好的錯誤結(jié)論。按照文獻6中給 出的原始數(shù)據(jù)計算出的 9 次擬合系數(shù)為：a0 = 0.000 000 56，a1 =0.999 88，a2 = 0.001 121 9，a3 = -0.169

11、 729 6，a4 = 0.004 017 2，a5 =0.005 361 59，a6 = 0.001 331 4，a7 = -0.000 566 5，a8 = 0.000 060 4，a9 = -0.000 002 14。擬合結(jié)果如圖 1 所示。另外，文獻6中的7 次擬合系數(shù)略有誤差，因此造成擬合效果較差?？梢姼唠A擬 合要求系數(shù)精度非常高。表 1 給出了文獻6中的數(shù)據(jù)對應(yīng)的311 次擬合的誤差平方和( 2)。(c) 7 次擬合(d) 9 次擬合* 測量值；理想正弦曲線；擬合曲線圖 2 加噪聲時擬合結(jié)果 (= 0.1 )正弦曲線。為了定量描述擬合曲線與理想正弦曲線 (真值曲線)的偏離程度，定義

12、擬合曲線與真值曲線之間的誤差平方和 2(2)2 =2= 0式中 擬合曲線上的點 不含有噪聲測量值(真值)。表 2 列出的是加入 = 0.01高斯白噪聲時擬合誤差的 1000 次平均結(jié)果，其中 2 為擬合曲線與測量值(包括測量噪聲)的誤差平 方和 2 為擬合曲線與正弦曲線真值的誤差平方和( 2 的多次平1.510.50-0.5-1-1.5均值相當(dāng)于均方誤差)?？梢娫谟性肼暤那闆r下，盡管正弦曲線擬合的誤差平方和隨擬合階數(shù)的增加而減小，但階數(shù)過高 時擬合曲線的均方誤差反而增加( 當(dāng)然，實際當(dāng)中由于并不知 道正弦曲線的真值 ，因此并不能利用 2 判斷擬合效果)。在表2 給定的條件下 7 次擬合效果最好

13、。01234567o 測量值圖 1 正弦曲線 9 次擬合結(jié)果表 2 有噪聲情況下的擬合誤差( = 0.01，1000 次平均)表 1 利用文獻 6 數(shù)據(jù)的擬合誤差測量數(shù)據(jù)的有效位數(shù)同樣會影響擬合精度，表 3 分別列出了對文獻 6 的測量數(shù)據(jù)只保留小數(shù)點后 3 位以及增加到 5 位的擬合誤差?？梢?7 階以下擬合精度受測量數(shù)據(jù)有效位數(shù) 影響不大，測量數(shù)據(jù)有效位數(shù)較少時，擬合階數(shù)太高沒有意義?？梢?，在不考慮測量噪聲(或測量噪聲很小) 的情況下，對于正弦曲線的擬合多項式的階數(shù)越高，擬合誤差平方和 2越小。 并不存在文獻6所得到的階數(shù)高了擬合效果反而不好的情況。表 3 測量數(shù)據(jù)有效位數(shù)對擬合誤差的影響

14、2有測量噪聲情況下擬合多項式階數(shù)對誤差的影響一般來說，擬合誤差平方和 2 隨擬合多項式階數(shù)的提高而減小。但應(yīng)當(dāng)注意最小二乘法的核心思想是擬合多項式階數(shù) 不變的情況下尋找使得 2 達到最小的多項式系數(shù)。在給定多 項式階數(shù)并且沒有關(guān)于測量噪聲的任何先驗知識的前提下， 可以認(rèn)為 2 越小擬合效果越好。但多項式階數(shù)不同的情況下，2 小不等于擬合效果就一定好。在有測量噪聲存在的情況下， 多項式階數(shù)越高濾除噪聲的性能越差。圖 2 為在文獻6所用 原始數(shù)據(jù)中加入 = 0.1高斯白噪聲時的擬合結(jié)果，可見有噪聲 的情況下，隨著擬合階數(shù)的增加，擬合曲線越來越明顯的偏離3正弦曲線擬合的外推誤差根據(jù)擬合曲線進行內(nèi)插一

15、般不會帶來太大的誤差，但是 外推的情況要復(fù)雜得多。如果是對直線或二次曲線進行擬合， 小范圍的外推可以得到滿意的結(jié)果。但對于正弦曲線根據(jù)擬 合多項式進行外推效果較差。圖 3 是根據(jù)文獻6中所用的一 個周期 17 個點的測量數(shù)據(jù)分別進行 5、7、9 和 11 階擬合然后 根據(jù)擬合曲線外推 6 點的仿真結(jié)果?？梢娂词?11 階擬合外推-擬合次數(shù)3579112(3 位有效數(shù))0.10594.30e-045.45e-073.12e-071.89e-072(5 位有效數(shù))0.10594.30e-044.23e-076.66e-113.29e-11擬合次數(shù)35791120.10604.3e-044.03e-

16、072.39e-091.43e-09擬合次數(shù)35791120.10730.00158.96e-047.07e-045.12e-0420.10640.00118.25e-040.00100.0013-/-：，：利用參數(shù)估計方法計算出正弦信號的頻率、幅度及相位等參數(shù)，從而獲得其表達式。在一些應(yīng)用當(dāng)中，正弦信號的頻率是 已知的，只有幅度和初始相位是未知的。針對這種頻率已知 的正弦波形可以采用下面基于傅利葉變換的擬合方法 (FT 擬 合)。設(shè)正弦信號可以表示為221100-1-1-2-2(3)= sin0 +，0< 0246802468(a) 5 次擬合(b) 7 次擬合的傅利葉變換為(只考慮頻

17、率正半軸)式中0 已知=在22sin/2=0 處的值0/2+ 0 /2(4)110=00(5)-1-10 =/22-2-2于是可根據(jù)0 分別得到幅度和相位的估計值0246802468= 2= 2(c) 9 次擬合(d) 11 次擬合(6)0表示外推值(7)0 可0圖 3正弦曲線外推結(jié)果式中0 0 在0 處的相角。在實際當(dāng)中=到第 6 點已經(jīng)出現(xiàn)明顯偏差。以用數(shù)值方法計算4多周期正弦曲線的擬合(8)0 =0= 1將幅度和相位的估計值帶入原表達式，得到擬合正弦曲線對于單個周期的正弦波形用 5 階或者 7 階多項式擬合誤差已經(jīng)很小。但是對于多個周期的正弦波形則需要用較高階 數(shù)來擬合。圖 4 所示的兩

18、個周期的正弦波形至少需要 9 階多 項式來擬合。從表 4 所列的 2 個周期正弦曲線擬合誤差可見 多周期正弦曲線擬合的誤差明顯高于(相同階數(shù)) 的單周期正 弦曲線情況 (即使排除波形周期數(shù)增加帶來的測量點數(shù)增加 對 2 的影響)。仿真結(jié)果表明隨著正弦波的周期數(shù)的進一步增 加，擬合誤差迅速增加。因此實際當(dāng)中多項式擬合方法不適 合多周期的正弦波形的擬合。= sin0 +(9)表 5 列出的為正弦信號頻率已知，A = 1，加 = 0.1 的測量噪聲時，利用 FT 擬合的結(jié)果和多項式擬和的 2 對比。 = 0.01時 FT 擬合 2 = 2.35e-04(同條件下多項式擬和結(jié)果見表 2)?？梢娫肼暠尘?/p>

19、中基于 FT 的擬和方法精度優(yōu)于多項式擬合。表 5 頻率已知利用 FT 擬合與多項式擬合誤差對比(= 0.1，單周期 17 點，1000 次平均結(jié)果)1.510.50-0.5-1-1.51.510.50-0.5-1-1.5在頻率已知情況下利用 FT 擬合正弦曲線，不僅算法簡單、擬合精度高，而且外推精度遠遠高于多項式擬合方法。 在測量過程中不僅存在高斯分布的白噪聲，有時還存在脈沖干擾。與多項式擬合相比，F(xiàn)T 擬合抗脈沖噪聲的能力更 強。圖 5 所示為在第 5 個測量點出現(xiàn)一個較大的噪聲時，F(xiàn)T 擬合和多項式擬合 (7 階) 結(jié)果，表 6 所示為 FT 擬合與多項式 擬合 (311 階) 受脈沖干

20、擾引起的誤差對比結(jié)果，可見多項式 擬合受影響較大。05(b) 5 次擬合1005(a) 3 次擬合101.510.50-0.5-1-1.51.510.50-0.5-1-1.505(c) 7 次擬合1005(d) 9 次擬合101100圖 4 兩個周期正弦波的擬合結(jié)果 (無測量噪聲)-1-10 12 3 4 5 (a) FT 擬合6 70 1 2 3 4 56 7 (b) 7 階多項式擬合表 4 2 個周期正弦曲線擬合誤差圖 5 脈沖噪聲對擬和效果的影響5頻率已知正弦曲線的擬合表 6脈沖干擾引起 FT 擬合與多項式擬合偏差對比對于正弦曲線的擬合，除了上面討論的最小二乘多項式擬合方法外，常用的方法

21、還有基于參數(shù)估計的擬合方法。即-<<-<<<-FT 擬合3 次擬合5 次擬合7 次擬合9 次擬合11 次擬合20.020.13630.03840.05880.07690.0879擬合次數(shù)357911215.47622.26030.07438.27e-043.91e-06-FT 擬合3 次擬合5 次擬合7 次擬合9 次擬合11 次擬合20.020.14480.05910.07920.09870.1119-基于 FT 的擬合方法在頻率較低( 即測量數(shù)據(jù)包含的正弦波周期數(shù)較少時)時應(yīng)盡量采用整數(shù)個周期的測量值，否則正 弦信號傅利葉變換的負頻率成分的頻譜泄漏對幅度和相位的

22、 估計精度影響較大。的多項式系數(shù)計算誤差所引起，因為高階多項式對系數(shù)誤差非常敏感。7結(jié)束語最小二乘多項式曲線擬合是一種常用的曲線擬合方法。 這種方法用于正弦曲線擬合時具有一定的局限性。對于一個 周期的正弦曲線擬合，采用最小二乘多項式曲線擬合可以取 得滿意的效果，在測量精度較高情況下，可以采用 7 階或 9 階 擬合，達到較高的擬合精度。測量數(shù)據(jù)有效位數(shù)對高階擬合 的效果影響較大，測量數(shù)據(jù)有效位數(shù)較少時，擬合階數(shù)不宜太 高。在有測量噪聲存在的情況下，擬合誤差平方和 2 不能完全 反映擬合效果，均方誤差能更好地反映擬合曲線與理想曲線 的偏差情況。有較明顯測量噪聲時，不宜采用 9 階以上多項 式擬合

23、，否則不但增加計算量，而且造成擬合均方誤差上升、 曲線平滑程度下降。采用最小二乘多項式擬合正弦曲線的外 推精度較差，不宜作多點外推。對于多周期的正弦曲線需要 較高階數(shù)的多項式進行擬合，例如 2 個周期的正弦波至少需 要 11 階多項式，才能取得較好的擬合效果。對于非正弦曲線， 根據(jù)曲線的具體類型不同，擬合多項式的階數(shù)越高不一定性 能越好。在實際當(dāng)中經(jīng)常遇到頻率已知的正弦曲線的擬合問 題。針對這種特定條件可以利用傅利葉變換估計其幅度和相 位實現(xiàn)曲線擬合，這種方法不但算法簡單、易于實現(xiàn)，而且精 度高、適應(yīng)性強。6非正弦曲線擬合階次的影響以上討論的都是正弦曲線的擬合情況。根據(jù)上面分析，對于正弦曲線，

24、在不考慮測量噪聲和計算誤差的情況下，擬合 多項式的階數(shù)越高，擬合曲線與實際測量值越接近。但如果 測量曲線本身并不是正弦曲線，而是一個多項式曲線，則如文 獻 8 所驗證的，并非擬合多項式的階數(shù)越高越好。以下面 5 次曲線為例110005= +332768根據(jù) 從 0 到 20，按上式計算的 值進行 39 階多項式擬合。擬合多項式系數(shù)及擬合誤差如表 7 所示，圖 6 為 5 階擬合 結(jié)果，11 次擬合的誤差為 1.6226e-013?？梢?5 次擬合效果最 好。由于曲線本身為 5 階多項式，當(dāng)用高于 5 階的多項式擬 合時，高于 5 階的系數(shù)并不為零，因此造成擬合誤差增加。100-10-20-30

25、-40-50-60-70參考文獻:051015201Steven M Kay. 統(tǒng)計信號處理基礎(chǔ)-估計與檢測理論M. 羅鵬飛,譯.北京:電子工業(yè)出版社,2003:182-225.IEEE Std 1057-1994, IEEE Standard for digitizing waveform re- cordersS.IEEE Standards Board,1994.Peter Handel.Properties of the IEEE-STD-1057 four-parameter sine wave fit algorithmJ.IEEE Transactions on Instrume

26、ntation and Measurement,2000,49(6):1189-1193.梁志國. 四參數(shù)正弦波曲線擬合的快速算法 J. 計測技術(shù),2006,26(1):4-7. 吳義華.正弦信號四參數(shù)的高精度估計算法J.中國科技大學(xué) 學(xué)報,2006,35(6):625-629. 司少玲,關(guān)永.三角函數(shù)曲線數(shù)據(jù)擬合最佳次數(shù)的確定J.計算 機工程與設(shè)計,2006,27(24):4660-4662.István Kollár,Jerome J Blair.Improved determination of the best fitting sine wave in ADC testingJ.IEEE Transactions on Instru- mentation and Measurement,2005,54(5):1978-1983