UWB雷達(dá)雜波物理建模與仿真分析

楊利民 許志勇 蘇衛(wèi)民 顧 紅

(南京理工大學(xué)電子工程與光電技術(shù)學(xué)院,江蘇 南京210049)

1.引 言

UWB雷達(dá)是指相對帶寬不小于 25%的雷達(dá)[1],相對帶寬定義為2(fH-fL)/(fH+fL),通常情況下UWB雷達(dá)具有很高的距離分辨力。雜波是指雷達(dá)接收機(jī)接收的不感興趣的回波[2],通常用σ0來描述。雖然UWB雷達(dá)分辨率高,具有良好的雜波抑制能力,但在雷達(dá)設(shè)計(jì)的分析階段,以及在優(yōu)化階段,當(dāng)無法獲取測試數(shù)據(jù)時(shí),雜波的建模和仿真依然具有重要的指導(dǎo)意義。

由于雜波具有隨機(jī)性,通常利用統(tǒng)計(jì)特性進(jìn)行描述[3-4]。采用統(tǒng)計(jì)模型的另一原因是傳統(tǒng)的目標(biāo)檢測理論基于雜波的統(tǒng)計(jì)特性來判斷是否存在目標(biāo)。但像UWB雷達(dá)這樣的高距離分辨率雷達(dá),該模型不再服從如Rayleigh分布之類的典型短拖尾[4]。此外,雜波統(tǒng)計(jì)模型物理特性不明顯,建立的雜波模型不能直接反映測量區(qū)域的地形變化特征。另常用雜波建模方法是反映地表面變化特征的物理建模。文獻(xiàn)[5]中視測量區(qū)由不同的平板組合而成(表面由許多小且方向不同的小平板組合),以簡化復(fù)雜地貌。每塊小平板的反射方向圖可視為等尺寸的線性天線,回波為入射波在平板上鏡面反射的一部分(法向方向反射最大)。

對于UWB雷達(dá),文獻(xiàn)[6]中描述了雜波物理模型,但是沒有考慮頻率的影響。目前采用平板模型對UWB雷達(dá)雜波建模的研究很少,雖然文獻(xiàn)[2]和[7]中描述了平板模型,但都是基于窄帶信號(hào)的。在我們的前期工作中[8]研究了寬帶雷達(dá)雜波平板模型,但該模型比較粗簡,未考慮諸多細(xì)節(jié)因素,比如入射波和反射波的極化、地表電磁特性(相對磁導(dǎo)率、介電常數(shù))、地面凹凸不平和植被覆蓋等,對UWB雷達(dá)雜波模型的影響。采用平板物理模型和“頻率分割子帶合成”法對UWB雷達(dá)地雜波進(jìn)行建模和仿真,并且給出了測量區(qū)域的雜波統(tǒng)計(jì)模型。其基本思路是首先根據(jù)地形起伏把測量區(qū)域劃分成若干個(gè)稱為平板的小區(qū)域;根據(jù)窄帶滿足的條件把UWB劃分成一系列子帶(即窄帶)。其次利用成熟的窄帶雷達(dá)雜波模型[9-13]建立子帶的平板雜波模型。最后,對平板內(nèi)的分辨單元的各個(gè)子帶貢獻(xiàn)利用如逆傅里葉變換(IFT)的方法進(jìn)行子帶合成,但子帶合成不是簡單的各個(gè)子帶的σ0相加,由于彼此間存在一定的相位關(guān)系,所以子帶合成事實(shí)上是矢量合成。以模擬各種UWB雷達(dá)雜波的細(xì)節(jié)部分[7],在子帶合成時(shí)還進(jìn)行了細(xì)節(jié)的補(bǔ)償。本文目的是通過仿真驗(yàn)證在包含不同的地貌和參數(shù)變化在內(nèi)的條件下的UWB雷達(dá)雜波的相關(guān)規(guī)律,為雷達(dá)分析和優(yōu)化提供參考,而不涉及具體的地貌匹配條件。

2.基本思路及雜波物理建模

2.1 基本思路

