基于GPS的電離層研究

1、中 國(guó) 地 質(zhì) 大 學(xué)本科生課程論文封面課程名稱：全球定位系統(tǒng)原理及應(yīng)用班 級(jí)： 姓 名： xxxx 學(xué) 號(hào)： 專 業(yè)： 地理信息科學(xué) 日期: 2015 年 1 月 19日評(píng) 語(yǔ)對(duì)課程論文的評(píng)語(yǔ):平時(shí)成績(jī)：課程論文成績(jī)：總 成 績(jī)：評(píng)閱人簽名：注：1、無(wú)評(píng)閱人簽名成績(jī)無(wú)效；2、必須用鋼筆或圓珠筆批閱，用鉛筆閱卷無(wú)效；3、如有平時(shí)成績(jī)，必須在上面評(píng)分表中標(biāo)出，并計(jì)算入總成績(jī)。 基于GPS的電離層研究 摘 要 本文主要介紹了電離層的基本特性，闡述了利用GPS觀測(cè)量計(jì)算電離層TEC的三種方法的基本原理：包括碼觀測(cè)解算電離層TEC的基本原理，載波相位觀測(cè)電離層TEC的基本原理以及碼觀測(cè)和載波相位觀測(cè)聯(lián)

2、合解算電離層TEC的基本原理。最后簡(jiǎn)單的介紹了基于基于GPS的三維電離層層析技術(shù)。電離層層析成像技術(shù)是計(jì)算機(jī)層析成像技術(shù)在電離層監(jiān)測(cè)中的一種新的應(yīng)用。該技術(shù)通過(guò)對(duì)電離層進(jìn)行分層研究，不僅克服了薄層假設(shè)電離層層析模型的局限性，也克服了經(jīng)驗(yàn)?zāi)Ｐ团c傳統(tǒng)地面探測(cè)手段的局限性，特別適合于監(jiān)測(cè)大尺度電離層電子密度垂直分布及其擾動(dòng)狀態(tài)。關(guān)鍵詞：GPS 電離層 TEC目 錄第一章 引言1第二章 電離層的基本特性22.1 電離層結(jié)構(gòu)22.2 電離層騷擾2第三章 計(jì)算TEC的基本理論33.1 TEC簡(jiǎn)介33.2 GPS觀測(cè)量及觀測(cè)方程33.2.1 偽距及碼觀測(cè)方程33.2.2 載波相位及其觀測(cè)方程43.3 碼觀測(cè)

3、解算電離層TEC的基本原理43.4 載波相位觀測(cè)解算電離層TEC的基本原理53.5 GPS碼與相位觀測(cè)聯(lián)合解算電離層TEC63.5.1 前因63.5.2 基本原理6第四章 基于GPS的三維電離層層析技術(shù)9總 結(jié)11參考文獻(xiàn)12 第一章 引言 電離層（Ionosphere）是地球大氣的一個(gè)電離區(qū)域。電離層（ionosphere）是受太陽(yáng)高能輻射以及宇宙線的激勵(lì)而電離的大氣高層。60千米以上的整個(gè)地球大氣層都處于部分電離或完全電離的狀態(tài)，電離層是部分電離的大氣區(qū)域，完全電離的大氣區(qū)域稱磁層。也有人把整個(gè)電離的大氣稱為電離層，這樣就把磁層看作電離層的一部分。除地球外，金星、火星和木星都有電離層。電離

4、層從離地面約50公里開(kāi)始一直伸展到約1000公里高度的地球高層大氣空域，其中存在相當(dāng)多的自由電子和離子，能使無(wú)線電波改變傳播速度，發(fā)生折射、反射和散射，產(chǎn)生極化面的旋轉(zhuǎn)并受到不同程度的吸收1。 人類對(duì)電離層的研究歷史已經(jīng)持續(xù)了一個(gè)多世紀(jì)，隨著科技的不斷進(jìn)步，電離層探測(cè)技術(shù)與方法不斷的發(fā)展。電離層形態(tài)監(jiān)測(cè)有助于加深認(rèn)識(shí)電離層活動(dòng)規(guī)律及其變化機(jī)制。長(zhǎng)期以來(lái)，人們對(duì)電離層形態(tài)的監(jiān)測(cè)與研究主要借助于電離層測(cè)高儀，雷達(dá)以及探空火箭等觀測(cè)手段，并以電離層峰值電子密度（NmF2)，峰值高度（hmF2）以及臨界頻率為研究對(duì)象。近二十年來(lái)，以美國(guó)全球定位系統(tǒng)（GPS）為代表的全球?qū)Ш叫l(wèi)星系統(tǒng)（GNSS）觀測(cè)技術(shù)

