在理想状况下(即忽略空气阻力)，对站台上的人：

$F=\Delta PS-F_f$

人的受力面积的数值可以近似为人的体积的数值：

$S=V$

根据伯努利效应，因为空气流速的增大，人前后产生的压强差：

$\Delta P=P_0-P$

人对地面的摩擦力：

$F_f=\mu mg$

其中 $\Delta P$为人所受空气压强差，S为人的受力面积， $F_f$为空气摩擦力。

总方程确定之后，接下来具体计算各个物理量。根据伯努利方程，计算因空气流速增大产生的压强：

$P=C-\frac{1}{2}\rho V^2-\rho gh$

$v=\frac{Q}{S}=\frac{Q}{dh}$

其中Q为空气流量，h为人的高度。

综合上述公式，可以得到：

$F=\left(P_0-C+\frac{1}{2}\rho \frac{Q^2}{d^2{h}^2}+\rho gh\right)V-\mu mg$

其中自变量为d、m、V，因变量为F。

运用控制变量法，对d、m、V三个变量分别进行定量分析。当m一定时，建立d、V与F的三维图像；当d一定时，建立m、V与F的三维图像；当V一定时，建立m、d与F的三维图像 [5] [6] ：

由 图1 分析可得，当人的体重一定时，人受到的力与人的体积成正比例关系，与人的距离成反比例关系。

由 图2 分析可得，当人的距离一定时，人受到的力与人的体积成正比例关系，与人的体重成反比例关系。由 图3 分析可得，当人的体积一定时，人受到的力与人的距离成反比例关系，与人的体重成反比例关系。

在实际情况下(即考虑空气阻力)，对站台上的人：

$F=\Delta PS-F_f-F_{空}$

$F_{空}=\frac{\text{1}}{\text{2}}C\rho S{V}_{相}^{\text{2}}$

式中C为空气阻力系数，该值通常是实验值，和物体的特征面积(迎风面积)，物体光滑程度和整体形状有关； $\rho$为空气密度，正常的干燥空气可取1.293 g/l，特殊条件下可以实地监测；S为物体迎风面积； $v_{相}$为物体与空气的相对运动速度。从上式可以得出，考虑空气阻力后，对比在理想情况下人受到的推力，只有人的体积这一自变量发生了变化，所以我们针对人的体积这一变量，利用MATLAB做出考虑空气阻力后的图像。

当d一定时，建立m、V与F的三维图像；当m一定时，建立d、V与F的三维图像。

若考虑空气阻力，利用MATLAB分析，与理想状况下结论相似。我们通过比较 图4 与 图2 ，以及比较 图5 与 图1 ，可知由于空气阻力的影响，计算得出的F较理想状态下数值偏小。这里就不赘述。