基于干擾保護比的同頻航向信標(biāo)臺最小保護間距計算*

程延松，李 銳，周自力，李 黎，王建強，王 軍

(1.中國民用航空局第二研究所，成都 610041；2.電子科技大學(xué) 通信抗干擾技術(shù)國家級重點實驗室，成都 611731)

?

基于干擾保護比的同頻航向信標(biāo)臺最小保護間距計算*

程延松**1，李 銳1，周自力1，李 黎2，王建強1，王 軍2

(1.中國民用航空局第二研究所，成都 610041；2.電子科技大學(xué) 通信抗干擾技術(shù)國家級重點實驗室，成都 611731)

針對航向信標(biāo)(LOC)臺站典型同頻場景，提出了最小保護間距的數(shù)值計算方法。利用 LOC臺站間距離曲線擬合得到滿足干擾保護比的保護間距函數(shù)，根據(jù)幾何函數(shù)確定LOC臺站空間間距，基于保護間距和空間間距的關(guān)系確定同頻LOC臺站的最小保護間距。所提計算方法可在滿足干擾保護比的前提下，獲得更加精確的同頻LOC臺站間距，該間距小于美國聯(lián)邦航空管理局(FAA)規(guī)定的相對保守的臺站間距，可有效提高LOC系統(tǒng)的頻率利用效率。

同頻航向信標(biāo)臺站；頻率利用效率；干擾保護比；需求保護間距；空間間距

1 引 言

儀表著陸系統(tǒng)(Instrument Landing System,ILS)是應(yīng)用最廣泛的飛機精密進近和著陸引導(dǎo)系統(tǒng)，由提供航向引導(dǎo)的航向信標(biāo)(Localizer,LOC)和提供垂直引導(dǎo)的下滑信標(biāo)(Glideslope,GS)確定下滑軌跡[1-2]。為了確保它們的頻率在其頻率保護服務(wù)區(qū)(Frequency Protected Service Volume,FPSV)中不受干擾，需要對ILS臺站的間距進行限制[3-5]。由于LOC的FPSV覆蓋區(qū)涵蓋了相關(guān)聯(lián)的GS的FPSV，故通常只考慮LOC的保護間距即可。LOC的保護間距包含LOC/LOC同頻、LOC/LOC鄰頻、LOC/甚高頻全向信標(biāo)(Very High Frequency Omnidirectional Range,VOR)鄰頻等情況。其中，LOC/LOC同頻場景下要求的保護間距最大，是本文的研究對象。

目前，美國聯(lián)邦航空管理局(Federal Aviation Administration,FAA)已有規(guī)范對LOC保護距離進行了規(guī)定，并提出了查表和計算法兩種方法[8]。其中，查表法較為直觀，但數(shù)據(jù)較為保守，具有較大的冗余；計算法的精度高于查表法，但文獻[8]中LOC/LOC站間距曲線圖不包含所有場景的數(shù)據(jù)，不能直接得到任意場景下的站間距。

針對上述表格法和計算法確定站間距的不足，本文提出一種新的綜合計算方法，通過天線輻射角度和天線方向圖確定同頻LOC的幾種典型場景。對每種場景，首先利用 LOC/LOC站間距曲線圖擬合得到滿足干擾保護比需求的保護間距函數(shù)，然后根據(jù)幾何函數(shù)確定LOC/LOC站間空間間距，最后利用保護間距和空間間距的關(guān)系確定同頻LOC的最小保護距離。計算表明，相比于FAA規(guī)定的相對保守的保護距離，本文所提計算方法可在滿足干擾保護要求的前提下，獲得更加精確的同頻LOC臺站間距離，該距離小于FAA規(guī)定的站間距離，從而可有效提高LOC系統(tǒng)的頻率利用效率。

2 同頻LOC模型

2.1 同頻干擾保護比

為確保同頻LOC安全工作，接收機處的期望信號功率與干擾信號功率的比值應(yīng)滿足信號干擾比(Desired to Undesired,D/U)要求[9]，特別地，需要在其FPSV的關(guān)鍵點提供D/U的保護。

在實際計算中，通常利用D/U值得到等效信號比(Equivalent Signal Ratio,ESR)，進而通過查表得到保護間距[8]。ESR是由期望信號和干擾信號的功率與天線增益差異引入的與D/U的差值，計算方法如下[7]：

ESR=D/U-PD+PU-AD+AU+GU-GD

(1)

