雷 雨，趙丹寧

（1．中國(guó)科學(xué)院國(guó)家授時(shí)中心，陜西 西安710600；2．中國(guó)科學(xué)院時(shí)間頻率基準(zhǔn)重點(diǎn)實(shí)驗(yàn)室，陜西 西安710600；3．中國(guó)科學(xué)院精密導(dǎo)航定位與定時(shí)技術(shù)重點(diǎn)實(shí)驗(yàn)室，陜西 西安710600；4．中國(guó)科學(xué)院，北京100049）

0 引 言

衛(wèi)星鐘差的精度直接影響著全球衛(wèi)星導(dǎo)航系統(tǒng)（GNSS）的導(dǎo)航定位結(jié)果，主要體現(xiàn)在以下兩個(gè)方面:首先，衛(wèi)星鐘差參數(shù)通過(guò)導(dǎo)航電文發(fā)播，任意時(shí)刻的鐘差是通過(guò)外推得到的，其預(yù)報(bào)精度直接決定了導(dǎo)航定位的精度［1］；其次，在實(shí)時(shí)精密單點(diǎn)定位（PPP）中，為了達(dá)到厘米級(jí)的定位精度，需采用預(yù)報(bào)鐘差參與解算［2］。因此，提高衛(wèi)星鐘差的預(yù)報(bào)精度具有重要意義。

較常用的鐘差預(yù)報(bào)模型為多項(xiàng)式模型和灰色模型，許多學(xué)者對(duì)其進(jìn)行了研究，并取得了顯著的效果［3－4］。但這兩種模型都存在一定的局限性:多項(xiàng)式模型的預(yù)報(bào)精度會(huì)隨著預(yù)報(bào)時(shí)間的延長(zhǎng)而降低［4］，此外，該方法易受多項(xiàng)式擬合階數(shù)等人為因素的影響；而灰色模型對(duì)鐘差數(shù)據(jù)類(lèi)型有較強(qiáng)的依賴(lài)性，要求鐘差數(shù)據(jù)平滑且呈指數(shù)規(guī)律變化，限制了其應(yīng)用范圍［5］。近年來(lái)，國(guó)內(nèi)外學(xué)者又發(fā)展了多種鐘差預(yù)報(bào)模型［6－7］，其中神經(jīng)網(wǎng)絡(luò)（NN）由于具有通過(guò)學(xué)習(xí)訓(xùn)練以任意精度逼近任意連續(xù)函數(shù)的能力，在鐘差預(yù)報(bào)中取得了良好的效果［8－10］。李孝輝等首次將誤差反向傳播（BP）神經(jīng)網(wǎng)絡(luò)應(yīng)用于原子鐘鐘差預(yù)報(bào)中，研究表明神經(jīng)網(wǎng)絡(luò)比AR模型的預(yù)報(bào)精度高，且鐘差異常值對(duì)BP神經(jīng)網(wǎng)絡(luò)模型預(yù)報(bào)性能的影響較??；郭承軍等將徑向基函數(shù)（RBF）神經(jīng)網(wǎng)絡(luò)應(yīng)用于衛(wèi)星鐘差預(yù)報(bào)中，研究表明:神經(jīng)網(wǎng)絡(luò)比灰色模型具有更好的預(yù)報(bào)精度和穩(wěn)定性；王威等采用RBF神經(jīng)網(wǎng)絡(luò)的變型—廣義回歸神經(jīng)網(wǎng)絡(luò)（GRNN）對(duì)衛(wèi)星鐘差進(jìn)行建模并預(yù)報(bào)，從輸入維數(shù)。光滑因子等幾個(gè)方面對(duì)預(yù)報(bào)性能的影響進(jìn)行了較為系統(tǒng)的研究，結(jié)果表明GRNN預(yù)報(bào)模型優(yōu)于傳統(tǒng)的二次多項(xiàng)式模型。應(yīng)用RBF神經(jīng)網(wǎng)絡(luò)進(jìn)行衛(wèi)星鐘差預(yù)報(bào)，由于神經(jīng)網(wǎng)絡(luò)易學(xué)習(xí)過(guò)度，從而影響網(wǎng)絡(luò)的泛化（預(yù)報(bào)）能力，因而采用交叉驗(yàn)證法（CV）優(yōu)化網(wǎng)絡(luò)，并利用數(shù)值試驗(yàn)討論了該方法的預(yù)報(bào)精度。

