孟慶松, 王 彬, 邵高平

(信息工程大學(xué)信息工程學(xué)院,河南 鄭州 450001)

0 引 言

線性調(diào)頻(linear frequency modulation,LFM)信號具有較好的自相關(guān)特性和較寬的多普勒容限,在水下經(jīng)常被用作水聲通信同步信標(biāo)、魚雷或主動聲納的探測信號。在水聲偵察和對抗過程中,識別出接收信號為LFM信號,有助于水聲通信信號的捕獲、信息恢復(fù)和水下目標(biāo)探測。

目前,對水聲LFM信號的識別方法多是假設(shè)在加性高斯白噪聲(additive Gaussian white noise, AGWN)環(huán)境下,比如基于魏格納分布的方法[1-2]、基于短時傅里葉變換(short-time Fourier transform,STFT)的方法[3-5]、基于自相關(guān)和循環(huán)譜[6-7]的方法及基于分?jǐn)?shù)階傅里葉變換(fractional Fourier transform,FRFT)[8-10]的方法。事實上,水聲信號傳播中受水聲多途效應(yīng)的影響,正確識別水聲信號必須選擇對水聲多途效應(yīng)穩(wěn)健的信號特征,文獻(xiàn)[11]討論了存在水聲多途信道和AGWN條件下的水聲信號識別問題,提出了對水聲多途效應(yīng)穩(wěn)健的水聲LFM循環(huán)譜特征。然而,在淺海水聲環(huán)境中,由于工業(yè)活動、地震風(fēng)暴及海洋生物等因素的影響,會存在大量的脈沖干擾,這種噪聲通常用α穩(wěn)定分布來描述[12-15]。α穩(wěn)定分布噪聲具有明顯的尖峰脈沖噪聲特性,不存在有限的二階及以上各階統(tǒng)計量[16],使得上述算法的檢測性能明顯下降。

近年來,在α穩(wěn)定分布噪聲下對LFM信號的處理問題逐漸受到人們的關(guān)注[17-19],但公開研究成果并不多見,文獻(xiàn)[17]利用了基于STFT的方法,通過設(shè)置損失函數(shù)對傳統(tǒng)的STFT進(jìn)行優(yōu)化,并進(jìn)行Radon變換,通過設(shè)定的閾值對LFM信號進(jìn)行檢測,但該方法在低信噪比(signal-to-noise ratio，SNR)條件下性能嚴(yán)重下降,且文中并未考慮信道條件。2016-2017年,文獻(xiàn)[18-19]在進(jìn)行α穩(wěn)定分布噪聲下的LFM信號參數(shù)估計時使用了基于穩(wěn)健S變換和最優(yōu)L-柯西加權(quán)的方法,取得了一定的改進(jìn)效果。

為了提高淺海水聲環(huán)境下LFM信號的識別性能,提出了基于非線性變換的離散分?jǐn)?shù)階傅里葉變換(discrete fractional Fourier transform,DFRFT)方法,用以檢測淺海水聲LFM信號。該方法首先對接收信號進(jìn)行非線性變換,然后進(jìn)行DFRFT,構(gòu)造出LFM信號識別特征量,最后利用支持向量機(jī)(support vector machine,SVM)作為分類器,對LFM信號進(jìn)行識別。仿真實驗結(jié)果表明,在混合信噪比(mixed signal-to-noise ratio,MSNR)為-15 dB時正確識別率高于94%。

1 信號模型

α穩(wěn)定分布噪聲下接收到的淺海水聲信號模型為

y(t)=s(t)?h(t)+w(t)

(1)

式中,符號“?”表示卷積運(yùn)算;y(t)為接收信號;s(t)為發(fā)送信號,其調(diào)制類型為集合{LFM,MPSK,MFSK,OFDM}的任意一種。如果發(fā)送信號為LFM信號,其表達(dá)式為

s(t)=Aexp(jπkt2+2jπf0t)

(2)

式中,A表示LFM信號幅度;k為LFM信號的調(diào)頻斜率;f0是初始頻率。

若發(fā)送信號為多相移鍵控(multiple phase shift keying，MPSK)信號,其表達(dá)式為

(3)

若發(fā)送信號為多頻移鍵控(multiple frequency shift keying,MFSK)信號,其表達(dá)式為

(4)

式中,fΔ為頻率偏移量;Δf=2fΔ為頻率間隔;sl表示傳輸符號。

