宮濤 楊建華 單振 王志樂 劉后廣

(中國礦業(yè)大學(xué)機(jī)電工程學(xué)院,江蘇省礦山機(jī)電裝備重點(diǎn)實(shí)驗(yàn)室,徐州 221116)

非線性調(diào)頻(NLFM)信號在雷達(dá)、通信、信號處理中應(yīng)用廣泛,該類信號所激勵(lì)下的非線性系統(tǒng)響應(yīng)有著豐富的信息,通過共振來增強(qiáng)NLFM 信號具有一定實(shí)際意義與價(jià)值.本文主要研究了受到不同類型NLFM信號所激勵(lì)的非線性系統(tǒng)共振現(xiàn)象,提出了實(shí)時(shí)尺度變換方法來處理高頻NLFM 信號,克服人為選擇造成信號輸出響應(yīng)較差的缺點(diǎn).同時(shí),提出實(shí)時(shí)譜放大因子作為共振評價(jià)指標(biāo),準(zhǔn)確評價(jià)NLFM 信號激勵(lì)下的系統(tǒng)共振響應(yīng),討論系統(tǒng)參數(shù)對系統(tǒng)最優(yōu)共振響應(yīng)的影響,參數(shù)選擇在合理的區(qū)間內(nèi)即可實(shí)現(xiàn)最優(yōu)共振響應(yīng),不僅實(shí)現(xiàn)信號特征的大幅增強(qiáng),還保持信號時(shí)頻特征的連續(xù)性.最后,將實(shí)時(shí)尺度變換方法和固定尺度變換方法進(jìn)行對比,說明所提方法處理NLFM 信號的優(yōu)越性.

1 引言

共振是系統(tǒng)響應(yīng)中一種典型的動(dòng)力學(xué)現(xiàn)象.通常來講,對于簡諧激勵(lì),當(dāng)外激勵(lì)頻率接近于振動(dòng)系統(tǒng)固有頻率時(shí)會(huì)發(fā)生共振現(xiàn)象,此時(shí)系統(tǒng)振動(dòng)振幅會(huì)急劇增強(qiáng).不同于傳統(tǒng)意義的共振,隨機(jī)共振(stochastic resonance,SR)[1]、振動(dòng)共振(vibrational resonance,VR)[2]得到廣泛研究,利用其獨(dú)特的動(dòng)力學(xué)特性可以增強(qiáng)微弱的信號.其中隨機(jī)共振是在微弱低頻信號、噪聲和非線性系統(tǒng)協(xié)同作用下使得系統(tǒng)輸出得到顯著增強(qiáng)的一種動(dòng)力學(xué)現(xiàn)象,最初由意大利學(xué)者Benzi 提出[3],之后在信號處理中得到廣泛應(yīng)用[4-6].隨機(jī)共振發(fā)生的基本條件是輸入的微弱信號滿足小參數(shù)條件,即信號頻率遠(yuǎn)小于1.而實(shí)際中的信號大部分不滿足小參數(shù)條件,為了實(shí)現(xiàn)大參數(shù)信號的隨機(jī)共振,一些尺度變換方法被引入,比如二次采樣[7]、頻移尺度變換[8]、歸一化尺度變換[9]、普通尺度變換[10].振動(dòng)共振與隨機(jī)共振類似,區(qū)別是輸入的噪聲被高頻輔助信號代替,尺度變換方法與隨機(jī)共振的方法一致.而單個(gè)信號激勵(lì)下的系統(tǒng)共振[11]是指不需要輔助信號,依靠系統(tǒng)參數(shù)變化實(shí)現(xiàn)系統(tǒng)輸出增強(qiáng)的一種現(xiàn)象,該現(xiàn)象有一定實(shí)際意義與價(jià)值,在信號增強(qiáng)方面,無需額外的高頻輔助信號、噪聲即可實(shí)現(xiàn)信號的增強(qiáng),操作方便,實(shí)現(xiàn)過程簡單.

