CATS(Create and Analyse Time Series):cats is a program to use Maximum Likelihood Estimation to fit a multi-parameter model to a time series (such as continuous GPS). The routine solves for all parameters simultaneously, but in two parts to increase speed; the linear part includes an offset and slope, the possibility of abrupt steps (as from earthquakes) and sinusoidal terms (for example annual and semiannual terms), while the non-linear part solves for several specific noise models and combinations.
在使用CATS软件做CORS站坐标时间序列的噪声分析时,想按照文献里的
这样的方法分析自己的数据,在测站较多,噪声模型组合也多的情况下,通过脚本实现循环,可以不用频繁手动操作。贴出脚本如下:
脚本里只有前五个噪声模型的组合方式,后面的也可以按照这个格式加进去。得到的其中一个结算结果如下:
通过执行上述脚本,可以实现多个测站多种噪声模型组合的自动运算,脚本运行成功后通过结果文件的相应的数据提取实现数据分析的目的。
[未完待续…后面等运算完贴上提取数据的脚本]
2019.10.8日更新
在运用CATS软件对各个测站进行噪声模型解算过后,需要对解文件中的数据进行提取,上述脚本中对各个测站的各个模型解算过后生成的解文件命名均是按照一定的规律来的,这给编写脚本提取解文件中的信息提供了极大的帮助,贴出提取时间序列的各个参数的运行脚本:
下面的脚本可以用来提取各个噪声模型解算结果文件中的最大似然估计值:
【注】在刚开始调试中,printf一行最后打的定向输出符号用的是“>”,导致结果文件*.dat中全部只有一行。后来通过查找资料发现,定向输出符号有“>”和“>>”之分,具体参考:https://blog.csdn.net/tq384998430/article/details/72820943
简单来说就是前者每条输出都会删除文件中原有的数据,而后者不会。这也是我的结果文件中始终只有一条数据的原因。
修改后生成的结果文件分别如下图:
以便于后面导入excel或者其他软件中进行分析。