孫文軍，芮國勝，王 林，田文飚

一種利用Duffing-VanderPol振子估計(jì)弱信號相位的方法*

孫文軍**，芮國勝，王 林，田文飚

（海軍航空工程學(xué)院電子信息工程系，山東煙臺264001）

根據(jù)Duffing-VanderPol振子大周期態(tài)時(shí)的最大LYaPunov指數(shù)與系統(tǒng)策動(dòng)力的關(guān)系，提出了一種基于該振子檢測微弱信號相位參量的方法。該方法利用系統(tǒng)策動(dòng)力與待測信號初始相位有關(guān)的特性，得到待測信號的初始相位，并將相位檢測過程中噪聲對系統(tǒng)的影響作為影響因子加以考慮。理論分析及數(shù)值仿真結(jié)果表明，該方法工程實(shí)現(xiàn)簡單，抗噪性能好。當(dāng)待測信號信噪比為-30 dB時(shí)，均方誤差僅為0.1，具有較高的估計(jì)精度。

微弱信號檢測；相位檢測；Duffing-VanderPol振子；LYaPunov指數(shù)

1 引 言

隨著現(xiàn)代數(shù)字通信技術(shù)的發(fā)展和各種應(yīng)用領(lǐng)域的需要，強(qiáng)噪聲背景下微弱信號的檢測已成為研究的重點(diǎn)。由于混沌振子對小信號的極端敏感性和對噪聲的免疫力，混沌理論在微弱信號檢測中得到了廣泛應(yīng)用[1-4]。在不同背景噪聲下，出現(xiàn)了各種基于混沌振子的微弱信號檢測理論和方法[5-6]，由于混沌振子檢測具有靈敏度高、抗噪聲能力強(qiáng)、設(shè)備成本低的特點(diǎn)，因而受到人們的廣泛關(guān)注。

許多研究者對這一領(lǐng)域進(jìn)行了深入研究。文獻(xiàn)[7]利用混沌振子法實(shí)現(xiàn)了對低信噪比環(huán)境下頻移鍵控信號解調(diào)；文獻(xiàn)[8]將多混沌振子陣列應(yīng)用于信號的檢測，實(shí)現(xiàn)了信號在不同初始相位下的檢測。近年來，弱信號檢測領(lǐng)域利用混沌振子的參數(shù)敏感特性取得了一些新的研究成果：文獻(xiàn)[9]通過對水聲信號的非線性特征進(jìn)行分析并對其進(jìn)行系統(tǒng)建模，將微弱水下目標(biāo)信號的檢測能力提升了約10 dB；文獻(xiàn)[10]將自相關(guān)和混沌理論相結(jié)合，進(jìn)一步增強(qiáng)了混沌振子的弱信號頻域檢測能力；文獻(xiàn)[11-12]基于主動(dòng)控制策略，通過設(shè)計(jì)比例微分控制器消除誤差系統(tǒng)中的非線性項(xiàng)，以獲得混沌系統(tǒng)的準(zhǔn)確同步。這些研究成果展示了混沌振子在弱信號檢測領(lǐng)域尤其是參數(shù)辨識等方面具有巨大的潛在研究價(jià)值，但仍存在一些問題：一方面，大多數(shù)文獻(xiàn)更多關(guān)注的是弱信號是否存在的物理檢測，而未對信號的相位檢測進(jìn)行論述；另一方面，基于數(shù)值分析理論的控制類方法需要設(shè)計(jì)復(fù)雜的代價(jià)函數(shù)，實(shí)現(xiàn)難度大，無法實(shí)現(xiàn)快速判決。

