摘 要： 針對甚高頻（VHF）云閃時差法雷電定位同源脈沖匹配難點,，提出了用FPGA實現(xiàn)相關(guān)性判斷的同源脈沖算法設(shè)計思路。為提升系統(tǒng)性能,、適應(yīng)云閃甚高頻信號的處理需求,，采用了改進的算法并用實測數(shù)據(jù)在MATLAB、FPGA上對算法進行驗證,，設(shè)計完成了云閃波形同源脈沖匹配算法在FPGA上的實現(xiàn),。系統(tǒng)通過了Modelsim的仿真驗證并在DE2平臺上完成了硬件測試。設(shè)計共消耗了3 499個邏輯單元,，最高處理速度可達98.07 MHz,，滿足了實際應(yīng)用的要求。
關(guān)鍵詞： 云閃,；時差法,；同源脈沖；FPGA
雷電是在雷暴天氣條件下,，發(fā)生于大氣中的一種長距離放電現(xiàn)象,，具有大電流、高電壓,、強電磁輻射等特征[1],。雷電主要包括云對地面的放電（地閃）和云與云之間的放電（云閃）[2]，地閃往往對地面民用工業(yè),、軍事設(shè)施有嚴(yán)重的破壞性[3],；云閃也嚴(yán)重影響民航飛機、航天飛行器等的正常飛行,。觀測表明,，云閃先于地閃5 min~15 min發(fā)生[4]，觀測云閃對地閃有預(yù)警的作用,，因此云閃監(jiān)測及定位研究具有重要意義,。
目前對于地閃的探測已經(jīng)有了比較完善的定位系統(tǒng)，如ADTD雷擊檢測儀[5],，但僅僅能檢測地閃,。云閃通?；谏醺哳l（VHF）頻段，已有的甚高頻云閃探測主要采用干涉法和時差法（TOA）建立雷電探測系統(tǒng),，時差法主要難點是同源脈沖匹配技術(shù),。OETZEL等利用VHF技術(shù)測量閃電輻射到達分離天線的時差確定閃電源位置的成果[6]，采用短基線避免判斷兩個脈沖是否屬于同一個輻射脈沖的問題,，但是短基線帶來的問題是信號到達兩天線的時差太小,，導(dǎo)致測向精度不高?；诖?，本文提出一種同源脈沖匹配算法——相關(guān)性判別法,該方法可以應(yīng)用于長基線的VHF云閃脈沖同源判別。
從數(shù)學(xué)上可以證明分子的模小于分母,，也即相關(guān)數(shù)Pxy的模不會大于1,。當(dāng)相關(guān)系數(shù)為0時相似度最差，即不相關(guān),。當(dāng)相關(guān)系數(shù)為1時,，誤差能量為0，說明兩信號相似度很好,，是線形相關(guān)的,。1980年，RUSTAN等發(fā)表了在KENNEDY航天中心觀測閃電放電過程的研究報告[7],，曾借助互相關(guān)分析,，通過圖形識別技術(shù)識別共源脈沖。據(jù)此,，可以將相關(guān)性判別應(yīng)用在VHF云閃波形相似度的判別中,。
2 算法的MATLAB驗證
2.1 算法驗證平臺
本文所采用的驗證平臺由天線、信道,、PC端數(shù)據(jù)采集卡以及DE2開發(fā)板四部分構(gòu)成,。如圖1所示。
天線單元由SG7200雙段車載天線在全頻段對云閃信號進行接收,，經(jīng)帶通濾波,、放大。為了前端觀測及處理方便,，選取美國GAGE公司產(chǎn)品CS21G8數(shù)據(jù)采集卡對云閃信號進行數(shù)據(jù)采集,可以實時觀測波形變化,。由于云閃定位系統(tǒng)采用FPGA實現(xiàn)，為了在FPGA中驗證算法的可行性,，將采集到的數(shù)據(jù)經(jīng)PC交由DE2開發(fā)板處理,。為了更真實地模擬實際系統(tǒng)，將需要進行相關(guān)性計算的兩組數(shù)據(jù)存在DE2上的SRAM中,。限于DE2上只有一片SRAM,，圖中將一片SRAM地址分為兩塊,，分別存放兩組擴大后的24位定點數(shù)據(jù)x(t)、y(t),，兩組數(shù)據(jù)在FPGA中完成同源匹配的計算,，若同源則輸出高電平給下一級定位系統(tǒng)。
2.2 MATLAB算法驗證
基于以上平臺,，2010年3月～5月,，項目組在北京歷經(jīng)兩個月的時間采集到大量雷電波形數(shù)據(jù)，選取其中編號為CH11707和CH11708的兩組數(shù)據(jù)在MATLAB中進行相關(guān)性算法驗證,，將兩組數(shù)據(jù)在MATLAB中作圖,，如圖2所示,。
圖2中尖峰脈沖為云閃輻射脈沖,，寬的部分為噪聲，由于信號到達兩天線的距離不同,，故有一定的時間差,，幅度衰減也不同。相關(guān)系數(shù)計算結(jié)果為-0.601 15,。對采集到的大量數(shù)據(jù)進行相關(guān)性計算,，對比結(jié)果可以發(fā)現(xiàn)，只有同一個云閃發(fā)出的同源脈沖相關(guān)系數(shù)結(jié)果是最大的,，并且與非同源計算出的相關(guān)系數(shù)差別很大,，如編號為CH11708和CH1203相關(guān)系數(shù)僅為0.001 974 7。
3 算法的FPGA實現(xiàn)
3.1 算法流程
基于以上理論分析,，可列出算法流程如圖3所示,。
(1)浮點定點化。采集卡采集到的數(shù)據(jù)為小數(shù)格式,，而FPGA中浮點數(shù)運算的電路比定點數(shù)復(fù)雜,，占用資源多，速度慢,?？紤]到算法與數(shù)據(jù)具體的幅值無關(guān)而與兩組波形的關(guān)系有關(guān)，因此,，可以將數(shù)據(jù)擴大適當(dāng)?shù)谋稊?shù)后截取為整數(shù)再輸入FPGA進行計算,。但是隨著處理數(shù)據(jù)位寬的增加，會導(dǎo)致硬件資源占用率提高,，降低運算效率,。綜合考慮，選擇擴大10 000倍,，以達到速度和精度上的折中,。
(2)定點浮點化,。由于最終計算結(jié)果是一個介于0～1之間的小數(shù)，因此除法計算前要將除數(shù)與被除數(shù)由定點數(shù)轉(zhuǎn)化成浮點數(shù),，再調(diào)用浮點數(shù)除法的IP核,，得到32位浮點格式的相關(guān)系數(shù)。
3.2 算法改進
按照上述設(shè)計思路,，在功能上實現(xiàn)了對云閃同源脈沖的判別,，但是算法最大時鐘頻率為52.70 MHz，由于云閃信號處于甚高頻頻段,，對系統(tǒng)處理速度上有較高要求,。為了更好地提升系統(tǒng)性能，分析發(fā)現(xiàn)限制時鐘頻率的最壞路徑位于開方模塊,，因此對開方模塊做如下改進：
(1)開方模塊：由于QUARTUS自帶的開方IP核將系統(tǒng)最大時鐘頻率限制在52.70 MHz,，為了改進算法的效率，采用循環(huán)冗余開方算法代替開方IP核,，循環(huán)冗余算法涉及到的運算有加法,、減法、移位,，加減法和移位較容易在硬件上實現(xiàn),，并且占用資源少，故可以提高算法的運算效率,。更改開方模塊后系統(tǒng)時鐘頻率可以提升至74.47 MHz,。此時系統(tǒng)時鐘頻率雖然有所提升，但還是被定點轉(zhuǎn)浮點模塊所限制,，故對定點轉(zhuǎn)浮點模塊進行改進,。
(2)定點轉(zhuǎn)浮點模塊：通常定點轉(zhuǎn)浮點的算法是通過查詢前導(dǎo)0的方式實現(xiàn)，由于定點轉(zhuǎn)浮點查詢前導(dǎo)0時采用了一系列的if,、else邏輯,，導(dǎo)致電路的關(guān)鍵路徑很長。因此考慮將if,、else邏輯替換成“與或”邏輯,，以此來縮短關(guān)鍵路徑、簡化電路,、達到速度提升的目的,。經(jīng)過更改定點轉(zhuǎn)浮點模塊，最大時鐘頻率可以提升至98.07 MHz,。
(3)優(yōu)化前后性能對比：改進算法后,，系統(tǒng)占用資源與之前對比結(jié)果如圖4所示。
由圖4可以看出,，經(jīng)過算法改進后,，最大時鐘頻率可以由原來的52.70 MHz提升至98.07 MHz,，并且邏輯單元占用數(shù)量上也明顯減少。
3.3 FPGA電路驗證
經(jīng)過對代碼和算法的優(yōu)化,，選取編號為CH11707和CH11708的兩組同源云閃輻射脈沖進行仿真,，具體仿真結(jié)果如圖5所示。
由圖5可以看到,，modelsim仿真的小數(shù)結(jié)果為：
Pxy=-0.601 144 9
MATLAB計算出的相關(guān)系數(shù)為-0.601 15,，modelsim仿真結(jié)果與MATLAB計算一致。
仿真成功后,，加入SRAM讀寫模塊,，在硬件上進行測試，通過嵌入式邏輯分析儀signaltap觀察實際的計算結(jié)果如圖6所示,。
由圖6可以看到,，實際硬件上計算的相關(guān)系數(shù)結(jié)果與modelsim仿真結(jié)果完全一致。
本文針對VHF云閃同源脈沖匹配的難題,，提出了相關(guān)性判別的匹配算法,，針對該思想完成了設(shè)計驗證以及FPGA硬件實現(xiàn),。在此基礎(chǔ)上,，利用電路和算法的優(yōu)化提高系統(tǒng)的性能，最大時鐘頻率由優(yōu)化前的52.70 MHz提升至98.07 MHz,，資源占用也有所減少,。所做的相關(guān)性算法的FPGA硬件實現(xiàn)可應(yīng)用于更廣泛的硬件實現(xiàn)相似性判斷，便于方便,、快速地識別信號間的相似性,。
參考文獻
[1] 佘會蓮，陳德生,，謝君.VHF閃電定位技術(shù)評述[J].氣象與環(huán)境科學(xué),，2008,，31(2).
[2] 張萍,，趙文光.基于云閃的雷電定位系統(tǒng)的誤差分析[J]. 測繪通報,，2009(1).
[3] 王懷斌,，郄秀書,，張義軍,，等.云閃和地閃的波形采集,、數(shù)據(jù)處理及其初步應(yīng)用[J].高原氣象,，2002(1).
[4] OETZEL G N,，PIERCE E T.Radio emissions from close lightning.Planetary Electrodynamics,，1 534-569，1969.
[5] 潘超.基于SoPC的雷電探測儀設(shè)計[J].微計算機信息,，2002.
[6] OETZEL G N,，PIERCE E T.VHF technique for locating lightning. Radio Science 4,，199-201，1969.
[7] RUSTAN P L,，UMAN M A,，CHIDERS D G，et al.Lightning source locations from VHF radiation data for a flash at kennedy space center[J].Geophys.Res.,，1985(C9)：4893-4903.