Ch3修改稿.ppt

例 23 求解微分方程 再输入命令 为方便求解, 先交换变量, 则有 结果为 整理后再交换变量, 得方程的通解为 即： 例24 求一阶微分方程 时的特解,并画出特解的图形. 程序为 的通解及 subs用法 举例说明： 相应的积分曲线为 例 25 求二阶微分方程 这里 程序为 练习 求解微分方程 命令为 上述方法可应用于求解微分方程组. 格式 且满足 例 26 求微分方程组 的特解, 并作出相应的图形. 应用 盐水的混合问题 一个圆柱形的容器, 内装 的均匀混合的盐水溶液, 如果纯水以每秒 的速度从容器顶部流入, 同时容器内 的混合的盐水又以每秒 的速度流出, 开始时, 容器 内盐的含量为 问经过时间 秒后, 容器内盐的含量 为多少? 该类问题具有很明显的微分方程模型的特征. 14升/秒，纯水 10.5升/秒，盐水 开始时：350升盐水 7公斤盐 假设: ⑴容器足够大, 在一定时间内不会溢出; ⑵注入后的纯水立刻会成为均匀的溶液. 模型建立 设在时刻 时, 容器内溶液的含盐量为 此时容 器内的溶液为 升，浓度为 在某一极短的时间 内，盐的改变量为： 因而盐量的改变量为 由此得到微分方程: 另一方面，由于 极短，所以在 时间内， 可以认为盐水浓度不发生变化， 流出的盐量 求解程序为 方程 解析解为： 图形为 ⑵数值解法 格式 区间 初始值 被积函数 例 27 求解微分方程: 这里 输入语句: 区间 初始值 画出积分曲线 结果为 例 28 用数值方法求解初值问题 并画出相 应的图形. 程序为 二元函数 将积分曲线与函数 的图形加以对比: 应用 Volterra问题求解 自然界中不同种群之间存在着一种非常有趣的既有依 存, 又有制约的生存方式: 种群甲靠丰富的自然资源生 长, 而种群乙靠捕食种群甲为生, 食用鱼和鲨鱼、美洲 兔和山猫、落叶松和蚜虫等都是这种生存方式的典型. 生态学上称种群甲为食饵（ Prey）, 种群乙为捕食者 （ Predator）, 二者共处组成食饵-捕食者系统（简称P-P 系统）. 近百年来许多数学家和生态学家对这一系统进 行了深入的研究, 建立了一系列的数学模型, 下面介绍 P-P系统最初的、最简单的一个模型—Volterra模型. 1.模型建立 捕食者（乙）数量为 设在时刻 时食饵（甲）数量为 甲独立生存的增长率为 即 乙的存在使甲的增长率减小, 减少量与 成正比, 即 ⑴ 乙独立生存的死亡率为 即 甲的存在使乙的死亡率减少, 减少量与 成正比, 即 ⑵ 上两式中的系数 分别表示捕食者掠取食饵的能力 食饵供养捕食者的能力. 方程⑴与⑵反应的是自然环境下食饵与捕食者之间相 互依存和相互制约的关系, 注意到, 该模型中并没有考虑 种群自身的阻滞增长关系. 模型分析 方程⑴和⑵无解析解! 2.模型求解 用MatLab中微分方程的数值解方法可以得到该问题的 数值解, 并用相轨线可以对模型作进一步的分析. 首先建立函数文件 再建立脚本文件 运行该文件: 鼠标控制输入字符 解：设T(x,y)表示金属板上点(x,y)的温度 则： 且设蚂蚁每向前走一步的长度为 米 现在假设蚂蚁从起点(11,51)处迈出第一步，往哪个 方向走呢？ 这个方向要能让温度下 降的最快，毫无疑问，就是 是在起点处沿梯度的方向 或梯度的反方向走. 下面求温度函数 在点 处 的梯度： 在matlab命令窗口中输入： 从而得到 在matlab命令窗口中继续输入： 可以看到： 0 负号表示什么呢？ 表示的方向如下图所示 因此，解题时要注意这一点，温度下降的方向应该是梯度的反方向 下面将矩形及第一步的路线画出来 顺时针画矩形 画第一步路线：必须确定第一步路线的两个端点， 第一个端点是起点，第二个端点的坐标确定如下： 起点 11 51 继续在matlab命令窗口中输入： x1 y1 局部放大后： 从第一步的终点出发，重复上述过程，可以 得到蚂蚁的爬行路线，程序编制如下： 蚂蚁爬行路线如下： 显然，蚂蚁沿着红线继续爬行，到达边界后，还会沿边界爬行，这不是本程序所能解决的。 4.多元函数的极值 求极值方法 格式1 功能 用单纯性方法求函数在指定处附近的极值. 格式2