1 RBF神經(jīng)網(wǎng)絡(luò)

1985年，Powell提出了多變量插值的徑向基函數(shù)方法。1988年，Moody和Darken提出了一種神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，即RBF神經(jīng)網(wǎng)絡(luò)，它能夠以任意精度逼近任意連續(xù)函數(shù)，與BP神經(jīng)網(wǎng)絡(luò)相比，RBF神經(jīng)網(wǎng)絡(luò)具有學(xué)習(xí)速度快、全局收斂等優(yōu)點(diǎn)。

1．1 RBF神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)

RBF神經(jīng)網(wǎng)絡(luò)是一種三層前饋式神經(jīng)網(wǎng)絡(luò)，包括一個(gè)輸入層、一個(gè)隱含層和一個(gè)輸出層。輸入層由信號(hào)源結(jié)點(diǎn)組成；第二層為隱含層，隱含層神經(jīng)元數(shù)目視所描述問(wèn)題的需要而定，隱含層神經(jīng)元的傳遞函數(shù)是對(duì)中心點(diǎn)徑向?qū)ΨQ(chēng)且衰減的非負(fù)非線(xiàn)性函數(shù)；第三層為輸出層，它對(duì)輸入模式的作用做出響應(yīng)。從輸入空間到隱含層空間的變換是非線(xiàn)性的，而從隱含層空間到輸出層空間的變換是線(xiàn)性的。

RBF網(wǎng)絡(luò)的基本思想是:以徑向基函數(shù)作為隱含層神經(jīng)元的基構(gòu)成隱含層空間，這樣就可以將輸入矢量直接映射到隱含層空間，而不像BP網(wǎng)絡(luò)那樣需要權(quán)連接，當(dāng)RBF的中心點(diǎn)確定以后，這種映射關(guān)系也就確定了。而隱含層空間到輸出空間的映射是線(xiàn)性的，即網(wǎng)絡(luò)的輸出是隱含層神經(jīng)元輸出的線(xiàn)性加權(quán)和，此處的權(quán)即為網(wǎng)絡(luò)可調(diào)參數(shù)。由此可見(jiàn)，從總體上看，網(wǎng)絡(luò)由輸入到輸出的映射是非線(xiàn)性的，而網(wǎng)絡(luò)輸出對(duì)可調(diào)參數(shù)而言卻又是線(xiàn)性的。這樣網(wǎng)絡(luò)的權(quán)就可以由線(xiàn)性方程直接解出，從而大大加快學(xué)習(xí)速度并避免局部極小問(wèn)題。

RBF網(wǎng)絡(luò)的關(guān)鍵在于隱含層節(jié)點(diǎn)徑向基函數(shù)的非線(xiàn)性逼近性能，徑向基函數(shù)表現(xiàn)為多種形式，常見(jiàn)的為高斯函數(shù)，其形式如下:

式中:x為輸入向量；ci為隱含層第i個(gè)節(jié)點(diǎn)的中心向量，與x同維；σi表示高斯函數(shù)的寬度；‖xci‖為向量x－ci的歐式范數(shù)；h表示隱含層神經(jīng)元的個(gè)數(shù)。

輸出層節(jié)點(diǎn)相應(yīng)的輸出可表示為

式中:＾yj為第j個(gè)輸出單元的輸出值；M表示輸出層神經(jīng)元的個(gè)數(shù)；wi，j為第i個(gè)隱單元到第j個(gè)輸出單元的權(quán)值。

1．2 RBF神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)算法

通常，RBF網(wǎng)絡(luò)的學(xué)習(xí)分為2個(gè)階段進(jìn)行，即隱含層的學(xué)習(xí)和輸出層的學(xué)習(xí)2步。

