数值分析常微分方程的差分方法精.ppt

《数值分析常微分方程的差分方法精.ppt》由会员分享，可在线阅读，更多相关《数值分析常微分方程的差分方法精.ppt（41页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数值分析常微分方程的差分方法第1页，本讲稿共41页问题的提出问题的提出 实际中，很多问题的数学模型都是微分方程。我们可实际中，很多问题的数学模型都是微分方程。我们可以研究它们的一些性质。但是，只有极少数特殊的方程有解以研究它们的一些性质。但是，只有极少数特殊的方程有解析解。对于绝大部分的微分方程是没有解析解的。析解。对于绝大部分的微分方程是没有解析解的。常微分方程作为微分方程的基本类型之一，在自然界与常微分方程作为微分方程的基本类型之一，在自然界与工程界有很广泛的应用。很多问题的数学表述都可以归结为工程界有很广泛的应用。很多问题的数学表述都可以归结为常微分方程的定解问题。很多偏微分方程问题，也

2、可以化为常微分方程的定解问题。很多偏微分方程问题，也可以化为常微分方程问题来近似求解。常微分方程问题来近似求解。第2页，本讲稿共41页常微分方程的定解问题常微分方程的定解问题 考虑考虑一阶一阶常微分方程的常微分方程的初值问题初值问题只要只要 f(x,y)在在a,b R1 上连续，且关于上连续，且关于 y 满足满足 Lipschitz 条件条件，即存在与即存在与 x,y 无关的常数无关的常数 L 使使对任意定义在对任意定义在 a,b 上的上的 y1(x)和和 y2(x)都成立，则上述问题都成立，则上述问题存存在唯一解在唯一解。第3页，本讲稿共41页差分方法差分方法要计算出解函数要计算出解函数 y

3、(x)在一系列节点在一系列节点 a=x0 x1 n）上产生的偏差均不超过）上产生的偏差均不超过，则称该方法是稳定的则称该方法是稳定的.稳定性问题比较复杂，为简化讨论，我们仅考察下列模型稳定性问题比较复杂，为简化讨论，我们仅考察下列模型方程方程 y=y，0 第31页，本讲稿共41页收敛性与稳定性收敛性与稳定性模型的欧拉格式为模型的欧拉格式为yn+1=（1+h）yn 模型的欧拉格式为模型的欧拉格式为则则n+1=（1+h）n要使要使|yn+1|yn|则则|1+h|1稳定条件稳定条件0 h-2/第32页，本讲稿共41页收敛性与稳定性收敛性与稳定性模型的隐式欧拉格式为模型的隐式欧拉格式为yn+1=yn+

4、hyn+1 解出解出恒成立恒成立总有总有结论结论恒稳定恒稳定|yn+1|yn|第33页，本讲稿共41页方程组与高阶方程的情形方程组与高阶方程的情形一阶方程组的一般形式一阶方程组的一般形式第34页，本讲稿共41页方程组与高阶方程的情形方程组与高阶方程的情形化高阶方程为一阶方程化高阶方程为一阶方程第35页，本讲稿共41页方程组与高阶方程的情形方程组与高阶方程的情形令令则有则有第36页，本讲稿共41页边值问题边值问题考虑常微分方程的边值问题：考虑常微分方程的边值问题：其中其中p(x)，q(x)和和f(x)均为均为a,b上给定的函数上给定的函数，为已知数为已知数。假定假定p(x)、q(x)及及f(x)

5、均为均为a,b上充分光滑的函数上充分光滑的函数，且且q(x)0，这时，边值问题存在连续可微的解，且唯一。这时，边值问题存在连续可微的解，且唯一。第37页，本讲稿共41页边值问题边值问题用差分法解边值问题的主要步骤是：用差分法解边值问题的主要步骤是：（1）将区间将区间a,b离散化离散化；（2 2）在这些节点上，将导数差商化，从而把微分方程）在这些节点上，将导数差商化，从而把微分方程 化为差分方程；化为差分方程；(3)解差分方程解差分方程实际上就是解线代数方程组。实际上就是解线代数方程组。将将a,b区间用节点区间用节点分成分成N等分，其中等分，其中x0=a与与xN=b 称为边界点，称为边界点，而而x1,x2,xN-1称为内点。称为内点。第38页，本讲稿共41页边值问题边值问题例例9.7 试用差分法解方程试用差分法解方程解解 将将0,1划分为四等分，即取划分为四等分，即取 ，得五个节点得五个节点差分方程为差分方程为第39页，本讲稿共41页边值问题边值问题将它改写成将它改写成 在每个内点列方程得在每个内点列方程得 由追赶法公式解得：由追赶法公式解得：y3=1.4855 y2=1.2802 y1=0.7753第40页，本讲稿共41页本章结束本章结束谢谢大家谢谢大家！第41页，本讲稿共41页