式中：PD和PU分別為期望信號載波功率(dBW)與干擾信號載波功率(dBW)；AD和AU分別為期望發(fā)射機天線增益(dBi)與干擾發(fā)射機天線增益(dBi)；GD和GU分別為期望發(fā)射機在考察位置相對于主瓣天線增益的相對天線增益(dBi)與干擾發(fā)射機相對天線增益(dBi)[8]。如圖 1所示，設(shè)A為期望LOC臺站，B為干擾LOC臺站，期望發(fā)射機的考察位置即是B臺站位置，干擾發(fā)射機的考察位置即是A臺站位置；關(guān)鍵點為相對于A臺站的FPSV場強最弱，同時相對于B臺站的FPSV距離最近的點；GD由A臺站到臨界點的連線與A臺站天線波束主瓣的夾角決定，GU由B臺站到臨界點的連線與B臺站天線波束主瓣的夾角Θ決定。

圖1 同頻LOC模型

GD和GU可通過天線方向圖得到。以常用的LOC標(biāo)準V環(huán)天線(Standard V-Ring Antenna)為例[8]，GD和GU如圖2所示。

圖2 天線輻射角度Θ和相關(guān)天線增益曲線圖Fig.2 Antenna radiation angle Θ and relative antenna gain curve

為了便于后續(xù)計算，將GU與角度Θ的關(guān)系利用天線方向圖進行擬合。針對圖2中的相對天線增益曲線，可分段擬合如表1所示。

表1 相對天線增益曲線擬合函數(shù)

Tab.1 The fitting function of relative antenna gain curve

Θ/(°)GU/dBi0～202～10-58×Θ+5410～90-739()×Θ-1503990～130-20130～17012×Θ-95170～17734×Θ-139177～180-6.5

計算ESR時，D/U和GD的數(shù)值是確定的，只需考慮GU的數(shù)值。利用上述GU與角度Θ的關(guān)系，即可通過式(1)得到ESR與Θ的函數(shù)關(guān)系。

2.2 保護距離與場景劃分

定義SD是通過文獻[8]中的需求保護間距曲線和對應(yīng)的ESR值得到的同頻需求保護間距。同頻需求保護間距即是可以建立同頻LOC臺站的間距，在同頻需求保護間距之外，不需要考慮同頻干擾的問題。

定義SP為干擾LOC臺站與期望LOC臺站間的空間間距。如圖1所示，空間間距為干擾LOC臺站到期望LOC臺站關(guān)鍵點的間距dU與期望LOC臺站到關(guān)鍵點的間距dD的和[8]：

SP=dU+dD。

(2)

如前文所述，需求保護距離SD可由ESR得到，而ESR與圖 2中的Θ有關(guān)。由圖 1可知，空間間距SP可基于幾何關(guān)系由Θ計算得到。顯然，當(dāng)空間間距SP變化時，Θ的值也隨之改變，進而導(dǎo)致需求保護間距SD的值發(fā)生變化。由于需求保護間距SD決定了同頻LOC能正常工作的最小間距，最小保護間距即是能滿足SD要求的SP的最小值。因此，確定最小保護間距的核心是圖 1中的角度Θ。

同頻LOC場景中，干擾LOC臺站可能位于期望LOC臺站的任何方向[6]，需要根據(jù)天線方向圖確定典型場景進行分析。以圖 2中LOC臺站的天線方向圖為例，0°～90°曲線單調(diào)遞減，90°～180°曲線單調(diào)遞增，故天線輻射強度在0°、90°、180°時分別處于最大、最小值，具有代表性；并且在實際LOC臺站建設(shè)中多處于正向、反向、平行正向和平行反向等方位，故只考慮正向、反向、平行正向和平行反向4種典型同頻LOC場景，如圖3所示。根據(jù)式(2)，下面給出對應(yīng)4類場景下的Θ與SP的函數(shù)，場景中的臺站功率均為13 dBW，天線增益均為12 dBi。

(a)干擾LOC站點在期望LOC站點左側(cè)

(b)干擾LOC站點在期望LOC站點右側(cè)

(c)干擾LOC站點與期望LOC站點平行同的

(d)干擾LOC站點與期望LOC站點平行反向

圖3 同頻LOC場景

Fig.3 The typical scenarios of cochannel LOC stations

