常微分方程matlab求解.pptxVIP

Outline;14.1 微分方程的基本概念;14.2 几种常用微分方程类型;4.可降阶的高阶微分方程 型的微分方程：微分方程 的右端仅含有自变量 x ，容易看出，只要把 作为新的未知函数，那么微分方程 即化为新未知函数 的一 阶微分方程，两边积分，就得到一个 阶的微分方程 同理可得 依此法继续进行，接连积分 n次，便得到方程 的含有 n 个任意常数的通解。 型的微分方程:方程的右端不显含未知函数 y。如果我们设 ，那么 因此，方程 就成为 ，这是一个关于变量 的一阶微分方程 ，设其通解为 ，又 因此又得到一个一阶微分方程 对它进行积分，便得到方程 的通解为 型的微分方程：方程中不显含自变??x ，为了求出它的解，我们令 ，并利用复合函数求导法则把 化为对 的导数，即 这样，方程 就成为 这是一个关于变量 的一阶微分方程，设它的通解为 分离变量并积分，便得方程 的通解为 ;14.3 高阶线性微分方程;14.4 一阶微分方程初值问题的数值解;2. Runge-Kutta法及其MATLAB实现 考虑微分方程 ，由Lagrange微分中值定理，存在 ，使得 于是，由 得 记 ，则称 为区间 上的平均斜率。这样，只要给 出了 的一种算法，就可以得到求解微分方程初值问题的一种计算公式。显然，显式Euler公式就是以 作为平均斜率 的近似。 经典四阶Runge-Kutta方法的迭代公式： ;14.5 一阶微分方程组和高阶微分方程的数值解;14.6 边值问题的数值解;非线性方程边值问题的打靶法： 考虑二阶常微分方程 的边值问题，边界条件为 。假定 该问题可以转换为下面的初值问题 则问题转化为求解 ，这是一个复杂的超越方程，可以考虑引入牛顿迭代法求 解参数 m。具体的迭代公式为： 式中 通过这些关系可以建立方程 具体计算中可以指定一个m 值，然后求解上面的初值问题，将结果代入上面的迭代公式中迭代一步，并将结果代入上式中重新计算，直至两次计算出来的 m值的误差在允许的范围内为止，最后将 m值代入初值问题 即可求解原始问题。 ;2.边值问题的MATLAB函数求解 MATLAB能求解的边值问题的一般形式如下 其中y 为状态变量向量 为方程中其他未知参数向量。该方程已知的边界值为 MATLAB提供了专门求解边值问题的bvp解算器bvpslover。要想求解一个常微分方程的边值问题，一般应该遵循以下几个步骤： （1）参数初始化 （2）微分方程和边值问题的MATLAB函数描述 （3）边值问题的求解 ; 谢谢大家！