分数阶微分方程的矩阵级数解

分数阶微分方程-课件

分数阶微分方程 第三讲分数阶微分方程基本理论 一、分数阶微分方程的出现背景及研究现状 1、出现背景 分数阶微积分是关于任意阶微分和积分的理论，它与整数阶微积分是统一的，是整数阶微积分的推广。 整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受，很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题，其无论在理论分析还是数值求解方面都已有较完善的理论。但当人们进入到复杂系统和复杂现象的研究时，经典整数阶微积分方程对这些系统的描述将遇到以下问题： （1）需要构造非线性方程，并引入一些人为的经验参数和与实际不符的假设条件； （2）因材料或外界条件的微小改变就需要构造新的模型； （3）这些非线性模型无论是理论求解还是数值求解都非常繁琐。 基于以上原因，人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程，其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势，因而成为复杂力学与物理过程数学建模的重要工具之一。 2、研究现状 在近三个世纪里，对分数阶微积分理论的研究主要在数学的纯理论领域里进行，似乎它只对数学家们有用。然而在近几十年来，分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。分数阶微积分理论也受到越来越多的国内外学者的广泛关注，特别是从实际问题抽象出来的分数阶微分方程成为很多数学工作者的研究热点。随着分数阶微分方程在越来越多的科学领域里出现，无论对分数阶微分方程的理论分析还是数值计算的研究都显得尤为迫切。然而由于分数阶微分是拟微分算子，它的保记忆性（非局部性）对现实问题进行了优美刻画的同时，也给我们的分析和计算造成很大困难。 在理论研究方面，几乎所有结果全都假定了满足李氏条件，而且证明方法也和经典微积分方程一样，换句话说，这些工作基本上可以说只是经典微积分方程理论的一个延拓。对分数阶微分方程的定性分析很少有系统性的结果，大多只是给出了一些非常特殊的方程的求解，且常用的求解方法都是具有局限性的。 在数值求解方面，现有分数阶方程数值算法还很不成熟，主要表现为： （1）在数值计算中一些挑战性难题仍未得到彻底解决，如长时间历程的计算和大空间域的计算等； （2）成熟的数值算法比较少，现在研究较多的算法主要集中在有限差分方法与有限单元法； （3）未形成成熟的数值计算软件，严重滞后于应用的需要。

一阶常微分方程的奇解

摘要.................................................... 错误!未定义书签。 1.何谓奇解.............................................. 错误!未定义书签。 2.奇解的产生............................................ 错误!未定义书签。 3.包络跟奇解的关系...................................... 错误!未定义书签。 4.理论上证明C-判别曲线与P-判别曲线方法................. 错误!未定义书签。 克莱罗微分方程 ..................................... 错误!未定义书签。 5.奇解的基本性质........................................ 错误!未定义书签。 定理１ ............................................. 错误!未定义书签。 定理２ ............................................. 错误!未定义书签。 定理３ ............................................. 错误!未定义书签。 6.小结.................................................. 错误!未定义书签。参考文献：.............................................. 错误!未定义书签。

二次微分方程的通解

教学目的：使学生掌握二阶常系数齐次线性微分方程的解法，了解二阶常系数非齐 次线性微分方程的解法 教学重点：二阶常系数齐次线性微分方程的解法 教学过程： 一、二阶常系数齐次线性微分方程 二阶常系数齐次线性微分方程 方程 y py qy 0 称为二阶常系数齐次线性微分方程 其中p 、q 均为常数 如果y 1、y 2是二阶常系数齐次线性微分方程的两个线性无关解 那么y C 1y 1C 2y 2就是它的通解 我们看看 能否适当选取r 使y e rx 满足二阶常系数齐次线性微分方程 为此将 y e rx 代入方程 y py qy 0 得 (r 2 pr q )e rx 0 由此可见 只要r 满足代数方程r 2 pr q 0 函数y e rx 就是微分方程的解 特征方程 方程r 2 pr q 0叫做微分方程y py qy 0的特征方程 特征方程 的两个根r 1、r 2可用公式 2 422,1q p p r -±+-= 求出 特征方程的根与通解的关系 (1)特征方程有两个不相等的实根r 1、r 2时 函数x r e y 11=、x r e y 22=是方程的两个线性无 关的解 这是因为

