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