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