[發(fā)明專(zhuān)利]基于硬判決的分布式真假?gòu)楊^識(shí)別方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201510230515.5
• 申請(qǐng)日: 2015-05-07
• 公開(kāi)（公告）號(hào): CN105259545B
• 公開(kāi)（公告）日: 2018-04-06
• 發(fā)明（設(shè)計(jì)）人: 張凱麗;何茜;何子述
• 申請(qǐng)（專(zhuān)利權(quán)）人: 電子科技大學(xué)
• 主分類(lèi)號(hào): G01S7/41
• 分類(lèi)號(hào): G01S7/41
• 代理公司: 電子科技大學(xué)專(zhuān)利中心51203
• 代理人: 張楊

1.一種基于硬判決的分布式真假?gòu)楊^識(shí)別方法，其特征在于包括以下步驟：

步驟1把MIMO雷達(dá)接收陣元?jiǎng)澐譃镚個(gè)接收子塊，各個(gè)子塊中的接收天線進(jìn)行聯(lián)合觀測(cè)，獲得k時(shí)刻局部觀測(cè)向量

步驟2計(jì)算如下所示的基于分布式處理的狀態(tài)空間模型中的參數(shù)：

x_k+1＝F_k+1x_k+G_k+1u_k+w_k+1

zkg=Hxk+ekg,g=1,...,G]]>

步驟2-1根據(jù)工程經(jīng)驗(yàn)和MIMO雷達(dá)估計(jì)性能得到狀態(tài)噪聲w_k的均值和協(xié)方差矩陣Σ_w；

步驟2-2根據(jù)彈頭目標(biāo)運(yùn)動(dòng)分析，求得系統(tǒng)輸入向量u_k；

步驟2-3由系統(tǒng)預(yù)處理得到F₀與G_k；

步驟2-4依據(jù)觀測(cè)向量確定觀測(cè)矩陣H；

步驟2-5計(jì)算出局部觀測(cè)噪聲的均值以及協(xié)方差矩陣

步驟2-5-1局部觀測(cè)噪聲建模為高斯白噪聲；

步驟2-5-2設(shè)定均值為0；

步驟2-5-3設(shè)定與接收子塊g有關(guān)，而與觀測(cè)時(shí)刻k無(wú)關(guān)，其對(duì)角線元素由狀態(tài)估計(jì)方差CRB計(jì)算得到；

步驟2-6根據(jù)彈頭目標(biāo)初始狀態(tài)向量x₀、MIMO雷達(dá)發(fā)射波形以及發(fā)射、接收天線位置，計(jì)算出步驟2-5-3中待求的狀態(tài)估計(jì)方差CRB；

步驟3計(jì)算局部LOUD檢測(cè)器，

其中n表示LOUD檢測(cè)器的判決時(shí)刻，L′是狀態(tài)轉(zhuǎn)移矩陣F_k+1的維數(shù)；H₀和H_l分別對(duì)應(yīng)真彈頭、假?gòu)楊^兩個(gè)假設(shè)；

步驟3-1給定卡爾曼濾波初值x_0|0和

步驟3-2計(jì)算出卡爾曼濾波過(guò)程的一步迭代結(jié)果；

步驟3-2-1 x_k+1|k＝F₀x_k|k+G_k+1u_k；

步驟3-2-2

步驟3-2-3

步驟3-2-4

步驟3-3若k＜n-2，令k＝k+1，重復(fù)步驟3-2；若k＝n-2，則停止迭代；

步驟3-4迭代停止后得到x_n-1|n-1與

步驟3-5計(jì)算條件概率密度函數(shù)

步驟3-5-1令F_n＝F₀；

步驟3-5-2 x_n|n-1＝F_nx_n-1|n-1+G_nu_n-1；

步驟3-5-3

步驟3-5-4求出的條件均值

步驟3-5-5求出的條件協(xié)方差矩陣

步驟3-5-6由求得其概率密度函數(shù)：

p(zng|Zn-1g;Fn)=exp{-12(zng-zn|n-1g)T(Σn|n-1g)-1(zng-zn|n-1g)}(2π)L2|Σn|n-1g|12]]>

步驟3-6計(jì)算檢測(cè)器分子部分的求和項(xiàng)

步驟3-6-1令i＝1,j＝1；

步驟3-6-2令F_n＝F₀+ΔFδ(i,j)，其中δ(i,j)是除了第i行第j列的元素為1其他元素均為0的L'維方陣；

步驟3-6-3利用步驟3-5-2至3-5-6求出

步驟3-6-4令F_n＝F₀-ΔFδ(i,j)；

步驟3-6-5利用步驟3-5-2至3-5-6求出

步驟3-6-6用偏導(dǎo)定義求出檢測(cè)器分子部分的求和項(xiàng)

Πng(i,j)=∂2p(zng|Zn-1g;Fn)∂Fn2(i,j)|Fn=F0=P(zng|Zn-1g;F0+ΔFδ(i,j))+p(zng|Zn-1g;F0-ΔFδ(i,j))-2p(zng|Zn-1g;F0)(ΔF)2]]>

步驟3-6-7若j＜L'，令j＝j(luò)+1，重復(fù)步驟3-6-2至3-6-6；若j＝L'且i＜L'，令i＝i+1，令j＝1，重復(fù)步驟3-6-2至3-6-6；若j＝L'且i＝L'，進(jìn)入下一步驟；

步驟3-7計(jì)算出局部處理器的LOUD檢驗(yàn)統(tǒng)計(jì)量，

ΓLOUDg(zng)=Σj=1L′Σi=1L′Πng(i,j)p(zng|Zn-1g;F0)]]>

步驟4求得虛警概率P_fa＝α下的檢測(cè)門(mén)限

步驟5將局部處理器的檢驗(yàn)統(tǒng)計(jì)量與門(mén)限比較，若則判定為假?gòu)楊^，即局部判決結(jié)果u_g為1，反之u_g為0；

步驟6重復(fù)利用步驟3～5得到G個(gè)局部處理器的判決結(jié)果u_g,g＝1,...,G；

步驟7融合中心收集到各個(gè)局部判決結(jié)果，利用K/G準(zhǔn)則進(jìn)行全局判決，

步驟8若全局判決u_f＝1，則判定為假?gòu)楊^；若u_f＝0，判定為真彈頭。