mathematica做恒星内部结构作业

基本代码

使用在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}}]

1. ListLinePlot[] 画连线图，ListPlot 画点图，可以检查点的疏密
2. ScalingFunctions -> {"Reverse"}，x轴反向，因为H-R图的温度轴是反的。ScalingFunctions -> {"Reverse"，"Reverse"} x,y轴全反向
3. PlotStyle -> {Red, Blue} 线的颜色
4. PlotLegends -> {"6M", "11M"} 图例
5. Frame -> True 边框
6. 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