若發(fā)送信號為正交頻分復(fù)用(orthogonal frequency division multiplexing,OFDM)信號,其表達(dá)式為

(5)

h(t)為水聲多途信道的沖激響應(yīng),其表達(dá)式為

(6)

式中,N為多徑數(shù)目,當(dāng)i=0,τ0=0時為直達(dá)聲道。τi(i≠0)為第i條聲線相對于第0條聲線的延時,Ai為第i條聲線的幅度,fdi為第i條聲線的多普勒頻移。

w(t)為與發(fā)送信號s(t)相獨立的噪聲信號,服從α穩(wěn)定分布。參數(shù)α為特征指數(shù),0<α≤2,脈沖性隨著α的增大而減小,當(dāng)α=2時脈沖噪聲退化為高斯噪聲。圖1為不同α值所對應(yīng)的標(biāo)準(zhǔn)對稱α穩(wěn)定(symmetricαstable,SαS)分布噪聲[16]的時域波形與時頻圖。

圖1 不同α值的標(biāo)準(zhǔn)SαS分布噪聲及時頻圖Fig.1 Standard SαS distributed noise and time-frequency diagram in different α

可以看出,隨著特征指數(shù)的減小,脈沖數(shù)量增多、幅度增強(qiáng)。當(dāng)α=2時,噪聲變成平穩(wěn)的高斯噪聲。從時頻圖1可以看出,α穩(wěn)定分布噪聲與高斯白噪聲(white Gaussian noise，WGN)具有明顯的區(qū)別。α穩(wěn)定分布噪聲在時頻域上會出現(xiàn)帶狀的頻域很寬的能量分布,且在時域上分布不均勻,而WGN在時頻域上的能量是均勻分布的。

2 基于非線性映射DFRFT的水聲LFM信號識別方法

本文針對淺海水聲多途脈沖噪聲信道環(huán)境,以典型水聲信號集{LFM,MPSK,MFSK,OFDM}為待識別信號集,研究LFM信號的識別問題,提出了基于非線性映射DFRFT的LFM信號識別算法,算法的處理框圖如圖2所示。

圖2 基于非線性映射DFRFT的LFM信號識別算法流程Fig.2 LFM signal recognition algorithm process based on DFRFT ofnonlinear mapping

2.1 對接收信號進(jìn)行非線性變換

α穩(wěn)定分布噪聲具有明顯的尖峰脈沖特性,可利用非線性函數(shù)對接收信號進(jìn)行非線性變換,對脈沖噪聲進(jìn)行抑制。首先對接收信號y(t)定義非線性變換[20]，即

(7)

將式(1)代入式(7),有

(8)

2.2 基于DFRFT的LFM信號特征提取

根據(jù)FRFT的定義式[21],對非線性變換后的接收信號進(jìn)行FRFT,有

(9)

式中,P為階數(shù);K(t,u)為核函數(shù),定義[16]為

KP(t,u)=

(10)

式中,n為整數(shù);u為變換域;β表示旋轉(zhuǎn)角度,β與P的關(guān)系為β=Pπ/2,將式(2)代入式(10)可得

(11)

從式(11)中可知,FRFT可以理解為旋轉(zhuǎn)角度為β的時頻面旋轉(zhuǎn),當(dāng)階數(shù)變化到一定范圍時,LFM信號在對應(yīng)階次的FRFT域上具有集聚效應(yīng),出現(xiàn)較高的峰值,當(dāng)階數(shù)為最佳階數(shù)時,k=-cotβ,此時變換域出現(xiàn)的峰值最大。因此,可利用基于FRFT的方法對LFM信號進(jìn)行檢測。

由式(2)和式(6)可得接收信號為

y(t)=w(t)+A0s(t)exp(j2πfd0t)+

(12)

根據(jù)FRFT的性質(zhì)[22],有

FP[y(t)]=Np(u)+A0Fp[s(t)exp(j2πfd0t)]+

Np(u)+A0Sp(u-fd0sinβ)·

(13)

式中,FP[y(t)]是階數(shù)為P時y(t)的FRFT結(jié)果;Np(u)為噪聲部分的FRFT,若發(fā)送信號為LFM信號且Np(u)的影響較小時,從式中可知,當(dāng)階數(shù)為最佳階數(shù)時,在FRFT域會形成一系列峰值。

