MATLAB第十三讲常微分方程初值问题数值解法.pptVIP

MATLAB第十三讲常微分方程初值问题数值解法演示文稿;MATLAB第十三讲常微分方程初值问题数值解法;9.1引言;很多物理系统与时间有关，从卫星运行轨道到单摆运动，从化学反应到物种竞争都是随时间的延续而不断变化的.解常微分方程是描述连续变化的数学语言，微分方程的求解就是确定满足给定方程的可微函数y(t)，研究它的数值方法是本章的主要目的.考虑一阶常微分方程的初值问题;定理1设f在区域D={(x,y)|a?x?b,y?R}上连续,关于y满足利普希茨条件，则对任意x0?[a,b],y0?R，常微分方程初值问题(1.1)式和(1.2)式当x?[a,b]时存在唯一的连续可微解y(x).;这个定理表明解对初值依赖的敏感性，它与右端函数f有关，当f的Lips.常数L比较小时，解对初值和右端函数相对不敏感，可视为好条件.若L较大则可认为坏条件，即病态问题.;虽然求解常微分方程有各种各样的解析方法，但解析方法只能用来求解一些特殊类型的方程，实际问题中归结出来的微分方程主要靠数值解法.;初值问题的数值解法有个基本特点，他们都采取“步进式”，即求解过程顺着节点排列的次序一步一步地向前推进.描述这类算法，只要给出用已知信息yn,yn-1,yn-2,?计算yn+1的递推公式.;9.2简单的数值方