微弱信號的相位檢測一直是弱信號檢測領(lǐng)域一個(gè)非常重要的問題。本文通過對Duffing-Vander_ Pol振子在大周期態(tài)時(shí)的LYaPunov指數(shù)進(jìn)行研究，得到其與系統(tǒng)策動(dòng)力的關(guān)系，為了能夠在較低信噪比條件下獲得較好的檢測性能，考慮了信號噪聲的影響，得到信號等效策動(dòng)力與信號相位的關(guān)系，最終得到LYaPunov指數(shù)與待檢測信號相位的關(guān)系。仿真實(shí)驗(yàn)結(jié)果表明，該方法具有較高的檢測精度，并且有較好的抗噪聲性能。

2 Duffing-Vanderpol系統(tǒng)信號檢測原理

考慮如下Duffing-VanderPol方程的弱信號檢測基本形式：

式中：ω為待檢測周期信號的頻率；k為阻尼比；-x（t）＋x3（t）為非線性恢復(fù)力；fcos（ωt）為內(nèi)置周期策動(dòng)力；rcos（ωt＋Δφ）是待測周期弱信號；n（t）為噪聲信號。傳統(tǒng)的信號檢測方式以Duffing-Vander_ Pol系統(tǒng)對噪聲免疫[13]為前提，將待測信號作為混沌系統(tǒng)周期策動(dòng)力的攝動(dòng)，在噪聲強(qiáng)烈的情況下，對系統(tǒng)狀態(tài)的改變影響很小，而一旦有特定頻率的待測信號，由于混沌系統(tǒng)對周期小信號的敏感性，即使幅值較小，也會(huì)使系統(tǒng)發(fā)生相變，在相位差滿足如下條件時(shí)，即π-arccos（r/（2f））≤Δφ≤π＋arccos（r/（2f）），通過判斷系統(tǒng)是否發(fā)生了相變來判定信號是否存在，從而達(dá)到在強(qiáng)背景噪聲下對微弱周期小信號進(jìn)行檢測的目的。上述信號檢測方法對信號相位只考慮了一定的檢測范圍，沒有具體的相位差值。可應(yīng)用多個(gè)內(nèi)置不同初始相位策動(dòng)力的混沌振子組成的陣列，對未知相位正弦信號進(jìn)行并行檢測，其檢測原理如圖1所示。

圖1 Duffing-VanderPol振子陣列檢測原理圖Fig.1 The detection sYstem comPosed of multiPle Duffing-VanderPol oscillators arraY

圖1 中，通過n個(gè)不同初始相位的振子檢測，對每個(gè)振子是否相變進(jìn)行判斷，綜合振子的狀態(tài)可以將待檢測信號的相位確定在一個(gè)較小范圍內(nèi)，從而達(dá)到弱信號相位檢測的目的。這種檢測方法使用的Duffing-VanderPol方程數(shù)量越多，檢測出的信號相位范圍越小，系統(tǒng)的靈敏度也越高，但是隨著方程數(shù)量的增多，檢測復(fù)雜度會(huì)大大增加，降低了檢測效率，并且此方法也沒有考慮到噪聲對系統(tǒng)相變判斷的影響。

3 本文相位檢測方法

本文綜合考慮Duffing-VanderPol系統(tǒng)LYaPunov指數(shù)與周期策動(dòng)力的關(guān)系以及周期策動(dòng)力與相位檢測的關(guān)系，提出基于Duffing-VanderPol振子的正弦信號相位檢測方法。

3.1 待檢測信號相位與幅值關(guān)系

由式（1）等號的右端可知，對混沌振子起相變作用的總周期驅(qū)動(dòng)信號為

式中：

由于待測信號的幅值r遠(yuǎn)遠(yuǎn)小于策動(dòng)力臨界值f，由式（4）可知，Δθ很小，所以對系統(tǒng)狀態(tài)的影響可以忽略不計(jì)，因此系統(tǒng)狀態(tài)是否改變的關(guān)鍵在于Fs與噪聲干擾n（t）的大小。由式（3）可知，系統(tǒng)是否會(huì)發(fā)生相變與待檢測信號和內(nèi)部驅(qū)動(dòng)信號的相位差有關(guān)，由于式（3）中cosΔφ與Fs存在非線性關(guān)系，不便于對相位進(jìn)行估計(jì)，所以對其進(jìn)行線性化可得