5、的興起以及全球分布的GNSS跟蹤站連續(xù)觀測(cè)的海量數(shù)據(jù)的累積，使得利用GNSS監(jiān)測(cè)電離層時(shí)空變化規(guī)律成為可能。電離層層析成像技術(shù)是計(jì)算機(jī)層析成像技術(shù)在電離層監(jiān)測(cè)中的一種新的應(yīng)用2。該技術(shù)通過(guò)對(duì)電離層進(jìn)行分層研究，不僅克服了薄層假設(shè)電離層層析模型的局限性，也克服了經(jīng)驗(yàn)?zāi)Ｐ团c傳統(tǒng)地面探測(cè)手段的局限性，特別適合于監(jiān)測(cè)大尺度電離層電子密度垂直分布及其擾動(dòng)狀態(tài)。近十多年來(lái)，基于GPS的電離層層析技術(shù)逐漸發(fā)展起來(lái)。借助于GPS的巨大優(yōu)勢(shì)，基于GPS的電離層層析技術(shù)可以實(shí)現(xiàn)三維甚至四維電離層結(jié)構(gòu)的重構(gòu)，從而克服了二維電離層層析模型的局限性，并逐漸形成了一種新的電離層監(jiān)測(cè)手段，在電離層形態(tài)與擾動(dòng)監(jiān)測(cè)研究方面具有

6、重大的科學(xué)意義和應(yīng)用價(jià)值，是目前大地測(cè)量，空間物理與無(wú)線電科學(xué)等領(lǐng)域中的一個(gè)研究熱點(diǎn)和重點(diǎn)。第二章 電離層的基本特性 2.1 電離層結(jié)構(gòu) 電離層結(jié)構(gòu)可用電離層特性參量電子密度、離子密度、電子溫度、離子溫度等的空間分布來(lái)表征3。但其研究主要是電子密度隨高度的分布.電子密度（或稱電子濃度）是指單位體積的自由電子數(shù).電子密度隨高度的變化與各高度上大氣成分、大氣密度以及太陽(yáng)輻射通量等因素有關(guān)。 電離層在垂直方向上呈分層結(jié)構(gòu),一般劃分為D層 、 E 層和F層,F層又分為F1 層和F2層.最大電子密度約為106厘米-3,大約位于300千米高度附近.除正規(guī)層次外,電離層區(qū)域還存在不均勻結(jié)構(gòu),如偶發(fā)E層（Es

7、）和擴(kuò)展F.偶發(fā)E層較常見(jiàn),是出現(xiàn)于E層區(qū)域的不均勻結(jié)構(gòu).厚度從幾百米至一二千米,水平延伸一般為0.10千米,高度大約在110千米處,最大電子密度可達(dá)106厘米-3.擴(kuò)展F是一種出現(xiàn)于F層的不均勻結(jié)構(gòu),在赤道地區(qū),常沿地磁方向延伸,分布于2501000千米或更高的電離層區(qū)域.電離層分層結(jié)構(gòu)只是電離層狀態(tài)的理想描述,實(shí)際上電離層總是隨緯度、經(jīng)度呈現(xiàn)復(fù)雜的空間變化,并且具有晝夜、季節(jié)、年、太陽(yáng)黑子周等變化.由于電離層各層的化學(xué)結(jié)構(gòu)、熱結(jié)構(gòu)不同,各層的形態(tài)變化也不盡相同。2.2 電離層騷擾 太陽(yáng)擾動(dòng)以及其他原因?qū)е聦?duì)電離層正常狀態(tài)的顯著偏離4。太陽(yáng)擾動(dòng)引起的電離層騷擾主要有電離層突然騷擾、電離層暴、

