GAMIT简介pjh.doc

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将没有站被查；站名前的负号