引 言
反應(yīng)速度更快、精度更高、目標(biāo)的自動錄取、使操作員能夠同時處理多批目標(biāo),是現(xiàn)代雷達所追求的技術(shù)指標(biāo)。雷達視頻信號的數(shù)字采集和檢測是達成這樣目標(biāo)的前提條件。眾所周知,雷達對目標(biāo)位置的測量,主要是依據(jù)目標(biāo)回波相對于發(fā)射電波的延遲時間以及雷達天線的方向來決定的。電波的傳播速度極快,傳播l海里的時間大約只有12.35μs,雷達發(fā)射電波的間隔一般在ms數(shù)量級。在此間隔時間內(nèi),理論上最多會有成百上千個目標(biāo)的回波。要在如此短的時間之內(nèi),對這么多的目標(biāo)回波信號進行數(shù)字采集和處理,快速的采集器件及高速計算機處理器必不可少。DSP(Digital Signal Processor)芯片,即數(shù)字信號處理器,就是這樣一種特別適用于進行這種實時信號處理的微處理器。本文以某雷達的技術(shù)參數(shù)為參考,用TMS320VC5402作為信號處理器,以AD9223為A/D變換器,采用雙端口存儲芯片IDT7203作數(shù)據(jù)暫存存儲器,介紹一種實現(xiàn)雷達回波信號的數(shù)字采集和處理方法。
1 主要器件特性
TMS320VC54x是TI公司開發(fā)的低功耗、高性能數(shù)字信號處理芯片,主要應(yīng)用于無線通信系統(tǒng)及雷達信號處
理系統(tǒng)等。本文所采用的TMS320VC5402芯片是此系列的一個典型產(chǎn)品。該芯片的主要特點有:速度快,指令周期10 ns,運算能力為100 MIPS;強大的尋址能力,最大可尋址外部存儲器1 M×16位,內(nèi)置16 K×16位的RAM,4K×16位的ROM;40位的算術(shù)邏輯運算單元(ALU),包括2個獨立的40位累加器和1個40位的桶形移位寄存器;1個17位×17位的硬件乘法器和1個40位的專用加法器,乘法器/加法器單元可在一個流水線狀態(tài)周期內(nèi)完成一次乘法累加運算(MA)。
AD9223是一種采用高速CMOS處理和新穎的4級流水線結(jié)構(gòu)的高性能、低噪聲、單電源、12位模/數(shù)變換器。適用于航海、通信、雷達、醫(yī)療等領(lǐng)域的數(shù)據(jù)采集系統(tǒng)。其主要特點有:高速度、高分辨力,采樣時鐘頻率可達3 Msps,分辨力為12位;單時鐘輸入,采用單時鐘輸入控制內(nèi)部所有的轉(zhuǎn)換周期;模擬信號輸入靈活,真正的差分輸入結(jié)構(gòu)允許模擬信號以單端或差分兩種分式輸入;靈活的參考電壓,適應(yīng)不同應(yīng)用的精度及溫度漂移的要求,參考電壓可選擇內(nèi)部或外部(1 V或2.5 V);采用4級流水線結(jié)構(gòu)及寬帶取樣保持放大器(SHA),使器件在每個時鐘周期都可捕獲輸入的采樣信號,數(shù)據(jù)的輸出延時為3個時鐘周期。
IDT7203是一種雙端口先進先出(FIFO)9位存儲緩沖器,存儲容量為2 KB。輸入和輸出具有各自的地址指針,每一個讀或?qū)懖僮?,對?yīng)的輸出及輸入地址指針自動加1。復(fù)位信號可將兩個地址指針全部清除。對外提供數(shù)據(jù)區(qū)空(empty)、半滿(half-full)及滿(full)信號,用來指示器件的狀態(tài)。最快讀寫速度為12 ns。
2 二級門限判決的工作原理
雷達的視頻信號中始終包含有目標(biāo)信號和噪聲信號,即x(t)=s(t)+n(t)。雷達信號處理的根本目的就是,從接收到的視頻信號x(t)中提取出有用的目標(biāo)信號s(t),濾除無用的雜波或干擾信號n(t)。在進行雷達信號處理時,預(yù)先并不知道目標(biāo)在一定條件下出現(xiàn)的概率,也很難確定一次漏警所造成的損失,故常采用奈曼一皮爾遜準則,即在允許一定虛警概率的條件下,使漏警概率達到最小。二級門限判決方法,是一種有效而實用的雷達信號處理方法。
2.1 二級門限判決原理
由最佳檢測理論知道,對于非相參高頻脈沖串,采用檢波后的積累(視頻積累)來改善雷達檢測性能是一種最佳方案,其原理如圖1所示。這種方案若采用模擬電路處理,比較復(fù)雜,實現(xiàn)起來也相當(dāng)困難;若采用數(shù)字化處理,直觀的辦法是將雷達接收機的視頻信號按距離(時間)和幅度進行量化處理變?yōu)閿?shù)字信號,然后存儲N個重復(fù)周期。對每一個重復(fù)周期按距離單元對每一個單元的信號幅度值設(shè)定一個門限值(r0),稱為第一門限。對每一個重復(fù)周期所量化的所有距離單元的視頻信號值x(ti),與其相對應(yīng)的門限值roi進行幅度值的比較。如果視頻信號值超出它的門限值,認為是“可能目標(biāo)”信號;否則,認為無信號。這個過程稱為第一級門限判決。對應(yīng)于每一個重復(fù)周期內(nèi)每一個距離單元所建立的第一門限值(roi)是各不相同的,而且,每一個重復(fù)周期所建立的也各不相同。這里將每一個重復(fù)周期所建立的第一門限值的集合稱為雜波均值估值表或叫雜波圖。對于“可能目標(biāo)”信號,并不能肯定它一定是目標(biāo)信號。因為有時突發(fā)的雜亂干擾信號或隨機噪聲信號,各次掃描的取樣是不相關(guān)的。它們可能會偶然一二次超過第一門限值,但連續(xù)多次超過第一門限值的可能性很?。欢鴮τ谀繕?biāo)來說,相鄰各次掃描均應(yīng)有回波信號,在N個取樣中連續(xù)超過第一門限值的概率較大。因此,需要按距離單元將超過各自第一門限值的量化脈沖進行計數(shù),如果在N個重復(fù)周期中有K個以上的量化脈沖超過各自的第一門限值,則判決為有信號。這個過程稱為第二級門限判決(K/N判決),而K/N稱為第二門限值,工作原王單如圖2所示
2.2 雜波均值估值表的建立
雜波均值估值表(雜波圖),是根據(jù)每個重復(fù)周期雷達視頻信號的量化值所建立起來的每個距離單元的雜波信號的平均值表。建立的方法為鄰近單元平均估值法。具體是,設(shè)某個距離單元ti上的視頻信號為x(ti),以ti單元為中心,前后共取n個參考距離單元,對這n個參考單元的視頻信號數(shù)值求平均值
作為這個距離單元的雜波平均估值。也叫作這個距離單元的第一門限值roi值的大小由雷達探測區(qū)域的雜波環(huán)境所決定。n值愈小,由于參考單元太少,會引起雜波平均估值的起伏變化愈大,將引起虛警概率的增加;而n值愈大時,雜波的平均估值起伏愈小,即第一門限值愈平穩(wěn),可得到恒虛警效果,但會使得資源開銷增大,增加實現(xiàn)的難度。對于氣象和海浪雜波來說,它們通常是連接成片的區(qū)域,比較適合采用鄰近單元平均估值法建立雜波均值表。這里n取16。由于以52.5 m為一個距離單元,相當(dāng)于對檢測點前后各420 m
的距離范圍取雜波的平均估值。
2.3 第二級門限值K/N的計算
為保證符合奈曼一皮爾遜準則,N值的選擇應(yīng)取決于天線波束掃過目標(biāo)的時間內(nèi)目標(biāo)的回波脈沖數(shù),是由雷達的工作參數(shù)決定的。對于某船用雷達來說,其部分工作參數(shù)為:重復(fù)頻率1 200 Hz,觸發(fā)脈沖寬度O.7μs,天線轉(zhuǎn)速20 r/min,天線波束寬度O.7°,作用距離40海里。N的取值為N=O.7×1 200/(20×6)=7。根據(jù)經(jīng)驗,K的取值為1.5√N。這里K=4,即K/N=4/7。
3 硬件電路設(shè)計
硬件電路工作原理如圖3所示。根據(jù)雷達的主要技術(shù)參數(shù),得知雷達的距離分辨力為150×0.7=105 m,在不降低原雷達性能的前提下,這里選擇2.857 142 857 MHz的采樣頻率(相當(dāng)于以52.5 m作為一個采樣距離單元)。存儲器以同樣的速率采集數(shù)據(jù)(大約每52.5 m采集一個雷達視頻信號的數(shù)值)。這樣的話,40海里的最大測距量程,對應(yīng)于雷達的每一次發(fā)射,電路應(yīng)當(dāng)采集1 412個距離單元的視頻信號量化值(這里取1500個距離單元),因此,雙端口RAM容量選用2 KB即可。由于訪問速度要快,選用了IDT7203。對于信號的量化位數(shù),即A/D轉(zhuǎn)換位數(shù)的選擇,考慮到量化噪聲的影響,所取的位數(shù)愈多影響愈小。為照顧DSP的處理能力以及A/D轉(zhuǎn)換器件的轉(zhuǎn)換速度,采用12位。
A/D轉(zhuǎn)換器AD9223的設(shè)置:采用片內(nèi)2.5 V作為參考電壓;VINB接參考電壓值,這樣VINA的輸入值最大可到5 V,最小為O V;運算放大器U3提供的視頻信號應(yīng)當(dāng)符合這個要求。
雷達的觸發(fā)脈沖作為系統(tǒng)工作的同步信號,從觸發(fā)脈沖的下降沿開始,為A/D轉(zhuǎn)換器提供2.857 142 857 MHz的采樣同步脈沖信號,A/D轉(zhuǎn)換器產(chǎn)生的12位精度的視頻信號數(shù)值以同樣的速率被打入到二片雙端口RAM(IDT7203)中。當(dāng)打入的數(shù)據(jù)達到1 504個時(由于A/D轉(zhuǎn)換器的數(shù)據(jù)輸出對視頻輸入信號有3個時鐘周期的延遲,因此,這里采用采集1 504個距離單元的數(shù)據(jù)),對DSP產(chǎn)生一個中斷信號。這些信號是由大規(guī)模可編程邏輯器件(Lattice M4A5-128/64)根據(jù)雷達觸發(fā)脈沖、晶體振蕩器以及DSP的I/O操作所產(chǎn)生的地址信號和讀寫信號共同作用而產(chǎn)生的。它們的時序關(guān)系如圖4所示。
DSP對中斷信號響應(yīng)后,執(zhí)行數(shù)據(jù)接收的中斷服務(wù)程序。該程序以極快的速度把雙端口RAM中的數(shù)據(jù)提取到自己的內(nèi)部存儲器中,這時DSP再運行濾波和目標(biāo)提取程序,最后將完全處理好的數(shù)據(jù)通過串行接口傳送出去。
雙端口RAM的讀信號RAMRD在執(zhí)行讀OFF00H地址的指令時,由DSP產(chǎn)生的地址信號(OFF00H)、WR信號以及10STRB信號的共同作用而產(chǎn)生。
4 軟件設(shè)計
軟件的主要任務(wù)是處理已經(jīng)被數(shù)字化了的視頻信號。為了實現(xiàn)快速處理,程序采用匯編語言編寫,保證在一個重復(fù)周期內(nèi)將數(shù)據(jù)處理完畢。程序運行的區(qū)域安排在TMS320VC5402片內(nèi)RAM的0000H~03FFH的1 KB內(nèi)存中。數(shù)據(jù)輸入?yún)^(qū)在0400H~09FFH,占1.5 KB內(nèi)存。雜波均值估值表(雜波圖)位于OA00H~OFFFH,占1.5 KB內(nèi)存。輸出數(shù)據(jù)區(qū)(處理結(jié)果)位于1000H~15FFH,占1.5 KB內(nèi)存。從1600H開始至3FFFH的10.5 KB存儲區(qū),被平均分為7個數(shù)據(jù)區(qū),作為7個重復(fù)周期的數(shù)據(jù)積累區(qū)。對于存儲雷達視頻信號量化數(shù)據(jù)的每一個存儲區(qū),其每一個地址單元的地址代表一個距離單元,其內(nèi)存儲的數(shù)據(jù)代表該距離單元處視頻信號的幅度量化值。程序響應(yīng)中斷后,工作過程為:數(shù)據(jù)的采集→建立雜波均值估值表(雜波圖)→第一門限判決→第二門限判決→數(shù)據(jù)輸出。方位信號的采集以及處理以后數(shù)據(jù)的輸出,此處不作討論。
4.1數(shù)據(jù)的采集
由于TMS320VC5402運行的速度極快,單指令周期為10 ns,而雙端口RAM的最快訪問速度為12 ns。為了保證數(shù)據(jù)的可靠讀取,把I/0操作設(shè)置為延遲2個執(zhí)行周期。DSP對于外部中斷0(INTO)響應(yīng)后進人中斷服務(wù),把雙端口RAM中的數(shù)據(jù)讀到DSP的內(nèi)置存儲器以0408H起始的1500個單元中。程序初始化時,將0400H起始的8個單元設(shè)置為0。這完全是為了計算雜波均值估值表(雜波圖)時的編程方便,其程序如下:
GETDATA:STM #0408H,AR6 ;初始地址
RPT#(1500—1) ;重復(fù)次數(shù)
PORTR OFF00H,*AR6+ ;讀數(shù)據(jù)
RETE &nbs
p; ;中斷返回
4.2 建立雜波均值估值表
根據(jù)3.2雜波均值估值表的建立方法,對每一個距離單元,前取8個,后取7個,共16個作參考單元,計算平均值,作該單元的雜波均值估值??偣灿嬎?500個距離單元的雜波均值估值,放置于OA00H~OFFFH的存儲區(qū)域。具體實現(xiàn)時,須先計算出第一個單元的雜波平均值估值,其后的逐個計算。程序代碼如下:
STM #O4ooH,AR5 ;讀入數(shù)據(jù)的首址
STM #OA00H,AR7 ;均值估值首址
RPTZ A,(16—1) ;循環(huán)16次
ADD*AR5+,A ;求和
STL A,一4,*AR7+ ;除以16,放入估值表
STM #0400H,AR6 ;計算余下的1 499個
STM#(1500—2),BRC
RPTB ZBEND一1
SUB*AR6+,A
ADD*AR5+,A
STL A,一4,*AR7+
ZBEND:RET
4.3 第一級門限判決
將讀入的雷達視頻信號值以距離單元為單位,逐個與雜波均值估值表中對應(yīng)距離單元的數(shù)值作比較:如果某單元雷達視頻信號的數(shù)值較其對應(yīng)的雜波均值估值大,認為此單元為“可能目標(biāo)”的回波信號,將它放置到數(shù)據(jù)積累區(qū)中(1600H起始)相應(yīng)的單元,同時也放置到輸出數(shù)據(jù)塊中(1000H起始);否則,這兩個單元都被清零。程序中用AR3作數(shù)據(jù)積累區(qū)的地址指針,其初值為1600H,由初始化程序設(shè)置。程序代碼如下:
STM #0408H,AR7 ;輸入數(shù)據(jù)區(qū)首址
STM #oA00H,AR6 ;均值估值表首址
STM #1000H,AR5 ;輸出數(shù)據(jù)區(qū)首址
STM#(1500—1),BRC: ;距離單元總數(shù)
RPTB PJlE一1 ;重復(fù)
LD*AR6+,A ;取雜波均值估值
STM*AR7,T ;取輸入的數(shù)據(jù)
ST #0,*AR3 ;積累區(qū)先清零
ST #O,*AR5 ;輸出區(qū)先清零
SUB*AR7+,A ;比較
STRCD*AR3+,ALT ;數(shù)據(jù)大,分別放置
STRCD*AR5+,ALT ;到積累區(qū)和輸出區(qū)
PJlE:LD AR3,B ;調(diào)整積累區(qū)指針
ADD #100H,B
AND #OFF00H,B
SUB #4000H,B
BC PJlEND,BLT
SUB #2A00H,B
PJIEND:ADD #4000H
,B
STML B,AR3 ;調(diào)整積累區(qū)指針結(jié)束
RET
4.4 第二級門限判決
第二門限判決就是在第一門限判決的基礎(chǔ)上,濾除那些偶然超過第一門限,但是又是在多次掃描中,孤立的、互不相關(guān)的雜波噪聲信號,就是對那些“可能目標(biāo)”的信號作進一步的處理,進一步降低檢測目標(biāo)的虛警概率。根據(jù)2.3節(jié)的討論,第二門限值(K/N)的取值為4/7。由于在第一門限判決處理時已經(jīng)建立了7次掃描的累積數(shù)據(jù)和初步的輸出數(shù)據(jù),它們分別位于數(shù)據(jù)累積區(qū)和數(shù)據(jù)輸出區(qū)。因此,在程序的具體實現(xiàn)時,就是根據(jù)最新的積累數(shù)據(jù),對量程范圍內(nèi)的每個距離單元逐個進行K/N判決,符合K/N判決準則的,輸出數(shù)據(jù)區(qū)的數(shù)值不變,而不符合準則的將該距離單元所在數(shù)據(jù)輸出區(qū)的數(shù)值清零。由于程序代碼較長,用圖5所示流程表示。
結(jié) 語
本文介紹的基于DSP的雷達視頻信號的數(shù)字采集和處理系統(tǒng)的設(shè)計,作為一種雷達視頻信號處理的方法,適用于多種雷達的信號處理,對于不同型號的雷達,由于其用途的不同,其技術(shù)參數(shù)也不相同。只要選擇合適的工作參數(shù),此設(shè)計對其他種類雷達的視頻信號仍然適用。