在信號處理領(lǐng)域,利用共振現(xiàn)象來增強(qiáng)不同類型的信號具有獨(dú)特的優(yōu)勢,相較于其他信號處理算法,如經(jīng)驗(yàn)?zāi)B(tài)分解[12]、變分模態(tài)分解[13]、小波變換[14]、信號濾波[15]等方法,該類方法不僅可以增強(qiáng)信號,而且可以有效抑制干擾噪聲.因此,基于共振的方法在信號處理中應(yīng)用廣泛.許多學(xué)者利用最優(yōu)共振現(xiàn)象來增強(qiáng)和處理各類信號,比如簡諧信號[16]、非周期二進(jìn)制信號[17]、M-ary 信號[18]、頻率調(diào)制(frequency modulation,FM)信號[19,20]等.其中FM 信號在雷達(dá)[21]、通信[22]、信號處理[23]中應(yīng)用廣泛,該類信號的特點(diǎn)是頻率受到調(diào)制,按照其瞬時(shí)頻率變化的規(guī)律主要分為線性調(diào)頻(linear frequency modulation,LFM)信號和非線性調(diào)頻(nonlinear frequency modulation,NLFM)信號,目前已有一些LFM 信號激勵(lì)下的非線性系統(tǒng)動(dòng)力學(xué)行為的研究,但是針對NLFM 信號的相關(guān)研究很少,考慮到NLFM 信號[24-26]廣泛存在,如何實(shí)現(xiàn)NLFM 信號最優(yōu)共振增強(qiáng)效果及NLFM 信號激勵(lì)下的非線性系統(tǒng)最優(yōu)共振響應(yīng)的研究是一個(gè)有實(shí)際意義與價(jià)值的問題.

目前關(guān)于系統(tǒng)共振的研究,主要針對周期信號,其信號頻率固定,通過尺度變換方法即可滿足共振的小參數(shù)要求,調(diào)整合適的系統(tǒng)參數(shù)可實(shí)現(xiàn)最優(yōu)共振響應(yīng).但是對于復(fù)雜NLFM 信號,其頻率隨時(shí)間非線性變化,此類信號的頻率變化復(fù)雜,信號在頻譜圖中表現(xiàn)為連續(xù)頻帶,而且之前所用的尺度變換方法沒有考慮頻率的變化性,采用傅里葉變換方法無法準(zhǔn)確描述信號特征,需要借助一些時(shí)頻分析方法來描述.因此,針對NLFM 信號激勵(lì)下的最優(yōu)共振響應(yīng)研究具有一定的難度與復(fù)雜性,采用的技術(shù)較周期信號的研究也更為復(fù)雜,采用固定尺度系數(shù)的尺度變換方法來處理NLFM 信號顯然不能滿足要求,無法實(shí)現(xiàn)最優(yōu)參數(shù)匹配,系統(tǒng)輸出響應(yīng)效果較差.同時(shí),NLFM 信號激勵(lì)下的系統(tǒng)共振響應(yīng)有一定的理論研究價(jià)值,對于實(shí)現(xiàn)NLFM 信號增強(qiáng)有一定的意義.基于此,本文研究了NLFM 信號激勵(lì)下的非線性系統(tǒng)最優(yōu)共振響應(yīng),提出了一種時(shí)變尺度系統(tǒng)共振方法有效實(shí)現(xiàn)復(fù)雜NLFM 信號的增強(qiáng).該方法的基本原理是利用實(shí)時(shí)變化的尺度系數(shù)來匹配瞬時(shí)頻率變化的NLFM 信號和系統(tǒng)參數(shù),以盡量保證在每一時(shí)刻信號的瞬時(shí)頻率滿足最優(yōu)小參數(shù)條件,通過調(diào)整系統(tǒng)參數(shù)和尺度系數(shù)實(shí)現(xiàn)最優(yōu)共振響應(yīng).