場景1：干擾LOC臺站在期望LOC臺站左側(cè)，如圖3(a)所示。圖中，A為期望LOC臺站，其前方的扇形區(qū)域為其信號的覆蓋區(qū)域，在方位角±10°的覆蓋區(qū)內(nèi)，引導(dǎo)間距不小于25 n mile(46.3 km)，在方位±10°～35°的覆蓋區(qū)內(nèi)，引導(dǎo)間距不小于17 n mile(31.5 km)[1];B為干擾LOC臺站;N為期望LOC臺站覆蓋區(qū)域的關(guān)鍵點;如前所述，Θ為關(guān)鍵點與B臺站的連線和B臺站主瓣所形成的夾角。此時有

(3)

場景2：干擾LOC臺站在期望LOC臺站右側(cè)，如圖3(b)所示。此時有

(4)

場景3：干擾LOC臺站與期望LOC臺站平行同向，如圖3(c)所示。此時有

(5)

場景4：干擾LOC臺站與期望LOC臺站平行反向，如圖3(d)所示。此時有

(6)

針對上述所有場景，我們的目標(biāo)是在干擾LOC臺站B對期望LOC臺站A的同頻干擾滿足需求保護間距要求SD的前提下，計算最小的空間間距SP。

3 計算方法

3.1 需求保護間距SD

為了便于后續(xù)計算，將文獻[8]中給出的需求保護間距曲線圖中距地高度1 905 m，引導(dǎo)間距分別為17 n mile和25 n mile的點全都提取出來,分別得到距地高度1 905 m、引導(dǎo)間距17 n mile和引導(dǎo)間距25 n mile的ESR與需求保護間距之間的曲線圖(在實際場景中，當(dāng)ESR小于0 dBi時，對應(yīng)需求保護間距數(shù)值很小，可以忽略，所以舍棄ESR小于0的部分數(shù)值)。為了便于計算，將曲線進行擬合[10]，分別得到引導(dǎo)間距17 n mile和引導(dǎo)間距25 n mile的ESR與需求保護間距SD之間的函數(shù)關(guān)系式如下：

(1)引導(dǎo)間距17 n mile的ESR與需求保護間距SD的二階擬合函數(shù)為

(7)

(2)引導(dǎo)間距25 n mile的ESR與需求保護間距SD的二階擬合函數(shù)為

(8)

計算表明，17 n mile的二階擬合函數(shù)的歸一化均方誤差(Normalized Mean Square Error,NMSE)數(shù)值0.038 7，25 n mile的二階擬合函數(shù)的NMSE數(shù)值為0.025 4。ESR與需求保護間距之間的實際曲線和擬合曲線如圖4所示，可以看出實際曲線和擬合曲線吻合很好，在后續(xù)的計算中直接采用擬合曲線。

圖4 ESR與需求保護間距之間的實際數(shù)值和擬合曲線(17 n mile和25 n mile)

Fig.4 The practical values and the corresponding fitting curve of required protection distance with respect to ESR

3.2 空間距離SP

干擾LOC臺站與期望LOC臺站間的空間間距SP可通過計算dU得到?；痉椒ㄊ牵豪檬?3)～(6)得到的Θ與SP的函數(shù)關(guān)系，根據(jù)Θ與ESR的對應(yīng)關(guān)系，得到與ESR對應(yīng)的SP值。以場景1為例，當(dāng)ESR為+13 dBi時，由式(1)可得到此時GU為-20 dBi，由圖 2可得此時Θ值為90°～130°，通過式(3)可計算得到此時SP的值為27～30 n mile。4種場景下的ESR與SP間的關(guān)系如表2～5所示。

表2 ESR與同頻LOC空間間距SP對應(yīng)表(場景1)

Tab.2 ESR and corresponding cochannel LOC space distance(Scenario 1)

ESR/dBiSP/nmile+1327～30+1432+1733+2064+2387+26157+26.5204→∞

表3 ESR與同頻LOC空間間距SP對應(yīng)表(場景2)

Tab.3 ESR and corresponding cochannel LOC space distance(Scenario 2)

ESR/dBiSP/nmile+830～31+1133+1440+1753+2075+2187+21.5108→∞

表4 ESR與同頻LOC空間間距SP對應(yīng)表(場景3)

Tab.4 ESR and corresponding cochannel LOC space distance(Scenario 3)

ESR/dBiSP/nmile+13+∞+14151+1755+2037+2334+2632+2932

表5 ESR與同頻LOC空間間距SP對應(yīng)表(場景4)

Tab.5 ESR and corresponding cochannel LOC space distance(Scenario 4)

