创新实验之脉冲微分方程的数值解实验报告技巧.doc

编号： 数学与计算科学学院 创新性实验结题报告书 实验题目： 脉冲微分方程的数值解 实验成绩（教师填写）： 学 院： 数学与计算科学学院 专 业： 数学与应用数学 学生姓名： 小组成员： 指导教师： 2014年 10月 17日 摘 要：脉冲微分方程在生态学、环境科学、电力工程及自动控制系统的实时仿真等科学与工程应用领域中有广泛的应用。但常见的脉冲微分方程中,只有极少数能够获得理论解的解析表达式,而许多实际问题不需要求解脉冲微分方程的解析解,只需要数值解,因此利用数值解法求解实际问题具有十分重要的理论意义及实用价值。 关键字：脉冲微分方程、数值解 一、实验目的及主要创新性点： 本次实验我们要研究的是脉冲时滞微分方程。脉冲时滞微分方程在生态学、环境科学、电力工程及自动控制系统的实时仿真等科学与工程应用领域中有广泛的应用。但常见的脉冲时滞微分方程中，只有极少数能够获得理论解的解析表达式，而许多实际问题不需要求解脉冲时滞微分方程的解析解，只需要数值解，因此利用数值解法求解实际问题具有十分重要的理论意义及实用价值。 实验要求： 1、了解脉冲微分方程的原理及其应用； 2、使用数学软件来求解脉冲时滞微分方程的数值解。 三、实验原理： 基本概念：(1)考虑微分系统， ， 其中，为开集，为n维欧式空间,; 2）对任意存在两个集合； 3）算子任意。 设是系统,满足的解。它具有如下特点：动点 起始于并沿曲线运动，在处，与集合相遇，这时算子将点作用到，其中。而点从出发继续沿解曲线运动，直到时刻再次遇到集合,点被作用到,其中。同前面一样，沿着系统的解曲线继续 运动，如果系统的解存在，则一直运动下去。将具有上述运动过程的综合起来称为脉冲微分方程系统，有所描述的曲线称为积分曲线，积分曲线所代表的函数称为脉冲微分系统的解，与集合相遇的时刻称为脉冲时刻。假设脉冲微分系统的解在脉冲时刻是左连续的，即上述的脉冲微分方程的解有以下几个特点： 微分方程的解是连续函数，若积分弧度与没有交集或者交其不动点； 微分方程的解是分段连续函数并且有有限个第一类间断点，积分弧和相交同时不是的不动点； 微分方程的解是分段连续函数并且第一类间断点的个数是可数的，积分弧与相交且这些点不是的不动点。 1.脉冲时滞微分方程初值问题的数值解法 考虑如下脉冲时滞微分方程初值问题: ， （1） 求初值问题(1)的数值解,就是寻求准确解在一系列离散节点上的近似值其中:脉冲时刻{,}满足,时,；是脉冲算子;,, 是正整数的集合；是时滞；是初始值函数. 显然,只有当初值问题(1)的解存在且唯一时,使用数值解法才有意义.当是连续且满足李普希茨条件时,李普希茨常数为,即对任意恒有时初值问题的解存在并且唯一. 2.基于Runge-Kutta法的脉冲时滞微分方程初值问题的数值解 龙格库塔法(Runge-Kutta)是一种微分方程近似求解的数值方法.龙格库塔法通过计算不同点上的函数值,并对这些函数值作线性组合,构造近似公式,再将近似公式与精确解的泰勒展开式进行比较,使前面的若干项相同,从而使近似公式的局部截断误差阶数尽可能地高.用个函数的值表示的Runge-Kutta方法称为级Runge-Kutta方法.截断误差为的称为阶Runge-Kutta方法。四级四阶龙格库塔法即其截断误差为,每步都要计算4次的值,其标准公式即古典Runge-Kutta法为: (3) 其中表示计算过程中选取的步长,通常采用等距节点.为了得到系统的精确 一般步长选的尽可能小.若,则将区间由右向左以步长h进行分割.设 则 代表网域,其中元素称为节点。(3)式中的表示点处的斜率；表示利用求得的点处的斜率；表示利用求得的点处的斜率；表示利用求得的点处的斜率。 （3）式是常微分方程的数值解的龙格库塔法.而所谓时滞,就是一种时间差,是时间延迟,是控制系统普遍存在的一种现象,与始点和终点都有关,即某一时刻系统(1)的数值解与系统(1)在时刻及时刻的数值解都有关。将(3)式应用于方程(1),得 (4) 即当计算的值时,需要用到的值,由于的值也在网格点上,因此将其代替时滞项的值进行迭代运算。 所谓脉冲,就是在短时间内突变,随后又迅速返回其初始值的物理量,