《常微分方程的常见解法》.pptVIP

定理 矩阵 （4.4.21） 是方程组（4.4.1）的基解矩阵，且 证明 当t=0时，由定义知 ，又因为 这表明， 是（4.4.1）的解矩阵，又因为 （4.4.22） 所以 是方程组（4.4.1）的基解矩阵. 由定理6可知，方程组（4.4.1）的通解为 这里c是一个常数向量. 若 是方程组（4.4.1）满足初始条件 的解，则由（4.4.22）知 即有 （4.4.23） * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 2.方程为常微分方程的充要条件 * * 3.全微分方程的积分 还可构造四阶(经典)Runge-Kutta方法 四阶(经典)Runge=Kutta方法有4阶精度 例 求初始值问题的数值解 利用四阶Runge=Kutta方法计算机编程 给出步长和初始值 循环计算各点上函数的近似值 显示结果 printlev1:=0: h:=0.1: x[0]:=0: y[0]:=0.5: f:=(x,y)-1+(y-x)^2; for n from 1 to 10 do x[n]:=h*n; k1:=f(x[n-1],y[n-1]); k2:=f(x[n-1]+h/2,y[n-1]+k1*h/2); k3:=f(x[n-1]+h/2,y[n-1]+k2*h/2); k4:=f(x[n-1]+h,y[n-1]+k3*h); y[n]:=y[n-1]+h*(k1+2*k2+2*k3+k4)/6; u[n]:=x[n]+1/(2-x[n]); print (x[n],y[n], u[n]); od: 运行结果 适应范围 与变化率有关的各种实际问题 应用三步曲 (1) 建模 即根据实际问题建立起适当的微分方程， 给出其定解条件. (2) 求解 求出所建立的微分方程的解 (3) 翻译 用所得结果来解释一些现象，或对问题的 解决提出建议或方法 建议: 模型要详略得当 在用微分方程解决实际问题的过程中一定要意识到实际问题是十分复杂的，微分方程只能是在一定程度上对问题的一种近似描述，只要结果的误差在一定范围内即可.任何模型都不可能把影响问题的所有因素都反映在微分方程中，或者要求所得结果 十分精确.一个好的微分方程模型是在实际问题的精确性和数学处理的可能性之间的一个平衡. 有一段时间,美国原子能委员会(现为核管理委员会)是这样处理浓缩放射性废物的,他们把这些废物装入密封性能很好的圆桶中,然后扔到水深300英尺的海里。 这种做法是否会造成放射性污染,很自然地引起了生态学家及社会各界的关注。原子能委员会一再保证,圆桶非常坚固,决不会破漏,这种做法是绝对安全的。然而一些工程师们却对此表示怀疑,他们认为圆桶在和海底相撞时有可能发生破裂。而原子能委员会有专家们则仍然坚持自己的看法。于是,双方展开了一场笔墨官司。 究竟谁的意见正确呢?看来只能让事实说话了。问题的关键在于圆桶到底能承受多大速度的碰撞,圆桶和海底碰撞时的速度有多大? 放射性废物的处理 大量破坏性实验,发现圆桶在40英尺／秒的冲撞下会发生破裂,剩下的问题就是计算圆桶沉入300英尺深的海底时,其末速度究竟有多大了。 美国原子能委员会使用的是55加仑的圆桶,装满放射性废物时的圆桶重量为W＝527.436磅,而在海水中受到的浮力B＝470.327磅。此外,下沉时圆桶还要受到海水的阻力,阻力D＝Ｃv,其中C为常数。工程师们做了大量实验,测得C＝0.08。现在,取一个垂直向下的坐标,并以海平面为坐标原点(ｙ＝0)。于是,根据牛顿第二定律建立圆桶下沉时应满足方程 质量·加速度=重力-浮力-摩擦阻力 模型及其解 o y mg B D 困难：无法知道下沉到海底的时间 积分和代入初始条件得： 最后再用数值计算可以得到水深300时的速度大小。 借助数值方法求出v(300)的近似值。计算结果表明, v(300)≈45.1英尺／秒＞40英尺／秒。 工程师们的猜测是正确的,他们打赢了这场官 司。现在,美国原子能委员会已改变了他们处理放 射性废物的方法,并明确规定禁止将放射性废物抛入海中。 一横截面积为常数A,高为H的水池内盛满了水,由池底一横截面积为B的小孔放水. 求在任意时刻的水面高度和将水放空所需的时间 . 例：水的流出时间 ：