ESR/dBiSP/nmile+1339→∞+1437+1734+2032+2332+2631+26.531→17

利用表2～5可分別得到4種場景的ESR與同頻LOC空間間距SP的曲線。

在已知ESR與需求保護間距SD的曲線圖和ESR與同頻LOC空間間距SP的曲線圖之后，可將4種場景的需求保護間距SD曲線和同頻LOC空間間距SP的曲線繪制于同一圖內(nèi)進行比較，如圖5所示。

(a)場景1

(b)場景2

(c)場景3

(d)場景4

圖5 ESR與需求保護間距和空間間距曲線圖

Fig.5 ESR,required protection distance and space distance curves

3.3 最小保護距離計算

在干擾LOC臺站B對期望LOC臺站A的同頻干擾滿足要求的前提下，最小空間間距SP_min為4種場景下需求保護間距SD和同頻LOC空間間距SP相等時的值。由于空間間距與ESR之間的函數(shù)關(guān)系難以得到解析的表達式，本文采用迭代的數(shù)值方法計算最小間距，下面介紹基本計算方法。

首先,由圖5得到SD和SP相等時SP值的范圍(SP_1,SP_2)，并確定SP在該范圍內(nèi)的單調(diào)性，取該范圍的平均值SP_mid。根據(jù)不同的場景，將該平均值代入式(3)～(6)，可求得對應(yīng)的角度Θ；

其次，根據(jù)獲得的Θ值，代入表1，可計算出此時GU的值;

第三，將GU代入式(1)即可得到ESR的值。通過ESR與需求保護間距函數(shù)式(7)或式(8)，可得到此時需求的保護間距SD_mid；

第四，將SD_mid與空間間距SP_mid進行比較：

(1)如果SP_mid與SD_mid恰好相等，SP_mid就是最小間距SP_min;

(2)若SP_mid大于SD_mid，在空間間距SP曲線單調(diào)遞增的情況下，用SP_mid替換邊界值SP_2，在空間間距曲線單調(diào)遞減的情況下，用SP_mid替換邊界值SP_1;

(3)若SP_mid小于SD_mid，在空間間距SP曲線單調(diào)遞增的情況下，用SP_mid替換邊界值SP_1，在空間間距曲線單調(diào)遞減的情況下，用SP_mid替換邊界值SP_2;

(4)縮小SP_min的范圍，并進行迭代運算。由于精度問題，SP_mid和SD_mid并不一定能取到相等值，故當(dāng)最終得到的SP_mid與SD_mid之間的差值小于誤差ε，即可認為此時的SP_mid為滿足同頻干擾保護的最小間距SP_min。

4 最小保護距離計算實例

采用上一節(jié)的計算方法，對第2節(jié)中給出的4種場景進行計算得到4種典型場景下的空間間距,如表6所示,計算時誤差ε=0.001。

表6 同頻LOC保護間距

Tab.6 The protection distance of cochannel LOC

同頻LOC模型最小保護間距/nmileFAA保守保護間距/nmile干擾LOC臺站在期望LOC臺站左側(cè)(場景1)7897干擾LOC臺站在期望LOC臺站右側(cè)(場景2)98135干擾LOC臺站與期望LOC臺站平行同向(場景3)63124干擾LOC臺站與期望LOC臺站平行反向(場景4)57124

從表6可以發(fā)現(xiàn)：不同場景下計算得到的最小保護間距均小于由FAA給出的保守保護間距；處于同一水平線的同頻LOC的保護距離要遠遠大于平行的同頻LOC的最小保護距離；在相同情況下，同向的同頻LOC的保護距離要大于反向的同頻LOC的最小保護距離。這是因為：由圖2，處于同一水平線時，LOC天線增益要遠大于平行時的增益；同時，LOC天線正向的增益也強于反向的增益。較大的天線增益導(dǎo)致期望LOC臺站與期望干擾LOC臺站之間更易產(chǎn)生干擾，故需要更遠的保護距離。

5 結(jié)束語

本文提出了一種更精確的同頻LOC站間距計算方法,利用同頻LOC的天線方向圖和天線輻射角度確定同頻場景，并利用擬合的方式將LOC站間曲線轉(zhuǎn)化為滿足干擾保護比需求的保護間距函數(shù)。進一步，由不同場景中的LOC臺站頻率保護區(qū)、關(guān)鍵點和天線輻射角度的幾何關(guān)系確定了LOC空間間距，利用保護間距與空間間距的相互關(guān)系確定同頻LOC最小保護間距。具體實例驗證結(jié)果表明，與FAA所規(guī)定的站間距相比，計算得到的同頻LOC最小保護間距在滿足干擾保護要求的前提下，數(shù)值更加小，從而可進一步提高LOC系統(tǒng)的頻率利用效率。