在討論雜波物理建模之前先介紹一下文中的基本思路:首先根據(jù)地形變化對雜波測量區(qū)域劃分成若干個(gè)小區(qū)域(稱為平板),再將UWB雷達(dá)信號(hào)分割成一系列子帶信號(hào),因而利用已有的成熟窄帶雜波模型通過子帶合成方式獲得 UWB雷達(dá)雜波模型,文中稱之為“頻率分割子帶合成”法[8]。在每個(gè)子帶內(nèi)采用平板物理模型得到低距離分辨率的窄帶雜波雷達(dá)截面積(RCS),然后通過加權(quán)IFT進(jìn)行子帶合成獲得UWB雷達(dá)雜波的數(shù)學(xué)模型。對 UWB雷達(dá)分辨單元面積歸一后獲得雜波散射系數(shù)(該徑向中每個(gè)樣點(diǎn)即為相應(yīng)UWB雷達(dá)分辨單元的雜波散射系數(shù)),帶寬越寬,它愈能反映地面散射細(xì)節(jié),即地面的一維距離像。最后利用最大似然參數(shù)估計(jì)(MLE)法估計(jì)其統(tǒng)計(jì)參數(shù),且利用柯爾莫諾夫-斯米爾諾夫檢驗(yàn)(Kolmogorov-Smirnov test)法優(yōu)選出雜波的統(tǒng)計(jì)分布。基本思路見圖1流程框圖。

圖1 全文基本思路框圖

2.2 平板模型

傳統(tǒng)的平板模型沒有考慮復(fù)雜地形,無論地表的形狀如何變化,都認(rèn)為是簡單的平板以不同的方向組合而成。對于平原,該模型具有比較好的適應(yīng),但對于人工建筑物或樹木等,采用傳統(tǒng)的平板模型顯然不再合適。廣義平板模型除了傳統(tǒng)的平板模型外,還包括柱面、二面角、邊緣和曲面等常見的幾何形狀。例如,利用柱面模型能較好描述樹木,二面角描述建筑物墻面和地面構(gòu)成的二面關(guān)系等。幾種常見廣義平板模型如圖2。

我們在前期工作[8]中通過引入廣義平板模型,構(gòu)建了UWB雷達(dá)雜波建模的路線框架,但未考慮諸多細(xì)節(jié)因素對UWB雷達(dá)雜波模型的影響。在此基礎(chǔ)上,文中對其中最簡單的平板模型進(jìn)行重點(diǎn)研究,加入一些實(shí)際細(xì)節(jié)因素的影響,且對原雜波模型進(jìn)行修正,通過仿真預(yù)測UWB雷達(dá)雜波的一些規(guī)律。關(guān)于廣義平板中其他更復(fù)雜的模型的相關(guān)研究將在后續(xù)工作中開展。傳統(tǒng)平板電磁波傳播示意圖見圖 3,其中入射角和反射角分別為θi和θs,且設(shè)θi=θs=θ,φ(即式(2)中的 φs)為反射波在 xoy平面上的投影與x軸之間的夾角。窄帶雷達(dá)雜波反射截面積的數(shù)學(xué)模型如下[7]

式中:

式中,R⊥(θ)為垂直方向上的菲涅耳散射系數(shù)

式中,μr和εr分別為地表面的相對磁導(dǎo)率和介電常數(shù)。如此可知式(1)不但考慮了地面的起伏,而且考慮了電磁波入射和反射的極性以及地表的電磁特性(介電常數(shù))。通常情況下,起伏的地表上同時(shí)考慮植被覆蓋更能反映實(shí)際地形情況。文中假設(shè)植被服從復(fù)高斯分布或其它弱拖尾的統(tǒng)計(jì)分布(例如雜草、農(nóng)作物等)。從而子帶的雜波平板散射模型[8]為

式中:r表示平板上距離門和雷達(dá)之間的距離(見式(6));ρ表示植被(例如雜草、農(nóng)作物等)引起的復(fù)高斯或其他弱拖尾的統(tǒng)計(jì)值;p(fi)和q(fi)分別為調(diào)節(jié)高低不平的地面和植被對雜波數(shù)據(jù)的貢獻(xiàn)大小的權(quán)值,并且假設(shè)式(5)中 γcpq的相位φ(f i)為平板相頻特性,假設(shè)服從[0,2π]的均勻分布(具體統(tǒng)計(jì)分布規(guī)律可通過雜波實(shí)測數(shù)據(jù)加以修正)。

