为网球发球飞行建模：对已发球的网球的路径进行建模，并创建模型结果的 3D 可旋转可视化。 包括升力和阻力。-matlab开发

用户定义网球的初始位置、速度和旋转。 单位是磅、英尺、秒。 网球场坐标为： x - 地面，垂直于网y - 地面，平行于网络z - 离地高度所采取的路径取决于 3 个力： 1) 重力（在 -z 方向恒定 32 ft/s^2 2) 拖动（与 v^2 成比例并在 -v-hat 方向上） 3）升力（与自旋和速度的叉积成正比） 阻力的比例常数与大气密度、球的横截面积、阻力系数和质量的倒数成正比： k_D = C_D*rho_atmosphere*area/(2*mass) 球的旋转使其升起。 常数有类似的形式： k_L = C_L*rho_atmosphere*area/(2*mass) 投放时旋转会增加错误的余地。 上旋（+y 分量）会产生向下的力，导致球下落并允许发球者以更大的初始角度击球。 后旋是通过将 y 分量设置为负值来创建的。 传播是通过 4/5 阶的 Runge-Kutta-Fehl