1）隱含層學(xué)習(xí)

隱含層學(xué)習(xí)表現(xiàn)在隱含層神經(jīng)元的數(shù)目h、中心向量c和高斯函數(shù)的寬度σ的確定上。其中，隱含層神經(jīng)元數(shù)目的確定是關(guān)鍵問(wèn)題，傳統(tǒng)的做法是使其與訓(xùn)練集樣本數(shù)目相同；而中心可根據(jù)訓(xùn)練樣本點(diǎn)在自變量空間的分布，選出或計(jì)算出有“代表性“的點(diǎn)作為中心點(diǎn)，再由中心周?chē)臉颖军c(diǎn)確定寬度σ，應(yīng)用較多的為聚類(lèi)分析法。

2）輸出層學(xué)習(xí)

輸出層學(xué)習(xí)體現(xiàn)在隱含層與輸出層間的連接權(quán)w的確定上，可以應(yīng)用最小二乘法對(duì)網(wǎng)絡(luò)輸出權(quán)值進(jìn)行訓(xùn)練，其學(xué)習(xí)訓(xùn)練的目標(biāo)是使總誤差達(dá)到最小。

1．3 RBF神經(jīng)網(wǎng)絡(luò)的MATLAB工具箱函數(shù)

1）RBF神經(jīng)網(wǎng)絡(luò)創(chuàng)建函數(shù)

MATLAB函數(shù)newrbe用于創(chuàng)建一個(gè)精確的RBF神經(jīng)網(wǎng)絡(luò)，其調(diào)用格式如下:

net＝newrbe（P，T，spread），

其中，P為網(wǎng)絡(luò)的輸入向量；T為網(wǎng)絡(luò)目標(biāo)向量（輸出向量）；spread為徑向基函數(shù)的寬度σ（默認(rèn)為1．0）；net為創(chuàng)建好的RBF網(wǎng)絡(luò)。

在MATLAB中，沒(méi)有專(zhuān)門(mén)的RBF神經(jīng)網(wǎng)絡(luò)訓(xùn)練和學(xué)習(xí)函數(shù)，這是因?yàn)樵谑褂煤瘮?shù)newrbe創(chuàng)建RBF網(wǎng)絡(luò)的過(guò)程中，newrbe函數(shù)就根據(jù)訓(xùn)練樣本P、T確定了網(wǎng)絡(luò)結(jié)構(gòu)，即隱含層神經(jīng)元數(shù)、權(quán)值和閾值等網(wǎng)絡(luò)參數(shù)。函數(shù)newrbe需要調(diào)節(jié)的參數(shù)只有一個(gè)spread，一般根據(jù)先驗(yàn)知識(shí)選取。newrbe函數(shù)創(chuàng)建網(wǎng)絡(luò)的過(guò)程為

①隱層神經(jīng)元數(shù)目等于輸入樣本數(shù)，其權(quán)值等于輸入向量的轉(zhuǎn)置:

式中:IW為隱層神經(jīng)元的權(quán)值；PT為輸入向量P的轉(zhuǎn)置。

②設(shè)置隱層神經(jīng)元的閾值為b

③以隱含層神經(jīng)元的輸出作為輸出層神經(jīng)元的輸入，確定輸出層神經(jīng)元的權(quán)值和閾值，使之滿(mǎn)足:

式中:W｛2，1｝、b｛2｝分別為輸出層神經(jīng)元的權(quán)值和閾值；A｛1｝為隱含層神經(jīng)元的輸出。

可以看出，上述過(guò)程只要進(jìn)行一次就可以得到一個(gè)零誤差的徑向基網(wǎng)絡(luò)，所以以newrbe函數(shù)創(chuàng)建網(wǎng)絡(luò)的速度是比較快的，但由于隱層神經(jīng)元數(shù)等于輸入樣本數(shù)，如果輸入向量數(shù)目很大時(shí)，將導(dǎo)致網(wǎng)絡(luò)的規(guī)模也很大。

2）RBF神經(jīng)網(wǎng)絡(luò)仿真函數(shù)