進(jìn)一步可得

通過上式可知，一旦信號幅度確定，系統(tǒng)狀態(tài)是否發(fā)生變化將由相位差值決定。因此，可以通過定量的系統(tǒng)相變判別方法來進(jìn)一步確定待測信號的相位差。

3.2 Lyapunov指數(shù)判別方法

由于利用系統(tǒng)相圖進(jìn)行判斷需要人工參與，且缺乏定量判別，并不適合工程應(yīng)用，所以實(shí)際中通常利用最大LYaPunov指數(shù)判定系統(tǒng)狀態(tài)。

對于映射系統(tǒng)xn＋1＝f（xn），考慮初值x0和其臨近值x0＋Δx0，其LYaPunov指數(shù)可表示為

式中：h為采樣間隔；n為采樣時(shí)間，且n＝mh。式（8）為計(jì)算連續(xù)系統(tǒng)的LYaPunov指數(shù)的方法。應(yīng)用LYaPunov指數(shù)對系統(tǒng)的狀態(tài)進(jìn)行判斷，為了能較為準(zhǔn)確，需要選取指數(shù)譜中穩(wěn)定的LYaPunov指數(shù)值，對于式（1）所示的Duffing-VanderPol系統(tǒng)來說，在Duffing-VanderPol方程阻尼比k＝0.5且h＝0.01、m＝30 000條件下，本文選取n＝300時(shí)的最大LYa_ Punov指數(shù)值，系統(tǒng)最大LYaPunov指數(shù)與系統(tǒng)的周期驅(qū)動(dòng)力強(qiáng)度F的關(guān)系如圖2所示。

圖2 系統(tǒng)最大LYaPunov指數(shù)與周期驅(qū)動(dòng)力強(qiáng)度F的關(guān)系Fig.2 The relationshiP betWeen the largest LYaPunov exPonent and Periodic driving force strength

圖2 中結(jié)果可以說明，Duffing-VanderPol振子的最大LYaPunov指數(shù)隨著系統(tǒng)的周期驅(qū)動(dòng)力強(qiáng)度F的增強(qiáng)近似呈單調(diào)遞減的關(guān)系。當(dāng)最大LYaPunov特性指數(shù)大于零時(shí)，說明系統(tǒng)處于混沌態(tài)；當(dāng)系統(tǒng)最大LYaPunov特性指數(shù)由大于零轉(zhuǎn)為小于零，則是系統(tǒng)從混沌態(tài)躍變到了周期態(tài)的標(biāo)志。

考慮實(shí)際計(jì)算LYaPunov指數(shù)時(shí)會(huì)受到輸入噪聲干擾，根據(jù)文獻(xiàn)[14]對噪聲干擾的分析，噪聲對LYaPunov指數(shù)的影響是由于輸入噪聲的方差導(dǎo)致相圖中增加了方差為的噪聲，且與系統(tǒng)采樣間隔h有關(guān)，即與成正比，根據(jù)式（8），LYa_ Punov指數(shù)為n次隨機(jī)變量的和，增加采樣時(shí)間n會(huì)相應(yīng)減少LYaPunov指數(shù)的方差，且隨機(jī)變量的噪聲符合正態(tài)分布，而不會(huì)隨機(jī)變量服從正態(tài)分布的特性，所以LYaPunov指數(shù)也服從正態(tài)分布。令LYaPunov指數(shù)的方差為，則其與/n成正比，所以噪聲對LYaPunov指數(shù)方差的影響為

通常Duffing-VanderPol振子處于大周期態(tài)時(shí)，周期驅(qū)動(dòng)力強(qiáng)度F與λ的關(guān)系可認(rèn)為是線性的，則噪聲對系統(tǒng)LYaPunov指數(shù)的影響可認(rèn)為對等效總策動(dòng)力的改變，即增加了噪聲的等效策動(dòng)力Fn，則

