模拟太阳系8大行星运行图matplotlib实现

　　1效果图
　　2说明
　　2。1我曾用matplotlib实现模拟太阳地球月亮的运行轨迹图，很受欢迎。
　　《python3的matplotlib的模拟太阳地球月亮运动示意图代码分析》
　　2。2这次再使用python3的语言，用matplotlib来实现太阳系的八大行星模拟示意图。
　　2。3熟悉python编程语言和思维，熟悉matplotlib作图。
　　3代码讲解：
　　3。1第1步：导入模块importnumpyasnpfrommatplotlibimportpyplotaspltfrommpltoolkits。mplot3dimportAxes3Dimportmatplotlib。animationasanimmation
　　3。2第2步：画布的相关设置fplt。figure（figsize（15，15））绘图的画布大小和定义初始化f。patch。setfacecolor（black）画布背景颜色设置axf。addsubplot（111，projection3d）设置3d坐标系这个只是图表中的背景颜色设置ax。setfacecolor（black）ORax。set（facecolorblack）三维坐标刻度设置ax。setxlim（〔65，65〕）ax。setylim（〔65，65〕）模型在坐标系中的范围ax。setzlim（〔3，3〕）
　　3。3第3步：太阳设置：形状，颜色和大小，位置固定ax。plot（〔0〕，〔0〕，〔0〕，markero，colorred，markersize70）
　　3。4第4步：8大行星的半径、轨道和角度设置常量的设置设置轨道半径rlist〔40，43，46，51，55，62，71，85〕设置行星角速度olist〔13，9，7，6，5，4，3，2〕设置环绕一周时间的范围以及运动间隔时间trangenp。arange（0，10。005，0。005）tlenlen（trange）行星坐标随时间变化方程z坐标zznp。zeros（tlen）x和y坐标xllist〔〕yllist〔〕foriinrange（len（rlist））：xllist。append（rlist〔i〕np。cos（olist〔i〕np。pitrange））yllist。append（rlist〔i〕np。sin（olist〔i〕np。pitrange））绘制行星的圆形轨迹foriinrange（len（rlist））：轨道线颜色和宽度设置ax。plot（xllist〔i〕，yllist〔i〕，zz，gray，linewidth0。2）
　　3。5第5步：8大行星颜色和大小设置ball1，ax。plot（〔〕，〔〕，〔〕，markero，coloryellow，markersize10）ball2，ax。plot（〔〕，〔〕，〔〕，markero，colorpink，markersize10）ball3，ax。plot（〔〕，〔〕，〔〕，markero，colorblue，markersize10）ball4，ax。plot（〔〕，〔〕，〔〕，markero，colororange，markersize10）ball5，ax。plot（〔〕，〔〕，〔〕，markero，colorbrown，markersize16）ball6，ax。plot（〔〕，〔〕，〔〕，markero，colorred，markersize16）ball7，ax。plot（〔〕，〔〕，〔〕，markero，colorgray，markersize14）ball8，ax。plot（〔〕，〔〕，〔〕，markero，colorgreen，markersize14）
　　3。6第6步：初始化行星的起始位置definit（）：returnball1，ball2，ball3，ball4，ball5，ball6，ball7，ball8
　　3。7第7步：数据更新，随时间t的变换设置新的行星位置坐标defupdate（data）：ball1。setdata（〔data〔0〕，data〔1〕〕）ball1。set3dproperties（data〔2〕）ball2。setdata（〔data〔3〕，data〔4〕〕）ball2。set3dproperties（data〔5〕）ball3。setdata（〔data〔6〕，data〔7〕〕）ball3。set3dproperties（data〔8〕）ball4。setdata（〔data〔9〕，data〔10〕〕）ball4。set3dproperties（data〔11〕）ball5。setdata（〔data〔12〕，data〔13〕〕）ball5。set3dproperties（data〔14〕）ball6。setdata（〔data〔15〕，data〔16〕〕）ball6。set3dproperties（data〔17〕）ball7。setdata（〔data〔18〕，data〔19〕〕）ball7。set3dproperties（data〔20〕）ball8。setdata（〔data〔21〕，data〔22〕〕）ball8。set3dproperties（data〔23〕）returnball1，ball2，ball3，ball4，ball5，ball6，ball7，ball8
　　3。8第8步：8大行星的数据更新，随时间t变换的行星坐标defdatagen（）：data〔〕fortiinrange（1，tlen）：ttrange〔ti〕行星与太阳坐标的关系方程xt1rlist〔0〕np。cos（olist〔0〕np。pit）yt1rlist〔0〕np。sin（olist〔0〕np。pit）xt2rlist〔1〕np。cos（olist〔1〕np。pit）yt2rlist〔1〕np。sin（olist〔1〕np。pit）xt3rlist〔2〕np。cos（olist〔2〕np。pit）yt3rlist〔2〕np。sin（olist〔2〕np。pit）xt4rlist〔3〕np。cos（olist〔3〕np。pit）yt4rlist〔3〕np。sin（olist〔3〕np。pit）xt5rlist〔4〕np。cos（olist〔4〕np。pit）yt5rlist〔4〕np。sin（olist〔4〕np。pit）xt6rlist〔5〕np。cos（olist〔5〕np。pit）yt6rlist〔5〕np。sin（olist〔5〕np。pit）xt7rlist〔6〕np。cos（olist〔6〕np。pit）yt7rlist〔6〕np。sin（olist〔6〕np。pit）xt8rlist〔7〕np。cos（olist〔7〕np。pit）yt8rlist〔7〕np。sin（olist〔7〕np。pit）zt0data。append（〔xt1，yt1，zt，xt2，yt2，zt，xt3，yt3，zt，xt4，yt4，zt，xt5，yt5，zt，xt6，yt6，zt，xt7，yt7，zt，xt8，yt8，zt〕）returndata
　　3。9第9步：动画挂起和展示anianimmation。FuncAnimation（f，update，framesdatagen（），initfuncinit，interval200）隐藏画轴plt。axis（off）plt。show（）
　　分享出来，在动画中学习计算机编程，学习python思维。