本文提出了時(shí)變尺度系統(tǒng)共振理論,同時(shí)定義了新的共振評價(jià)指標(biāo)-實(shí)時(shí)譜放大因子,可有效評價(jià)NLFM 信號的共振增強(qiáng)效果;處理了兩類典型的NLFM 信號,研究了相關(guān)參數(shù)對最優(yōu)共振響應(yīng)的影響,在保證信號時(shí)頻特征連續(xù)性的同時(shí),實(shí)現(xiàn)了信號的共振增強(qiáng),與固定尺度變換方法相比,本文所提方法克服了尺度系數(shù)選擇人為性造成信號增強(qiáng)效果較差的缺點(diǎn).

2 時(shí)變尺度系統(tǒng)共振理論

一種時(shí)變雙穩(wěn)態(tài)系統(tǒng)共振模型的動(dòng)力學(xué)方程表示如下:

其中a(t),b(t)分別為時(shí)變系統(tǒng)參數(shù) (a(t),b(t)?1),系統(tǒng)為雙穩(wěn)態(tài)系統(tǒng).x(t) 為系統(tǒng)的輸出,s(t) 為輸入的NLFM 信號.NLFM 信號s(t) 的瞬時(shí)頻率f(t)隨時(shí)間不斷變化,為實(shí)現(xiàn)尺度系數(shù)與信號s(t)及系統(tǒng)參數(shù)a(t),b(t)相匹配,引入時(shí)變尺度系數(shù)β(t),表達(dá)式為

其中β0為尺度系數(shù)常數(shù),f1(t)為s(t)的瞬時(shí)頻率f(t)的離散化表示,β(t) 滿足條件β(t)?1.同時(shí),系統(tǒng)參數(shù)a(t),b(t)與時(shí)變尺度系數(shù)β(t)滿足以下關(guān)系:

其中a1?1,b1?1,(5)式滿足系統(tǒng)共振小參數(shù)條件.(5)式相較(1)式,NLFM 信號s(t)幅值縮小到1/β(t),將(5)式的NLFM 信號幅值乘以β(t),變?yōu)?/p>

若系統(tǒng)共振發(fā)生在(6)式中,則系統(tǒng)共振也發(fā)生在(7)式中:

(7)式為一個(gè)大參數(shù)共振模型,可檢測任意大頻率的信號.

譜放大因子是經(jīng)典的共振性能評價(jià)指標(biāo),針對NLFM 信號,信號的瞬時(shí)頻率非線性變化,利用譜放大因子難以刻畫NLFM 信號的共振效果.因此,本文定義了一種實(shí)時(shí)譜放大因子η來刻畫NLFM信號激勵(lì)下的系統(tǒng)共振效果,如下式:

其中η為系統(tǒng)輸出的特征頻帶的放大程度,N為總數(shù)據(jù)點(diǎn)數(shù),M為頻率曲線f(t) 兩側(cè)的頻帶寬度m對應(yīng)數(shù)據(jù)點(diǎn)數(shù),X(·) 和S(·) 分別為輸入信號和輸出信號的離散短時(shí)傅里葉變換譜值,p(f(ti)) 為ti時(shí)刻信號的瞬時(shí)頻率處點(diǎn)的位置.核心是在STFT 頻譜上定義譜放大因子,將輸出信號STFT 頻譜圖中信號時(shí)頻一定范圍內(nèi)的能量之和與輸入NLFM 信號的時(shí)頻圖中脊線能量之比作為評價(jià)指標(biāo),綜合考慮NLFM 信號的時(shí)頻特征,描述系統(tǒng)的共振響應(yīng)更準(zhǔn)確、可靠.若η＞ 1,說明NLFM 信號被有效增強(qiáng).

3 幾種不同NLFM 信號增強(qiáng)仿真分析

3.1 分段NLFM 信號