需要說明的是，與本文分析中采用的典型場景相比，在一般場景的計算中，由于臺站方位角和朝向的不確定，故首先需要由干擾臺站到期望臺站覆蓋區(qū)域邊界的最近點確定關(guān)鍵點的位置，并將臺站方位角和朝向加入到關(guān)鍵點到干擾LOC臺站主瓣夾角的計算中，最終確定干擾發(fā)射機相對天線增益，即可采用本文的方法計算最小間距。

[1] 王世錦，王湛.機載雷達與通信導(dǎo)航設(shè)備[M].北京：科學(xué)出版社,2009. WANG Shijin，WANG Zhan.Airborne radar and communication navigation equipment[M].Beijing:Science Press,2009.(in Chinese)

[2] 馬存寶.民機通信導(dǎo)航與雷達[M].西安：西北工業(yè)大學(xué)出版社,2004:365-374. MA Cunbao.Civil aircraft navigation and radar communication[M].Xi′an:Northwestern Polytechnic University Press,2004:365-374.(in Chinese)

[3] 陳康.儀表著陸系統(tǒng)航向信標(biāo)的分析與研究[J].信息通信,2011(6):9-11. CHEN Kang.Research and analysis on the ground course beacon subsystem of ILS[J].Information & Communications,2011(6):9-1.(in Chinese)

[4] MH/T 4006.1-1998,航空無線電導(dǎo)航設(shè)備[S]. MH/T 4006.1-1998,Aeronautical radio navigation aids [S].(in Chinese)

[5] ICAO Inernational Standards and Recommendation Practices Annex 10，Aeronutical Telecommunications Volume 1：Radios Navigation Aids Sixth Edition[S].

[6] CCAR-123-2004-98TM,平行跑道同時儀表運行管理規(guī)定[S]. CCAR-123-2004-98TM,Regulations for the management of parallel tracks instrument operation simultaneously[S].(in Chinese)

[7] 陳亞青，羅亮，詹建明.窄距平行跑道運行新技術(shù)探討[J] .中國民用航空,2013,12(12)：39-40. CHEN Yaqing，LUO Liang，ZHAN Jianming.Discussion on the new technology of the narrow range parallel track[J].China Civil Aviation,2013,12(12)：39-40.(in Chinese)

[8] FAA order 6050.32b,Spectrum management regulations and procedures manual[S].

[9] 中國民用航空總局政策法規(guī)司.附件10[G]//國際民用航空條約匯編.北京：中國民航出版社,2005. Civil Aviation Adminstration of China.APPENDIX 10[G]//International Civil Aviation Convention.Beijing:China Civil Aviation Press,2005.(in Chinese)

[10] 石振東.誤差和曲線擬合[M].哈爾濱:哈爾濱工程大學(xué)出版社,2010. SHI Zhendong.Error and curve fitting[M].Harbin:Harbin Enigineering University Press,2010.(in Chinese)

程延松(1970—)，男，河南魯山人，2000年于四川大學(xué)獲理學(xué)博士學(xué)位，現(xiàn)為中國民用航空局第二研究所高級工程師，主要從事航空通信導(dǎo)航監(jiān)視技術(shù)等領(lǐng)域的研究工作；

CHENG Yansong was born in Lushan,Henan Province,in 1970.He received the Ph.D. degree from Sichuan University in 2000.He is now a senior engineer.His research concerns aviation communication,navigation and surveillance technology.

Email:ys.cheng@163.com

李 銳(1968—)，女，四川成都人，1999年于中國科學(xué)院成都計算機應(yīng)用研究所獲工學(xué)碩士學(xué)位，現(xiàn)為中國民用航空局第二研究所高級工程師，主要從事航空通信導(dǎo)航監(jiān)視技術(shù)等領(lǐng)域的研究工作；

LI Rui was born in Chengdu,Sichuan Province,in 1968.She received the M.S.degree from Chengdu Institute of Computer Application,Chinese Academy of Sciences in 1999.She is now a senior engineer.Her research concerns aviation communication,navigation and surveillance technology.

