地基雷達探測臨近空間高超聲速目標優(yōu)化部署方法*

?

地基雷達探測臨近空間高超聲速目標優(yōu)化部署方法*

肖松，譚賢四，王紅，汪毅

(空軍預(yù)警學(xué)院，湖北 武漢430019)

摘要：為提高地基雷達對臨近空間高超聲速目標的探測能力，探討了一種實用高效的地基雷達優(yōu)化部署方法，對地基雷達探測臨近空間高超聲速目標的難點和數(shù)學(xué)模型進行了分析，提出了地基雷達探測臨近空間高超聲速目標的部署原則和量化指標，建立了地基雷達優(yōu)化部署模型，并進行了信息素引導(dǎo)性控制的蟻群算法設(shè)計。仿真結(jié)果表明，該方法可實現(xiàn)多種程式地基雷達的優(yōu)化部署，提高了地基雷達優(yōu)化部署效率和可操作性，為地基雷達探測臨近空間高超聲速目標的優(yōu)化部署提供了一種新思路。

關(guān)鍵詞：優(yōu)化部署；地基雷達；臨近空間；高超聲速；蟻群算法；探測

臨近空間高超聲速飛行器具有高高空、超高速、高機動和遠距離精確打擊等特性[1-2]，可用于天地往返運輸系統(tǒng)及遠程精確打擊系統(tǒng)[3]，這給防空預(yù)警系統(tǒng)提出了新的嚴峻挑戰(zhàn)[4]。如何優(yōu)化部署地基雷達，提高雷達探測此類目標的效能，滿足最低保障能力指標和構(gòu)成最大保障能力探測網(wǎng)絡(luò)，實現(xiàn)雷達資源優(yōu)化配置，是當前亟待解決的問題。

針對地基雷達探測常規(guī)目標優(yōu)化部署方法的研究較多[5]，但是，針對臨近空間高超聲速目標的地基雷達優(yōu)化部署方法還未見報道。文獻[6]在充分考慮雷達部署的各種限制條件基礎(chǔ)上，采用遺傳算法，通過染色體群世代更迭，使部署接近最優(yōu)(次優(yōu))解；文獻[7]建立了環(huán)形、直線和扇形3種最優(yōu)部署模型，采用完全枚舉法、黃金分割法

和坐標輪換法求出最優(yōu)解；文獻[8]針對防區(qū)內(nèi)分散于不同位置且重要度不同的責(zé)任區(qū)，基于NSGA-Ⅱ算法進行多目標優(yōu)化，實現(xiàn)全方位全縱深的雷達部署；文獻[9]通過建立目標優(yōu)化函數(shù)和利用結(jié)構(gòu)演變更新青蛙位置，實現(xiàn)組網(wǎng)雷達的優(yōu)化部署。上述方法都是基于一定原理，對目標進行搜索，通常適用于大量部署方案中的尋優(yōu)過程，針對性強，而難以滿足雷達探測臨近空間高超聲速目標時，發(fā)現(xiàn)距離遠、連續(xù)穩(wěn)定跟蹤時間長和銜接高度特殊等需求。鑒于此，結(jié)合臨近空間高超聲速目標特性和地基雷達部署原則，擬提出地基雷達探測臨近空間高超聲速目標的約束條件并建立優(yōu)化部署數(shù)學(xué)模型，設(shè)計一種基于改進蟻群算法的優(yōu)化部署方法，為臨近空間預(yù)警系統(tǒng)中地基雷達優(yōu)化部署提供理論支撐。

1問題描述

1.1　地基雷達探測臨近空間高超聲速目標

雷達天線架設(shè)高度以及目標飛行高度與雷達直視距離之間的關(guān)系如圖1所示。

圖1　雷達直視距離與目標高度關(guān)系Fig.1　Relationship of radar direct-viewdistance and target height

從圖1中可以看出，理想情況下，對于飛行高度為20km的目標，地基雷達天線架高為1.0km時，地基雷達直視距離只有約700km。所以，對于飛行距離可以達到幾千公里的臨近空間高超聲速目標，必須采用多部雷達進行組網(wǎng)探測。

另外，臨近空間高超聲速飛行器在高超聲速飛行時，飛行器與大氣強烈作用，在飛行器周圍會形成等離子體包覆流場[10]，而包覆流場會對雷達電磁波產(chǎn)生衰減[11-12]。文獻[11]把單位距離的衰減系數(shù)公式簡化為

