GAMIT簡介pjh.doc

GAMIT簡介 上海天文臺(tái)程宗頤一， GAMIT的特點(diǎn)用雙差做基本觀測(cè)量。1， 優(yōu)點(diǎn) ： 消除了站鐘和星鐘的主要誤差； 軌道誤差對(duì)測(cè)站相對(duì)位置不敏感，有利于精密定位；2， 缺點(diǎn) ： 有效觀測(cè)減少, 各站有效觀測(cè)因其在網(wǎng)內(nèi)的位置的變化而變化； 不能用于單點(diǎn)定位；3, 在X-Windows 支持下的 UNIX 操作系統(tǒng) 或 LINUX 操作系統(tǒng)下運(yùn)行。二， GAMIT的功能1, 估計(jì)衛(wèi)星狀態(tài)矢量（座標(biāo)和速度）；2, 估計(jì)太陽光壓參數(shù)（39個(gè)）；3, 估計(jì)y-bias；4, 估計(jì) erp 參數(shù)；5, 估計(jì)站坐標(biāo)，可得站坐標(biāo)的時(shí)間序列，確定測(cè)站位置隨時(shí)間而變的規(guī)律，從而可研 究板塊運(yùn)動(dòng)和地殼形變以及其它地球物理現(xiàn)象；6, 估計(jì)對(duì)流層引起的天頂延遲，從而能開展GPS氣象學(xué)方面的應(yīng)用和研究。三， GAMIT軟件的使用簡介。1， 建立以（積日）天命名的工作目錄，rinex,eph目錄和 tables目錄(建立兩個(gè)目錄積日工作目錄與tables目錄即可。積日目錄：比如，2003322是指2003年第322天的工作目錄，里面存放o-files（相位與偽距觀測(cè)文件），n-files（導(dǎo)航文件，衛(wèi)星狀態(tài)和星鐘改正信息），以及該工作日的sp3-file （GPS衛(wèi)星精密星歷，從IGS網(wǎng)站下載）2, 在 rinex 目錄中再建以天命名的觀測(cè)文件目錄，并在觀測(cè)文件目錄中準(zhǔn)備每天的 rinex o-files（相位和偽距），n-file(衛(wèi)星狀態(tài)和星鐘改正信息)，m-file(氣象參數(shù))3, 2 在 tables目錄中建立1 測(cè)站坐標(biāo)初始文件 （L 文件lfile.），地心球坐標(biāo);2 測(cè)站信息文件 （ ），接收機(jī)和天線型號(hào)信息，天線高；3 測(cè)段信息控制文件 （ sestbl. ）4 測(cè)站信息控制文件 （ sittbl. ）5 星號(hào)對(duì)照表 （ svnav.dat ）, 衛(wèi)星質(zhì)量，yaw rate,6 接收機(jī)及其天線型號(hào)對(duì)照表 （ rcvant.dat ）7 天線相位中心改正表 （ antmod.dat ）8 地球形狀參數(shù)表 （ gdetic.dat ）9 太陽表 ( soltab. )10 月亮表 ( luntab. )11 章動(dòng)表 ( nutabl. )12 跳秒表 ( leap.sec )13 周跳的自動(dòng)探測(cè)和修改命令表 ( autcln.cmd )14 海潮表 ( scherneck_grid, scherneck_stations.oct )4, 在 eph 目錄中建立GPS衛(wèi)星精密星歷 (igs sp3 or sp1 ) 15 極移表 （ pole. ）16 TAIUT1 表 （ ut1. ） 執(zhí)行 erp 命令即可產(chǎn)生這兩個(gè)文件。 17 壞衛(wèi)星信息文件 ( svs_exclude.dat )53, 選擇計(jì)算方按案和設(shè)置先驗(yàn)參數(shù), 輸入正確的天線高1 編輯sestbl. 文件選擇計(jì)算方按案和設(shè)置與計(jì)算方按案相應(yīng)的參數(shù)；2 編輯sittbl.文件，對(duì)每一個(gè)站設(shè)置先驗(yàn)約束；3 編輯 文件，以求每一個(gè)站有正確的天線高;4 把上述所有文件軟連接到工作目錄：。進(jìn)入工作目錄，比如，cd 2003322 , 用linksday year doy project name, linksday 2003 322 lake. 工程名（4個(gè)字符）在 文件中第一行，后4個(gè)字符是軌道碼。64, 按如下計(jì)算流程計(jì)算 rinex o-filesrinex n-file t-fileg-filey-filesh_sp3fitigs sp3 or sp1 fileut1., pole., soltab.,luntab., nuntab. t-filesh_check_j-filerinex n-filemakejSsvs_j-filesh_check_sessl-file, j-file, rinex o-file,rinex n-filex-filek-file makexi-fileb*.bat fileb*ddd filesd-file, t-file, j-file,x or c-files, sestbl., sittbl. fixdrvq-file,o-file, h-file c-filenew l-file, new g-file etc.csh b*.bat75, 批處理模塊b*.bat的計(jì)算流程ut1., pole., leap.sec,soltab., luntab.,nuntab., g-filet-filearcy?t.dddyawtabt-file, y?.dddl-file, x or c-files, i-file,j-file, t-file, rcvant.dat,svnav.dat, gdetic.dat,antmod.dat,grid.oct, stations.octP-filesc-files,modelsh_auteditAutcln.cmd.prefitautcln.cmdautcln.sumautcln.outc-filesautcln.cmd.prefitc-files, d-fileautclncfmrgm-filec-filesl- , g- , q- , o- , h- , t-fileetc.solvec-files, m-filerun scandd do civewIf rerun autcln ?nrms 0.4 ? Choice 3 yes no choice 2 choice 1 no yesAdjust parameters in autcln.cmd StopDo one-way by civew See next pagesh_auteditAutcln.cmd.postfitn-fileautcln.sumautcln.cmd autclnc-filesc-filesAutcln.cmd.postfitc-file , n-filec-filescfmrgm-filet-, g-, l-, q-,o-, h-file,etc.c-files, m-filesolvevscan.outv?1.ddd.worstscan.ddm-file,c-filesscandd8) makexp，sh_sp3fit , sh_check_sess , makej , makex , fixdrv 等模塊的運(yùn)行。A， makexp 的運(yùn)行。目的： 產(chǎn)生 文件； 為makej 和 makex產(chǎn)生輸入文件(1) sh_makexp -yr -doy -nav -sinfo (2) makexp makexp.inpyeardoysessionauto2350.00n60 00 00 1440makexp.inpB， sh_sp3fit的運(yùn)行 目的： 使用偽觀測(cè)（從觀測(cè)所得位置和速度）獲得精確的衛(wèi)星初始狀態(tài)。Sh_sp3fit f -d -o -I -r -t 產(chǎn)生歷表文件的其它方法 bccheck （產(chǎn)生新的導(dǎo)航文件 ） sh_bcfit bctot.inpC, sh_check_sess 的運(yùn)行Sh_check_sess -sess 278 -type gfile -file Sh_check_sess -sess 278 -type jfile -file D, makej 的運(yùn)行 makej auto2280.00n jpgga0.228 makej E, makex的運(yùn)行， makex makex F, fixdrv的運(yùn)行， fixdrv fixdrv四, GAMIT 文件介紹 （見說明書第3章）1， 與站名有關(guān)的文件（x-,c-,k-,p-,v-,w-,z-, etc.）2， 與解有關(guān)的文件 (d-,autcln.out,autcln.sum,a-, b-,g-,h-,I-, j-, m-, n-,o-,q-,t-,u-,v-,y- etc.)3, 與特定的實(shí)驗(yàn)有關(guān)的文件(,,l-, etc.)4, 全局文件(gdetic.dat, svnav.dat, antmod.dat, rcvant.dat, ut1., pole., leap.sec, nutabl., luntab., soltab., stations.oct, grid.oct. etc.)五， 所需文件的準(zhǔn)備。1， L 文件的準(zhǔn)備1) 由 rinex 觀測(cè)文件準(zhǔn)備L文件ls ./rinex/?/*o .grep APPRO *o xyz.rnxrx2apr xyz.rnx cut c5- xyz.rnx.apr tmpcoordinate ratio 這里 C1 是這個(gè)周跳參數(shù)在實(shí)驗(yàn)取到整數(shù)值的過程中所獲得的最小chi*2; C2 是下一個(gè)最小chi*2, min 是該命令設(shè)置的 min chi*2, 而 ratio 是命令中的 ratio chi*2. 該命令的一般形式為： remove_bias_cond chi*2 ratio min chi*2 max gap 實(shí)際的GAP 必須小于命令中所給的 max gap . 該命令的例子為： remove_bias_cond 10.0 3.0 1800 12 3.0 3600 該命令的最后一個(gè)參數(shù) 是可選擇的。假如它被選擇，則由前式所算得的 X 須減去 1 large gap scale*atan(gap/min data)， 其中，gap是在數(shù)據(jù)里的 gap, min data 是空隙兩邊資料數(shù)較小的那段資料的資料數(shù)。Latge gap scale的省缺值是 5.0 3, phs_fit_tol 該命令設(shè)置判別在用預(yù)擬合的鐘擬合非差相位資料時(shí)是否出現(xiàn)周跳的判據(jù)。其參數(shù)的物理意義及其省缺值為：（1） 第一次通通過時(shí)平均相位殘差相對(duì)于距離解的偏離，缺省值為1000（2） 第一次通通過時(shí)最壞相位殘差相對(duì)于距離解的偏離，缺省值為 500（3） 第二次通通過時(shí)平均相位殘差相對(duì)于距離解的偏離，缺省值為 200（4） 第二次通通過時(shí)最壞相位殘差相對(duì)于距離解的偏離，缺省值為 100第3與第4個(gè)值通常要比前兩個(gè)值小很多。實(shí)際插標(biāo)處的偏離通過設(shè)置 status_report, pass1_slips pass2_slips 看到。 該命令的一個(gè)例子如下： phs_fit_tol 5000 4999 500 499 以上3條命令最為重要，使用 autcln 時(shí)應(yīng)特別注意。 其它重要的命令有：4, rng_jump_tol n-allan sd min jump (sec) 這個(gè)命令給出在偽距的 OC 中能被看成周跳的不連續(xù)性的大小。 其第1個(gè)值 n-allan sd 是由 ALLAN 標(biāo)準(zhǔn)偏差所給出的總的穩(wěn)定性的 multiplier; 其第2 個(gè)值，以微秒為單位，是將被考慮為周跳的最小值。那就是說， 當(dāng)這兩個(gè)值都被超過時(shí)，一個(gè)鐘的跳躍出現(xiàn)了。這兩個(gè)參數(shù)的省缺值是 100.0 和 0.95， 但對(duì)于較壞的資料可使用 20 和 0.1。5， clk_reset_tol jump difference (sec) 這個(gè)命令給出鐘重新設(shè)置一個(gè)毫秒的跳躍的容忍程度。 一般設(shè)置為10 sec。 但是假如先驗(yàn)鐘的多項(xiàng)式（由 FIXDRV）與資料匹配得不是很好（由于壞的先驗(yàn) 坐標(biāo)或衛(wèi)星軌道），這個(gè)值需要增加。對(duì)于質(zhì)量低的很多石英鐘，該值可設(shè)為 100sec。6， ion_jump_tol receiver code/all max gap (sec) multiplier min dion (cyc) max dion (cyl) 這個(gè)命令設(shè)置電離層延遲觀測(cè)量 LG 跳躍的容忍程度。也就是 LG 跳躍的可容 忍程度在 min dion 和 max dion之間. 超出這個(gè)范圍就是周跳。這個(gè)命令的一個(gè) 例子為： ion_jump_tol all 30.0 6.0 2.0 5.0 這樣的設(shè)置即使對(duì)于極區(qū)和赤道地區(qū)的接收機(jī)也非常適合。對(duì)于老的接收機(jī)和電 離層比較平靜的地區(qū)，這個(gè)命令應(yīng)如下設(shè)置 ion_jump_tol shao 240.0 4.0 0.8 5.07， allow_one_bg yes/no. 當(dāng)一個(gè)接收機(jī)在同一時(shí)刻對(duì)所有衛(wèi)星都出現(xiàn)周跳時(shí)，如 設(shè)置 yes， 則允許 autcln 選擇一個(gè)衛(wèi)星的信號(hào)為基準(zhǔn)，而對(duì)其它衛(wèi)星的信號(hào)進(jìn) 行修改。一般應(yīng)如此設(shè)置。8， remove_first_bia yes/no。autcln 總是在每一個(gè)非差序列的開始插標(biāo)，這對(duì) solve 沒有影響。但如設(shè)置yes，在生成 c文件時(shí) 這個(gè)額外的標(biāo)志將被移去。而生成單 差文件則不能去掉這個(gè)標(biāo)志，因此當(dāng)要求生成單差文件時(shí)，該命令須設(shè)置為 no。9， edit_site_sv site code PRN start epock stop epock 該命令設(shè)置一個(gè)特定的歷元范圍 start epock stop epock，在這個(gè)范圍內(nèi)的某 個(gè)站 site的某個(gè)衛(wèi)星 PRN的資料不能用于決定鐘的性能，也不能用于形成雙差。 這個(gè)命令對(duì)于去掉壞衛(wèi)星資料是很有用的。 all對(duì)于特定的衛(wèi)星可刪去所有的 站的資料；0 可刪去對(duì)于特定的站的所指歷元范圍內(nèi)的所有衛(wèi)星的資料。10， min_elevation min clean elevation (deg) min output elev (deg) 該命令設(shè)置資料取舍的最小仰角，小于該仰角的資料將不會(huì)被使用。最小輸出仰 角是指在生成 c文件時(shí)的資料取舍的最小仰角。這是可選的。 假如設(shè)置 use_gamit_elev = yes， 則這個(gè)命令被忽略。11， trim_oneway_tol min_dtl_bias min-good-bias min_dtr_end min_good_end min_dtl_bias 是指兩個(gè)標(biāo)志之間要求的最小時(shí)間間隔 min-good-bias 在非差資料中兩個(gè)標(biāo)志之間的標(biāo)志被移去的這兩個(gè)標(biāo)志之間要 求的最小時(shí)間間隔 min_dtr_end 最后一個(gè)標(biāo)志后的資料數(shù)在總的資料數(shù)中的北例 min_good_end 最后一個(gè)標(biāo)志后所要求的最小資料數(shù) 該命令是用來刪去小段非差資料的。12， dd_return-size max WL max LC max LG 該命令用來設(shè)置修改周跳所用的資料數(shù)。 max WL 用非差資料中的 WL 估計(jì)的最大資料數(shù)，省缺值是100； max LC 用雙差資料中的 LC 估計(jì)的最大資料數(shù)，省缺值是50； max LG 用雙差資料中的 LG 估計(jì)的最大資料數(shù)，省缺值是10。 假如資料曲線較為彎曲，比如對(duì)于電離層活動(dòng)明顯地區(qū)的LG曲線，這些值應(yīng) 該被減少。但他們不能小于5。是可選擇的。 它設(shè)置一個(gè)時(shí)間 間隔。在這個(gè)時(shí)間間隔上，非差的L1、L2 距離資料將被形成 WL 和LG 進(jìn) 行批處理。如按省缺值，則該命令為： dd_return-size 100 50 10 它的意思是在修改周跳時(shí)，對(duì)于30秒采樣間隔、非差的WL及雙差的LC 和 LG ，該周跳的每一邊應(yīng)分別有50、25、5分鐘的資料。者可以減少這些值2 4倍。 當(dāng)在 sestbl. 文件中設(shè)置 autcln postfit = R or Y; autcln command file = autcln.cmd ; use N-file = Yes ; 則在 autcln.cmd中應(yīng)包含 : 13， apply_phs_clk max iter non-int converged% over shoot 這個(gè)命令的作用是產(chǎn)生去除了星鐘和站鐘的影響后的非差(one-way)相位殘 差。這是通過迭代完成的。max iter給出了最大迭代次數(shù)（省缺值是30）； 而non-int給出了在估計(jì)鐘的零點(diǎn)差時(shí)非整數(shù)標(biāo)志(bias)參數(shù)被使用之前的迭代 次數(shù)（省缺值是3）；converged% 給出了非差相位殘差的RMS在兩次迭代 之間的相對(duì)其平均值的變化的百分比（省缺值是0.1）。小于這個(gè)百分比， 則迭代已經(jīng)收斂；over shoot 則用來加快收斂進(jìn)程，其省缺值是1.5。當(dāng)企圖 尋找問題時(shí)， 這條命令是有用的 。 14， use_posifit 。 該命令啟動(dòng)后擬合編輯。其所需之參數(shù)全部或部 分來自被 solve更新了的 m文件.15, postfit_edit start iter nsigma max restore max rms 該命令允許使用 n-sigma判據(jù)編輯相位殘差。 當(dāng)這個(gè)命令啟動(dòng)時(shí)， apply_phs_clk命令也必然啟動(dòng)。start iter 表明在編輯開始之前alliy_phs_clk 的迭代次數(shù)。 若干次迭代應(yīng)該在編輯之前完成，一般是9次。這個(gè)nsigma 是每一個(gè)站的 LC 相位資料的RMS 的彌散度。 其省缺值是4.0。 殘差小于 nsigma，且由于標(biāo)志(bias)靠近而被插標(biāo)的 LC 的殘差，如其小于max restore ， 則它們能被重新恢復(fù)過來。 如果一個(gè)站的RMS大于 max rms，它的資料將 被全部剔掉。 16, pf_remove_bf yes/no 該命令允許去掉后擬合編輯后的非差資料中的周跳標(biāo)志 17， summary autcln.post.sum 其它較為重要的命令有： 18， use_gamit no 該命令告訴 autcln 忽略由接收機(jī)所插入的表示失鎖的任 何標(biāo)志。19， max_scan_edit 該命令設(shè)置在 autcln決定在刪除所有的資料之前對(duì)特定的站 和衛(wèi)星的資料所能加的標(biāo)志數(shù).20, remove_ms-jump yes/no 該命令決定當(dāng)生成C 文件時(shí) 是否要把鐘中的毫秒 的跳去掉。21， scan_site all/none/list of site 讓使用者指明哪些站在雙差形成之前應(yīng)該被檢 查一遍。 假如在雙差形成之前發(fā)現(xiàn)未插標(biāo)的周跳，這個(gè)命令應(yīng)該被使用。 設(shè) 置 ALL， 則所有的站都要被檢查， NONE將沒有站被查；站名前的負(fù)號(hào)