式(5)中由于γcpq中相位和ρ都為隨機(jī)項(xiàng),所以(,r)也具有統(tǒng)計(jì)特性。對于UWB雷達(dá),具有高分辨能力,強(qiáng)散射體貢獻(xiàn)未被平均掉,導(dǎo)致雜波統(tǒng)計(jì)上表現(xiàn)出長拖尾特性,雜波仿真結(jié)果(第3節(jié))證實(shí)了這一特征。但是對于測量區(qū)域茂密的植被,雜波主要由q(f i)ρ決定。綜上可看出,式(5)較文獻(xiàn)[8]模型更能反映地形變化以及植被覆蓋的一般情況。

2.3 子帶合成

一系列子帶雜波的貢獻(xiàn)通過子帶合成的方法實(shí)現(xiàn)超寬帶雷達(dá)雜波建模。設(shè)超寬帶雷達(dá)帶寬為BW,從而距離分辨率為δR=c/(2BW)(c為光速)。對位于rk(k=1,2,…),長為L的平板,覆蓋的距離門個(gè)數(shù)M由L/δR取整獲得,則M個(gè)距離門和雷達(dá)之間的距離為

則中心位于r k(k=1,2,…)處的平板的第m個(gè)分辨單元的雜波模型(雷達(dá)目標(biāo)截面積)為

上式事實(shí)上是逆傅里葉變換,其中I為子帶個(gè)數(shù);f i為第i個(gè)子窄帶的中心頻率;rk為平板中心斜距;W(f i,rk,m)為細(xì)節(jié)補(bǔ)償因子,實(shí)際上反映UWB雷達(dá)雜波RCS的頻譜結(jié)構(gòu),該結(jié)構(gòu)受諸多因素(如多徑、頻率和散射點(diǎn)分布等)的影響,通常是非白的,即頻譜的能量可能聚集在較少的幾個(gè)子帶內(nèi),即從定性上分析,W(fi,rk,m)在統(tǒng)計(jì)上應(yīng)具有長拖尾特征。文中暫時(shí)采用典型的 Lognormal長拖尾分布經(jīng)過模歸一化后進(jìn)行FFT得到W(fi,rk,m),用以補(bǔ)償子帶合成前子帶間雜波細(xì)節(jié)的光滑化效應(yīng)[7]。但在實(shí)際中具體采用何種長拖尾分布模型,還需進(jìn)一步地研究,并通過實(shí)測數(shù)據(jù)加以修正。

由式(7)計(jì)算σ0[14]為

式中,A為分辨單元的面積。

3.實(shí)驗(yàn)及結(jié)果分析

仿真條件為:測量區(qū)為長15 m、寬1 m的地面;式(6)中r k=5 km;設(shè)圖 3 中 θi=θs=θ=15°,φ=180°;頻率分割的子帶相對帶寬為1%;入射波和反射波均為水平極化(記為H-H);式(7)中W(f i,r)服從(0,10)的Lognormal分布。

文中UWB信號(hào)采用具有代表性的線性調(diào)頻(LFM)信號(hào),由于篇幅所限,在此未對其他 UWB信號(hào)如時(shí)域沖擊波形(impulse波形)等一一進(jìn)行討論。仿真分析的第一種情況是地貌分別為裸地(p(f)=0.8,q(f)=0.1)和植被覆蓋(p(f)=0.1,q(f)=0.8)且中心頻率為6 GHz,相對帶寬為25%(即5.25～6.75 GHz)時(shí)的雜波統(tǒng)計(jì)特性比較。另一種情況是中心頻率分別為6 GHz(5.25～6.75 GHz)和10 GHz(9.25～10.75 GHz),且地貌為裸地(p(f)=0.8,q(f)=0.1)時(shí)對雜波統(tǒng)計(jì)特性的影響。地貌和中心頻率變化時(shí)第一組仿真雜波數(shù)據(jù)分別如圖4和圖5所示(共進(jìn)行了15次統(tǒng)計(jì)獨(dú)立重復(fù)仿真)