通過線性最小二乘擬合的方法，可以得到大周期態(tài)時(shí)周期驅(qū)動(dòng)力強(qiáng)度與LYaPunov指數(shù)的關(guān)系公式為

進(jìn)一步推導(dǎo)可得

式中：λp為混沌態(tài)到大尺度周期態(tài)的最大LYaPunov指數(shù)閾值；α、b為擬合參數(shù)。在圖1所示的條件下，根據(jù)對F的統(tǒng)計(jì)結(jié)果可得，λp＝-0.084 5，與其對應(yīng)的等效總策動(dòng)力幅值為Fp＝0.837 0，K＝0.897 9，擬合參數(shù)分別為α＝0.165 5，b＝0.823 1，其擬合誤差如圖3所示。

圖3 周期驅(qū)動(dòng)力強(qiáng)度擬合的標(biāo)準(zhǔn)差Fig.3 The adaPtive mean-square error of Periodic driving force

由圖可見，在Duffing-VanderPol系統(tǒng)處于穩(wěn)定的大尺度周期狀態(tài)時(shí)，對系統(tǒng)等效總策動(dòng)力擬合的誤差很小，擬合函數(shù)具有很高的精度，據(jù)此可以利用其估計(jì)待測微弱信號的相位。

3.3 根據(jù)Lyapunov指數(shù)估計(jì)微弱信號的相位

由式（7）和式（8），可以得到大尺度周期態(tài)時(shí)LYaPunov指數(shù)與待檢測信號相位是線性關(guān)系，所以，為了提高待檢測相位的估計(jì)精度，對內(nèi)置正弦策動(dòng)力f設(shè)置等間隔分布的多個(gè)值，對其加和平均進(jìn)行相位估計(jì)，設(shè)f分別為f1，f2，…，fN，且fi-r≥Fp，這樣可以保證Duffing-VanderPol振子處于大尺度周期態(tài)，每一個(gè)fi對應(yīng)一個(gè)最大LYaPunov指數(shù)λi，可得相位估計(jì)公式

上述表達(dá)式對于相位的估計(jì)可以提高一定的精度，并且有利于工程實(shí)現(xiàn)。根據(jù)以上分析，對信號相位檢測的步驟歸納如下：

（1）設(shè)置內(nèi)置策動(dòng)力：根據(jù)圖1所示系統(tǒng)最大LYaPunov指數(shù)與等效策動(dòng)力的關(guān)系，設(shè)置內(nèi)置驅(qū)動(dòng)力f1，令f1＝Fp＋r，將式（1）的f替換成f1，得到系統(tǒng)的檢測方程；

（2）計(jì)算最大LYaPunov指數(shù)：對第一步中的檢測方程，應(yīng)用歐拉-丸山方法進(jìn)行數(shù)值求解，積分步長為0.01，求得方程的數(shù)值解，根據(jù)3.2節(jié)的方法，進(jìn)一步得到方程的最大LYaPunov指數(shù)λ1；

（3）重新設(shè)置內(nèi)置驅(qū)動(dòng)力：令fi＝fi-1＋0.001，重復(fù)步驟2，得到在不同內(nèi)置策動(dòng)力下多個(gè)最大LYa_ Punov指數(shù)λi；

（4）估計(jì)弱信號相位：根據(jù)公式（16）計(jì)算得到弱信號相位Δφ的估計(jì)。

4 仿真分析

設(shè)置待檢測正弦信號的幅值為0.001 V，加入高斯白噪聲，信噪比為-5 dB，分別設(shè)置不同的信號相位，，對Δφ進(jìn)行估計(jì)，并與設(shè)置的信號相位進(jìn)行對比，結(jié)果如圖4所示。

圖4 微弱信號相位檢測結(jié)果Fig.4 The Phase detection results of Weak signals