函數(shù)sim用于利用已經(jīng)創(chuàng)建好的RBF神經(jīng)網(wǎng)絡(luò)進(jìn)行仿真預(yù)測(cè)，其調(diào)用格式為

＾y＝sim（net，＾P），

其中，net為創(chuàng)建好的RBF網(wǎng)絡(luò)；＾P為網(wǎng)絡(luò)的輸入向量；＾y為網(wǎng)絡(luò)輸出向量。

2 交叉驗(yàn)證法

神經(jīng)網(wǎng)絡(luò)對(duì)訓(xùn)練樣本的逼近稱(chēng)為擬合，而對(duì)非訓(xùn)練樣本的逼近稱(chēng)為泛化（預(yù)報(bào)）。顯然，衡量一個(gè)網(wǎng)絡(luò)訓(xùn)練好壞的標(biāo)準(zhǔn)不僅是要考察網(wǎng)絡(luò)的擬合能力，更要檢驗(yàn)其泛化能力，即訓(xùn)練神經(jīng)網(wǎng)絡(luò)的最終目的是確保訓(xùn)練好的網(wǎng)絡(luò)同時(shí)具有良好的擬合能力及泛化能力。一般情況下，訓(xùn)練能力差時(shí)網(wǎng)絡(luò)的泛化能力也弱，并且在一定程度上隨著訓(xùn)練能力的提高泛化能力也會(huì)隨之提高，但是這種趨勢(shì)有一個(gè)極限，達(dá)到此極限后，隨著訓(xùn)練能力的提高，泛化能力反而下降，即出現(xiàn)所謂過(guò)適性或過(guò)擬合，也稱(chēng)為過(guò)學(xué)習(xí)，此時(shí)網(wǎng)絡(luò)學(xué)習(xí)了過(guò)多的樣本細(xì)節(jié)而不能反映樣本所蘊(yùn)含的內(nèi)在規(guī)律。欠學(xué)習(xí)。過(guò)學(xué)習(xí)的網(wǎng)絡(luò)都屬于不成熟網(wǎng)絡(luò)，泛化能力會(huì)很差，實(shí)際中網(wǎng)絡(luò)往往是過(guò)擬合。從上述意義上講，網(wǎng)絡(luò)訓(xùn)練過(guò)程中存在兩方面的問(wèn)題:

1）神經(jīng)網(wǎng)絡(luò)對(duì)訓(xùn)練樣本有較強(qiáng)的依賴(lài)性，如何從有限的訓(xùn)練樣本中挖掘更多的信息是網(wǎng)絡(luò)訓(xùn)練好壞的關(guān)鍵；

2）在網(wǎng)絡(luò)訓(xùn)練過(guò)程中要著重考察網(wǎng)絡(luò)的泛化能力，避免發(fā)生網(wǎng)絡(luò)過(guò)學(xué)習(xí)的狀態(tài)。

解決上述問(wèn)題的一個(gè)有效方法是采用交叉驗(yàn)證法優(yōu)化網(wǎng)絡(luò)。交叉驗(yàn)證法的基本思想是在某種意義下將學(xué)習(xí)集進(jìn)行分組，一部分做為訓(xùn)練集，另一部分做為驗(yàn)證集，首先用訓(xùn)練集對(duì)神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練，再利用驗(yàn)證集來(lái)測(cè)試訓(xùn)練得到的模型，以此來(lái)做為評(píng)價(jià)網(wǎng)絡(luò)性能的指標(biāo)。交叉驗(yàn)證法的一種改進(jìn)方法是 K 重交叉驗(yàn)證法（K－fold Cross－Validation）［11］，其工作方式如下:假定有個(gè) N 樣本，將其隨機(jī)均勻分成K份，輪流選擇其中K－1份作為訓(xùn)練樣本對(duì)網(wǎng)絡(luò)進(jìn)行訓(xùn)練，剩余的1份作為驗(yàn)證（其樣本個(gè)數(shù)為N／K）校本以檢驗(yàn)網(wǎng)絡(luò)的泛化性能。網(wǎng)絡(luò)需要訓(xùn)練K次，也就可以創(chuàng)建K個(gè)網(wǎng)絡(luò)，選擇泛化誤差最小時(shí)的網(wǎng)絡(luò)模型為最優(yōu)模型，用均方根誤差（RMSE）來(lái)衡量泛化誤差:

式中:RMSE（j）表示第j次訓(xùn)練的網(wǎng)絡(luò)泛化誤差；j＝1，2，…，K；yi為理論值；＾yi為網(wǎng)絡(luò)輸出值。

3 基于RBF神經(jīng)網(wǎng)絡(luò)的鐘差預(yù)報(bào)

3．1 利用RBF神經(jīng)網(wǎng)絡(luò)進(jìn)行鐘差預(yù)報(bào)的過(guò)程

利用RBF神經(jīng)網(wǎng)絡(luò)方法進(jìn)行鐘差預(yù)報(bào)的過(guò)程為:

1）確定網(wǎng)絡(luò)的輸入，輸出量及相應(yīng)的維數(shù)。輸入量為已知的鐘差值，輸出量為鐘差預(yù)報(bào)值。

2）根據(jù)已知數(shù)據(jù)，進(jìn)行學(xué)習(xí)集與測(cè)試集的分配。

3）數(shù)據(jù)的歸一化，由于各歷元鐘差數(shù)據(jù)相差懸殊，即輸入變量差異較大，為了使網(wǎng)絡(luò)訓(xùn)練一開(kāi)始就給各輸入分量以同等重要的地位，需要將輸入數(shù)據(jù)變換到同一范圍中。選用由MATLAB提供的premnmx函數(shù)，將輸入數(shù)據(jù)和與之對(duì)應(yīng)的輸出數(shù)據(jù)進(jìn)行歸一化處理，使處理后數(shù)據(jù)全部落在（－1，1）的范圍內(nèi)，這樣有利于提高神經(jīng)網(wǎng)絡(luò)的收斂速度，節(jié)省網(wǎng)絡(luò)訓(xùn)練耗時(shí)，并且可以防止發(fā)生突變的鐘差數(shù)據(jù)對(duì)網(wǎng)絡(luò)訓(xùn)練的影響。

premnmx函數(shù)的模型為

式中:p為需要?dú)w一化的數(shù)據(jù)；pn為歸一化后的數(shù)據(jù)；minp、maxp分別表示p中的最大值和最小值。

4）網(wǎng)絡(luò)創(chuàng)建和訓(xùn)練。采用MATLAB神經(jīng)網(wǎng)絡(luò)工具箱函數(shù)newrbe建立RBF神經(jīng)網(wǎng)絡(luò)，并根據(jù)學(xué)習(xí)集采用K重交叉驗(yàn)證法優(yōu)化網(wǎng)絡(luò)模型。

5）網(wǎng)絡(luò)仿真測(cè)試。利用訓(xùn)練好的RBF神經(jīng)網(wǎng)絡(luò)進(jìn)行鐘差預(yù)報(bào)。

6）預(yù)報(bào)結(jié)果的逆歸一化處理，并對(duì)還原結(jié)果進(jìn)行檢驗(yàn)。

3．2 試驗(yàn)結(jié)果與分析

為了分析RBF神經(jīng)網(wǎng)絡(luò)模型預(yù)報(bào)衛(wèi)星鐘差的精度，選取3顆鐘差變化較典型的衛(wèi)星，分別是PRN03、PRN05和PRN08，其中PRN03呈遞增變化，PRN05呈遞減變化，而PRN08既呈遞增也呈遞減變化，具體的鐘差變化如圖4～6所示。采用的數(shù)據(jù)為2010年11月28日至30日的IGS事后GPS精密鐘差，采樣間隔為5min．由于灰色模型的鐘差預(yù)報(bào)效果優(yōu)于二次多項(xiàng)式模型［3］，為了驗(yàn)證RBF網(wǎng)絡(luò)的效果，采用灰色模型GM（1，1）對(duì)衛(wèi)星鐘差預(yù)報(bào)，并將兩者加以比較。其中，3．2．1節(jié)給出了基于交叉驗(yàn)證法的網(wǎng)絡(luò)訓(xùn)練結(jié)果；3．2．2節(jié)將傳統(tǒng)灰色系統(tǒng)模型與RBF神經(jīng)網(wǎng)絡(luò)模型的預(yù)報(bào)精度進(jìn)行對(duì)比分析。