圖4給出地貌分別是裸地((a)圖)和植被覆蓋((b)圖)時(shí)的雜波散射系數(shù)與平板斜距的關(guān)系,從圖中可以看出在植被條件下雜波和裸地相比起伏要小些,是由于植被的存在,削弱了強(qiáng)散射體的貢獻(xiàn)。圖5給出了中心頻率分別為10 GHz((a)圖)和6 GHz((b)圖),絕對帶寬均為1.5 GHz的雜波散射系數(shù)與平板斜距的關(guān)系,從圖中可知在絕對帶寬和測量區(qū)域尺寸相同的情況下,中心頻率低的雜波散射系數(shù)總體要大于中心頻率高的,這是因?yàn)樵诮^對帶寬相同時(shí),中心頻率低,則對應(yīng)的相對帶寬更寬,強(qiáng)散射的反射特性將更加凸顯,同時(shí)由于低頻段平板散射方向圖的主瓣更寬,旁瓣更高,所以中心頻率低的雜波散射系數(shù)大。

利用MLE估計(jì)的參數(shù)[15]和統(tǒng)計(jì)分布擬合情況(柯爾莫諾夫-斯米爾諾夫(K-S)檢驗(yàn)統(tǒng)計(jì)量D)[16]如表1和表2所示。

表1 不同地貌條件下統(tǒng)計(jì)分布參數(shù)估計(jì)和統(tǒng)計(jì)量

表2 不同頻段條件下統(tǒng)計(jì)分布參數(shù)估計(jì)和統(tǒng)計(jì)量

不同的地貌變化和波段變化時(shí)雜波仿真數(shù)據(jù)直方圖及其擬合曲線分別如圖6和圖7所示(15次統(tǒng)計(jì)獨(dú)立重復(fù)仿真)。

從表1看出當(dāng)?shù)孛矠槁愕貢r(shí)雜波服從Gen.Gamma分布,而當(dāng)?shù)孛灿兄脖桓采w時(shí)則服從Gamma分布。對于表2,當(dāng)絕對帶寬為1.5 GHz,中心頻率為 6 GHz(相對帶寬為 25%)時(shí)雜波服從Weibull分布,而中心頻率為10 GHz(相對帶寬為15%)時(shí)雜波服從Dagum分布。圖6表示其他條件相同,而地貌變化(其上圖表示裸地、下圖表示植被覆蓋)時(shí)雜波散射系數(shù)的直方圖及其統(tǒng)計(jì)分布擬合曲線;圖7表示其他條件相同,而中心頻率變化(其上圖表示6 GHz、下圖表示10 GHz)時(shí)雜波散射系數(shù)的直方圖及其統(tǒng)計(jì)分布擬合曲線。從圖6和圖7均可發(fā)現(xiàn),雜波散射系數(shù)PDF擬合曲線均具有“低重心、重拖尾”的特點(diǎn)。表現(xiàn)出重拖尾是由于UWB具有高距離分辨能力,強(qiáng)散射體可能超過一個(gè)或多個(gè)距離門,表現(xiàn)出強(qiáng)反射特性,不像窄帶信號(hào)由于存在大量的弱散射體可能抵消強(qiáng)散射體的貢獻(xiàn),所以對于UWB雷達(dá)雜波PDF曲線均具有重拖尾;低重心是指在σ0近0處具有最大的概率,和窄帶信號(hào)相比,雜波統(tǒng)計(jì)分布重心向低散射系數(shù)方向移動(dòng),其原因是因?yàn)楦呔嚯x分辨率,極強(qiáng)散射體可能遮蔽一定數(shù)量的且比其弱的散射體,形成極弱散射或0散射(完全遮蔽),所以PDF曲線重心同時(shí)往0散射系數(shù)方向移動(dòng)。由于UWB信號(hào)具有高距離分辨,在距離單元中的未被遮蔽的弱散射體數(shù)目也很少,表現(xiàn)出重心對應(yīng)的概率較窄帶信號(hào)很小。由于雜波的“低重心、長拖尾”的特性,和窄帶雷達(dá)相比,UWB雷達(dá)具有更大的動(dòng)態(tài)范圍,從而易造成雷達(dá)接收機(jī)的飽和甚至過載,導(dǎo)致系統(tǒng)的損壞。同時(shí)從圖7看出,在絕對帶寬相同,而中心頻率變小時(shí),由于導(dǎo)致相對帶寬的變寬,其統(tǒng)計(jì)分布的拖尾比中心頻率大的統(tǒng)計(jì)分布拖尾更長。

圖6 不同地貌條件下雜波統(tǒng)計(jì)分布模型擬合

圖7 不同頻段條件下雜波統(tǒng)計(jì)分布模型擬合

