微分方程求解

二次微分方程的通解

第六节 二阶常系数齐次线性微分方程 教学目的：使学生掌握二阶常系数齐次线性微分方程的解法，了解二阶常系数非齐 次线性微分方程的解法 教学重点：二阶常系数齐次线性微分方程的解法 教学过程： 一、二阶常系数齐次线性微分方程 二阶常系数齐次线性微分方程 方程 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 2 1y y i x e x -= βα 故e x cos x 、y 2e x sin x 也是方程解 可以验证 y 1e x cos x 、y 2e x sin x 是方程的线性无关解

微分方程公式运用表

微分方程公式运用表 一、 一阶微分方程 判断特征： (,)dy f x y dx = 类型一：()()dy g x h y dx =（可分离变量的方程） 解法（分离变量法）： ()()dy g x dx h y =，然后两边同时积分。 类型二：()()dy P x y Q x dx +=（一阶线性方程） 解法（常数变易法）：()()(())P x dx P x dx y e C Q x e dx -??=+? 类型三： (,)(,)dy f x y f tx ty dx ==（一阶齐次性方程） 解法（换元法）：y u x =?令类型一 类型四：P()y=Q(x)y n dy x dx +（伯努利方程） 解法（同除法）：1()()n n dy y P x y Q x dx --+=?类型二 二、 可降阶的高阶微分方程 类型一：()()n y f x = 解法（多次积分法）：(1)()()n du u y f x f x dx -=? =?令多次积分求 类型二：''(,')y f x y = 解法：'(,)dp p y f x p dx =?=?令一阶微分方程 类型三：''(,')y f y y = 解法：'(,)dp dp dy dp p y p f y p dx dy dx dy =?==??令类型二 三、线性微分方程 类型一：''()'()0y P x y Q x y ++=（二阶线性齐次微分方程） 解法：找出方程的两个任意线性不相关特解：12(),()y x y x

则：1122()()()y x c y x c y x =+ 类型二：''()'()()y P x y Q x y f x ++=（二阶线性非齐次微分方程） 解法：先找出对应的齐次微分方程的通解：31122()()()y x c y x c y x =+ 再找出非齐次方程的任意特解()p y x ，则：1122()()()()p y x y x c y x c y x =++ 类型三：'''0y py q ++=（二阶线性常系数齐次微分方程） 解法（特征方程法）：2 1,20p q λλλ++=?= （一）122121240x x p q y c e c e λλλλ?=->?≠?=+ （二）12120()x y c c x e λλλλ?=?==?=+ （三）12120,(cos sin )x i i y e c x c x αλαβλαβββ?

二次微分方程的通解

教学目的：使学生掌握二阶常系数齐次线性微分方程的解法，了解二阶常系数非齐 次线性微分方程的解法 教学重点：二阶常系数齐次线性微分方程的解法 教学过程： 一、二阶常系数齐次线性微分方程 二阶常系数齐次线性微分方程 方程 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 是方程的线性无关解 因此方程的通解为

微分方程几种求解方法

第五章 控制系统仿真 §5.2 微分方程求解方法 以一个自由振动系统实例为例进行讨论。 如下图1所示弹簧-阻尼系统，参数如下： M=5 kg, b=1 N.s/m, k=2 N/m, F=1N F 图1 弹簧-阻尼系统 假设初始条件为：00=t 时，将m 拉向右方，忽略小车的摩擦阻力，m x 0)0(= s m x /0)0(=? 求系统的响应。 )用常微分方程的数值求解函数求解包括ode45、 ode23、ode113、ode15s 、ode23s 等。 wffc1.m myfun1.m 一、常微分方程的数值求解函数ode45求解 解：系统方程为 F kx x b x m =++??? 这是一个单变量二阶常微分方程。

将上式写成一个一阶方程组的形式，这是函数ode45调用规定的格式。 令： x x =)1( （位移） )1()2(? ?==x x x （速度） 上式可表示成： ??????--=??????=??? ???????)1(*4.0)2(*2.02.0)2()2()2()1(x x x x x x x && 下面就可以进行程序的编制。 %写出函数文件myfun1.m function xdot=myfun1(t,x) xdot=[x(2);0.2-0.2*x(2)-0.4*x(1)]; % 主程序wffc1.m t=[0 30]; x0=[0;0]; [tt,yy]=ode45(@myfun1,t,x0); plot(tt,yy(:,1),':b',tt,yy(:,2),'-r') hold on plot(tt,0.2-0.2*yy(:,2)-0.4*yy(:,1),'-k') legend('位移','速度',’加速度’)

偏微分方程组解法

偏微分方程组解法 某厚度为10cm 平壁原温度为20C ?，现其两侧面分别维持在20C ?和120C ?，试求经过8秒后平壁温度分布，并分析温度分布随时间的变化直至温度分布稳定为止。 22x t a t ??=??τ 式中a 为导温系数，/s m c 2；2=a 。 解： 模型转化为标准形式： 2 21x t t a ??=??τ 初始条件为： ()200,=x t 边界条件为： ()120,0=τt ，()20,1.0=τt 函数： pdefun.m %偏微分方程(一维动态传热) function [c,f,s]=pdefun(x,t,u,dudx) c=1/2e-4;f=dudx;s=0; icbun.m %偏微分方程初始条件(一维动态传热) function u0=icbun(x) u0=20; bcfun.m %偏微分方程边界条件(一维动态传热) function [pl,ql,pr,qr]=bcfun(xl,ul,xr,ur,t) pl=ul-120;ql=0;pr=ur-20;qr=0; 命令： x=linspace(0,10,20)*1e-2; t=linspace(0,15,16); sol=pdepe(0,pdefun,icfun,bcfun,x,t); mesh(x,t,sol(:,:,1)) %温度与时间和空间位置的关系图 %画1、2、4、6、8、15s 时刻温度分布图

plot(x,sol(2,:,1)) 1s时刻，（因为本题sol第一行为0时刻） hold on plot(x,sol(3,:,1)) plot(x,sol(5,:,1)) plot(x,sol(7,:,1)) plot(x,sol(9,:,1)) plot(x,sol(16,:,1)) 计算结果： %第8秒时温度分布 x sol(9,:,1) 经过8秒时的温度分布为： x/cm 0 0.5263 1.0526 1.5789 2.1053 2.6316 3.1579 t/C ?120.0000 112.5520 105.1653 97.8994 90.8100 83.9477 77.3562 x/cm 3.6842 4.2105 4.7368 5.2632 5.7895 6.3158 6.8421 t/C ?71.0714 65.1202 59.5200 54.2784 49.3930 44.8518 40.6338 x/cm 7.3684 7.8947 8.4211 8.9474 9.4737 10.0000 t/C ?36.7095 33.0419 29.5877 26.2982 23.1207 20.0000 或者求第8秒时，x=0,2,4,,6,8,10cm处的温度 [uout,duoutdx]=pdeval(0,x,sol(9,:,:),[0,2,4,6,8,10]*1e-2) 120.0000 92.2279 67.5007 47.5765 32.3511 20.0000

各种类型的微分方程及其相应解法

各种类型的微分方程及其相应解法 专业班级：交土01班 姓名：高云 学号：1201110102 微分方程的类型有很多种，解题时先判断微分方程是哪种类型，可以帮助我们更快解题，所以我们有必要归纳整理一下各类型（主要是一阶和二阶）的微分方程及其相应解法。 一、一阶微分方程的解法 1.可分离变量的方程 dx x f dy y g )()(=,或)()(y g x f dx dy = 其特点是可以把变量x 和y 只分别在等式的两边，解法关键是把变量分离后两边积分。 例1.求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项,得dx y dy x y )1()1(2-=- 设,01,012≠-≠-x y 分离变量得 dx x dy y y 1112-=- 两端积分??-=-dx x dy y y 1112得 ||ln |1|ln |1|ln 2 112C x y +-=- 于是 2212)1(1-±=-x C y 记,21C C ±=则得到题设方程的通解 .)1(122-=-x C y 2.齐次方程 （1）)(x y f dx dy = （2） )(c by ax f dx dy ++=（a ，b 均不等于0） 例2求解微分方程.2222xy y dy y xy x dx -=+- 解 原方程变形为=+--=2222y xy x xy y dx dy ,1222?? ? ??+--??? ??x y x y x y x y 令,x y u =则,dx du x u dx dy +=方程化为,1222u u u u dx du x u +--=+ 分离变量得?? ????-+--??? ??--112212121u u u u ,x dx du = 两边积分得 ,ln ln ln 2 1)2ln(23)1ln(C x u u u +=----

二阶常系数齐次线性微分方程的通解证明教学提纲

二阶常系数齐次线性微分方程的通解证明

二阶常系数齐次线性微分方程的通解证明 来源：文都教育 在考研数学中，微分方程是一个重要的章节，每年必考，其中的二阶常系数齐次线性微分方程是一个基本的组成部分，它也是求解二阶常系数非齐次线性微分方程的基础，但很多同学对其求解公式不是十分理解，做题时也感到有些困惑，为了帮助大家对其通解公式有更深的理解和更牢固的掌握，文都网校的蔡老师下面对它们进行一些分析和简捷的证明，供考研的朋友们学习参考。 一、二阶常系数齐次线性微分方程的通解分析 通解公式：设0y py qy '''++=，,p q 为常数，特征方程02=++q p λλ的特征根为 12,λλ，则 1）当12λλ≠且为实数时，通解为1212x x y C e C e λλ=+； 2）当12λλ=且为实数时，通解为1112x x y C e C xe λλ=+； 3）当12,i λλαβ=±时，通解为12(cos sin )x y e C x C x αββ=+； 证：若02=++q p λλ的特征根为12,λλ，则1212(),p q λλλλ=-+ =，将其代入方程0y py qy '''++=中得1212()y py qy y y y λλλλ''''''++=-++= 212212()()()0y y y y y y y y λλλλλλ'''''''=---=---=， 令2z y y λ'=-，则11110x dz z z z z c e dx λλλ'-=? =?=，于是121x y y c e λλ'-=，由一阶微分方程的通解公式得 221212()()()1212[][]dx dx x x x y e c e e dx C e c e dx C λλλλλλ----??=+=+?? (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 提供

常微分方程数值解

第四章常微分方程数值解 [课时安排]6学时 [教学课型]理论课 [教学目的和要求] 了解常微分方程初值问题数值解法的一些基本概念，如单步法和多步法，显式和隐式，方法的阶数，整体截断误差和局部截断误差的区别和关系等；掌握一阶常微分方程初值问题的一些常用的数值计算方法，例如欧拉（Euler）方法、改进的欧拉方法、龙贝－库塔（Runge-Kutta）方法、阿达姆斯（Adams）方法等，要注意各方法的特点及有关的理论分析；掌握构造常微分方程数值解的数值积分的构造方法和泰勒展开的构造方法的基本思想，并能具体应用它们导出一些常用的数值计算公式及评估截断误差；熟练掌握龙格－库塔（Ｒ－Ｋ）方法的基本思想，公式的推导，Ｒ－Ｋ公式中系数的确定，特别是能应用“标准四阶Ｒ－Ｋ公式”解题；掌握数值方法的收敛性和稳定性的概念，并能确定给定方法的绝对稳定性区域。[教学重点与难点] 重点：欧拉方法，改进的欧拉方法，龙贝－库塔方法。 难点：R—K方法，预估-校正公式。 [教学内容与过程] 4.1 引言 本章讨论常微分方程初值问题 (4.1.1) 的数值解法，这也是科学与工程计算经常遇到的问题，由于只有很特殊的方程能用解析方法求解，而用计算机求解常微分方程的初值问题都要采用数值方法.通常我们假定(4.1.1)中 f(x,y)对y满足Lipschitz条件，即存在常数L＞0，使对，有 (4.1.2) 则初值问题(4.1.1)的解存在唯一. 假定(4.1.1)的精确解为，求它的数值解就是要在区间上的一组离散点 上求的近似.通常取 ，h称为步长，求(4.1.1)的数值解是按节点的顺序逐步 推进求得.首先，要对方程做离散逼近，求出数值解的公式，再研究公式的局部截

常微分方程解题方法总结.doc

常微分方程解题方法总结 来源：文都教育 复习过半， 课本上的知识点相信大部分考生已经学习过一遍 . 接下来， 如何将零散的知 识点有机地结合起来， 而不容易遗忘是大多数考生面临的问题 . 为了加强记忆， 使知识自成 体系，建议将知识点进行分类系统总结 . 著名数学家华罗庚的读书方法值得借鉴， 他强调读 书要“由薄到厚、由厚到薄”，对同学们的复习尤为重要 . 以常微分方程为例， 本部分内容涉及可分离变量、 一阶齐次、 一阶非齐次、 全微分方程、 高阶线性微分方程等内容， 在看完这部分内容会发现要掌握的解题方法太多， 遇到具体的题 目不知该如何下手， 这种情况往往是因为没有很好地总结和归纳解题方法 . 下面以表格的形 式将常微分方程中的解题方法加以总结，一目了然，便于记忆和查询 . 常微分方程 通解公式或解法 ( 名称、形式 ) 当 g( y) 0 时，得到 dy f (x)dx ， g( y) 可分离变量的方程 dy f ( x) g( y) 两边积分即可得到结果； dx 当 g( 0 ) 0 时，则 y( x) 0 也是方程的 解 . 解法：令 u y xdu udx ，代入 ，则 dy 齐次微分方程 dy g( y ) x dx x u g (u) 化为可分离变量方程 得到 x du dx 一 阶 线 性 微 分 方 程 P ( x)dx P ( x) dx dy Q(x) y ( e Q( x)dx C )e P( x) y dx

伯努利方程 解法：令 u y1 n，有 du (1 n) y n dy ， dy P( x) y Q( x) y n（n≠0,1）代入得到du (1 n) P(x)u (1 n)Q(x) dx dx 求解特征方程：2 pq 三种情况： 二阶常系数齐次线性微分方程 y p x y q x y0 二阶常系数非齐次线性微分方程 y p x y q x y f ( x) (1)两个不等实根：1, 2 通解： y c1 e 1x c2 e 2x (2) 两个相等实根：1 2 通解： y c1 c2 x e x (3) 一对共轭复根：i , 通解： y e x c1 cos x c2 sin x 通解为 y p x y q x y 0 的通解与 y p x y q x y f ( x) 的特解之和. 常见的 f (x) 有两种情况： x （ 1）f ( x)e P m ( x) 若不是特征方程的根，令特解 y Q m ( x)e x；若是特征方程的单根，令特 解 y xQ m ( x)e x；若是特征方程的重根， 令特解 y*x2Q m (x)e x； （2）f (x) e x[ P m ( x) cos x p n ( x)sin x]

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 提供了多种求解器solver ，对于不同的ODE 问题，采用不同的solver.

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

高阶线性微分方程常用解法简介 摘要：本文主要介绍高阶线性微分方程求解方法，主要的内容有高阶线性微分方程求解的常 用方法如。 关键词：高阶线性微分方程 求解方法 在微分方程的理论中，线性微分方程是非常值得重视的一部分内容，这不仅 因为线性微分方程的一般理论已被研究的十分清楚，而且线性微分方程是研究非线性微分方程的基础，它在物理、力学和工程技术、自然科学中也有着广泛应用。下面对高阶线性微分方程解法做一些简单介绍. 讨论如下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 为任意常数. 如果特征方程有复根，则因方程的系数是实常数，复根将称对共轭的出现.

解微分方程方法

MATLAB解微分方程(2011-07-15 17:35:25) 转载▼ 分类：matlab学习标签： 教育

先说明一下最常用的ode45调用方式，和相应的函数文件定义格式。 [t,x]=ode45(odefun,tspan,x0); 其中，Fun就是导函数，tspan为求解的时间区间（或时间序列，如果采用时间序列，则必须单调），x0为初值。 这时，函数文件可以采用如下方式定义 function dx=odefun(t,x) 对于上面的小例子，可以用如下的程序求解。

2.终值问题 tspan可以是递增序列，也可以为递减序列，若为递减则可求解终值问题。 [t,x]=ode45(@zhongzhiode,[3,0],[1;0;2]);plot(t,x) function dx=zhongzhiode(t,x) dx=[2*x(2)^2-2; -x(1)+2*x(2)*x(3)-1; -2*x(2)+2*x(3)^2-4]; 结果如下 3.odeset options = odeset('name1',value1,'name2',value2,...) [t,x]=solver(@fun,tspan,x0,options) 通过odeset设置options 第一，通过求解选项的设置可以改善求解精度，使得原本可能不收敛的问题收敛。options=odeset('RelTol',1e-10);

第二，求解形如M(t,x)x'=f(t,x)的方程。 例如，方程 x'=-0.2x+yz+0.3xy y'=2xy-5yz-2y^2 x+y+z-2=0 可以变形为 [1 0 0][x'] [-0.2x+yz+0.3xy] [0 1 0][y']=[2xy-5yz-2y^2 ] [0 0 1][z'] [x+y+z-2 ] 这样就可以用如下的代码求解该方程 function mydae M=[1 0 0;0 1 0;0 0 0]; options=odeset('Mass',M); x0=[1.6,0.3,0.1]; [t,x]=ode15s(@daedot,[0,1.5],x0,options);plot(t,x) function dx=daedot(t,x) dx=[ -0.2*x(1)+x(2)*x(3)+0.3*x(1)*x(2); 2*x(1)*x(2)-5*x(2)*x(3)-2*x(2)*x(2); x(1)+x(2)+x(3)-2]; 4.带附加参数的ode45

几类三阶常微分方程的通解公式【开题报告】

毕业论文开题报告 数学与应用数学 几类三阶常微分方程的通解公式 一、选题的背景、意义 常微分方程是指包含一个自变量和它的未知函数以及未知函数的微分的等式。微分方程差不多是和微积分同时产生的，它的形成和发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关。20世纪30年代中期法国数学家勒雷和绍尔建立了LeraySchauder度理论[1]。他们的方法用于研究线性微分、积分、泛函数方程时，取得了巨大成功。 常微分方程在很多学科领域内有着重要的作用，自动控制、各种电子学装置的设计、弹道的计算、飞机和导弹飞行的稳定性的研究、化学反应过程稳定性的研究等等，这些问题都可以归结为高阶微分方程的模型[1，2]，或者化为研究解的性质的问题。很多物理与技术问题都可以化归为微分方程的求解问题。牛顿研究天体力学和机械力学的时候，利用了微分方程这个工具，从理论上得到了行星运动规律。后来，法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星的位置。这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量。 微分方程的理论逐步完善的时候，利用它就可以精确地表述事物变化所遵循的基本规律，只要列出相应的微分方程，就会有解方程的方法[3-5]。微分方程也就成了最有生命力的数学分支。常微分方程是数学分析或基础数学的一个组成部分，在整个数学大厦中占据着重要位置。 有关三阶常微分方程的求解研究已经取得了较为丰富的结果，下面对研究三阶常微分方程的通解详见文献[6-10]。 二、研究的基本内容与拟解决的主要问题 本文主要是对三阶常微分方程通解的研究，具体研究的基本内容与拟解决的主要问题如下： 问题1 如果已知三阶线性微分方程 ()()()() +++= y P x y Q x y R x y f x ''''''

二次微分方程的通解

二次微分方程的通解 Company number：【0089WT-8898YT-W8CCB-BUUT-202108】

第六节 二阶常系数齐次线性微分方程 教学目的：使学生掌握二阶常系数齐次线性微分方程的解法，了解二阶常系数非齐次线 性微分方程的解法 教学重点：二阶常系数齐次线性微分方程的解法 教学过程： 一、二阶常系数齐次线性微分方程 二阶常系数齐次线性微分方程 方程 ypyqy 0 称为二阶常系数齐次线性微分方程 其中p 、q 均为常数 如果y 1、y 2是二阶常系数齐次线性微分方程的两个线性无关解 那么yC 1y 1C 2y 2就是它的通解 我们看看 能否适当选取r 使ye rx 满足二阶常系数齐次线性微分方程 为此将ye rx 代入方程 ypyqy 0 得 (r 2prq )e rx 0 由此可见 只要r 满足代数方程r 2prq 0 函数ye rx 就是微分方程的解 特征方程 方程r 2prq 0叫做微分方程ypyqy 0的特征方程 特征方程的两个根r 1、r 2可用公式 求出 特征方程的根与通解的关系 (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 )(21212 1-==不是常数 因此方程的通解为 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=是方程的解 又

微分方程公式运用表

错误！未找到引用源。 微分方程公式运用表 一、 一阶微分方程 判断特征： (,)dy f x y dx = 类型一：()()dy g x h y dx =（可分离变量的方程） 解法（分离变量法）： ()()dy g x dx h y =，然后两边同时积分。 类型二：()()dy P x y Q x dx +=（一阶线性方程） 解法（常数变易法）：()()(())P x dx P x dx y e C Q x e dx -??=+? 类型三： (,)(,)dy f x y f tx ty dx ==（一阶齐次性方程） 解法（换元法）：y u x =?令类型一 类型四：P()y=Q(x)y n dy x dx +（伯努利方程） 解法（同除法）：1()()n n dy y P x y Q x dx --+=?类型二 二、 可降阶的高阶微分方程 类型一：()()n y f x = 解法（多次积分法）：(1)()()n du u y f x f x dx -=? =?令多次积分求 类型二：''(,')y f x y = 解法：'(,)dp p y f x p dx =?=?令一阶微分方程 类型三：''(,')y f y y = 解法：'(,)dp dp dy dp p y p f y p dx dy dx dy =?==??令类型二 三、线性微分方程 类型一：''()'()0y P x y Q x y ++=（二阶线性齐次微分方程） 解法：找出方程的两个任意线性不相关特解：12(),()y x y x

则：1122()()()y x c y x c y x =+ 类型二：''()'()()y P x y Q x y f x ++=（二阶线性非齐次微分方程） 解法：先找出对应的齐次微分方程的通解：31122()()()y x c y x c y x =+ 再找出非齐次方程的任意特解()p y x ，则：1122()()()()p y x y x c y x c y x =++ 类型三：'''0y py q ++=（二阶线性常系数齐次微分方程） 解法（特征方程法）：22 1,2402p p q p q λλλ-±-++=?= （一）122121240x x p q y c e c e λλλλ?=->?≠?=+ （二）12120()x y c c x e λλλλ?=?==?=+ （三）12120,(cos sin )x i i y e c x c x αλαβλαβββ?

matlab常微分方程和常微分方程组的求解

常微分方程和常微分方程组的求解 一、实验目的： 熟悉Matlab 软件中关于求解常微分方程和常微分方程组的各种命令，掌握利用Matlab 软件进行常微分方程和常微分方程组的求解。 二、相关知识 在MATLAB 中，由函数dsolve()解决常微分方程（组）的求解问题，其具体格式如下： X=dsolve(‘eqn1’,’eqn2’,…) 函数dsolve 用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解。 例1：求解常微分方程1dy dx x y = +的MATLAB 程序为：dsolve('Dy=1/(x+y)','x')， 注意，系统缺省的自变量为t ，因此这里要把自变量写明。 结果为：-lambertw(-C1*exp(-x-1))-x-1 其中：Y=lambertw(X)表示函数关系Y*exp(Y)=X 。 例2：求解常微分方程 2 '''0yy y -=的MATLAB 程序为：Y2=dsolve('y*D2y-Dy^2=0’,’x’) 结果为： Y2 =[ exp((x+C2)/C1)] [ C2] 我们看到有两个解，其中一个是常数。 例3：求常微分方程组253t t dx x y e dt dy x y e dt ?++=??? ?--=??通解的MATLAB 程序为： [X,Y]=dsolve('Dx+5*x+y=exp(t),Dy-x-3*y=exp(2*t)','t') 例4：求常微分方程组020 210cos ,224,0 t t t dx dy x t x dt dt dx dy y e y dt dt =-=?+-==??? ?++==??通解的MATLAB 程序 为：

求下列微分方程的通解

第一章 绪 论 例1-1 求下列微分方程2 3x dx dy =的通解，并分别求满足下列条件的特解。 （1）通过点)1,2(； （2）与直线x y =相切； （3）与直线13+-=x y 正交。 解 直接积分得方程的通解为C x y +=3。 （1）将代入通解中1,2==y x 得7-=C ，则通过点)1,2(解为73-=x y 。 （2）与直线x y =相切的解满足在切点处斜率相同，有132=x ，即得3 1± =x ，切 点坐标为)3 1, 3 1( 和)3 1,31(- - 。同（1）的解法，与直线x y =相切的解为 3 323 + =x y 和3 32 3- =x y 。 （3）与直线13+-=x y 正交的解在正交点处斜率满足3 132 =x ，即得3 1± =x ，正交 点坐标为)0,31 (和)2,3 1(- 。同（1）的解法所求方程的解为27 553 +=x y 和27 13 -=x y 。 评注：求方程满足某条件的特解，关键要找到所求积分曲线经过的某一特定点的坐标，代入通解中确定出任意常数即可得特解。 例1-2 求与曲线族x Ce y =正交的曲线族。 解 因为曲线族x Ce y =满足的微分方程为y y ='，所以与曲线族x Ce y =正交的曲线族满足的微分方程为y y 1- ='，解之得C x y +-=22 ，这就是所求曲线族方程。 评注：首先对已给定的曲线族求得其满足的微分方程，其次借助于正交性得到所求曲线族满足的微分方程，再求解此微分方程。有时直接给出一个微分方程，要求求得与此微分方程的积分曲线族正交（或夹角为某一固定值）的曲线族。 例1-3 求一曲线方程，使曲线上任一点平分过该点的法线在两坐标轴之间的线段。 解 设所求的曲线为)(x y y =，过曲线上任一点),(y x 的法线方程为