基于Cornell的自適應(yīng)電離層閃爍強度的模型研究

易首成, 紀元法, 孫希延, 韋照川, 賈茜子

(桂林電子科技大學信息與通信學院, 桂林 541004)

電離層的電子密度分布不均勻體的百米尺度不均勻結(jié)構(gòu)會造成全球?qū)Ш叫l(wèi)星系統(tǒng)(global navigation satellite system, GNSS)信號的繞射、聚焦和散射,形成電離層閃爍[1-3]。根據(jù)電離層閃爍對信號的不同影響,可將電離層閃爍可分為幅度閃爍和相位閃爍,在地面接收端,不同傳播路徑上的無線電波信號相互干渉,造成電波信號振幅的隨機增強或削弱,稱為幅度閃爍,不規(guī)則結(jié)構(gòu)的折射效應(yīng)或衍射效應(yīng),造成接收到的無線電波信號相位發(fā)生快速起伏,稱為相位閃爍[4-5]。低緯地區(qū)以幅度閃爍為主,而極區(qū)以相位閃爍為主,幅度閃爍不一定伴隨有相位閃爍發(fā)生,而相位閃爍必然對應(yīng)有幅度閃爍[6],幅度和相位的閃爍,致使信號接收機的測量誤差增大、頻繁周跳,嚴重時還會出現(xiàn)頻繁、持續(xù)的衛(wèi)星信號失鎖現(xiàn)象[7]。由于電離層閃爍效應(yīng)具有突發(fā)性、偶發(fā)性和區(qū)域性等特點[8],并且受電離層閃爍影響的GNSS信號很難復(fù)現(xiàn),因此,建立能夠準確描述電離層閃爍特性的仿真模型尤為重要。

目前,能夠近似描述電離層閃爍信號模型有基于多重相位屏技術(shù)提出全球電離層閃爍模型(global ionospheric scintillation model,GISM)[9],幅度和相位具有一定的相關(guān)性的概率密度統(tǒng)計模型AJ-Stanford模型[10],基于赤道地區(qū)的大量實測數(shù)據(jù)的統(tǒng)計特性提出Cornell模型[11],以及可模擬出對接收機環(huán)路有嚴重影響的Canonical模型[12],該電離層閃爍模型被廣泛應(yīng)用于接收機算法研究和性能的評估中。徐昌元[13]針對Cornell模型對閃爍頻率限定不嚴格,提出基于Cornell模型三參數(shù)電離層閃爍信號仿真模型。安盼盼[14]基于Cornell模型提出了多頻點電離層閃爍模型。盧丹等[15]基于Cornell模型提出輸入?yún)?shù)為電子總含量(total electron content,TEC)和導(dǎo)航數(shù)據(jù)的電離層閃爍模型。但在使用Cornell模型結(jié)合衛(wèi)星信號模擬器生成的閃爍信號與實際設(shè)計的輸入?yún)?shù)存在誤差,未考慮幅度和相位的更新速率對模擬的閃爍信號的影響。因此,使用Cornell模型生成不同數(shù)據(jù)更新周期和閃爍強度的電離層閃爍信號,并周期的估算閃爍信號幅度閃爍指數(shù)S4的大小。在相同的數(shù)據(jù)更新周期下,隨著電離層閃爍的增強,估算的幅度閃爍指數(shù)S4與模型輸入的指數(shù)S4誤差越大;在同一閃爍條件下,數(shù)據(jù)更新周期越大,兩者的誤差越大,并隨著電離層閃爍的增強而增強。

由于研究電離層閃爍,重要的是電離層閃爍仿真模型生成的閃爍信號能夠描述實際的閃爍變化,針對上述問題,提出基于Cornell模型的自適應(yīng)指數(shù)S4的電離層閃爍模型,稱為AS4-Cornell模型;通過仿真以及使用全球定位系統(tǒng)(global positioning system,GPS)模擬器模擬基于Cornell模型和AS4-Cornell模型的電離層閃爍信號,使用GNSS信號接收機提取載噪比,估算兩個模型的電離層幅度閃爍指數(shù)S4,通過數(shù)據(jù)分析,驗證AS4-Cornell模型的正確性以及對電離層幅度閃爍指數(shù)S4的改善性能。改善后的模型有助于低緯地區(qū)電離層閃爍模型的研究和建立有一定的參考意義。

