一種外輻射源雷達(dá)同頻干擾的消除方法
【專利摘要】本發(fā)明屬于雷達(dá)技術(shù)領(lǐng)域，公開了一種外輻射源雷達(dá)同頻干擾的消除方法。該方法包括：獲取參考天線接收的參考信號(hào)及陣元天線接收的M個(gè)回波信號(hào)；利用參考信號(hào)消除M個(gè)回波信號(hào)中的直達(dá)波信號(hào)以及多徑信號(hào)得到M個(gè)第一回波信號(hào)，構(gòu)造得到回波信號(hào)矩陣；對(duì)回波信號(hào)矩陣進(jìn)行獨(dú)立分量分析得到3個(gè)獨(dú)立分量信號(hào)；利用3個(gè)獨(dú)立分量信號(hào)對(duì)回波信號(hào)矩陣進(jìn)行3次時(shí)域干擾相消，得到3個(gè)第二回波信號(hào)；對(duì)每個(gè)第二回波信號(hào)與參考信號(hào)作距離?多普勒二維相關(guān)得到3個(gè)相關(guān)矩陣，利用相關(guān)矩陣進(jìn)行目標(biāo)檢測(cè)，其中能夠獲得目標(biāo)信息的相關(guān)矩陣對(duì)應(yīng)的第二回波信號(hào)即為消除同頻干擾后的回波信號(hào)。本發(fā)明能夠有效消除主瓣內(nèi)的同頻干擾，從而獲得較好的目標(biāo)檢測(cè)結(jié)果。
【專利說明】
一種外輻射源雷達(dá)同頻干擾的消除方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及雷達(dá)技術(shù)領(lǐng)域，尤其涉及一種外輻射源雷達(dá)同頻干擾的消除方法。
【背景技術(shù)】
[0002] 外輻射源雷達(dá)借助第三方輻射源信號(hào)(如調(diào)頻廣播信號(hào)、模擬電視信號(hào)、通信信號(hào) 等)作為照射源照射目標(biāo)進(jìn)行探測(cè)，其探測(cè)過程屬于被動(dòng)探測(cè)。因此對(duì)外輻射雷達(dá)而言，相 比干擾和噪聲，其目標(biāo)回波信號(hào)十分微弱，所以對(duì)干擾與噪聲的抑制程度會(huì)直接影響外輻 射雷達(dá)的探測(cè)性能。
[0003] 外輻射源雷達(dá)體系下的干擾包括相干干擾和非相干干擾，其中，相干干擾指的是 直達(dá)波干擾與多徑干擾，非相干干擾指的是來自輻射源的同頻源的同頻干擾。由于外輻射 源雷達(dá)所接收的回波信號(hào)中干擾信號(hào)的功率遠(yuǎn)強(qiáng)于目標(biāo)回波信號(hào)的功率，因此直接根據(jù)所 接收的回波信號(hào)很難檢測(cè)到目標(biāo)，所以需要先濾除其中的干擾信號(hào)。目前，對(duì)于相干干擾， 可采用時(shí)域干擾相消算法來消除。對(duì)于同頻干擾，則可采用空域?yàn)V波的方法來濾除，具體來 說，可通過自適應(yīng)波束形成技術(shù)在干擾來波方向形成方向圖零點(diǎn)以濾除同頻干擾。然而，存 在同頻干擾與目標(biāo)回波的來波方向位于同一主瓣內(nèi)的情況，此種情況下，現(xiàn)有的空域?yàn)V波 算法即不再適用。對(duì)于主瓣內(nèi)的同頻干擾，目前還沒有行之有效的消除方法。

【發(fā)明內(nèi)容】