由于工程計算中為離散形式,計算時需要使用DFRFT方法,這里采用基于分解型的DFRFT[21],即

(14)

分別以步長ΔP和Δm對參數(shù)P和m進(jìn)行遍歷,其中ΔP=0.005,Δm=1得到接收信號的FRFT信息。

將非線性變換后的信號帶入DFRFT式(14)中,可得

(15)

現(xiàn)有的基于FRFT對LFM信號檢測的方法多通過使用閾值的方法來判定有無LFM信號[17],在低SNR和α穩(wěn)定分布噪聲下的檢測性能不穩(wěn)定,對閾值的依賴性較高。對此,本文提出了基于方差比的檢測識別方法,可有效提高檢測性能。

從式(13)可知,若發(fā)送信號為LFM信號且噪聲分量在分?jǐn)?shù)域的影響較小,當(dāng)階數(shù)為最佳階數(shù)時,在FRFT域會形成一系列峰值;當(dāng)階數(shù)為1時,FRFT退化為普通的傅里葉變換,接收信號在頻域表現(xiàn)為一個帶寬較大的信號。因此可將兩者的歸一化方差比作為特征量進(jìn)行LFM信號的識別。

(16)

(17)

式中,max(·)表示求最大值。

(18)

式中,max(·)表示求最大值。

通過分別計算這2種階數(shù)時信號在FRFT域上的方差,得到方差的比值并作為特征量,即

(19)

(20)

識別特征量定義為

(21)

2.3 通過SVM對接收信號的特征量進(jìn)行分類識別

SVM是一種有監(jiān)督的機(jī)器學(xué)習(xí)算法,可用于分類問題。他使用一種稱為核函數(shù)的技術(shù)找到2種分類之間的最佳邊界[23]。

不同的海域和不同的SNR條件會對LFM信號的識別特征量產(chǎn)生一定影響。與通過閾值進(jìn)行判決相比,使用SVM進(jìn)行分類可在不同條件下均找到最佳分類邊界,且可以通過不斷增加訓(xùn)練信號改進(jìn)最佳邊界。近年來,也有人提出一些實用的方法對SVM的性能進(jìn)行改進(jìn),如對數(shù)據(jù)處理方式進(jìn)行調(diào)整,大樣本類多取樣,小樣本類少取樣等[24-26]。

通過Matlab仿真軟件產(chǎn)生典型海洋水聲信號聲源波形數(shù)據(jù)。首先,產(chǎn)生不同調(diào)制方式的水聲通信信號(比如調(diào)制方式為MPSK、MFSK、OFDM等)、不同參數(shù)的LFM信號、不同種類的艦船輻射噪聲(比如不同螺旋槳數(shù)、不同螺旋槳葉片數(shù)的艦船)等。

根據(jù)全球海洋信息數(shù)據(jù)(比如Argo全球海洋信息數(shù)據(jù)集[27]),通過經(jīng)驗公式可以得到不同海域的聲剖面信息,即

c=1 449.2+4.6T-0.055T2+

(1.34-0.010T)(S-35)+0.016z

(22)

式中,c為聲速，m/s;T為溫度，℃;S為鹽度，‰;z為深度，m。

然后,利用Bellhop仿真軟件,得到不同深度、不同傳輸距離條件下的聲信道沖激響應(yīng)。各種聲源波形與不同信道的沖激響應(yīng)卷積運(yùn)算,然后再疊加上不同強(qiáng)度的噪聲,得到接收傳感器的接收波形,將此仿真信號波形作為SVM訓(xùn)練的樣本數(shù)據(jù)。

其中,疊加的噪聲為不同特征指數(shù)的α穩(wěn)定分布噪聲,用MSNR來刻畫其噪聲強(qiáng)度,即

(23)

選用SVM為分類器,根據(jù)待檢測信號所發(fā)生的海域,選取該海域或最臨近海域的樣本數(shù)據(jù),計算識別特征量,并對SVM分類器進(jìn)行訓(xùn)練。

然后,計算待測信號的識別特征量,并將其送入訓(xùn)練好的SVM分類器中進(jìn)行LFM信號的識別,即可完成在淺海水聲條件下對LFM信號的檢測。

3 仿真實驗及分析