分段NLFM 信號較為典型,在機(jī)械設(shè)備信號處理中十分常見,由軸承、齒輪等零件經(jīng)過復(fù)雜路徑產(chǎn)生的信號為FM 信號,按照轉(zhuǎn)速變化的特征,會(huì)產(chǎn)生反映零件運(yùn)行狀態(tài)的LFM 信號和NLFM信號.其中與分段NLFM 信號頻率特征相對應(yīng)的信號也很常見,比如文獻(xiàn)[27,28]中這類轉(zhuǎn)速工況,礦井提升機(jī)運(yùn)行的典型轉(zhuǎn)速工況為速度曲線線性上升-平穩(wěn)-線性下降,反映出提升過程是線性加速-穩(wěn)速-線性減速,這樣的轉(zhuǎn)速是人為控制的,關(guān)鍵部位的軸承故障特征信號均為此類分段NLFM 信號;流水線上的工業(yè)機(jī)器人的零件運(yùn)行工況所產(chǎn)生的信號也為此類NLFM 信號,在故障診斷領(lǐng)域該類信號很常見,應(yīng)用廣泛,此類信號具有一定的代表意義.

一種典型的分段NLFM 信號s1(t)的頻率f(t)隨時(shí)間變化的分段函數(shù)為

其中f0為起始頻率,γ為頻率變化率,t1,t2分別為頻率規(guī)律發(fā)生變化的時(shí)刻.對應(yīng)余弦函數(shù)的NLFM信號s1(t)描述為

其中A為信號幅值,φ為初始相位,s1(t)的頻率變化規(guī)律為線性上升-不變-線性下降.

圖1 給出了該NLFM 信號的相關(guān)信息,圖1(a),(b)分別為信號的時(shí)域波形圖和STFT 頻譜圖,在STFT 頻譜圖中可以觀察到信號的頻率特征曲線在0—3 s 內(nèi)從0 Hz 增大到150 Hz,然后保持至7 s,最后降到0 Hz.

針對圖1 中的信號,圖2 分別給出了不同尺度系數(shù)下的系統(tǒng)輸出實(shí)時(shí)譜放大因子η隨系統(tǒng)參數(shù)a1的變化曲線,選取的固定尺度系數(shù)有β=β0×10,β0× 75,β0× 150,時(shí)變尺度系數(shù)有β=β0×f1(t),在固定尺度系數(shù)方法下的系統(tǒng)輸出實(shí)時(shí)譜放大因子η隨著尺度系數(shù)的增大而增強(qiáng),在最大值處實(shí)現(xiàn)了NLFM 信號最優(yōu)系統(tǒng)共振響應(yīng).同時(shí),在時(shí)變尺度系數(shù)方法下的系統(tǒng)輸出實(shí)時(shí)譜放大因子最大值ηmax比β=β0× 150 時(shí)的最大值略小,但是區(qū)別不大.

圖1 輸入的NLFM 信號 (a)時(shí)域波形圖;(b) STFT 頻譜圖 (仿真參數(shù)A=0.1,f0=0,γ=50,φ=0,t1=3,t2=10)Fig.1.The input NLFM signal:(a) The time domain waveform;(b) the STFT spectrum.The simulation parameters are A=0.1,f0=0,γ=50,φ=0,t1=3,t2=10.

圖2 不同尺度系數(shù)β=β0 × 10,β0 × 75,β0 × 150,β=β0 ×f1(t) 下的系統(tǒng)參數(shù)a1 與實(shí)時(shí)譜放大因子η 的關(guān)系圖(仿真參數(shù)為b1=1,β0=100,m=2)Fig.2.Curves of η-a1 are obtained under different scale coefficient β=β0 × 10,β0 × 75,β0 × 150,β=β0 ×f1(t).The simulation parameters are b1=1,β0=100,m=2.

為了探究系統(tǒng)參數(shù)對最優(yōu)共振響應(yīng)的影響,圖3 給出了實(shí)時(shí)譜放大因子η與系統(tǒng)參數(shù)a1,b1的偽三維關(guān)系圖.其中η的大小表示共振響應(yīng)的強(qiáng)弱程度,η值越大表示共振響應(yīng)越強(qiáng),在圖形中可以發(fā)現(xiàn)明顯的共振區(qū)域,在這個(gè)共振區(qū)域內(nèi)均可發(fā)生系統(tǒng)最優(yōu)共振響應(yīng),實(shí)現(xiàn)NLFM 信號增強(qiáng).

