数值计算方法[第6章 常微分方程初值问题的数值解法.pdfVIP

6章 常微分方程初值问题的数值解法 /* Numerical Methods for Ordinary Differential Equations */ 1 问题的提出及基本理论 $考虑一阶常微分方程的初值问题/* Initial-Value Problem */: dy = f ( x , y ) x ˛ [a , b ] dx y (a ) = y 0 定义 微分方程和定解条件共同组成了定解问题。由给定 的初始条件和微分方程所给出的定解问题称为初值问题。 让我们来看看下面两 个应用的例子吧！ 人口接近饱和 例: 人口发展模型 dP (t) P (t ) r(1- )P (t ) K dt K P (t ) P 0 0 t P (t ) t 其中 是所研究地区人口的上界， 是 时刻所研究地区 K 的人口数， 是该地区人口的自然增长率，当 越接近 r P (t ) K 的时候，该地区的人口就越接近饱和，从而人口的自然增长 率会随之而下降。该Logistic模型的解析解为 K K - P ( ) , 0 P t C r t t 1 Ce- ( - 0 ) P 0 例: Lorenz 吸引子(混沌学中的微分方程组) 在所有的常微分方程中，被研究的最多的毫无疑问就是Lorenz方程。该方程 是1963年由气象学家爱德华·罗伦兹(Edward Lorenz)提出的。为了预报天 气，他用计算机求解仿真地球大气的13个方程式。为了更细致地考察结果， 他把一个中间解取出，提高精度再送回。而当他喝了杯咖啡以后回来再看时 竟大吃一惊：本来很小的差异，结果却偏离了十万八千里！计算机没有毛 病，于是，洛伦兹 Lorenz ）认定，他发现了新的现象：“对初始值的极端不 稳定性”，即：“混沌 ”，又称“蝴蝶效应”。 蝴蝶效应 纽约的一只蝴蝶翅膀一拍，风和日丽的杭州就刮起台风来了？！ NY HZ 简化的地球大气环流模型 让我们作个实验吧！ y (t) b 0 y y (t ) 1 2 1 d dt y 2 (t) 0 s s y 2 (t )