基本代码
使用在LINUX环境下的Mathematica,注意路径
- 载入数据(依次,设置目录,载入数据文件并选择数据范围)
SetDirectory["/home/yiwen/Desktop/Data"];
data6 = Import["20171218-6.data", "Table"][[30 ;; 800, {38, 39}]];
data05 = Import["20171218-11.data", "Table"][[30 ;; 400, {38, 39}]];
Tips:MESA的数据是data文件,可以用Excel或者WPS打开,然后点数据、分列、空格符号(空格)来分列,检查数据,这样比较清晰。但是不要保存,用Excel和WPS保存文件会造成精度的损失。(你会画出锯齿状的图的!)
- 画图
HW1 = ListLinePlot[{data6, data05}, ScalingFunctions -> {"Reverse"},
PlotStyle -> {Red, Blue}, PlotLegends -> {"6M", "11M"},
Frame -> True,
FrameLabel -> {{HoldForm[log_L], None}, {HoldForm[log_Teff], None}}]
- ListLinePlot[] 画连线图,ListPlot 画点图,可以检查点的疏密
- ScalingFunctions -> {"Reverse"},x轴反向,因为H-R图的温度轴是反的。ScalingFunctions -> {"Reverse","Reverse"} x,y轴全反向
- PlotStyle -> {Red, Blue} 线的颜色
- PlotLegends -> {"6M", "11M"} 图例
- Frame -> True 边框
- FrameLabel -> {{HoldForm[log_L], None}, {HoldForm[log_Teff], None}} 边框上的标题(当坐标轴名)
- 输出(保存为jpg格式,命名为HW1)
Export["H-R.jpg", HW1]
##其它图
data62 = Import["20171218-6.data", "Table"][[30 ;; 400, {46, 47}]];
data112 = Import["20171218-11.data", "Table"][[7 ;; 400, {46, 47}]];
HW2 = ListLinePlot[{data62, data112}, PlotStyle -> {Red, Blue},
PlotLegends -> {"6M", "11M"}, Frame -> True,
FrameLabel -> {{HoldForm[log_Teff], None}, {HoldForm[log_Rho],
None}}]
Export["rho-Teff.jpg", HW2]
data6cT = Import["20171218-6.data", "Table"][[30 ;; 800, {3, 53}]];
data6oT = Import["20171218-6.data", "Table"][[30 ;; 800, {3, 54}]];
data11cT = Import["20171218-11.data", "Table"][[30 ;; 490, {3, 53}]];
data11oT = Import["20171218-11.data", "Table"][[30 ;; 490, {3, 54}]];
HW3 = ListLinePlot[{data6cT, data6oT, data11cT, data11oT},
PlotStyle -> {Red, Blue, Green, Yellow},
PlotLegends -> {"C-6M", "O-6M", "C-11M", "O-11M"}, Frame -> True,
FrameLabel -> {{HoldForm[c12 & o16], None}, {HoldForm[T], None}}]
Export["fengdu.jpg", HW3]
一些资源
https://www.zhihu.com/question/40525831
mathematica Linux、Win、Mac版及破解网页下载
链接:https://pan.baidu.com/s/1slmguDZ 密码:ijo1