圖4 中，1～17監(jiān)測點(diǎn)分別是初始相位為0～π的微弱正弦信號，由圖可見，應(yīng)用本文方法可以對弱信號相位進(jìn)行較好的估計(jì)。取其中較有代表性的9點(diǎn)對其估計(jì)的均方誤差做統(tǒng)計(jì)，結(jié)果如表1所示。

表1 不同初始相位的檢測誤差Tab.1 The detection error of Weak signals With different initial Phases

通過觀察各估計(jì)點(diǎn)的均方誤差可以發(fā)現(xiàn)，除了相位差在0和π附近的偏差較大，其他相位估計(jì)的均方誤差均在0.01以下，這是由于相位差在0和π處的反余弦函數(shù)斜率接近于0，函數(shù)曲線靠近漸近線，擬合參數(shù)誤差會(huì)對估計(jì)結(jié)果造成較大影響。

進(jìn)一步對檢測系統(tǒng)抗噪性能進(jìn)行研究，對信號設(shè)置9個(gè)不同的初始相位值，在不同的信噪比條件下，得到的相位檢測誤差如圖5所示。

圖5 不同信噪比的相位檢測誤差Fig.5 The Phase detection error of Weak signals under different SNR conditions

從結(jié)果中可以看出，該檢測方法對噪聲有一定的免疫力，在信噪比-30 dB條件下，均方誤差在0.1以下，仍可以有效檢測出信號的初始相位，這是由于本文方法在擬合策動(dòng)力與最大LYaPunov指數(shù)的關(guān)系時(shí)，將噪聲對LYaPunov指數(shù)的影響作為擬合項(xiàng)加以考慮，這樣使在較低信噪比條件下的弱信號相位檢測精度得到相應(yīng)提高。

5 結(jié)束語

本文研究了一種基于Duffing-VanderPol振子的正弦信號相位檢測方法，該方法基于系統(tǒng)處于大周期態(tài)時(shí)的最大LYaPunov指數(shù)對參數(shù)仍有極端敏感的特性，通過這種特性進(jìn)行定量的描述，找到最大LYaPunov指數(shù)與弱信號相位之間的函數(shù)關(guān)系，并在實(shí)現(xiàn)相位檢測時(shí)考慮噪聲的影響，得到信號相位的估計(jì)式。仿真實(shí)驗(yàn)表明，該檢測方法精度較高，易于實(shí)現(xiàn)，并且具有很強(qiáng)的抗噪聲能力。

文中算法提出的依據(jù)是基于混沌振子相變和最大LYaPunov指數(shù)準(zhǔn)則的，當(dāng)混沌振子沒有發(fā)生可靠相變或者最大LYaPunov指數(shù)無法實(shí)現(xiàn)相變檢測時(shí)，本文方法同樣也不再適用。

[1] MURALI K，DITTO W L.Realization of the fundamental nor gate using a chaotic circuit[J].PhYsical RevieW E，2013，68（1）：1-4.

[2] FIERENS P I，IBANEZ S A，Perazzo R P J.A memorY de_ vice sustained bY noise[J].PhYsics Letters A，2014，374（22）：2207-2209.

[3] 芮國勝，張洋，苗俊，等.聯(lián)合增益遞推的Duffing系統(tǒng)弱信號檢測算法[J].電子學(xué)報(bào)，2012，40（6）：1269-1273. RUI Guosheng，ZHANG Yang，MIAO Jun，et al.A Weak signal detection method bY duffing sYstem With the gain [J].Acta Electronica Sinica，2012，40（6）：1269-1273.（in Chinese）

[4] 張淑清，師榮艷，李盼，等.基于混沌關(guān)聯(lián)積分的暫態(tài)電能質(zhì)量擾動(dòng)分類[J].儀器儀表學(xué)報(bào)，2015（1）：160-166. ZHANG Shuqing，SHI RongYan，LI Pan，et al.Transient PoWer qualitY disturbance classification based on chaos correlation integral[J].Chinese Journal of Scientific In_ strument，2015（1）：160-166.（in Chinese）

