ode45解微分方程

ode45求解微分方程

一、ode45求解一阶常微分方程

求解一阶常微分方程：23*t

y

y t *=+ ，求解区间[1 4]，初值y0=-2。 MATLAB 程序：

%%%%%% ode45求解一阶常微分方程

%%%%%% t^2*Dy=y+3*t

%%%%%% 求解区间[1 4], 初值y0=-2

odefun=@(t,y)(y+3*t)/t^2; %定义函数

tspan=[1 4]; %求解区间

y0=-2; %初值

[t,y]=ode45(odefun,tspan,y0);

plot(t,y) %作图

title('t^2y''=y+3t,y(1)=-2,1

legend('t^2y''=y+3t')

xlabel('t')

ylabel('y')

% 精确解

% dsolve('t^2*Dy=y+3*t','y(1)=-2')

% ans =一阶求解结果图

% (3*Ei(1)-2*exp(1))/exp(1/t)-(3*Ei(1/t))/exp(1/t) Ei(x)=exp(t)/t 从负无穷到x 的积分

图1

二、ode45求解二阶常微分方程

求解二阶常微分方程：

sin

y y y t

++=

，初始条件为y(0)=5，y'(0)=6，求

解区间为[0 20]。

MATLAB程序：

%%%%%%% y''+y'+y=sin(t)

%%%%%%% 初始条件为y(0)=5,y'(0)=6.

%%%%%%% 先降阶为一阶微分方程组

%%%%%%% y'=z

%%%%%%% z'=-z-y+sin(t)

%%%%%%% 然后用ode45解方程

%%%%%%% [t,y]=ode45(@weifen,[0 20],[5 6])

%%%%%%% plot(t,y)就可以画出y和y'的图像

[t,y]=ode45(@weifen,[0 20],[5 6]);

%%%%%[0 20]为t的区间，[5 6]为y(0)和y'(0)的初值

plot(t,y)

function dy=weifen(t,x)

dy=zeros(2,1);

%y=x(1)

%z=x(2)

dy(1)=x(2);

dy(2)=sin(t)-x(2)-x(1);

end

图2

注：求解二阶微分方程只需把二阶微分方程转化为一阶微分方程，再进行求解即可。

二阶微分方程解法知识讲解

二阶微分方程解法

第六节 二阶常系数齐次线性微分方程 教学目的：使学生掌握二阶常系数齐次线性微分方程的解法，了解二阶常系数非齐 次线性微分方程的解法 教学重点：二阶常系数齐次线性微分方程的解法 教学过程： 一、二阶常系数齐次线性微分方程 二阶常系数齐次线性微分方程: 方程 y ''+py '+qy =0 称为二阶常系数齐次线性微分方程, 其中p 、q 均为常数. 如果y 1、y 2是二阶常系数齐次线性微分方程的两个线性无关解, 那么y =C 1y 1+C 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 )(21212 1-==不是常数. 因此方程的通解为 x r x r e C e C y 2121+=. (2)特征方程有两个相等的实根r 1=r 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, 2=α±i β时, 函数y =e (α+i β)x 、y =e (α-i β)x 是微分方程的两个线性无关的复数形式的解. 函数y =e αx cos βx 、y =e αx sin βx 是微分方程的两个线性无关的实数形式的解. 函数y 1=e (α+i β)x 和y 2=e (α-i β)x 都是方程的解, 而由欧拉公式, 得 y 1=e (α+i β)x =e αx (cos βx +i sin βx ), y 2=e (α-i β)x =e αx (cos βx -i sin βx ), y 1+y 2=2e αx cos βx , )(2 1cos 21y y x e x +=βα, y 1-y 2=2ie αx sin βx , )(21sin 21y y i x e x -=βα. 故e αx cos βx 、y 2=e αx sin βx 也是方程解. 可以验证, y 1=e αx cos βx 、y 2=e αx sin βx 是方程的线性无关解. 因此方程的通解为

二阶线性微分方程的解法

二阶常系数线性微分方程 一、二阶常系数线形微分方程的概念 形如 )(x f qy y p y =+'+'' (1) 的方程称为二阶常系数线性微分方程.其中p 、q 均为实数,)(x f 为已知的连续函数. 如果0)(≡x f ,则方程式 (1)变成 0=+'+''qy y p y (2) 我们把方程(2)叫做二阶常系数齐次线性方程,把方程式(1)叫做二阶常 系数非齐次线性方程. 本节我们将讨论其解法. 二、二阶常系数齐次线性微分方程 1．解的叠加性 定理1 如果函数1y 与2y 是式(2)的两个解, 则2211y C y C y +=也是 式(2)的解,其中21,C C 是任意常数. 证明 因为1y 与2y 是方程(2)的解,所以有 0111 =+'+''qy y p y 0222 =+'+''qy y p y 将2211y C y C y +=代入方程(2)的左边,得 )()()(22112211221 1y C y C q y C y C p y C y C ++'+'+''+'' =0)()(2222111 1=+'+''++'+''qy y p y C qy y p y C 所以2211y C y C y +=是方程(2)的解. 定理1说明齐次线性方程的解具有叠加性. 叠加起来的解从形式看含有21,C C 两个任意常数,但它不一定是方程式(2)的通解. 2.线性相关、线性无关的概念

设,,,,21n y y y 为定义在区间I 内的n 个函数,若存在不全为零的常数 ,,,,21n k k k 使得当在该区间内有02211≡+++n n y k y k y k , 则称这n 个函数在区间I 内线性相关,否则称线性无关. 例如 x x 22sin ,cos ,1在实数范围内是线性相关的,因为 0sin cos 12 2≡--x x 又如2,,1x x 在任何区间(a,b)内是线性无关的,因为在该区间内要使 02321≡++x k x k k 必须0321===k k k . 对两个函数的情形,若=21y y 常数, 则1y ,2y 线性相关,若≠2 1y y 常数, 则1y ,2y 线性无关. 3.二阶常系数齐次微分方程的解法 定理 2 如果1y 与2y 是方程式(2)的两个线性无关的特解,则 212211,(C C y C y C y +=为任意常数)是方程式(2)的通解. 例如, 0=+''y y 是二阶齐次线性方程,x y x y cos ,sin 21==是它的 两个解,且≠=x y y tan 2 1常数,即1y ,2y 线性无关, 所以 x C x C y C y C y cos sin 212211+=+= ( 21,C C 是任意常数)是方程0=+''y y 的通解. 由于指数函数rx e y =(r 为常数)和它的各阶导数都只差一个常数因子, 根据指数函数的这个特点,我们用rx e y =来试着看能否选取适当的常数r , 使rx e y =满足方程(2).

二阶常微分方程解

第七节 二阶常系数线性微分方程 在上节我们已经讨论了二阶线性微分方程解的结构，二阶线性微分方程的求解问题，关键在于如何求二阶齐次方程的通解和非齐次方程的一个特解。本节讨论二阶线性方程的一个特殊类型，即二阶常系数线 性微分方程及其求解方法。先讨论二阶常系数线性齐 §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)

