常系数线性微分方程的解的结构分析

常系数线性微分方程的解的结构分析
常系数线性微分方程的解的结构分析

常系数线性微分方程的解的结构分析

【 摘要】在参考和总结了许多场系数线性微分方程的解法的基础上,本文总结了一些常系数微分方程的解的解法,并针对一类常系数线性微分方程的已有结论给予证明,以解给予一些结论证明思路,以及一些实例,并向高阶推广。 【关键词 】常系数 线性 微分方程 结构

一阶常系数齐次线性微分方程

0=+ax dt dx

, (1.1)

的求解 上式可以改写为

adt x

dx

-= , (1.2) 于是变量x 和t 被分离,再将两边积分得

c at x +-=ln , (1.3) 这里的c 为常数。又由对数的定义,上式可以变为 at

ce

x -= , (1.4)

其中c= , 因为x=0也是方程的解,因此c 可以是任意常数。

这里首先是将变量分离,然后再两边积分,从而求出方程的解。这便要方程式可以分离变量的,也就是变量分离方程。 一阶常系数微分方程

)()(x Q y x P dx

dy

+= , (2.1) 其中P (x ),Q(x)在考虑的区间上式连续函数,若Q (x )=0 ,上式就变为

y

x P dx dy

)(= , (2.2)

上式为一阶齐次线性微分方程。还是变量分离方程我们可以参考上面变量分离方程的解法,先进行变量分离得到 dx x P y

dy

)(= , (2.3) 两边同时积分,得到

?

=dx

x p ce y )( , (2.4)

这里c 是常数。

若Q (x )≠ 0 , 那么上式就变成了 一阶非齐次线性微分方程。 我们知道一阶齐次线性微分方程是一阶常微分方程的一种特殊情况,那么可以设想将一阶

齐次线性微分方程的解

?

=dx

x p ce y )( , (2.5)

中的常数c 变易成为待定的函数c (x ),令 ?=dx

x p e x c y )()( , (2.6)

微分之,就可以得到

?+?=dx x p dx

x p e x P x c e dx

x dc dx dy )()()()()( , (2.7)

以(2.7),(2.6)代入2.1,得到 )()()()()()()()()(x Q e x c x p e x P x c e dx

x dc dx x p dx x p dx

x p +?=?+?,(2.8) 即 ?=-dx x p e x Q dx

x dc )()()