1 電離層閃爍模型仿真

1.1 Cornell模型原理

Cornell電離層閃爍模型生成的閃爍信號的數(shù)學模型為

(1)

概率密度分布服從均值為零的高斯分布,閃爍引起的信號強度變化δI=δA2。Cornell模型生成電離層閃爍序列的實現(xiàn)框圖如圖1[8]所示。

n(t)為復(fù)高斯噪聲;τ0為去相關(guān)時間;為歸一化前的時變分量 為歸一化前的直接分量;為歸一化前的閃爍序列; z(t)為歸一化后的閃爍序列

1.2 AS4-Cornell模型原理

針對Cornell模型在同一閃爍條件下,數(shù)據(jù)更新周期越長,估計和模型輸入的幅度閃爍指數(shù)S4的誤差越大,并隨著電離層閃爍的增強而增大的問題,提出AS4-Cornell模型。其實現(xiàn)框圖如圖2所示,在Cornell模型的基礎(chǔ)上,增加復(fù)高斯白噪聲權(quán)值分配和權(quán)值更新,三者構(gòu)成一個閉環(huán)系統(tǒng)。

nr(t)為更新后的復(fù)高斯信號;MAX_ΔS4為最大S4偏差; zf(t)為閃爍信號;GW(n)為權(quán)值

1.2.1 BP神經(jīng)網(wǎng)絡(luò)PID原理

BP(back propagation)神經(jīng)網(wǎng)絡(luò)由信號的正向傳播和誤差的反向傳播構(gòu)成,即誤差的反向傳播算法[16]。采用的BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)如圖3所示,它分為輸入層、隱含層和輸出層[17-18]。輸出層的神經(jīng)節(jié)點有4個,其中,r(k)為期望幅度閃爍指數(shù)S4,y(k)為實際估算的幅度閃爍指數(shù)S4,e(k)為期望與實際的幅度閃爍指數(shù)S4的偏差;隱含層有5個神經(jīng)節(jié)點,輸出層有kp、ki和kd3個節(jié)點,wij(k)為輸入層到隱含層的調(diào)節(jié)權(quán)值,wlj(k)為隱含層到輸出層的調(diào)節(jié)權(quán)值?；贐P神經(jīng)網(wǎng)絡(luò)的PID控制結(jié)構(gòu)如圖4所示,采用的增量式PID調(diào)節(jié)控制系統(tǒng)的輸出,整個系統(tǒng)的工作原理為,BP神經(jīng)網(wǎng)絡(luò)通過系統(tǒng)反饋的誤差,自學習調(diào)節(jié)加權(quán)系數(shù),輸出最優(yōu)的PID參數(shù),最后通過PID參數(shù)的調(diào)節(jié)使系統(tǒng)達到最優(yōu)的輸出。

j為輸入層節(jié)點;i為隱藏層節(jié)點;l為輸出層節(jié)點;

r(t)為輸入;e(t)為誤差;u(t)為被控參數(shù);y(t)為輸出

1.2.2 AS4-Cornell模型

通過Cornell模型和上述BP神經(jīng)網(wǎng)絡(luò)PID算法建立AS4-Cornell模型。首先根據(jù)幅度閃爍指數(shù)S4的估算周期對閃爍信號和復(fù)高斯白噪聲進行分段;然后估算每段閃爍信號的幅度閃爍指數(shù)S4以及與模型輸入幅度閃爍指數(shù)S4的偏差;最后采用BP神經(jīng)網(wǎng)絡(luò)PID算法根據(jù)偏差計算每段復(fù)高斯白噪聲的權(quán)值,并對復(fù)高斯白噪聲進行調(diào)整,使模型最終產(chǎn)生的閃爍信號滿足閃爍指數(shù)S4的指標,具體如下。