8、極蓋吸收、極光帶吸收等.人為因素如核爆炸、大功率發(fā)射機(jī)對(duì)電離層加熱也能引起電離層騷擾.電離層騷擾對(duì)無(wú)線電波傳播會(huì)產(chǎn)生嚴(yán)重影響.電離層突然騷擾.太陽(yáng)色球在耀斑爆發(fā)期間發(fā)出強(qiáng)烈的紫外線和射線輻射 ,使 D 層的電子密度突然增大 ,對(duì)通訊造成嚴(yán)重影響,甚至中斷.突然騷擾持續(xù)時(shí)間一般為幾分鐘至幾小時(shí).電離層暴.F2層狀態(tài)的異常變化.極蓋吸收.在強(qiáng)烈的太陽(yáng)耀斑爆發(fā)時(shí),由太陽(yáng)噴射出來(lái)的高能質(zhì)子流沿地磁力線沉降到極蓋區(qū)上空,使 D 層的電離急劇增大,以至通過(guò)該區(qū)的無(wú)線電波被強(qiáng)烈吸收,常造成無(wú)線電通訊中斷. 極光帶吸收.來(lái)自太陽(yáng)擾動(dòng)區(qū)的高能電子和質(zhì)子沉降到極區(qū)上空,使極光帶低電離層的電離增強(qiáng),以至通過(guò)該區(qū)域的

9、電磁波被強(qiáng)烈吸收.極光帶吸收甚至使電波訊號(hào)中斷。 第三章 計(jì)算TEC的基本理論 3.1 TEC簡(jiǎn)介 電離層電子濃度總含量(TEC)又稱電離層電子濃度柱含量、積分含量等，是一個(gè)非常重要的電離層參量，對(duì)電離層物理的理論研究及電離層電波傳播的應(yīng)用研究均具有十分重要的意義。理論上，TEC的空間分布及時(shí)間變化，反映了電離層的主要特性，因此通過(guò)探測(cè)與分析電離層TEC參量，可以研究電離層不同時(shí)空尺度的分布與變化特性，如電離層擾動(dòng)，電離層的周日、逐日變化，電離層年度變化，以及電離層的長(zhǎng)期變化等5。應(yīng)用中，電離層的TEC與穿透電離層傳播的無(wú)線電波時(shí)間延遲與相位延遲密切相關(guān)，因此可用于在衛(wèi)星定位、導(dǎo)航等空間應(yīng)用工

10、程中的電波傳播修正。 電離層TEC探測(cè)手段以衛(wèi)星信標(biāo)測(cè)量為主，如微分多普勒方法、法拉第旋轉(zhuǎn)方法等。隨著全球定位系統(tǒng)（GPS）的使用，采用GPS雙頻信標(biāo)的測(cè)量獲取電離層TEC參量成為當(dāng)前最為重要的和廣泛采用的方法。利用局域與全球的GPS臺(tái)網(wǎng)觀測(cè)，通過(guò)實(shí)時(shí)處理與傳輸，并采用不同的TEC mapping方法，可以進(jìn)行大范圍電離層TEC的現(xiàn)報(bào)，這在當(dāng)前空間天氣研究中具有特別意義3。3.2 GPS觀測(cè)量及觀測(cè)方程3.2.1 偽距及碼觀測(cè)方程 偽距是GPS接收機(jī)產(chǎn)生的偽隨機(jī)信號(hào)與衛(wèi)星發(fā)射的信號(hào)偏移的時(shí)間量乘以光速。當(dāng)接收機(jī)鐘與衛(wèi)星鐘完全同步時(shí)，測(cè)定的就是信號(hào)傳播的時(shí)間延遲量。因信號(hào)傳播受到多種誤差源的影響

11、，所測(cè)得的距離是一種包含多種誤差的觀測(cè)值，故稱為偽距3。偽距精度一般為碼元寬帶的1%，故對(duì)于C/A碼，精度約為2.9m，對(duì)于P碼，精度約為0.29m。 碼的觀測(cè)方程：（1式）P為碼偽距觀測(cè)值；T，R分別代表衛(wèi)星和接收機(jī)；下標(biāo)k=1,2表示兩個(gè)不同的載頻L1和L2；為信號(hào)發(fā)射時(shí)刻衛(wèi)星至接收時(shí)刻的接收機(jī)天線的幾何距離，稱為站星距；c為真空中的光速；rR為接收機(jī)相對(duì)論效應(yīng)；I為電離層折射誤差；Tr為對(duì)流層折射誤差；Bq為衛(wèi)星和接收機(jī)的儀器偏差；M為碼觀測(cè)中的多路徑效應(yīng)；Ba接收機(jī)天線相位中心偏差；為碼量測(cè)噪聲。設(shè)衛(wèi)星在信號(hào)發(fā)射時(shí)刻的坐標(biāo)為，接收機(jī)天線在信號(hào)接收時(shí)刻的坐標(biāo)為，于是站星距可記為：|=(2

12、式)3.2.2 載波相位及其觀測(cè)方程 載波相位觀測(cè)值是GPS接收機(jī)接收到的，具有多普勒頻移的載波信號(hào)與接收機(jī)本身產(chǎn)生的參考載波信號(hào)之間的相位差3。在基于GPS的三維電離層層析中，必須應(yīng)用高精度載波相位觀測(cè)和碼觀測(cè)聯(lián)合解算電離層TEC。載波相位觀測(cè)文件中提供的是衛(wèi)星信號(hào)的相位與接收機(jī)本機(jī)振蕩相位之差，即載波拍頻相位。通常，載波相位觀測(cè)方程記為：（10式）式中：k為相應(yīng)載波波長(zhǎng)；Nk為對(duì)應(yīng)頻率fk的相位觀測(cè)常數(shù)，下標(biāo)k=1,2,表示載波相位L1和L2；bq為衛(wèi)星和接收機(jī)的儀器偏差；m為多路徑效應(yīng)對(duì)相位觀測(cè)的影響；e為相位觀測(cè)值的噪聲；其他符號(hào)的含義與1式相同。 3.3 碼觀測(cè)解算電離層TEC的基本

