简单微分方程求解及导弹追踪问题.ppt

实验 8 简单微分方程求解及导弹追踪问题 内容提要 在实际问题中，我们要要研究的往往是变动的量以及它们之间的关系，由于这些量是时刻变化的，因此它们之间的关系不能用简单的代数方程，而要用微分方程来表达，从而寻求微分方程的解成为一个十分重要的任务，这里我们首先给出利用Mathematia解微分方程的方法，以及求微分方程数值解的方法，最后给出“导弹追踪及其模拟”的演示。 简单微分方程求解及导弹追踪问题 实验步骤 微分方程精确解 在Mathematica中对于简单的微分方程可以求得其精确，但能精确求解的微分方程类型是非常有限的，下面按照我们学过的微分方程的类型来考察。 可分离变量的微分方程 简单微分方程求解及导弹追踪问题 例 1求解微分方程=ey，并作出其积分曲线。 解 打开Mathematica后键入： 并运行，即得到微分方程得解为 下面我们在同一坐标系中作出这个微分方程的积分曲线（对应C=-3,-2,-1,0,1,2,3）,我们键入 运行即得积分曲线图37（我们将在下册利用多元函数等值线同样可以作出微分方程的积分曲线）. 简单微分方程求解及导弹追踪问题 线性微分方程 例 2 求解微分方程 ，并作出其积分曲线。 解 打开Mathematica后键入 并运行，即得到微分方程的解为 y= 简单微分方程求解及导弹追踪问题 为在同一坐标系中作出这个微分方程的积分曲线（对应），只须键入： 运行即可得到4条积分曲线（图38）。 简单微分方程求解及导弹追踪问题 通过实际操作可知，除了上述两类方程外，用精确求解都比较困难，为此必须考虑数值解问题。 微分方程数值解 在区间[1，4]上求出微分方程 ,的数值解，并作出数值解的图形。 解 打开Mathematica后键入： 简单微分方程求解及导弹追踪问题 运行后即得数值解为： 并已储存在变量sol中，为作出数值解的图形，我们首先把刚才得到的结果取代为y[x]后定义为函数f[x]，最后在[1,4]区间上作出f[x]的图形，即键入： 简单微分方程求解及导弹追踪问题 运行即可（图39）。这里要注意的是：1）要使得到的数值能运算，必须使用命令”Evalute”;2）得到的数值解只在求解范围内才意义（这里即为区间[1，4]）。 * *