為驗證本文方法的有效性,設(shè)計仿真實驗進(jìn)行驗證。仿真實驗條件為:噪聲為加性標(biāo)準(zhǔn)SαS 分布噪聲,噪聲特征指數(shù)為α=1.5。根據(jù)Argo全球海洋信息數(shù)據(jù)集,得到東海某海域的聲剖面信息,如圖3所示。

圖3 東海某海域聲剖面圖Fig.3 Acoustic profile of sea area in the East China Sea

將發(fā)送端、接收端的深度分別設(shè)置為200 m、300 m,發(fā)送端、接收端的水平距離為2 000 m,通過Bellhop仿真得到傳輸信道的沖激響應(yīng),如圖4所示。

圖4 東海某海域沖激響應(yīng)圖Fig.4 Impulse response diagram of sea area in the East China Sea

3.1 非線性變換改進(jìn)性能分析

首先,驗證α穩(wěn)定分布噪聲下非線性變換對LFM信號檢測性能改進(jìn)的有效性。仿真實驗條件為:發(fā)送LFM信號的初始頻率為13 kHz,帶寬為5 kHz,采樣頻率50 kHz,采樣時間為50 ms;信道采用上述Bellhop仿真得到的東海某海域信道;噪聲為加性標(biāo)準(zhǔn)SαS分布噪聲,噪聲特征指數(shù)為α=1.5。圖5(a)是LFM信號未經(jīng)過非線性變換后的DFRFT(MSNR=-5 dB),圖5(b)為LFM信號經(jīng)過非線性變換后的DFRFT(MSNR=-5 dB),圖5(c)是LFM信號未經(jīng)過非線性變換后的DFRFT(MSNR=-15 dB),圖5(d)是LFM信號經(jīng)過非線性變換后的DFRFT(MSNR為-15 dB)。

圖5 非線性變換對DFRFT結(jié)果的影響Fig.5 Influence of nonlinear transformation on DFRFT results

對比圖5可以看出,由于α穩(wěn)定分布噪聲的影響,傳統(tǒng)的基于FRFT的方法性能嚴(yán)重下降,特征遭到嚴(yán)重破壞,而非線性變換可以很好地抑制α穩(wěn)定分布噪聲。

3.2 提取特征信息并通過SVM對接收信號的特征量進(jìn)行分類識別

首先,在東海某海域仿真信道和穩(wěn)定分布噪聲條件下,MSNR在-15～10 dB范圍內(nèi),選取12個不同參數(shù)的仿真信號,包括FSK、PSK、OFDM和LFM信號,其中LFM的信號初始頻率為13 kHz,帶寬為2 kHz、3 kHz、4 kHz、5 kHz、6 kHz、7 kHz,信號長度均為50 ms,每間隔5 dB對每種信號進(jìn)行50次蒙特卡羅仿真實驗,計算特征量歸一化方差比的值并取平均,得到不同類型信號在不同MSNR條件下的方差比。

圖6 不同類型信號在不同MSNR條件下的方差比Fig.6 Variance ratio of different types of signals under different MSNR corditions

從仿真數(shù)據(jù)可知,利用方差比作為特征量對LFM信號進(jìn)行識別是可行的。

根據(jù)Argo全球海洋信息數(shù)據(jù)集,得到東海某海域和南海某海域的聲剖面信息。通過Bellhop仿真得到相應(yīng)的水聲信號仿真波形。對仿真接收信號計算識別特征量,并對SVM分類器進(jìn)行訓(xùn)練。

訓(xùn)練完成后,對LFM的識別率進(jìn)行測試。在MSNR -15 ～10 dB范圍內(nèi),每間隔1 dB產(chǎn)生500組隨機(jī)的信號,包括PSK、FSK、OFDM及帶寬為2～7 kHz的LFM信號。將仿真得到的信道沖激響應(yīng)與信號卷積,并選擇噪聲條件為加性標(biāo)準(zhǔn)SαS分布噪聲,特征指數(shù)為α=1.5,得到待測信號。

將待測信號送入訓(xùn)練好的SVM分類器中進(jìn)行LFM信號的識別,得到LFM的識別率隨MSNR變化的曲線。此外,將文獻(xiàn)[17]所述的基于Robust STFT的方法在信號長度為50 ms情況下進(jìn)行仿真,結(jié)果如圖7(a)所示,仿真結(jié)果表明,本文提出的方法有了較大改進(jìn)。選取南海某海域的仿真信道,在同樣條件下得到如圖7(b)所示的結(jié)果。