圖3 實(shí)時(shí)譜放大因子η 與系統(tǒng)參數(shù)a1,b1 的關(guān)系圖 (仿真參數(shù)為m=2,β0=100)Fig.3.Contour plot of the real time spectral amplification factor η in the a1 -b1 plane.The simulation parameters are m=2,β0=100.

在圖2 中,當(dāng)a1=0.32,b1=1 時(shí),時(shí)變尺度方法對應(yīng)的譜放大因子η取得最大值ηmax=6.4,與其對應(yīng)的最優(yōu)系統(tǒng)共振響應(yīng)如圖4 所示.圖4(a),(b)分別為響應(yīng)的時(shí)域波形圖和STFT 頻譜圖,信號幅值得到明顯增強(qiáng)的同時(shí)保證了信號的連續(xù)性.

圖4 b1=1 時(shí)系統(tǒng)的共振輸出結(jié)果 (a) 時(shí)域波形圖;(b) STFT 頻譜圖(仿真參數(shù)為m=2,β0=100)Fig.4.The output response of system resonance when b1=1:(a) The time domain waveform;(b) the STFT spectrum.The simulation parameters are m=2,β0=100.

為說明時(shí)變尺度方法的必要性,將固定尺度系數(shù)變尺度方法和時(shí)變尺度系數(shù)變尺度方法的最優(yōu)輸出響應(yīng)進(jìn)行對比.本文選擇固定尺度系數(shù)為β=β0× 10,β0× 75,β0× 150,目的是模擬尺度系數(shù)較小、適中、較大的情況,尺度變換后信號的頻率變?yōu)閒(t)/β,在圖1 中f(t)的范圍為0 ≤f(t) ≤ 150,采用固定尺度系數(shù)變換后頻率范圍分別為0 ≤f(t)/(β0×10) ≤ (15/β0),0 ≤f(t)/(β0×75) ≤ (2/β0),0 ≤f(t)/(β0×150) ≤ (1/β0),尺度系數(shù)基數(shù)β0固定(文中取100),可以發(fā)現(xiàn)固定尺度系數(shù)為β0×10時(shí),尺度變換后頻率達(dá)到0.15,此時(shí)雖然滿足系統(tǒng)共振小參數(shù)要求,但是沒有達(dá)到最優(yōu)輸出效果(見圖5),說明尺度變換后的頻率需要進(jìn)一步降低;當(dāng)尺度系數(shù)逐漸增大,尺度變換后的頻率逐漸減小,此時(shí)針對整個(gè)信號均可以滿足共振的小參數(shù)要求,實(shí)現(xiàn)最優(yōu)系統(tǒng)共振響應(yīng).

圖5 不同的固定尺度系數(shù)β 下的系統(tǒng)共振輸出結(jié)果 (a) β=β0 × 10 下系統(tǒng)輸出STFT 譜圖;(b) β=β0 ×10 下系統(tǒng)輸出時(shí)域波形圖;(c) β=β0 × 75 下系統(tǒng)輸出STFT 譜圖;(d) β=β0 × 75 下時(shí)域波形圖;(e) β=β0 ×150 下系統(tǒng)輸出STFT 譜圖;(f) β=β0 ×150 下時(shí)域波形圖 (仿真參數(shù)為β0=100)Fig.5.The output response of system resonance under different fixed scale coefficients β:(a) The output STFT spectrum under β=β0 × 10;(b) the output time domain waveform under β=β0 × 10;(c) the output STFT spectrum under β=β0 × 75;(d) the output time domain waveform under β=β0 × 75;(e) the output STFT spectrum under β=β0 × 150;(f) the output time domain waveform under β=β0 × 150.The simulation parameters are β0=100.