函数x r e y 11=、x r e y 22=是方程的解 又x r r x r x r e e e y y )(212121-==不是常数 因此方程的通解为 x r x r e C e C y 2121+= (2)特征方程有两个相等的实根r 1r 2时 函数x r e y 11=、x r xe y 12=是二阶常系数齐次线性微分方程的两个线性无关的解 这是因为 x r e y 11=是方程的解 又 x r x r x r x r x r x r qxe e xr p e xr r xe q xe p xe 111111)1()2()()()(1211++++=+'+'' 0 )()2(121111=++++=q pr r xe p r e x r x r 所以x r xe y 12=也是方程的解 且x e xe y y x r x r ==1112不是常数 因此方程的通解为 x r x r xe C e C y 1121+= (3)特征方程有一对共轭复根r 1, 2i 时 函数y e ( i )x 、y e (i )x 是微分方程的 两个线性无关的复数形式的解 函数y e x cos x 、y e x sin x 是微分方程的两个线性无关 的实数形式的解 函数y 1e ( i )x 和y 2e (i )x 都是方程的解 而由欧拉公式 得 y 1e ( i )x e x (cos x i sin x ) y 2e ( i )x e x (cos x i sin x ) y 1y 22e x cos x ) (2 1cos 21y y x e x +=βα y 1y 22ie x sin x ) (21sin 21y y i x e x -=βα 故e x cos x 、y 2e x sin x 也是方程解 可以验证 y 1e x cos x 、y 2e x sin x 是方程的线性无关解 因此方程的通解为

MATLAB求解常微分方程数值解

利用MATLAB求解常微分方程数值解

目录 1. 内容简介 (1) 2. Euler Method(欧拉法)求解 (1) 2.1. 显式Euler法和隐式Euler法 (2) 2.2. 梯形公式和改进Euler法 (3) 2.3. Euler法实用性 (4) 3. Runge-Kutta Method(龙格库塔法)求解 (5) 3.1. Runge-Kutta基本原理 (5) 3.2. MATLAB中使用Runge-Kutta法的函数 (7) 4. 使用MATLAB求解常微分方程 (7) 4.1. 使用ode45函数求解非刚性常微分方程 (8) 4.2. 刚性常微分方程 (9) 5. 总结 (9) 参考文献 (11) 附录 (12) 1. 显式Euler法数值求解 (12) 2. 改进Euler法数值求解 (12) 3. 四阶四级Runge-Kutta法数值求解 (13) 4.使用ode45求解 (14)

1.内容简介 把《高等工程数学》看了一遍，增加对数学内容的了解，对其中数值解法比较感兴趣，这大概是因为在其它各方面的学习和研究中经常会遇到数值解法的问题。理解模型然后列出微分方程，却对着方程无从下手，无法得出精确结果实在是让人难受的一件事情。 实际问题中更多遇到的是利用数值法求解偏微分方程问题，但考虑到先从常微分方程下手更为简单有效率，所以本文只研究常微分方程的数值解法。把一个工程实际问题弄出精确结果远比弄清楚各种细枝末节更有意思，因此文章中不追求非常严格地证明，而是偏向如何利用工具实际求解出常微分方程的数值解，力求将课程上所学的知识真正地运用到实际方程的求解中去，在以后遇到微分方程的时候能够熟练运用MATLAB得到能够在工程上运用的结果。 文中求解过程中用到MATLAB进行数值求解，主要目的是弄清楚各个函数本质上是如何对常微分方程进行求解的，对各种方法进行MATLAB编程求解，并将求得的数值解与精确解对比，其中源程序在附录中。最后考察MATLAB中各个函数的适用范围，当遇到实际工程问题时能够正确地得到问题的数值解。 2.Euler Method(欧拉法)求解 Euler法求解常微分方程主要包括3种形式，即显式Euler法、隐式Euler法、梯形公式法，本节内容分别介绍这3种方法的具体内容，并在最后对3种方法精度进行对比，讨论Euler法的实用性。 本节考虑实际初值问题 使用解析法，对方程两边同乘以得到下式

二阶线性常微分方程的幂级数解法

二阶线性常微分方程的幂级数解法 从微分方程学中知道，在满足某些条件下，可以用幂级数来表示一个函数。因此，自然想到，能否用幂级数来表示微分方程的解呢？ 例1、求方程 ''0y xy -=的通解 解：设2012n n y a a x a x a x =+++++…… 为方程的解，这里(0,1,2,,,)i a i n =……是待定常系数，将它对x 微分两次，有 ''212312132(1)(1)n n n n y a a x n n a x n na x --+=?+?++-+++ 将y ，'y 的表达式代入方程，并比较的同次幂的系数，得到 x -∞<<∞2210a ?=，30320,a a ?-= 41430,a a ?-= 52540,a a ?-= 或一般的可推得 32356(31)3k a a k k = ?????-? ， 1 3134673(31) k a a k k += ??????+ ， 320k a += 其中1a ，2a 是任意的，因而代入设的解中可得： 36347 01[1][] 2323562356(31)33434673(31) n x x x x x y a a x n n n n =+++++++++?????????-????????+ 这个幂级数的收敛半径是无限大的，因而级数的和（其中包括两个任意常数0a 及1a ）便是所要求的通解。

例6 求方程'''240y xy y --=的满足初值条件(0)0y =及'(0)1y =的解。 解 设级 2012n n y a a x a x a x =+++++……为方程的解。首先，利用初值 条件，可以得到 00a =， 11a =， 因而 2323'2123''223123232(1)n n n n n n y x a x a x a x y a x a x na x y a a x n n a x --=+++++=+++++=+?++-+ 将y ，'y ，''y 的表达式带入原方程，合并x 的各同次幂的项，并令各项系数等于零，得到 21422 0,1,0,,,1 n n a a a a a n -==== - 因而 567891111 ,0,,0,,2!63!4! a a a a a = ===== 最后得 21111 (1)!! k a k k k += ?=- , 20k a =, 对一切正整数k 成立。 将i a (0,1,2,)i = 的值代回2012n n y a a x a x a x =+++++……就得到 521 3 2!! k x x y x x k +=+++++ 2 422 (1),2!! k x x x x x xe k =++++ += 这就是方程的满足所给初值条件的解。 是否所有方程都能按以上方式求出其幂级数解？或者说究竟方程应该满足什么条件才能保证它的解可用幂级数来表示呢？级数的

总结一阶常微分方程奇解的求法

总结一阶微分方程奇解的求法 摘要：利用有关奇解的存在定理，总结出求一阶微分方程奇解的几种方法，并通过一些具体的例题说明这几种方法的应用 Using relevant theorems to develop several methods of finding singular solution of ordinary differential equation. In addition, illustrate the application of these methods through the concrete examples. 关键词：常微分方程 奇解 c-判别式 p-判别式 方法一：利用c-判别式求奇解 设一阶微分方程0, ,=?? ? ?? dx dy y x F ① 可求出方程①的通解为()0,,=c y x φ ② 如果()()???==0 ,,0,,' c y x c y x c φφ ③ 是微分方程①的解，且对③式满足：()()02 '2 '≠+y x φφ ④ 则③是微分方程①的奇解，且是通解②的包络。 例1：方程() 2 2 2 x x y dy dx dy dx + -= 的奇解 解：首先，本具题意求出该微分方程的通解为2 2 2 c cx y x ++= 与4 2 x y = 其中c 为任意常数 当时2 2 2 c cx y x ++= ， ()y c cx x c y x -++= 2 2 2 ,,φ 其相应的c -判别式为 ? ??=+=-++02022x 2 c x y c cx 易得到： ? ??=-=2 2c y c x

代入原微分方程，可知? ??=-=2 2c y c x 不是原微分方程的解； 当4 2 x y = 时，易求出2 ,1''x y x ==φφ，则有()()02 '2 '≠+y x φφ 故4 2 x y = 为原微分方程的奇解 例2：试求微分方程() () y y dy dx 9 42 2 1= -的奇解 解：首先，根据题意求出微分方程的通解为：()()0322=---y y c x 其中c 为任意常数 再由相应的c-判别式： ()()()? ??=--=---020 322c x y y c x 易求出：? ??==0y c x 或 ???==3y c x 当???==0y c x 时，代入原微分方程成立； 所以? ??==0y c x 为原微分方程的解 且有()02'=--=c x x φ；()()93232 '-=---=y y y y φ 满足(Φ‘ x )2 +(Φ‘ y )2≠0 易验证???==3y c x 不是原微分方程的解 故x=c, y=0 是元微分方程的奇解。 方法二：利用p-判别法求奇解 在微分方程①中，设y ′=p,则此方程的p-判别式为： ()()?????==0,,0 ,,' p y x F p y x F p ⑤ 消去p 之后得到的函数y=?(x)是微分方程①身为解，

一阶常微分方程解法总结

页脚内容1 第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程： ①、形如 )()(y g x f dx dy = 当0)(≠y g 时，得到 dx x f y g dy )()(=，两边积分即可得到结果； 当0)(0=ηg 时，则0)(η=x y 也是方程的解。 例1.1、xy dx dy = 解：当0≠y 时，有xdx y dy =，两边积分得到)(2ln 2为常数C C x y += 所以)(11212 C x e C C e C y ±==为非零常数且 0=y 显然是原方程的解； 综上所述，原方程的解为)(1212 为常数C e C y x = ②、形如0)()()()(=+dy y Q x P dx y N x M 当0)()(≠y N x P 时，可有dy y N y Q dx x P x M ) ()()()(=，两边积分可得结果； 当0)(0=y N 时，0y y =为原方程的解，当0(0=） x P 时，0x x =为原方程的解。 例1.2、0)1()1(22=-+-dy x y dx y x

页脚内容2 解：当0)1)(1(22≠--y x 时，有dx x x dy y y 1 122-=-两边积分得到 )0(ln 1ln 1ln 22≠=-+-C C y x ，所以有)0()1)(1(22≠=--C C y x ； 当0)1)(1(22=--y x 时，也是原方程的解； 综上所述，原方程的解为)()1)(1(22为常数C C y x =--。 ⑵可化为变量可分离方程的方程： ①、形如)(x y g dx dy = 解法：令x y u = ，则udx xdu dy +=，代入得到)(u g u dx du x =+为变量可分离方程，得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x x y f =。 ②、形如)0(),(≠+=ab by ax G dx dy 解法：令by ax u +=，则b du adx dy +=，代入得到)(1u G b a dx du b =+为变量可分离方程，得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x by ax f =+。 ③、形如 )(222111c y b x a c y b x a f dx dy ++++= 解法：01、02211 =b a b a ，转化为)(by ax G dx dy +=，下同①； 02、0221 1 ≠b a b a ，???=++=++00222111c y b x a c y b x a 的解为),(00y x ，令???-=-=00y y v x x u

二阶常微分方程解

第七节 二阶常系数线性微分方程 在上节我们已经讨论了二阶线性微分方程解的结构，二阶线性微分方程的求解问题，关键在于如何求二阶齐次方程的通解和非齐次方程的一个特解。本节讨论二阶线性方程的一个特殊类型，即二阶常系数线 性微分方程及其求解方法。先讨论二阶常系数线性齐 §7.1 二阶常系数线性齐次方程及其求 22dx y d ＋p dx dy ＋qy ＝ 0 (7.1) 其中p 、q 是常数，由上节定理二知，要求方程(7.1)的通解，只要求出其任意两个线性无关的特解y 1,y ２ 我们先分析方程(7.1)可能具有什么形式的特解， 从方程的形式上来看，它的特点是22 dx y d ，dx dy ，y 各乘 以常数因子后相加等于零，如果能找到一个函数y ，

其22dx y d ，dx dy ，y 之间只相差一个常数因子，这样的函 数有可能是方程(7.1)的特解，在初等函数中，指数函数e rx y ＝e rx (其中r 为待定常数) 将y ＝e rx ，dx dy ＝re rx ，22dx y d ＝r 2e rx 代入方程 (7.1) 得 r 2e rx ＋pre rx ＋qe rx ＝ 0 或 e rx （r 2＋pr ＋q ）＝ 因为e rx ≠ 0 r 2 ＋pr ＋q ＝ 由此可见，若 r r 2＋pr ＋q ＝ 0 (7.2) 的根，那么e rx 就是方程(7.1)的特解，于是方程(7.1)的求解问题，就转化为求代数方程(7.2)的根问题。称(7.2)式为微分方程(7.1) 特征方程(7.2)是一个以r 为未知函数的一元二次代数方程。特征方程的两个根r １，r 2，称为特征根，由代数知识，特征根r 1，r 2 有三种可能的情况，下面 (1)若特证方程(7.2)有两个不相等的实根r １，r 2，此时e r １x ，e r2x 是方程(7.1)

用MATLAB解常微分方程

实验四 求微分方程的解 一、问题背景与实验目的 实际应用问题通过数学建模所归纳而得到的方程，绝大多数都是微分方程，真正能得到代数方程的机会很少．另一方面，能够求解的微分方程也是十分有限的，特别是高阶方程和偏微分方程（组）．这就要求我们必须研究微分方程（组）的解法，既要研究微分方程（组）的解析解法（精确解），更要研究微分方程（组）的数值解法（近似解）． 对微分方程（组）的解析解法(精确解)，Matlab 有专门的函数可以用，本实验将作一定的介绍． 本实验将主要研究微分方程(组)的数值解法（近似解），重点介绍 Euler 折线法． 二、相关函数（命令）及简介 1．dsolve('equ1','equ2',…)：Matlab 求微分方程的解析解．equ1、equ2、…为方程（或条件）．写方程（或条件）时用 Dy 表示y 关于自变量的一阶导数，用用 D2y 表示 y 关于自变量的二阶导数，依此类推． 2．simplify(s)：对表达式 s 使用 maple 的化简规则进行化简． 例如： syms x simplify(sin(x)^2 + cos(x)^2) ans=1 3．[r,how]=simple(s)：由于 Matlab 提供了多种化简规则，simple 命令就是对表达式 s 用各种规则进行化简，然后用 r 返回最简形式，how 返回形成这种形式所用的规则． 例如： syms x [r,how]=simple(cos(x)^2-sin(x)^2) r = cos(2*x) how = combine 4．[T,Y] = solver(odefun,tspan,y 0) 求微分方程的数值解． 说明： (1) 其中的 solver 为命令 ode45、ode23、ode113、ode15s 、ode23s 、ode23t 、ode23tb 之一． (2) odefun 是显式常微分方程：?????==0 0)() ,(y t y y t f dt dy (3) 在积分区间 tspan =],[0f t t 上，从0t 到f t ，用初始条件0y 求解．

高阶线性微分方程常用解法介绍

高阶线性微分方程常用解法简介 关键词：高阶线性微分方程 求解方法 在微分方程的理论中，线性微分方程是非常值得重视的一部分内容，这不仅 因为线性微分方程的一般理论已被研究的十分清楚，而且线性微分方程是研究非线性微分方程的基础，它在物理、力学和工程技术、自然科学中也有着广泛应用。下面对高阶线性微分方程解法做一些简单介绍. 讨论如下n 阶线性微分方程：1111()()()()n n n n n n d x d x dx a t a t a t x f t dt dt dt ---++++= （1），其中()i a t （i=1,2,3,,n ）及f(t)都是区间a t b ≤≤上的连续函数，如果 ()0f t ≡，则方程（1）变为 1111()()()0n n n n n n d x d x dx a t a t a t x dt dt dt ---++++= （2），称为n 阶齐次线性微分方程，而称一般方程（1）为n 阶非齐次线性微分方程，简称非齐次线性微分方程，并且把方程（2）叫做对应于方程（1）的齐次线性微分方程. 1.欧拉待定指数函数法 此方法又叫特征根法，用于求常系数齐次线性微分方程的基本解组。形如 111121[]0,(3),n n n n n n n d x d x dx L x a a a x dt dt dt ---≡++++=其中a a a 为常数，称为n 阶常系数齐次线性微分方程。 111111111111[]()()()n t n t t t t n n n n n n n t t n n n n n n n d e d e de L e a a a e dt dt dt a a a e F e F a a a n λλλλλλλλλλλλλλλλ---------≡++++=++++≡≡++++其中=0(4)是的次多项式. ()F λ为特征方程，它的根为特征根. 1.1特征根是单根的情形 设12,,,n λλλ是特征方程111()0n n n n F a a a λλλλ--≡++++=的n 个彼此不相等的根，则应相应地方程（3）有如下n 个解：12,,,.n t t t e e e λλλ（5）我们指出这n 个解在区间a t b ≤≤上线性无关，从而组成方程的基本解组. 如果(1,2,,)i i n λ=均为实数，则（5）是方程（3）的n 个线性无关的实值 解，而方程（3）的通解可表示为1212,n t t t n x c e c e c e λλλ=+++其中12,,,n c c c 为任意常数. 如果特征方程有复根，则因方程的系数是实常数，复根将称对共轭的出现.设1i λαβ=+是一特征根，则2i λαβ=-也是特征根，因而于这对共轭复根

一阶常微分方程解法总结

第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程： ①、形如 )()(y g x f dx dy = 当0)(≠y g 时，得到 dx x f y g dy )() (=，两边积分即可得到结果； 当0)(0=ηg 时，则0)(η=x y 也是方程的解。 例1.1、 xy dx dy = 解：当0≠y 时，有xdx y dy =，两边积分得到)(2ln 2为常数C C x y += 所以)(112 12 C x e C C e C y ±==为非零常数且 0=y 显然是原方程的解； 综上所述，原方程的解为)(12 12 为常数C e C y x = ②、形如0)()()()(=+dy y Q x P dx y N x M 当0)()(≠y N x P 时，可有 dy y N y Q dx x P x M ) () ()()(=，两边积分可得结果； 当0)(0=y N 时，0y y =为原方程的解，当0(0=） x P 时，0x x =为原方程的解。 例1.2、0)1()1(22=-+-dy x y dx y x 解：当0)1)(1(22≠--y x 时，有 dx x x dy y y 1 122-=-两边积分得到 )0(ln 1ln 1ln 22≠=-+-C C y x ，所以有)0()1)(1(22≠=--C C y x ； 当0)1)(1(22=--y x 时，也是原方程的解； 综上所述，原方程的解为)()1)(1(22为常数C C y x =--。 ⑵可化为变量可分离方程的方程：

①、形如 )(x y g dx dy = 解法：令x y u =，则udx xdu dy +=，代入得到)(u g u dx du x =+为变量可分离方程，得 到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x x y f =。 ②、形如)0(),(≠+=ab by ax G dx dy 解法：令by ax u +=，则b du adx dy +=，代入得到)(1u G b a dx du b =+为变量可分离方程， 得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x by ax f =+。 ③、形如 )(2 221 11c y b x a c y b x a f dx dy ++++= 解法：0 1、 02 2 11=b a b a ，转化为 )(by ax G dx dy +=，下同①； 02、 022 1 1≠b a b a ，???=++=++00 222111 c y b x a c y b x a 的解为),(00y x ，令???-=-=0 0y y v x x u 得到，)()( )(221 12211u v g u v b a u v b a f v b u a v b u a f du dv =++=++=，下同②； 还有几类：xy u dy xy xg dx xy yf ==+,0)()( 以上都可以化为变量可分离方程。 例2.1、 2 5--+-=y x y x dx dy 解：令2--=y x u ，则du dx dy -=，代入得到u u dx du 7 1+= - ，有dx udu 7-= 所以)(72 2 为常数C C x u +-=，把u 代入得到)(72 22 为常数） （C C x y x =+--。 例2.2、 1 212+-+-=y x y x dx dy 解：由???=+-=+-012012y x y x 得到?????=-=3131y x ，令?? ???-=+=3131y v x u ，有???==du dx dv dy ，代入得到

二阶常微分方程的解法及其应用.

目录 1 引言 (1) 2 二阶常系数常微分方程的几种解法 (1) 2.1 特征方程法 (1) 2.1.1 特征根是两个实根的情形 (2) 2.1.2 特征根有重根的情形 (2) 2.2 常数变异法 (4) 2.3 拉普拉斯变化法 (5) 3 常微分方程的简单应用 (6) 3.1 特征方程法 (7) 3.2 常数变异法 (9) 3.3 拉普拉斯变化法 (10) 4 总结及意义 (11) 参考文献 (12)

二阶常微分方程的解法及其应用 摘要：本文通过对特征方程法、常数变易法、拉普拉斯变换法这三种二阶常系数常微分方程解法进行介绍，特别是其中的特征方程法分为特征根是两个实根的情形和特征根有重根的情形这两种情况，分别使用特征值法、常数变异法以及拉普拉斯变换法来求动力学方程，现今对于二阶常微分方程解法的研究已经取得了不少成就，尤其在二阶常系数线性微分方程的求解问题方面卓有成效。应用常微分方程理论已经取得了很大的成就，但是，它的现有理论也还远远不能满足需要，还有待于进一步的发展，使这门学科的理论更加完善。 关键词：二阶常微分方程;特征分析法；常数变异法；拉普拉斯变换

METHODS FOR TWO ORDER ORDINARY DIFFERENTIAL EQUATION AND ITS APPLICATION Abstract:This paper introduces the solution of the characteristic equation method, the method of variation of parameters, the Laplasse transform method the three kind of two order ordinary differential equations with constant coefficients, especially the characteristic equation method which is characteristic of the root is the two of two real roots and characteristics of root root, branch and don't use eigenvalue method, method of variation of constants and Laplasse transform method to obtain the dynamic equation, the current studies on solution of ordinary differential equations of order two has made many achievements, especially in the aspect of solving the problem of two order linear differential equation with constant coefficients very fruitful. Application of the theory of ordinary differential equations has made great achievements, however, the existing theory it is still far from meeting the need, needs further development, to make the discipline theory more perfect. Keywords:second ord er ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform 1 引言 数学发展的历史告诉我们，300年来数学分析是数学的首要分支，而微分方程

分数阶微积分发展现状及展望教学文稿

分数阶微积分发展现状及展望 在数学领域中，大体分为五种研究方向：基础数学，应用数学，计 算数学，概率论与数理统计，统计学与控制论。这五个方向对数学在当 代的发展都有不可或缺的作用。从研究内容来讲，方程、算子、群论、 图论、代数、几何等等都是数学领域重要的研究对象。作为基础数学专 业分数阶微分方程方向的博士生，本文将从分数阶微分方程的发展的历 史及现状、本人对分数阶微分方程未来发展的看法来介绍分数阶微分的 基本知识。 (一)、发展历史及现状 牛顿和莱布尼兹发明的微积分是现代数学与古典数学的分水岭。分数阶微积分是关于任意阶微分和积分的理论，它与整数阶微积分是统一的，是整数阶微积分的推广。整数阶微积分作为描述经典物理及相关学科理论的解析数学工具已为人们普遍接受，很多问题的数学模型最终都可以归结为整数阶微分方程的定解问题，其无论在理论分析还是数值求解方面都已有了比较完善的理论。但当人们进入到复杂系统和复杂现象的研究时，经典整数阶微积分方程对这些系统的描述将遇到一些问题，如：需要构造非线性方程，并引入一些人为的经验参数和与实际不符的假设条件；因材料或外界条件的微小改变就需要构造新的模型等等。基于以上原因，人们迫切期待着有一种可用的数学工具和可依据的基本原理来对这些复杂系统进行建模。分数阶微积分方程非常适合于刻画具有记忆和遗传性质的材料和过程，其对复杂系统的描述具有建模简单、参数物理意义清楚、描述准确等优势，因而成为复杂力学与物理过程数学建模的重要工具之一。 对大多数研究人员和工程师而言，分数阶微积分也许还是比较陌生的，但它实际上早在300多年前就被提出。1695年9月，洛必达 （L’Hospital）在给莱布尼兹的著名信件中就写到“对于简单的线性函数 f(x)=x，如果函数导数次数为分数而不是整数那会怎样”。这是公认的第一次提及分数阶微分。1832年，刘维尔（Liouville）成功的应用了自己提出的分数阶导数的定义，解决了势理论问题。之后刘维尔发表的一系列文 章使他成为分数阶微积分理论的实际级创始人。1974年，Oldham与Spanier出版了第一本关于分数阶微积分理论的专著。 在近三个世纪里，对分数阶微积分理论的研究主要在数学的纯理论领域里进行，但是从近几十年，分数阶微分方程越来越多的被用来描述光学和热学系统、流变学及材料和力学系统、信号处理和系统识别、控制和机器人及其他应用领域中的问题。分数阶微积分理论也受到越来越多的国内

二阶常微分方程解

二阶常微分方程解

————————————————————————————————作者: ————————————————————————————————日期:

第七节 二阶常系数线性微分方程 的解法 在上节我们已经讨论了二阶线性微分方程解的结构，二阶线性微分方程的求解问题,关键在于如何求二阶齐次方程的通解和非齐次方程的一个特解。本节讨论二阶线性方程的一个特殊类型，即二阶常系数线性微分方程及其求解方法。先讨论二阶常系数线性齐次方程的求解方法。 §7.1 二阶常系数线性齐次方程及其求解方法 设给定一常系数二阶线性齐次方程为 ?? 22 dx y d +p dx dy ＋ｑy=０ (7.1) 其中p 、q 是常数,由上节定理二知,要求方程(7．1)的通解，只要求出其任意两个线性无关的特解y 1,ｙ2就可以了，下面讨论这样两个特解的求法。 我们先分析方程(7.1)可能具有什么形式的特解, 从方程的形式上来看，它的特点是22dx y d ，dx dy ,y 各乘以 常数因子后相加等于零，如果能找到一个函数ｙ,其

22dx y d ，dx dy ，ｙ之间只相差一个常数因子，这样的函数有可能是方程(７.1)的特解,在初等函数中,指数函数e rx ,符合上述要求，于是我们令 y=e r ｘ (其中r 为待定常数）来试解 将y ＝e rx ,dx dy =re r ｘ,22dx y d ＝r ２e r ｘ 代入方程(7.1) 得 r 2e ｒx ＋ｐre rx ＋ｑｅrx ＝０ 或 e r ｘ(r ２+pr+q ）＝0 因为e rx ≠０,故得 ? r 2 +pr ＋q=0 由此可见，若r 是二次方程 ?? r 2＋pr ＋q=０ (7．2) 的根，那么e r ｘ就是方程(7.１)的特解,于是方程（７.1)的求解问题，就转化为求代数方程(7.2）的根问题。称（7.2)式为微分方程(７.１)的特征方程。 特征方程(7.2)是一个以r 为未知函数的一元二次代数方程。特征方程的两个根r 1，r 2，称为特征根,由代数知识，特征根r １，r 2有三种可能的情况,下面我们分别进行讨论。 (1)若特证方程（7．2）有两个不相等的实根r 1， r 2，此时e r １ｘ ,e r2x 是方程（７.1)的两个特解。

Matlab求解微分方程(组)及偏微分方程(组)

第四讲 Matlab 求解微分方程（组） 理论介绍：Matlab 求解微分方程（组）命令 求解实例：Matlab 求解微分方程（组）实例 实际应用问题通过数学建模所归纳得到的方程，绝大多数都是微分方程，真正能得到代数方程的机会很少.另一方面，能够求解的微分方程也是十分有限的，特别是高阶方程和偏微分方程（组）.这就要求我们必须研究微分方程（组）的解法：解析解法和数值解法. 一．相关函数、命令及简介 1.在Matlab 中，用大写字母D 表示导数，Dy 表示y 关于自变量的一阶导数，D2y 表示y 关于自变量的二阶导数，依此类推.函数dsolve 用来解决常微分方程（组）的求解问题，调用格式为： X=dsolve(‘eqn1’,’eqn2’,…) 函数dsolve 用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解. 注意，系统缺省的自变量为t 2.函数dsolve 求解的是常微分方程的精确解法，也称为常微分方程的符号解.但是，有大量的常微分方程虽然从理论上讲，其解是存在的，但我们却无法求出其解析解，此时，我们需要寻求方程的数值解，在求常微分方程数值解方面，MATLAB 具有丰富的函数，我们将其统称为solver ，其一般格式为： [T,Y]=solver(odefun,tspan,y0) 说明：(1)solver 为命令ode45、ode23、ode113、ode15s 、ode23s 、ode23t 、ode23tb 、ode15i 之一. (2)odefun 是显示微分方程'(,)y f t y =在积分区间tspan 0[,]f t t =上从0t 到f t 用初始条件0y 求解. (3)如果要获得微分方程问题在其他指定时间点012,,, ,f t t t t 上的解，则令 tspan 012[,,,]f t t t t =(要求是单调的). (4)因为没有一种算法可以有效的解决所有的ODE 问题，为此，Matlab 提供

分数阶微分方程数值解的一种逼近方法.

分数阶微分方程数值解的一种逼近方法 By:Pankaj Kumar, Om Prakash Agrawal 摘要 本文提出了一类分数阶微分方程（FDEs）的数值解方案.在这种方法中,FDEs 被Caputo型分数阶导数所表现. Caputo型分数阶导数的属性可以让一个分数阶微分方程减弱为一个Volterra型积分方程. 这样做了之后,许多研究Volterra 型积分方程的数值方法也可以应用于寻找FDEs的数值解. 本文总时间被划分为一组小区间,在两个连续区间中，用二次多项式逼近未知函数. 这些近似被替换成转化的Volterra型积分方程由此获得一组方程. 这些方程的解提供了FDE的解. 这种方法被应用于解决两种类型的FDEs,线性和非线性. 用这里给出的方法得到的解能与解析解和其他方法的数值解较好的吻合. 同时结果说明这种数值方法是稳定的. 1.引言 本文讨论分数阶微分方程的数值解. 分数阶导数和分数阶积分近年来收到了广泛的关注. 在许多实际应用中，分数阶导数和分数阶积分为考虑的系统提供了更加精确地模型. 比如，分数阶导数已经被成功地运用到模拟许多粘性材料的依赖频率的阻尼行为.1980年之前，Bagley 和Torvik提出了这个领域已经被研究的工作的一个回顾，并且说明了半阶导数模型可以非常好地描述阻尼材料的频率以来. 另一些学者说明了分数阶导数和分数阶积分在电化学过程，电解质极化，有色噪声，粘性材料和混沌领域的应用. Mainardi，Rossikhin和Shitikova 提出了分数阶导数和分数阶积分在一般固体力学，特定粘弹性阻尼模型中的应用的调查. Magin提出了分数阶微积分在生物工程的三个关键部分的回顾. 分数阶导数和分数阶积分在其他领域的应用以及相关的数学工具和技巧还可以在许多其他文献上找到. 系统模型中分数阶导数的引进大多会导致分数阶微分方程的出现. 对某些特定的分数阶微分方程在通常系统条件下的解，已经有几种方法被找到. 这些方法包括，拉普拉斯变换，傅里叶变换，模态综合法和特征向量展开法，数值法以

一阶常微分方程的奇解

摘要 (2) 1.何谓奇解 (2) 2.奇解的产生 (3) 3.包络跟奇解的关系 (4) 4.理论上证明C-判别曲线与P-判别曲线方法 (5) 4.1 克莱罗微分方程 (9) 5.奇解的基本性质 (12) 5.1 定理１ (12) 5.2 定理２ (14) 5.3 定理３ (14) 6.小结 (14) 参考文献： (15)

一阶常微分方程的奇解 摘要 在常微分方程中，我们知道方程的解可以有多种，现在我们来讨论求奇解的方法。我们看到某些微分方程，会存在一些特殊的积分曲线，他并不属于这方程的积分曲线族，但是，在这些特殊的积分曲线上的每一点处，都有积分曲线族中的一条曲线和他在此处相切。从而我们引出了积分曲线族的包络，而为了求微分方程的奇解，，我们应先求出他的通解，然后求通解的包络。 关键词：奇解，包络，C-判别式，P-判别式 1.何谓奇解 设一阶隐式方程),,(,y y x F =0有一特解

)(:x y ψ=Γ，j x ∈ 如果对每一点Γ∈P ，在P 点的任何一个领域，方程),,(,y y x F =0都有一个不同于Γ的解在P 点与Γ相切，则称Γ是微分方程的),,(,y y x F =0的奇解 定义：如果一个一阶微分方程的一个特解的积分曲线上的每一点都至少和这个微分方程的不同的积分曲线相切，并且这相切的积分曲线在切点的任何邻域都不重合，则称这个特解为这个微分方程的奇解 2.奇解的产生 先看一个例子，求方程 033=-?? ? ??y dx dy （1） 或与它等价的方程 3y dx dy = 的解。 经分离变量后，可得（1）的通解 3)(27 1c x y += 容易看出，y=0也是原方程的一个解。现在来研究这个解y=0有什么特殊的地方。由图我们看到，在解y=0上的每一 点)0,(0x 处相切，这种特殊的积分曲线y=0 称为奇积分曲线，他所对应的解就是奇 解，这就是奇解的产生。 我们现在给出曲线族包络的定义 某些微分方程，存在一些特殊的积分 曲线，会存在一些特殊的积分曲线，他并 不属于这方程的积分曲线族，但是，在这些特殊的积分曲线上的每一点处，都有积分曲线族中的一条曲线和他在此处相切。在几何学里，这些特殊的积分曲线称为上述积分曲线族的包络，在微分方程里，这些特殊的积分曲线所对应的解酒称为方程的奇解。