[發(fā)明專利]一種基于杜芬振子逆向相變的弱信號(hào)檢測(cè)方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201710251003.6
• 申請(qǐng)日: 2017-04-18
• 公開（公告）號(hào): CN107066740B
• 公開（公告）日: 2020-04-14
• 發(fā)明（設(shè)計(jì)）人: 劉海波;姜慧;唐波;譚堅(jiān)文;馬大瑋;宋傳座;趙有華;張紀(jì)忠
• 申請(qǐng)（專利權(quán)）人: 中國(guó)人民解放軍重慶通信學(xué)院;重慶大學(xué)
• 主分類號(hào): G06F30/20
• 分類號(hào): G06F30/20
• 代理公司: 北京天奇智新知識(shí)產(chǎn)權(quán)代理有限公司 11340
• 代理人: 馬冬新

1.一種基于杜芬振子逆向相變的弱信號(hào)檢測(cè)方法，其特征在于：所述檢測(cè)方法包括：

(1)構(gòu)建檢測(cè)陣列，所述檢測(cè)陣列包括兩組1×N₁個(gè)杜芬振子，第1組第n個(gè)杜芬振子和第2組的第n個(gè)杜芬振子初始相軌跡點(diǎn)、參考信號(hào)初始相位均相同，調(diào)整兩組1×N₁檢測(cè)振子的參考信號(hào)頻率等于待測(cè)信號(hào)頻率，使得1×N₁個(gè)杜芬振子參考信號(hào)初始相位等間隔分布于[0,2π]，所述間隔為2π/N₁，所述檢測(cè)陣列中第(m,n)杜芬振子的檢測(cè)方程為：

(2)根據(jù)步驟(1)中檢測(cè)陣列，輸入待測(cè)信號(hào)為γ_m＞＞εr(t)，將第一組振子參考信號(hào)及第二組振子參考信號(hào)與待測(cè)信號(hào)中的小信號(hào)矢量合成，調(diào)整第一組杜芬振子的參考信號(hào)幅度γ₁及第二組杜芬振子的參考信號(hào)幅度γ₂，使第一組N₁個(gè)杜芬振子中N₂個(gè)杜芬振子發(fā)生逆向相變，第二組N₁個(gè)杜芬振子中N₁-N₂'個(gè)杜芬振子發(fā)生逆向相變，根據(jù)矢量合成關(guān)系，計(jì)算得待測(cè)信號(hào)幅度值為：

(3)重復(fù)步驟(2)，測(cè)量L次，第l次的N₂和N₂'的記錄值分別為N₂^(l)和N₂'^(l)，待測(cè)信號(hào)幅度值的第l次估計(jì)結(jié)果為：

對(duì)待測(cè)信號(hào)幅度值測(cè)量結(jié)果求平均值，所述平均值為待測(cè)信號(hào)幅度值：

其中，m為振子陣列組號(hào)，1≤m≤2，n為每組振子的編號(hào)，1≤n≤N₁；N₁為每組振子包含的振子數(shù)，且N₁≥3；N₂為第一組振子中發(fā)生相變的振子數(shù)，N₂＜N₁/2；N₂'為第二組振子中未發(fā)生逆向相變的振子數(shù)，N₂'＜N₁/2；Ω為參考信號(hào)的角頻率；Φ_n＝2π(n-1)/N₁為第n個(gè)振子的參考信號(hào)初始相位；γ_m為第m組振子參考信號(hào)幅度值，ε為增益控制參數(shù)；n(t)為輸入待測(cè)信號(hào)中的噪聲；L為獨(dú)立測(cè)量的次數(shù)，l為獨(dú)立測(cè)量次數(shù)編號(hào)，1≤l≤L，第l次測(cè)量相應(yīng)參數(shù)分別記為：N₂記為N₂^(l)，N₂'為記N₂'^(l)，γ_m記為γ_m^(l)。

2.根據(jù)權(quán)利要求1所述的基于杜芬振子逆向相變的弱信號(hào)檢測(cè)方法，其特征在于：所述N₂＝N₂'，l次測(cè)量中N₂^(l)和N₂'^(l)統(tǒng)計(jì)相等。

3.根據(jù)權(quán)利要求1所述的基于杜芬振子逆向相變的弱信號(hào)檢測(cè)方法，其特征在于：所述步驟(2)還包括設(shè)置振子初始狀態(tài)為穩(wěn)定大周期態(tài)，包括利用數(shù)值仿真方法確定穩(wěn)定大周期態(tài)下步驟(1)所述杜芬振子相軌跡點(diǎn)(x,y)和系統(tǒng)策動(dòng)力相位Φ的對(duì)應(yīng)關(guān)系，根據(jù)所述對(duì)應(yīng)關(guān)系設(shè)置杜芬振子的初始相軌跡點(diǎn)位置(x_0n,y_0n)和參考信號(hào)初始相位Φ_n；

其中，(x_0n,y_0n)和Φ_n分別為第n個(gè)振子的初始相軌跡點(diǎn)和初始相位。