根據(jù)幅度閃爍指數(shù)S4的估算周期把閃爍信號z(t)平均分成N段。

(2)

式(2)中:N為整數(shù);zN(t)為周期分段后的閃爍信號;Ta為閃爍信號持續(xù)時間;T為估算周期。

估算每段閃爍信號zN(t)的幅度閃爍指數(shù)S4,S4與信號強度I的關(guān)系可表示為

(3)

式(3)中:〈·〉為均值運算。

計算估算指數(shù)S4N與模型輸入的指數(shù)S4偏差eN,當所有的偏差eN低于最大允許偏差MAX_ΔS4,輸出最終的閃爍信號zf(t),否則以估算指數(shù)S4N模型輸入指數(shù)S4以及兩者誤差eN作為BP神經(jīng)網(wǎng)絡(luò)PID算法輸入層的參數(shù),迭代求出適合系統(tǒng)的最優(yōu)的PID參數(shù)kpN、kiN和kdN,然后由增量式PID原理更新權(quán)值GWN(n),可表示為

GWN(n)=GWN(n-1)+kpN[eN(n)-eN(n-1)]+kiNeN(n)+kdN[eN(n)-2eN(n-1)+eN(n-2)]

(4)

式(4)中:GWN(0)=1。

更新復(fù)高斯白噪聲nr(t),并作為Cornell模型的輸入,重新生成閃爍信號,它與分段后復(fù)高斯白噪聲n(t)={n1(t),n2(t),…,nN(t)}和更新后的權(quán)值GWN(n)關(guān)系為

nr(t)=n(t)GWN(n) ={n1(t)GW1(n),n2(t)GW2(n),…,nN(t)GWN(n)}

(5)

2 模型仿真分析

仿真閃爍信號的配置參數(shù):去相關(guān)時間τ0為1,仿真時間Ta為300 s,發(fā)生弱閃爍時S4取0.3,發(fā)生中等閃爍時S4取0.6,發(fā)生強閃爍時S4取0.8。通過分析Cornell模型和AS4-Cornell模型的時域特性和頻率特性,以確保建立的Cornell仿真模型和AS4-Cornell仿真模型的正確性,同時在不同數(shù)據(jù)更新周期下,估算幅度閃爍指數(shù)S4。

2.1 時域特性分析

Cornell模型和AS4-Cornell模型的時域特性分析主要包括閃爍信號的幅度和相位的概率密度分布,圖5(a)和圖5(b)分別給出電離層發(fā)生弱閃爍、中閃爍和強閃爍時,Cornell模型幅度和相位概率密度分布仿真結(jié)果,圖5(c)和圖5(d)分別給出了AS4-Cornell模型的幅度和相位概率密度分布仿真結(jié)果。

圖5 時域特性仿真結(jié)果Fig.5 Simulation results of time domain characteristics

從圖5(a)和圖5(c)可以看出,Cornell模型和AS4-Cornell模型在不同的閃爍程度下的幅度閃爍序列概率密度分布服從萊斯分布;從圖5(b)和圖5(c)可以看出,相位閃爍序列概率密度分布服從均值為0的高斯分布,仿真的閃爍信號的幅度和相位的閃爍序列概率密度分布與理論相符,并且當去相關(guān)時間保持不變時,隨著電離層閃爍強度的增強,幅度和相位的概率密度分布圖的中心點密度逐漸減小。

2.2 頻域特性分析

Cornell模型和AS4-Cornell模型的頻域特性主要從閃爍信號的幅度和相位的頻譜分布展開分析,并使用最小二乘法線性擬合頻率在0.7～3 Hz的幅度譜和0.8～10 Hz相位譜,其中p為最小二乘法擬合的冪律指數(shù)[19]。圖6和圖7分別給出電離層發(fā)生弱閃爍、中閃爍和強閃爍時,Cornell模型的幅度和相位的頻譜分布的仿真結(jié)果。圖8和圖9分別給出AS4-Cornell模型的幅度和相位的頻譜分布的仿真結(jié)果。