(1)

(2)

式中, ρ為反射系數(shù)。

(3)

(4)

其中，s為

(5)

式中，p=fp/f為fp相對f的歸一化值，q=v/2πf為等離子體碰撞頻率的歸一化。

總回波功率是由電磁波在空氣與等離子體界面上的反射功率加上反射目標回波功率等組成的，假設(shè)等離子體參數(shù)不變，等離子體厚度為10cm，此時雙程衰減為20cm，考慮fp=108~1012Hz，v/2π=107~1011Hz，C波段回波總衰減如圖2所示。

(a)回波總衰減(a) Echo total attenuation

由圖2可得，C波段條件下，回波總衰減最嚴重為22dB左右，總衰減大于5dB的范圍約為總感興趣范圍的1/16(大致處于高度20～60km，速度8～15Ma的范圍)。同理可得P波段條件下，回波總衰減小于8.5dB，總衰減大于2dB的范圍約為總感興趣范圍的1/8(大致處于高度20~55km，速度7.6~10.3Ma的范圍)。

由于等離子體對電磁波的衰減，使得雷達探測臨近空間高超聲速目標的難度進一步加大。所以，要通過對指定區(qū)域多部雷達裝備進行優(yōu)化部署，充分發(fā)揮雷達組網(wǎng)系統(tǒng)對臨近空間高超聲速目標的探測效能。

1.2　優(yōu)化部署數(shù)學(xué)模型

地基雷達優(yōu)化部署問題的數(shù)學(xué)模型為

(6)

式中，F(xiàn)i(X) 為目標函數(shù)， fi(x)為約束函數(shù)，X=(x1,x2,…,xn)T為n維決策變量空間，Rm為m維解空間。不存在X使得所有子目標函數(shù)Fi(X)同時取最大值，即目標函數(shù)相互之間是沖突的。

地基雷達探測臨近空間高超聲速目標優(yōu)化部署問題可以轉(zhuǎn)化為數(shù)學(xué)上的多邊形對多邊形的覆蓋問題，用數(shù)學(xué)語言可以描述如下：針對來襲目標的可能發(fā)射地點和攻擊對象，確定目標可能經(jīng)過的區(qū)域Q，其面積為S，對于不同類型雷達Ri(i=1,2,…,n，表示不同雷達程式其探測范圍可能為圓形、圓環(huán)、扇形等)，假設(shè)第i種程式的雷達對h高度層的覆蓋范圍為Sih，第j部雷達的中心記為(xij,yij)。該問題轉(zhuǎn)換為求解一種優(yōu)化方案，使得雷達對多邊形的覆蓋率α達到某一數(shù)值α0，即α≥α0，雷達探測面積重疊率β小于或等于某一數(shù)值β0，即β≤β0，同時，使用雷達的數(shù)量m(m≥m0)最小，所需雷達的最少值m0為理想值，即區(qū)域所需覆蓋的面積除以單個平臺覆蓋面積。要使部署方案效益最大，即問題解目標為

(7)

式中，ξh為h高度層的權(quán)重，約束條件為α≥α0，β≤β0，(xj,yj)∈Q。在工程應(yīng)用中，各目標之間常常難以很好兼顧，所以，采用改進的蟻群算法，通過建立多個優(yōu)化目標和評價指標，力求尋找較為均衡的優(yōu)化方案，實現(xiàn)對臨近空間高超聲速目標的有效探測。

2地基雷達探測臨近空間高超聲速目標優(yōu)化部署

2.1　量化指標

地基雷達探測臨近空間高超聲速目標的優(yōu)化部署指標主要有首次發(fā)現(xiàn)距離、銜接高度、探測范圍重疊度、責(zé)任區(qū)覆蓋系數(shù)和頻率重疊系數(shù)等。

1)首次發(fā)現(xiàn)距離

首次發(fā)現(xiàn)距離Ls指雷達網(wǎng)能有效探測到來襲目標時，目標與保護對象之間的最遠距離。表達式為

(8)

式中，Lp為雷達網(wǎng)在p高度層對目標的首次發(fā)現(xiàn)距離。由于臨近空間高超聲速目標飛行速度極快，為了能給攔截系統(tǒng)提供更充足的準備時間，要求首次發(fā)現(xiàn)距離越遠越好，如果沒有提出具體要求，Ls應(yīng)能滿足攔截武器在預(yù)定攔截窗口實施攔截的要求。

