matlab代码尝试绘制脉冲星周期跃变数据图

使用matlab代码尝试绘制脉冲星周期跃变的数据图

获得自转模型文件.par和到达时间文件.tim

使用pulsar software软件tempo2画图，命令为：

在图像中拟合得到周期跃变现象，使用Ctrl+J保存得到glitch1.txt为周期跃变的残差图数据

使用Ctrl+R分段选择跃变前后数据，得到regions.dat数据，之后使用tempo glitch插件作图

得到result.dat即为g1.dat以及插件绘制的图像g1.ps

可以自己使用matlab根据数据手动绘制图像，需要的数据是跃变历元glitch epoch和残差数据glitch1.txt以及F0，F1值g1.dat，采用如下代码：

这个代码一次绘制了三幅图：残差图，ΔF0图和F1变化图。

画图时要注意：

[1] 残差数据和F0数据范围不要差太多

[2] F1演化这部分，线性拟合需要看是否有线性趋势，如果没有就不用拟合

[3] 大跃变，无法拟合跃变后数据，不需要画残差图