(, 积分后得到

c (x )=c dx e x Q dx x p +??

-)()( , (2.9)

这里c 是任意常数,将上式代入(2.6)得到方程(2.1)的通解 ))(()()(c dx e x Q e y dx

x p dx

x p +?

?

=?-

(2.91)

在上面的一阶线性微分方程中,是将一阶齐次线性微分方程中的通解中的常数c 变成c(x) ,常数变易法一阶非齐次线性微分方程的解, 感觉这个方法之所以用x 的未知函数u(x)替换任意常数C,是因为C 是任意的,C 与x 形成函数关系,要确定C,需要由初始条件确定,一个x,确定一个C,也就形成一对一或多对多的映射,也就是函数关系,而这里的C 是任意的,也就可以用一个未知的,也就是任意的函数u(x)来代替,进而求得非齐次线性微分方程的解。这种将常数变异为待定函数的方法,我们通常称为常数变易法。常数变易法实质也是一种变量变换的方法,通过变换(2.6可将方程(2.1)化为变量分离方程。

二阶常系数线性微分方程

(1)二阶常系数线性齐次方程

022=++qy dx

dy

p dx y d (3.1)

其中p 、q 是常数,我们知道,要求方程(3.1)的通解,只要求出其任意两个线性无关的特

解y 1,y 2就可以了,下面讨论这样两个特解的求法。

我们先分析方程(3.1)可能具有什么形式的特解,从方程的形式上来看,它的特点是22dx

y

d ,

dx dy ,y 各乘以常数因子后相加等于零,如果能找到一个函数y ,其22dx

y d ,dx dy ,y 之间只相差一个常数因子,这样的函数有可能是方程(3.1)的特解,在初等函数中,指数函数e rx

符合上述要求,于是我们令

y =e

rx

(其中r 为待定常数)来试解

将y =e rx

,dx dy =re rx ,22dx

y d =r 2e rx

代入方程(3

得 r 2e rx +pre rx +qe rx

或 e rx

(r 2

+pr +q )=

因为e rx ≠0,故得

r 2

+pr +q =

由此可见,若r 是二次方程

r 2

+pr +q =0 (3

的根,那么e rx

就是方程(3.1)的特解,于是方程(3.1)的求解问题,就转化为求代数方程(3.2)

的根问题。称(3.2)式为微分方程(3.1)的特征方程。

特征方程(3.2)是一个以r 为未知函数的一元二次代数方程。特征方程的两个根r 1,r 2,称

为特征根,由代数知识,特征根r 1,r 2有三种可能的情况,下面我们分别进行讨论。

(1)若特证方程(3.2)有两个不相等的实根r 1,r 2,此时x

r e 1,x

r e 2是方程(3.1)的两个特解。

因为x r x r e e 2

1 =x

r r e )(21-≠常数

所以e r1x

,e r2x

为线性无关函数,由解的结构定理知,方程(3.1)的通解为

x r x r e c e c y 2121+=

(2)若特征方程(3.2)有两个相等的实根r 1=r 2,此时p 2

-4q =0,即

有2

21p r r -==,这样只能得到方程(3.1)的一个特解x

r e y 11=,因此,我们还要设法找出另一个满足

≠12y y 常数,的特解y 2,故12y y 应是x 的某个函数,设u y y

=1

2 ,其中u =u(x)为待定函数,即

x

r ue uy y

112==

对y 2求一阶,二阶导数得

x r x r x r e u r dx

du

ue r e dx du dx dy 111)(112+=+= x r e dx

u d dx du r u r dx y d 1)2(2212

12

22++= 将它们代入方程(7.1)得

0)()2(11112212

1=+++++x r x r x r que e u r dx

du

p e dx u d dx du r u r

0])()2([1121122=+++++x

r e u q pr r dx

du p r dx u d 因为01≠x

r e ,且因r 1是特征方程的根,故有2

1r +pr 1+q =0,又因2

1p

r -

=故有2r 1+p =0,于是上式成为

02

2=dx u

d 显然满足02

2=dx u

d 的函数很多,我们取其中最简单的一个

u(x)=

则y 2=xe

是方程(3.1)的另一个特解,且y 1,y 2是两个线性无关的函数,所以方程(3.1)

的通解是

x

r x r x r e x c c xe c e c y

111)(2121+=+=

(3)若特征方程(3.2)有一对共轭复根 r 1=α+i β,r 2=α-i β

此时方程(3.1)有两个特解

y 1=e

(α+i β)x

2

=e

(α-i β)x

y =c 1e

(α+i β)x

+c 2e

(α-i β)

其中c 1,c 2为任意常数,但是这种复数形式的解,在应用上不方便。在实际问题中,常常需

要实数形式的通解,为此利用欧拉公式

e ix

=cosx +isinx ,e

-ix

=cosx -

2

1(e ix +e -ix

)=

21(y 1+y 2)=21e αx (e i βx +e -i βx )=e αx

cos β

i 21(y 1-y 2)=i 21e αx (e i βx -e -i βx )=e αx

sin β

由上节定理一知,2

1(y 1+y 2),2

1(y 1-y 2)是方程(3.1)的两个特解,也即e αx cos βx ,e αx

sin βx 是方程(3.1)的两个特解:且它们线性无关,我们已知,方程(3.1)的通解为

y =c 1e αx

cos βx +c 2e αx

sin β

i

21(e ix -e -ix

)=

或 y =e αx

(

1c cos βx +2c sin β

其中c 1,c 2为任意常数,至此我们已找到了实数形式的通解,其中α,β分别是特征方程

(3.2)复数根的实部和虚部。

综上所述,求二阶常系数线性齐次方程(3.1)的通解,只须先求出其特征方程(3.2)的根,再

根据他的三种情况确定其通解,现列表如下

特征方程r 2

+pr +q =0的根

微分方程02

2=++qy dx dy

p dx y d 的通解

有二个不相等的实根r 1,r 2

x r x r e c e c y 2121+=

有二重根r 1=r 2

y =(c 1+c 2x)x

r e 1

有一对共轭复根

β

αβαi r i r -=+=21

y =e αx

(c 1cos βx +c 2sin βx)

在解决二阶常系数线性微分方程的解的问题中,我们用到了一个转换,也就是令rx

e y =,那我们为什么会想到它呢?为什么要用它不用其它的呢?首先观察微分方程

022=++qy dx dy p dx y d ,由4个部分构成,2

2dx y

d ,p 与dx dy 相乘,q 与y 相乘,然后相加等于0,我们可以观察到他们的特点就是2

2dx y d ,dx dy

,y 都是与常数因子相乘,然后相加等于0,如

果能找到一个函数y ,其22dx

y d ,dx dy

,y 之间只相差一个常数因子,综合我们以前所学过的

知识,一个函数的一阶,二阶倒数和该函数只相差一个常数,这样的函数有可能是方程(3.1)的特解,在初等函数中,我们所学过的指数函数e rx

,就符合上述要求,于是我们令rx

e y =,

从而进行下去,得到了结果。

求二阶常系数线性非齐次方程

+p +qy =f(x) (3.3

的通解,只要先求出其对应的齐次方程的通解,再求出其一个特解,而 后相加就得到非齐

次方程的通解,而且对应的齐次方程的通解的解法,前面已经解决,因此下面要解决的问题

是求方程(3.3)的一个特解。

方程(3.1)的特解形式,与方程右边的f(x)有关,这里只就f(x)的两种常见的形式进行讨论。 一.f(x)=p n (x)e

αx

,其中p n (x)是n 次多项式,我们先讨论当a 0=时,即当

f(x)=p n (x )时方程

)(22x p qy dx

dy

p dx y d n =++) (3.4

的一个特解。

(1)如果q≠0,我们总可以求得一n 次多项式满足此方程,事实上,可设特解y ~=Q n (x)

=a 0x n

+a 1x

n -1

+…+a n

,其中a 0,a 1,…a n 是待定常数,将y ~

及其导数代入方程(3.4),得方程左右两边都是n 次多项式,比较两边x 的同次幂系数,就可确定常数a 0,a 1,…a n 。

例1. 求32222-=++

x y dx

dy

dx y d 的一个特解。 解 自由项f(x)=x 2

-3是一个二次多项式,又q =2≠0,则可设方程的特解为

2120~a x a x a y ++=

求导数 1

02~a x a y +='

02~a y =''

代入方程有20a x 2

+(20a +21a )x +(20a +a 1+2a 2)=x 2

-3比较同次幂系数

解得

所以特解y ~

=21x 2-21x -4

7

(2)如果q =0,而p≠0,由于多项式求导一次,其次数要降低一次,此时y ~

=Q n (x)不能满足方程,但它可以被一个(n +1)次多项式所满足,此时我们可设

y ~=xQ n (x)=a 0x n +1+a 1x n +…+a n

代入方程(7.4),比较两边系数,就可确定常数a 0,a 1,…a n 。

例2. 求方程234

222+=+x dx

dy

dx y d 的一个特解。 解 自由项 f(x)=3x 2

+2是一个二次多项式,又q =0,p =4≠0,故设特解

x a x a x a y 2130~++=

求导数

212023~a x a x a y ++='

1

026~a x a y +=''

代入方程得

120a x 2

+(8a 1+6a 0)x +(2a 1+4a 2)=3x 2

+2,比较两边同次幂的系数

解得

所求方程的特解x x x y 32

19

16341~

23+-=

(3)如果p =0,q =0,则方程变为2~

y =p n (x),此时特解是一个(n +2)次多项式,可设 y ~=x 2Q n (x),代入方程求得,也可直接通过两次积分求得。

下面讨论当a ≠0时,即当f(x)=p n (x)e αx

时方程

)(~

2x p qy dx

dy p y n =++e αx (3

的一个特解的求法,方程(3.4)与方程(3.5)相比,只是其自由项中多了一个指数函数因子e αx ,如果能通过变量代换将因子e αx

去掉,使得(3.5)化成(3.4)式的形式,问题即可解决,

为此设y =ue αx ,其中u =u(x)是待定函数,对y =ue αx

,求导得

=dx dy ax ax aue dx

du e + 求二阶导数ax ax ax ue a dx

du ae dx u d e dx y d 222

222++=

代入方程(3.5)得

ax ax ax

que au dx

du pe u a dx du a dx u d e +++++][]2[222=p n (x)e αx

消去e αx

)()()2(2

22x p u q pa a dx

du p a dx u d n =+++++ (3

由于(3.6)式与(3.4)形式一致,于是按(3.4)的结论有:

(1)如果02

≠++q pa a ,即a 不是特征方程r 2

+pr +q =0的根,则可设(3.6)的特解u =Q

n

(x),从而可设(3.5)的特解为

y ~=Q n (x)e αx

(2) 如果02

=++q pa a ,而02≠+p a ,即α是特征方程r 2+pr +q =0的单根,则可

设(3.6)的特解u =xQ n (x),从而可设(3.5)的特解为

y ~=xQ n (x)ax e

如果02

=++q pa r ,且02=+p a ,此时a 是特征方程r 2

+pr +q =0的重根,则

可设(3.6)的特解u =x 2

Q n (x),从而可设(3.5)的特解为

y ~=x 2Q n (x)ax e

例3. 求下列方程具有什么样形式的特解

(1)x e y dx

dy

dx y d 32265

=++

(2)x e y dx

dy

dx y d 222365

-=++

(3)x e a x y dx

dy

a

dx y d -+-=++)3(222

解(1)因a =3不是特征方程r 2

+5r +6=0的根,故方程具有形如

x e a y 3

0~=的特解。

(2)因a =-2是特征方程r 2

+5r +6=0的单根,故方程具有形如

x e a x a x y 210)(~-+=的特解。

(3)因a =-1是特征方程r 2

+2r +1=0的二重根,所以方程具有形如

x e a x a x a x y -++=)(~21202的特解。

结论

在一阶二阶微分方程中,有一个函数因子rx

e y =,都是通过它的转换,巧妙的解决了该类微分方程的求解,其实对于n 阶微分方程,首先观察微分方程

0)(.......)1(11=+++--x t a dt x

d a dt x d n n n n

n ,n n dt x d ,)1(a 与11--n n dt x d 相乘,.........,)(t a 与x 相乘,然后相加等于0,我们可以观察到他们的特点就是n n dt x d ,1

1--n n dt x

d ,.......,x 都是与常数因子相乘,然后相加等于0,如果能找到一个函数y ,其n n dt x d ,1

1--n n dt x

d ,.....,x 之间也只相差一

个常数因子,综合我们以前所学过的知识,一个函数的一阶,二阶倒数和该函数只相差一个常数,在初等函数中,我们所学过的指数函数e rx

,就符合上述要求,于是我们令rx e y ,

从而进行下去,也能得出结果。 .

变系数线性常微分方程的求解

变系数线性常微分方程的求解 张慧敏,数学计算机科学学院 摘要:众所周知,所有的常系数一阶、二阶微分方程都是可解的,而变系数 二阶线性微分方程却很难解,至今还没有一个普遍方法。幂级数解法是一个非常有效的方法,本文重点讨论二阶变系数线性常微分方程的解法,从幂级数解法、降阶法、特殊函数法等方面探究了二阶微分方程的解法,简单的介绍了几种高阶微分方程的解法,并讨论了悬链线方程等历史名题。 关键词:变系数线性常微分方程;特殊函数;悬链线方程;幂级数解法 Solving linear ordinary differential equations with variable coefficients Huimin Zhang , School of Mathematics and Computer Science Abstract:As we know, all of ordinary differential equations of first, second order differential equations with constant coefficients are solvable. However, the linear differential equations of second order with variable coefficients are very difficult to solve. So far there is not a universal method. The method of power-series solution is a very efficient method. This article focuses on solving linear ordinary differential equations of second order with variable coefficients, and exploring the solution of in terms of power-series solution, the method of reducing orders, the method of special functions. Also, this paper applies the above methods to solve several linear differential equations of higher order and especially discusses the famous catenary equation. Key words:Linear ordinary differential equations with variable coefficients; Special Functions; catenary equation; Power Series Solution.

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

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

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

附录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)

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

附录A 线性常微分方程 本课程的研究内容与常微分方程理论有非常密切的联系,因此在本附录里,我们将对线性常微分方程的知识——包括解的存在性、解的结构和求解方法做一些回顾和总结。 把包含未知函数和它的j 阶导数()j y (的方程称为常微分方程。线性常微分方程的标准形式 ()(1)110()()'()()n n n y p x y p x y p x y f x --++++=L (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 不恒等于零,则上式等价于 而()'ln 'y y y =,从而(A.3)的通解为 ()d ()p x x y x Ce -?= ( A.4) 对于非齐次一阶线性常微分方程(A.2),在其两端同乘以函数()d p x x e ?

注意到上面等式的左端 因此有 两端积分 其中C 是任意常数。进一步有 综上有如下结论 定理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 -?加上函数()d ()d *()()d p x x p x x y x e e f x x -??=?。容易验证,*()y x 是方程(A.1)的一个特解。这符合线性方程解的结构规律。 例1 求解一阶常微分方程 解 此时()2()1p x f x =-=,,由(A.5)式,解为 其中C 是任意常数。 A.2 二阶线性常微分方程 将具有以下形式的方程 "()'()()y p x y q x y f x x I ++=∈,, (A.6) 称为二阶线性常微分方程,其中(),(),()p x q x f x 都是变量x 的已知连续函数。称 "()'()0y p x y q x y x I ++=∈,, (A.7) 为与(A.6)相伴的齐次方程. A .2.1 二阶线性微分方程解的结构 首先讨论齐次方程(A.7)解的结构。

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

附录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 -?加上函数

一阶线性微分方程组

第4章 一阶线性微分方程组 一 内容提要 1. 基本概念 一阶微分方程组:形如 ??? ????? ???===) ,,,,( ),,,,(),,,,(2121222111 n n n n n y y y x f dx dy y y y x f dx dy y y y x f dx dy ΛΛΛΛΛ (3.1) 的方程组,(其中n y y y ,,,21Λ是关于x 的未知函数)叫做一阶微分方程组。 若存在一组函数)(,),(),(21x y x y x y n Λ使得在[a,b]上有恒等式 ),,2,1))((,),(),(,() (21n i x y x y x y x f dx x dy n i i ΛΛ==成立,则 )(,),(),(21x y x y x y n Λ称为一阶微分方程组(3.1)的一个解 含有n 任意常数n C C C ,,,21Λ的解 ?????? ?===) ,,,,( ),,,,(),,,,(21321222111n n n n C C C x y C C C x y C C C x y ΛΛΛΛΛ??? 称为(3.1)通解。如果通解满方程组 ???????=Φ=Φ=Φ0 ),,,,,,,,( 0),,,,,,,,(0),,,,,,,,(21212121221211n n n n n n n C C C y y y x C C C y y y x C C C y y y x ΛΛΛΛΛΛΛΛ 则称这个方程组为(3.1)的通积分。 满足初始条件,)(,,)(,)(0020021001n n y x y y x y y x y ===Λ的解,叫做初值问题的解。

(整理)常系数线性微分方程的解法

常系数线性微分方程的解法 摘要:本文对常系数线性方程的各种解法进行分析和综合,举出了每个方法的例题,以便更好的掌握对常系数线性微分方程的求解. 关键词:特征根法;常数变易法;待定系数法 Method for solving the system of differential equation with Constant Coefficients Linear Abstract: Based on the linear equations with constant coefficients of analysis and synthesis method, the method of each sample name, in order to better grasp of the linear differential equation with constant coefficients of the solution. Key Words: Characteristic root ;Variation law ;The undetermined coefficient method 前言:常系数性微分方程因形式简单,应用广泛,解的性质及结构已研究的十分清楚,在常微分方程中占有十分突出的地位。它的求解是我们必须掌握的重要内容之一,只是由于各种教材涉及的解法较多,较杂,我们一般不易掌握,即使掌握了各种解法,在具体应用时应采用哪种方法比较适宜,我们往往感到困难。本文通过对一般教材中涉及的常系数线性微分方程的主要解法进行分析和比较,让我们能更好的解常系数线性微分方程。 1.预备知识 复值函数与复值解 如果对于区间a t b ≤≤中的每一实数t ,有复值()()()z t t i t ?ψ=+与它对应,其中()t ?和()t ψ是在区间a t b ≤≤上定义的实函数,1i =-是虚数单位,我们就说在区间a t b ≤≤上给定了一个复值函数()z t .如果实函数()t ?,()t ψ当t 趋于 0t 时有极限,我们就称复值函数()z t 当t 趋于0t 时有极限,并且定义

线性方程组解的判定与解的结构

***学院数学分析课程论文 线性方程组解的判定与解的结构 院系数学与统计学院 专业数学与应用数学(师范) 姓名******* 年级 2009级 学号200906034*** 指导教师 ** 2011年6月

线性方程组解的判定与解的结构 姓名****** (重庆三峡学院数学与计算机科学学院09级数本?班) 摘 要:线性方程组是否有解,用系数矩阵和增广矩阵的秩来刻画.在方程组有解且有 多个解的情况下,解的结构就是了解解与解之间的关系. 关键词:矩阵; 秩; 线性方程组; 解 引言 通过系数矩阵和增广矩阵的秩是否相同来给出判定线性方程组的解的判别条件.在了解了线性方程组的判别条件之后,我们进一步讨论解的结构.对于齐次线性方程组,解的线性组合还是方程组的解.在线性方程组有无穷个解时可用有限多个解表示出来.另外以下还涉及到线性方程组通解的表达方式. 1 基本性质 下面我们分析一个线性方程组的问题,导出线性方程组有解的判别条件. 对于线性方程组 1111221121122222 1122n n n n s s sn n s a x a x a x b a x a x a x b a x a x a x b ++???+=??++???+=???????++???+=? (1) 引入向量 112111s αααα??????=?????????,122222s αααα??????=?????????,…12n n n sn αααα??????=????????? ,12s b b b β?? ?? ??=??????? ?? 方程(1)可以表示为 1122n n x x x αααβ++???+= 性质 线性方程组⑴有解的充分必要条件为向量β可以表成向量组α1,α2,…,αn 的线性组合. 定理1 线性方程组⑴有解的充分必要条件为它的系数矩阵

常系数线性微分方程的解法

常系数线性微分方程的解法 摘 要:本文主要介绍了常系数线性微分方程的解法.着重讨论利用代数运算和微分运算来求常系数齐次线性微分方程和非齐次线性微分方程的通解. 关键词:复值函数与复值解;欧拉方程;比较系数法;拉普拉斯变换法 The Solution of Linear Differential Equation with Constant Coefficients Abstract :The solutions of linear differential equation with constant coefficients are introduced in this article. And using the algebraic operation and differential operation to solv the general solution of homogeneous linear differential equation and nonhomogeneous linear differential equation are discussed emphatically. Key Words :complex flnction and complex answer; euler equation;the method of coefficients comparison; the method of laplace transformation. 前言 为了让我们更多的认识和计算常系数线性微分方程,本文通过对复值函数和复值解以及常系数线性微分方程和欧拉函数的简单介绍,进而简单讨论了常系数线性微分方程的解法,以此来帮助我们解决常系数线性微分方程的解. 1. 预备知识 1.1复值函数与复值解 如果对于区间a t b ≤≤中的每一个实数t ,有复数()()()z t t i t ?ψ=+与它对应,其中 ()t ?和()t ψ是在区间a t b ≤≤上定义的实函数,i =是虚数单位,我们就说在区间 a t b ≤≤上给定了一个复值函数()z t .如果实函数()t ?,()t ψ当t 趋于0t 时有极限,我们 就称复值函数()z t 当t 趋于0t 时有极限,并且定义 lim ()lim ()lim ()t t t t t t z t t t ?ψ→→→=+. 如果0 0lim ()()t t z t z t →=,我们就称()z t 在0t 连续.显然,()z t 在0t 连续相当于()t ?,()t ψ在0 t 连续.当()z t 在区间a t b ≤≤上每点都连续时,就称()z t 在区间a t b ≤≤上连续.如果极

阶线性微分方程解的结构

阶线性微分方程解的结 构 公司内部编号:(GOOD-TMMT-MMUT-UUPTY-UUYY-DTTI-

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

假设()y x 不恒等于零,则上式等价于 而()'ln 'y y y =,从而()的通解为 ()d ()p x x y x Ce -?= ( ) 对于非齐次一阶线性常微分方程(),在其两端同乘以函数()d p x x e ? 注意到上面等式的左端 因此有 两端积分 其中C 是任意常数。进一步有 综上有如下结论 定理 假设()()p x f x I 和在上连续,则一阶线性非齐次常微分方程()的通解具有如下形式 ()d ()d ()d ()()d p x x p x x p x x y x Ce e e f x x --???=+?‘ () 其中C 是任意常数。 观察()式和()式,我们发现一阶线性非齐次常微分方程()的解等于一阶线性齐次常微分方程()的通解()d p x x Ce -?加上函数 ()d ()d *()()d p x x p x x y x e e f x x -??=? 。容易验证,*()y x 是方程()的一个特解。这符合线性方程解的结构规律。 例1 求解一阶常微分方程

第三章 一阶线性微分方程组 第四讲 常系数线性微分方程组的解法(1)

第四讲 常系数线性微分方程组的解法(4课时) 一、目的与要求: 理解常系数线性微分方程组的特征方程式, 特征根, 特征向量的概念, 掌 握常系数线性微分方程组的基本解组的求法. 二、重点:常系数线性微分方程组的基本解组的求法. 三、难点:常系数线性微分方程组的特征方程式, 特征根, 特征向量的概念. 四、教学方法:讲练结合法、启发式与提问式相结合教学法. 五、教学手段:传统板书与多媒体课件辅助教学相结合. 六、教学过程: 1 新课引入 由定理3.6我们已知道,求线性齐次方程组(3.8)的通解问题,归结到求其基本解组. 但是对于一般的方程组(3.8),如何求出基本解组,至今尚无一般方法. 然而对于常系数线性齐次方程组 dY AY dx = (3.20) 其中A 是n n ?实常数矩阵,借助于线性代数中的约当(Jordan)标准型理论或矩阵指数,可以使这一问题得到彻底解决. 本节将介绍前一种方法,因为它比较直观. 由线性代数知识可知,对于任一n n ?矩阵A ,恒存在非奇异的n n ?矩阵T ,使矩阵 1T AT -成为约当标准型. 为此,对方程组(3.20)引入非奇异线性变换 Y TZ = (3.21) 其中()(,1,2,,),ij T t i j n == det 0T ≠,将方程组(3.20)化为 1dZ T ATZ dx -= (3.22) 我们知道,约当标准型1 T AT -的形式与矩阵A 的特征方程 11121212221 2 det()0n n n n nn a a a a a a A E a a a λ λλλ ---= =-

的根的情况有关. 上述方程也称为常系数齐次方程组(3.20)的特征方程式.它的根称为矩阵 A 的特征根. 下面分两种情况讨论. (一) 矩阵A 的特征根均是单根的情形. 设特征根为12,,,,n λλλ 这时 12 1 00 n T AT λλλ-????? ?=?????? 方程组(3.20)变为 11122 200n n n dz dx z dz z dx z dz dx λλλ?????????????? ????????= ???????????????? ?????? (3.23) 易见方程组(3.23)有n 个解 1110(),00x Z x e λ????????=???????? 220010(),,()0001n x x n Z x e Z x e λλ???????????? ????==???????????????? 把这n 个解代回变换(3.21)之中,便得到方程组(3.20)的n 个解 12()i i i i x x i i ni t t Y x e e T t λλ?? ????==?????? (1,2,,)i n =

线性方程组解的结构

线性方程组解的结构 11111221n n b a x a x a x =++???+ 22112222n n b a x a x a x =++???+ 33113223n n b a x a x a x =++???+ ………………………………… 1122n n n nn n b a x a x a x =++???+ 表示从变量12 ,n x x x ???到变量12,n b b b ???的线性变换,其中ij a 是常数。确 定了线性变换,它的系数所构成的矩阵(系数矩阵)也就确定,线性变换根矩阵是一一对应的关系。 上式可以表示为以向量x 为未知元的向量方程: Ax=b 线性方程组如果是有解的,称它是相容的,否则称为不相容。 一、 定理4:N 元线性方程组Ax=b (1) 无解的充要条件是R(A)

(2) 若R(A)=R(B),则进一步把B 化成最简型,而对于齐次线性 方程组,则把系数矩阵A 化成最简型。 (3) 设R(A)=R(B)=r ,把行最简型中r 个非0行的非0首个元素所对应的未知数取做非自由未知数,其他的元素取做自由未知数。带入原方程,就可以得到一个关于自由为未知量的表达式。 三、 齐次线性方程组求解步骤:Ax=0 (1) 根据R(A)与n (变量个数)来判断解的结构: A. R(A)

常系数线性微分方程的解的结构分析

常系数线性微分方程的解的结构分析 【 摘要】在参考和总结了许多场系数线性微分方程的解法的基础上,本文总结了一些常系数微分方程的解的解法,并针对一类常系数线性微分方程的已有结论给予证明,以解给予一些结论证明思路,以及一些实例,并向高阶推广。 【关键词 】常系数 线性 微分方程 结构 一阶常系数齐次线性微分方程 0=+ax dt dx , (1.1) 的求解 上式可以改写为 adt x dx -= , (1.2) 于是变量x 和t 被分离,再将两边积分得 c at x +-=ln , (1.3) 这里的c 为常数。又由对数的定义,上式可以变为 at ce x -= , (1.4) 其中c= , 因为x=0也是方程的解,因此c 可以是任意常数。 这里首先是将变量分离,然后再两边积分,从而求出方程的解。这便要方程式可以分离变量的,也就是变量分离方程。 一阶常系数微分方程 )()(x Q y x P dx dy += , (2.1) 其中P (x ),Q(x)在考虑的区间上式连续函数,若Q (x )=0 ,上式就变为 y x P dx dy )(= , (2.2) 上式为一阶齐次线性微分方程。还是变量分离方程我们可以参考上面变量分离方程的解法,先进行变量分离得到 dx x P y dy )(= , (2.3) 两边同时积分,得到 ? =dx x p ce y )( , (2.4) 这里c 是常数。 若Q (x )≠ 0 , 那么上式就变成了 一阶非齐次线性微分方程。 我们知道一阶齐次线性微分方程是一阶常微分方程的一种特殊情况,那么可以设想将一阶

齐次线性微分方程的解 ? =dx x p ce y )( , (2.5) 中的常数c 变易成为待定的函数c (x ),令 ?=dx x p e x c y )()( , (2.6) 微分之,就可以得到 ?+?=dx x p dx x p e x P x c e dx x dc dx dy )()()()()( , (2.7) 以(2.7),(2.6)代入2.1,得到 )()()()()()()()()(x Q e x c x p e x P x c e dx x dc dx x p dx x p dx x p +?=?+?,(2.8) 即 ?=-dx x p e x Q dx x dc )()() (, 积分后得到 c (x )=c dx e x Q dx x p +?? -)()( , (2.9) 这里c 是任意常数,将上式代入(2.6)得到方程(2.1)的通解 ))(()()(c dx e x Q e y dx x p dx x p +? ? =?- (2.91) 在上面的一阶线性微分方程中,是将一阶齐次线性微分方程中的通解中的常数c 变成c(x) ,常数变易法一阶非齐次线性微分方程的解, 感觉这个方法之所以用x 的未知函数u(x)替换任意常数C,是因为C 是任意的,C 与x 形成函数关系,要确定C,需要由初始条件确定,一个x,确定一个C,也就形成一对一或多对多的映射,也就是函数关系,而这里的C 是任意的,也就可以用一个未知的,也就是任意的函数u(x)来代替,进而求得非齐次线性微分方程的解。这种将常数变异为待定函数的方法,我们通常称为常数变易法。常数变易法实质也是一种变量变换的方法,通过变换(2.6可将方程(2.1)化为变量分离方程。 二阶常系数线性微分方程 (1)二阶常系数线性齐次方程 022=++qy dx dy p dx y d (3.1) 其中p 、q 是常数,我们知道,要求方程(3.1)的通解,只要求出其任意两个线性无关的特 解y 1,y 2就可以了,下面讨论这样两个特解的求法。 我们先分析方程(3.1)可能具有什么形式的特解,从方程的形式上来看,它的特点是22dx y d ,

(整理)二阶常系数线性微分方程的解法版.

第八章 8.4讲 第四节 二阶常系数线性微分方程 一、二阶常系数线形微分方程的概念 形如 )(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 2 2 sin ,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 为常数)和它的各阶导数都只差一个常数因子,

齐次和非齐次线性方程组的解法整理

线性方程组解的结构(解法) 一、齐次线性方程组的解法 【定义】r (A )= r 时,若()r A n ≤,则存在齐次线性方程组的同解方程组; 若()r A n >,则齐次线性方程组无解。 1、求AX = 0(A 为m n ?矩阵)通解的三步骤 (1)?? →A C 行 (行最简形);写出同解方程组CX =0. (2) 求出CX =0的基础解系,,,n r -12ξξξ; (3) 写出通解n r n r k k k --=+++1122X ξξξ其中k 1,k 2,…, k n-r 为任意常数.

常系数线性方程组基解矩阵的计算

常系数线性方程组基解矩阵的计算

常系数线性方程组基解矩阵的计算 董治军 (巢湖学院数学系,安徽巢湖238000) 摘要:微分方程组在工程技术中的应用时非常广泛的,不少问题都归结于它的求解问题,基解矩阵的存在和具体寻求是不同的两回事,一般齐次线性微分方程组的基解矩阵是无法通过积分得到的,但当系数矩阵是常数矩阵时,可以通过方法求出基解矩阵,这时可利用矩阵指数exp A t,给出基解矩阵的一般形式,本文针对应用最广泛的常系数线性微分方程组,结合微分方程,线性代数等知识,讨论常系数齐次线性微分方程的基解矩阵的几个一般的计算方法. 关键词;常系数奇次线性微分方程组;基解矩阵;矩阵指数 Calculation of Basic solution Matrix of

Linear Homogeneous System with Constant Coefficients Zhijun Dong (Department of Mathematics, Chaohu College Anhui, Chaohu) Abstract: Differential equations application in engineering technology is very extensive, when many problems are attributable to its solving problem, base solution matrix existence and specific seek is different things, general homogeneous linear differential equations is not the base solution matrix by integral get, but when coefficient matrix is constant matrix, can pass out the base solution matrix method, then are available matrix exponential t, the general form base solution matrix, the paper discusses the most widely used differential equations with constant coefficients, combined with differential equations, linear algebra, discuss knowledge of homogeneous linear differential equation with constant coefficients of base solution matrix several general calculation method. Keyword: linear homogeneous system with constant coefficients; matrix of basic solutions; matrix exponent 引言: 线性微分方程组的求解历来是常微分方程的重点,根据线性微分方程组的解的结构理论,求解线性微分方程组的关键在于求出对应齐次线性微分方程组的基解矩阵,本文主要讨论齐次线性微分方程组 X ’=AX ★ 的基解矩阵的计算问题,这里A 是n n ?常数矩阵. 一.矩阵指数exp A 的定义和性质: 1.矩阵范数的定义和性质 定义:对于n n ?矩阵A =ij a ???? n ×n 和n 维向量X =()1,...,T n X X 定义A 的范数为A =,1 n ij i j a =∑ ,X =1 n i i x =∑ 设A ,B 是n ×n 矩阵,x ,y 是n 维向量,易得下面两个性质:

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

附录A 线性常微分方程 本课程的研究内容与常微分方程理论有非常密切的联系,因此在本附录里,我们将对线性常微分方程的知识——包括解的存在性、解的结构和求解方法做一些回顾和总结。 把包含未知函数和它的j 阶导数()j y (的方程称为常微分方程。线性常微分方程的标准形式 ()(1)110()()'()()n n n y p x y p x y p x y f x --++++=L (A.1) I 上满A.1)的) 当f ) 假设()y x 不恒等于零,则上式等价于 而()'ln 'y y y =,从而(A.3)的通解为 ()d ()p x x y x Ce -?= ( A.4) 对于非齐次一阶线性常微分方程(A.2),在其两端同乘以函数()d p x x e ?

注意到上面等式的左端 因此有 两端积分 其中C 是任意常数。进一步有 综上有如下结论 定理 A.1 假设()()p x f x I 和在上连续,则一阶线性非齐次常微分方程(A.1)的通解 阶函数y ) 称为二阶线性常微分方程,其中(),(),()p x q x f x 都是变量x 的已知连续函数。称 "()'()0y p x y q x y x I ++=∈,, (A.7) 为与(A.6)相伴的齐次方程. A .2.1 二阶线性微分方程解的结构 首先讨论齐次方程(A.7)解的结构。

定理 A.2 如果函数12()()y x y x 与是线性齐次方程(A.7)的两个解,则函数1122()()y c y x c y x =+仍为该方程的解,其中12,c c 是任意的常数。 定理1 说明齐次线性常微分方程(A.7)的解如果存在的话,一定有无穷多个。为了说明齐次线性常微分方程(A.7)通解的结构,首先给出函数线性无关的定义。 定义A.1设函数12(),(),,()n y x y x y x L 是定义在区间I 上的n 个函数,如果存在n 个 不全为零的常数12,,n k k k L , ,使得1122()()()0n n k y x k y x k y x ++=L 在区间I 上恒成立, (,a (A.7)。 则y 其基础解系。 关于二阶线性非齐次常微分方程(A.6)的通解,有如下结论 定理A.5 若函*()y x 是方程(A.6)的一个特解,()Y x 是方程(A.6)相伴的齐次方程的通解,则()()*()y x y x Y x =+是二阶线性非齐次常微分方程(A.6)的通解。 从定理A.4,A.5可以得到求解二阶线性非齐次常微分方程(A.6)的通解的一般步骤: (1)求解与(A.6)相伴的齐次方程(A.7)的线性无关的两个特解12()()y x y x 与,得该

相关文档
最新文档