2)威力銜接高度

臨近空間高超聲速目標能夠在臨近空間進行跳躍飛行，為了保證雷達網(wǎng)在臨近空間各高度層都具有較好的威力銜接，沒有具體規(guī)定時，應(yīng)以來襲目標飛行的最低高度來進行約束。銜接高度示意圖如圖3所示。

圖3　銜接高度示意圖Fig.3　Connect height graph

3)探測范圍重疊度

相鄰雷達探測范圍重疊度Ch是指相鄰雷達在h高度層上重疊的最大跨度與平均探測距離的比值。重疊度示意圖如圖4所示。

圖4　重疊度示意圖Fig.4　Graph degree of overlap

則相鄰雷達探測范圍重疊度可表示為

(9)

式中，r1、r2分別為雷達1和2的探測距離，chmax為h高度層上相鄰雷達重疊的最大寬度。相鄰雷達在進行任務(wù)交接時，需滿足3點快速航跡起始的條件，即chmax≥3vT(v為目標速度，T為雷達探測周期)。

4)責(zé)任區(qū)覆蓋系數(shù)

責(zé)任區(qū)覆蓋系數(shù)?是指雷達網(wǎng)覆蓋區(qū)域與責(zé)任區(qū)域面積之比。在戰(zhàn)場環(huán)境下，可將責(zé)任區(qū)分為警戒區(qū)域、重點區(qū)域和核心區(qū)域。為了提高雷達網(wǎng)對臨近空間高超聲速目標的連續(xù)探測能力，要求擴大警戒區(qū)域，保障重點區(qū)域，重點保護核心區(qū)域。

(10)

5)頻率重疊系數(shù)

頻率重疊系數(shù)φ是指相鄰雷達的頻域重疊大小和空域重疊大小的關(guān)系。為了避免同頻干擾，相同頻率雷達部署時要保持相對距離，并在雷達網(wǎng)內(nèi)配置不同頻率的雷達。

(11)

式中，Sih∩Sjh為第i部雷達和第j部雷達在第h高度層的探測重疊面積；fih∩fjh為第i部雷達和第j部雷達在第h高度層上頻率重疊大??；Sih∪Sjh為第i部雷達和第j部雷達在第h高度層的探測面積總和；fih∪fjh為第i部雷達和第j部雷達在第h高度層上頻帶寬度總和。

2.2　部署模型

以上量化指標可以作為優(yōu)化部署方案選取的主要依據(jù)，地基雷達探測臨近空間高超聲速目標優(yōu)化部署的目標函數(shù)可以表示為多個高度層的優(yōu)化組合函數(shù)，即

(12)

Fs=f(Ls,H,Ch,?,φ,ρ)

(13)

式中，Ls、H、Ch、?、φ和ρ分別為首次發(fā)現(xiàn)目標距離、銜接高度、探測范圍重疊度、責(zé)任區(qū)覆蓋系數(shù)、頻率重疊系數(shù)和其他約束條件。綜合考慮以上因素和雷達網(wǎng)部署原則，以及探測臨近空間高超聲速目標的特殊需求，優(yōu)化部署的目標函數(shù)為

(14)

式中，ξh為各個高度層優(yōu)化指數(shù)所占比重，F(xiàn)sh為h高度層上雷達網(wǎng)連續(xù)覆蓋面積。約束條件為

(15)

式中，L0、H0、C0和ρ0分別為發(fā)現(xiàn)目標距離、銜接高度、探測范圍重疊度和其他約束條件的常數(shù)值，?→1、φ→0分別表示責(zé)任區(qū)覆蓋系數(shù)越大越好、頻率重疊系數(shù)越小越好。另外，還有一些約束條件無法進行數(shù)學(xué)描述，如氣候條件、地形條件，因此，優(yōu)化部署方案最后還需要根據(jù)作戰(zhàn)實際情況進行微調(diào)。

3基于改進蟻群算法的優(yōu)化部署算法

3.1　優(yōu)化部署步驟

根據(jù)作戰(zhàn)需要、目標特性以及雷達性能等因素，對雷達網(wǎng)進行優(yōu)化部署，采用固定部署和機動部署相結(jié)合的方式，主要分為以下幾步：

1)劃分防御范圍