通過仿真結(jié)果分析進(jìn)一步發(fā)現(xiàn),在其他條件相同,而地貌變化時(shí),雜波統(tǒng)計(jì)分布也隨著變化,如圖6(上、下圖)。根據(jù)PDF曲線擬合得知當(dāng)?shù)孛矠槁愕貢r(shí)服從Gen.Gamma分布,但當(dāng)?shù)孛灿兄脖桓采w時(shí),則服從Gamma分布。同樣,對于其他條件相同而只有中心頻率變化時(shí),雜波統(tǒng)計(jì)分布也發(fā)生變化,見圖7(上、下圖),在中心頻率為6 GHz時(shí)服從Weibull分布,但當(dāng)中心頻率變?yōu)?0 GHz時(shí)則服從Dagum分布。由此發(fā)現(xiàn)對于UWB雷達(dá),無法用單一的統(tǒng)計(jì)模型對其建模,這和窄帶雷達(dá)雜波模型(通常地雜波服從具有弱拖尾的Rayleigh分布)具有普適性有很大不同。因此,如果依然對 UWB雜波采用統(tǒng)計(jì)建模則帶來很大的困難,這也是文中采用物理建模的原因之一。此外,由于傳統(tǒng)的目標(biāo)檢測是采用萊曼-皮爾遜(Neyman-Pearson)準(zhǔn)則,該準(zhǔn)則是基于統(tǒng)計(jì)模型的,所以雜波統(tǒng)計(jì)模型隨參數(shù)變化而變化的特性為UWB雷達(dá)目標(biāo)檢測帶來挑戰(zhàn)。

圖8和圖9分別表示地貌和波段的變化時(shí)統(tǒng)計(jì)分布曲線比較。圖8是地貌分別為裸地和植被時(shí)PDF曲線局部放大圖(由于植被條件下拖尾比裸地的要長得多,局部放大以便更清楚看出兩者PDF曲線的關(guān)系),根據(jù)計(jì)算得出擬合曲線在聯(lián)合概率密度函數(shù)CDF值為0.9處對應(yīng)的散射系數(shù)為1.3408×104(裸地)和5.6899×104(植被)。圖9是中心頻率分別為6 GHz和 10 GHz的 PDF曲線,在CDF值為0.9時(shí),計(jì)算出散射系數(shù)分別1.259×104(6 GHz)和1.3618×103(10 GHz)。由圖9可知當(dāng)相同的絕對帶寬但中心頻率變化時(shí),相同CDF值對應(yīng)的散射系數(shù)位置也發(fā)生移動(dòng),對于較小的中心頻率,當(dāng)CDF值為0.9時(shí)對應(yīng)的雜波散射系數(shù)位置更遠(yuǎn)離0散射系數(shù)的位置,這是由于較小的中心頻率在相同的絕對帶寬條件下有較大的相對帶寬,從而具有更長的拖尾,且PDF曲線重心沿PDF軸下移所致。但由此會(huì)導(dǎo)致虛警概率增大,為雷達(dá)目標(biāo)檢測、辨識(shí)和數(shù)據(jù)處理帶來困難,更有可能在強(qiáng)散射體占優(yōu)的分辨單元內(nèi)由于很小的信雜比(SCR)而導(dǎo)致目標(biāo)被“埋沒”,造成距離盲區(qū)。

4.結(jié) 論

文中采用物理平板模型和“頻率分割子帶合成”法對UWB雷達(dá)雜波進(jìn)行建模與仿真。首先根據(jù)地形把地面劃成若干個(gè)小區(qū)域,同時(shí)把超寬帶分割成一系列子帶;然后分別計(jì)算子帶的雜波貢獻(xiàn);最后利用子帶合成得到超寬帶雷達(dá)雜波的模型。通過統(tǒng)計(jì)模型建模只能適應(yīng)具體地形、季節(jié)和某時(shí)間的雜波情況;文中構(gòu)建的UWB雷達(dá)雜波模型不但反映了地面凹凸不平,植被覆蓋的現(xiàn)象,而且考慮了地表面的電磁特性(如介電常數(shù))和電磁波入射與反射的極化特性,具有高保真性。仿真結(jié)果表明:對于UWB雷達(dá)雜波,在其他條件相同的情況下,地貌或者波段(或中心頻率)的變化,雜波的統(tǒng)計(jì)分布也跟隨變化,因此,對于UWB雷達(dá)雜波不能用固定的統(tǒng)計(jì)分布對其建模,同時(shí)由于長拖尾的雜波容易導(dǎo)致虛警概率增大,這給UWB雷達(dá)目標(biāo)檢測帶來挑戰(zhàn);UWB雷達(dá)雜波均具有長拖尾的統(tǒng)計(jì)分布特征,在接近0散射系數(shù)σ0處概率密度最大,即“低重心,長拖尾”。且隨著相對帶寬的增大,拖尾加重,則UWB雷達(dá)具有更大的動(dòng)態(tài)范圍,從而容易出現(xiàn)接收機(jī)的飽和甚至過載現(xiàn)象。經(jīng)仿真分析可知:仿真結(jié)果與UWB雷達(dá)雜波定性分析的規(guī)律基本符合。后續(xù)工作中將進(jìn)一步研究針對其他UWB信號(hào)(如時(shí)域impulse波形等)的雜波建模和仿真。