3．2．1 交叉驗(yàn)證法試驗(yàn)分析

算例中將28日和29日2天的數(shù)據(jù)作為學(xué)習(xí)集用于RBF神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)訓(xùn)練，將30日的數(shù)據(jù)作為測(cè)試集以檢驗(yàn)?zāi)Ｐ托Ч?，采?重交叉驗(yàn)證法將學(xué)習(xí)樣本隨機(jī)分成4份，其中3份作為訓(xùn)練樣本，剩余1份作為驗(yàn)證樣本，由于是用MATLAB函數(shù)newrbe實(shí)現(xiàn)RBF神經(jīng)網(wǎng)絡(luò)的鐘差建模，因此每次訓(xùn)練時(shí)需要調(diào)節(jié)的參數(shù)只一個(gè)spread，即徑向基函數(shù)的寬度σ，在（0．1，5）的范圍內(nèi)選取，步長(zhǎng)取0．1，對(duì)于每次訓(xùn)練，驗(yàn)證誤差（泛化誤差）最小時(shí)所對(duì)應(yīng)的值即為最優(yōu)值。只給出PRN08衛(wèi)星的鐘差在訓(xùn)練過(guò)程中訓(xùn)練誤差（擬合誤差）和驗(yàn)證誤差的變化情況，如圖1所示，其中橫坐標(biāo)代表徑向基函數(shù)的寬度σ（spread），縱坐標(biāo)代表均方根誤差（RMSE）。

由圖1可以看出，每次訓(xùn)練的結(jié)果是訓(xùn)練誤差幾乎相同，都不超過(guò)1ps，但驗(yàn)證誤差卻不盡相同，這說(shuō)明網(wǎng)絡(luò)的擬合能力與其泛化能力在某種程度上并沒(méi)有必然的聯(lián)系。為了更好、更直觀的反映網(wǎng)絡(luò)訓(xùn)練效果，將圖1中的驗(yàn)證誤差圖（圖1b）的部分區(qū)域放大作圖，放大后如圖2（a）、圖2（b）所示。另外，表1給出了每次訓(xùn)練的最小驗(yàn)證誤差。

結(jié)合圖1、2，分析表1，可以更清楚地看出:

1）spread對(duì)網(wǎng)絡(luò)性能有重要的影響。實(shí)際應(yīng)用中，對(duì)于特定的衛(wèi)星鐘差數(shù)據(jù)，spread值的調(diào)節(jié)過(guò)程有一定的規(guī)律，對(duì)于PRN08衛(wèi)星鐘差數(shù)據(jù)，在一定的范圍內(nèi)，spread值越大網(wǎng)絡(luò)性能越好，另外，每次訓(xùn)練得到的最優(yōu)spread值并不相同；

2）隨著訓(xùn)練次數(shù)的增加，擬合誤差幾乎沒(méi)有變化，而驗(yàn)證誤差則是先減小后增大，存在最小點(diǎn)，對(duì)應(yīng)的最小值為0．211 4ns，出現(xiàn)在第2次訓(xùn)練過(guò)程中，對(duì)應(yīng)的擬合誤差和spread值分別為0．62 ps、2．4；

3）交叉驗(yàn)證法通過(guò)多次訓(xùn)練網(wǎng)絡(luò)能夠優(yōu)化網(wǎng)絡(luò)模型，比只訓(xùn)練一次的泛化誤差更小。

表1 驗(yàn)證誤差／ns