13、原理 由上文中的碼的觀測(cè)方程1式和（相速度）給出的觀測(cè)量，在忽略上標(biāo)i,j的前提下，通過(guò)頻率間求差，可得： （3式）（4式）在3式和4式中，1,2為兩個(gè)頻率上碼測(cè)量的噪聲之差；e1,2為兩個(gè)頻率上載波相位測(cè)量之差。假如用和分別表示碼觀測(cè)中衛(wèi)星和接收機(jī)的差分儀器偏差，及，用分別表示L1和L2兩個(gè)載波相位上的衛(wèi)星和接收機(jī)的差分儀器偏差，即，則3式和4式可分別表示為： （5式） （6式）在GPS碼觀測(cè)中，多路徑效應(yīng)對(duì)觀測(cè)量的影響通常與衛(wèi)星觀測(cè)的截止高度角有關(guān)。在實(shí)際的電離層TEC計(jì)算中，為了減少多路徑效應(yīng)的影響，通常選取一個(gè)合適的截止高度角。對(duì)于雙頻GPS觀測(cè)來(lái)說(shuō)，L1頻率上的電離層延遲量可以表示為

14、： （7式）將7式代入5式和6式中，可得到電離層TEC的表達(dá)式如下： （8式） （9式） 3.4 載波相位觀測(cè)解算電離層TEC的基本原理 利用GPS載波相位測(cè)量值也可以推導(dǎo)出電離層TEC，其形式如下： （11式）上式中，11=L1，22=L2，于是，11式可以表示為： （12式）上式中，L1-L2即為L(zhǎng)4組合，TEC表示由該組合導(dǎo)出的電離層TEC。 由載波相位測(cè)量的L4組合雖然可以求得GPS衛(wèi)星信號(hào)傳播路徑上的電離層TEC,但由于其中含有未知的相位初值，因此，利用載波相位測(cè)量的L4組合只能獲得電離層TEC的相對(duì)變化。通常，根據(jù)載波相位L4組合求得的電離層TEC稱為相對(duì)TEC。對(duì)應(yīng)的，根據(jù)碼觀測(cè)

15、形成的P4組合得到的電離層TEC稱為絕對(duì)TEC。然而，根據(jù)碼觀測(cè)的P4組合得到的絕對(duì)TEC的精度通常不是很高，即使使用精碼（P碼或Y碼），其TEC的測(cè)量精度一般也只能達(dá)到的量級(jí)，即1TECU（TECU為電離層總電子含量的常用單位）。 3.5 GPS碼與相位觀測(cè)聯(lián)合解算電離層TEC3.5.1 前因由以上的分析可得，由載波相位觀測(cè)量的L4組合得到的電離層TEC的相對(duì)精度很高，但其只能反應(yīng)電離層TEC的相對(duì)變化情況；而用碼觀測(cè)的P4組合雖然能得到電離層的絕對(duì)值，但它的精度較低。所以我們需要第三種方法，才能得到高精度的絕對(duì)TEC?？紤]到由載波相位觀測(cè)量的L4組合和用碼觀測(cè)的P4組合這兩種組合計(jì)算電離層

