專利名稱：判斷涉及固定或移動發(fā)射機、反射器和傳感器的傳播信道的脈沖響應的方法和裝置的制作方法
判斷涉及固定或移動發(fā)射機、反射器和傳感器的傳播信道 的脈沖響應的方法和裝置本發(fā)明涉及一種借助于組合式空間/延遲-距離/多普勒運動學模糊函數而能夠 在電磁波或聲波的傳播信道涉及固定或移動發(fā)射器、傳感器和反射器時評估傳播信道脈沖 響應的特定參數的方法，用于探測和確定發(fā)射器和反射器的位置和運動學參數。這種方法用于電磁學和聲學很多領域，用于探測、傳輸、定位、導航中，以便改善對 傳播介質的認識，從而改善對有用信號的處理、無線電通信和廣播無線電網絡的工程設計、 航空控制、海防控制等?？梢詫⑵鋺糜趯嵤┍景l(fā)明的接收系統(tǒng)的自我定位。例如，將其用于電磁學或聲學中固定或移動的發(fā)射器、反射器和傳感器。常規(guī)電磁學或聲學反射器的傳播、發(fā)聲、探測和定位測量系統(tǒng)時常是有源的，一般 使用-利用定向天線(例如盤式反射器)的機械掃略或利用從相位和幅度加權傳感器 陣列形成的射束的電子掃略，以及-對于射束的每個瞄準位置(或空間小區(qū))〇已知信號的發(fā)射，已知信號可以是連續(xù)的、脈沖的或已知脈沖串的形式，〇對于窄帶信號而言，基于在天線輸出處觀察到且在時間與頻率上發(fā)生偏移的信 號與被發(fā)射信號的相關性，對于寬帶信號而言，基于在天線輸出處觀察到的延遲的、在頻率 上發(fā)生偏移的、在時間和頻率上被壓縮的信號與被發(fā)射信號的相關性，對距離/速度模糊 函數的計算，-對于每個距離/多普勒小區(qū)〇針對給定誤報警概率的閾值設置，〇模糊函數與閾值的比較。這種操作的主要特別特征是空間分析(射束掃略)和距離/速度或延遲/多普勒 分析之間的解耦。無論射束中是否有發(fā)射器或反射器，這種解耦都需要針對每個射束位置 系統(tǒng)實施距離/速度分析。此外，瞄準射束隱含地預示著自由空間(沒有多路徑)中的傳播，對于電子掃略而 言，需要控制天線方向圖(輻射元件和網絡的型號、傳感器陣列的校準等)。而且，對于給定的空間/距離/多普勒小區(qū)而言，調節(jié)閾值需要事先基于無反射器 的空間/距離/多普勒小區(qū)的觀測估計噪聲水平，這可能是難以實施的，且在計算上成本高 昂。此外，在存在干擾時，必需要針對被掃略射束的每個位置在接收中插入通過空間濾波排 除干擾的技術，于是變得抗干擾且相對于干擾魯棒。然而，由于被掃射束(在波長數量方 面)具有與天線或網絡的孔徑成反比的一定角寬度，針對給定空間小區(qū)濾去干擾可能會伴 隨著濾去相同這一小區(qū)中存在的反射回波。因此，必需要在處理操作中并入防止濾去待測 量反射器的增強魯棒性的技術，代價是可能會損失濾去干擾的性能并增大實施的復雜性。測量傳播或電磁或聲學探測的問題涉及到在一定期間0≤k≤K-I內探測發(fā)射信號8(1^^)的存在并從被發(fā)射信號的知識和在傳感器上接收的信號的該矢量x(kTe)頻移偏移和平移版本的觀測估計信道矢量hs (涉及多傳感器接收并對應于自由空間中傳播的反射器或發(fā)射器位置的指引矢量)、延遲IoTe (為簡單起見假設為取樣周期的倍數，但這決不 是強制的或限制性的)和多普勒頻移△ fo = mo/KTe (頻率分辨率為Ι/KTe，為簡單起見假 設多普勒頻移是這個分辨率的倍數)。對于電磁學或聲學中的常規(guī)傳播測量或探測應用而言，常規(guī)接收器假設自由空間 中的傳播，即~ ^ ei如S,并以對應于所用網絡形成的射束的波瓣寬度(通常是以“3dB” 的波瓣寬度)分辨率逐個方向或逐個矢量s地掃描空間，其中hs是與反射器方向相關聯(lián)的 信道的脈沖響應矢量，68和s分別對應于發(fā)射器或反射器的相位和指引矢量。這定義了 本領域的技術人員通常使用的空間小區(qū)的上述概念。此外，利用等于Ι/Be的分辨率估計延 遲loTe，標稱地這是所發(fā)射信號的等價頻帶Be的倒數的函數(一般正比于該倒數)，該倒 數還定義了距離分辨率。這定義了上文所述并被本領域的技術人員通常使用的距離小區(qū)的 概念。最后，多普勒頻移的估計分辨率是(通常正比于)個體觀測時間的倒數的函數，即所 發(fā)射信號的時長KTe的倒數的函數。這定義了上文所述并被本領域的技術人員通常使用的 多普勒小區(qū)的概念?，F有技術定義了不同的接收結構。通常，最佳探測器的接收結構取決于源自要探 測的發(fā)射器和反射器的信號傳播信道上先驗可用的信息和總噪聲，總噪聲包括接收器的熱 噪聲和潛在干擾[1]。電磁[2]或聲學[4]探測中使用的常規(guī)接收器通過射束的電子或機 械掃略來掃描空間并實施位于距離/多普勒分析上游且與之解耦的空間分析，通常隱含地 或明示地假設，針對射束的每個位置和掃描的每個距離/多普勒小區(qū)_在自由空間中傳播，-掃描方向已知，-總體高斯噪聲是圓形的且未知，-相對于背景噪聲，源自發(fā)射器或反射器的信號是微弱的，-源自發(fā)射器或反射器的信號相位未知。僅在這些假設之下，這些接收器才是最佳的。本發(fā)明的目的是利用能夠至少克服 上述系統(tǒng)缺點的接收結構取代以上常規(guī)結構。這種結構尤其在于實現空間分析的耦合，以 及組合過程中的延遲_距離/多普勒_動力學分析。這種結構實施延遲_距離空間變量的 耦合/組合處理。應當指出，本發(fā)明的方法能夠確定脈沖響應的參數。這些參數例如是涉 及無線電電場的空間、時間和頻率結構(到達角分布、角度、時間和多普勒擴散等)的參數。 這些參數也可以是地面上、空中等處的障礙物漫射的特征，甚至是漫射體的運動學參數。本 發(fā)明涉及一種在包括N個傳感器的系統(tǒng)中探測和/或定位固定或移動的發(fā)射器和反射器的 方法，傳感器接收源自固定或移動的所述發(fā)射器或反射器的信號，其特征在于至少包括以 下步驟〇通過耦合空間分析和延遲_距離/多普勒_運動學分析確定模糊函數，〇確定與對應于所發(fā)射信號的復包絡的已知信號s(kTe)和濾波器的輸出w(l，m)
之間的相關對應的至少一個充分統(tǒng)計量 (I m, K),其中ι是時域/延遲-距離假設的符 號表示，且其中m是頻率/多普勒-運動學假設的符號表示，通過將針對假設(l，m)的參數 的統(tǒng)計值與閾值比較確定假設(l，m)的參數值。
多普勒/運動學參數引用m作為變量的索引。應當指出，針對統(tǒng)計值 (I m> K), 針對信號s(kTe)和濾波器w(l，m)的以上表示是以形式的意義來解釋的完全一般性地， 參數1、K和m的每個都能夠代表與發(fā)射器、反射器、甚至傳感器(如果其自身是移動的)的 位置和運動學量相關聯(lián)的若干變量。然而，對于窄帶信號以及小到注意確保減小發(fā)射器、反 射器和傳感器位置和運動學變化的過程相干積分時間而言，回到上述簡化情況將時域假 設1簡化到關于延遲值的假設，將頻率假設m簡化到關于多普勒頻移值的假設，延遲值和多 普勒頻移值分別對應于上述分辨率的倍數(整數或其他)。為了確保在大多數復雜情況下 該方法有效，例如在傳感器、發(fā)射器或反射器的運動學性質在實施中的信號相干積分時間 期間誘發(fā)延遲或多普勒參數的不固定性時(例如與快速相關的不固定性)，1的定義不僅考 慮到源自反射器的信號延遲，而且考慮到源自反射器的信號的一階參數，例如相對時域斜 變(或延遲在積分時間上的導數，其自身取決于信號的時域包絡和作為本發(fā)明主題的處理 操作的參數化)，不僅考慮到源自反射器的信號的多普勒頻移m，而且考慮到源自反射器的 信號的一階參數，例如信號的多普勒斜變(或多普勒頻移在積分時間上的導數，其自身受 限于信號的頻譜包絡和作為本發(fā)明主題的處理操作的參數化)。在更為復雜的情況下，可以 使用二階相對角度、時域和多普勒參數以及 對發(fā)射器、反射器或傳感器的運動學性質建模 的任何其他類型的變量。在將1表示的時域參數簡化為延遲或距離的情況下，或在將m表示的運動學參數 簡化為多普勒或相對速度時，該方法至少包括以下步驟對于每種延遲假設1和每種多普勒假設m _利用所發(fā)射的信號的K個已知樣本s(kTe)以及在頻率上偏移的K個觀測值 xm((k+l)Te) = χ ((k+1) Te) exp [-j2 π m (k+1) /K], 0 彡 k 彡 K_l，估計矢量Zi jra^Ure) (5)，-從頻率上偏移的所述K 個觀測值 xm ((k+1) Te) = χ ((k+1) Te) exp [-j2 Jim (k+1)/ κ]估計矩陣及 κ-1 ；-計算所述充分統(tǒng)計量t(l, m, K) ^ (1/ %) rxmS(lTe)^ Rxm(ITe)-1 rXmS(lTe).其中‘是由所述已知信號的
平均功率的最小二乘方估計定義的，-針對給定誤報警概率設置所述探測閾值，-將所述統(tǒng)計值 (1，m, K)與所述固定閾值比較，如果超過所述閾值，要求探測距 離/多普勒小區(qū)(l，m)中的發(fā)射器或反射器；-從1估計所述發(fā)射器或反射器的相對延遲和距離，-從m估計所述發(fā)射器或反射器的相對多普勒和速度；_從由》s = fjszasGT^給出的源自所述發(fā)射器或反射器的信號天線上的指引 矢量》估計所述發(fā)射器或反射器的空間小區(qū)。在更復雜的情況下，該方法例如包括以下步驟3A)將時域參數引入諸如由1符號表示的延遲-位置的定義、引入諸如由m符號表 示的多普勒_運動學參數的頻率參數的定義，并引入每個假設1和m上的所述已知信號s 的定義，信號的變化和變換由所述發(fā)射器、接收器或反射器的運動引起，由傳播信道的波動引起，這取決于延遲-位置假設、相對多普勒-速度假設、處理積分時間、信號的等價帶寬、 當在處理積分時間期間信號是脈沖式時的信號的等價時長及其周期、延遲隨這一時間的波 動、多普勒效應隨這一時間的波動，3B)根據由1符號表示的時域假設和由m符號表示的頻率假設使用從s的變換獲 得的信號S1,m。對于信號受到高于擾的發(fā)射器或反射器的探測而言，該方法使用信號X’ m((k+l) Te)和仿樣而不是輸入信號xm((k+l)Te和原始仿樣，信號x’ m((k+l)Te))是從信號 xffl((k+l)Te)獲得的，仿樣是借助于涉及增強空時濾波器的操作從變換的仿樣S1,m獲得的， 所述操作如下6A)對于每個干擾發(fā)射器通過使準則L(U)I x((lo+l)Te) - Siattsint j 2最小化估計從發(fā)射器到接
厶 ο 收器的傳播濾波器的脈沖響應 gint，仿樣 Sint = (sint((l0+l)Te), sint((l0+l+l)Te),..., 、((10+1+卜1)1^))1或5化“1 ” = (sintm, J((Vl)Te), Sintnbl“ ((10+1+1) Te),…，Sintffl, 廣((lfl+L-DTe))1是通過原始干擾信號的總先驗知識，或通過原始干擾信號的部分先驗 知識，或借助本領域技術人員已知的任何方式通過原始干擾信號的解調/均等獲得的，并 可能被轉換成Sin^16B)對于固定的或運動學已知且根據它們的頻率和時間穩(wěn)定性選擇的一個或多個 干擾發(fā)射器而言，向仿樣Sint中引入不同的頻率假設，并采用仿樣中存在的重現，以便在大的積分長 度K上精確估計相對于干擾發(fā)射器的頻率和時間偏離，針對用于接收和取樣用于實施本發(fā) 明的信號的系統(tǒng)的混頻器和時鐘的本地振蕩器估計相對于干擾發(fā)射器的漂移。在步驟6B期間，該方法例如包括以下步驟更新與干擾發(fā)射器相關聯(lián)的傳播濾波器gint (步驟6A)，以及6C/對于每個干擾發(fā)射器，由與其關聯(lián)的傳播濾波器gint執(zhí)行原始Sint或變換的Sintnbl干擾信號的卷積，從總 的所接收信號減去這一操作獲得的信號，6D/補償頻率和時間漂移是-在前面操作之后從xm((k+l)Te)獲得的信號，以便產生信號χ’m((k+l)Te)，s然 后保持不變(S’ = S)_或者在該方法中使用的仿樣信號，s于是可以變?yōu)镾’ ^nbl可以變?yōu)镾nbl’6E/以X，m和S，實施上述步驟6F/以χ，m和sm, ！，實施上述步驟。使用觀測矢量？歷((A+J)Te) - [Xb( (^-i) Te)，Xm{ (J^1) Te)卞],例如觀測
矢量由觀測矢量及其共軛復數構成，且從7/ (O+7) Te)而非xm((k+l)Te)實施上述步驟。例如，也可能利用經由僅噪聲參考樣本獲得的擴展總僅噪聲矩陣&7> (""e)， 而非擴展觀測相關矩陣的估計i^w(/7V)，來執(zhí)行總僅噪聲參考和上述步驟，其中<formula>formula see original document page 11</formula>
由于對于直線發(fā)射信號和可能非圓干擾而言，反射器的指引矢量是已知的，因此 該方法的特征在于其至少包括以下步驟計算充分統(tǒng)計量<formula>formula see original document page 11</formula>
構造矩陣<formula>formula see original document page 11</formula>構造統(tǒng)計量
<formula>formula see original document page 11</formula>(20)
hm^ RxmilTeT1i M [M^ R^m(ITe)-1M 廠1 M^ ^(/^Γ1 ^ms(ITe)
-tnc(l’m,K)針對給定誤報警概率設置探測閾值，將統(tǒng)計值W, S, Λ)與閾值比較，如果超過閾值，則探測空間小區(qū)s和距 離/多普勒小區(qū)(l，m)中的發(fā)射器或反射器，從1估計發(fā)射器或反射器的相對延遲和距離，從m估計發(fā)射器或反射器的相對多普勒和速度/運動學性質，從s估計發(fā)射器或反射器的角位置。本發(fā)明還涉及一種包括適于執(zhí)行本發(fā)明方法 的步驟的諸如處理器的模塊的裝置。作為本發(fā)明主題的空間/延遲_距離/多普勒_運動學耦合使得能夠直接訪問發(fā) 射器和反射器的位置而無需掃描所有空間，因此無需掃描沒有發(fā)射器和反射器的方向。這 樣能夠放寬對傳感器網絡模式的控制的限制(例如，對于輻射元件之間的耦合而言，校準 不再是必要的)，完美地適合于多路徑傳播。最后，明智的耦合不僅允許獨立于輸入處存在 的噪聲水平進行閾值設置，而且允許相對于干擾進行固有保護，而沒有拒絕對測量有用的 信號的風險，且無需插入處理操作以提高魯棒性。閱讀通過例示而非限制方式給出的簡化示范性實施例的以下描述，本發(fā)明的其他 優(yōu)點和特征將變得更明顯，附圖提供了 〇

圖1，根據本發(fā)明的空間/延遲_距離/多普勒_運動學模糊函數的工作原理 圖，在這種非限制性簡化情況下將多普勒-運動學參數簡化成直接鏈接到發(fā)射器或反射器 的相對速率的多普勒參數。〇圖2，考慮干擾的方法的示范性實施，〇圖3，空間/延遲-距離/多普勒_運動學模糊函數的非圓擴展的工作原理圖， 在這種非限制性簡化情況下將多普勒_運動學參數簡化成直接鏈接到發(fā)射器或反射器的相對速率的多普勒參數，以及〇圖4，在s已知時的空間/延遲-距離/多普勒_運動學模糊函數的工作原理圖，在這種非限制性簡化情況下將多普勒_運動學參數簡化成直接鏈接到發(fā)射器或反射器 的相對速率的多普勒參數?？紤]具有N個窄帶(NB)數字傳感器Ci的天線，數字傳感器從源自已知時機的電 磁或聲發(fā)射系統(tǒng)的連續(xù)、脈沖或重復脈沖信號的發(fā)射器或反射器接收貢獻(非限制性范 例預先解調的通信發(fā)射器，或其信號的特定部分是先驗已知的-同步詞和代碼-信道發(fā)聲 器的發(fā)射部分)，向其添加由干擾和背景噪聲構成的總噪聲。給定這些假設，假設在傳感器 Ci輸出處采樣和觀測到的信號的復包絡的矢量x(kTe)由下式給出<formula>formula see original document page 12</formula>(1)其中Te是數字傳感器的采樣周期，s (t)是所發(fā)射信號的復包絡，為接收器所知并 是KTe的時長，bT (kTe)是采樣的總噪聲矢量，不與源自反射器的信號相關，hs是與源自反 射器的信號相關聯(lián)的信道的脈沖響應矢量，IoTe是源自反射器的信號的傳輸延遲，具體而 言包含關于反射器、使用的發(fā)射系統(tǒng)和傳感器之間的距離的信息、為簡單起見假設等于多 個Te的延遲，作為非限制性范例給出的假設，Afo = mo/KTe是反射器(可能是移動的)導 致的多普勒頻移，為簡單起見假設等于處理操作KTe的積分時間倒數的倍數，作為非限制 性范例給出的假設，其包含反射器相對于發(fā)射系統(tǒng)和所用傳感器位置的相對速度的信息。 嚴格來說，模型(1)表示對于例如自由空間中傳播出現的沒有時域擴散的傳播信道(衛(wèi)星 電信、空中無線電通信等)或具有信號衰落的信道，對于具有信號衰落的信道而言，信道的 相干頻帶大于信號的帶寬，更多地稱為術語“平坦衰落”，這影響到城區(qū)中的特定無線電通 信或探測環(huán)境。具體而言，對于自由空間中的傳播，As ^ e^h,其中Φ8和s分別對應于源 自發(fā)射器或反射器的信號的相位和指引矢量。模型(1)還表示，將發(fā)射器或反射器的時域 和頻率參數簡化為延遲和多普勒頻移，僅僅對于特定運動學性質(勻速直線運動以及相對 于過程的相干積分時間的有限速度)的發(fā)射器或反射器，本說明書才是嚴格完整的。然而，該方法還適用于時域上擴散的信道。在這些狀況下，(1)的有用部分實際上 對應于源自發(fā)射器或反射器的信號的路徑或傳播模式的貢獻，其他路徑集成在總噪聲矢量 中。此外，如前面強調的，該方法還適用于表現出的運動學性質比均勻且相對于過程 積分時間而言速度有限的直線運動更復雜的發(fā)射器、反射器或傳感器。在這些狀況下，(1) 的有用部分實際上對應于經歷了這種運動學性質誘發(fā)的變形的信號的貢獻以及必需要考 慮與如前所述的一階參數(1階)互補的有用信號的觀測值(例如時域斜變和多普勒斜 變)，以便報告高的相對速度或運動，這時在作為本發(fā)明主題的過程相干積分階段期間發(fā)射 器、反射器或傳感器的加速度有限；甚至與二階參數(2階)互補，以便報告作為本發(fā)明主題 的過程的相干積分期間發(fā)射器、反射器或傳感器的顯著加速/減速階段；最后在將該方法 應用于通過這種方式通過s的變換獲得的信號S1,m之前對1表示的時域假設和m表示的頻 率假設誘發(fā)的信號s的變形進行建模。注意，上述參數列表不是限制性的，可以引入其他參 數以對發(fā)射系統(tǒng)、傳感器或反射器的運動對源自后者的信號的影響建模。二階觀測統(tǒng)計量在前述假設的條件下，在說明書下文中考慮的二階統(tǒng)計量對應于在下式定義的傳感器輸出x(kTe)處采樣和觀測的信號復包絡矢量的第一相關矩陣<formula>formula see original document page 13</formula>這通常取決于時間，因為所發(fā)射信號可能是不固定的(例如，用于通信信號的數
字調制……)。在以上表達中，卞代表共軛轉置，對好e) -ElbTikTe) ATUTfe)卞]是矢
量bT(kTe)的第一相關矩陣， (ATe) A EfkW^I2]是將由用于自由空間傳播的全向傳感 器接收的源自發(fā)射器或反射器的信號的瞬時功率。例如，在包括幾個接收天線Cp……Cn和處理器P的接收器中實施下文所述的各 實施例，處理器被編程控制來執(zhí)行各步驟。結合圖1、2、3和4給出各示范性實施例。本發(fā)明的方法尤其依賴于在組合過程中利用空間分析和延遲-距離/多普勒-運 動學分析的耦合來確定模糊函數。本發(fā)明的示范性實施例圖1所示的本發(fā)明方法的第一變體實施例包括找到時域/延遲-距離假設 /=/;和頻率/多普勒-運動學假設W=扎，從而從觀測矢量Xm((k+i)ig = x((k+i)Te) exp [-j2 Ji m (k+1) /K] (0彡k彡K_l)最佳地探測到已知樣本s (kTe)，0彡k彡K-1，假設R (k) 和hs是未知的，并且對傳播信道和源自發(fā)射器或反射器的信號電平不做假設。一種方法包括找到/二&和講=；^，使得從最小二乘方的意義上，從對觀測值 Xm ((k+1) Te) (0彡k彡K-1)的線性濾波操作最佳地估計已知樣本S (kTe)，0彡k彡K-I。換
言之，該方法試圖找到使最小二乘方基準最小化的，
Λ1 κ~1Λ2
t(l,m,K) ^ — Σ I s(kTe) - w{hm^xm{{i&t)Te)\(3)
K 4 = 0其中圖1中引用的A (7, /77)，是使準則(3)最小化且由下式定義的空間濾波器in) ^ Rxm(ITe)'1 rXmS(lTe)(4)其中矢量i5α 7b)和矩陣左xmU Te)由下式給出
<formula>formula see original document page 13</formula>將(4)、(5)和(6)代入(3)，由下式給出最小二乘方準則S {1, m, Λ)的新表達<formula>formula see original document page 13</formula>(7)其中圖1中標示為4的0 Qj m>幻使得0 S C (I, m> JO < 1,由下式給出
<formula>formula see original document page 13</formula>
且彥S由下式定義<formula>formula see original document page 14</formula>從(7)導出，/0和而0分別是使(S)給出的充分統(tǒng)計量0 (J’m，幻最大化的參數1和 m0在實踐中，為探測閾值設置值β (圖1中的5)，使得在沒有源自發(fā)射器或反射器的 信號時超過該閾值的統(tǒng)計量(8)等于被稱為誤報警概率(FAP)的特定值。因此，如果統(tǒng)計 量(8)超過對(l，m)的閾值β，針對距離小區(qū)1和多普勒小區(qū)m探測發(fā)射器或反射器。該接收器不要求對傳播信道、干擾或所發(fā)射信號做任何假設，因為后者是部分或 全部已知的(例如，在完成在先解調時)。該接收器不要求掃略空間。此外，它獨立于觀測矢 量的功率，這意味著探測閾值獨立于輸入處存在的噪聲水平，因為在K大于幾個單位或幾 十單位時，統(tǒng)計量概率(？ (I m, Λ)定律大致遵循高斯定律且不再受到干擾定律的影響(無
論它們的水平如何)。類似性比率方面的最佳性條件從最佳探測時刻IoTe和最佳多普勒頻移Δ f。= mo/KTe開始，并假設具有兩個假 設HO和Hl的情形，HO 僅在xmo((k+lo)Te)中存在總噪聲，Hl 在xmo ((k+lo) Te)中存在總 噪聲和有用信號，可以如下表達Hl :xm。((k+l。)Te) s(kTe)hs+bTm。((k+l。)Te)(IOa)HO :xm。((k+l。)Te) bTm。((k+l。)Te)(IOb)其中bTm。((k+l。)Te) = bT ((k+l0) Te) exp [-j2 π m0 (k+l0) /K],在這一環(huán)境下，根據探測的統(tǒng)計理論(Neyman Pearson之后)[1]，基于所發(fā)射信 號期間的觀測值Xm。( (k+lo) Te)探測源自發(fā)射器或反射器的信號s (kTe)的最佳策略包括 將類似性比率(RV)，L(X) (lo, mo)與閾值進行比較，類似性比率被定義為
_]職^ P[W(明0“仏-口叫(11)
V[xmomio)Te\ 0 ^ 1, / H0]其中p[xm。((k+l。)Te)，0彡k彡K-l/Hi] (i = 0,1)是在假設Hi的條件下，矢量 [Xffl0(I0Te)1, Xmo((1+10)Te)τ, ... , xm。((K+l。-l)Te)T]T 的條件概率密度。然后能夠表明，如果 符合以下條件，在一般化類似性比率的意義上，由(8)表征的接收器是對源自發(fā)射器或反 射器的信號進行最佳探測的問題的充分統(tǒng)計量Al 樣本 bTmo((k+lo)Te)彼此不相關，0 彡 k 彡 K-1 ；A2 矩陣 R((k+lo)Te)不取決于 k，A3 矩陣 R((k+lo)Te)未知，A4:矢量hs未知，A5 矢量 bTmo((k+lo)Te)是高斯型的，0 彡 k 彡 K-1A6 矢量bTmo((k+lo)Te)在二階是圓形的，0彡k彡K-1統(tǒng)計量(8)還可以被寫成標示為4的形式<formula>formula see original document page 14</formula>
其中;會i</，(“十是用于觀測矢量xm((k+l)Te)的濾波器
w (I歷)的輸出。于是，從(I2)導出，在歸一化因子之內，充分統(tǒng)計量6 (又馮幻，4對應于 已知信號s(kTe)和濾波器& (1,范)的輸出ym((k+l)Te)之間的相關3的結果，如圖3所示 在最小二乘方的意義上，濾波器A (1, Λ7)是在K—c 時獲得的在K個觀測矢量上 對濾波器w(/，m) ^及JT1Gws(Zre)的估計，其使得s(kTe)和輸出對皿((如乃Te)之
間的平均二次誤差在無限觀測時界上的時域平均最小化，其中Κχ < E[Xm((k+l)Te) xm{{k+l)Te)^>, rXmS(lTe) ^ <E[xm({Hl)Te) s(kTef] >
且其中 <. > 是在無限觀測時界上以k為單位的時域平均操作。只要1和m分別與Io和 mo保持遠離，濾波器w(l，m)就保持接近零，這樣產生也接近零的統(tǒng)計值(？(人/H, Λ)， 由于所發(fā)射信號的時長有限，在估計噪聲之內。在1和m分別接近Io和mo時， C(l, m’K)么^00 C(l, m, /Q增大并在1 = 1。和m = m。時達到其最大值。在這種情況下， 所發(fā)射信號s(kTe)和觀測矢量Xm。((k+lo)Te)的有用部分同相，假設總噪聲和所發(fā)射信號 之間統(tǒng)計獨立，濾波器w (lo，mo)對應于以下空間適配濾波器(SAF)
m0)厶 Rx-1 rXmoS(l0Te) =[π8/(1 +π8JT1 hs (13)S卩，對應于使輸出處的信號與干擾加噪聲比(SINR)最大化的空間濾波器，其中Tis = < E[\s(kTe)\2] > and R = < E[bj{kTe) bj{kTe)^] >。因此，C(l。，m。，K)被
寫為
<formula>formula see original document page 15</formula>其中SINRy是在SAF，w(lo，mo)的輸出處的SINR，被定義為
SINRy = %s hs^R~l hs(15)函數C(l0，m0，K)是SINRy的遞增函數，在輸出SINRy處信號與干擾加噪聲比值高 時接近一。只要干擾的數量保持小于傳感器的數量，SAF (13)就拒絕干擾，同時使輸出處的 信噪比(SNR)最大化。在這些條件下，所發(fā)射信號和SAF估計的輸出之間的相關操作從有 限數量的樣本K產生最大相關，樣本K的最小值因此應當在輸出處的SINR減小時增大。于是，對于每個被掃描距離/多普勒小區(qū)而言，計算關聯(lián)的SAF使得能夠使竊聽路 徑輸出處的SINR最大化，而對于常規(guī)接收器而言，空間濾波器對于同一個空間小區(qū)的所有 距離/多普勒小區(qū)都是恒定的。說明時域假設或提供探測的距離小區(qū)給出關于從發(fā)射器到接收器的相對延遲/距離 (如果發(fā)射器上有同步信息可用，也關于絕對延遲/距離)的信息或關于反射器相對于發(fā)射器和接收器的距離的信息。頻率假設或提供探測的多普勒小區(qū)給出關于發(fā)射器相對于接收器的相對速度或 反射器相對于發(fā)射器和接收器的位置和運動的速度的信息。與探測相關聯(lián)的空間小區(qū)被定義為^ =左Xww(Zo^)W=Zixws(Zre),能夠想象
提供探測以及對發(fā)射器和反射器的同時表征，而對它們的位置沒有先驗了解，且沒有空間 掃略。提出的處理操作并非唯一針對使用固定或旋轉定向天線，無論是在被操作發(fā)射系 統(tǒng)的水平還是在作為本發(fā)明主題的傳感器和方法的水平上，這樣能夠實現逐扇區(qū)地搜索并 針對噪聲(改善靈敏度)提供額外的天線增益和對干擾的補充空間濾波。即使發(fā)射系統(tǒng)甚至在同一頻帶中同時發(fā)射不同信號(例如“彩色”發(fā)射)，也可以 想到本發(fā)明提出的處理。對該方法進行空時擴展以在存在干擾和多傳播路徑且存在接收系統(tǒng)故障時改善 性能在有大量干擾時，例如在存在源自與多個傳播路徑相關聯(lián)的多個發(fā)射器的信號的 情況下，要處理的信號源自反射器時，或者在使用的接收系統(tǒng)呈現出一定的不穩(wěn)定或靈敏 度降低故障時，能夠增加該方法中自由度的數量，因此改善其分離能力、其區(qū)分低速發(fā)射器 和反射器的能力以及其靈敏度。上文所述的步驟使用對所接收總信號的時域預過濾，對于 每個干擾發(fā)射器，所述步驟包括評估干擾發(fā)射器的總體貢獻_這種貢獻源自由干擾源、由 與其相關聯(lián)的傳播濾波器發(fā)射的信號的卷積，這種信號和這種濾波器是預先估計的_，然 后，從輸入信號減去干擾發(fā)射器導致的總貢獻。在下文中通過例示和非限制性方式給出了 簡化的示范性實施例，附圖2表示與作為本擴展主題的空時預濾波相關聯(lián)的前面所述的一 般方法的實施。例如，強干擾對應于值比要探測的接收器大10到40dB的干擾。為了對信號進行預過濾，必需要估計干擾源發(fā)射的信號Sint (kTe)以及與其相關聯(lián) 的信號傳播信道gint (ITe)。圖2示出了 Sint的不同構造可能。為此，根據以下情況可以實施不同方法 在干擾源例如是(以非限制性方式)信道發(fā)聲器、預先調制的無線電通信或無 線電廣播發(fā)射器、雷達或有源聲納時，在原來發(fā)射的干擾信號Sint完全已知時 從干擾發(fā)射器向接收器傳播信號誘發(fā)該方法通過本領域的技術人員已知的任 何估計方法評估濾波器，例如(以非限制性方式)，方式是搜索在時刻Itl使下面的時域最 小二乘方準貝站int’(/o，L)最小化的傳播濾波器gint = (gintl，gint2，...，gintL)的系數gintQ = gint (I0Te)，gintl = gint((l0+l).Te), ···，gintL = gint((l。+L).Te)，根據在適當選擇的天線上 所取的信號X并根據在L個時域樣本上所取的仿樣Sint= (sint((l0+l)Te), sint((l0+l+l) Te)，. . .，Sint ((lQ+l+L-l)Te))T 評估所述系數。
<formula>formula see original document page 16</formula>適當恢復的濾波器的脈沖響應gint包含從發(fā)射器到接收器的所有傳播路徑。在部分知道原始發(fā)射的干擾信號Sint時，在干擾源例如是標準化的或已知技術規(guī)范主題的通信發(fā)射器時通過本領域的技術人員已知的方法估計用于接收到的干擾信號 的平衡濾波器；如果存在信號中先驗已知的序列，通過對該序列進行時間和頻率同步，并應 用本領域技術人員已知的估計gint的任何方法，例如(非限制性方式)使全部或部分工作 于源自發(fā)射器的信號Sint的已知部分的先前準則&nt，(]0> U最小化；或者如果存在信號 中先驗已知的序列，通過對該序列進行時間和頻率同步，然后對整個干擾信號Sint進行解調 /均衡，并應用本領域技術人員已知的估計gint的任何方法，例如(非限制性方式)通過使 用所恢復干擾信號Sint的部分來使先前準則SirZ Uc L)最小化；或者，最后通過盲符號 同步技術、涉及所接收信號χ的載波恢復/解調/均等并利用源自干擾源Sint的信號上的 可用先驗知識以恢復之，然后應用本領域技術人員已知的估計gint的任何方法，例如(非限 制性方式)通過使用所恢復干擾信號Sint的部分來使先前準則U Uo, L)最小化。注意，原則上，將前述步驟的不管哪一個應用于在例如專用于對前向信號解調 (這是圖2的非限制性范例)的適當選擇的特定天線上接收的信號X，或在天線網絡上接收 的信號，可能與用于將該方法用于有用信號的信號相同在這種情況下，通過針對xGldl) Te)，gintl和Sint ((10+1) Te)，考慮到維度為NX 1的信號矢量(N為天線數量)并考慮N維復矢 量的空間中的矢量范數，要最小化的準則U {lo, L)的前述表達在形式上保持相同。在 這種情況下，獲得干擾源和傳感器之間的傳播信道的脈沖響應的空時(不再僅僅是時間) 估計。從前述步驟中的一個或其他步驟，除了源自干擾源Sint (kTe)的信號之外，還導出 干擾發(fā)射器和接收器之間的傳播濾波器gint。
此外，通過將多普勒假設引入仿樣Sint并將信號Sint的重復和重現作為基準，能夠 精確地估計振蕩器、混頻器和傳感器時鐘的相對漂移。步驟6B)在所選的干擾發(fā)射器是固定的，或位置和運動學性質已知且比傳感器更 穩(wěn)定時，在應用作為本發(fā)明主題的方法時使用估計的漂移以增加積分時間(即增大參數 K)，因此改善該方法的靈敏度、其多普勒分辨率以及其區(qū)分低速反射器或發(fā)射器的能力。最 后，將該方法的這種擴展重復到幾個干擾發(fā)射器能夠在下文所述的信號校正之前進一步細 化估計。步驟6C)通過發(fā)射器_接收器傳播濾波器的前述估計gint對來自干擾發(fā)射器的信 號^㈣進行卷積，然后從天線上接收的總信號減去這種操作得到的貢獻或結果，對有用信 號、天線網絡上的總信號分量進行預過濾。通過這種方式，通過(時域)處理消除了干擾發(fā) 射器(前向路徑+多路徑)造成的能量貢獻的大部分，在這種預過濾獲得的信號中大大增 強了源自發(fā)射器或反射器的有用信號的信號干擾比和相對貢獻。在完成這種預過濾時，然 后將本發(fā)明的方法應用于從輸入信號X獲得的信號X’，因此提高了總靈敏度和總體性能水 平。步驟6D)通過估計并補償接收器中的時鐘漂移和本地振蕩器漂移，通過這種擴展 還大大增大了在實施該方法時能夠使用的樣本數量的值K，因此增大了多普勒分辨率、區(qū)分 低速反射器的能力以及該方法的處理增益和靈敏度。在總噪聲非圓時擴展該方法圖3示出了另一示范性實施例，在干擾呈現出非圓性質時，在前述方法的擴展時實施的步驟利用了這種性質。盡管所發(fā)射信號具有非圓性質，也可以想到使用這種接收器。 最好的情況對應于所發(fā)射信號是直線性的，換言之其復包絡S (kTe)具有實數值的情況。尤 其對于調幅(AM、ASK)或2狀態(tài)調相(BPSK)發(fā)射信號就是這種情況。在這些狀況下，利用 如下定義的統(tǒng)計量(圖3中的10)替換充分統(tǒng)計量(8)是有利的
^ncih m, K) ^ (1/ %) hms^ hm^TeTl ^ms(ITe)(17)其中石^[xmmi)Te)T, ^((W^)1"]1^其中&和是 由(5)和(6)定義的，其中xffl((k+l)Te) mmmt)Te),(圖3中的6)替換，圖3為充分統(tǒng)計 量(16)的工作原理圖。將統(tǒng)計量(17)擴展到所發(fā)射信號，在退旋預處理之后，所發(fā)射信號變?yōu)闇手本€性 的。尤其對于η/2-BPSK或OQPSK型的所發(fā)射信號，甚至特定的連續(xù)相調制(CPM)，即調頻 的且具有連續(xù)相，例如MSK和GMSK信號而言，就是這種情況。在這些狀況下，必需要實施充 分統(tǒng)計量(17)，但要通過用L個時域去旋的空時矢量替換空間矢量Sm((HZ)A)來實施如果L· 為奇數，4 [3^w((HiL-iy2)re)T，…，Sidm((k - (I-l)/2)re)T]T，
以及如果 L 為偶數， [xdm((k+L/2)Te)T,..., xdmiik-L/2)Te)T]Tf其中切《((Λ+/)7ν)厶[xdmmi)Te)T,XdmmOTe)^ 'xdm{(M)Te)厶 ^(Λ+/)^)
exPH27tm(^>^’姒且τ是所發(fā)射信號的符號時間。在已知或假設反射器的指引矢量時擴展該方法標稱擴展在傳播發(fā)生于自由空間中時，有用信道矢量、被表達為/^ Aeic^s其中小3和8分
別對應于源自發(fā)射器或反射器的信號的相位和指引矢量。在考慮空間的空間掃略時(對于 常規(guī)探測結構就是這種情況)，可以認為在掃描發(fā)射器或反射器的空間小區(qū)時矢量s是已 知的。在這些狀況下，采取假設Al到Α6，但用單個未知相位Φ8的假設Α4’替換假設Α4， 常規(guī)的探測結構僅對源自低功率發(fā)射器或反射器的信號而言證明是理想的，對于源自大功 率發(fā)射器或反射器的信號變?yōu)椴荒敲蠢硐?。在這些狀況下，假設以非常低的FAP作為目標， 針對所有信號電平優(yōu)化探測性能水平可以被證明是有用的。在特定狀況下，如圖4中的13、12、14所示，用于在空間小區(qū)S、距離小區(qū)1和多普 勒小區(qū)m中的時長上探測已知信號s(kTe)的充分統(tǒng)計量由下式給出
F(/, m, s,K) = (/ %) —---x--(18)
(s^ Rxm(ITe)-1 s) (I-CihnuK)) 其中0 (l，m，K)，l2，由⑶定義。表達式(1了)是在矢量s已知時⑶的擴展。對
于低的SINRy值(SINRy << 1)，從(16)中導出 (1，m, K) 1和統(tǒng)計量(17)可以近似為
<formula>formula see original document page 19</formula>
其對應于自適應信道形成時常規(guī)使用的結構IRxm(ITe)-1 rXmS(lTe) I2, 14,
但由S"T為;^(/7^廠、項進行標準化，對于獨立于輸入噪聲或干擾功率的給定FAP而言， 這允許對探測閾值β，15進行調節(jié)。另一方面，對于不再低的SINRy值而言，具有已 知s的最佳探測器采用前面解釋的統(tǒng)計量0 (人鞏幻，因為在這種情況下空間濾波器
|^，m)會^^(/？ —1；^^/ 收斂得快于濾波器&(/，m)會吳;^^^廠、,如[3]中所示。
在圖4中示出了充分統(tǒng)計量(17)的工作原理圖。僅總噪聲參考是可用的，從經由僅噪聲參考樣品獲得的僅總噪聲矩陣的估計 ^7^/7 ，而非觀測相關矩陣的估計實施前述步驟。非圓擴展在干擾表現出非圓性質時，即使在假設s已知時，也可以使用能夠利用這種性質 的接收器。盡管可以構思實施這種接收器而不顧所發(fā)射信號的非圓性質，最佳情況還對應 于所發(fā)射信號為直線性的情況。在這些狀況下，利用如下定義的統(tǒng)計量/H, s, AO替換充
分統(tǒng)計量(17)是有利的
<formula>formula see original document page 19</formula>
<formula>formula see original document page 19</formula>
其中M是如下定義的矩陣(Nx2)<formula>formula see original document page 19</formula>(21)針對在去旋預處理之后變成準直線性的所發(fā)射信號擴展統(tǒng)計量(20)。在這些狀況 下，必需要實施充分統(tǒng)計量(20)，但要通過前面段落中定義的用L個時域獲取去旋的空時 矢量，對((A+/)re)，替換空間矢量來實施。具有僅噪聲參考的非圓擴展對于一定數量的應用而言，通過接收沒有任何有用信號的二次數據，僅噪聲加干 擾的參考是可用的。例如，在這種二次數據能夠源自與要探測的有用信號的存在相關聯(lián) 的那些距離_多普勒小區(qū)相鄰的距離_多普勒小區(qū)時，就是這種情況。在這些狀況下， 可以通過任何相關矩陣估計器從這種二次數據估計擴展的總噪聲的相關矩陣
hTm(lTe)>"從那U后，胃具W僅總、噪考白勺非圓《 ^1。如果s是卩白勺，S禾中
探測器可以對應于接收器(17)，或者如果s是已知的，對應于接收器(20)，但要通過用經由
僅噪聲參考樣本獲得的公而非擴展觀測相關矩陣的估計^^//Τ^，替換 >力 (/7^來對應，其中<formula>formula see original document page 20</formula>以同樣的方式，還可以根據同樣的方法構造具有僅總噪聲參考的接收器，該參考既使用基本數據(信號加總噪聲)又使用二次數據。針對寬帶源的擴展在接收的源(源自發(fā)射器或反射器的有用信號、干擾)對于使用的天線網絡是寬 帶時，前述處理性能水平降低，可以實施以下兩種變體實施例，在比例Byftl或π DBe/c超過 幾個百分比時發(fā)生這種情況，其中D是網絡的孔徑，Be是信號的等價帶寬(視為接收頻帶中 包括的Be)，f0為信號的載波，c是傳播速度第一種實施例包括將接收頻帶整理成幾個對于網絡而言窄的子帶，以及逐個子 帶實施前述處理操作。這種策略的優(yōu)點是其并行安排相對簡單的處理操作。第二實施例包括實施接收空間_時間結構。在對于傳感器網絡而言信號不再是窄帶的情況下(π DB/c超過幾個 百分比)前述空間矢量Xm((k+i)ig被具有L個時域獲取的空時矢量替換
XmM(^Te) [xm((k+(L-l)/2)Te)T,..xm((k - (I-l)/2)re)T]T 如果 L 為奇數，而
Xmfst((k+0Te) ~ [xdmiik + Lf2)Te)T,…，xdm((k - I/2)re)T]T如果 L 為偶數。在與載波相比信號也是寬帶的情況下(Byftl超過幾個百分比)以及發(fā)射器或反射 器的運動學性質要求考慮在仿樣信號s上誘發(fā)的形變的情況下向前述處理操作中不僅引 入延遲ITe和多普勒m/ (KTe)變量，而且引入諸如積分時間K上的延遲壓縮和多普勒壓縮的 額外一階變量，諸如前述變量的時域漂移的二階變量等，通過1表示的時域變量(延遲-位 置等)和m表示的頻率變量(多普勒_運動學性質等)寫入對應于原始信號s的變換的信 號S1,m，并使用信號S1,m來實施如前述權利要求中主張的方法的步驟并將本發(fā)明的方法應用 于該信號S1, m。這種方式提供了不向距離分辨率中弓丨入任何損失的優(yōu)點。補充處理操作可以利用前面描述的結構之一構思補充處理操作。具體而言，可以基于與距 離/多普勒小區(qū)相關聯(lián)的觀測值xm((k+i)Te)，或或基于在作 為K的倍數或約數的間隔上觀測值與所發(fā)射信號S(kig的相關結果，換言之，基于矢量 bcmsU.K'Te)’r除以κ或κ除以K’(這是參考測角的情形)構思目標在于從角度上定位 給定距離/多普勒小區(qū)的所有發(fā)射器的高分辨率測角技術。在S已知的情況下，可以在將觀測值的空間濾波限制到空間小區(qū)s之后實施這些 技術。本發(fā)明的方法例如用于位置已知的接收系統(tǒng)，以及面對同步、位置或運動學信息 在本發(fā)明使用的信號發(fā)射或反射上可用的發(fā)射器或反射器，這樣能夠從本發(fā)明基于作為本 發(fā)明主題的方法產生的相對信息推斷出絕對的時域信息(延遲/距離)和絕對頻率信息 (速度/絕對多普勒效應)。根據一個示范性實施，在彼此同步和通信或與操作中心同步和通信的幾個接收系 統(tǒng)上重復該方法的使用，以產生關于發(fā)射器和反射器的同步、位置或運動學信息，以便獲得發(fā)射器和反射器的絕對時域和頻率信息?？梢栽谝韵聭弥恢惺褂帽景l(fā)明的方法傳播信道測量(發(fā)聲)、無線電通信和 無線電廣播網絡的計量、通信發(fā)射器的解調和測角、有源和無源雷達、有源聲納、敵友詢問， 敵友詢問更多地被稱為“鑒別朋友和敵人(IFF) ”。根據變體實施例，在該方法中，例如，分別用未多普勒化(m = 0)且未時域偏移(1 =0)的觀測矢量^(1^)，3 ( ^431；)或 ^ (λ?！?，替換多普勒化和時域偏移的觀測矢 量Xm ((k+1) Te)，xmmJ)Te\ xffl,st ((k+1) Te)或省，其中參考信號是多普勒化的 且時域偏移的，即被 sm((k+1)Te) = s(kTe)exp[j2Jim(k+l)/K].替換。本發(fā)明的方法可以包括基于延遲_距離/多普勒_運動學或空間/延遲-距離/ 多普勒_運動學假設的常規(guī)測角步驟(高分辨率等)或參考測角步驟，包括以下方式實施 常規(guī)或高分辨率測角技術參基于與能量距離/多普勒(1，m)相關聯(lián)的觀測值xm((k+l)Te)，iw((W)7V)，xm, st((k+i)Te)或 或者基于在作為K的約數的時段上的觀測值Xm ((k+1) Te)，XmmDTe), xffl, st((k+l)Te)或與所發(fā)射信號 s(kTe)相關的結果， 在指引矢量s已知的情況下，在將觀測值的可能空間濾波限制到空間小區(qū)s之后。本發(fā)明的裝置和方法使得尤其能夠測量傳播，探測和/或定位包括N個傳感器的 系統(tǒng)中的固定或移動發(fā)射器或反射器，傳感器從發(fā)射器或反射器接收信號，發(fā)射器或反射 器包括至少一個模塊，使其能夠執(zhí)行該方法的步驟。根據變體實施例，發(fā)射器或反射器的指引矢量是已知的，該方法包括以下步驟。 從所發(fā)射信號的K個已知樣本s(kTe)和K個頻率偏移觀測值xm((k+1)Te)=
X ((k+1) Te) exp [-j2 JI m (k+1) /K]估計矢量、 燦(/7>)，0 彡 k 彡 K_1， 從 K 個頻率偏移觀測值 Xm ((k+1) Te) = χ ((k+1) Te) exp [-j2 Jim (k+1)/K]估計矩 陣、 (/re)，o 彡 k 彡 κ-i， 計算充分統(tǒng)計量d(l, m, K) ^ (1/ %s) rXms(lTe)^ Rxm(ITe)-1、脅(嘰)其中&是由⑶定義的， 計算充分統(tǒng)計量
0" n-n^、丨,^V)-1。m抓)丨2 Hh 叫 s，/Q = (W %) —-----λ-
(ST ^(^e)"1 S) (1 - C(/, m, K)) 針對給定誤報警概率設置探測閾值， 將統(tǒng)計值f(/，m, S,尺)與設置的探測閾值比較，如果超過該閾值，要求探測距離/ 多普勒小區(qū)(1，m)和空間小區(qū)s中的反射器， 從1估計發(fā)射器或反射器的相對延遲和距離，
從m估計發(fā)射器或反射器的多普勒和相對速度/運動學性質， 從s估計反射器的角位置。參考文獻[1]H. L. VAN TREES, "Detection, Estimation and ModulationTheory", Part I， II，III，Wiley and Sons,1971[2]J. DARRICAUT, "Physique et Theorie du radar", Volumes 1，2，3，Editions Sodipe,1993,[3]R. A. M0NZING0,T. W MILLER,"Introduction to adaptivearrays", John Wiley and sons,1981[4]ff. S. BURDIC,"Underwater Acoustic System Analysis",Prentice Hall Sig. Proc. Series 2nd edition, 19991.
權利要求
一種用于在包括N個傳感器的系統(tǒng)中確定涉及固定或移動的發(fā)射器、反射器或傳感器的傳播信道脈沖響應的一個或多個參數、探測和確定這些發(fā)射器和反射器的位置和運動學參數的方法，所述傳感器接收來自所述反射器上的反射或所述發(fā)射器的信號，所述方法至少包括以下步驟○利用空間分析和延遲-距離/多普勒-運動學分析的耦合確定模糊函數，○確定對應于已知信號s(kTe)和濾波器w(l，m)輸出之間的相關性的至少一個充分統(tǒng)計量所述已知信號對應于所發(fā)射信號的復包絡且經過了與關于所述發(fā)射器、接收器和反射器的位置和運動學參數和關于傳播的波動的假設相關聯(lián)的可能變換，其中l(wèi)是一組時域假設的符號表示，其中m是一組頻率假設的符號表示，通過將用于對(l，m)的統(tǒng)計值與探測閾值比較確定所述對(l，m)的值。FPA00001087378300011.tif
2.根據權利要求1所述的方法，在將1符號表示的時域參數簡化為延遲或距離的情況 下，且在將m符號表示的運動學參數簡化為多普勒或相對速度時，其特征在于至少包括以 下步驟對于每種延遲假設1和每種多普勒假設m 其中m是對變量Xm和..Gmi做索引的多普勒/運動學參數，-利用所發(fā)射的信號的K個已知樣本s (kTe)以及在頻率上偏移的K個觀測值<formula>formula see original document page 2</formula>，估計矢量<formula>formula see original document page 2</formula> (5)；-從頻率上偏移的所述K個觀測值<formula>formula see original document page 2</formula>估計 矩陣及K-I ；-計算所述充分統(tǒng)計量<formula>formula see original document page 2</formula>其中&是由所述已知信號的平均功率的最小二乘方估計定義的，-針對給定誤報警概率設置所述探測閾值，-將所述統(tǒng)計值 (1, m, K)與所述固定閾值比較，如果超過所述閾值，要求探測距離/ 多普勒小區(qū)(l，m)中的發(fā)射器或反射器；-從1估計所述發(fā)射器或反射器的相對延遲和距離， -從m估計所述發(fā)射器或反射器的相對多普勒和速度；-從由》=給出的源自所述發(fā)射器或反射器的信號天線上的指引矢量it 估計所述發(fā)射器或反射器的空間小區(qū)。
3.根據權利要求2所述的方法，其特征在于，在最復雜的情況下，3A)將時域參數引入諸如由1符號表示的延遲-位置的定義、引入諸如由m符號表示 的多普勒-運動學參數的頻率參數的定義，并引入每個假設1和m上的所述已知信號s的 定義，信號的變化和變換由所述發(fā)射器、接收器或反射器的運動引起，由傳播信道的波動引 起，這取決于延遲-位置假設、相對多普勒-速度假設、處理積分時間、信號的等價帶寬、當 在處理積分時間期間信號是脈沖式時的信號的等價時長及其周期、延遲隨這一時間的波 動、多普勒效應隨這一時間的波動，3B)根據由1符號表示的時域假設和由m符號表示的頻率假設使用從s的變換獲得的f曰號S1 mo
4.根據權利要求3所述的方法，其特征在于所述延遲的波動是由時域斜變建模到一階 的波動。
5.根據權利要求3所述的方法，其特征在于所述多普勒波動由多普勒斜變建模到一階。
6.根據權利要求2和3之一所述的方法，其特征在于對于信號受到高干擾的發(fā)射器或 反射器的探測而言，使用信號X’ m((k+l)Te)和仿樣S’而不是輸入信號xm((k+l)Te和原始 仿樣s，信號x’m((k+l)Te))是從信號1(&+1)1；)獲得的，仿樣S’是借助于涉及增強空-時 濾波器的操作從原始仿樣s或變換的仿樣S1,m獲得的，所述操作如下6A)對于每個干擾發(fā)射器通過使準則‘《仏，L) ^-.Y |對(/。+07>) - ginttSi nl I 最小化估計從所述發(fā)身寸器至丨J所述L O接收器的傳播濾波器的脈沖響應gint，仿樣Sint = (sint((l0+l)Te), sint((l0+l+l)Te),..., 、((1。+1+卜1)1^))1或5她 1,1 “ = (sintm,/'((Vl)Te), Sintnbl“ ((10+1+1) Te), ···，Sintffl, 廣((lfl+L-DTe))1是通過原始干擾信號的總先驗知識、或通過原始干擾信號的部分先驗 知識、或借助本領域技術人員已知的任何方式通過對所述原始干擾信號的解調/均等獲得 的，并經過可能的變換而成Sin^1 ；6B)對于固定的或運動學已知且根據它們的頻率和時間穩(wěn)定性選擇的一個或多個干擾 發(fā)射器而言，向所述仿樣Sint中引入不同的頻率假設，并采用所述仿樣中存在的重現，以便在大的積 分長度K上精確估計相對于所述干擾發(fā)射器的頻率和時間偏離，針對用于接收和取樣用于 實施本發(fā)明的信號的系統(tǒng)的混頻器和時鐘的本地振蕩器估計相對于所述干擾發(fā)射器的漂 移。
7.根據權利要求6所述的方法，其特征在于，在步驟6B期間 更新與所述干擾發(fā)射器相關聯(lián)的傳播濾波器gint (步驟6A)，以及 6C/對于每個干擾發(fā)射器，由與其關聯(lián)的傳播濾波器gint執(zhí)行原始Sint或變換的Sintnbl干擾信號的卷積，從總的所 接收信號減去這一卷積操作獲得的信號， 6D/補償頻率和時間漂移是-在前面操作之后從xm((k+l)Te)獲得的信號上執(zhí)行，以便產生信號χ’ m((k+l)Te),s 然后保持不變(S’ = S)-或者在該方法中使用的仿樣信號上執(zhí)行，在權利要求2中限定的s于是可以變?yōu)閟’、 在權利要求3中定義的Snul可以變?yōu)镾ma’ ；6E/以χ’ m和S’實施根據權利要求2所述的步驟 6F/以χ’ m和Snul’實施根據權利要求3所述的步驟。
8.根據權利要求2、3和6之一所述的方法，其特征在于使用觀測矢量 Xm((k+l)Te) ^ [xm{{k+DTe)T, xm((k+l)Te)hT ,所述觀測矢量由觀測矢量及其共軛復數構成，且其特征在于從仏)而不是xm((k+l)Te)實施權利要求2和3所述的步驟。
9.根據權利要求8所述的方法，其特征在于還有僅總噪聲參考可用，且其特征在于利用經由僅噪聲參考樣本獲得的擴展僅總噪聲矩陣的估計》δ7> .(/7^，取 代擴展的觀測相關矩陣的估計來實施權利要求2和3所述的步驟，其中b4 [bm((k+l)re)T, bm((k+l)Te)^]T.。
10.根據權利要求2、3、4、5和9之一所述的方法，其中在去旋之后所 發(fā)射信號是準直線性的，所述干擾不是圓形的，或者所述傳播信道在時間 上是擴散的，其特征在于，使用的空時矢量對于準直線性調制而言是去旋 的，或對于具有L個時域獲取的直線性調制是非去旋的，如果L是奇數，XdmM(^Te)Hx dm((kHL-l)/2)Te)T,..., χ dm((k- (L-l)/2)Te)T]T,如 果 L 是偶數，χ dm,st{^+l)Te)-[^dm{{k +"2)Γβ)Τ，…，3M(k -L/2)Te)T]T , 其中艾 dm((k+I)Te) ^ [xdm_Te)\ xdm((k+l)Te)^T,xdm((k+l)Te)^xd((k+l)Te) Qxp[-}2nm(k+l)/K], xd(tTΓ t/Tx(0，了 是符號時長。
11.根據權利要求2、3、4、8和9所述的方法，對于所述方法而言，所述發(fā)射器或反射器 的指引矢量是已知的，且其包括以下步驟 從所發(fā)射信號的K個已知樣本s(kTe)和K個頻率偏移觀測值Xm((k+l)Te)= x((k+lTe)exp[-j23im(k+l)/K]估計矢量？ XwS(ZT1e), 0 < k< K- 從 K 個頻率偏移觀測值 xm((k+l)Te) = x((k+l)Te)exp[-j2Jim(k+l)/K]估計矩陣RMlTe) 計算所述充分統(tǒng)計量'C(IfMiK) ^ (1 /%) 的)卞 Iw(Zre)-1 &其中b是由⑶定義的， 計算所述充分統(tǒng)計量=(!/%)丨 f ‘肌) (叫丨2 (s卞 RxmQTeT1 s) (1 — ￡(/, m，K)) 針對給定誤報警概率設置所述探測閾值， 將所述統(tǒng)計值:K/, W，5, TQ與設置的探測閾值比較，如果超過該閾值，要求探測距離 /多普勒小區(qū)(1，m)和空間小區(qū)S中的反射器， 從1估計所述發(fā)射器或反射器的相對延遲和距離， 從m估計所述發(fā)射器或反射器的相對多普勒和速度/運動學性質， 從s估計所述反射器的角位置。
12.根據權利要求2、3、4、7、8和9之一所述的方法，其特征在于有僅總噪聲參考可用， 且其特征在于利用經由僅噪聲參考樣本獲得的僅總噪聲矩陣的估計力取代觀測 相關矩陣的估計i ； 來實施權利要求2、3或7所述的步驟。
13.根據權利要求2、3、4、7、8和9之一所述的方法，其中對于直線發(fā)射信號和潛在非圓 干擾而言，所述反射器的指引矢量是已知的，其特征在于其至少包括以下步驟計算所述充分統(tǒng)計量<formula>formula see original document page 5</formula>構造矩陣<formula>formula see original document page 5</formula>構造統(tǒng)計量<formula>formula see original document page 5</formula>針對給定誤報警概率設置所述探測閾值，將所述統(tǒng)計值fnc(7，m, s, k)與所述閾值比較，如果超過所述閾值，則探測空間 小區(qū)S和距離/多普勒小區(qū)(l，m)中的發(fā)射器或反射器， 從1估計所述發(fā)射器或反射器的相對延遲和距離， 從m估計所述發(fā)射器或反射器的相對多普勒和速度/運動學性質， 從s估計所述發(fā)射器或反射器的角位置。
14.根據權利要求2、3、4、8、9、10和11之一所述的方法，其特征在于還有僅總噪 聲參考可用，且其特征在于利用經由僅噪聲參考樣本獲得的擴展僅總噪聲矩陣的估計Rbtm(lTe)，取代擴展的觀測相關矩陣的估計Rxm(lTe)來實施權利要求13所述的步驟，其中bm((k+l)Te)△[bmiik+^Te)1 ,bm{{k+riTe^]Tο
15.根據前述權利要求之一所述的方法，其特征在于分別用未多普勒化(m= 0)且未時 域偏移(1 = 0)的觀測矢量x(kTe)，7 U7^)，Xst(kTe)或文d,MTe)，替換多普勒化和時 域偏移的觀測矢量 Xffl ((k+1) Te)，X m{ U&i) Te)，Xffl, st ((k+1) Te)或文 st ((I&l) Te)，其中 參考信號S (kTe)是多普勒化的且時域偏移的，即被Sm ((k+1) Te) = S (kTe) exp [j2 π m(k+1) / K]替換。
16.根據前述權利要求之一所述的方法，其特征在于，在源自所述發(fā)射器、所述反射器 或所述干擾的信號對于所用天線陣列而言是寬帶時，一旦比例Byftl或π DBe/c超過幾個百 分比就發(fā)生這種情況，其中D是所述陣列的孔徑，Be是信號的等價帶寬(視為包括在接收頻 帶中)，fo是信號的載波，c是傳播速度，實施以下變體實施例之一_將所述接收帶構造為幾個子帶，對于陣列而言是窄的；一次一個子帶地實施作為前 述權利要求主題的處理操作， -實施接收空間-時間結構， 在所述信號對于傳感器陣列不再是窄帶的情況下(n DByc超過幾個百分比) 利用具有L個時域獲取的空時矢量替換前述空間矢量xm((k+l) Te)如果L為奇數，xm,st((k+l)Te)^ [xm((k+(L-l)/2)Te)T,...,xm((k- (L-l)/2)Te)T]T，如果 L 為偶數，xm,st((k+l)Te) ^ [xdm((k+ L/2)Te)T，...，xdm((k-L/2)Te)T]T; 在與載波相比信號也是寬帶的情況下(Byftl超過幾個百分比)以及所述發(fā)射器或反射器的運動學性質要求考慮在仿樣信號S上誘發(fā)形變的情況下向前述處理操作中不僅引 入延遲和多普勒變量，而且引入諸如積分時間K上的延遲壓縮和多普勒壓縮的額外1階變 量，諸如前述變量的時域導數的2階變量等，通過1符號表示的時域變量(延遲_位置等) 和m符號表示的頻率變量(多普勒-運動學性質等)寫入對應于原始信號s的變換的信號 Sl,m，并使用信號S1,m來實施如前述權利要求中主張的方法的步驟。
17.根據前述權利要求之一所述的方法，其特征在于其包括基于延遲-距離/多普 勒_運動學或空間/延遲_距離/多普勒_運動學假設的常規(guī)測角步驟(高分辨率等)或 參考測角步驟，包括通過以下方式實施常規(guī)或高分辨率測角技術 從與能量距離/多普勒(1，m)相關聯(lián)的觀測值<formula>formula see original document page 6</formula>實施，或 從作為K約數的時段上的觀測值<formula>formula see original document page 6</formula>或與所發(fā)射信號s(kTe)的相關結果實施，或 在指引矢量s已知的情況下，在將觀測值的可能空間濾波限制到空間小區(qū)S之后實施。
18.—種在包括N個傳感器的系統(tǒng)中測量傳播，探測和/或定位固定或移動的發(fā)射器或 反射器的裝置，所述傳感器從所述發(fā)射器或反射器接收信號，其特征在于其包括適于執(zhí)行 根據權利要求1到17之一所述的方法步驟的至少一個模塊。
全文摘要
用于在包括N個傳感器的系統(tǒng)中測量涉及固定或移動發(fā)射器和反射器的傳播信道的脈沖響應特定參數的方法，以便探測和確定關于這些發(fā)射器和反射器的位置和運動學性質的參數，或自動定位實施本發(fā)明的接收系統(tǒng)，傳感器從所述發(fā)射器或從所述反射器的反射接收信號，所述方法至少包括以下步驟確定耦合空間分析和延遲-距離/多普勒-運動學分析的模糊函數，確定與對應于所發(fā)射信號的復包絡的已知信號s(kTe)和濾波器w(l，m)的輸出之間的相關相對應的至少一個充分統(tǒng)計量C(l，m，K)，其中l(wèi)對應于時域假設(延遲-距離，時域斜變等)，m對應于頻率假設(多普勒-運動學、多普勒斜變等)，通過將對(l，m)的統(tǒng)計值與閾值比較確定對(l，m)的值。
文檔編號G01S7/292GK101821646SQ200880110838
公開日2010年9月1日 申請日期2008年8月8日 優(yōu)先權日2007年8月8日
發(fā)明者F·德拉沃, F·皮彭, P·舍瓦利耶 申請人:塔萊斯公司