為了驗(yàn)證交叉驗(yàn)證法的有效性，突出其優(yōu)點(diǎn)，將其同常規(guī)方法（將學(xué)習(xí)樣本全部用來(lái)訓(xùn)練網(wǎng)絡(luò)）進(jìn)行對(duì)比，圖3示出了兩種方法對(duì)同一測(cè)試樣本的測(cè)試誤差。對(duì)于常規(guī)方法，測(cè)試誤差RMSE為17．30ns，而對(duì)于交叉驗(yàn)證法，測(cè)試誤差RMSE僅為6．03ns，顯然，交叉驗(yàn)證法可以明顯提高網(wǎng)絡(luò)的泛化能力。

圖3 兩種方法的測(cè)試誤差

3．2．2 灰色系統(tǒng)模型與RBF神經(jīng)網(wǎng)絡(luò)模型預(yù)報(bào)精度對(duì)比分析

分別利用灰色系統(tǒng)模型GM（1，1）和基于交叉驗(yàn)證法的RBF神經(jīng)網(wǎng)絡(luò)模型對(duì)3顆衛(wèi)星的鐘差作24h的短期預(yù)報(bào)并進(jìn)行誤差比較，如圖4～6所示。

由圖4～6可以看出，傳統(tǒng)灰色模型對(duì)不同衛(wèi)星鐘差在短期預(yù)報(bào)誤差異常波動(dòng)很大，對(duì)鐘差數(shù)據(jù)類(lèi)型有較強(qiáng)的依賴(lài)性，因此對(duì)于衛(wèi)星鐘差短期預(yù)報(bào)而言，RBF神經(jīng)網(wǎng)絡(luò)模型的預(yù)報(bào)精度和穩(wěn)定性要明顯優(yōu)于傳統(tǒng)灰色模型，可以在實(shí)際應(yīng)用中來(lái)代替灰色模型。

為了進(jìn)一步研究RBF神經(jīng)網(wǎng)絡(luò)模型在衛(wèi)星鐘差短期預(yù)報(bào)中的性能，對(duì)上述預(yù)報(bào)結(jié)果作誤差統(tǒng)計(jì)特性分析，如表2所示，其中MAX代表最大誤差，MEAN代表平均誤差，STD代表標(biāo)準(zhǔn)差。

表2 預(yù)報(bào)誤差統(tǒng)計(jì)／ns

由表2可以更清楚地看出，對(duì)于衛(wèi)星鐘差短期預(yù)報(bào)來(lái)講，RBF神經(jīng)網(wǎng)絡(luò)模型較傳統(tǒng)灰色模型有更高的預(yù)報(bào)精度和更好的穩(wěn)定性，然而RBF神經(jīng)網(wǎng)絡(luò)模型也有一定的局限性，主要是對(duì)訓(xùn)練樣本有一定的要求，體現(xiàn)在以下2個(gè)方面:

1）對(duì)訓(xùn)練樣本的數(shù)量有一定要求，如果訓(xùn)練樣本過(guò)少，網(wǎng)絡(luò)的學(xué)習(xí)效果就不是很好，相應(yīng)地就會(huì)影響網(wǎng)絡(luò)的泛化能力，但訓(xùn)練樣本也不能太多，過(guò)多的樣本將會(huì)導(dǎo)致過(guò)長(zhǎng)的網(wǎng)絡(luò)訓(xùn)練耗時(shí)；

2）要求訓(xùn)練樣本數(shù)據(jù)具有一定的代表性，即能夠表征系統(tǒng)的變化特性，否則網(wǎng)絡(luò)的泛化能力也就不強(qiáng)。

4 結(jié) 論

