飞行器轨迹仿真大作业.docxVIP

飞行器轨迹仿真大作业[文档副标题]1. 积分方法欧拉积分二阶龙格库塔四阶龙格库塔2. 计算例子实现自由落体运动的轨迹求解。，阻力系数，小球半径, ，重力加速度取常数，大气密度采用插值计算。求落地速度和落地时间，观察速度的变化规律。大气密度变化规律：[高度/M] 0.0 1000 3000 5000 7000 9000 11000 15000 20000[密度/kg/m^3] 1.225 1.112 0.909 0.736 0.589 0.466 0.364 0.194 0.0883. 图形3.1 用Euler法画出的图图1 速度随时间的变化图图2 高度随时间的变化3.2 用Runge-Kutta法画的图形图3 高度随时间变化图图4 速度随时间变化图4. 结果分析表一 不同积分法对应的落地时间和落地速度Euler落地时间RungeKutta落地时间Euler落地速度RungeKutta落地速度59.8659.85-345.984-345.949由迭代结果可知，小球的高度，速度变化规律是：开始加速度为负，并逐渐趋近于0，速度由100m/s逐渐减小到0，此时高度逐渐上升但上升的速度逐渐变慢，到达最高点后速度由正变负，绝对值不断变大，小球加速下降。图5 有无空气阻力时速度时间图像对比图6 有无空气阻力时