圖5 為不同固定尺度系數(shù)β下的最優(yōu)系統(tǒng)共振響應(yīng)結(jié)果,當(dāng)β=β0×10 時(shí),系統(tǒng)最優(yōu)輸出響應(yīng)變差,在圖5(a)中的信號幅值增強(qiáng)效果不顯著,圖5(b) 的時(shí)域波形圖中信號波形產(chǎn)生很大失真.當(dāng)β分別為β0× 75,β0× 150 時(shí),系統(tǒng)共振達(dá)到最優(yōu)輸出響應(yīng),與時(shí)變尺度方法下的最優(yōu)輸出響應(yīng)區(qū)別不大.但是,固定尺度系數(shù)下系統(tǒng)最優(yōu)響應(yīng)依賴于人為選擇的尺度系數(shù),一旦參數(shù)選擇不理想,共振輸出響應(yīng)會(huì)變得很差,而時(shí)變尺度方法將系數(shù)定為NLFM 信號的離散化頻率f1(t),克服了人為選擇造成的系數(shù)最優(yōu)輸出響應(yīng)不穩(wěn)定的缺點(diǎn).

3.2 調(diào)幅-調(diào)頻NLFM 信號

在信號處理中,由正弦信號頻率調(diào)制和幅值調(diào)制的信號也是一種典型的NLFM 信號,該類信號特點(diǎn)是瞬時(shí)頻率隨時(shí)間波動(dòng),變化規(guī)律復(fù)雜,研究這類信號最優(yōu)系統(tǒng)共振響應(yīng)及增強(qiáng)具有重要意義.因此,根據(jù)文獻(xiàn)[29],本文構(gòu)造了一種調(diào)幅-調(diào)頻的NLFM 仿真信號,如 (11) 式所示:

該仿真信號受到正弦幅值、雙正弦頻率調(diào)制,如圖6所示.圖6(a),(b) 分別為該信號的時(shí)域波形圖和STFT 頻譜圖,不同時(shí)刻信號的幅值發(fā)生變化,頻率也發(fā)生變化,頻率分布在10—75 Hz.黑色實(shí)線為NLFM 信號通過Hilbert 變換相位提取的瞬時(shí)頻率曲線,與時(shí)頻特征曲線相一致.對于此類NLFM信號而言,固定尺度變換方法不適用.

圖6 輸入的NLFM 信號 (a) 時(shí)域波形圖;(b) STFT 時(shí)頻圖Fig.6.The input NLFM signal:(a) The time domain waveform;(b) the STFT spectrum.

圖7 為a1=0.01 時(shí)不同尺度系數(shù)下實(shí)時(shí)譜放大因子η與系統(tǒng)參數(shù)b1的關(guān)系圖,在不同的尺度系數(shù)下,隨著b1的增大,η表現(xiàn)出先增大后減小的規(guī)律.對于時(shí)變尺度系數(shù)β=β0×f1(t),當(dāng)b1=0.0001,a1=0.01 時(shí),實(shí)時(shí)譜放大因子最大為101.1.對于固定尺度系數(shù)的情況,僅從數(shù)值上來看,其譜放大因子既可能比對應(yīng)的實(shí)時(shí)尺度系數(shù)下共振輸出結(jié)果小,也可能比其大.

圖7 不同尺度系數(shù)下的實(shí)時(shí)譜放大因子η 與系統(tǒng)參數(shù)b1 關(guān)系圖 (仿真參數(shù)為a1=0.01,m=5,β0=1000)Fig.7.Curves of η-b1 are obtained under different scale coefficients.The simulation parameters are a1=0.01,m=5,β0=1000.

在時(shí)變尺度方法處理下,其對應(yīng)的最優(yōu)輸出信號如圖8 所示.圖8(a),(b) 分別為最優(yōu)系統(tǒng)共振響應(yīng)時(shí)域波形圖和STFT 頻譜圖,信號幅值大幅增強(qiáng),同時(shí)保證信號時(shí)頻特征的連續(xù)性.需要說明的是在圖8(b)中短時(shí)傅里葉譜能量計(jì)算范圍,黑色實(shí)線包圍區(qū)間表示響應(yīng)信號的能量計(jì)算區(qū)間,該區(qū)間內(nèi)所有STFT 能量值之和作為輸出譜能量,輸入信號譜能量為NLFM 信號頻率曲線上STFT能量值之和.