臨近空間高超聲速飛行器的發(fā)射點所在區(qū)域相對固定，根據(jù)其作戰(zhàn)意圖，可以預(yù)測其落點區(qū)域，如圖5所示。

圖5　防御區(qū)域示意圖Fig.5　Defense area graph

在圖5中，目標發(fā)射區(qū)域為以O(shè)點為圓心的區(qū)域，目標可以在敵國陸地、島嶼或者艦船上發(fā)射，保衛(wèi)目標為P點所在的圓形區(qū)域，核心區(qū)域為來襲目標最有可能經(jīng)過的區(qū)域(圖中頂點為J、C、E、L)，重點區(qū)域為在目標機動范圍內(nèi)可能攻擊周邊其他目標的區(qū)域(圖中頂點為I、B、F、M)，防御區(qū)域為由于環(huán)境、裝備誤差等，造成對保衛(wèi)區(qū)域預(yù)測存在偏差，或者在重點區(qū)域附近存在需要保衛(wèi)的目標，對重點區(qū)域進行擴大而形成的防御范圍(圖中頂點為O、A、G)。

2)預(yù)估雷達數(shù)目

根據(jù)作戰(zhàn)需要、防御范圍和雷達性能等，確定雷達類型及數(shù)量。結(jié)合圖1和圖3，確定相鄰雷達威力銜接高度H以及相鄰雷達在h高度層上的最大探測距離，得到相鄰雷達之間的最大布站距離，并估算完全覆蓋防御區(qū)域的雷達數(shù)量。

3)地基雷達部署的初始化，設(shè)置一種部署方案作為初始方案。

4)根據(jù)約束條件計算目標函數(shù)。

5)調(diào)整部署方案，對于高機動雷達，可以根據(jù)機動部署原則在優(yōu)化過程中進行考慮。

6)重復(fù)步驟4和步驟5。

7)根據(jù)優(yōu)化函數(shù)確定可行的部署方案，并選取最優(yōu)或次優(yōu)的部署方案，搜索結(jié)束。

通常在得到優(yōu)化部署方案后，還需要根據(jù)實際情況進行人工微調(diào)。

3.2　改進蟻群算法

優(yōu)化算法是雷達網(wǎng)優(yōu)化部署需要考慮的一個問題，力求節(jié)省尋優(yōu)步驟，快速搜索到最佳方案。改進的蟻群算法將雷達優(yōu)化部署問題轉(zhuǎn)換為蟻群網(wǎng)絡(luò)，通過對初始化信息素和每次搜索結(jié)束后的信息素進行優(yōu)化，提高蟻群搜索效率。改進蟻群算法流程圖如圖6所示。

圖6　改進蟻群算法流程圖Fig.6　Modified ant colony algorithm flow chart

如圖6所示，初始化包括確定防御區(qū)域、核心區(qū)域范圍和雷達相關(guān)參數(shù)等。搜索路徑上的每個點如果滿足所有約束條件，則可形成一個部署方案，當m只螞蟻隨機放在n(m≥n)個部署點上時，為了使螞蟻快速從大量路徑中找到較好的路徑，k時刻待選部署點的概率P(k)為

P(k)=τ[s(k), j(k)]α·φ[s(k), j(k)]β

(16)

式中，τ、φ分別為信息素矩陣和啟發(fā)因子矩陣，α、β分別為信息素重要程度參數(shù)和啟發(fā)式因子重要程度參數(shù)，s(k)、j(k)分別表示k時刻已訪問部署點和待訪問部署點。

(17)

式中，ε、φ分別為信息素蒸發(fā)系數(shù)和增加強度系數(shù)，λ(k)為k時刻部署點與前一點的距離，ξ為當前部署點在防御區(qū)域中重要程度的參數(shù)。對處于核心區(qū)域的螞蟻給出啟發(fā)性的信息素，用以引導(dǎo)螞蟻選擇路徑，剔除處于防御區(qū)域邊緣的螞蟻，因為，如果雷達中心位于防御區(qū)域邊緣上，則其一半的探測范圍在防御區(qū)域之外，對探測資源是一種浪費。另外，螞蟻每次搜索結(jié)束后，對其搜索路徑進行分析，比較好的部署方案則留下更多的信息素，減少下一次無效的搜索。

4仿真分析

4.1　仿真環(huán)境

