第四章 高阶微分方程
- 格式:ppt
- 大小:2.15 MB
- 文档页数:66
第四章 高阶微分方程§4.1 线性微分方程的一般理论习题4.11.设)(t x 和)(t y 是区间[]b a ,上的连续函数,证明:若在区间[]b a ,上有≠)()(t y t x 常数或≠)()(t x t y 常数,则)(t x 和)(t y 在区间[]b a ,上线性无关.(提示:用反证法) 证明 )(t x 和)(t y 是区间[]b a ,上线性相关,则存在不全为0的常数21,c c 使得0)()(21≡+t y c t x c ,[]b a t ,∈,若)0(,021≠≠c c 或得12)()(c c t y t x -≡(或21)()(c c t x t y -≡)[]b a t ,∈∀成立。
与假设矛盾,故)(t x 和)(t y 在区间[]b a ,上线性无关.2.证明非齐次线性方程的叠加原理:设)(1t x ,)(2t x 分别是非齐次线性方程)()()(1111t f x t a dt xd t a dt x d n n n n n =+++-- (1) )()()(2111t f x t a dtxd t a dt x d n n n nn =+++-- (2) 的解,则)()(21t x t x +是方程)()()()(21111t f t f x t a dtxd t a dt x d n n n n n +=+++-- (3) 的解.证明 因为)(1t x ,)(2t x 分别是方程(1)、(2)的解,所以)()()(1111111t f x t a dt x d t a dt x d n n n n n =+++-- , )()()(2212112t f x t a dtx d t a dt x d n n n nn =+++-- , 二式相加得,)()())(()()()(21211211121t f t f x x t a dt x x d t a dt x x d n n n n n +=++++++-- ,即)()(21t x t x +是方程(3)的解.3.(1).试验证022=-x dt x d 的基本解组为tt e e -,,并求方程t x dtx d cos 22=-的通解。
第四章高阶微分方程[教学目标]1. 理解高阶线性微分方程的一般理论,n阶齐次(非齐次)线性微分方程解的性质与结构,熟练掌握n阶常系数齐次线性微分方程的待定指数函数解法。
2.掌握n阶非齐次线性微分方程的常数变易法,理解n阶常系数非齐次线性微分方程特解的待定系数法和Laplce变换法。
3.熟练欧拉方程与高阶方程的降阶法和幂级数解法。
4.掌握高阶方程的应用。
[教学重难点]重点是线性微分方程解的性质与结构,高阶方程的各种解法。
难点是待定系数法求特解。
[教学方法] 讲授,实践。
[教学时间] 16学时[教学内容]线性微分方程的一般理论,齐次(非齐次)线性微分方程解的性质与结构,非齐次线性微分方程的常数变量易法;常系数线性方程与欧拉方程的解法,非齐线性方程的比较系数法与拉氏变换法;高阶方程的降阶法和幂级数解法及高阶方程的应用。
[考核目标]1.理解高阶线性微分方程的一般理论,能够求解高阶常系数线性微分方程。
2.掌握n阶非齐次线性微分方程的常数变易法。
3.n阶常系数非齐次线性微分方程特解的待定系数法和Laplce变换法。
4.熟练高阶方程的降阶法和幂级数解法及高阶方程的应用。
§4.1线性微分方程的一般理论4.1.1引言讨论n阶线性微分方程1111()()()()n n n n n n d x d x dxa t a t a t x f t dt dt dt---++++= (4.1) 其中()(1,2,,)i a t i n = 及()f t 都是区间a t b ≤≤上的连续函数 如果()0f t ≡,则方程(4.1)变为:1111()()()0n n n n n n d x d x dxa t a t a t x dt dt dt---++++= (4.2) 称它为n 阶齐线性微分方程,而称一般的方程(4.1)为n 阶非齐线性微分方程,并且通常把方程(4.2)叫对应于方程(4.1)的齐线性方程。
定理1 如果()(1,2,,)i a t i n = 及()f t 都是区间a t b ≤≤上的连续函数,则对于任一[]0,t a b ∈ (1)(1)000,,,n x x x - ,方程(4.1)存在唯一解()x t ϕ=,定义于区间a t b ≤≤上,且满足初始条件:1(1)(1)0000001()()(),,,n n n d t d t t x x x dt dtϕϕϕ---=== (4.3) 从这个定理可以看出,初始条件唯一地确定了方程(4.1)的解,而且这个解在所有()(1,2,,)i a t i n = 及()f t 连续的整个区间a t b ≤≤上有定义。
第四章 高阶微分方程一.一般的概念和性质 一般的n 阶线性微分方程具有如下的形式:)()()()(1111t f x t a dt dx t a dt xd t a dtx d n n n n n n=++++--- ,(4.1)0)()()(1111=++++---x t a dtdx t a dtxdt a dtx d n n n n nn .(4.2)其中),,2,1()(n i t a i =和)(t f 都是某区间],[b a 上的连续函数. (4.2)称为齐线性(微分)方程, (4.1)称为非齐线性(微分)方程. (4.2)也称为(4.1)的对应齐方程.1.函数组的线性相关与线性无关. 区间],[b a 上的k 个函数)(),(),(21t x t x t x k 称为是线性相关的, 如果存在不全为零的常数k c c c ,,21, 使得在],[b a 上恒成立0)()()(2211≡+++t x c t x c t x c k k .如果这样的常数不存在, )(),(),(21t x t x t x k 称为是线性无关的.2.(伏)郎斯基行列式. 如果函数)(),(),(21t x t x t x k 还有直到1-k 阶的导数, 行列式)()()()()()()()()()](),(),([)()1()1(2)1(1212121t x t x t x t x t x t x t x t x t x t x t x t x W t W k kk k k k k ---'''==称为这些函数的伏郎斯基行列式或郎斯基行列式.典型例题:已知123()4,()4,()xy x y x x y x e x ==-=+是某一三阶齐线性方程的解, 试求 )(),(21x y x y 和)(3x y 的伏朗斯基行列式123[,,]()W y y y x . (见模拟试题)3.基本解组. 齐线性方程(4.2)的n 个线性无关解称为(4.2)的一个基本解组。
第四章高阶微分方程教学安排说明章节题目:§4.1 线性微分方程的一般理论;§4.2 常系数线性方程的解法§4.3 高阶方程的降阶和幂级数解法学时分配:共16学时。
§4.1 线性微分方程的一般理论(4学时)§4.2 常系数线性方程的解法(8时§4.3 高阶方程的降阶和幂级数解法(4学时)本章教学目的与要求:1. 理解高阶线性微分方程的一般理论,n阶齐次(非齐次)线性微分方程解的性质与结构,熟练掌握n阶常系数齐次线性微分方程的待定指数函数解法。
2.掌握n阶非齐次线性微分方程的常数变易法,理解n阶常系数非齐次线性微分方程特解的待定系数法和Laplce变换法。
3.熟练欧拉方程与高阶方程的降阶法和幂级数解法。
4.掌握高阶方程的应用。
教学重点与难点:重点是线性微分方程解的性质与结构,高阶方程的各种解法。
难点是待定系数法求特解。
教学内容线性微分方程的一般理论,齐次(非齐次)线性微分方程解的性质与结构,非齐次线性微分方程的常数变量易法;常系数线性方程与欧拉方程的解法,非齐线性方程的比较系数法与拉氏变换法;高阶方程的降阶法和幂级数解法及高阶方程的应用。
考核目标1.理解高阶线性微分方程的一般理论,能够求解高阶常系数线性微分方程。
2.掌握n阶非齐次线性微分方程的常数变易法。
3.n阶常系数非齐次线性微分方程特解的待定系数法和Laplce变换法。
4.熟练高阶方程的降阶法和幂级数解法及高阶方程的应用。
课 堂 教 学 方 案课程名称:§4.1 线性微分方程的一般理论授课时数:4学时授课类型:理论课 教学方法与手段:讲授法教学目的与要求:理解高阶线性微分方程的一般理论,n 阶齐次(非齐次)线性微分方程解的性质与结构,熟练掌握n 阶常系数齐次线性微分方程的待定指数函数解法。
教学重点、难点:线性微分方程解的性质与结构教学内容§4.1 线性微分方程的一般理论4.1.1 引言一般形式:1111()()()()n n n n n n d x d xdxa t a t a t x f t dt dtdt---++++= (4.1) 其中(),1,2,,i a t i n =及()f t 都是a t b ≤≤上的连续函数。
第四章 高阶微分方程例4-1 设()()t x t x t x n n 11),(,,+ 是n 阶非齐次线性方程()()()()()t f x t a x t a x n n n =+++- 11 (1)的在区间[]b a ,上的1+n 个线性无关的解,则方程(1)在区间[]b a ,上的任何解()t x ,都可以表示为()t x ()()()t x C t x C t x C n n n n 1111+++++=其中 111=++++n n C C C ,反过来,若()()t x t x t x n n 11),(,,+ 是方程(1)在区间[]b a ,上的1+n 个线性无关的解,则()()()t x C t x C t x C n n n n 1111+++++必为方程(1)在区间[]b a ,上的解,其中111=++++n n C C C 。
证 构造函数11211,,,+++---n n n n x x x x x x ,显然,它们是方程(1)所对应的齐次 方程()()()()011=+++-x t a x t a x n n n (2)的n 个解,并且在区间[]b a ,上线性无关。
若不然,假设存在一组不全为零的数n C C C ,,,21 使得[]b a t ,∈()()()01122111=-++-+-+++n n n n n x x C x x C x x C , 即()01212211=+++-++++n n n n x C C C x C x C x C ,与()()t x t x t x n n 11),(,,+ 是区间[]b a ,上的1+n 个线性无关的解相矛盾。
因此,由定理4.7,方程(1)的任何解()t x 都可以表示为()t x ()()()11122111+++++-++-+-=n n n n n n x x x C x x C x x C()12122111+----++++=n n n n x C C C x C x C x C 。
第四章高阶微分方程本章教学目的及要求1、使学生掌握n阶线性齐次方程和非齐次方程解的性质。
熟练掌握通解结构。
2、要求学生熟练掌握常系数线性齐次方程的基本解组的特征根法和常系数线性非齐次方程特解的待定系数法。
熟练掌握常数变易法,掌握较简单的二阶变系数方程的解法。
本章重点线性方程的性质及其解的结构,解常系数线性方程的特征根法和待定系数法及视常数为变数法。
难点函数组线性相关性及判别,一般理论中定理证明二阶变系数线性方程求解。
课时安排:讲授16学时,习题课4学时第四章高阶微分方程本章主要介绍n阶线性微分方程的一般理论和求解方法.我们将把n阶线性微分方程化成等价的一阶线性微分方程组,这样可以把第3章的主要结论自然地应用到本章内容中.即把n阶线性微分方程作为一阶线性微分方程组的特例加以处理,以避免理论推导上的重复.§4.1 线性微分方程的一般理论4.1.1 引言n阶线性微分方程在自然科学与工程技术中有着极其广泛的应用.在介绍线性方程的一般理论之前,先让我们来研究两个实际例子.例1 弹簧振动.图 4-1设一质量为m 的物体B 被系于挂在顶板上一弹簧的末端,(我们将假设弹簧的质量与这一物体的质量比较起来是小的可以忽略不计的),现在来求该物体在外力扰动时的运动微分方程式.当物体B 不受外力扰动时,重力被作用于物体B 上的弹簧的弹力所平衡而处于静止位置,把物体B 的静止位置取为坐标轴x 的原点0,向下方向取为正向,如图4-1的(a ).若有一外力f(t)沿垂直方向作用在物体B 上,那么物体B 将离开静止位置0,如图4-1的(b ),记x=x(t)表物体B 在t 时刻关于静止位置0的位移,于是22,dtx d dx dx 分别表示物体B 的速度和加速度. 由牛顿第二定律F=ma, m 是物体B 的质量,a =22dtx d 是物体B 位移的加速度,而F 是作用于物体B 上的合外力.这时,合外力F 由如下几部分构成.(1)弹簧的恢复力f 1,依虎克定律,弹簧恢复力f 1与物体B 的位移x 成正比,即f 1=-cx式中比例常数c (>0)叫作弹性系数,根据所取的坐标系,恢复力f 1的方向与位移x 的方向相反,所以上式右端添一负号.(2)空气的阻力f 2,当速度不太大时,空气阻力f 2可取为与物体B 位移的速度成正比,亦即f 2=—μdt dx 式中比例常数μ(>0)叫作阻尼系数,式中右边的负号,是由于阻力f 2的方向与物体B 的速度dtdx 的方向相反. (3)外力f(t).因此,我们得到F=-cx-μdt dx +f(t) 从而我们得物体B 在外力f(t)作用下的运动微分方程式).(22t f cs dt dx dtx d m =++μ (4.1) 我们将在本章第4节,详细叙述方程(4.1)所描述的弹簧振动性质.由于方程(4.1)是描述物体B 在外力f(t)经常作用下的运动,所以方程(4.1)亦称为阻尼强迫振动.我们将在本章第4节,详细叙述方程(4.1)所描述的弹簧振动性质.由于方程(4.1)是描述物体B 在外力f(t)经常作用下的运动,所以方程(4.1)亦称为阻尼强迫振动.例2 电振荡在很多无线电设备(如收音机和电视机)中,我们经常见到如图4-2的回路.它由四个元件组成,即电源(设其电动势为E)图4-2电阻R,电感L 以及电容器C.为了简单起见,电容器的电容量我们也用C 表示,它所储藏的电荷量为q.这时电容器的两个极板分别带着等量但符号相反的电荷,极板间的电位差等于 1c E q C =. 此外,当电路中流过交流电时,电容器极板上的电量以及它们的正负符号均随时间发生变化.根据电流定义,这时有dq i dt=. 根据基尔霍夫(Kirchhoff)第二定律,在闭合回路中全部元件的电压的代数和等于零,即0di q E Ri Ldt C ---= 整理后可得di q LRi E dt C ++= (4.2) 考虑到dq i dt=,上式可写成 22d q dq q L R E dt dt C++= (4.3) 于是,得到了关于电荷量q 的方程.如果在式(4.2)两端对t 求导数,并假设E 是常量(直流电压),则可得关于电流的方程220d i di iL Rdt dt C++= (4.4)实验表明,在一定条件下,上述回路中的电流会产生周期振荡,因此我们把上述回路称为电振荡回路.不难看出,方程(4.1),(4.3)和(4.4)都具有一个明显的特点,就是在这些方程中,未知函数及其导数是一次式,因此这些方程称为线性微分方程.又由于出现在上述方程中的导数的最高阶数为2,故我们称上述方程为二阶线性微分方程.一般的n阶线性微分方程可以写成如下形状:y(n)+p1(x)y(n-1) +…+p n-1(x)y′+p n(x)y=f(x) (4.5) 方程(4.5)的初始条件记为y(x0)=y0,y′(x0)=y0′,…,y (n-1)(x0)=y0 (n-1)(4.6) n阶线性微分方程与第三章讲过的一阶线性微分方程组有着密切的关系,即可以把前者化成后者,而且二者是等价的,这样就可以把前者作为后者的特例加以处理.在方程(4.5)中,令y′=y1,y″=y2,…,y(n-1)=y n-1,(4.5)就可以化成一阶方程组⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧+----====-----)()()()(1111112211x f y x p y x p y x p dxdy y dx dy y dx dy y dx dy n n n n n n(4.7) (4.7)可以写成向量形式()()d x Y x dx=+Y A F (4.8) 其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=--)()()(01000010)(121x p x p p x p x A n n n ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=)(000)(x f x F , ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=-121n y y y y Y 方程组(4.8)的初始条件可记为Y (x 0)=Y 0其中⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=--)1(000010100')()()(n n y y y x y x y x y Y 引理 4.1 方程(4.5)与方程组(4.7)是等价的,即若()y x ϕ=是方程(4.5)在区间I 上的解,则()y x ϕ=,1'()y x ϕ=,…(1)1()n n y x ϕ--=是方程组(4.7)在区间I 上的解;反之,若()y x ϕ=,1'()y x ϕ=,…(1)1()n n y x ϕ--=是方程组(4.7)在区间I 上的解,则()y x ϕ=是方程(4.5)在区间I 上的解.证明 设()y x ϕ=是方程(4.5)在区间I 上的解.令(1)121()'(),()''(),,()()n n x x x x x x ϕϕϕϕϕϕ--=== (4.9)则有⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧+----===---)()()()()()()()()()()(1111)1(211x f x p x x p x x p dx x d x dx x d x dx x d n n n n ϕϕϕϕϕϕϕϕ (4.10)在区间I 上恒成立.这表明, ()y x ϕ=,11()y x ϕ=,…,11()n n y x ϕ--=是方程组(4.7)在区间I 上的解.反之,设()y x ϕ=,11()y x ϕ=,…,11()n n y x ϕ--=是方程组(4.7)在区间I 上的解.于是(4.10)式在区间I 上恒成立.由(4.10)的前n-1个等式.可以看出,函数()x ϕ,1()x ϕ,…,1()n x ϕ-满足关系式(4.9),将它们代入到(4.10)的最后一个等式就有()(1)11()()()()'()()()()n n n n x p x x p x x p x x f x ϕϕϕϕ--++++=在区间I 上恒成立,这就表明()y x ϕ=是方程(4.5)在区间I 上的解.证毕.由引理4.1和第三章的定理3.1′,我们立即可以得到下面的 定理 4.1 如果方程(4.5)的系数p k (x)(k=1,2,…,n)及其右端函数f(x)在区间I 上有定义且连续,则对于I 上的任一x 0及任意给定的y 0,y 0′,…,y 0(n-1),方程(4.5)的满足初始条件(4.6)的解在I 上存在且唯一.在下面的讨论中,总假设(4.5)的系数p k (x)(k=1,2,…,n)及其右端函数f(x)在区间I 上连续,从而,方程(4.5)的满足初始条件(4.6)的解在整个区间I 上总存在且唯一.如果在(4.5)中,f(x)在区间I 上恒等于零,(4.5)变成y (n)+p 1(x)y n-1+…+p n-1(x)y ′+p n (x)y=0 (4.11)方程(4.11)称为n 阶线性齐次微分方程 (或简称n 阶齐次方程),与此相应,(4.5)称为n 阶线性非齐次微分方程(或简称n 阶非齐次方程).有时,为了叙述上的方便,还称(4.11)为(4.5)的对应的齐次方程.4.1.2 齐次线性微分方程的解的性质与结构由引理 4.1,齐次方程(4.11)等价于下面的一阶线性齐次微分方程组Y x A dxdY )(= (4.12) 这里A(x)和Y 与(4.8)中的相同.于是由第三章的定理3.2可知,齐次方程(4.11)的所有解也构成一个线性空间.为了研究这个线性空间的性质,进而搞清楚(4.11)的解的结构,我们需要下面的定义和引理.定义4.1 函数组12(),(),,()n x x x ϕϕϕ称为在区间I 上线性相关,如果存在一组不全为零的常数α1,α2,…,αn ,使得1122()()()0n n x x x αϕαϕαϕ+++= (4.13)在区间I 上恒成立. 反之,如果只当α1=α2=…=αn =0时,才能使(4.13)在I 上成立,则称函数组12(),(),,()n x x x ϕϕϕ在I 上线性无关.引理4.2 一组n-1阶可微的数值函数12(),(),,()n x x x ϕϕϕ在I 上线性相关的充要条件是向量函数组⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-)()(')()1(111x x x n ϕϕϕ , ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-)()(')()1(222x x x n ϕϕϕ ,…, ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-)()(')()1(x x x n n n n ϕϕϕ (4.14)在I 上线性相关.证明 若12(),(),,()n x x x ϕϕϕ在I 上线性相关,则存在一组不全为零的常数α1,α2,…,αn ,使得1122()()()0n n x x x αϕαϕαϕ+++= (4.15)0在I 上恒成立.将(4.15)0式对x 逐次微分n-1次,得'''1122()()()0n n x x x αϕαϕαϕ+++= (4.15)1………………………………………(1)(1)(1)1122()()()0n n n n n x x x αϕαϕαϕ---+++= (4.15)n-1联合(4.15)0,(4.15)1,…,(4.15)n-1就得到向量函数组(4.14)是线性相关的.反之,若向量函数组(4.14)在I 上线性相关,则存在不全为零的常数α1,α2,…,αn ,使得(4.15)0,(4.15)1,…,(4.15)n-1各式在I 上恒成立,由(4.15)0表明12(),(),,()n x x x ϕϕϕ在I 上线性相关.证毕.由引理 4.2,为了建立函数组线性相关与线性无关的判别法则,自然需要引入下面的定义.定义4.2 设函数组12(),(),,()n x x x ϕϕϕ中每一个函数()k x ϕ (k=1,2,…,n)均有n-1阶导数,我们称行列式)()()()(')(')(')()()()()1()1(2)1(12121x x x x x x x x x x W n n n n n n ---=ϕϕϕϕϕϕϕϕϕ为已知函数组的朗斯基(Wronski)行列式.有了以上的准备工作,我们现在可以清楚地看到,齐次方程(4.11)的一般理论完全可以归结为第三章中一阶线性齐次微分方程组的一般理论来加以处理.由§3.3中关于齐次方程组的有关定理,可以自然地得到下面的关于齐次方程(4.11)的一系列定理.定理4.2 齐次方程(4.11)的n 个解12(),(),,()n x x x ϕϕϕ在其定义区间I 上线性无关(相关)的充要条件是在I 上存在点x 0,使得它们的朗斯基行列式W(x 0)≠0(W(x 0)=0).定理4.3 如果12(),(),,()n x x x ϕϕϕ是方程(4.11)的n 个线性无关解,则1122()()()n n y C x C x C x ϕϕϕ=+++ (4.16)是方程(4.11)的通解,其中C 1,C 2,…,C n 为n 个任意常数.通常称定理4.3为方程(4.11)的基本定理.定义4.3 方程(4.11)的定义在区间I 上的n 个线性无关解称为(4.11)的基本解组.由定义4.3,方程(4.11)的基本定理又可叙述为:方程(4.11)的通解为它的基本解组的线性组合.例3 易于验证函数y 1=cos x, y 2=sin x 是方程y ″+y=0的解.并且由它们构成的朗斯基行列式01cos sin sin cos ''2121≠=-=x x xx y y y y在(-∞,+∞)上恒成立,因此,这两个函数是已知方程的两个线性无关解,即是一基本解组,故该方程的通解可写为y(x)=C 1cos x+C 2sin x其中,C 1,C 2是任意常数.不难看出,对于任意的非零常数k 1和k 2函数组y 1=k 1cos x, y 2=k 2sin x都是已知方程的基本解组.基本定理表明,齐次方程(4.11)的所有解的集合是一个n 维线性空间.进一步,我们还有定理4.4 n 阶齐次方程(4.11)的线性无关解的个数不超过n 个. 定理4.5 n 阶齐次方程(4.11)总存在定义在区间I 上的基本解组. 最后,齐次方程(4.11)的解与它的系数之间有如下关系.定理4.6 设12(),(),,()n x x x ϕϕϕ是方程(4.11)的任意n 个解,W(x)是它们的朗斯基行列式,则对区间I 上的任一x 0有10-()0W(x)=W(x )exx p t dt⎰ (4.17)上述关系式称为刘维尔(Liouvill e )公式.由公式(4.17)可以再次看出齐次方程(4.11)的朗斯基行列式的两个重要性质:1.方程(4.11)解的朗斯基行列式W(x)在区间I 上某一点为零,则在整个区间I 上恒等于零.2. 方程(4.11)解的朗斯斯行列式W(x)在区间I 上某一点不等于零,则在整个区间I 上恒不为零.下面给出刘维尔公式的一个简单应用:对于二阶线性齐次方程y ″+p(x)y ′+q(x)y=0如果已知它的一个非零特解y 1,依刘维尔公式(4.17),可用积分的方法求出与y 1线性无关的另一特解,从而可求出它的通解.设y 是已知二阶齐次方程一个解,根据公式(4.17)有⎰=-dsx p C y y y y )(11e '' 或⎰=--ds x p C yy y y )(11e '' 为了积分上面这个一阶线性方程,用211y 乘上式两端,整理后可得 ⎰=⎪⎪⎭⎫ ⎝⎛-dxd p y C y y x )(211e d d 由此可得⎰+⎰=-*21)(1e C dx y C y y dx x p 易见⎰⎰=-dx y y y dxx p )(211e 1是已知方程的另一个解,即C *=0, C=1所对应的解.此外,由0e '')(11≠⎰=-dxx p C y y y y所以,所求得的解y 与已知解y 1是线性无关解.从而,可得已知方程的通解⎰⎰+=-x y Cy C*y y dxx p d e 1)(2111 (4.18) 其中C *和C 是任意常数.例4 求方程(1-x 2)y ″-2xy ′+2y=0的通解.解 容易看出,已知方程有特解y 1=x.此处212)(x xx p --=, 根据公式(4.18),立刻可以求得通解⎪⎭⎫⎝⎛--++=⎪⎭⎫⎝⎛-++-+=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛++-++=⎥⎦⎤⎢⎣⎡-+=⎥⎥⎦⎤⎢⎢⎣⎡⎰+=⎰⎰⎰-111ln 211ln 211112111211)1(e 1**2*22*1221*12x x x C C x x x Cx C dx x x x C C x x x dx C C x dx y C C y y dx x x习题4.1.21. 试讨论下列各函数组在他们的定义区间上是线性相关的还是线性无关的?(1) sin 2,cos ,sin t t t (2) ,tan x x (3) 223,2,24x x x x x -+++ (4) 2,,t t t e te t e2.设在方程()()0y p x y q x y '''++= 中,()p x 在某区间I 上连续且恒不为零,试证它的任意两个线性无关的解的朗斯基行列式是区间I 上的严格单调函数.3.试证明:二阶线性齐次方程的任意两个线性的解组的朗斯基行列式之比是一个不为零的常数..4.已知方程(1)0x y xy y '''--+=的一个解1y x =,是求其通5.已知方程211(1ln )0x y y y x x'''-+-=的一个解1ln y x =,试求其通解. 6.在方程中()()0y p x y q x y '''++=,当系数满足什么条件时,其基本解组的朗斯基行列式等于常数. 7.设1()y x 是n 阶线性齐次方程()(1)11()()()0n n n n y a x y a x y a x y --'+++=的一个非零解. 试证明:利用线性变换1()y y x z =可将已知方程化为1n -阶的齐次方程.8.求方程3260x y x y xy y ''''''-+-=的通解,一只他的两个特解1y x =,22y x =. 9.设111213212223313233()()()()()()()()()()a x a x a x x a x a x a x a x a x a x ∆=, 求证:111213212223313233111213111213212223212223313233313233()()()() ()()()()()()()()()()()() ()()()()()()()()()()()()a x a x a x d x a x a x a x dxa x a x a x a x a x a x a x a x a x a x a x a x a x a x a x a x a x a x a x a x a x '''∆='''++'''并试证对n 阶行列式也有类似结果.4.1.3 非齐次线性微分方程与常数变易法由于n 阶非齐次方程(4.5)等价于一阶非齐次方程组(4.7),于是由第三章的定理3.10,我们有下面的定理4.7 n 阶线性非齐次方程(4.5)的通解等于它的对应齐次方程的通解与它本身的一个特解之和.由此可见,求(4.5)的通解问题,就归结为求(4.5)的一个特解和对应齐次方程的一个基本解组的问题了.和一阶非齐次线性微分方程组一样,对于非齐次方程(4.5),也能够由对应齐次方程的一个基本解组求出它本身的一个特解,即常数变易法.具体作法如下.设y 1,y 2,…,y n 是(4.5)的对应齐次方程的n 个线性无关解,则函数y=C 1y 1+C 2y 2+…+C n y n是(4.5)的通解,其中C 1,C 2,…,C n 是任意常数.现在设一组函数C 1(x),C 2(x),…,C n (x),使n n y x C y x C y x C y )()()(~2211+++= (4.19)成为非齐次方程(4.5)的解.由非齐次方程(4.5)与一阶非齐次方程组(4.7)的等价关系和第三章的(3.18)式,可知,'''12(),(),,()n C x C x C x 满足下面的非齐次方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---)()()()(')(')(')()()()1()1(2)1(12121x y x y x y x y x y x y x y x y x y n n n n n n ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡)(')(')('21x C x C x C n =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡)(00x f (4.20) 它是关于变量'()i C x (i=1,2,…,n)的线性代数方程组,由于它的系数行列式恰是齐次方程的n 个线性无关解的朗斯基行列式,故它恒不为零,因此,上述方程组关于'()i C x 有唯一解.解出后再积分,并代入到(4.19)中,便得到(4.5)的一个特解.例5 求非齐次方程xy y cos 1"=+ 的通解.解 由例3知y 1=cosx,y 2=sinx 是对应齐次方程的线性无关解,故它的通解为y=C 1cosx+C 2sinx现在求已知方程形如y 1=C 1(x)cosx+C 2(x)sinx的一个特解.由关系式(4.20), '1()C x ,'2()C x 满足方程组⎥⎦⎤⎢⎣⎡-cos sin sin cos x x x ⎥⎦⎤⎢⎣⎡)(')('21x C x C =⎥⎥⎦⎤⎢⎢⎣⎡x cos 10或写成纯量方程组⎪⎩⎪⎨⎧=+-=+x x x C x x C x x C x x C cos 1cos )('sin )('0sin )('cos )('2121 解上述方程组,得xxx C cos sin )('1-= , 1)('2=x C 积分得C 1(x)=ln |cosx |, C 2(x)=x故已知方程的通解为y=C 1cos x+C 2sin x+cos x ln |cos x |+x sin x习题4.1.31.设()p x ,()q x ,()f x 在[0,1]上连续,试证明,方程()()()'''++=y p x y q x y f x 满足条件(0)(1)0==y y 的解唯一的充要条件是:方程()()0'''++=y p x y q x y 只有零解满足条件:(0)(1)0==y y .2.用常数变易法求方程21xx e y y e ''-=-的通解.3.在方程32()y y y f x '''++=中,()f x 在[,)a +∞上连续,且lim ()0x f x →+∞=,试证明,已知方程的任一解()y x 均有:lim ()0x y x →+∞=.§4.2 常系数线性微分方程的解法本节先讨论常系数线性齐次方程()(1)11'0n n n n y a y a y a y --++++= (4.21)的求解问题,这里12,,,n a a a 为实常数.由定理4.3,我们知道(4.21)的求解问题归结为求其基本解组即可.虽然对于一般的线性齐次微分方程,人们至今没有找到一个求其基本解组的一般方法,但是对于方程(4.21),这一问题已彻底解决.其中,一个自然的作法是把(4.21)化成与之等价的一阶线性常系数齐次微分方程组,然后按3.5节的有关解法及引理4.1和引理4.2,就可以求得(4.21)的基本解组.但是这样的推导过程并不十分简洁,因此我们这里将对方程(4.21)采用下面的待定指数函数法求解.首先,研究一个简单的一阶方程y ′+ay=0 (4.22)其中a 是常数,不难求出它有特解y=e -ax .比较(4.21)与(4.22),我们可以猜想方程(4.21)也有形如y=e λx (4.23)的解,其中λ是待定常数.将(4.23)代入(4.21)中得到111()0n n x n n a a a e λλλλ--++++= (4.24)因为0x e λ≠,所以有111()0n n n n P a a a λλλλ--=++++=, (4.25)我们称(4.25)为方程(4.21)的特征方程,它的根称为特征根.这样,y =e λx 是方程(4.21)的解,当且仅当λ是特征方程(4.25)的根. 下面分两种情形讨论. 4.2.1 特征根都是单根.定理4.8 若特征方程(4.25)有n 个互异根λ1,λ2,…,λn ,则x n x x n y , y y λλλe ,,e e 2121=== (4.26)是方程(4.21)的一个基本解组.证明 显然,x i iy λe = (i=1,2,…,n)分别是(4.21)的解.它们的朗斯基行列式),,(,0)(e 111e e e ee e e e e e )(1)(1121121)(11211212121212121+∞-∞∈≠-===∏≤<≤++---++---ni j j ixn nn n nx x n n x n x n x n x xxx xx x W n n n n n λλλλλλλλλλλλλλλλλλλλλλλλλλλλλ从而(4.26)是方程(4.21)的一个基本解组.上述行列式为著名的范德蒙(Vandermond)行列式.例1 求方程y ″-5y ′=0的通解.解 特征方程为λ2-5λ=0特征根为λ1=0,λ2=5,故所求通解为y=C 1+C 2e 5x其中C 1,C 2为任意常数.例2 求方程y ″-5y ′+6y =0的通解及满足初始条件:当x =0时,y =1, y ′=2的特解.解 特征方程为2560λλ-+=。