資訊專欄INFORMATION COLUMN
摘要：用畫一個平面的太陽系得到一些朋友的欣賞，然后有同學提出了繪制三維太陽系的要求。從畫圖的角度來說，三維太陽系其實并不難，問題在于八大行星對黃道面的傾斜太小，所以盡管畫個三維的圖，但就觀感而言，無非是把二維的嵌入到三維空間罷了。
用Python畫一個平面的太陽系得到一些朋友的欣賞，然后有同學提出了繪制三維太陽系的要求。
從Python畫圖的角度來說，三維太陽系其實并不難，問題在于八大行星對黃道面的傾斜太小，所以盡管畫個三維的圖，但就觀感而言，無非是把二維的嵌入到三維空間罷了。
來點小行星
代碼如下
from os import cpu_countimport numpy as npfrom numpy.random import randimport matplotlib.pyplot as pltfrom matplotlib import animationau,G,RE,ME = 1.48e11,6.67e-11,1.48e11,5.965e24m = np.array([3.32e5,0.055,0.815,1,0.107,317.8])*ME*Gr = np.array([0,0.387,0.723,1,1.524,5.203])*REv = np.array([0,47.89,35.03,29.79,24.13,13.06])*1000theta = rand(len(m))*np.pi*2cTheta,sTheta = np.cos(theta), np.sin(theta)xyz = r*np.array([cTheta, sTheta, 0*r]) #位置三分量，因為參數太多，所以把這三個分量寫在了一起uvw = v*np.array([-sTheta, cTheta, 0*v]) #速度三分量N_ast = 100m_ast = rand(N_ast)*1e20r_ast = (rand(N_ast)*3.5+1.6)*REv_ast = np.sqrt(G*3.32e5*ME/r_ast) #小行星速度sqrt(GM/R)theta = rand(N_ast)*np.pi*2phi = (rand(N_ast)-0.5)*0.3 #給一個隨機的小傾角cTheta,sTheta = np.cos(theta), np.sin(theta)cPhi,sPhi = np.cos(phi),np.sin(phi)xyza = r_ast*np.array([cTheta*cPhi, sTheta*cPhi, sPhi])uvwa = v_ast*np.array([-sTheta*cPhi, cTheta*cPhi, sPhi])name = "solar.gif"fig = plt.figure(figsize=(10,10))ax = fig.add_subplot(projection="3d")ax.grid()ax.set_xlim3d([-5.5*RE,5.5*RE])ax.set_ylim3d([-5.5*RE,5.5*RE])ax.set_zlim3d([-5.5*RE,5.5*RE])traces = [ax.plot([],[],[],"-", lw=0.5)[0] for _ in range(len(m))]pts = [ax.plot([],[],[],marker="o")[0] for _ in range(len(m))]pt_asts = [ax.plot([],[],[],marker=".")[0] for _ in range(N_ast)]N = 500dt = 3600*50ts = np.arange(0,N*dt,dt)xyzs,xyzas = [],[]for _ in ts: xyz_ij = (xyz.reshape(3,1,len(m))-xyz.reshape(3,len(m),1)) r_ij = np.sqrt(np.sum(xyz_ij**2,0)) xyza_ij = (xyz.reshape(3,1,len(m))-xyza.reshape(3,N_ast,1)) ra_ij = np.sqrt(np.sum(xyza_ij**2,0)) for j in range(len(m)): for i in range(len(m)): if i!=j : uvw[:,i] += m[j]*xyz_ij[:,i,j]*dt/r_ij[i,j]**3 for i in range(N_ast): uvwa[:,i] += m[j]*xyza_ij[:,i,j]*dt/ra_ij[i,j]**3 xyz += uvw*dt xyza += uvwa*dt xyzs.append(xyz.tolist()) xyzas.append(xyza.tolist())xyzs = np.array(xyzs).transpose(2,1,0)xyzas = np.array(xyzas).transpose(2,1,0)def animate(n): for i in range(len(m)): xyz = xyzs[i] traces[i].set_data(xyz[0,:n],xyz[1,:n]) traces[i].set_3d_properties(xyz[2,:n]) pts[i].set_data(xyz[0,n],xyz[1,n]) pts[i].set_3d_properties(xyz[2,n]) for i in range(N_ast): pt_asts[i].set_data(xyzas[i,0,n],xyzas[i,1,n]) pt_asts[i].set_3d_properties(xyzas[i,2,n]) return traces+pts+pt_astsani = animation.FuncAnimation(fig, animate, range(N), interval=10, blit=True)plt.show()ani.save(name)文章版權歸作者所有，未經允許請勿轉載,若此文章存在違規行為，您可以聯系管理員刪除。
轉載請注明本文地址：http://m.specialneedsforspecialkids.com/yun/122575.html
文章目錄 日地月三體日地火太陽系 你們要的3D太陽系 圖片上傳之后不知為何幀率降低了許多。。。 日地月三體 所謂三體，就是三個物體在重力作用下的運動。由于三點共面，所以三個質點僅在重力作用下的運動軌跡也必然無法逃離平面。 三體運動所遵循的規律就是古老而經典的萬有引力 F ...
摘要：看看這些驚人的純實驗，也許你自己也可以嘗試一下。項目鏈接疊疊高游戲你可以不用來編寫一個游戲。這個純粹用實現的疊疊高游戲看上去很簡單，但是很有趣，而且圖形也很漂亮。項目鏈接鬼影漸變效果按鈕令人驚訝的是它是只用編寫的。 翻譯：瘋狂的技術宅原文：https://1stwebdesigner.com/12... 本文首發微信公眾號：jingchengyideng歡迎關注，每天都給你推送新鮮...
摘要：看看這些驚人的純實驗，也許你自己也可以嘗試一下。項目鏈接疊疊高游戲你可以不用來編寫一個游戲。這個純粹用實現的疊疊高游戲看上去很簡單，但是很有趣，而且圖形也很漂亮。項目鏈接鬼影漸變效果按鈕令人驚訝的是它是只用編寫的。 翻譯：瘋狂的技術宅原文：https://1stwebdesigner.com/12... 本文首發微信公眾號：jingchengyideng歡迎關注，每天都給你推送新鮮...
摘要：看看這些驚人的純實驗，也許你自己也可以嘗試一下。項目鏈接疊疊高游戲你可以不用來編寫一個游戲。這個純粹用實現的疊疊高游戲看上去很簡單，但是很有趣，而且圖形也很漂亮。項目鏈接鬼影漸變效果按鈕令人驚訝的是它是只用編寫的。 翻譯：瘋狂的技術宅原文：https://1stwebdesigner.com/12... 本文首發微信公眾號：jingchengyideng歡迎關注，每天都給你推送新鮮...
摘要：紋理貼圖的應用以及實現一個太陽系的自轉公轉點擊查看演示地址中的紋理紋理貼圖是通過將圖像應用到對象的一個或多個面，來為對象添加細節的一種方法。 紋理貼圖的應用以及實現一個太陽系的自轉公轉 點擊查看demo演示 demo地址：https://nsytsqdtn.github.io/d... three.js中的紋理 紋理貼圖是通過將圖像應用到對象的一個或多個面，來為3D對象添加細節的一種...
閱讀 1040·2021-11-22 13:53
閱讀 1594·2021-11-17 09:33
閱讀 2396·2021-10-14 09:43
閱讀 2855·2021-09-01 11:41
閱讀 2275·2021-09-01 10:44
閱讀 2916·2021-08-31 09:39
閱讀 1453·2019-08-30 15:44
閱讀 1864·2019-08-30 13:02
