《1、常微分方程的常见解法》由会员分享,可在线阅读,更多相关《1、常微分方程的常见解法(98页珍藏版)》请在金锄头文库上搜索。
1、一、 向量场设一阶微分方程 满足解的存在唯一性定理的条件。平面的一个区域的右端函数在中有定义,那么,过中任一点 有且仅有的一个解 ,满足从几何方面看,解 就是通过点 的一条常微分方程的解法介绍曲线(称为积分曲线),且 就是该曲线上的点 处的切线斜率,特别在 切线斜率解,但我们知道它的解曲线在区域D中任意点 的切线斜率是 。就是 尽管我们不一定能求出方程 的如果我们在区域D内每一点 处,都画上一个就得到一个方向场,将这个方向场称为由微分方程所确定的向量场。的值为斜率中心在 以点的线段,我们它所确定的向量场中的一条曲线,该曲线所经过的从几何上看,方程 的一个解 就是位于每一点都与向量场在这一点的方
2、向相切。行进的曲线,因此,求方程满足初始值的这样的一条曲线。的解,就是求通过点形象的说,解就是始终沿着向量场中的方向向量场对于求解微分方程的近似解和研究微分方程的几何性质极为重要,因为,可根据向量场的走向来近似求积分曲线,同时也可根据向量场本身的性质来研究解的性质。例1.3.1 在区域 内画出方程 的向量场和几条积分曲线。解:用计算各点的斜率的方法手工在网格点上画出向量场的方向可以得到向量场,但手工绘图误差较大。我们可以用Maple 软件包来完成。点的向量相重合。 L在每点均与向量场的向量相切。 在L上任一点,L的切线与所确定的向量场在该定理1.3L为的积分曲线的充要条件是:曲线Maple指令
3、:DEtoolsphaseportrait # 画向量场及积分曲线 (diff(y(x),x)=-y(x),y(x), # 定义微分方程 x=-2.2, # 指定x范围 y(-2)=2,y(-2)=1,y(-2)=-2, # 给出3个初始值 dirgrid=17,17, # 定义网格密度 arrows=LINE, # 定义线段类型 axes=NORMAL); # 定义坐标系类型在MATLAB的向量场命令为 quiver(x,y,px,py) 回车后Maple就在 条积分曲线,见下图 的图形,并给出了过点的网格点上画出了向量场的三 所谓图解法就是不用微分方程解的具体表达式, 直接根据右端函数的结
4、构和向量场作出积分曲线 的大致图形。 图解法只是定性的,只反映积分曲线的一部分主要特征。 该方法的思想却十分重要。因为能够用初等方法求解的方程极少,用图解法来分析积分曲线的性 态对了解该方程所反映的实际现象的变化规律就 有很重要的指导意义。二、 积分曲线的图解法三、一阶常微分方程的解法1线性方程2 变量可分离方程3 全微分方程4 变量替换法5 一阶隐式方程6 近似解法7 一阶微分方程的应用初值问题 的解为 初值问题 的解为 Bernoulli方程求出此方程通解后,令解法:伯努利方程的标准形式:除方程两边 , 得换回原变量即得伯努利方程的通解。例 湖泊的污染设一个化工厂每立方米的废水中含有3.0
5、8kg盐酸,这些废水流入一个湖泊中,废水流入的速率20 立方米每小时. 开始湖中有水400000立方米. 河水 中流入不含盐酸的水是1000立方米每小时, 湖泊中混合均匀的水的流出的速率是1000立方米每小时,求该厂排污1年时, 湖泊水中盐酸的含量。解: 设t时刻湖泊中所含盐酸的数量为考虑 内湖泊中盐酸的变化。因此有该方程有积分因子两边同乘以后,整理得积分得利用初始条件得当 , 得 变量可分离方程的求解方程(2.2.1)两边同除以 这样对上式两边积分得到齐次函数: 函数称为m次齐次函数, 如果齐次方程: 形如的方程称为齐次方程。 引入一个新变量化为变量可分离方程求解思想:求解。齐次方程可化为齐
6、次方程的方程形如的方程可化为齐次方程.其中都是常数.1. 当时, 此方程就是齐次方程 . 2. 当时, 并且(1)此时二元方程组有惟一解引入新变量此时, 方程可化为齐次方程:(2) 若则存在实数使得:或者有不妨是前者, 则方程可变为令则4. 对特殊方程令则例 求方程 的通解。 解:解方程组 得 令 代入原方程可得到齐次方程令 得还原后得原方程通解为变量分离后积分例:雪球融化问题设雪球在融化时体积的变化率与表面积成比例,且融化过程中它始终为球体,该雪球在开始时的半径为6cm ,经过2小时后,其半径缩小为3cm。求雪球的体积随时间变化的关系。 解:设t时刻雪球的体积为 ,表面积为 ,球体与表面积的
7、关系为 变量可分离方程的应用由题得引入新常数 再利用题中的条件得分离变量积分得方程得通解为再利用条件 确定出常数C和r代入关系式得 t的取值在 之间。 中连续且有连续的一阶偏导数,则 定理2.1 设函数 和 在一个矩形区域是全微分方程的充要条件为:(2.3.3)方程为全微分方程的充要条件例:验证方程是全微分方程,并求它的通解。由于 3.全微分方程的积分解:当一个方程是全微分方程时,我们有三种解法 .(1) 线积分法:或由公式(2.3.4)得 : 故通解为其中为任意常数所以方程为全微分方程。(2)偏积分法 的通解.例:求方程由于 解: 假设所求全微分函数为 ,则有 求 而 即从而即例:验证方程是
8、全微分方程,并求它满足初始条件: 的解。 解: 所以方程为全微分方程。 由于 由于 (3)凑微分法方程的通解为: 利用条件 得 最后得所求初值问题得解为:根据二元函数微分的经验,原方程可写为四、微分方程的近似解法 用一些函数去近似微分方程的解 在一些点上计算方程解的近似值 逐次迭代法 Taylor级数法 Euler折线法 Runge-Kutta法能得到解析解的方程:线性方程、变量可分离的方程、 全微 分方程以及能通过各种方法化为这些类 型的方程. 绝大部分方程无法求得解析解,一些近似 解法也对实际问题的解决有很大帮助,我 们需要讨论在得不到解析解时寻求近似解 的方法。对初始值问题构造迭代序列
9、该序列一致收敛到解,故迭代一定次 数后就可以作为一个近似1、 逐次迭代法解:该初值问题近似解的迭代序列 如下例 求初值问题的近似解迭代的误差 (|x|1+(y-x)2; f2:=(x,y)-2*(x-y)+2*(y-x)*(1+(y-x)2); for n from 0 to 9 do xn+1:=h*(n+1); yn+1:=yn+h*f1(xn,yn); zn+1:=zn+h*f1(xn,zn)+h2*f2(xn,zn)/2; un+1:=xn+1+1/(2-xn+1); print (xn+1,yn+1,zn+1,un+1); od: 可以改变步长和增加分点来观察计算精度的变化情况对于常
10、微分方程的边值问题的解即- (1)Runge-Kutta(龙格 - 库塔)法Runge-Kutta方法的导出有上使用微分中值定理,在区间-(2)引入记号的近似值K。就可得到相应的-(3)Runge-Kutta方法即(3)式只要使用适当的方法求出y(x)上平均斜率在区间K可以认为是在区间上的平均斜率。低阶Runge-Kutta方法如下图即则(4)式化为即Euler方法Euler方法也称为一阶Runge-Kutta方法由于-(4)(由(4)式)令则(3)式化为-(5)称为二阶Runge-Kutta法高阶Runge-Kutta方法未知令令)(2111-nnxyKx预测处的斜率如果以取则-(6)(6)
11、式称为三阶Runge-Kutta方法还可构造四阶(经典)Runge-Kutta方法四阶(经典)Runge=Kutta方法有4阶精度例 求初始值问题的数值解 利用四阶Runge=Kutta方法计算机编程给出步长和初始值循环计算各点上函数的近似值显示结果printlev1:=0: h:=0.1: x0:=0: y0:=0.5: f:=(x,y)-1+(y-x)2; for n from 1 to 10 do xn:=h*n; k1:=f(xn-1,yn-1); k2:=f(xn-1+h/2,yn-1+k1*h/2); k3:=f(xn-1+h/2,yn-1+k2*h/2); k4:=f(xn-1+
12、h,yn-1+k3*h); yn:=yn-1+h*(k1+2*k2+2*k3+k4)/6; un:=xn+1/(2-xn); print (xn,yn, un); od:运行结果 适应范围与变化率有关的各种实际问题 应用三步曲(1) 建模即根据实际问题建立起适当的微分 方程, 给出其定解条件. (2) 求解求出所建立的微分方程的解(3) 翻译用所得结果来解释一些现象,或对 问题的 解决提出建议或方法建议: 模型要详略得当在用微分方程解决实际问题的过程中一 定要意识到实际问题是十分复杂的,微分 方程只能是在一定程度上对问题的一种近 似描述,只要结果的误差在一定范围内即 可.任何模型都不可能把影响
13、问题的所有 因素都反映在微分方程中,或者要求所得 结果 十分精确.一个好的微分方程模型是在实际 问题的精确性和数学处理的可能性之间的 一个平衡.有一段时间,美国原子能委员会(现为核管理委员会)是 这样处理浓缩放射性废物的,他们把这些废物装入密封性能 很好的圆桶中,然后扔到水深300英尺的海里。这种做法是否会造成放射性污染,很自然地引起了生态 学家及社会各界的关注。原子能委员会一再保证,圆桶非常 坚固,决不会破漏,这种做法是绝对安全的。然而一些工程 师们却对此表示怀疑,他们认为圆桶在和海底相撞时有可能 发生破裂。而原子能委员会有专家们则仍然坚持自己的看 法。于是,双方展开了一场笔墨官司。究竟谁的
14、意见正确呢?看来只能让事实说话了。问题的 关键在于圆桶到底能承受多大速度的碰撞,圆桶和海底碰撞 时的速度有多大?放射性废物的处理大量破坏性实验,发现圆桶在40英尺秒的冲撞下 会发生破裂,剩下的问题就是计算圆桶沉入300英尺深 的海底时,其末速度究竟有多大了。美国原子能委员会使用的是55加仑的圆桶,装满放 射性废物时的圆桶重量为W527.436磅,而在海水中 受到的浮力B470.327磅。此外,下沉时圆桶还要受 到海水的阻力,阻力Dv,其中C为常数。工程师们 做了大量实验,测得C0.08。现在,取一个垂直向下 的坐标,并以海平面为坐标原点(0)。于是,根据 牛顿第二定律建立圆桶下沉时应满足方程质
15、量加速度=重力-浮力-摩擦阻力 模型及其解oymgBD困难:无法知道下沉到海底的时间积分和代入初始条件得:最后再用数值计算可以得到水深300时的速度 大小。借助数值方法求出v(300)的近似值。计算结果 表明,v(300)45.1英尺秒40英尺秒。工程师们的猜测是正确的,他们打赢了这场官 司。现在,美国原子能委员会已改变了他们处理放 射性废物的方法,并明确规定禁止将放射性废物抛 入海中。 一横截面积为常数A,高为H的水池内 盛满了水,由池底一横截面积为B的小 孔放水. 求在任意时刻的水面高度和 将水放空所需的时间 .例:例:水的流出时间水的流出时间: 有高为1米的半球形容器, 水从它的底部 小孔流出, 小孔横截面积为1平方厘米(如图). 开始时容器内盛满了水, 求水从小孔流出过程中 容器里水面的高度h(水面与孔口中心间的距离) 随时间t的变化规律.例 解: 由力学知识得,水从孔口 流出的流量为流量系数孔口截面面积重力加速度值得进一步探讨的问题: 不同的形状设在微小的时间间隔水面的高度由h 降至