假設(shè)：來襲臨近空間高超聲速目標距離保衛(wèi)目標8000km，巡航段速度近似為勻速v≈15Ma，飛行高度h≈30km，防御區(qū)域大小為3500×6000km2，雷達程式、部署數(shù)量估算及威力(以30km高度為例)如表1所示。

表1　雷達程式、部署數(shù)量及威力

在工程實踐中，地基雷達存在頂空盲區(qū)，但本文探討的是骨干雷達優(yōu)化部署方式，其盲區(qū)可以通過相應(yīng)的補盲雷達來彌補，而雷達探測臨近空間目標不受地物反射干擾，氣候環(huán)境相對較穩(wěn)定，且防御區(qū)域遠離敵方發(fā)射區(qū)域，難以實施干擾，故這里假設(shè)雷達探測區(qū)域近似為圓形。

4.2　結(jié)果分析

不同迭代次數(shù)尋優(yōu)結(jié)果如圖7所示(由于篇幅原因，圖中橫縱坐標的真實比例難以很好顯示)。

(a)第20次尋優(yōu)結(jié)果(a)The 20th. optimizing result

(b) 第40次尋優(yōu)結(jié)果(b)The 40th. optimizing result

(c) 第60次尋優(yōu)結(jié)果(c)The 60th. optimizing result圖7　不同迭代次數(shù)仿真結(jié)果Fig.7　Simulation results of different iteration frequency

從圖7可以看出，優(yōu)化部署是蟻群根據(jù)約束條件，不斷尋找和判斷最佳部署點的過程，尋優(yōu)路線及長度如圖8所示。

圖8　尋優(yōu)路線及長度Fig.8　Optimizing path and length

圖8中，左圖為螞蟻在一次尋優(yōu)過程中歷經(jīng)的路線，右圖為各次最佳尋優(yōu)路線長度和尋優(yōu)路線平均長度。從圖中可以看出，隨著迭代次數(shù)的增加，最佳尋優(yōu)路線長度在不斷減小，說明通過信息素的引導(dǎo)作用，蟻群在不斷接近最優(yōu)解。量化指標在計算過程中的變化情況如圖9所示。

(a) 平均探測范圍重疊度和首次發(fā)現(xiàn)距離(a) Overlap degree of average detection rangeand first find distance

(b) 責(zé)任區(qū)覆蓋系數(shù)和頻率重疊系數(shù)(b) Obligation area covering factor and frequencyoverlap factor圖9　量化指標變化情況Fig.9　Variation of quantize index

從仿真結(jié)果來看，隨著迭代次數(shù)的增加，平均探測范圍重疊度和頻率重疊系數(shù)在逐漸減小，首次發(fā)現(xiàn)距離和責(zé)任區(qū)覆蓋系數(shù)在逐漸增大。當?shù)螖?shù)超過5時，責(zé)任區(qū)覆蓋系數(shù)增大明顯，首次發(fā)現(xiàn)距離和頻率重疊系數(shù)接近平均值，這說明通過信息素的有效引導(dǎo)，大大縮小了蟻群搜索尋優(yōu)的解空間，提高了優(yōu)化部署的實施效率。

5結(jié)論

針對地基雷達探測臨近空間高超聲速目標的優(yōu)化部署問題，提出地基雷達量化指標，建立優(yōu)化部署模型，并改進蟻群算法中信息素的引導(dǎo)性。通過仿真得出以下結(jié)論：

1)新方法針對性強，能夠?qū)崿F(xiàn)多種程式地基雷達探測臨近空間高超聲速目標的優(yōu)化部署，為地基雷達優(yōu)化部署提供了一種新思路。

2)通過對信息素進行引導(dǎo)性的控制，能夠有效提高蟻群算法的尋優(yōu)速度，減少無效搜索次數(shù)。

臨近空間高超聲速目標與常規(guī)目標同時存在時，如何對地基雷達以及異源傳感器進行優(yōu)化部署，將是下一步的研究重點。

參考文獻(References)

[1]Xu Y F, Jiang B, Gao Z F, et al. Fault tolerant control for near space vehicle: a survey and some new results[J]. Journal of Systems Engineering and Electronics, 2011, 22(1):88-94.

[2]Huang W, Ma L, Wang Z G, et al. A parametric study on the aerodynamic characteristics of a hypersonic wave rider vehicle[J]. Acta Astronautica, 2011, 69(3-4):135-140.