圖8 時(shí)變尺度系數(shù)下最優(yōu)共振輸出結(jié)果 (a) 時(shí)域波形圖;(b) STFT 頻譜圖 (仿真參數(shù)為m=5,β0=1000)Fig.8.The output response of optimal resonance:(a) The time domain waveform;(b) the STFT spectrum.The simulation parameters are m=5,β0=1000.

雖然在圖7 實(shí)時(shí)譜放大因子η與系統(tǒng)參數(shù)b1的關(guān)系圖中,固定的尺度系數(shù)所導(dǎo)致的共振結(jié)果η(β=β0×40,β0×75)有的比時(shí)變尺度系數(shù)共振結(jié)果的η大,最優(yōu)共振結(jié)果如圖9 所示,但是固定尺度變換方法處理后的NLFM 信號的連續(xù)性很差,能量發(fā)散現(xiàn)象嚴(yán)重,遠(yuǎn)沒有時(shí)變尺度共振效果好.

圖9 不同固定尺度系數(shù)β 下共振響應(yīng) (a) β=β0 × 10 下輸出STFT 譜圖;(b) β=β0 × 10 下輸出時(shí)域波形圖;(c) β=β0 ×40 下輸出STFT 譜圖;(d) β=β0 × 40 下輸出時(shí)域波形圖;(e) β=β0 × 75 下輸出STFT 譜圖;(f) β=β0 × 75 下輸出時(shí)域波形圖Fig.9.The resonance response under different fixed scale coefficient β:(a) The output STFT spectrum under β=β0 × 10;(b) the output time domain waveform under β=β0 × 10;(c) the output STFT spectrum under β=β0 × 40;(d) the output time domain waveform under β=β0 × 40;(e) the output STFT spectrum under β=β0 × 75;(f) the output time domain waveform under β=β0 × 75.

4 結(jié)論

1) 針對傳統(tǒng)的尺度變換方法不能匹配NLFM信號的時(shí)變頻率特征的問題,提出了隨信號瞬時(shí)頻率變化而變化的時(shí)變尺度共振方法,以實(shí)現(xiàn)每一時(shí)刻的最優(yōu)共振響應(yīng),進(jìn)而增強(qiáng)NLFM 信號.

2) 針對傳統(tǒng)的共振效果評價(jià)指標(biāo)-譜放大因子在NLFM 信號激勵(lì)下的系統(tǒng)共振響應(yīng)評價(jià)方面存在一定局限性,定義了新的指標(biāo)-實(shí)時(shí)譜放大因子,可準(zhǔn)確評價(jià)NLFM 信號激勵(lì)下系統(tǒng)共振響應(yīng),同時(shí)保證了系統(tǒng)輸出信號時(shí)頻特征連續(xù)性.

3) 研究了非線性系統(tǒng)參數(shù)對共振輸出響應(yīng)的影響,指出最優(yōu)共振響應(yīng)發(fā)生在明顯的共振帶上,選擇位于共振帶的系統(tǒng)參數(shù)即可實(shí)現(xiàn)最優(yōu)非線性系統(tǒng)共振輸出,實(shí)現(xiàn)對NLFM 信號的增強(qiáng).

本文實(shí)現(xiàn)了NLFM 信號的最優(yōu)非線性系統(tǒng)共振響應(yīng),為NLFM 信號增強(qiáng)及檢測提供參考,在信號增強(qiáng)方面有一定的實(shí)際意義與價(jià)值.同時(shí),給出了NLFM 信號激勵(lì)下系統(tǒng)最優(yōu)輸出響應(yīng)的相關(guān)影響規(guī)律,對研究復(fù)雜信號激勵(lì)下的系統(tǒng)動(dòng)力學(xué)行為有一定的參考價(jià)值.