实验五_常微分方程及追击问题.pptVIP

实验五 常微分方程及追击问题 东南大学 数学系 * 实验五 常微分方程及追击问题 * 1、求解常微分方程的精确解 在微分方程的学习中，对于一些特殊类型的微分方程可以求得其精确解，在Mathematica中对于这些类型，利用内部函数同样可以求得。 例1 求解微分方程，并作出其积分曲线。 解：打开Mathematica并输入语句： 下面在同一坐标系中做出这个微分方程的6条积分曲线，如图1（其中分别对应）。 例2 求微分方程满足初始条件的特解。 解：只要输入语句： 运行即可得到结果 微分方程的数值解 例3 求在区间上初值问题的数值解，并作出数值解的图形。 解：输入命令： 此时返回的解放在一个表中，函数y[x]不能调用进行运算，要绘出数值解的图形，先要用Evaluate”命令把它转化为可运算的。为此先把刚才得到的结果取代为y[x]后定义为f[x]，然后在区间[0，10]上作出f[x]即近似解的图形。输入命令如下：运行后便得到近似解的图形。 3、追击问题 下面看一个实际问题：我缉私艇雷达发现，正东1 海里处一艘走私船正以常速向北方向逃窜，缉私艇立即以的速度追赶，借助于雷达，缉私艇航行的方向始终对准走私船。试求缉私艇的航行曲线方程，并问走私船航行多远时被我缉私艇追上。 （一）建立微分方程模型： 首先如图建立坐标系。 此为可降阶的二阶微分方程，可解得缉私艇的航行曲线方程为： 当时，，故走私船航行海里时被缉私艇追上。 两边对x求导，则实际问题化为求解微分方程： 运行后得到： 显然第一个解为所求的方程，为了求得当时函数得值，可由以下命令得到： 在这个问题上，就是一步步地以时间间隔为来模拟缉私艇追踪敌艇的实际过程。 （二）仿真方法：即模仿真实事件的行为和过程。 如图，当时，敌艇在处，缉私艇在处，方向指向敌艇，即。认为缉私挺的运动曲线是由水平方向和垂直方向组合而成。 当时，敌艇运行到，缉私艇运行到点，则 运行结果的最后一项为： 从此结果可以看出，当时迭代结束，此时追赶的时间为小时，说明敌艇航行了0.667497海里被缉私艇追上，这和前面的结果基本上一致。 运行后可得： 即方程的通解为。 对于一阶微分方程，除了可分离变量的微分方程，线性微分方程也可用Mathematica求出精确解，请同学们自己举例求解。除了这两类方程外，由实际操作可知，用Mathematica精确求解都比较困难，为此必须要考虑求微分方程的数值解问题。 运行后输出的结果为： 它表示得到了插值函数InterpolatingFunction[domain, table]类型的近似解，近似解的定义域domain为[0,10] 设缉私艇的航行曲线方程为， 在时刻t位于，则有 又，消去t，得。 或者用mathematica软件来求解上述微分方程，输入命令： 当时，敌艇运行到， 而缉私艇运行到点，则 方向由点指向点，即； 方向由点指向点，即； 于是，可以得到仿真的迭代格式为： 利用mathematica可以对该实际问题进行模拟，不妨设海里/小时，则输入程序： 在程序中，设定时间间隔为小时，并且在且迭代结束。 另外，可以在程序中将改变进行模拟，而从运行结果来看不会改变敌艇被追上时的航行路程，改变的只是被追击上的时间。 运行后也可得，即走私船航行海里时被缉私艇追上。