[3]曾開春,向錦武. 高超聲速飛行器飛行動力學(xué)特性不確定分析[J]. 航空學(xué)報, 2013, 34(4): 798-808.

ZENG Kaichun, XIANG Jinwu. Uncertainly analysis of flight dynamic characteristics for hypersonic vehicles[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(4):798-808. (in Chinese)

[4]肖松,譚賢四,王紅,等. 臨近空間高超聲速目標斷續(xù)點跡航跡起始方法[J].華中科技大學(xué)學(xué)報(自然科學(xué)版), 2014, 42(3):52-57.

XIAO Song, TAN Xiansi, WANG Hong, et al. Feasible track initiation method for near space hypersonic target[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2014, 42(3):52-57. (in Chinese)

[5]馬梁,王雪松,李永禎,等. 雷達部署對進動目標微多普勒頻率的可測性分析[J].國防科技大學(xué)學(xué)報, 2011, 33(2):54-59.

MA Liang, WANG Xuesong, LI Yongzhen, et al. Influence of radar location on the measurability of micro-doppler detecting for precession target[J]. Journal of National University of Defense Technology, 2011, 33(2):54-59. (in Chinese)

[6]張遠,方青,曲成華.基于遺傳算法的組網(wǎng)雷達優(yōu)化部署[J].雷達科學(xué)與技術(shù), 2014, 12(1):76-80.

ZHANG Yuan, FANG Qing, QU Chenghua. Optimal deployment of netting radars based on genetic algorithms[J]. Radar Science and Technology, 2014, 12(1):76-80. (in Chinese)

[7]Gao S. Optimal deployment problems of radar network[J]. Research Journal of Applied Sciences, Engineering and Technology, 2013, 6(10): 1879-1883.

[8]劉彥君,黃金才,王江. 有源干擾條件下基于NSGA-Ⅱ的雷達網(wǎng)優(yōu)化部署方法[J]. 指揮控制與仿真, 2014, 36(1):36-40.

LIU Yanjun, HUANG Jincai, WANG Jiang. Optimal deployment of radar network based on NSGA-Ⅱ under active jamming[J]. Command Control & Simulation, 2014, 36(1):36-40. (in Chinese)

[9]Cui Y J, Cha H, Tian B. Cultural shuffled frog leaping algorithm and its applications for the radar network deployment[J]. Applied Mechanics and Materials, 2014, 624:516-519.

[10]樂嘉陵. 再入物理[M]. 北京:國防工業(yè)出版社, 2005.

YUE Jialing. Re-entry physical[M].Beijing:National Defense Industry Press, 2005. (in Chinese)

[11]Garg P, Dodiyal A K. Reducing RF blackout during re-entry of the reusable launch vehicle[C]//Aerospace Conference, IEEE, 2009:1-15.

[12]Pippert G F. On the structure of wake turbulence deduced from field radar measurements[C]//AIAA Conference on Physics of Entry into Planetary Atmosphere, 1963:6320446.

Deployment optimization method for ground-based radar to detect near space hypersonic target

XIAOSong,TANXiansi,WANGHong,WANGYi

(Air Force Early Warning Academy, Wuhan 430019, China)

Abstract：In order to improve detection ability of ground-based radar to near space hypersonic target and explore a practical and efficient deployment optimization method, the difficulties and mathematical models of the ground-based radar to detect near space hypersonic target were analyzed and the deployment principle and quantitative indicators of it were put forward. Deployment optimization model of ground-based radar was established and the ant colony algorithm controlled by introductory pheromone was designed. The simulation results show that this method can achieve a variety of programs of ground-based radar deployment optimization, improve the efficiency and operability of the ground-based radar deployment optimization, and provide a new way of thinking for ground-based radar to detect near space hypersonic target.

Key words：deployment optimization; ground-based radar; near space; hypersonic; ant colony algorithm; detection

中圖分類號：TN95

文獻標志碼：A

文章編號：1001-2486(2015)02-121-07

收稿日期：2014-08-11基金項目：國家自然科學(xué)基金資助項目(61271451)

作者簡介：肖松(1985—)，男，湖北武漢人，博士研究生，E-mail：tracyoioi@163.com；譚賢四(通信作者)，男，教授，博士，博士生導(dǎo)師， E-mail: tanxs-hust@163.com

doi:10.11887/j.cn.201502023

http://journal.nudt.edu.cn