[1] JAMES D.Taylor.Introduction to ultra-wideband radar system[M].CRC Press,Inc.,1995:1-2.

[2] SKOLNIK M I.Radar handbook[M].3nd ed.McGraw-Hill,Inc.,2008:14.

[3] ANASTASSOPOULOS V,LAMPROPOULOS G A,DROSOPOULOS A,et al.High resolution radar clutter statistics[J].IEEE Trans on AES,1999,35(1):43-60.

[4] GRECO M S,GINI F.Statistical analysis of high-resolution SAR ground clutter data[J].IEEE Trans on GRS,2007,45(3):556-575.

[5] GARCIA-FERNANDEZ A F,GRAJAL J,YESTEOJEDA O A.Analysis of ISAR images of a helicopter by a facet model[C]//IEEE International Conference on Radar.USA,2008:32-37.

[6] JACHSON JA,MOSES R L.Clutter model for VHF SA R imagery[J].Proceedings of the SPIE,2004,5427.271-282.

[7] RUCK G T,BARRICK D E,STUART W D,et al.Radar Cross Section Handbook Volume 2[M].New York:Plenum Press,1970:700-702.

[8] XU Zhiyong.Subband Systhesis mehtod for simulation of wideband radar clutter[C]//IET International Radar Conference.Guilin,China,April,2009:20-22.

[9] BILLINGSLEY B,FARINA A,GINI F,et al.Statistical analysis of measureed radar ground clutter data[J].IEEE Trans on AES,1999 35(2):579-592.

[10] MA Xiaoyan,FANG Xueli,CHEN Song.Radar clutter recognition based on feature extraction byαtruncation set[C]//Proc.IEEE,IRC'.Beijing,China,2001:436-439.

[11] BILLINGSLEY J B.Low angle radar land clutter:measurements and empirical models[M].USA,William Andrew Publishing,Inc.,2001:143-242.

[12] 張長隆,王 展,周良柱.線性調(diào)頻脈沖壓縮雷達(dá)雜波統(tǒng)計(jì)模型分析[J].電波科學(xué)學(xué)報(bào),2004,19(2):240-244.

ZHANG Changlong,WANG Zhan,ZOU Liangzhu.Clutter statistical analysis for linear f requency modulated pulse compression radar[J].Chinese Journal of Radio Science,2004,19(2):240-244.(in Chinese)

[13] 楊 博,周一宇,黃知濤.空基運(yùn)動(dòng)平臺(tái)照射源的無源雷達(dá)雜波建模[J].電波科學(xué)學(xué)報(bào),2009,24(6):1099-1103.

YANG Bo,ZHOU Yiyu,HUANG Zhitao.Modeling of the ground clutter on passive radar based on airborn illuminator[J].Chinese Journal of Radio Science,2009,24(6):1099-1103.(in Chinese)

[14] SKOLNIK M I.Introduction to Radar Systems(3rd edition)[M].Beijing,China,Publishing House of E-lectronics Industry,2006:403-404.

[15] LIANG Jing,LIANG Qilian,SAMN S W.Foliage clutter modeling using the UWB radar[C]//IEEE International Conference on Communications.Beijing,China,19-23 M ay 2008:1937-1941.

[16] STEPHENS M A.EDF Statistics for goodness of fit and some comparisons[J].Journal of the American Statistical Association 69,1974:730-737.