16、TEC的特點(diǎn)以及兩者時(shí)間變化趨勢(shì)的一致性，如果將兩種方法聯(lián)合起來(lái)解算電離層TEC，那么理論上，能夠得到高精度的TEC。3.5.2 基本原理對(duì)于觀測(cè)歷元t，載波相位得到的與碼觀測(cè)得到的之間應(yīng)該有一個(gè)差值，該差值可表示為： （13式）令0=,則上式可變?yōu)椋?（14式）理論上，TEC為一常量。如果GPS相位觀測(cè)中沒(méi)有周跳發(fā)生或?qū)χ芴M(jìn)行了精確的校正后，對(duì)每一衛(wèi)星和接收機(jī)時(shí)，只要其觀測(cè)值可以利用，在每個(gè)歷元上都可以得到一個(gè)TEC，將TEC對(duì)時(shí)間進(jìn)行平滑處理，則可以得到一個(gè)精度更高的TEC數(shù)據(jù)。在觀測(cè)歷元N，TEC遞歸計(jì)算公式可以表示為： = =（15式）通常，衛(wèi)星和接收機(jī)的儀器偏差在幾天之內(nèi)是穩(wěn)定的，

17、因此，TEC在一天之內(nèi)也可以看作一個(gè)常量。綜合13式和15式，即可獲得高精度的絕對(duì)TEC。具體方式如下： （16式）將13式代入15式可得 = （17式）如果GPS的相位觀測(cè)中存在周跳現(xiàn)象，在修復(fù)了周跳以后，17式中最后一項(xiàng)中的模糊度N1和N2均為常量，顧及衛(wèi)星和接收機(jī)的儀器偏差在一天之內(nèi)也可以視作常量，則上市可變化為：（18式）綜合18式和16式可得： （19式） 上式中，P1和P2為每個(gè)歷元的碼觀測(cè)；L1和L2為每個(gè)歷元的相位觀測(cè)；0=9.52，為GPS碼觀測(cè)中衛(wèi)星和接收機(jī)的差分儀器偏差。由大量數(shù)據(jù)進(jìn)行處理后可以得到結(jié)論：在不出現(xiàn)硬件變更的情況下，接收機(jī)與衛(wèi)星發(fā)射器儀器偏差的日變化較為穩(wěn)定

18、，在連續(xù)兩天中，接收機(jī)的儀器偏差變化在1ns以內(nèi)，衛(wèi)星儀器偏差的變化在0.5ns以內(nèi)，因此，19式可以轉(zhuǎn)換為： （20式）根據(jù)上式，求出儀器偏差后，利用碼和相位觀測(cè)的聯(lián)合即可求出高精度的絕對(duì)TEC。 第四章 基于GPS的三維電離層層析技術(shù) 如前所述，電離層斜距TEC為電離層電子密度沿信號(hào)傳播路徑的線積分，一般表示為： TEC= （21式）上式中，Ne為GPS衛(wèi)星與接收機(jī)間信號(hào)傳播路徑上電離層電子密度；r是由經(jīng)度，緯度和高度組成的位置向量；l為GPS衛(wèi)星與接收機(jī)之間的衛(wèi)星信號(hào)傳播路徑?；贕PS的電離層層析問(wèn)題就是根據(jù)反演區(qū)域內(nèi)一系列衛(wèi)星信號(hào)傳播路徑上的TEC信息反演該區(qū)域內(nèi)電離層電子密度的時(shí)空

19、分布。由于GPS衛(wèi)星傳輸?shù)氖歉哳l信號(hào)，因此，在電離層CT系統(tǒng)中，GPS衛(wèi)星信號(hào)的傳播路徑可以近似的看作直線3。由21式可以看出，電離層電子密度與電離層TEC之間是非線性的。在實(shí)際反演過(guò)程中，為了反演的方便，通常利用離散反演方法將待反演的電離層空間離散化，為此，事先需要選取合適的基函數(shù)bj,模型化待反演的電離層電子密度，于是有 （22式）上式中，(j=1,2,3.J)代表基函數(shù)的系數(shù)，J表示所選擇的基函數(shù)的數(shù)量，其大小與反演區(qū)域內(nèi)的像素總素相等。對(duì)每一條射線路徑的電離層TEC測(cè)量值，將22式代入21式，可得：= （i=1,2,.,m）(23式)其中，m代表電離層斜距TEC觀測(cè)值總數(shù)或GPS射線路

