使用matlab代碼嘗試繪制脈沖星周期躍變的數(shù)據(jù)圖

獲得自轉(zhuǎn)模型文件.par和到達時間文件.tim

使用pulsar software軟件tempo2畫圖，命令為：

在圖像中擬合得到周期躍變現(xiàn)象，使用Ctrl+J保存得到glitch1.txt為周期躍變的殘差圖數(shù)據(jù)

使用Ctrl+R分段選擇躍變前后數(shù)據(jù)，得到regions.dat數(shù)據(jù)，之后使用tempo glitch插件作圖

得到result.dat即為g1.dat以及插件繪制的圖像g1.ps

可以自己使用matlab根據(jù)數(shù)據(jù)手動繪制圖像，需要的數(shù)據(jù)是躍變歷元glitch epoch和殘差數(shù)據(jù)glitch1.txt以及F0，F(xiàn)1值g1.dat，采用如下代碼：

這個代碼一次繪制了三幅圖：殘差圖，ΔF0圖和F1變化圖。

畫圖時要注意：

[1] 殘差數(shù)據(jù)和F0數(shù)據(jù)范圍不要差太多

[2] F1演化這部分，線性擬合需要看是否有線性趨勢，如果沒有就不用擬合

[3] 大躍變，無法擬合躍變后數(shù)據(jù)，不需要畫殘差圖