圖7 在不同MSNR條件下對LFM信號的識別率曲線Fig.7 Recognition rate curve of LFM signal under different MSNR conditions

由圖可知,與文獻(xiàn)[17]所述方法相比,本文提出的方法有了較大改進(jìn)。

為了驗證在不同α值的穩(wěn)定噪聲條件下識別率的變化情況,選取了0.8、1.2、1.5、1.8、2作為α值進(jìn)行仿真實驗。實驗采用上述東海某海域信道條件,對每個α值,MSNR在-15～10 dB范圍內(nèi),每間隔1 dB產(chǎn)生100組隨機(jī)的帶寬為2～7 kHz的LFM信號進(jìn)行仿真實驗,結(jié)果如圖8所示。

圖8 在不同α值條件下對LFM信號的識別率曲線Fig.8 Recognition rate curve of LFM signal under different α

仿真實驗說明,本文所述方法在不同的α值環(huán)境中性能較穩(wěn)定,且均有較高的識別率。

4 結(jié) 論

本文提出了基于非線性映射DFRFT的LFM信號識別方法。首先,通過非線性變換對α穩(wěn)定分布噪聲進(jìn)行抑制,然后進(jìn)行FRFT并通過提取的歸一化方差比特征量識別出LFM信號。仿真表明,該方法在MSNR和α穩(wěn)定分布噪聲條件下能夠有效識別LFM信號。

參考文獻(xiàn)：

[1] 梁國龍,生雪莉. 基于WVD-HT的寬帶調(diào)頻信號檢測技術(shù)研究[J]. 電子學(xué)報, 2011, 32(12):1941-1944.

LIANG G L, SHENG X L. Research on wideband FM signal detection technology based on WVD-HT[J]. Atca Electronica Sinica, 2011, 32(12): 1941-1944.

[2] WANG H W, FAN X Y, CHEN Y, et al. Wigner-hough transform based on slice’s entropy and its application to multi-LFM signal detection[J]. Journal of Systems Engineering and Electronics, 2017, 28(4):634-642.

[3] DING Y F, SUN L, ZHANG H Y, et al. A multi-component LFM signal parameters estimation method using STFT and Zoom-FRFT[C]∥Proc.of the 8th IEEE International Conference on Communication Software and Networks,2016:112-117.

[4] DUAN Y, WANG J Z, SU S Y, et al. Detection of LFM signals in low SNR based on STFT and wavelet denoising[C]∥Proc.of the International Conference on Audio, Language and Image Processing, 2014:921-925.

[5] MA X, LIU D, SHAN Y. Intra-pulse modulation recognition using short-time ramanujan Fourier transform spectrogram[J]. Eurasip Journal on Advances in Signal Processing, 2017,2017(1):42.

[6] 崔化超,周勝增,蔣小勇. 水聲脈沖信號實時檢測估計軟件的設(shè)計與實現(xiàn)[J]. 聲學(xué)技術(shù), 2015, 34(4): 74-77.

CUI H C, ZHOU S Z, JIANG X Y. Design and realization of real-time detection and estimation software for underwater acoustic pulse signal[J].Technical Acoustics, 2015,34(4):74-77.

[7] 黃冠欽,楊坤德,馬遠(yuǎn)良. 微弱信號的自適應(yīng)相干累積檢測[J].聲學(xué)技術(shù), 2011,30(3): 198-200.

HUANG G Q, YANG K D, MA Y L. Adaptive coherent accumulation detection of weak signal[J]. Technical Acoustics, 2011, 30(3): 198-200.

[8] 王曉燕,方世良,朱志峰. 基于ST-FRFT的非合作水聲脈沖信號檢測方法[J]. 信號處理, 2011, 21(8): 1271-1278.

WANG X Y, FANG S L, ZHU Z F. Non-cooperative underwater acoustic pulse signal detection method based on ST-FRFT[J]. Journal of Signal Processing, 2011, 21(8): 1271-1278.

[9] CHEN Y L, GUO L X, GONG Z X. The concise fractional Fourier transform and its application in detection and parameter estimation of the linear frequency-modulated signal[J].Chinese Journal of Acustica, 2017(1): 70-86.

[10] ZHANG Z C, WANG H Y, LIU Z G, et al. An new LFM detection algorithm based on ADTR-FrFT for active sonar system[C]∥Proc.of the IEEE International Conference on Signal Processing, Communications and Computing, 2017:1-5.

