学问思辨行

为满足从CDDIS下载GNSS等大地测量数据的需求，开发了本软件，软件界面如图所示：软件相关说明及下载请关注公众号，后续相关文章仅在公众号发布，敬请关注：

可以登陆从CDDIS下载数据，如图：可以下载的数据类型：长文件名或短文件名观测文件；星历文件仅支持brdm；图片所示分析中心的产品及天线文件和dcb文件；下载的数据会自动解压转换为相应类型；已知bug：下载可能会进入假死状态，即出现未响应的状况，不用担心，他正在工作；用户名或者密码错误时，可能会在日志文件中提示下载某些文件失败，请检查你的登陆信息是否正确；不美观的status bar；下载文件的情况可能会跟你想的不一样：下载erp及snx，不仅仅下载周解文件，回吧选定日

1、首先新建一个空项目，在此命名为GAMP；2、将下载的GAMP压缩包解压，将GAMP -> mannual_GAMP -> GAMP_src -> Windows -> gamp_c 下的gamp_src文件夹复制到所建项目路径下GAMP -> GAMP，将测试数据也一并复制过来，GAMP -> mannual_GAMP -> Examples下的20...

转自：GPS观测值的线性组合

我们在使用消电离层组合模型的时候，可能会发现，为什么很多材料上会写着完全不同的两种消电离层组合，其实把不同的写在一起就会发现：PIF=f12f12−f22P1−f22f12−f22P2P_{IF} = \frac{f_1^2}{f_1^2-f_2^2} P_1 - \frac{f_2^2}{f_1^2-f_2^2} P_2PIF​=f12​−f22​f12​​P1​−f12​−f22​f22...

MW组合对相位观测值进行宽巷组合，即宽巷相位观测值为：ϕw=ϕ1−ϕ2\phi_w = \phi_1 - \phi_2ϕw​=ϕ1​−ϕ2​因为整周模糊度NNN与ϕ\phiϕ表示的含义相同，故：Nw=N1−N2N_w = N_1 - N_2Nw​=N1​−N2​根据ϕ=fwt\phi=f_wtϕ=fw​t，其中等式右边分别为频率和时间，可得：fwt=f1t−f2tf_w t ...

GF组合相比于MW组合探测周跳，GF更为简单。Geometry_free组合，顾名思义，消去了与几何相关的项（与频率无关项）。给出如下长度形式的相位观测方程，相关符号含义见MW组合博文107-周跳探测之MW：{L1=ρ+λ1N1−I1L2=ρ+λ2N2−f12f22I1\begin{cases}L_1 = \rho + \lambda_1 N_1 - I_1 \\L_2 = \rho +...

静态定位模式下以2019年cut02440.17o作为观测数据对比两个软件，结果如下：gamp收敛速度快一些，但rtklib收敛后似乎更稳定。

转自：JhonnyK-SP3精密星历简介  IGS精密星历采用sp3格式，其存储方式为ASCII文本文件，内容包括表头信息以及文件体，文件体中每隔15 min给出1个卫星的位置，有时还给出卫星的速度。它的特点就是提供卫星精确的轨道位置。采样率为15分钟，实际解算中可以进行精密钟差的估计或内插，以 提高其可使用的历元数。命名规则：  常用的sp3格式的命名规则为：tttwwwwd.sp3其...

  将crx2rnx.exe与*.*d文件放在同一路径下，运行cmd，进入该路径。  在cmd下，直接输入f:可直接进入f盘，然后输入cd rinex进入相应文件夹，在此为rinex文件夹。输入cd..可返回上一级目录。如图：  若转换单个d文件，可执行：crx2rnx 019b0010.19d即可得到转换后的o文件。  若批量转换，可通过循环执行：for %a in (*.*d)...

一、卫星轨道六要素1、升交点赤经Ω\OmegaΩ春分点：太阳沿黄道从天赤道以南向北通过天赤道的点，称为春分点。升交点：卫星由南向北运行时，与地球赤道面的交点。升交点赤经为地心与春分点连线 和 地心与升交点连线之间的夹角。2、轨道倾角iii卫星运动的轨道面和赤道面之间的夹角。升交点赤经和轨道倾角 确定轨道在惯性空间中的位置。3、近地点张角ω\omegaω地心与升交点之间的连线...

一、无线电波频率二、多址技术

关于RINEX文件中的N文件，摘录说明如下：格式说明采用的是Fortran 程序设计语言中的格式说明方式, 一个格式说明项通常具有如下形式: [r]fw.[m]其中:r : 重复因子, 表示后面的内容将重复的次数, 该部分是可选的;f: 数据类型符, 在RINEX 格式的说明中, 用到了如下数据类型:    X, 空格    A, 字符型    I , 整型    F, 单精度...

关于RINEX(2.10)文件中的O文件，摘录说明如下：  对于观测值，每一行表示每颗卫星的观测值，每一行表示五个观测类型，如果类型过多，则换行，因此在写程序读取时要注意这一点。...

四级移位寄存器，初始状态为1，模二相加规则：0⊕0=0;0⊕1=1;1⊕0=1;1⊕1=0;0 \oplus 0 = 0;\\0 \oplus 1 = 1;\\1 \oplus 0 = 1;\\1 \oplus 1 = 0;0⊕0=0;0⊕1=1;1⊕0=1;1⊕1=0;状态1：1,1,1,11,1,1,11,1,1,1将第三位和第四位模二相加，得到0，存至第一位，前三位移动...

1 C/A码C/A码是由10级寄存器在1.023Mbps的频率(即每秒生成1.023×1061.023 \times 10^{6}1.023×106个码元)下产生的。10级寄存器一个周期产生的码元个数为：210−1=10232^{10}-1=1023210−1=1023那么在1.023Mbps的频率下，其周期为0.001秒即1ms，码元宽度即为：0.0011023×c≈293.0522...

对于卫星定位而言，卫星发射时刻的时间归化虽然道理浅显易懂，但很容易就会陷入思维误区，博主就是活生生的例子。设定ttt是求解得到的卫星发射时刻(各种改正都已完成)，现要归化到以toet_{oe}toe​为基准的时刻，就存在一个加减604800的问题，理由如下：如图所示，对于第一种情况，toet_{oe}toe​在N+1周的开始，ttt在N周的末尾，这种条件下必有t−toe>302...

一、MATLAB相关函数介绍1、 uigetfile[fname, fpath] = uigetfile('*.*N', 'File N');功能：打开一个扩展名为“*N”的文件，获取文件名与其绝对路径； File N为弹出的对话框的名称。2、strcatfilename = strcat(fpath, fname);功能：将两个字符串连接在一起，这样就获取了N文件的路径。3、f...

涉及到时间系统时，总是被搞得眼花缭乱，看到一篇文档介绍的挺好，在此做了一部分摘录，具体请参考：空间大地测量（二）-陈义-豆丁空间大地测量（二）-陈义-百度摘录如下：...