[5] 李月，石要武，馬海濤，等.湮沒在色噪聲背景下微弱方波信號的混沌檢測方法[J].電子學(xué)報(bào)，2004，32（1）：87-90. LI Yue，SHI YaoWu，MA Haitao，et al.Chaotic detection method for Weak square Wave signal submerged in colored noise[J].Acta Electronica Sinica，2004，32（1）：87-90.（in Chinese）

[6] 王文波，汪祥莉.噪聲模態(tài)單元預(yù)判的經(jīng)驗(yàn)?zāi)B(tài)分解脈沖星信號消噪[J].物理學(xué)報(bào)，2013，13（7）：1-12. WANG Wenbo，WANG Xiangli.EmPirical mode decomPo_ sition Pulsar signal de-noising method based on Predic_ ting of noise mode cell[J].Acta PhYsica Sinica，2013，13（7）：1-12.（in Chinese）

[7] WANG J，REN X，ZHANG S W，et al.AdaPtive bistable stochastic resonance aided sPectrum sensing[J].IEEE Transactions on Wireless Communication，2014，3（7）：4014-4024.

[8] LEANDRODOS S F，GUERRA A，BATISTELA N J，et al. Multi-objective cuckoo search algorithm based on duffing oscillator aPPlied to hYsteresis Parameters estimation[J]. IEEE Transactions on Magnetics，2013，49（5）：1745-1748.

[9] 劉海波，吳德偉，金偉，等.Duffing振子微弱信號檢測方法研究[J].物理學(xué)報(bào)，2013，62（5）：42-47. LIU Haibo，WU DeWei，JIN Wei，et al.StudY on Weak sig_ nal detection method With duffing oscillators[J].Acta PhYsica Sinica，2013，62（5）：42-47.（in Chinese）

[10] 田晶晶，李世武，蘇建，等.基于經(jīng)驗(yàn)?zāi)B(tài)分解的載貨汽車載荷動(dòng)態(tài)檢測策略研究[J].振動(dòng)與沖擊，2013，32（4）：173-178. TIAN Jingjing，LI ShiWu，SU Jian，et al.DYnamic detec_ tion strategY of vehicle load based on EMD[J].Journal of Vibration and Shock，2013，32（4）：173-178.（in Chinese）

[11] 孫自強(qiáng)，陳長征，谷艷玲，等.基于混沌和取樣積分技術(shù)的大型風(fēng)電增速箱早期故障診斷[J].振動(dòng)與沖擊，2013，32（9）：113-117. SUN Ziqiang，CHEN Changzheng，GU Yanling，et al.In_ ciPient fault diagnosis of large scale Wind turbine gearbox based on chaos theorY and samPling integral technologY [J].Journal of Vibration and Shock，2013，32（9）：113-117.（in Chinese）

[12] 徐一紅，祝長生，趙耀培.基于稀疏分解和混沌理論的微弱信號檢測[J].電訊技術(shù)，2015，55（11）：1194-1199. XU Yihong，ZHU Changsheng，ZHAO YaoPei.Weak signal detection based on sParse decomPosition and cha_ otic theorY[J].Telecommunication Engineering，2015，55（11）：1194-1199.（in Chinese）

[13] LANDAUER R.Reversible comPuting and PhYsical magne_ tometer[J].PhYsics TodaY，2014，45（3）：100-103.

[14] 馬知恩，周義倉.常微分方程定性與穩(wěn)定性方法[M].北京：科學(xué)出版社，2001：55-60. MA Zhien，ZHOU Yicang.StabilitY and qualitY analYsis of nonlinear differential equation[M].Beijing：Science Press，2001：25-30.（in Chinese）

