石 敏
(1. 水聲對抗技術(shù)重點實驗室, 廣東 湛江 524022;2. 中國人民解放軍91388部隊, 廣東 湛江 524022)
混沌臨界值對混沌系統(tǒng)檢測性能的影響分析
石 敏1,2
(1. 水聲對抗技術(shù)重點實驗室, 廣東 湛江 524022;2. 中國人民解放軍91388部隊, 廣東 湛江 524022)
利用混沌系統(tǒng)進行微弱信號檢測時,確定系統(tǒng)由混沌態(tài)轉(zhuǎn)化為大尺度周期態(tài)的混沌臨界閾值至關(guān)重要。仿真表明:選取不同精度的混沌臨界閾值將對混沌系統(tǒng)檢測性能產(chǎn)生影響,選取的臨界閾值精度越高,能夠檢測到信號的信噪比越低,但同時其對噪聲檢測的虛警率也越高。因此,實際檢測中需根據(jù)檢測要求來選取合適精度的臨界閾值,以使檢測性能達到檢測要求。
微弱信號檢測;混沌系統(tǒng)檢測性能;混沌臨界值
利用混沌系統(tǒng)進行微弱信號檢測的基本原理是混沌系統(tǒng)具有對噪聲及與其內(nèi)置信號頻差較大信號免疫,而對與混沌系統(tǒng)內(nèi)置信號頻差較小信號敏感的特點[1–5]。
當混沌系統(tǒng)處于由混沌狀態(tài)轉(zhuǎn)化為大尺度周期狀態(tài)的臨界狀態(tài)時,若輸入噪聲則系統(tǒng)仍保持混沌狀態(tài),而輸入與之同頻率的微弱信號,則系統(tǒng)迅速進入大尺度周期狀態(tài),據(jù)此可以檢測待測信號中是否含有微弱信號[6–10]。
在混沌微弱信號檢測中,如何確定混沌臨界值至關(guān)重要[11]。仿真中發(fā)現(xiàn):選取不同精度的臨界值將影響混沌系統(tǒng)的檢測性能,本文主要研究不同精度的臨界值將如何影響系統(tǒng)的檢測性能。
考慮可用于任意頻率微弱信號檢測的 duffing 混沌檢測系統(tǒng),其 duffing 方程為:
式中,ω 和γ分別為混沌系統(tǒng)內(nèi)置信號頻率和幅度;k為阻尼系數(shù),一般取 k=0.5;(x–x3)為非線性項和分別為 x 的一階和二階微分。
隨著γ值的變化,系統(tǒng)的相軌跡將發(fā)生變化,歷經(jīng)同宿軌道、倍周期分叉軌道、混沌狀態(tài)及大尺度周期狀態(tài)等,當γ等于混沌臨界閾值 γd時,系統(tǒng)進入混沌臨界狀態(tài)。此時若輸入純噪聲,系統(tǒng)仍保持混沌狀態(tài),若輸入含有與內(nèi)置信號頻率相同的信號時,系統(tǒng)迅速從混沌狀態(tài)轉(zhuǎn)化為大尺度周期狀態(tài),據(jù)此可檢測微弱信號是否存在。
Lyapunov 指數(shù)可用于定量確定系統(tǒng)的狀態(tài)和混沌臨界值:系統(tǒng)的最大 Lyapunov 特性指數(shù)大于 0,說明系統(tǒng)處于混沌狀態(tài);當系統(tǒng)的最大 Lyapunov 特性指數(shù)小于 0 時,說明系統(tǒng)處于大尺度周期狀態(tài);當 Lyapunov 特性指數(shù)由大于 0 轉(zhuǎn)為小于 0,則說明系統(tǒng)從混沌態(tài)躍變到了周期態(tài),最大 Lyapunov 特性指數(shù)符號轉(zhuǎn)變的那一刻所對應(yīng)的系統(tǒng)內(nèi)置信號幅值即為系統(tǒng)臨界閾值。文獻[12]中給出了 Lyapunov 指數(shù)的具體計算方法及利用該指數(shù)計算系統(tǒng)臨界閾值的方法。
表1 11 個采樣點對應(yīng)的 Lyapunov 指數(shù)值Tab. 1 Lyapunov characteristic exponents corresponding to 11 sample points
圖1 系統(tǒng)對應(yīng)的時域圖和相軌跡圖Fig. 1 Time-domain diagramand phase contrail diagramof chaotic system
由表1 可知,小數(shù)點后取 3 位的混沌臨界閾值為γd=0.825,由圖1 中 γ=0.825 和 γ=0.826 對應(yīng)的時域圖和相軌跡圖也可以判定 0.825 是混沌臨界閾值。
若要取更高精度的臨界閾值,則可將 [0.825, 0.826] 以 0.000 1 為步長進行等分,同樣計算各采樣點對應(yīng)的 Lyapunov 指數(shù),根據(jù) Lyapunov 指數(shù)值得到更高精度的臨界閾值,得到其精度取小數(shù)點后 4 位對應(yīng)的臨界閾值為 0.825 8。該過程一直持續(xù)下去,直到得到預(yù)先設(shè)置的精度。這里以小數(shù)點后 8 位為設(shè)置精度,得到各采樣點對應(yīng)的 Lyapunov 指數(shù)如表2 所示。
由表2 可知,小數(shù)點后取 8 位對應(yīng)的混沌臨界閾值為 γd=0.825 843 00。
表2 12 個采樣點對應(yīng)的 Lyapunov 指數(shù)值Tab. 2 Lyapunov characteristic exponents corresponding to 12 sample points
由第 2 節(jié)可知式(1)中的 duffing 混沌系統(tǒng)的臨界閾值按不同精度可取為 0.82,0.825,0.825 8 和 0.825 843 00。
(6)收益轉(zhuǎn)讓:支付期內(nèi),無論發(fā)明人、創(chuàng)建人或作者、參與者是否在慕尼黑工業(yè)大學任職,相關(guān)費用都由慕尼黑工業(yè)大學支付。相關(guān)人員如果離開慕尼黑工業(yè)大學,不能帶走其收益部分,而是要轉(zhuǎn)給所在部門。如果調(diào)入慕尼黑工業(yè)大學另一所屬機構(gòu),其收益可隨之轉(zhuǎn)入。
3.1 只輸入高斯白噪聲時對混沌系統(tǒng)檢測性能的影響
將 duffing 系統(tǒng)處于混沌臨界狀態(tài),即將式(1)中的γ分別取為 0.82,0.825,0.825 8,0.825 843 00時,在系統(tǒng)中加入高斯白噪聲,式(1)變?yōu)槿缦滦问剑?/p>
式中:γd為混沌臨界閾值;n(t)為高斯白噪聲。
理論上,高斯白噪聲輸入到處于混沌臨界狀態(tài)的混沌系統(tǒng)中,系統(tǒng)仍保持混沌狀態(tài)不變,即系統(tǒng)對應(yīng)的最大 Lyapunov 指數(shù)應(yīng)保持為正數(shù)。但仿真實驗中發(fā)現(xiàn)在不同精度的臨界閾值下輸入一定功率的白噪聲也會引起系統(tǒng)的狀態(tài)變化,即系統(tǒng)對應(yīng)的最大 Lyapunov指數(shù)可能會變?yōu)樨摂?shù)。
以輸入噪聲功率為 10–6為例,分別取 γd=0.82, 0.825,0.825 8,0.825 843 00,每個臨界閾值下進行100 次仿真實驗,得到噪聲功率為 10–6時,各仿真中系統(tǒng)對應(yīng)的最大 Lyapunov 指數(shù)如圖2 所示。
圖2 在不同精度臨界閾值下,系統(tǒng)輸入功率為 10–6的高斯白噪聲后,系統(tǒng)對應(yīng)的最大 Lyapunov 指數(shù)Fig. 2 Maximal Lyapunov exponent of chaotic systemafter inputting gauss noise with power 10–6into the systemunder different precision critical value
定義輸入噪聲時,系統(tǒng)的虛警率為:
由圖2 可知,γd=0.82,0.825,0.825 8,0.825 843 00時,對應(yīng)的虛警率分別為 0%,0%,7%,62%。隨著混沌臨界閾值的不斷精確,系統(tǒng)對噪聲越來越敏感,其虛警率增加,原因如下:
1)噪聲的頻帶很寬,其包含外力驅(qū)動力頻率的窄帶可能導(dǎo)致系統(tǒng)由混沌臨界狀態(tài)跳變到大尺度周期狀態(tài);
3.2 輸入含噪信號時對檢測性能的影響
在式(1)的 duffing 系統(tǒng)處于臨界狀態(tài)時,加入含噪信號,其變?yōu)椋?/p>
式中:γd和 n(t)與式(2)的意義相同;a 為待測信號幅度。
理論上,含噪信號輸入到處于混沌臨界狀態(tài)的混沌系統(tǒng)中,系統(tǒng)應(yīng)從混沌狀態(tài)迅速轉(zhuǎn)化為大尺度周期狀態(tài),即系統(tǒng)對應(yīng)的最大 Lyapunov 指數(shù)變?yōu)樨摂?shù)。但仿真實驗中發(fā)現(xiàn)在不同精度的臨界閾值下輸入一定信噪比的含噪信號也可能不會引起系統(tǒng)的狀態(tài)變化,即系統(tǒng)對應(yīng)的最大 Lyapunov 指數(shù)仍保持為正數(shù)。
以輸入噪聲為 10–6,輸入信號的信噪比為–10 dB為例,分別取 γd=0.82,0.825,0.825 8,0.825 843 00,每個臨界閾值下進行 100 次仿真實驗,得到噪聲功率為 10–6,信噪比為–10 dB 時,各仿真中系統(tǒng)對應(yīng)的最大 Lyapunov 指數(shù)如圖3 所示。
其中信噪比的定義為:
定義輸入含噪信號時,系統(tǒng)的檢測率為:
由圖3 可得,γd=0.82,0.825,0.825 8,0.825 843 00時,對應(yīng)的檢測率分別為 0%,0%,100%,100%??梢?,隨著臨界閾值的不斷精確,系統(tǒng)對信號的靈敏度增加,這主要是因為隨著臨界閾值的不斷精確,系統(tǒng)越來越接近周期狀態(tài),其有序度增加,微小的信號加入到系統(tǒng)中就可以引起系統(tǒng)狀態(tài)的改變。
圖3 在不同精度臨界閾值下,系統(tǒng)輸入信噪比為–10 dB,噪聲功率為 10–6的含噪信號后,系統(tǒng)對應(yīng)的最大 Lyapunov 指數(shù)Fig. 3 Maximal Lyapunov exponent of chaotic systemafter inputting signal with SNR–10 dB and noise with power 10–6into the systemunder different precision critical value
3.3 不同情況下的檢測性能
取噪聲功率分別為 10–6,10–4和 10–2,信噪比分別為–10 dB,–20 dB,–30 dB,–35 dB 和–40 dB,得到的總檢測結(jié)果如表3 所示。
表3 混沌系統(tǒng)在不同情況下的檢測性能Tab. 3 Detection performance of chaotic systemunder different cases
從表可見,不同精度的臨界值將影響混沌系統(tǒng)的檢測結(jié)果,當選取的臨界值精度越高時,其能檢測到信號的信噪比越低,但其虛警率也越高,因此實際檢測時需根據(jù)檢測的要求來選取合適精度的臨界閾值。從表中還可看到噪聲強度對混沌系統(tǒng)檢測性能也會產(chǎn)生影響,這將在下一階段中詳細研究。
在混沌系統(tǒng)微弱信號檢測中,一般利用 Lyapunov指數(shù)定性確定混沌系統(tǒng)中由混沌態(tài)轉(zhuǎn)化為大尺度周期態(tài)的臨界閾值,而選取具有不同精度的臨界閾值將影響系統(tǒng)的檢測性能,臨界閾值的精度越高,能夠檢測到的信噪比越低,但同時虛警率也越高,因此實際檢測中需根據(jù)對檢測性能的要求選取合適精度的臨界閾值進行信號檢測。
[1]聶春燕. 混沌系統(tǒng)與弱信號檢測[M]. 北京: 清華大學出版社, 2009.
[2]ZHANGm, LIU Z Z, CAO Y, et al. Research and simulation on weak signal detection based on duffing oscillator and damping ratio perturbation[J]. Journal ofmeasurement Science and Instrumentation, 2011, 2(2): 161–163.
[3]張剛, 胡韜, 王穎. 基于Melnikov函數(shù)Duffing系統(tǒng)微弱信號檢測[J]. 電子測量技術(shù), 2015, 38(1): 109–112. ZHANG Gang, HU Tao, WANG Ying. The weak signal detection based on Duffing systemandmelnikov function[J]. Electronicmeasurement Technology, 2015, 38(1): 109–112.
[4]陳軍. 基于混沌理論的檢測系統(tǒng)應(yīng)用研究綜述[J]. 甘肅高師學報, 2013, 18(2): 21–25. CHEN Jun. Application research and discussing of detective systems based on chaos theory[J]. Journal of Gansu Normal Colleges, 2013, 18(2): 21–25.
[5]陳新國, 王潔蕓. 混沌振子在不同初值下檢測弱信號的性能分析[J]. 儀器儀表學報, 2012, 33(12): 2857–2862. CHEN Xin-guo, WANG Jie-yun. Performance analysis of weak signal detection based on chaotic oscillator under different initial condition[J]. Chinese Journal of Scientific Instrument, 2012, 33(12): 2857–2862.
[6]朱斌. 基于混沌理論的微弱信號檢測[J]. 電子科技, 2010, 23(2): 65–67. ZHU Bin. Weak signal detection based on the chaos theory[J]. Electronic Science & Technology, 2010, 23(2): 65–67.
[7]石敏, 徐襲. 基于混沌系統(tǒng)的水下目標輻射噪聲線譜檢測[J].艦船科學技術(shù), 2013, 35(5): 30–33. SHImin, XU Xi. Detection of underwater target radiated noise line-spectra based on chaos system[J]. Ship Science and Technology, 2013, 35(5): 30–33.
[8]朱來普, 張陸勇, 謝文鳳, 等. 基于Duffing混沌振子的微弱信號檢測研究[J]. 無線電工程, 2012, 42(1): 17–20. ZHU Lai-pu, ZHANG Lu-yong, XIE Wen-feng, et al. Research of weak signal detection based on Duffing chaotic oscillator[J]. Radio Engineering, 2012, 42(1): 17–20.
[9]李琳, 劉春剛, 石碩, 等. 基于混沌振子和Lyapunov指數(shù)的微弱信號檢測方法[J]. 黑龍江大學自然科學學報, 2012, 29(4): 556–560. LI Lin, LIU Chun-gang, SHI Shuo, et al. Amethod of weak signal detection based on chaotic oscillator and Lyapunov exponent[J]. Journal of Natural Science of Heilongjiang University, 2012, 29(4): 556–560.
[10]冉莉, 王民. 一種強噪聲背景下微弱信號檢測方法研究[J].信息技術(shù), 2012(2): 41–45. RAN Li, WANGmin. Research on weak signal detectionmethod under the strong noise background[J]. Information Technology, 2012(2): 41–45.
[11]JIN T, ZHANG H. Statistical approach to weak signal detection and estimation using Duffing chaotic oscillators[J]. Science China Information Sciences, 2011, 54(11): 2324–2337.
[12]張賓. Lyapunov特性指數(shù)的算法研究及其在弱信號混沌檢測中的應(yīng)用[D]. 長春: 吉林大學, 2004.
Influence analysis of critical value on detection performance of chaotic system
SHImin1,2(1. Science and Technology on Underwater Acoustic Antagonizing Laboratory, Zhanjiang 524022, China; 2. No. 91388 Unit of PLA, Zhanjiang 524022, China)
Determining the critical value of chaotic systemwhich brings the systemfromchaotic state to large-scale periodic state is very important in detecting weak signal by chaotic system. Simulation shows that choosing critical value of chaotic systemwith different precision will influence the detection performance of chaotic system. The higher precision of critical value of chaotic system, the signal with lower SNR (SNR, signal to noise ratio) can be detected by chaotic system, while the higher false alarmrate is derived. Therefore in actual signal detection, suitable precision of critical value of chaotic systemshould be chosen based on detection requirement in order to bring detection performance achieve the detection requirement.
weak signal detection;detection performance of chaotic system;critical value of chaotic system
TP391
A
1672–7619(2016)10–0137–05
10.3404/j.issn.1672-7619.2016.010.028
2016–03–02;
2016–05–10
石敏(1979–),女,博士,工程師,主要從事水聲信號處理方面的研究。