二阶常微分方程解

二阶常微分方程解

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

第七节 二阶常系数线性微分方程 的解法 在上节我们已经讨论了二阶线性微分方程解的结构，二阶线性微分方程的求解问题,关键在于如何求二阶齐次方程的通解和非齐次方程的一个特解。本节讨论二阶线性方程的一个特殊类型，即二阶常系数线性微分方程及其求解方法。先讨论二阶常系数线性齐次方程的求解方法。 §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)的两个特解。

二阶线性微分方程解的结构

二阶线性微分方程解的结构

————————————————————————————————作者: ————————————————————————————————日期： ?

附录A 线性常微分方程 本课程的研究内容与常微分方程理论有非常密切的联系,因此在本附录里,我们将对线性常微分方程的知识——包括解的存在性、解的结构和求解方法做一些回顾和总结。 把包含未知函数和它的j 阶导数()j y (的方程称为常微分方程。线性常微分方程的标准形式 ()(1)110()()'()()n n n y p x y p x y p x y f x --++++= (A.1) 其中n 称为方程的阶数,()j p x 和()f x 是给定的函数。可微函数()y y x =在区间 I 上满足方程（Ａ．１），则称其为常微分方程(Ａ．1）在 I 上的一个解。，()f x 称为方程（Ａ.1)的自由项,当自由项()0f x ≡时方程（A.1）称为是齐次方程,否则称为非齐次方程。一般来说常微分方程的解是不唯一的，我们将方程的全部解构成的集合称为解集合,解集合中全部元素的一个通项表达式称为方程的通解，而某个给定的解称为方程的特解。 在本附录里,我们重点介绍一阶和二阶常微分方程的相关知识。 A.１ 一阶线性常微分方程 一阶线性常微分方程表示为 '()()y p x y f x x I +=∈，. (A.2) 当()0f x ≡,方程退化为 '()0y p x y +=, （A.3) 假设()y x 不恒等于零,则上式等价于

'()y p x y =- 而()'ln 'y y y =，从而（A.3)的通解为 ()d ()p x x y x Ce -?= ( A.4） 对于非齐次一阶线性常微分方程(A ．２）,在其两端同乘以函数()d p x x e ? ()d ()d ()d '()()p x x p x x p x x e y p x e y e f x ???+= 注意到上面等式的左端 ()d ()d ()d ''()p x x p x x p x x e y p x e y e y ?????+= ??? ‘ 因此有 ()d ()d '()p x x p x x e y e f x ????= ??? ‘ 两端积分 ()d ()d ()d p x x p x x e y C e f x x ??=+?‘ 其中C 是任意常数。进一步有 ()d ()d ()d p x x p x x y e C e f x x -????=+ ??? ?‘ 综上有如下结论 定理A.１ 假设()()p x f x I 和在上连续，则一阶线性非齐次常微分方程（A.1)的通解具有如下形式 ()d ()d ()d ()()d p x x p x x p x x y x Ce e e f x x --???=+?‘ （A.5)

二阶线性微分方程解的结构

附录A 线性常微分方程 本课程的研究内容与常微分方程理论有非常密切的联系，因此在本附录里，我们将对线性常微分方程的知识——包括解的存在性、解的结构和求解方法做一些回顾和总结。 把包含未知函数和它的j 阶导数()j y (的方程称为常微分方程。线性常微分方程的标准形式 ()(1)110()()'()()n n n y p x y p x y p x y f x --++++= （A.1） 其中n 称为方程的阶数，()j p x 和()f x 是给定的函数。可微函数()y y x =在区间 I 上满足方程（A.1），则称其为常微分方程（A.1）在 I 上的一个解。，()f x 称为方程（A.1）的自由项，当自由项()0f x ≡时方程（A.1）称为是齐次方程，否则称为非齐次方程。一般来说常微分方程的解是不唯一的，我们将方程的全部解构成的集合称为解集合，解集合中全部元素的一个通项表达式称为方程的通解，而某个给定的解称为方程的特解。 在本附录里，我们重点介绍一阶和二阶常微分方程的相关知识。 A.1 一阶线性常微分方程 一阶线性常微分方程表示为 '()()y p x y f x x I +=∈，. （A.2） 当()0f x ≡，方程退化为 '()0y p x y +=, （A.3） 假设()y x 不恒等于零，则上式等价于 '()y p x y =-

而()'ln 'y y y =，从而（A.3）的通解为 ()d ()p x x y x Ce -?= （ A.4） 对于非齐次一阶线性常微分方程（A.2），在其两端同乘以函数()d p x x e ? ()d ()d ()d '()()p x x p x x p x x e y p x e y e f x ???+= 注意到上面等式的左端 ()d ()d ()d ''()p x x p x x p x x e y p x e y e y ?????+= ??? ‘ 因此有 ()d ()d '()p x x p x x e y e f x ????= ??? ‘ 两端积分 ()d ()d ()d p x x p x x e y C e f x x ??=+?‘ 其中C 是任意常数。进一步有 ()d ()d ()d p x x p x x y e C e f x x -??? ?=+ ??? ?‘ 综上有如下结论 定理A.1 假设()()p x f x I 和在上连续，则一阶线性非齐次常微分方程（A.1）的通解具有如下形式 ()d ()d ()d ()()d p x x p x x p x x y x Ce e e f x x --???=+?‘ （A.5） 其中C 是任意常数。 观察（A.4）式和（A.5）式，我们发现一阶线性非齐次常微分方程（A.1） 的解等于一阶线性齐次常微分方程（A.2）的通解()d p x x Ce -?加上函数

二次微分方程的通解

第六节 二阶常系数齐次线性微分方程 教学目的：使学生掌握二阶常系数齐次线性微分方程的解法，了解二阶常系数非齐 次线性微分方程的解法 教学重点：二阶常系数齐次线性微分方程的解法 教学过程： 一、二阶常系数齐次线性微分方程 二阶常系数齐次线性微分方程: 方程 y ''+py '+qy =0 称为二阶常系数齐次线性微分方程, 其中p 、q 均为常数. 如果y 1、y 2是二阶常系数齐次线性微分方程的两个线性无关解, 那么y =C 1y 1+C 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 1=r 2时, 函数x r e y 11=、x r xe y 12=是二阶常系数齐次线性微分

二阶常微分方程的降阶解法

郑州航空工业管理学院 毕业论文（设计） 2015届数学与应用数学专业1111062班级 题目二阶常微分方程的降阶解法 姓名贾静静学号111106213 指导教师程春蕊职称讲师 2015年4月5号

二阶常微分方程的降阶解法 摘要 常微分方程是数学领域的一个非常重要的课题，并在实践中广泛于解决问题，分析模型。常微分方程在微分理论中占据首要位置,普遍应用在工程应用、科学研究以及物理学方面,不少应用范例都归结为二阶线性常微分方程的求解问题。而正常情况下，常系数微分方程依据线性常微分方程的日常理论是可以求解的.不过对于变系数二阶线性常微分方程的求解却有一定程度的困难,迄今为止还没有一个行之有效的普遍方法。 本文主要考虑了二阶常系数线性微分方程的降阶法。关于二阶常系数线性微分方程的求解问题,首先,我们给出二阶齐次常系数线性微分方程的特征方程，并求解出特征方程的两个特征根；其次，利用积分因子乘以微分方程和导数的运算，将二阶常系数线性微分方程化为一阶微分形式；最后，将一阶微分形式两边同时积分，求解一阶线性微分方程，可求得二阶常系数线性微分方程的一个特解或通解。关于二阶变系数齐次线性微分方程的求解问题，化为恰当方程通过降阶法求解二阶齐次变系数微分方程的通解。对于非齐次线性微分方程，只需再运用常数变易法求出它的一个特解，问题也就相应地解决了。 关键词 二阶常微分方程；降阶法；特征根；常数变易法；一阶微分形式

Order reduction method of second order ordinary differential equations Jingjing Jia Chunrui Cheng 111106213 Abstract Ordinary differential equation is a very important topic in the field of mathematics, it has been widely used in solving the problem and analyzing model in practice . Ordinary differential equations in the theory of differential occupied first place, it has been widely used in engineering application and scientific research as well as physics, many application examples are attributed to second order linear ordinary differential equation solving problem. And under normal circumstances,ordinary coefficient differential equation on the basis of the linear often daily theory of differential equations is can be solved. But for the solution for variable coefficient second order linear ordinary differential equations have a certain degree of difficulty, so far we haven't a well-established general method. This paper mainly introduces the method of reduction of order two order linear differential equation with constant coefficients.On the problem of solving the linear differential equation with two order constant coefficients,first, we give homogeneous ordinary coefficient linear differential equation of the characteristic equation and solve the two characteristic roots of characteristic equation;secondly,we should use the integral factor times differential equation and derivative operation and turn two order constant

二阶微分方程解法

第六节二阶常系数齐次线性微分方程 教学目的：使学生掌握二阶常系数齐次线性微分方程的解法，了解二阶常系数非齐次线性微分方程的解法 教学重点：二阶常系数齐次线性微分方程的解法 教学过程： 一、二阶常系数齐次线性微分方程 二阶常系数齐次线性微分方程：方程 y'H；Py qy -0 称为二阶常系数齐次线性微分方程.其中p、q均为常数， 如果y i、y是二阶常系数齐次线性微分方程的两个线性无关解.那么y=C1y1?C2y2就是它的通解 我们看看.能否适当选取r使yw rx满足二阶常系数齐次线性微分方程.为此将y=e rx代入方程 y py qy=0 得 (r2Pr q)e rx=0 . 由此可见.只要r满足代数方程r2Pr q -0 .函数y=e rx就是微分方程的解， 特征方程：方程r2pr q -0叫做微分方程y…卩y'qy=0的特征方程，特征方程的两个根r i、r2可用公式 _p +±J p2 _4q r 1,2 = 2 求出 特征方程的根与通解的关系： (1) 特征方程有两个不相等的实根 的解 这是因为. 函数y i =e l^ix、 y e「i X y2=e r2x是方程的解?又-e rx -e(ri -r2)x不是常数，y2 e2 因此方程的通解为 y =C i e rix C2e r2x. (2) 特征方程有两个相等的实根r i=Γ2时?函数y i =e rix、y2 =xe rix是二阶常系数齐次线性微分 r i、—时.函数y1 =e"x、y2 =e r2χ是方程的两个线性无关

方程的两个线性无关的解 这是因为.y 1=e r ιx 是方程的解.又 (Xe rI X ) - P(Xe rlX ) ? q(xe l ^ιx ) = (2r 1 ? xr 12)e l ^ιx ? p(1 ■ xr 1)e rlx ■ qxe r ι = e l ^ιx (2r 1 亠 P)亠Xe rI X (r 12 亠 pr 1 亠q) =O y 护 所以y 2 =xe r M 也是方程的解 且竺=空一 X 不是常数. y 1 e r1x 因此方程的通解为 y =C I e r1X C 2xe r1x ⑶特征方程有一对共轭复根 r 1,2二〉L 时.函数以y=e c J)X 是微分方程的两个线性无 关的复数形式的解，函数y=e ：X CoS x > y w ：X Sin X 是微分方程的两个线性无关的实数形式的解 函数y^e cJ )X 和y 2=e (H)X 都是方程的解.而由欧拉公式.得 y 1 =e (Cr fr *=eE(cos βx4isinPx). y 2=e (C ?i ^*=e c x (cosβx4sinPx). y 1 4y 2=2e c x cosPx . e" cos Px=1 (y 1 +y 2). 2 y 1-y 2=2ie c x si nβ×.e c x sin βx=?y 1-y 2), 2i 故e^cosPx 、y 2=e 0?nβ×也是方程解， 可以验证$1 =e"cos ：x 、y2=e "sin ：X 是方程的线性无关解， 因此方程的通解为 y=e c x (C 1cosβx 4C 2sinPx ), 求二阶常系数齐次线性微分方程 / p/ q^0的通解的步骤为： 第一步 写出微分方程的特征方程 2 r Pr q=0 第二步求出特征方程的两个根 「1、匕， 第三步 根据特征方程的两个根的不同情况 .写出微分方程的通解. 例1求微分方程y …~2,?3y=0的通解 解所给微分方程的特征方程为 2 r -2r-3=0 .即(r 1)(r-3^0 . 其根「1—1「2 =3是两个不相等的实根.因此所求通解为 yQef C 2e 3x . 例2求方程y=2y?yp 满足初始条件y ∣x=o=4、y]χ=o=…2的特解

常系数二阶微分方程的齐次通解

常系数二阶微分方程的齐次通解

————————————————————————————————作者：————————————————————————————————日期：

附录2 常系数二阶微分方程的齐次通解 常系数二阶齐次微分方程 0=+2+2022y dt dy dt y d ωα 设其中α、ω0都是正实数。 要使二阶微分方程有确定的解，必须知道两个初始条件：初始值y (0)和一阶导数的初始值0 =t dt dy 。 这里只讨论齐次通解在一些典型的系数值下的特点，不求出解中的待定常数。目的在于避免过多的数学式子，突出对有普遍意义的特征的认识。 尝试St e y =（S 为实的或复的常数）是否能为方程的解。 代入方程可得恒等式： 0=)+2+(202S S S e St ωα 由此得到决定常数S 的特征方程： 0=+2+202ωαS S 该一元二次代数方程的根为： 202-±-=ωααS 因常数项的值不同，解的形式不同： 1.自由振荡情况(无阻尼情况)（0=α） 此时，S 是一对共轭虚数： 01j =ωS 02-j =ωS 齐次通解为： t t e K e K t y 00-j 2j 1+=)(ωω 变为常用的三角函数式 )+sin(=)(0θωt K t y 这是一个等幅正弦振荡，ω0 是自由振荡角频率或谐振角频率。K 和θ 是由初始条件决定的常数。 2.欠阻尼情况（ 0<<0ωα ） 此时，S 是一对共轭复数： d 1j +-=ωαS d 2j --=ωαS 齐次通解为： )+sin(=)(d -θωαt Ke t y t 这是一个衰减振荡。其中，220-=αωωd （正实数）是衰减振荡角频率。 振幅按指数函数t e α-衰减，故称α为衰减系数。 K 和θ 是由初始条件决定的常数。 这种情况下，系统开始会有正弦振荡，但随时间而衰减，过一段时间后就消失。 3.过阻尼情况（0>ωα）

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

目录 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年来数学分析是数学的首要分支，而微分方程

(完整版)专题一(二阶常微分方程解法)

二阶微分方程： 时为非齐次 时为齐次，0)(0)()()()(22≠≡=++x f x f x f y x Q dx dy x P dx y d 二阶常系数齐次线性微分方程及其解法： 2 122,)(2,,(*)0)(1,0(*)r r y y y r r q pr r q p qy y p y 式的两个根、求出的系数； 式中的系数及常数项恰好是，，其中、写出特征方程：求解步骤： 为常数； ，其中?'''=++?=+'+''式的通解：出的不同情况，按下表写、根据(*),321r r 二阶常系数非齐次线性微分方程 型 为常数； 型，为常数 ，]sin )(cos )([)()()(,)(x x P x x P e x f x P e x f q p x f qy y p y n l x m x ωωλλλ+===+'+'' 二阶常系数非齐次线性微分方程的一般形式是 ''+'+=y py qy f x () （1） 其中p q ,是常数。 方程（1）的通解为对应的齐次方程 0=+'+''qy y p y （2） 的通解Y 和方程（1）的一个特解*y 之和。即 *y Y y +=.我们已解决了求二阶常系数齐 次线性方程通解的问题，所以，我们只需讨论求二阶常系数非齐次线性微分方程的特解* y 的方法。 下面我们只介绍当方程（1）中的)(x f 为如下两种常见形式时求其特解*y 的方法。 一、 f x e P x x m ()()=?λ型 由于方程（1）右端函数f x ()是指数函数e x λ?与m 次多项式P x m ()的乘积，而指数

函数与多项式的乘积的导数仍是这类函数，因此，我们推测： 方程（1）的特解应为 y e Q x x *?=λ()( Q x ()是某个次数待定的多项式 ) y e Q x e Q x x x *??'=+'λλλ()() y e Q x Q x Q x x *?"=?+'+''λλλ[()()()]22 代入方程（1），得 e Q x p Q x p q Q x e P x x x m λλλλλ???''++'+++≡?[()()()()()]()22 消去e x λ?，得 ''++'+++≡Q x p Q x p q Q x P x m ()()()()()()22λλλ （3） 讨论 01、如果λ不是特征方程 r pr q 20++=的根。 即 02≠++q p λλ 由于P x m ()是一个m 次的多项式，欲使（3）的两端恒等，那未Q x ()必为一个m 次多项式，设为 Q x b x b x b x b m m m m m ()=++++--0111Λ 将之代入（3），比较恒等式两端x 的同次幂的系数，就得到以b b b b m m 01 1,,,,Λ-为未知数的m +1个线性方程的联立方程组，解此方程组可得到这m +1个待定的系数，并得到特解 y e Q x x m *?=λ() 02、如果λ是特征方程 r pr q 20++=的单根。 即 λλ20++=p q ，但 20λ+≠p 欲使（3）式的两端恒等，那么'Q x ()必是一个m 次多项式。 因此，可令 Q x x Q x m ()()=? 并且用同样的方法来确定)(x Q 的系数b b b b m m 0 11,,,,Λ-。 03、如果λ是特征方程 r pr q 20++=的二重根。 即 λλ20++=p q ，且 20λ+=p 。 欲使（3）式的两端恒等，那么''Q x ()必是一个m 次多项式 因此， 可令 Q x x Q x m ()()=?2 并且用同样的方法来确定)(x Q 的系数b b b b m m 011,,,,Λ-。

阶常微分方程解

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

y ＝e rx (其中r 为待定常数)来试解 将y ＝e rx ，dx dy ＝re rx ，22 dx y d ＝r 2e rx 代入方程 得 r 2e rx ＋pre rx ＋qe rx ＝0 或 e rx （r 2 ＋pr ＋q ）＝0 因为e rx ≠0，故得 r 2＋pr ＋q ＝0 由此可见，若r 是二次方程 r 2＋pr ＋q ＝0 的根，那么e rx 就是方程的特解，于是方程的求解问题，就转化为求代数方程的根问题。称式为微分方程的特征方程。 特征方程是一个以r 为未知函数的一元二次代数方程。特征方程的两个根r １，r 2，称为特征根，由代数知识，特征根r 1，r 2有三种可能的情况，下面我们分别进行讨论。 (1)若特证方程有两个不相等的实根r １，r 2，此时e r １x ， e r2x 是方程的两个特解。 因为 x r x r 21e e ＝e x )r r (21 ≠常数 所以e r1x ，e r2x 为线性无关函数，由解的结构定理知，方程的通解为 y ＝C 1e r1x ＋C 2e r2x

二阶线性微分方程及其解法

n 阶微分方程的一般形式为： () (,,',",,)0n F x y y y y =L ， 一般情况下，求n 阶微分方程的解是困难的. 作为基础知识，本节仅讨论二阶常系数线性微分方程的求解方法. 一、 二阶线性微分方程解的结构 如果二阶微分方程)',,(''y y x F y =的未知函数及其导数都是一次项的，称为二阶线性微分方程. 二阶线性微分方程的一般形式为 ).()(')(''x f y x q y x p y =++ （） 如果0)(≡x f ，则方程（）成为 .0)(')(''=++y x q y x p y （） 方程（）称为二阶齐次线性微分方程，相应地，方程（）称为二阶非齐次线性微分方程. 定理 齐次线性微分方程解的叠加性定理. 设1y 和2y 是二阶齐次线性微分方程（）的两个解，则 2211y c y c y += 也是微分方程（）的解，其中21,c c 为任意常数. 证： 将2211y c y c y +=代入方程（）的左端，可得 ))(()')((')'(221122112211y c y c x q y c y c x p y c y c +++++ ))(()'')(()''''(221122112211y c y c x q y c y c x p y c y c +++++= =+++))(')(''(1111y x q y x p y c ))(')(''(2222y x q y x p y c ++ =0， 所以，2211y c y c y +=也是微分方程（）的解.□ 定理表明，二阶齐次线性微分方程的解可叠加. 如果我们已知二阶齐次线性微分方程的两个解1y 和2y ，很容易得到含有任意常数21,c c 的解，2211y c y c y +=. 如果解1y 和2y 有一定关系，那么，解2211y c y c y +=中的任意常数21,c c 可以合并成一个任意常数. 因此，依据本章第一节的论述，它并不是二阶齐次线性微分方程的通解. 那么，二阶齐次线性微分方程的两个解1y 和2y 要满足哪些条件才能使解2211y c y c y +=成为二阶齐次线性微分方程的

二阶非齐次线性微分方程的解法

目 录 待定系数法 常数变异法 幂级数法 特征根法 升阶法 降阶法 关键词：微分方程，特解，通解， 二阶齐次线性微分方程 常系数微分方程 待定系数法 解决常系数齐次线性微分方程[]21220, (1) d x dx L x a a x dt dt ≡++= 12,. a a 这里是常数 特征方程212()0F a a λλλ=++= (1.1) (1)特征根是单根的情形 设 12,,,n λλλ是特征方程的 (1.1)的2个彼此不相等的根，则相应的方程 (1)有如 下2个解： 12,t t e e λλ (1.2) 如果(1,2)i i λ=均为实数，则 (1.2)是方程 (1)的2个线性无关的实值解，而方程 (1)的通解可表示为 1212t t x c e c e λλ=+ 如果方程有复根，则因方程的系数是实系数，复根将成对共轭出现。设 i λαβ=+是一特征根，则i λαβ=-也是特征根，因而与这对共轭复根对应，方程 (1)有两个复值解 (i)t (cos t sin ),t e e i t αβαββ+=+

(i)t (cos t sin ).t e e i t αβαββ-=- 它们的实部和虚部也是方程的解。这样一来，对应于特征方程的一对共轭复根 i λαβ=±，我们可求得方程 (1)的两个实值解 cos ,sin .t t e t e t ααββ （2）特征根有重跟的情形 若10λ=特征方程的 k 重零根，对应于方程 (1)的k 个线性无关的解21 1,t,t ,k t -。 若这个 k 重零根10, λ≠设特征根为12,,,,m λλλ其重数为 1212,,,k (k 2)m m k k k k ++ =。方程 (1)的解为 11112222111,t ,t ;,t , t ; ;,t , t ;m m m m t t k t t t k t t t k t e e e e e e e e e λλλλλλλλλ--- 对于特征方程有复重根的情况，譬如假设i λαβ=+是k 重特征根，则i λαβ=-也是k 重特征根，可以得到方程 (1)的2k 个实值解 2121cos ,cos ,cos ,,cos ,sin ,sin ,sin , ,sin .t t t k t t t t k t e t te t t e t t e t e t te t t e t t e t ααααααααββββββββ-- 例1 求方程 220d x x dt -=的通解。 解 特征方程 210λ-=的根为121,1λλ==-有两个实根，均是单根，故方程的通 解为 12,t t x c e c e -=+ 这里12,c c 是任意常数。 例2 求解方程 220d x x dt +=的通解。 解 特征方程 210λ+=的根为12,i i λλ==-有两个复根， 均是单根，故方程的通解 为 12sin cos ,x c t c t =+