[15] 聶春燕，石要武.基于互相關(guān)檢測和混沌理論的弱信號檢測方法研究[J].儀器儀表學(xué)報(bào)，2001，22（1）：32-35. NIE ChunYan，SHI YaoWu.The research of Weak signal detection based on cross-correlation and chaos theorY [J].Chinese Journal of Scientific Instrument，2001，22（1）：32-35.（in Chinese）

孫文軍（1987—），男，山東聊城人，2012年于海軍航空工程學(xué)院獲碩士學(xué)位，現(xiàn)為博士研究生，主要研究方向?yàn)榛煦缛跣盘枡z測和現(xiàn)代通信系統(tǒng)；

SUN Wenjun Was born in Liaocheng，Shan_ dong Province，in 1987.He received the M.S. degree from Naval Aeronautical and Astronautical UniversitY in 2012.He is currentlY Working toWard the Ph.D.de_ gree.His research concerns Weak signals detection based on cha_ otic theorY and modern Wireless communication sYstem.

Email：djresearch@126.com

芮國勝（1968—），男，江蘇南京人，2001年于哈爾濱工程大學(xué)獲博士學(xué)位，現(xiàn)為海軍航空工程學(xué)院教授、博士生導(dǎo)師，主要研究方向?yàn)楝F(xiàn)代通信系統(tǒng)、非線性濾波理論、小波理論與應(yīng)用等；

RUI Guosheng Was born in Nanjing，Jiangsu Province，in 1968.He received the Ph.D.degree from Harbin Engineering U_ niversitY in 2001.He is noW a Professor and also the Ph.D.su_ Pervisor.His research concerns modern communication sYstem，nonlinear filter theorY，Wavelet theorY and its aPPlication，etc.

王 林（1985—），男，山東煙臺人，博士，講師；

WANG Lin Was born in Yantai，Shandong Province，in 1985.He is noW a lecturer With the Ph.D.degree.

田文飚（1987—），男，江西南昌人，博士，講師。

TIAN Wenbiao Was born in Nanchang，Jiangxi Province，in 1987.He is noW a lecturer With the Ph.D.degree.

Estimation of Weak Signal Phase by Using Duffing-Vanderpol Oscillator

SUN Wenjun，RUI Guosheng，WANG Lin，TIAN Wenbiao
（DePartment of Electronic and Information Engineering，Naval Aeronautical and Astronautical UniversitY，Yantai 264001，China）

When Duffing-VanderPol oscillator is in large-scale Periodic motion，there are some nexus be_ tWeen the largest LYaPunov exPonent and sYstem driving force.According to this conclusion，an algorithm for detecting Weak signals’Phase based on Duffing-VanderPol oscillator is ProPosed.Through finding the relationshiP betWeen sYstem driving force and the Phase of Weak signal，the Weak signals’Phase is detec_ ted.The noise imPact on sYstem is also considered.Theoretical analYsis and simulation result shoW that this algorithm is easY to be realized in engineering With strong restrain abilitY to noise and high estimation accu_ racY.When the signal-to-noise ratio（SNR）of to-be-detected signals is equal to-30 dB，the meansquare error（MSE）is onlY 0.1.

Weak signal detection；Phase detection；Duffing-VanderPol oscillator；LYaPunov exPonent

The National Natural Science Foundation of China（No.41476089）

TN911

A

1001-893X（2016）01-0014-06

10.3969/j.issn.1001-893x.2016.01.003

孫文軍，芮國勝，王林，等.一種利用Duffing-VanderPol振子估計(jì)弱信號相位的方法[J].電訊技術(shù)，2016，56（1）：14-19.[SUN Wenjun，RUI Guosheng，WANG Lin，et al.Estimation of Weak signal Phase bY using Duffing-VanderPol oscillator[J].Telecommunication En_ gineering，2016，56（1）：14-19.]

2015-09-15；

2015-12-02 Received date：2015-09-15；Revised date：2015-12-02

國家自然科學(xué)基金資助項(xiàng)目（41476089）

**通信作者：djresearch@126.com Corresponding author：djresearch@126.com