通過(guò)以上對(duì)RBF神經(jīng)網(wǎng)絡(luò)在衛(wèi)星鐘差預(yù)報(bào)中的分析可以看出，RBF神經(jīng)網(wǎng)絡(luò)模型較傳統(tǒng)的灰色系統(tǒng)模型在鐘差預(yù)報(bào)方面有一定的優(yōu)勢(shì)，模型精度和穩(wěn)定性均有所提高。在MATLAB環(huán)境下對(duì)RBF神經(jīng)網(wǎng)絡(luò)建模，對(duì)參數(shù)要求較低，需要調(diào)節(jié)的只有一個(gè)徑向基函數(shù)的寬度參數(shù)，容易實(shí)現(xiàn)。但模型本身也有一定的局限性，對(duì)樣本數(shù)據(jù)依賴(lài)性強(qiáng)，且容易學(xué)習(xí)過(guò)度，對(duì)于后者，可采用交叉驗(yàn)證法訓(xùn)練網(wǎng)絡(luò)避免，試驗(yàn)也表明該方法可以明顯提高網(wǎng)絡(luò)的泛化能力。另外，模型只適用于短期預(yù)報(bào)，如果要用于長(zhǎng)期預(yù)報(bào)，則需要不斷增加新的樣本數(shù)據(jù)對(duì)模型進(jìn)行再訓(xùn)練。

僅僅對(duì)RBF神經(jīng)網(wǎng)絡(luò)在衛(wèi)星鐘差預(yù)報(bào)中應(yīng)用的初步研究，所用方法尚需改進(jìn)和完善，如何對(duì)模型進(jìn)行優(yōu)化以提高預(yù)報(bào)精度，以及如何將模型用于長(zhǎng)期預(yù)報(bào)是下一步的研究重點(diǎn)。

［1］ 吳海濤，李孝輝，盧曉春，等．衛(wèi)星導(dǎo)航系統(tǒng)時(shí)間基礎(chǔ)［M］．北京:科學(xué)出版社，2011．

［2］ 李 黎，匡翠林，朱建軍．基于IGU預(yù)報(bào)軌道實(shí)時(shí)估計(jì)精密衛(wèi)星鐘差［J］．大地測(cè)量與地球動(dòng)力學(xué)，2011，31（2）:111－116．

［3］ 崔先強(qiáng)，焦文海．灰色系統(tǒng)模型在衛(wèi)星鐘差預(yù)報(bào)中的應(yīng)用［J］．武漢大學(xué)學(xué)報(bào)·信息科學(xué)版，2005，30（5）:447－450．

［4］ 鄭作亞，盧秀山．幾種GPS衛(wèi)星鐘差預(yù)報(bào)方法比較及精度分析［J］．山東科技大學(xué)·自然科學(xué)版，2008，27（4）:6－11，15．

［5］ 朱陵鳳，吳曉平，李 超，等．灰色模型在衛(wèi)星鐘差預(yù)報(bào)鐘的缺陷分析［J］．宇航計(jì)測(cè)技術(shù)，2007，27（4）:34，42－44．

［6］ 朱祥維，肖 華，雍少為，等．衛(wèi)星鐘差預(yù)報(bào)的Kalman算法及其性能分析［J］．宇航學(xué)報(bào)，2008，29（3）:966－970，1052．

［7］ 黃觀文，張 勤，許國(guó)昌．基于頻譜分析的IGS精密星歷衛(wèi)星鐘差精度分析研究［J］．武漢大學(xué)學(xué)報(bào)·信息科學(xué)版，2008，33（5）:496－499．

［8］ 李孝輝，柯熙政．原子鐘信號(hào)的神經(jīng)網(wǎng)絡(luò)模型［J］．陜西天文臺(tái)臺(tái)刊，2001，23（2）:110－115．

［9］ 郭承軍，滕云龍．基于小波分析和神經(jīng)網(wǎng)絡(luò)的衛(wèi)星鐘差預(yù)報(bào)性能分析［J］．天文學(xué)報(bào)，2010，51（4）:395－403．

［10］ 王 威，榮志飛，范國(guó)清．基于GRNN的導(dǎo)航衛(wèi)星鐘差仿真方法研究［J］．系統(tǒng)仿真學(xué)報(bào)，2011，23（11）:2339－2345．

［11］ 丁常富，王 亮．基于交叉驗(yàn)證法的BP神經(jīng)網(wǎng)絡(luò)在汽輪機(jī)故障診斷中的應(yīng)用［J］．電力科學(xué)與工程，2008，24（3）:31－34．