圖6 Cornell模型幅度頻譜Fig.6 Cornell model amplitude spectrum

圖7 Cornell模型相位頻譜Fig.7 Cornell model phase spectrum

圖8 AS4-Cornell模型幅度頻譜Fig.8 AS4-Cornell model amplitude and spectrum

圖9 AS4-Cornell模型相位頻譜Fig.9 AS4-Cornell model phase spectrum

從圖6～圖9中可以看出,Cornell模型和AS4-Cornell模型在不同的閃爍程度下的幅度和相位的頻譜均服從一定的冪律分布。因此,閃爍信號的幅度和相位的頻域特性也與理論相符。

2.3 估算幅度閃爍指數(shù)S4

在不同電離層閃爍強度和數(shù)據(jù)更新周期為1、10、100 ms的實驗條件下,仿真生成1 h的基于Cornell模型的電離層閃爍信號,然后每分鐘估算一次幅度閃爍指數(shù)S4,結(jié)果如圖10所示。

圖10 Cornell閃爍信號的幅度閃爍指數(shù)S4估算結(jié)果Fig.10 Estimation results of amplitude flicker index S4 of Cornell flicker signal

從圖10中可以看出,數(shù)據(jù)更新周期為1 ms時,估算的幅度閃爍指數(shù)S4與模型輸入的指數(shù)S4在不同閃爍條件下的最大誤差分別為0.003 2、0.015 4和0.029 7;數(shù)據(jù)更新周期為10 ms時,在不同閃爍條件下的最大誤差分別為0.013 2、0.028 1、0.055 6;數(shù)據(jù)更新周期為100 ms時,在不同閃爍條件下的最大誤差分別為0.056 3、0.153 8、0.174 3;在同一閃爍條件下,數(shù)據(jù)更新周期越大,兩者的誤差越大,并隨著電離層閃爍的增強而增強;后續(xù)模型在GPS模擬器的應(yīng)用中,這種波動更加明顯,極大影響了實際情況下的電離層閃爍信號的模擬。

分析可知,當數(shù)據(jù)更新周期為100 ms時,估算幅度閃爍指數(shù)S4與模型輸入幅度閃爍指數(shù)S4偏差最大,因此測試AS4-Cornell的改善性能和有效性,只需要在不同電離層閃爍強度、數(shù)據(jù)更新周期為100 ms和最大允許偏差MAX_ΔS4為0.001的條件下,仿真生成1 h的電離層閃爍信號,每分鐘估算一次幅度閃爍指數(shù)S4,AS4-Cornell閃爍信號的結(jié)果如圖11所示。

圖11 AS4-Cornell閃爍信號的幅度閃爍指數(shù)S4估算結(jié)果Fig.11 Results of amplitude scintillation index S4 estimation of AS4-Cornell scintillation signal

圖11所示的估算幅度閃爍指數(shù)S4的結(jié)果表明,在不同的電離層閃爍強度下, AS4-Cornell模型產(chǎn)生的電離層閃爍信號能夠很好地反映模型輸入的電離層閃爍指數(shù)S4的強度,并且與模型輸入幅度閃爍指數(shù)S4的誤差在設(shè)置的最大允許偏差0.001之內(nèi),相對于Cornell模型有了極大的改進。

3 基于GPS信號模擬器的電離層閃爍模型驗證

3.1 模擬器原理

當電離層發(fā)生閃爍時,將會對GNSS信號的幅度造成乘性衰減干擾,對相位造成加性偏移干擾。因此在干擾下,GNSS接收機接收到的單顆GPS衛(wèi)星信號數(shù)學模型為

r(t)=A0δAcos[2π(fIF+fd)t+φ0+ δφ]c(t)d(t)+n(t)

(6)

式(6)中:c(t)為擴頻碼;d(t)為導(dǎo)航電文;fIF為中頻頻率;fd為射頻頻率;δA為電離層閃爍引起的幅度衰減指數(shù);δφ為電離層閃爍引起的相位偏移誤差。