20、徑總數(shù)。假如以Aij代表23式中的積分項(xiàng)，則有： （24式）考慮到GPS測(cè)量中噪聲的影響，將24式代入23式，則每條GPS信號(hào)傳播路徑上夫人電離層TEC測(cè)量數(shù)據(jù)可以表達(dá)為： （25式）上式中，i為第i條GPS射線傳播路徑上的觀測(cè)噪聲。為了簡(jiǎn)化計(jì)算，本文采用像素指標(biāo)函數(shù)作為基函數(shù)，即： 1,r像素 b(r)= 0，其他 （26式） 利用像素指標(biāo)函數(shù)作為基函數(shù)后，待反演的電離層空間被離散化為一系列小的三維像素。將26式代入24式后，有： （27式）上式中，表示第i條GPS射線在第j個(gè)像素內(nèi)的截距。當(dāng)GPS射線穿過(guò)該像素時(shí)，Aij=0,反之，當(dāng)像素內(nèi)沒(méi)有任何GPS射線穿過(guò)時(shí)，Aij=0。將25式以矩

21、陣的形式表示，則可得： （28式）式中，y為m條GPS信號(hào)傳播途徑上電離層斜距TEC觀測(cè)值所構(gòu)成的列向量；A為GPS射線在對(duì)應(yīng)像素內(nèi)的截距構(gòu)成的mJ維向量，e為m條GPS射線上的觀測(cè)噪聲和離散誤差組成的m維列向量。 總 結(jié)本文主要介紹了電離層的基本特性，闡述了利用GPS觀測(cè)量計(jì)算電離層TEC的三種方法的基本原理：包括碼觀測(cè)解算電離層TEC的基本原理，載波相位觀測(cè)電離層TEC的基本原理以及碼觀測(cè)和載波相位觀測(cè)聯(lián)合解算電離層TEC的基本原理。最后簡(jiǎn)單的介紹了基于基于GPS的三維電離層層析技術(shù)。利用雙頻GPS測(cè)量電離層TEC是當(dāng)前最為重要的和廣泛采用的方法。而近十多年來(lái)，基于GPS的電離層層析技術(shù)也

22、逐漸發(fā)展起來(lái)。借助于GPS的巨大優(yōu)勢(shì)，基于GPS的電離層層析技術(shù)可以實(shí)現(xiàn)三維甚至四維電離層結(jié)構(gòu)的重構(gòu)，從而克服了二維電離層層析模型的局限性，并逐漸形成了一種新的電離層監(jiān)測(cè)手段，在電離層形態(tài)與擾動(dòng)監(jiān)測(cè)研究方面具有重大的科學(xué)意義和應(yīng)用價(jià)值，是目前大地測(cè)量，空間物理與無(wú)線電科學(xué)等領(lǐng)域中的一個(gè)研究熱點(diǎn)和重點(diǎn)。本文主要參照了基于GNSS的電離層層析算法及其應(yīng)用一書(shū)和其他文獻(xiàn)，由于作者知識(shí)的淺薄與能力的欠缺，本文只是對(duì)有關(guān)文獻(xiàn)資料做了一個(gè)初步的整理，總結(jié)與歸納，存在很多不足之處。參考文獻(xiàn)1 熊年祿，唐存琛，李行健.電離層物理概論.武漢大學(xué)出版社：1999年.35-37.2熊波，萬(wàn)衛(wèi)星，劉立波，等.武漢電離

23、層TEC和NmF2及板厚的季節(jié)變化.見(jiàn)：空間科學(xué)學(xué)報(bào),2007,27（2）：25-131.3聞德保.基于GNSS的電離層層析算法及其應(yīng)用.測(cè)繪出版社：2013年.24-37.4萬(wàn)衛(wèi)星，寧百齊，袁洪.電離層擾動(dòng)的GPS探測(cè).見(jiàn)：空間科學(xué)學(xué)報(bào),1998,18(3):243-2515蔡昌盛，高井祥，李征航.利用GPS監(jiān)測(cè)電離層總電子含量的季節(jié)性變化.見(jiàn)：武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2006年第5期（5）:451-453.6楊東凱，樊江濱，張波.GNSS應(yīng)用與方法.電子工業(yè)出版社:2011年.317-3197BUST G S,COKER C,COCO D.IRI data ingestion and i

24、onospheric tomography.Advances in Space Research,2001,27:157-165.8LIU Zhizhao.Ionospheric tomgraphic modeling and applications using global positioning system(GPS)measurements.Calgary:The University of Calgary.9MANNUCCI A J,WILSON S D,YUAN D N.A global mapping technique for GPS-derived ionospheric total electron content measurements.Radio Science,1998