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

飞行器轨迹仿真大作业 [文档副标题] 积分方法 欧拉积分 二阶龙格库塔  x ? x ? h ? f (x , t ) n?1 n n n x ? x n?1 h (k ? k ) n 2 1 2 k ? f (t , x ) n n k ? f (t n h, x n hk ) 1 四阶龙格库塔 x ? x ? h (K ? 2K ? 2K ? K ) i?1 K ? 1 i f (t i 6 1 2 3 4 , x ) i K ? f (t ? h , x ? h K ) 2 i 2 i 2 1 K ? f (t ? h , x ? h K ) 3 i 2 i 2 2 计算例子 K ? f (t ? h , x ? hK ) 4 i i 3 实现自由落体运动的轨迹求解。 ? h? ? v ? ? D(v, h) ?v? ? ?g ? ? m v ? 100m / s, h 0 0 ? 10000m, m ? 5kg ，阻力系数 C d ? 0.02 ，小球半径 r ? 0.1m , g ? 9.81m / s2 ，重力加速度取常数，大气密度采用插值计算。求落地速度和落 0 地时间，观察速度的变化规律。 大气密度变化规律： [高度/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.088 图形 用 Euler 法画出的图 100 50 0 -50 v -100 度 速 -150 -200 -250 -300 -350  0 10 20 30 时间t  40 50 60 图 1 速度随时间的变化图 11000 10000 9000 8000 7000 h 度 高 6000 5000 4000 3000 2000  0 10 20 30 时间t  40 50 60 图 2 高度随时间的变化 用 Runge-Kutta 法画的图形 11000 10000 9000 8000 7000 h 度 高 6000 5000 4000 3000 2000  0 10 20 30 时间t  40 50 60 图 3 高度随时间变化图 100 50 0 -50 v -100 度 速 -150 -200 -250 -300 -350  0 10 20 30 时间t  40 50 60 图 4 速度随时间变化图 4. 结果分析 Euler 落地时间 Euler 落地时间 RungeKutta 落地时间 Euler 落地速度 RungeKutta 落地速度 59.86 59.85 -345.984 -345.949 由迭代结果可知，小球的高度，速度变化规律是：开始加速度为负，并逐渐趋近于 0，速度由 100m/s 逐渐减小到 0，此时高度逐渐上升但上升的速度逐渐变 慢，到达最高点后速度由正变负，绝对值不断变大，小球加速下降。 无空气阻力有空气阻力 无空气阻力 有空气阻力 0 -100 v 度 -200 速 -300 -400 -500  0 10 20 30 时间t  40 50 60 图 5 有无空气阻力时速度时间图像对比 无空气阻力有空气阻力 无空气阻力 有空气阻力 10000 9000 8000 7000 h 度 6000 高 5000 4000 3000 2000 1000  0 10 20 30 时间t  40 50 60 图 6 有无空气阻力时高度时间对比 从图中可以看出，随着时间的增加，有阻力时速度和高度和无阻力时的速度和高度差值越来越大。 : 5. 代码 代码 1： Euler 方法 #includeiostream #includefstream #includemath.h #includeiomanip using namespace std; void Euler(double h,double* x,double* f,double D,double p); void Fct(double h,double* x,double* f,double D,double p); double Rho(double* x); int main() { double t0=0; double tf=60; double h=0.01; double x[2]; double f[2]; double D; double p; x[0]=10000; x[1]=100; double t; ofstream outFile(D:\\Matlab\\Kevin\\prime.txt); for(t=