[0004] 為此，本發(fā)明的實(shí)施例提供一種外輻射源雷達(dá)同頻干擾的消除方法，能夠在目標(biāo) 回波通道中混有時(shí)域和空域都不可分的同頻干擾時(shí)，有效消除主瓣內(nèi)的同頻干擾，從而獲 得較好的目標(biāo)檢測(cè)結(jié)果。
[0005] 為達(dá)到上述目的，本發(fā)明的實(shí)施例采用如下技術(shù)方案：
[0006] 提供一種外輻射源雷達(dá)同頻干擾的消除方法，所述外輻射源雷達(dá)設(shè)置有參考天線 和接收天線，所述參考天線指向所述外輻射源雷達(dá)的輻射源，所述接收天線指向所述外輻 射源雷達(dá)的觀測(cè)區(qū)；所述接收天線為包括M個(gè)陣元天線的均勻線陣，所述M個(gè)陣元天線的相 位中心間距為半波長，M為整數(shù)，M多3;所述方法包括如下步驟：
[0007] S1:獲取所述參考天線接收的參考信號(hào)，以及，獲取所述M個(gè)陣元天線中每個(gè)陣元 天線接收的回波信號(hào)得到M個(gè)回波信號(hào)；
[0008] S2:利用所述參考信號(hào)，消除M個(gè)回波信號(hào)中的直達(dá)波信號(hào)以及多徑信號(hào)，得到M個(gè) 第一回波信號(hào)，進(jìn)而利用所述M個(gè)第一回波信號(hào)構(gòu)造得到回波信號(hào)矩陣；
[0009] S3:對(duì)所述回波信號(hào)矩陣進(jìn)行獨(dú)立分量分析，得到3個(gè)獨(dú)立分量信號(hào)，其中1個(gè)獨(dú)立 分量信號(hào)為同頻干擾信號(hào)；
[0010] S4:利用所述3個(gè)獨(dú)立分量信號(hào)對(duì)所述回波信號(hào)矩陣進(jìn)行時(shí)域干擾相消，得到3個(gè) 第二回波信號(hào)；其中，每次時(shí)域干擾相消分別消除所述3個(gè)獨(dú)立分量信號(hào)中的1個(gè)獨(dú)立分量 信號(hào)，得到對(duì)應(yīng)的1個(gè)第二回波信號(hào)，每次時(shí)域干擾相消所消除的獨(dú)立分量信號(hào)不同；
[0011] S5:對(duì)每個(gè)第二回波信號(hào)與所述參考信號(hào)作距離-多普勒二維相關(guān)，得到3個(gè)相關(guān) 矩陣;利用每個(gè)相關(guān)矩陣進(jìn)行目標(biāo)檢測(cè)，其中能夠獲得目標(biāo)信息的相關(guān)矩陣對(duì)應(yīng)的第二回 波信號(hào)即為所求的消除同頻干擾后的回波信號(hào)。
[0012] 基于上述方案，通過對(duì)回波信號(hào)矩陣進(jìn)行獨(dú)立分量分析，可得到3個(gè)獨(dú)立分量信 號(hào)，并且其中一個(gè)獨(dú)立分量信號(hào)為同頻干擾信號(hào)。由于同頻干擾信號(hào)的強(qiáng)度遠(yuǎn)大于目標(biāo)回 波信號(hào)及噪聲，因此消除同頻干擾信號(hào)后的回波信號(hào)的主要成分即為目標(biāo)回波信號(hào)。這樣， 利用分別消除了 3個(gè)獨(dú)立分量信號(hào)得到的3個(gè)第二回波信號(hào)與參考信號(hào)作距離-多普勒二維 相關(guān)得到3個(gè)相關(guān)矩陣，進(jìn)而利用每個(gè)相關(guān)矩陣進(jìn)行目標(biāo)檢測(cè)，其中能夠獲得目標(biāo)信息的相 關(guān)矩陣對(duì)應(yīng)的第二回波信號(hào)即為所求的消除同頻干擾后的回波信號(hào)。綜上，本發(fā)明實(shí)施例 提供的外輻射源雷達(dá)同頻干擾的消除方法能夠在目標(biāo)回波通道中混有時(shí)域和空域都不可 分的同頻干擾時(shí)，有效消除主瓣內(nèi)的同頻干擾，從而獲得較好的目標(biāo)檢測(cè)結(jié)果。
【附圖說明】
[0013] 為了更清楚地說明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見地，下面描述中的附圖僅僅是本 發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以 根據(jù)這些附圖獲得其他的附圖。
[0014]圖1為本發(fā)明實(shí)施例中外輻射源雷達(dá)基地配置示意圖；
[0015] 圖2為本發(fā)明實(shí)施例提供的一種外輻射源雷達(dá)同頻干擾的消除方法的流程示意 圖；
[0016] 圖3為直接利用雷達(dá)接收天線接收的回波信號(hào)進(jìn)行目標(biāo)檢測(cè)的檢測(cè)結(jié)果圖；其中， 圖3(a)為多普勒維的檢測(cè)結(jié)果圖，圖中橫坐標(biāo)為多普勒頻移單元(單位:kHz)，縱坐標(biāo)為幅 度(單位:dB);圖3(b)為時(shí)延維的檢測(cè)結(jié)果圖，圖中橫坐標(biāo)為回波時(shí)延單元(單位:km)，縱坐 標(biāo)為幅度(單位:dB);
[0017] 圖4為利用消除獨(dú)立分量信號(hào)71后的回波信號(hào)進(jìn)行目標(biāo)檢測(cè)的檢測(cè)結(jié)果圖；其中， 圖4(a)為多普勒維的檢測(cè)結(jié)果圖，圖中橫坐標(biāo)為多普勒頻移單元(單位:kHz)，縱坐標(biāo)為幅 度(單位:dB);圖4(b)為時(shí)延維的檢測(cè)結(jié)果圖，圖中橫坐標(biāo)為回波時(shí)延單元(單位:km)，縱坐 標(biāo)為幅度(單位:dB);
[0018] 圖5為利用消除獨(dú)立分量信號(hào)72后的回波信號(hào)進(jìn)行目標(biāo)檢測(cè)的檢測(cè)結(jié)果圖；其中， 圖5(a)為多普勒維的檢測(cè)結(jié)果圖，圖中橫坐標(biāo)為多普勒頻移單元(單位:kHz)，縱坐標(biāo)為幅 度(單位:dB);圖5(b)為時(shí)延維的檢測(cè)結(jié)果圖，圖中橫坐標(biāo)為回波時(shí)延單元(單位:km)，縱坐 標(biāo)為幅度(單位:dB);
[0019] 圖6為利用消除獨(dú)立分量信號(hào)y3后的回波信號(hào)進(jìn)行目標(biāo)檢測(cè)的檢測(cè)結(jié)果圖；其中， 圖6(a)為多普勒維的檢測(cè)結(jié)果圖，圖中橫坐標(biāo)為多普勒頻移單元(單位:kHz)，縱坐標(biāo)為幅 度(單位:dB);圖6(b)為時(shí)延維的檢測(cè)結(jié)果圖，圖中橫坐標(biāo)為回波時(shí)延單元(單位:km)，縱坐 標(biāo)為幅度(單位:dB)。
【具體實(shí)施方式】
[0020] 下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完 整地描述，顯然，所描述的實(shí)施例僅僅是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例?；?本發(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他 實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。
[0021] 另外，為了便于清楚描述本發(fā)明實(shí)施例的技術(shù)方案，在本發(fā)明的實(shí)施例中，采用了 "第一"、"第二"等字樣對(duì)功能和作用基本相同的相同項(xiàng)或相似項(xiàng)進(jìn)行區(qū)分，本領(lǐng)域技術(shù)人 員可以理解"第一"、"第二"、等字樣并不對(duì)數(shù)量和執(zhí)行次序進(jìn)行限定。
[0022] 為便于理解，以下將結(jié)合圖1所示的外輻射源雷達(dá)的應(yīng)用場(chǎng)景，對(duì)外輻射源雷達(dá)的 工作原理進(jìn)行簡(jiǎn)要介紹如下：
[0023]參見圖1，在外輻射源雷達(dá)接收站的遠(yuǎn)場(chǎng)設(shè)置有第三方輻射源，第三方輻射源作為 發(fā)射站輻射電磁波信號(hào)。當(dāng)輻射源輻射的電磁波信號(hào)照射目標(biāo)上時(shí)，目標(biāo)會(huì)將電磁波信號(hào) 反射出去。外輻射源雷達(dá)通過設(shè)置的接收天線即可接收目標(biāo)反射的電磁波信號(hào)(通常稱為 目標(biāo)回波信號(hào)）。同時(shí)，輻射源輻射的部分電磁波信號(hào)會(huì)直接照射在外輻射源雷達(dá)的接收站 上，形成直達(dá)波信號(hào)。外輻射源雷達(dá)通過設(shè)置的參考天線即可接收直達(dá)波信號(hào)，進(jìn)而以直達(dá) 波信號(hào)作為參考信號(hào)，利用所接收的目標(biāo)回波信號(hào)估計(jì)目標(biāo)的位置及速度等信息，實(shí)現(xiàn)對(duì) 目標(biāo)的探測(cè)、跟蹤及定位。
[0024] 然而，如圖1所示，作為發(fā)射站的第三方輻射源周圍往往存在一些建筑物，當(dāng)輻射 源輻射的電磁波信號(hào)照射在周邊的建筑物上，即會(huì)產(chǎn)生多徑信號(hào)，部分多徑信號(hào)即會(huì)隨目 標(biāo)回波信號(hào)一同被外輻射源雷達(dá)接收。同時(shí)，空間中可能存在輻射源的同頻源，當(dāng)同頻源發(fā) 射的電磁波信號(hào)照射在外輻射源雷達(dá)的接收站上，同頻信號(hào)也會(huì)一同被接收天線接收。此 外，在接收目標(biāo)回波信號(hào)的同時(shí)，還可能會(huì)接收到輻射源直接輻射的直達(dá)波信號(hào)。綜上，外 輻射源雷達(dá)接收天線所接收的回波信號(hào)中除包含目標(biāo)回波信號(hào)之外，還可能包含多徑信 號(hào)、同頻信號(hào)以及直達(dá)波信號(hào)。其中，直達(dá)波信號(hào)與多徑信號(hào)為相干干擾信號(hào)，同頻信號(hào)為 非相干干擾信號(hào)。目前，對(duì)于非相干干擾信號(hào)中的同頻信號(hào)，可采用自適應(yīng)波束形成技術(shù)在 干擾來波方向形成方向圖零點(diǎn)予濾除。然而，存在同頻干擾與目標(biāo)回波的來波方向位于同 一主瓣內(nèi)的情況，此種情況下，現(xiàn)有方法即不再適用。
[0025] 為此，本發(fā)明的實(shí)施例提供一種外輻射源雷達(dá)同頻干擾的消除方法，以在目標(biāo)回 波通道中混有時(shí)域和空域都不可分的同頻干擾時(shí)，消除主瓣內(nèi)的同頻干擾，以獲得較好的 目標(biāo)檢測(cè)結(jié)果。其中，所述外輻射源雷達(dá)設(shè)置有參考天線和接收天線，參考天線指向外輻射 源雷達(dá)的輻射源，接收天線指向外輻射源雷達(dá)的觀測(cè)區(qū)。接收天線為包括M個(gè)陣元天線的均 勻線陣，M個(gè)陣元天線的相位中心間距為半波長，M為整數(shù)，M多3。
[0026] 如圖2所示，本發(fā)明實(shí)施例提供的外輻射源雷達(dá)同頻干擾的消除方法具體包括以 下步驟：
[0027] S1:獲取參考天線接收的參考信號(hào)，以及，獲取M個(gè)陣元天線中每個(gè)陣元天線接收 的回波信號(hào)得到M個(gè)回波信號(hào)。
[0028]其中，本領(lǐng)域普通技術(shù)人員可以理解，參考天線接收的參考信號(hào)中包含直達(dá)波信 號(hào)及多徑信號(hào)，接收天線的每個(gè)陣元天線接收的回波信號(hào)中包含目標(biāo)回波信號(hào)以及直達(dá)波 信號(hào)、同頻信號(hào)及多徑信號(hào)構(gòu)成的干擾信號(hào)。
[0029] S2:利用參考信號(hào)，消除M個(gè)回波信號(hào)中的直達(dá)波信號(hào)以及多徑信號(hào)，得到M個(gè)第一 回波信號(hào)，進(jìn)而利用M個(gè)第一回波信號(hào)構(gòu)造得到回波信號(hào)矩陣。
[0030] 步驟S2的理論依據(jù)在于：回波信號(hào)中的目標(biāo)回波信號(hào)相對(duì)于直達(dá)波干擾信號(hào)及多 徑干擾信號(hào)具有多普勒頻移，因此，可認(rèn)為目標(biāo)回波信號(hào)與直達(dá)波干擾信號(hào)及多徑干擾信 號(hào)是正交的?；诖耍纯衫脜⒖夹盘?hào)（包括直達(dá)波干擾信號(hào)以及多徑干擾信號(hào)）張成的 信號(hào)子空間構(gòu)造正交投影矩陣，并將回波信號(hào)投影在所構(gòu)造的正交投影矩陣上，以獲得回 波信號(hào)中的直達(dá)波干擾和多徑干擾成分，進(jìn)而減去回波信號(hào)中的直達(dá)波干擾和多徑干擾成 分，即可消除回波信號(hào)中直達(dá)波干擾和多徑干擾成分。
[0031 ] 具體來說，可采用擴(kuò)展相消批處理算法（英文：Batch Version of Extensive Cancellation Algorithm，簡(jiǎn)稱:ECA-B)消除M個(gè)回波信號(hào)中的直達(dá)波信號(hào)以及多徑信號(hào)以 獲得M個(gè)第一回波信號(hào)，進(jìn)而利用M個(gè)第一回波信號(hào)構(gòu)造得到回波信號(hào)矩陣。
[0032] 即，步驟S2具體可以包括如下子步驟：
[0033] S2a:分別對(duì)參考信號(hào)以及M個(gè)回波信號(hào)進(jìn)行分段，得到參考信號(hào)的B段分段信號(hào)以 及每個(gè)回波信號(hào)的B段分段信號(hào)。
[0034]其中，B為正整數(shù)，N為參考信號(hào)及回波信號(hào)的信號(hào)長度。
[0035]優(yōu)選的，本發(fā)明實(shí)施例中，參考信號(hào)及M個(gè)回波信號(hào)的每段分段信號(hào)的長度不小于 1000?？梢岳斫猓瑢⑿盘?hào)劃分為B段后，每段分段信號(hào)的長度為
若要使每段分段信號(hào)的長 度不小于1000,則對(duì)應(yīng)的_
，因此需要將參考信號(hào)及回波信號(hào)均勻劃分
。其中，表示向下取整。
[0036] S2b:利用參考信號(hào)的第i段分段信號(hào)
構(gòu)造與其對(duì)應(yīng)的子空間矩陣
[0037] 其中，i取1到B之間的所有整數(shù)值，
表示參考信號(hào)的第i段分段信號(hào)，
依次表示將參考信號(hào)的第i段分段信號(hào)
延遲1個(gè)時(shí)延單元、2個(gè)時(shí)延單元、???（？-：〇個(gè)時(shí)延單元后的信號(hào)，p = 1+1，1表示回波 信號(hào)中多徑干擾的時(shí)延單元個(gè)數(shù)，1多1。需要說明的是，在實(shí)際應(yīng)用中，可根據(jù)經(jīng)驗(yàn)對(duì)參數(shù)P 進(jìn)行賦值。具體而言，當(dāng)外輻射源雷達(dá)的接收站位于較為空曠的地方時(shí)，對(duì)P賦以相對(duì)較小 的值;而當(dāng)外輻射源雷達(dá)的接收站位于建筑物較多的地方時(shí)，對(duì)P賦以相對(duì)較大的值。
[0038] 示例性的，假設(shè)參考信號(hào)的第b段分段信號(hào)
可以表示為：
，則利用參考信號(hào)的第b段分段信
構(gòu)造的與其對(duì)應(yīng)的子空間矩陣可表示為：
[0040]
表示第b段分段信號(hào)
）的起始數(shù)據(jù)
丨表示第b段
分段信號(hào)
的結(jié)尾數(shù)據(jù)。
[0041 ] S2c:利用子空間矩陣^構(gòu)造第一正交投影矩陣叫=& ，進(jìn)而利 用第一正交投影矩陣對(duì)第j個(gè)回波信號(hào)的第i段分段信號(hào)
進(jìn)行時(shí)域干擾相消，得 到消除直達(dá)波信號(hào)及多徑信號(hào)后的第i段分段信號(hào)
[0042]其中，j取1到M之間的所有整數(shù)值，H共輒轉(zhuǎn)置，-1表示求逆。第j個(gè)回波信號(hào)的第i 段分段信號(hào)
_可表不為：
[0044] 式中，
表示第j個(gè)回波信號(hào)的第i段分段信號(hào)+
的起始數(shù)據(jù)，
表示第j個(gè)回波信號(hào)的第i段分段信號(hào)
的結(jié)尾數(shù)據(jù)。
[0045] S2d:對(duì)M個(gè)回波信號(hào)的消除直達(dá)波信號(hào)以及多徑信號(hào)后的B段分段信號(hào)進(jìn)行拼接， 得到M個(gè)第一回波信號(hào)。
[0046] S2e:利用M個(gè)第一回波信號(hào)構(gòu)造得到回波信號(hào)矩陣KS5,Hif。
[0047] 其中，*^*1 _，足《r2，…，式!^分別表不M個(gè)第一回波信號(hào)中的第1個(gè)第一回波信號(hào)、第2個(gè) 第一回波信號(hào)、…第M個(gè)第一回波信號(hào)。
[0048] S3:對(duì)回波信號(hào)矩陣進(jìn)行獨(dú)立分量分析，得到3個(gè)獨(dú)立分量信號(hào)。
[0049] 其中，一個(gè)獨(dú)立分量信號(hào)為同頻干擾信號(hào)。
[0050] 步驟3的理論依據(jù)在于:經(jīng)過步驟S2,可消除M個(gè)回波信號(hào)中的直達(dá)波干擾和多徑 干擾，這樣一來，回波信號(hào)中即剩余同頻干擾信號(hào)、目標(biāo)回波信號(hào)以及噪聲，即M個(gè)回波信號(hào) 即可表示為 =&+民+?，式中，j取1到M之間的所有整數(shù)值，Sd表示回波信號(hào)中的同頻 干擾，S t表示回波信號(hào)中的目標(biāo)回波信號(hào)，n表示回波信號(hào)中的噪聲。同時(shí)，同頻干擾Sd、目標(biāo) 回波信號(hào)St以及噪聲n彼此相互獨(dú)立，因此，可將回波信號(hào)分離為3個(gè)相互獨(dú)立的分量。并 且，由于同頻干擾信號(hào)的強(qiáng)度遠(yuǎn)強(qiáng)于噪聲和目標(biāo)回波信號(hào)，因此經(jīng)過獨(dú)立分量分析之后，同 頻干擾信號(hào)的分離效果最好，進(jìn)而可以認(rèn)為分離回波信號(hào)得到的3個(gè)分量信號(hào)中的某一獨(dú) 立分量即是回波信號(hào)中的同頻干擾信號(hào)。
[0051 ] 具體而言，可采用快速獨(dú)立分量分析（英文：Fast Independent Component Analysis，簡(jiǎn)稱:FastICA)算法對(duì)回波信號(hào)矩陣進(jìn)行獨(dú)立分量分析，即，步驟S3具體可以包 括如下子步驟：
[0052] S3a:計(jì)算回波信號(hào)矩陣的協(xié)方差矩陣，并對(duì)協(xié)方差矩陣進(jìn)行特征分解得到特征值 及特征向量，利用特征值以及特征向量構(gòu)造得到白化向量V =『1/2ET。
[0053] 其中，E = (ei，…，en)，D = diag(di，…，dn)，ei，…，en為協(xié)方差矩陣的特征向量， cU，…，dn為協(xié)方差矩陣的特征值，（?廣1/2表示對(duì)（?）中的元素開平方后對(duì)矩陣求逆，T表示 轉(zhuǎn)置。
[0054] 以下，對(duì)白化向量V = D4/2ET的得出過程分析如下：
[0055] 容易理解，對(duì)回波信號(hào)矩陣進(jìn)行白化預(yù)處理實(shí)質(zhì)上即是對(duì)回波信號(hào)進(jìn)行去相關(guān)和 縮放，得到白化后的回波信號(hào)矩陣，使得白化后的回波信號(hào)矩陣的協(xié)方差矩陣為單位矩陣。 換言之，白化預(yù)處理過程即是尋找白化向量V，使得白化后的回波信號(hào)矩陣滿 足:S "丨=/，1為單位矩陣。 *> SUf 5ilU 》
[0056] 為此，可以選擇線性變化V =『1/2ET，這樣一來，白化后的回波信號(hào)矩陣&的協(xié)方 差矩陣為：
[0058] 即，V = D~2ET滿足使的要J :，因此可選擇線性變化V =『V2ET作為 白化向量。
[0059] S3b:利用白化向量V對(duì)回波信號(hào)矩陣進(jìn)行白化預(yù)處理，得到白化信號(hào)<".= V&。
[0060] 其中，Ssur表示回波信號(hào)矩陣。
[0061 ] S3c :將待估計(jì)的獨(dú)立成分個(gè)數(shù)設(shè)置為3，
作為 Fast ICA中不動(dòng)點(diǎn)迭代的基本公式，迭代獲得3個(gè)權(quán)矢量wpl、Wp4PWp3。
[0062] 以下，對(duì)迭代公式
的得出過程分析如下：
[0063] 獨(dú)立分量分析即是尋找線性分咼矩陣1=[￥1，'\￥2，￥3]1，使得)/=|^.>， :.>，」/=認(rèn)/義_,,|,. 的3個(gè)行向量71、72和73相互獨(dú)立?？赏ㄟ^負(fù)熵估計(jì)函數(shù)尋找線性分量矩陣W，具體而言，使用 一般形式的非二次函數(shù)的期望來估計(jì)負(fù)熵可得到如下所示的負(fù)熵估計(jì)函數(shù)：
[0064] J(y)^ki(E{G1(y)})2+k2(E{G 2(y)}-E{G2(v)})2
[0065] 其中，ki和k2為正整數(shù)，G1和G2為非二次函數(shù)，G1為奇函數(shù)，G 2為偶函數(shù)，v為零均值 單位方差的高斯變量。
[0066] 對(duì)上式所示的負(fù)熵估計(jì)函數(shù)進(jìn)行進(jìn)一步簡(jiǎn)化，可得：
[0067] J(y)^ki(E{G(y)})2+k2(E{G(y)}-E{G(v)})2
[0068] ^(ki+k2)(E{G(y)})2-2k2E{G(y)}E{G(v)}+k 2(E{G(v)})2
[0069] ^(ki+k2)[E{G(y)}-E{G(v)}]2+2kiE{G(y)}E{G(v)}-ki(E{G(v)}) 2
[0070] 進(jìn)一步的，省略約等于號(hào)右側(cè)的后兩項(xiàng)，使用一個(gè)非二次函數(shù)G來表示負(fù)熵估計(jì)函 數(shù)J(y)，有：
[0071] J(y)^(ki+k2)[E{G(y)}-E{G(v)}]2
[0072] 根據(jù)上式可以得到：
[0073] J(y)cx[E{G(y)}-E{G(v)}]2
[0074]式中，G可以是任意的實(shí)際非二次函數(shù)，常用的G有G(y)=y4,或G(y)=-eXp(-y 2/ 2)，為簡(jiǎn)便起見，本發(fā)明實(shí)施例中采用G(y)=y4,其導(dǎo)數(shù)g = 4y3。
[0075 ]由于 r = [ , y2, y3 f =' [ h\，w 2，w'3 f 則=v/ ，將 y .代入最終簡(jiǎn)化后 的負(fù)熵估計(jì)函數(shù)，有：
[0077] 對(duì)h忒的傷熵估i+通教T(v)龍?zhí)荻鹊玫剑?br>[0079] 式中，wT表示線性分離矩陣W的行向量，
，函數(shù)g為負(fù)熵 估計(jì)函數(shù)J(y)中非二次函數(shù)G的導(dǎo)數(shù)，g = 4y3。
[0080] 令y =1，使用不定點(diǎn)迭代算法對(duì)上式求解，可以得到：
[0082]對(duì)w進(jìn)行歸一化處理，同時(shí)，由于非多項(xiàng)式矩的收斂性較差，不利于迭代，因此將上 式調(diào)整為：
[0084]接下來，需求解式中0的取值。首先，
'求F關(guān)于w的梯
，.有：

[0096]具體的，結(jié)合上式，使用FastICA算法獲取3個(gè)權(quán)矢量wpi、wP2和w P3的具體步驟如下： [0097] 步驟1:初始化獨(dú)立成分個(gè)數(shù)m = 3，并令p = 1。
[0098]步驟2:選擇具有單位范數(shù)的初始化向量wP (可以隨機(jī)選?。?。
[0101] 步驟5:對(duì)wp進(jìn)行標(biāo)準(zhǔn)化:wp-wp/I I wp I I。
[0102] 步驟6:若wp尚未收斂，返回步驟3。
[0103]步驟7:若wp收斂，保存wp，并令p加1，若p<m，則返回步驟2。
[0104] 至此，即得到3個(gè)權(quán)矢量Wpi、wP2和wP3。
[0105] S3d:將3個(gè)權(quán)矢量wpi、wP2和wP3分別與白化信號(hào)1§"相乘，得到3個(gè)獨(dú)立分量信號(hào)yi、 sur y2 和 y3。
[0107]需要說明的是，由于獨(dú)立分量分析存在不確定性，即分量幅度和順序不確定，因此 無法確定3個(gè)獨(dú)立分量yi、y2和y3中哪個(gè)獨(dú)立分量為同頻干擾信號(hào)。本發(fā)明實(shí)施例中，可通過 分別消除3個(gè)獨(dú)立分量信號(hào)yi、y2和y3,由于3個(gè)獨(dú)立分量信號(hào)yi、y2和y3中有一個(gè)獨(dú)立分量信 號(hào)為同頻干擾信號(hào)，并且同頻干擾信號(hào)的強(qiáng)度遠(yuǎn)大于目標(biāo)回波信號(hào)及噪聲，因此消除同頻 干擾信號(hào)后回波信號(hào)的主要成分即為目標(biāo)回波信號(hào)。這樣，利用消除同頻干擾信號(hào)的回波 信號(hào)即可檢測(cè)得到目標(biāo)的目標(biāo)信息。因此，可利用消除3個(gè)獨(dú)立分量信號(hào) yi、y.y3后的回波 信號(hào)進(jìn)行目標(biāo)檢測(cè)，進(jìn)而根據(jù)檢測(cè)結(jié)果確定消除同頻干擾的回波信號(hào)，具體過程請(qǐng)參考下 面的步驟S4和S5。
[0108] S4:利用3個(gè)獨(dú)立分量信號(hào)對(duì)回波信號(hào)矩陣進(jìn)行時(shí)域干擾相消，得到3個(gè)第二回波 信號(hào)。
[0109] 其中，每次時(shí)域干擾相消分別消除所述3個(gè)獨(dú)立分量信號(hào)中的1個(gè)獨(dú)立分量信號(hào)， 得到對(duì)應(yīng)的1個(gè)第二回波信號(hào)，每次時(shí)域干擾相消所消除的獨(dú)立分量信號(hào)不同。
[0110]優(yōu)選的，步驟S4具體可以包括如下子步驟：
[0111] S4a:利用3個(gè)獨(dú)立分量信號(hào)中的第k個(gè)獨(dú)立分量信號(hào)yk，構(gòu)造得到第二正交投影矩 陣 wk=(ykHyk)-VkH〇
[0112]其中，k分別取 1、2以及3。即，wi = (yiHyi)-ViH，W2 = (y2Hy2)-V2H，W3 = (y3Hy3)-V3H， yi、y2和y3為3個(gè)獨(dú)立分量信號(hào)。
[0113] S4b:利用第二正交投影矩陣概對(duì)回波信號(hào)矩陣進(jìn)行時(shí)域干擾相消，以消除第k個(gè) 獨(dú)立分量信號(hào)yk，獲得第二回波信號(hào)-叫^",.。
[0114] SP，之=& _ h.A,"，心.=Ss",.- >從,".，義= 陣。
[0115] 由前述分析可知，3個(gè)獨(dú)立分量信號(hào)yi、y2和y3中有一個(gè)獨(dú)立分量信號(hào)為同頻干擾 信號(hào)，因此分別消除獨(dú)立分量信號(hào)yi、ydPy 3得到3個(gè)第二回波信號(hào)義、匕以及兒,中有一 個(gè)信號(hào)為消除同頻干擾信號(hào)的回波信號(hào)。但由于獨(dú)立分量分析存在不確定性，即分量幅度 和順序不確定，因此無法確定3個(gè)獨(dú)立分量 yi、y4Py3中哪個(gè)獨(dú)立分量為同頻干擾信號(hào)，進(jìn)而 也就無法確定第二回波信號(hào)<",、Si和^"_中哪個(gè)信號(hào)為消除同頻干擾信號(hào)的回波信號(hào)。
[0116] S5:對(duì)每個(gè)第二回波信號(hào)與參考信號(hào)作距離-多普勒二維相關(guān)，得到3個(gè)相關(guān)矩陣， 利用每個(gè)相關(guān)矩陣進(jìn)行目標(biāo)檢測(cè)，其中能夠獲得目標(biāo)信息的相關(guān)矩陣對(duì)應(yīng)的第二回波信號(hào) 即為所求的消除同頻干擾后的回波信號(hào)。
[0117] 優(yōu)選的，步驟S5中，對(duì)每個(gè)第二回波信號(hào)與參考信號(hào)作距離-多普勒二維相關(guān)，得 到 3個(gè)相關(guān)矩陣，具體可以包括：
[0118] 利用預(yù)存的計(jì)算式：
，對(duì)每個(gè)第二回波信號(hào)與 參考信號(hào)作距離-多普勒二維相關(guān)，得到3個(gè)相關(guān)矩陣|1[l，p]、|2[l，p^&| 3[l，p]。
[0119] 其中，k分別取1、2以及表示第k個(gè)第二回波信號(hào)，Sref[n]表示參考信號(hào)， Sref[n-1]表示將參考信號(hào)Sref[n]延遲1個(gè)時(shí)延單元后的信號(hào)，n表示采樣次數(shù)，1表示時(shí)延單 元，*表示共輒運(yùn)算，P表示多普勒頻移單元，N表示信號(hào)長度。
[0120] 按照步驟S5,利用分別消除了 3個(gè)獨(dú)立分量信號(hào)得到的3個(gè)第二回波信號(hào)與參考信 號(hào)作距離-多普勒二維相關(guān)，得到3個(gè)相關(guān)矩陣，進(jìn)而利用每個(gè)相關(guān)矩陣進(jìn)行目標(biāo)檢測(cè)，由于 3個(gè)第二回波信號(hào)中有一個(gè)信號(hào)為消除同頻干擾后的回波信號(hào)，因此利用該信號(hào)與參考信 號(hào)作距離-多普勒二維相關(guān)得到的相關(guān)矩陣可檢測(cè)得到目標(biāo)信息。也就是說，能夠獲得目標(biāo) 信息的相關(guān)矩陣對(duì)應(yīng)的第二回波信號(hào)即為所求的消除同頻干擾后的回波信號(hào)。
[0121] 基于本發(fā)明實(shí)施例上述方案，通過對(duì)回波信號(hào)矩陣進(jìn)行獨(dú)立分量分析，可獲得3個(gè) 獨(dú)立分量信號(hào)，并且其中一個(gè)獨(dú)立分量信號(hào)為同頻干擾信號(hào)。由于同頻干擾信號(hào)的強(qiáng)度遠(yuǎn) 大于目標(biāo)回波信號(hào)及噪聲，因此消除同頻干擾信號(hào)后的回波信號(hào)的主要成分即為目標(biāo)回波 信號(hào)。這樣，利用分別消除了 3個(gè)獨(dú)立分量信號(hào)得到的3個(gè)第二回波信號(hào)與參考信號(hào)作距離- 多普勒二維相關(guān)得到3個(gè)相關(guān)矩陣，進(jìn)而利用每個(gè)相關(guān)矩陣進(jìn)行目標(biāo)檢測(cè)，其中能夠獲得目 標(biāo)信息的相關(guān)矩陣對(duì)應(yīng)的第二回波信號(hào)即為所求的消除同頻干擾后的回波信號(hào)。綜上，本 發(fā)明實(shí)施例提供的外輻射源雷達(dá)同頻干擾的消除方法能夠在目標(biāo)回波通道中混有時(shí)域和 空域都不可分的同頻干擾時(shí)，有效消除主瓣內(nèi)的同頻干擾，從而獲得較好的目標(biāo)檢測(cè)結(jié)果。
[0122] 以下，通過仿真實(shí)驗(yàn)進(jìn)一步說明本發(fā)明的上述效果：
[0123] 1)仿真實(shí)驗(yàn)場(chǎng)景
[0124] 具體實(shí)驗(yàn)場(chǎng)景可參考圖1，其中，輻射源的高度為200m，初始基線距離為13Km，發(fā)射 信號(hào)的頻率為88MHz、帶寬為100kHz ;雷達(dá)的采樣頻率為200kHz，積累時(shí)間為1秒；同頻源發(fā) 射的同頻信號(hào)的載頻88MHz、帶寬為100kHz;系統(tǒng)的信噪比為-12dB，干信比為20dB。同時(shí)，設(shè) 置2個(gè)待檢測(cè)目標(biāo)，其中目標(biāo)1的位置為：回波時(shí)延單元42，多普勒頻移單元43;目標(biāo)2的位置 為：回波時(shí)延單元65，多普勒頻移單元70。
[0125] 2)仿真內(nèi)容
[0126] ①直接利用雷達(dá)接收天線接收的回波信號(hào)進(jìn)行目標(biāo)檢測(cè)，并繪制檢測(cè)結(jié)果圖。
[0127] ②按照本發(fā)明實(shí)施例提供的同頻信號(hào)的消除方法獲取回波信號(hào)的3個(gè)獨(dú)立分量信 號(hào)yi、y2和y 3,進(jìn)而分別消除3個(gè)獨(dú)立分量信號(hào)71、72和73，根據(jù)消除獨(dú)立分量信號(hào)后的回波信 號(hào)進(jìn)行目標(biāo)檢測(cè)，并繪制檢測(cè)結(jié)果圖。
[0128] 3)仿真結(jié)果分析：
[0129] 圖3所示為直接利用雷達(dá)接收天線接收的回波信號(hào)進(jìn)行目標(biāo)檢測(cè)的檢測(cè)結(jié)果圖。 其中，圖3(a)為多普勒維的檢測(cè)結(jié)果圖，圖3(b)為時(shí)延維的檢測(cè)結(jié)果圖。觀察圖3(a)和圖3 (b)，容易看出，主瓣含有同頻干擾，目標(biāo)的尖峰被雜波所掩蓋，無法有效地檢測(cè)到目標(biāo)。
[0130] 圖4為利用消除獨(dú)立分量信號(hào)yi后的回波信號(hào)進(jìn)行目標(biāo)檢測(cè)的檢測(cè)結(jié)果圖。其中， 圖4(a)為多普勒維的檢測(cè)結(jié)果圖，圖4(b)為時(shí)延維的檢測(cè)結(jié)果圖。觀察圖4(a)和圖4(b)發(fā) 現(xiàn)，在回波時(shí)延單元42、多普勒頻移單元43以及回波時(shí)延單元65、多普勒頻移單元70處分別 有一個(gè)明顯尖峰，且尖峰高于雜波平臺(tái)10db左右。對(duì)比仿真實(shí)驗(yàn)場(chǎng)景設(shè)置中目標(biāo)的坐標(biāo)信 息發(fā)現(xiàn)，圖中所示的2個(gè)峰值正是仿真實(shí)驗(yàn)中設(shè)置的兩個(gè)目標(biāo)的回波的峰值。可見，根據(jù)消 除獨(dú)立分量yi后的回波信號(hào)能夠檢測(cè)得到目標(biāo)，由此可知，獨(dú)立分量信號(hào)7:為同頻干擾信 號(hào)，而消除獨(dú)立分量信號(hào)7:后的回波信號(hào)即為所求的消除同頻干擾信號(hào)后的回波信號(hào)。
[0131] 圖5為利用消除獨(dú)立分量信號(hào)y2后的回波信號(hào)進(jìn)行目標(biāo)檢測(cè)的檢測(cè)結(jié)果圖。其中， 圖5(a)為多普勒維的檢測(cè)結(jié)果圖，圖5(b)為時(shí)延維的檢測(cè)結(jié)果圖。觀察圖5(a)和圖5(b)發(fā) 現(xiàn)，目標(biāo)峰值被雜波掩蓋，無法有效檢測(cè)到目標(biāo)。由此可知，獨(dú)立分量信號(hào)y 2并非同頻干擾 信號(hào)，回波信號(hào)經(jīng)過ECA算法消除72后消除的不是同頻干擾成分，而是噪聲或目標(biāo)回波成 分。
[0132] 圖6為利用消除獨(dú)立分量信號(hào)y3后的回波信號(hào)進(jìn)行目標(biāo)檢測(cè)的檢測(cè)結(jié)果圖。其中， 圖6(a)為多普勒維的檢測(cè)結(jié)果圖，圖6(b)為時(shí)延維的檢測(cè)結(jié)果圖。觀察圖6(a)和圖6(b)發(fā) 現(xiàn)，目標(biāo)峰值被雜波掩蓋，無法有效檢測(cè)到目標(biāo)。由此可知，獨(dú)立分量信號(hào)y 3并非同頻干擾 信號(hào)，回波信號(hào)經(jīng)過ECA算法消除y3后消除的不是同頻干擾成分，而是噪聲或目標(biāo)回波成 分。
[0133] 本領(lǐng)域普通技術(shù)人員可以理解:實(shí)現(xiàn)上述方法實(shí)施例的全部或部分步驟可以通過 程序指令相關(guān)的硬件來完成，前述的程序可以存儲(chǔ)于一計(jì)算機(jī)可讀取存儲(chǔ)介質(zhì)中，該程序 在執(zhí)行時(shí)，執(zhí)行包括上述方法實(shí)施例的步驟;而前述的存儲(chǔ)介質(zhì)包括:ROM、RAM、磁碟或者光 盤等各種可以存儲(chǔ)程序代碼的介質(zhì)。
[0134]以上所述，僅為本發(fā)明的【具體實(shí)施方式】，但本發(fā)明的保護(hù)范圍并不局限于此，任何 熟悉本技術(shù)領(lǐng)域的技術(shù)人員在本發(fā)明揭露的技術(shù)范圍內(nèi)，可輕易想到變化或替換，都應(yīng)涵 蓋在本發(fā)明的保護(hù)范圍之內(nèi)。因此，本發(fā)明的保護(hù)范圍應(yīng)以所述權(quán)利要求的保護(hù)范圍為準(zhǔn)。
【主權(quán)項(xiàng)】
1. 一種外輻射源雷達(dá)同頻干擾的消除方法，其特征在于，所述外輻射源雷達(dá)設(shè)置有參 考天線和接收天線，所述參考天線指向所述外輻射源雷達(dá)的輻射源，所述接收天線指向所 述外輻射源雷達(dá)的觀測(cè)區(qū)；所述接收天線為包括M個(gè)陣元天線的均勻線陣，所述M個(gè)陣元天 線的相位中心間距為半波長，M為整數(shù)，M多3;所述方法包括如下步驟： S1:獲取所述參考天線接收的參考信號(hào)，以及，獲取所述M個(gè)陣元天線中每個(gè)陣元天線 接收的回波信號(hào)得到M個(gè)回波信號(hào)； S2:利用所述參考信號(hào)，消除M個(gè)回波信號(hào)中的直達(dá)波信號(hào)以及多徑信號(hào)，得到M個(gè)第一 回波信號(hào)，進(jìn)而利用所述M個(gè)第一回波信號(hào)構(gòu)造得到回波信號(hào)矩陣； S3:對(duì)所述回波信號(hào)矩陣進(jìn)行獨(dú)立分量分析，得到3個(gè)獨(dú)立分量信號(hào)，其中1個(gè)獨(dú)立分量 信號(hào)為同頻干擾信號(hào)； S4:利用所述3個(gè)獨(dú)立分量信號(hào)對(duì)所述回波信號(hào)矩陣進(jìn)行3次時(shí)域干擾相消，得到3個(gè)第 二回波信號(hào)；其中，每次時(shí)域干擾相消分別消除所述3個(gè)獨(dú)立分量信號(hào)中的1個(gè)獨(dú)立分量信 號(hào)，得到對(duì)應(yīng)的1個(gè)第二回波信號(hào)，每次時(shí)域干擾相消所消除的獨(dú)立分量信號(hào)不同； S5:對(duì)每個(gè)第二回波信號(hào)與所述參考信號(hào)作距離-多普勒二維相關(guān)，得到3個(gè)相關(guān)矩陣； 利用每個(gè)相關(guān)矩陣進(jìn)行目標(biāo)檢測(cè)，其中能夠獲得目標(biāo)信息的相關(guān)矩陣對(duì)應(yīng)的第二回波信號(hào) 即為所求的消除同頻干擾后的回波信號(hào)。2. 根據(jù)權(quán)利要求1所述的方法，其特征在于，步驟S2具體包括如下子步驟：S2a:分別對(duì)所述參考信號(hào)以及M個(gè)回波信號(hào)進(jìn)行分段，得到所述參考信號(hào)的B段分段信 號(hào)以及每個(gè)回波信號(hào)的B段分段信號(hào);其中，B為正整數(shù)，N為所述參考信號(hào)及所述回波信號(hào) 的信號(hào)長度； S2b:利用所述參考信號(hào)的第i段分段信號(hào) 構(gòu)造與其對(duì)應(yīng)的子空間矩陣丨其中，i取1到B之間的所有整數(shù) 值，表示參考信號(hào)的第i段分段信號(hào)分別表示將參考信號(hào)的第i段分段信號(hào)延遲1個(gè)時(shí)延單元、2個(gè)時(shí)延單元、???(p-1)個(gè) 時(shí)延單元后的信號(hào)，P = l+l，l表示回波信號(hào)中多徑干擾的時(shí)延單元個(gè)數(shù)，1彡1; S2c:利用所述子空間矩陣~構(gòu)造第一正交投影矩陣叫，進(jìn)而利 用所述第一正交投影矩陣對(duì)第j個(gè)回波信號(hào)的第i段分段信號(hào)進(jìn)行時(shí)域干擾相消，得到消除直達(dá)波信號(hào)及多徑信號(hào)后的第i段分段信號(hào) 其中，j取1到M之間的所有整數(shù)值； S2d:對(duì)所述M個(gè)回波信號(hào)的消除直達(dá)波信號(hào)以及多徑信號(hào)后的B段分段信號(hào)進(jìn)行拼接， 得到M個(gè)第一回波信號(hào)； 326:利用所述1個(gè)第一回波信號(hào)構(gòu)造得到回波信號(hào)矩陣1=隊(duì)<>^，"％心,, 1/，其 中，Lq 分別表不M個(gè)第一回波信號(hào)中的第1個(gè)第一回波信號(hào)、第2個(gè)第一回波信 號(hào)、…第M個(gè)第一回波信號(hào)。3. 根據(jù)權(quán)利要求1或2所述的方法，其特征在于，步驟S3具體包括如下子步驟： S3a:計(jì)算所述回波信號(hào)矩陣的協(xié)方差矩陣，并對(duì)所述協(xié)方差矩陣進(jìn)行特征分解得到特 征值及特征向量，利用所述特征值以及所述特征向量構(gòu)造得到白化向量V =『1/2ET;其中，E = (ei，"_，en)，D = diag(di,…，dn)，ei,…，en為所述協(xié)方差矩陣的特征向量，di，"_，d n為所述 協(xié)方差矩陣的特征值，（?）_1/2表示對(duì)（?）中的元素開平方后對(duì)矩陣求逆，T表示轉(zhuǎn)置； S3b:利用所述白化向量V對(duì)所述回波信號(hào)矩陣進(jìn)行白化預(yù)處理，得到白化信號(hào) =VSs,ir;其中，Ssur表示回波信號(hào)矩陣； S3c:將待估計(jì)的獨(dú)立成分個(gè)數(shù)設(shè)置為3，作為快速獨(dú) 立分量分析FastICA算法中不動(dòng)點(diǎn)迭代的基本公式，迭代獲得3個(gè)權(quán)矢量wPi、wP2和w P3; S3d:將所述3個(gè)權(quán)矢量wpl、wp2和wp3分別與所述白化信號(hào)I相乘，得到3個(gè)獨(dú)立分量信號(hào) yi、y2 和 y3〇4. 根據(jù)權(quán)利要求1-3任一項(xiàng)所述的方法，其特征在于，步驟S4具體包括如下子步驟： S4a:利用所述3個(gè)獨(dú)立分量信號(hào)中的第k個(gè)獨(dú)立分量信號(hào)yk，構(gòu)造得到第二正交投影矩 陣wk=(ykHyk)-VkH;其中，k分別取1、2以及3; S4b:利用所述第二正交投影矩陣wk對(duì)所述回波信號(hào)矩陣進(jìn)行時(shí)域干擾相消，以消除所 述第k個(gè)獨(dú)立分量信號(hào)yk，獲得第二回波信號(hào)，式中，S表示回波信號(hào)矩 陣。5. 根據(jù)權(quán)利要求1-4任一項(xiàng)所述的方法，其特征在于，步驟S5中，所述對(duì)每個(gè)第二回波 信號(hào)與所述參考信號(hào)作距離-多普勒二維相關(guān)，得到3個(gè)相關(guān)矩陣，包括： 利用預(yù)存的計(jì)算式：，對(duì)每個(gè)第二回波信號(hào)與所述 參考信號(hào)作距離-多普勒二維相關(guān)，得到3個(gè)相關(guān)矩陣|1[14]、|2[14]以及|3[14];其中沽 分別取1、2以及3，〇]表示第k個(gè)第二回波信號(hào)，Sref [n]表示參考信號(hào)，Sref [n-1 ]表示將 參考信號(hào)Sref [n ]延遲1個(gè)時(shí)延單元后的信號(hào)，n表示采樣次數(shù)，1表示時(shí)延單元，*表示共輒運(yùn) 算，p表示多普勒頻移單元，N表示信號(hào)長度。
【文檔編號(hào)】G01S7/41GK106054143SQ201610542988
【公開日】2016年10月26日
【申請(qǐng)日】2016年7月12日
【發(fā)明人】王俊, 郭帥, 鄧亞琦
【申請(qǐng)人】西安電子科技大學(xué)