陣，同時(shí)構(gòu)建與該參數(shù)矩陣同型的零矩陣作為第一累加矩陣；1η、1 12、φι參數(shù)矩陣中的每一 參數(shù)數(shù)組與第一累加矩陣中每一行對(duì)應(yīng);a為1ιι、1?2、φι參數(shù)矩陣中參數(shù)數(shù)組的數(shù)量；
[0092] 并且，將121、122、炳參數(shù)數(shù)組中的每一數(shù)組逐行排列，生成13*3的121、122、9 2參數(shù)矩 陣，同時(shí)構(gòu)建與該參數(shù)矩陣同型的零矩陣作為第二累加矩陣；121、122、φ2參數(shù)矩陣中的每一 參數(shù)數(shù)組與第二累加矩陣中每一行對(duì)應(yīng);b為1 21、122、參數(shù)矩陣中參數(shù)數(shù)組的數(shù)量。
[0093] S235.將第一雷達(dá)獲取的目標(biāo)徑向長度信息中的任一時(shí)間點(diǎn)、及1η、112、φ1:參數(shù)矩 陣中的每一參數(shù)數(shù)組代入公式6,得到對(duì)應(yīng)于每一參數(shù)數(shù)組的a個(gè)運(yùn)算結(jié)果;計(jì)算上述a個(gè)運(yùn) 算結(jié)果與該時(shí)間點(diǎn)對(duì)應(yīng)的目標(biāo)徑向長度的差值;將得到的對(duì)應(yīng)于每一參數(shù)數(shù)組的a個(gè)差值 與第一閾值進(jìn)行比較，標(biāo)記差值小于第一閾值的參數(shù)數(shù)組，并在該參數(shù)數(shù)組對(duì)應(yīng)的第一累 加矩陣的行加1;
[0094] 并且，將第二雷達(dá)獲取的目標(biāo)徑向長度信息中的任一時(shí)間點(diǎn)、及121、122、φ2參數(shù)矩 陣中的每一參數(shù)數(shù)組代入公式7，得到對(duì)應(yīng)于每一參數(shù)數(shù)組的b個(gè)運(yùn)算結(jié)果;計(jì)算上述b個(gè)運(yùn) 算結(jié)果與該時(shí)間點(diǎn)對(duì)應(yīng)的目標(biāo)徑向長度的差值;將得到的對(duì)應(yīng)于每一參數(shù)數(shù)組的b個(gè)差值 與第二閾值進(jìn)行比較，標(biāo)記差值小于第二閾值的參數(shù)數(shù)組，并在該參數(shù)數(shù)組對(duì)應(yīng)的第二累 加矩陣的行加1。
[0095] S236.將第一雷達(dá)獲取的目標(biāo)徑向長度信息中的P個(gè)時(shí)間點(diǎn)按照步驟S235依次執(zhí) 行;選取第一累加矩陣中最大值所在行對(duì)應(yīng)的參數(shù)數(shù)組作為第一雷達(dá)目標(biāo)徑向長度曲線的 參數(shù)；
[0096]并且，將第二雷達(dá)獲取的目標(biāo)徑向長度信息中的Q個(gè)時(shí)間點(diǎn)按照步驟S235依次執(zhí) 行;選取第二累加矩陣中最大值所在行對(duì)應(yīng)的參數(shù)數(shù)組作為第二雷達(dá)目標(biāo)徑向長度曲線的 參數(shù)。
[0097] S3.利用第一雷達(dá)目標(biāo)徑向長度曲線、及第二雷達(dá)目標(biāo)徑向長度曲線計(jì)算目標(biāo)進(jìn) 動(dòng)角、及第一雷達(dá)與第二雷達(dá)的進(jìn)動(dòng)軸視線角差值，所述進(jìn)動(dòng)軸視線角為雷達(dá)視線與目標(biāo) 進(jìn)動(dòng)軸的夾角。
[0098] 在本發(fā)明優(yōu)選實(shí)施例中，步驟S3具體包括：
[0099] S31.利用第一雷達(dá)目標(biāo)徑向長度曲線、第二雷達(dá)目標(biāo)徑向長度曲線及公式1計(jì)算 目標(biāo)進(jìn)動(dòng)角。
[0100]
[0101 ]其中，Χι為第一雷達(dá)在第一時(shí)間段Τι內(nèi)，目標(biāo)徑向長度的最大值與最小值的和;Υι 為第一雷達(dá)在第一時(shí)間段!^內(nèi)，目標(biāo)徑向長度的最大值與最小值的差;X2為第二雷達(dá)在第二 時(shí)間段Τ2內(nèi)，目標(biāo)徑向長度的最大值與最小值的和;Υ2為第二雷達(dá)在第二時(shí)間段Τ2內(nèi)，目標(biāo) 徑向長度的最大值與最小值的差。
[0102] S32.利用第一雷達(dá)目標(biāo)徑向長度曲線、第二雷達(dá)目標(biāo)徑向長度曲線、目標(biāo)進(jìn)動(dòng)角 及公式2計(jì)算第一雷達(dá)與第二雷達(dá)的進(jìn)動(dòng)軸視線角差值。
[0103]
[0104] 其中，cu為第一雷達(dá)的進(jìn)動(dòng)軸視線角，且〇1從0°到90°按照第一步進(jìn)值設(shè)置為與 ai對(duì)應(yīng)的第二雷達(dá)的進(jìn)動(dòng)軸視線角；i > 1、且i ΕΝ。
[0105] 在本發(fā)明的優(yōu)選實(shí)施例中，若第一步進(jìn)值為5°，則〇1為0°、5°、10°、15°、20°、 25°……90° ;若第一步進(jìn)值為10°，則cti為10°、20°、30°……90°。
[0106] 公式1、公式2的推導(dǎo)過程如下：
[0107]假設(shè)在彈道中段某短時(shí)間段T內(nèi)，雷達(dá)進(jìn)動(dòng)軸視線角近似不變，通過步驟S2,在時(shí) 間段T1內(nèi)第一雷達(dá)能夠得到一組徑向長度Lem，在時(shí)間段T2ft第二雷達(dá)能夠得到一組徑向 長度Lem，即便TjPT 2都小于一個(gè)進(jìn)動(dòng)周期，都能夠估計(jì)徑向長度曲線，從而獲得一個(gè)周期T 內(nèi)徑向長度的最大值及最小值。對(duì)于第一雷達(dá)，根據(jù)公式12,可知徑向長度的最大值及最小 值出現(xiàn)在：⑶8(碼+私)=- 1和) =1 ，對(duì)應(yīng)的徑向長度分別為：
[0108] Len(ti) =Lcos(a-0)-Rsin(a-9)
[0109] Len(t2) =Lcos(a+0)-Rsin(a+0)
[011 0]則在?\時(shí)間段內(nèi)的徑向長度最大值及最小值為：
[0111] max(Lem) =Lcos(a-0)-Rsin(a-6)
[0112] min(Lem) =Lcos(a+0)-Rsin(a+6)
[0113] 令：
[0114] Xi=max(Leni)+min(Lem)
[0115] Yi=max(Leni)-min(Lem)
[0116] 可得：
[0117] Xi = 2LcosacosB-2RcosBsina
[0118] Yi = 2LsinasinB+2RcosasinB 公式 15
[0119] 由此：
[0136] 將h沭兩式聯(lián)立可得：
[0137]
[0138] 在步驟S3中，第一雷達(dá)及第二雷達(dá)在較短的觀測(cè)時(shí)間內(nèi)，即可準(zhǔn)確計(jì)算目標(biāo)進(jìn)動(dòng) 角及進(jìn)動(dòng)軸視線角差值，在此基礎(chǔ)上，本發(fā)明能夠?qū)崿F(xiàn)短時(shí)間對(duì)目標(biāo)結(jié)構(gòu)參數(shù)的提取;與現(xiàn) 有技術(shù)相比，大幅減小了提取進(jìn)動(dòng)特征所需的觀測(cè)時(shí)間。同時(shí)，通過公式1，本發(fā)明實(shí)現(xiàn)了以 車父尚精度對(duì)目標(biāo)進(jìn)動(dòng)角的提取。
[0139] S4.基于目標(biāo)進(jìn)動(dòng)角、進(jìn)動(dòng)軸視線角差值、第一雷達(dá)目標(biāo)徑向長度曲線及第二雷達(dá) 目標(biāo)徑向長度曲線計(jì)算目標(biāo)長度及目標(biāo)底面半徑。
[0140] 在本發(fā)明優(yōu)選實(shí)施例中，步驟S4具體包括：
[0141] S41.基于目標(biāo)進(jìn)動(dòng)角、進(jìn)動(dòng)軸視線角差值及公式3、公式4計(jì)算對(duì)應(yīng)于每一個(gè)第一 雷達(dá)進(jìn)動(dòng)軸視線角<^的目標(biāo)樣本長度U及目標(biāo)樣本底面半徑心；
[0145] 其中，Τη…以為!^內(nèi)的N個(gè)時(shí)間點(diǎn)，Τ21···Τ2μ為T 2內(nèi)的Μ個(gè)時(shí)間點(diǎn);Μ>1，Ν>1，ΜΕΝ， NEN;Len(Tii)為與Τιι，··Τιν對(duì)應(yīng)的目標(biāo)徑向長度，Len(T2i)Τ2Μ)為與Τ2?··· T2M對(duì)應(yīng)的目標(biāo)徑向長度；ω為目標(biāo)進(jìn)動(dòng)角速度；φ?、啊2分別為第一雷達(dá)、第二雷達(dá)視角下 的目標(biāo)進(jìn)動(dòng)初相角；Tir·· tin為與Tir·· Tin對(duì)應(yīng)的第一雷達(dá)視線角，Τ2?···Τ2Μ為與Τ2?···Τ2Μ對(duì)應(yīng) 的第二雷達(dá)視線角;且
[0146
[0147] 公式3、公式4由公式12得到。
[0148] 在本發(fā)明優(yōu)選實(shí)施例中，利用最小二乘法求得對(duì)應(yīng)于每一個(gè)第一雷達(dá)進(jìn)動(dòng)軸視線 角<^的目標(biāo)樣本長度U及目標(biāo)樣本底面半徑心。
[0149] S42.根據(jù)目標(biāo)樣本長度U及目標(biāo)樣本底面半徑&提取目標(biāo)長度及目標(biāo)底面半徑。
[0150]在本發(fā)明優(yōu)選實(shí)施例中，步驟S42具體包括：
[0151 ] S421.根據(jù)公式5計(jì)算每一組目標(biāo)樣本長度h及目標(biāo)樣本底面半徑1^的誤差值e (〇i)；
[0152]
[0153] S422.比較所有誤差值的大小，并提取最小的誤差值對(duì)應(yīng)的目標(biāo)樣本長度作為目 標(biāo)長度，提取最小的誤差值對(duì)應(yīng)的目標(biāo)樣本底面半徑作為目標(biāo)底面半徑。
[0154]在步驟S4中，不需要了解進(jìn)動(dòng)目標(biāo)在各種姿態(tài)下的RCS特性或形狀，即可以較高的 精度提取其結(jié)構(gòu)參數(shù)。
[0155] 在本發(fā)明優(yōu)選實(shí)施例中，本發(fā)明提供的進(jìn)動(dòng)目標(biāo)結(jié)構(gòu)參數(shù)提取方法，在步驟S1之 前，還包括：
[0156] S0.從U部雷達(dá)中任取兩部合為一組，形成V組雷達(dá);根據(jù)所述V組雷達(dá)中任一組中 的兩部雷達(dá)提取目標(biāo)長度k及目標(biāo)底面半徑心；
[0157] 在步驟S4之后，還包括：
[0158] S5.根據(jù)目標(biāo)長度b及目標(biāo)底面半徑&計(jì)算最佳目標(biāo)長度及最佳目標(biāo)底面半徑;其 中，U>2,且UEN; ; 1 < j < V，且jEN。
[0159]在本發(fā)明優(yōu)選實(shí)施例中，步驟S5具體為：
[0160]根據(jù)公式8計(jì)算k的平均值L#，將k的平均值1^#作為最佳目標(biāo)長度；
[0161] 根據(jù)公式9計(jì)算心的平均值1?#，將心的平均值妒作為最佳目標(biāo)底面半徑；
[0162]
[0163]
[0164] 在本發(fā)明優(yōu)選實(shí)施例中，步驟S5具體為：
[0165] 根據(jù)公式10計(jì)算k的加權(quán)平均值ΙΛ將k的加權(quán)平均值f作為最佳目標(biāo)長度；
[0166] 根據(jù)公式11計(jì)算Rj的加權(quán)平均值R'將Rj的加權(quán)平均值浐作為最佳目標(biāo)底面半徑；
[0167]
[0168]
[0169] 其中，Wj為V組雷達(dá)中第j組的權(quán)值，Wj根據(jù)第j組中兩部雷達(dá)中的每一雷達(dá)的性能 參數(shù)如測(cè)量精度、探測(cè)距離、穩(wěn)定性參數(shù)、可靠性參數(shù)等確定。
[0170] 在本發(fā)明優(yōu)選實(shí)施例中，基于V組雷達(dá)中每一組得到的目標(biāo)進(jìn)動(dòng)角，計(jì)算上述目標(biāo) 進(jìn)動(dòng)角的平均值或加權(quán)平均值，作為最佳目標(biāo)進(jìn)動(dòng)角。
[0171]由此，實(shí)現(xiàn)了通過多于兩部的雷達(dá)提取目標(biāo)結(jié)構(gòu)參數(shù)及目標(biāo)進(jìn)動(dòng)角的技術(shù)效果， 進(jìn)一步提高了目標(biāo)結(jié)構(gòu)參數(shù)及目標(biāo)進(jìn)動(dòng)角的提取精度。
[0172] 圖3是本發(fā)明的進(jìn)動(dòng)目標(biāo)結(jié)構(gòu)參數(shù)提取方法的第二流程圖；