[11] SHI H, XIONG J, ZHOU C, et al. A new recognition and classification algorithm of underwater acoustic signals based on multi-domain features combination[C]∥Proc.of the Ocean Acoustics, 2016: 1-7.

[12] HOW J, HATZINAKOS D. Impulsive noise modeling with stable distributions in fading environments[C]∥Proc.of the Statistical Signal and Array Processing, 1996: 140-140.

[13] LAGUNA S G, LOPEZ G M. On the use of alpha-stable distributions in noise modeling for PLC[J]. IEEE Trans.on Power Delivery, 2015, 30(4):1863-1870.

[14] CHEN Z, GENG X, YIN F. A harmonic suppression method based on fractional lower order statistics for power system[J]. IEEE Trans.on Industrial Electronics,2016,63(6):3745-3755.

[15] MAHMOOD A, CHITRE M, Armand M A. Maximum-likelihood detection performance of uncoded OFDM in impulsive noise[C]∥Proc.of the IEEE Global Communications Conference, 2013: 3516-3521

[16] QIU T S, ZHANG X. Statistical signal processing: non-gaussion signal processing and applications[M].Beijing: Electronic Industry Press, 2004: 169-172.

[17] 陳文武,蔡征宇,陳如山,等.脈沖噪聲下基于Robust STFT的信號檢測與參數(shù)估計[J].南京理工大學(xué)學(xué)報,2012,36(2):328-331.

CHEN W W, CAI Z Y, CHEN R S, et al. Signal detection and parameter estimation based on robust STFT under impulse noise[J]. Journal of Nanjing University of Science and Technology, 2012, 36(2): 328-331.

[18] 金艷,高舵,姬紅兵.α穩(wěn)定分布噪聲下基于穩(wěn)健S變換的LFM信號參數(shù)估計[J].系統(tǒng)工程與電子技術(shù),2017, 39(4):693-699.

JIN Y, GAO D, JI H B. Parameter estimation of LFM signal based on robust S transform inαstable distribution noise[J]. Systems Engineering and Electronics, 2017, 39(4):693-699.

[19] 金艷,胡碧昕,姬紅兵.α穩(wěn)定分布噪聲下基于最優(yōu)L-柯西加權(quán)的LFM信號參數(shù)估計[J].系統(tǒng)工程與電子技術(shù),2016,38(7): 1488-1495.

JIN Y, HU B X, JI H B. Parameter estimation of LFM signal based on optimal L-Cauchy weighted method inαstable distribution noise[J]. Systems Engineering and Electronics, 2016, 38(7):1488-1495.

[20] 張俊林.α穩(wěn)定分布噪聲下通信信號調(diào)制分析技術(shù)研究[D].鄭州: 信息工程大學(xué), 2005.

ZHENG J L. Research on communication signal modulation analysis technology under alpha stable distributed noise[D]. Zhengzhou:Information Engineering University, 2005.

[21] OZAKTAS H M, ARIKAN O, KUTAY M A, et al. Digital computation of the fractional Fourier transform[J]. IEEE Trans.on Signal Processing, 1996, 44(9): 2141-2150.

[22] 陶然,鄧兵,王越. 分?jǐn)?shù)階傅里葉變換及其應(yīng)用[M].北京: 清華大學(xué)出版社, 2009: 13-25.

TAO R, DENG B, WANG Y. Fractional Fourier transform and its application[M]. Beijing:Tsinghua University Press, 2009: 13-25.

[23] VAPNIK V. The nature of statistical learning theory[M].New York: Springer, 1995:142-165.

[24] ZHAO H, CHEN H, NGUYEN T. Stratified over-sampling bagging method for random forests on imblanced data[M]. Switzerland: Springer International Publishing, 2016:46-60.

[25] VLUYMANS S, TARRAGOD S, SAEYS Y. Fuzzy rough classifiers for class imbalanced multi-instance data[J].Pattern Recognition, 2016, 53(C):36-45.

[26] LIU ALEXANDER, GHOSH J, MARTIN C. Generative oversampling for mining imbalanced datasets[C]∥Proc.of the International Conference on Data Mining, 2007:66-72.

[27] China Argo Real-time Data Center. The International Argo Project[EB/OL].[2017-8-08].http:∥www.argo.org.cn/data/argo.php/.