### Python 实现自由落体模拟 自由落体是指物体仅受重力作用而下落的过程。在理想情况下，忽略空气阻力时，可以通过牛顿第二定律来描述运动过程[^2]。然而，在实际应用中，空气阻力会对物体的加速度产生显著影响。 #### 落地时间计算 假设初始高度为 \( h_0 \)，初速度为零，则不考虑空气阻力的情况下，落地时间可以使用以下公式计算： \[ t = \sqrt{\frac{2h_0}{g}} \] 其中 \( g \) 是重力加速度（约为9.8 m/s²）。如果引入空气阻力，情况变得更加复杂，因为阻力与速度成正比或平方关系。此时需通过数值方法求解微分方程。 #### 加入空气阻力的功能 空气阻力通常可以用以下形式表示： \[ F_d = -k v^n \] 其中： - \( k \) 表示比例系数， - \( n \) 取值一般为1（线性模型）或2（二次模型）， - \( v \) 是瞬时速度。 当 \( n=1 \) 时，称为斯托克斯阻力；\( n=2 \) 则更接近真实环境中的高速流动条件。 以下是基于Python的一个简单实现方案，采用欧拉法求数值解： ```python import numpy as np import matplotlib.pyplot as plt def free_fall_with_air_resistance(h0, mass, drag_coefficient, time_step=0.01, gravity=9.8): """ Simulate a falling object with air resistance. Parameters: h0 : float Initial height of the object (meters). mass : float Mass of the object (kg). drag_coefficient : float Drag coefficient for air resistance. time_step : float Time step size used in simulation (seconds). gravity : float Gravitational acceleration (m/s^2). Returns: times : list[float] List of simulated times. heights : list[float] Corresponding heights at each time point. velocities : list[float] Velocities corresponding to each time point. """ # Initialize variables velocity = 0.0 # initial velocity is zero height = h0 # start from given height times = [0.0] heights = [height] velocities = [velocity] while height >= 0: # continue until ground level reached force_gravity = -mass * gravity # gravitational force downwards (-ve direction) force_drag = -(drag_coefficient / mass) * abs(velocity)**1 * velocity # linear model net_force = force_gravity + force_drag # total forces acting on body acceleration = net_force # since F = ma -> a = F/m and here we assume unit mass # Update state using Euler method velocity += acceleration * time_step # update velocity based on current acceleration height += velocity * time_step # update position according to new velocity # Record data points times.append(times[-1] + time_step) heights.append(max(height, 0)) velocities.append(velocity) return times[:-1], heights[:-1], velocities[:-1] # Example usage if __name__ == "__main__": H0 = 100 # meters MASS = 0.5 # kg DRAG_COEFFICIENT = 0.1 # arbitrary value; adjust accordingly ts, hs, vs = free_fall_with_air_resistance(H0, MASS, DRAG_COEFFICIENT) # Plot results fig, ax1 = plt.subplots() color = 'tab:red' ax1.set_xlabel('Time (s)') ax1.set_ylabel('Height (m)', color=color) ax1.plot(ts, hs, label="Height", color=color) ax1.tick_params(axis='y', labelcolor=color) ax2 = ax1.twinx() color = 'tab:blue' ax2.set_ylabel('Velocity (m/s)', color=color) ax2.plot(ts, vs, label="Velocity", color=color) ax2.tick_params(axis='y', labelcolor=color) fig.tight_layout() plt.title("Free Fall Simulation With Air Resistance") plt.show() ``` 此脚本定义了一个函数 `free_fall_with_air_resistance` 来模拟带空气阻力的自由落体运动，并绘制随时间变化的高度和速度曲线图。注意这里选择了简单的线性拖曳项作为例子展示[^3]。