考虑空气阻力的自由落体运动python画图简单的

好的，以下是一个考虑空气阻力的自由落体运动的 Python 代码示例，使用 Matplotlib 库绘制：

import matplotlib.pyplot as plt
import numpy as np

# 定义常量
g = 9.8  # 重力加速度
c = 0.1  # 空气阻力系数
m = 1  # 物体质量
t_max = 10  # 模拟时间
dt = 0.01  # 时间步长

# 定义初始条件
v = 0  # 初始速度为0
y = 100  # 初始高度为100

# 定义数组保存计算结果
t_array = np.arange(0, t_max, dt)
v_array = np.zeros_like(t_array)
y_array = np.zeros_like(t_array)

# 计算运动轨迹
for i in range(len(t_array)):
    v_array[i] = v
    y_array[i] = y
    a = g - c/m*v**2  # 计算加速度，考虑空气阻力
    v += a*dt  # 计算速度
    y -= v*dt  # 计算高度

# 绘制图像
plt.plot(t_array, y_array)
plt.xlabel('Time (s)')
plt.ylabel('Height (m)')
plt.title('Projectile Motion with Air Resistance')
plt.show()

运行这段代码，将会得到一个高度随时间变化的图像，表示考虑空气阻力的自由落体运动的运动轨迹。