利用execl的数据库进行仿真

Shark丶

仿真数据：来自于execl数据库中的空间坐标（X,Y,Z，T）
仿真目的：通过execl的空间坐标数据，利用python实现类似于飞行器飞行、起降的仿真
求大神给个思路，如何将execl中的大量数据联动到python开发平台内，然后通过类似于描点的方式模拟出飞机的起降，飞行路线。然后以3D的形式展现出来。

Shark丶

类似于下面这段代码，是模拟日-地-月的运动轨迹，但是飞机的飞行路线要从execl数据库中导入，所以有些没有头绪。求大神解惑

1. import numpy as np
   
2. import matplotlib as mpl
   
3. mpl.use("TkAgg")
   
4. from matplotlib import pyplot as plt
   
5. from mpl_toolkits.mplot3d import Axes3D
   
6. import matplotlib.animation as animmation
   
7. 
   
8. r1 = 10 #地球环绕轨迹的半径
   
9. r2 = 1 #月球环绕轨迹的半径
   
10. omega1 = 2 * np.pi #地球运动的角速度
    
11. omega2 = 24 * np.pi #地球自转速度
    
12. phi = 5.1454 * np.pi / 180 #地球自转轴角
    
13. 
    
14. def update(data): #随时间t变换设置新的地球位置的坐标，也即是随animmation.FuncAnimation函数的帧数参数变化每一帧的不同坐标位置的图片对象参数
    
15.     global line1, line2 , line3
    
16.     line1.set_data([data[0], data[1]])
    
17.     line1.set_3d_properties(data[2])
    
18.     line2.set_data([data[3], data[4]])
    
19.     line2.set_3d_properties(data[5])
    
20.     line3.set_data([data[6], data[7]])
    
21.     line3.set_3d_properties(data[8])
    
22.     return line1,line2,line3,
    
23. 
    
24. def init(): #地球的起始位置
    
25.     global line1, line2, line3
    
26.     ti = 0
    
27.     t = t_drange[np.mod(ti, t_dlen)]
    
28.     xt1 = x0 + r1 * np.cos(omega1 * t)
    
29.     yt1 = y0 + r1 * np.sin(omega1 * t)
    
30.     zt1 = z0 + 0
    
31.     xt2 = xt1 + r2 * np.sin(omega2 * t)
    
32.     yt2 = yt1 + r2 * np.cos(omega2 * t)/(np.cos(phi) * (1 + np.tan(phi) ** 2))
    
33.     zt2 = zt1 + (yt2 - yt1) * np.tan(phi)
    
34.     xt21 = xt1 + r2 * np.sin(2 * np.pi * t_range)
    
35.     yt21 = yt1 + r2 * np.cos(2 * np.pi * t_range)/(np.cos(phi) * (1 + np.tan(phi) ** 2))
    
36.     zt21 = zt1 + (yt21 - yt1) * np.tan(phi)
    
37.     line1, = ax.plot([xt1], [yt1], [zt1], marker='o', color='blue',markersize=8)
    
38.     line2, = ax.plot([xt2], [yt2], [zt2], marker='o', color='orange',markersize=4)
    
39.     line3, = ax.plot(xt21, yt21, zt21, color='purple')
    
40.     return line1,line2,line3
    
41. 
    
42. def data_gen(): #随时间t变换的地球和月球坐标，也即是animmation.FuncAnimation函数的帧数参数
    
43.     global x0,y0,z0,t_dlen
    
44.     data = []
    
45.     for ti in range(1,t_dlen):
    
46.         t = t_drange[ti]
    
47.         xt1 = x0 + r1 * np.cos(omega1 * t)
    
48.         yt1 = y0 + r1 * np.sin(omega1 * t)
    
49.         zt1 = z0
    
50.         xt2 = xt1 + r2 * np.sin(omega2 * t)
    
51.         yt2 = yt1 + r2 * np.cos(omega2 * t)/(np.cos(phi) * (1 + np.tan(phi) ** 2))
    
52.         zt2 = zt1 + (yt2 - yt1) * np.tan(phi)
    
53.         xt21 = xt1 + r2 * np.sin(2 * np.pi * t_range)
    
54.         yt21 = yt1 + r2 * np.cos(2 * np.pi * t_range)/(np.cos(phi) * (1 + np.tan(phi) ** 2))
    
55.         zt21 = zt1 + (yt21 - yt1) * np.tan(phi)
    
56.         data.append([xt1, yt1, zt1, xt2, yt2, zt2, xt21, yt21, zt21])
    
57.     return data
    
58. 
    
59. t_range = np.arange(0, 1 + 0.005, 0.005) #设置环绕一周时间的范围以及运动间隔时间
    
60. t_drange = np.arange(0, 1, 0.005 )
    
61. t_len = len(t_range)
    
62. t_dlen = len(t_drange)
    
63. #sun's coordination
    
64. x0 = 0
    
65. y0 = 0
    
66. z0 = 0
    
67. #地球轨道
    
68. x1 = x0  + r1 * np.cos(omega1 * t_range)
    
69. y1 = y0 + r1 * np.sin(omega1 * t_range)
    
70. z1 = z0 + np.zeros(t_len)
    
71. 
    
72. #月球轨道
    
73. x2 = x1 + r2 * np.sin(omega2 * t_range)
    
74. y2 = y1 + r2 * np.cos(omega2 * t_range)/(np.cos(phi) * (1 + np.tan(phi) ** 2))
    
75. z2 = z1 + (y2 - y1) * np.tan(phi)
    
76. 
    
77. 
    
78. f = plt.figure(figsize=(6,6)) #绘图的画布
    
79. ax  = f.add_subplot(111,projection='3d') #设置3d坐标系
    
80. ax.set_title("Sun-Earth-Moon Model") #设置图像标题（太阳-地球-月亮转动模型）
    
81. 
    
82. ax.plot([0], [0], [0], marker='o', color= 'red', markersize=16) #绘制太阳的各种属性
    
83. ax.plot(x1, y1, z1, 'r') #绘制地球图像
    
84. ax.plot(x2, y2, z2, 'b') #绘制月球图像
    
85. ax.set_xlim([-(r1 + 2), (r1 + 2)]) #太阳用动模型在坐标系中的范围
    
86. ax.set_ylim([-(r1 + 2), (r1 + 2)]) #地球用动模型在坐标系中的范围
    
87. ax.set_zlim([-5, 5]) #月球用动模型在坐标系中的范围
    
88. 
    
89. #红色球体代表太阳，蓝色球体代表地球，橙色球体代表月亮
    
90. line1, = ax.plot([], [], [], marker='o', color='blue',markersize=8,animated = True)
    
91. line2, = ax.plot([], [], [], marker='o', color='orange',markersize=4,animated = True)
    
92. line3, = ax.plot([], [], [], color='purple',animated = True)
    
93. 
    
94. # 将上述函数对象传如animmation.FuncAnimation函数以生成连读的地球运动模型
    
95. ani = animmation.FuncAnimation(f, update, frames = data_gen(), init_func = init,interval = 20)
    
96. 
    
97. plt.show()
    

复制代码

Torwnexial

我也不太懂，但刚从知乎上看见一篇python中的excel用法：https://zhuanlan.zhihu.com/p/38875272  希望能帮到你