使用電離層閃爍模型產(chǎn)生成電離層閃爍信號,提取閃爍信號中的幅度和相位序列,并應(yīng)用與模擬器的數(shù)字中頻信號生成單元,電離層閃爍信號硬件實現(xiàn)框圖如圖12所示。

圖12 電離層閃爍信號硬件實現(xiàn)框圖Fig.12 Block diagram of hardware implementation of ionospheric scintillation signal

3.2 結(jié)果與分析

為驗證AS4-Cornell模型對幅度閃爍指數(shù)S4的改進性能,模擬的閃爍信號配置為閃爍衛(wèi)星號2,持續(xù)時間為15 min,閃爍序列更新周期為100 ms,去相關(guān)時間為1,載噪比為50 Hz,幅度閃爍指數(shù)S4分別為0.3、0.6、0.8;首先使用GPS模擬器分別模擬基于Cornell模型和AS4-Cornell模型的電離層閃爍的GPS信號,然后分析不同電離層閃爍強度下的GPS信號的星座圖[20],結(jié)果如圖13～圖15所示。

I為同相信號;Q為正交信號

從圖13中可以看出,電離層發(fā)生弱閃爍時,信號的幅度和相位擾動小,星座圖呈環(huán)形分布;電離層發(fā)生中閃爍和強閃爍時,信號的幅度和相位擾動增大,星座圖呈散點云狀,隨著閃爍強度的增強,幅度和相位的擾動劇烈,導(dǎo)致接收機的鑒相器工作點偏離0附近,I支路的幅度值出現(xiàn)180度跳變和波動值接近于0的情況增多,信號發(fā)生深度衰落。

載噪比變化如圖14所示。可以看出,隨著電離層閃爍的增強,載噪比的波動范圍變大。當發(fā)生電離層閃爍時,載噪比為50 Hz,當發(fā)生電離層弱閃爍時,載噪比有3～6 Hz的波動;當發(fā)生電離層中閃爍時,載噪比有6 ～22 Hz的波動;當發(fā)生電離層強閃爍時,載噪比有7 ～26 Hz的波動。

CNO為載噪比;N為歷元

每分鐘估算的S4指數(shù)如圖15所示?？梢钥闯?1 min和17 min為未模擬電離層閃爍階段,電離層幅度閃爍指數(shù)S4為0,在2～16 min分別模擬電離層發(fā)生弱閃爍、中閃爍以及強閃爍,兩種模型的幅度閃爍指數(shù)S4的估算結(jié)果在預(yù)設(shè)值附近,但基于AS4-Cornell模型模擬生成的閃爍信號的電離層幅度閃爍指數(shù)S4與模擬器預(yù)設(shè)值的最大偏差為0.09,從而驗證了改進模型的有效性?？紤]到閃爍信號的幅度閃爍序列和相位閃爍序列是基于仿真的閉環(huán)條件下生成的,而不是基于實際硬件環(huán)境生成,因此指數(shù)S4的偏差未達到最大允許偏差的要求。

圖15 PRN2 S4指數(shù)估算結(jié)果Fig.15 Estimated results of PRN2 S4 index

4 結(jié)論

為了改善傳統(tǒng)的Cornell電離層閃爍模型的幅度閃爍指數(shù)S4與模型輸入S4在周期估算時存在偏差的問題,提出AS4-Cornell電離層閃爍模型,并通過仿真和實際閃爍信號的模擬得出以下結(jié)論。

(1)仿真結(jié)果表明,AS4-Cornell模型生成的閃爍信號的幅度閃爍序列時域特性服從萊斯分布,相位閃爍序列時域特性服從均值為零的高斯分布;幅度和相位的頻域特性符合冪律分布,與理論相符。估算的S4指數(shù)與模型輸入的S4指數(shù)誤差可以限定在設(shè)定的最大偏差內(nèi)。

(2)對GPS模擬器模擬生成的閃爍信號進行S4指數(shù)的估算,結(jié)果表明,相比Cornell模型,基于AS4-Cornell模型生成的閃爍信號更符合設(shè)置的幅度閃爍強度。