[CATS软件学习笔记]CATS分析多个测站多种噪声模型组合的脚本

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或者其他软件中进行分析。

评论 15
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值