周自力(1971—)，男，重慶人，1993年于上海交通大學(xué)獲雙學(xué)士學(xué)位，現(xiàn)為中國民用航空局第二研究所高級工程師，主要從事航空通信導(dǎo)航監(jiān)視技術(shù)研究工作；

ZHOU Zili was born in Chongqing,in 1971.He received the dual B.S. degree from Shanghai Jiaotong University in 1993.He is now a senior engineer.His research concerns aviation communication,navigation and surveillance technology.

李 黎(1994—)，男，重慶人，2016年于電子科技大學(xué)獲學(xué)士學(xué)位，現(xiàn)為電子科技大學(xué)通信抗干擾重點實驗室碩士研究生，主要從事無線通信算法及其實現(xiàn)的研究工作；

LI Li was born in Chongqing,in 1994.He received the B.S. degree from University of Electronic Science and Technology of China in 2016.He is now a graduate student.His research concerns wireless communication algorithm and its implementation.

王建強(1977—)，男，山西沁縣人，2001年于四川大學(xué)獲工學(xué)學(xué)士學(xué)位，現(xiàn)為中國民用航空局第二研究所工程師，主要從事計算機應(yīng)用、航空無線電導(dǎo)航頻率工程等方面的研究；

WANG Jianqiang was born in Qinxian,Shanxi Province,in 1977.He received the B.S. degree from Sichuan University in 2001.He is now an engineer.His research concerns computer application,aeronautical navigational aid frequency engineering.

王 軍(1974—)，男，四川人，2009年于電子科技大學(xué)獲工學(xué)博士學(xué)位，現(xiàn)為電子科技大學(xué)通信抗干擾技術(shù)國家級重點實驗室教授、博士生導(dǎo)師，主要從事無線與移動通信技術(shù)方面的研究。

WANG Jun was born in Sichuan Province,in 1974.He received the Ph.D. degree from University of Electroic Science and Technology of China in 2009.He is now a professor and also the Ph.D. supervisor.His research interests include signal processing for wireless communications,air-to-ground communication systems,and cognitive radio.

Interference Protection Ratio Based Computation of Minimum Cochannel Localizer Station Protection Distance

CHENG Yansong1,LI Rui1,ZHOU Zili1,LI Li2,WANG Jianqiang1,WANG Jun2

(1.The Second Research Institute of Civil Aviation Administration of China,Chengdu 610041,China;2.National Key Laboratory of Science and Technology on Communications,University of Electronic Science and Technology of China,Chengdu 611731,China)

A numerical computation method for minimum protection distance is proposed for cochannel localizer(LOC) stations under several typical scenarios.For this proposed method,the function of required protection distance is fitted based on the distance curves between LOC stations,and the space distance curves corresponding to these scenarios are obtained from the geometrical function.The minimum LOC cochannel distance is determined based on the relationship between the required protection distance and the space distance.Compared with the relatively conservative protection distance provided by the Federal Aviation Administration(FAA) regulation,the proposed method can obtain more accurate and smaller LOC cochannel distances under the fundamental premise of safety.Therefore,the frequency utilization efficiency of LOC systems can be improved.

cochannel localizer station;frequency utilization efficiency;interference protection ratio;demand protection distance;physical distance

10.3969/j.issn.1001-893x.2016.07.004

程延松，李銳，周自力，等.基于干擾保護比的同頻航向信標(biāo)臺最小保護間距計算[J].電訊技術(shù)，2016,56(7):737-743.[CHENG Yansong,LI Rui,ZHOU Zili,et al.Interference protection ratio based computation of minimum cochannel localizer station protection distance[J].Telecommunication Engineering,2016,56(7):737-743.]

2016-01-15；

2016-03-15 Received date:2016-01-15；Revised date：2016-03-15

國家自然科學(xué)基金資助項目(61471099)；國家高技術(shù)研究發(fā)展計劃(863計劃)項目(2015AA01A705)；民航安全能力建設(shè)項目(TMSA1510)；高等學(xué)校博士學(xué)科點專項科研基金資助課題(博導(dǎo)類)(20130185110005)

Foundation Item:The National Natural Science Foundation of China(No.61471099);The National High-tech R&D Program of China(863 Program)(2015AA01A705);The Civil Aviation Safety Capacity Build Project(TMSA1510)；The Specialized Research Fund for the Doctoral Program of Higher Education of China(20130185110005)

V249

A

1001-893X(2016)07-0737-07

**通信作者：ys.cheng@163.com Corresponding author：ys.cheng@163.com