常微分matla实b验报告.docVIP

数学与信息科学系实训报告 常微分方程初值问题数值解法及其MATLAB实现 姓名： 学号： 专业：信息与计算科学 年级：2010级 指导教师：朱耀生 完成时间：2013年11月25—2013年12月6 实验目的: (1)研究满足给定方程的可微函数的数值方法 (2)培养matlab编程和上机调试能力 实验基本原理和内容: 根据给定的初始条件，确定常微分方程惟一解的问题叫常微分方程初值问题。大多数实际问题难以求得解析解，必须将微分问题离散化，用数值方法求其近似解。　一阶常微分方程的初值问题的提法是，求出函数，使满足条件????????????　　(1)利用数值方法解问题 (1)时，通常假定解存在且惟一，解函数及右端函数具有所需的光滑程度。数值解法的基本思想是:先取自变量一系列离散点,把微分问题(1)离散化，求出离散问题的数值解,并以此作为微分问题解的近似。例如取步长0,以剖分区间[,]，令=+,把微分方程离散化成一个差分方程。以()表微分方程初值问题的解,以表差分问题的解，就是近似解的误差，称为全局误差。因此，设计各种离散化模型,求出近似解,估计误差以及研究数值方法的稳定性和收敛性等构成了数值解法的基本内容function [x,y]=meuler(df,xspan,y0,h) %用途：改进欧拉公式解常微分方程y=f(x,y)