三章阶线性微分方程组三阶线性非齐次方程组的般理论
三阶常系数线性非齐次微分方程通解的降阶法
三阶常系数线性非齐次微分方程通解的降阶法降阶法是高阶线性微分方程的一种解法，它可以解决三阶常系数非齐次微分方程。
下面我们来分析一下它是如何解决三阶常系数非齐次微分方程的。
1. 定义降阶法降阶法是一种用于解决三阶常系数非齐次微分方程的算法，它将三阶常系数非齐次微分方程转化为一组互相关的线性一阶方程组。
2. 三阶常系数非齐次微分方程三阶常系数非齐次微分方程是在三阶线性常系数微分方程的基础上，涉及右端非齐次项，则称为三阶常系数非齐次微分方程，它的一般形式为：$$y^{'''}+a_2y''+a_1y'+a_0y=g(x)$$3. 降阶法的基本思想降阶法的基本思想是将三阶general equation降低到一组互相关的线性一阶方程组，通过求解这个方程组来解决三阶general equation，换言之，就是将三阶微分方程转化为三个一阶微分方程。
4. 降阶法的具体步骤（1）令$u_1=y$ 、$u_2=y'$和$u_3=y''$ ，引入三个新变量。
（2）将三阶常系数非齐次微分方程变换为三个一阶微分方程：$\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad u_1'=u_2$$\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad u_2'=u_3$$\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad u_3'=g(x)-a_2u_2-a_1u_1-a_0u_0$（3）解上述方程组，即可求得三阶常系数非齐次微分方程原方程的通解。
5. 降阶法的优缺点（1）优点：相比于其他解法，降阶法计算量较小，易于推导和实现。
（2）缺点：当微分方程非常复杂时，降阶法可能会出现运算失真或者不稳定的现象，影响最终结果的准确性。
三阶常系数线性非齐次微分方程的常数变易解法
若 方程（ ７ ） 对应 的齐次方程 有复特征根 知， 方程（ ７ ） 有特 解
， 则 由文 献【 ４ 】 可
＝ ｅ － ， （ ） ．
是方程（ ４ ） 化成 Ｃ （ ） ＝ ： ｅ ， （ ） ，
（ ４ ）
（ ２ ５ ＋ 】 ｍｌ ｈ
（ｉ）若实根 ａ是三重 特征根 ， 则Ｐ （ ） － Ｐ （ ￣ ） － － ｏ。 于
Ｄ （ ）
ｉ ｒ 恤ｌ 【 — Ｆ 一・
ｓ ｉ ｎ｜ Ｅ
ｆ 。
于是
） ｓ ｉ
，
，
， ． （ 撕 ３ 口 ｋ
（ ９ ）
积 分三 次， 即 得到ｃ （ ） ＝ Ｉ ｛ ｌ 【 Ｉ ｅ － － ｙ （ ） 出 】 ｌ ． 从而 方程
（ １ ） 的特解 为
三阶常 系数 线性非齐次微分方程 的常数变 易解法
胡 爱 莲
（ 喀什师范学院 数学系 ， 新疆 喀什 ８ ４ ４ ０ ０ ８ ）
摘 要： 将常数变 易法 应用 于三阶常系数线性非齐次微分方 程 ， 对一般 非齐次 自由项形式 ， 给 出了方 程的特解 公式 ， 进 而求得 了通解 ． 关键词 ： 常数变易法 ； 三阶常 系数线性非 齐次微分方程 ； 通解 ； 特解
ｙ ＝ ｅ ａ Ｉ ｆ ｛ ｆ ［ ｅ －
，
ｆ ‘
一
） 出 】 也 ｝ 出 ．
（ ８ ）
数变易法， 可设方程（ １ ） 有特解 ｙ ＝ ｃ （ ） ， 代人（ １ ） 得：
Ｃ 肿 （ ） ＋ （ ３ ａ ＋ 口 １ ）Ｃ （ ） ＋ （ ３ ０ ｌ ａ ＋ ａ ａ ） Ｃ （ ）
系数 法求得 （ １ ） 的特 解． 但 当方 程（ Ｉ ） 的非 齐次项是 其他 函
常微分方程复习提要全文
式
dyi (x) dx
fi (x, y1(x),
, yn (x)), (i 1.2
n)
则称 y1(x), , yn (x) 为微分方程组(3.1)在区间 [a,b] 的一个解。
通解及通积分：
含有n个任意常数 c1, cn 的方程组（3.1）的解
y1 1(x, c1, cn )
yn
n (x, c1,
齐次方程组的解组线性相关性的判别法：
推论3.3 方程组(3.8)的n个解在其定义区间I上线性 无关的充要条件是它们的朗斯基行列式W(x)在I上任一点
不为零.
解组
线性相关 W ( x0 )=0 线性无关 W ( x0 ) 0
我们把一阶线性齐次方程组(3.8)的n个线性无关解 称为它的基本解组。其对应的矩阵称为基本解矩阵。
（其中F为已知的函数）
定义（P3） ：微分方程中出现的未知函数的 最高阶导数的阶数(或微分的阶数)称为微分方程的 阶数.
定义（P4） ：如果一个微分方程关于未知函数 及其各阶导数都是一次的，则称它为线性微分方程， 否则称之为非线性微分方程.
定义（P4）: 设函数 y x在区间I上连续，且有
dy1
dx
a11( x) y1
a12 ( x) y2
dy2 dx
a21( x) y1
a22 ( x) y2
dyn dx
an1( x) y1
an2 ( x) y2
a1n ( x) yn f1( x),
a2n ( x) yn f2 ( x), （3.6）
ann ( x) yn fn ( x).
解法：两边除以yn ，得 yn dy p( x) y1n f ( x) dx
令z y1n ，则 dz (1 n) yn dy ,代入方程
常微分方程第五章微分方程组总结
一．线性微分方程组的一般理论1. 线性微分方程组一般形式为：1111122112211222221122()()()(),()()()(), 1 ,()()()(),n n n n nn n nn n n x a t x a t x a t x f t x a t x a t x a t x f t x a t x a t x a t x f t '=++++⎧⎪'=++++⎪⎨⋅⋅⎪⎪'=++++⎩（） 记：111212122212111222()()()()()()()()()()()()(), , ()n n n n nn n n n a t a t a t a t a t a t A t a t a t a t f t x x f t x x f t x x f t x x ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦'⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥'⎢⎥⎢⎥⎢⎥'===⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥'⎣⎦⎣⎦⎣⎦非齐次线性方程组表示为：()() x A t x f t '=+齐次线性方程组表示为：()x A t x '=2.齐次线性方程组的一般理论(1)定理 (叠加原理) 如果12(),(),,()n x t x t x t ⋯是齐次方程组()x A t x '=的k 个解，则它们的线性组合1212()()()n n c x t c x t c x t ++⋯+也是齐次方程组的解，这里12,,,n c c c ⋯是任意常数(2)向量函数线性相关性定义在区间],[b a 上的函数12(),(),,()n x t x t x t ⋯，如果存在不全为零的常数k c c c ,,,21⋯使得1212()()()0n n c x t c x t c x t ++⋯+≡在],[b a 上恒成立，我们称这些向量函数是线性相关的，否则称这些向量函数线性无关。
三阶非齐次常系数微分方程解的表达式
教学单位学生学号本科毕业论文（设计）题目学生姓名专业名称指导教师年月日三阶非齐次常系数线性微分方程解的表达式胡青（111114110）（湖北工程学院数学与统计学院湖北孝感 432000）摘要：利用常数变易法求解二阶常系数非齐次线性微分方程，这种方法被许多作者研究过。
本文利用常数变易法求出了三阶非齐次常系数微分方程解的表达式，利用解的表达式，可以很方便地求出三阶非齐次常系数微分方程的解。
关键词：常数变易法；三阶；非齐次；常系数微分方程；解的表达式The Solution of third-order Non-homogeneous Ordinary Differential Equation with Constant C oefficientHuQing (111114110)(Hubei Engineering College, Xiaogan, Hubei, 432000)Abstract:Uses the method of constant variation to second-order constant coefficient non-homogeneous linear differential equation, this method by many authors studied. In this paper, the constant variation method is used to derive the third order nonhomogeneous the expressions of the constant coefficient differential equations, using the expression of the solution can be easily calculated third-order differential equation with constant .Key words:constant variation method; The third order. Non-homogeneous; Constant coefficient differential equations; Expression of solution1. 预备知识定义1.1：方程0...)(111=++++≡--n n n n a a a F λλλλ称为常系数齐次线性方程0...)1(1)2(2)1(1)(=+++++---x a y a y a y a y n n n n n 的特征方程。
线性微分方程的基本理论
( k 1) x1 (t )
x2 ( t ) x2 ( t )
( k 1) x2 (t )
xk ( t ) xk ( t )
( k 1) xk (t )
t 0 t 0
显然对所有的 t, 恒有 W [ x1 ( t ), x2 ( t )] 0,
但 x1 ( t ), x2 ( t ) ( , ) 上线性无关. 在 事实上, 假设存在恒等式
c1 x1 (t ) c2 x2 (t ) 0,
则当 t 0 时, 有c2 0, 当 t 0 时, 有c1 0, 故x1 ( t ), x2 ( t ) ( , )上线性无关. 在
n n
a1 ( t )
d
n1
x
dt
n1
an1 ( t )
dx dt
an ( t ) x
L[e
t
n
]?
n1
[ a1 (t )
a2 (t )
n 2
an1 (t ) an (t )]e .
t
性质3.1 L(cx ) cL( x ) c 为常数. 性质3.2
或
在(a, b) 上有定义,
则在（a, b）上线性无关的充要条件为
x1 ( t ) x2 ( t )
或
x2 ( t ) x1 ( t )
在（a, b）上不恒为常数.
8
例3: sin t , cos t在任何区间上都线性无关. 2 2 cos t ,1 sin t 在任何区间上都线性相关. 注3：函数组的线性相关与线性无关是 依赖于所取的区间。 例4: 函数x1 ( t ) t , x2 ( t ) t 在（ ， ） 上 是线性无关, 而在 (0, )和 ( , 0) 上是线性相关的.
4.1--线性微分方程的一般理论
(1) ( n1) ( n1) ( t0 ) x0 , ( t0 ) x0 , , ( t0 ) x0 .
二、齐次线性方程的解的结构与性质
dn x d n 1 x dx an ( t ) x 0 n a1 ( t ) n 1 a n 1 ( t ) dt dt dt ( 2)
定理5 n 阶齐次线性方程(2)必存在n个线性无关的解.
定理6 (通解结构定理) 若 x1(t), x2(t),, xn(t)是方程(2)
的n个线性无关的解, 则方程(2)的通解可表为： x = c1x1(t)+ c2x2(t)+ +cnxn (t) (3) (其中 c1, c2, ,cn 是任意常数.)
2t
定理 设非齐次方程(1)的右端 f(t) 是几个函数之和,
(n 阶非齐次线性方程解的叠加原理)
作业
P131: 1, 3.(1),(2), 6,7
*例 已知 y = x 及 y = sinx 为某二阶齐线性方程的 解, 求该方程 .
解 设方程为 y P ( x ) y Q( x ) y 0 , y x 为其解 , 有 P ( x ) xQ( x ) 0 ,
例1 验证x1=cos t, x2=sin t 为方程 x+x = 0 的基本 解组, 并求其通解. 解 ∵ x1+x1 = cos t +cos t =0, x1 为方程的解.
∵ x2+x2 = sin t +sin t =0, x2 为方程的解.
由它们构成的朗斯基行列式
故 x1, x2 是方程的一个基本解组,
3第三章 线性微分方程
第三章 线性微分方程线性微分方程是常微分方程理论的重要组成部分，它在自然科学和工程技术方面有着极其广泛的应用，很多实际问题都可以用线性微分方程来处理.在第一章中，我们学习了一阶线性微分方程的处理方法，本章介绍高阶线性微分方程的概念以及常系数线性微分方程的解法.3.1 线性微分方程的一般理论在第一章，我们介绍了一阶线性微分方程)()(x q y x p y =+'的解法，这里介绍n 阶线性微分方程的概念，以及解的存在唯一性定理.首先来研究下面一个实际的例子.例1 弹簧振动设一质量为m 的物体A 悬挂在一上端固定的弹簧的末端（假设弹簧的质量相对于物体A 的质量可以忽略）试求该物体在外力)(t f 作用下的所满足的微分方程.当物体A 不受外力时，重力与弹簧的拉力平衡时的位置选为坐标轴x 的原点O ，向下的方向取为x 轴的正向.设t 时刻，物体A 的位移为)(t x ，速度为)(t v ，加速度为)(t a ，则22)(,)(dtx d t a dt dx t v ==. 由牛顿第二定律知，ma F =，其中m 是物体A 的质量，a 是加速度，F 是合外力.下面对物体A 所受到的力进行分析，由以下几个部分构成.（1）弹簧的拉力1f ，设弹簧的弹性系数为k ，在t 时刻，物体的位移为)(t x ，依据胡克定律知kx f -=1（2）空气的阻力2f ，当速度不太大时，空气的阻力与物体的速度成正比，设比例常数为)0(>μ，因为阻力的方向与速度的方向相反，所以dtdx v f μμ-=-=2 （3）外力)(t f因此，我们得到合外力)(t f dtdx kx F +--=μ 代入ma F =得物体A 所满足的微分方程为)(22t f kx dt dx dtx d m =++μ （3.1） 那么物体A 的运动规律方程)(t x 就是上述微分方程（3.1）的解，如何解该方程是我们本章学习的重点.方程（3.1）和我们第一章中所学的一阶线性微分方程一样，都是线性微分方程，因为方程（3.1）中的导数的最高阶数为2，所以（3.1）是二阶线性微分方程.一般的n 阶线性微分方程具有如下的形式：)0)((),()()()()(01)1(1)(0≠=+'+++--x a x y x a y x a y x a y x a n n n n ϕ因为0)(0≠x a ，所以上式可化为)()()()(1)1(1)(x f y x p y x p y x p y n n n n =+'+++-- （3.2） 其中)()()();,,2,1(,)()()(00x a x x f n i x a x a x p i i ϕ=== 方程（3.2）的初值条件为)1(00)1(0000)(,,)(,)(--='='=n n y x y y x y y x y （3.3） 方程（3.2）在什么条件下存在满足初值条件（3.3）的解呢？有解的话，其解是否唯一？存在区间又是什么呢？为了解决这些问题，我们先给出一般的n 阶方程),,,,()1()(-'=n n y y y x f y满足初值条件)1(00)1(0000)(,,)(,)(--='='=n n y x y y x y y x y 的解的存在唯一性定理.定理3.1 如果函数),,,,()1(-'n y y y x f 在闭区域1)1(0)1(100001,,,,:---+≤-≤'-'≤-≤-n n n n b y y b y y b y y a x x R上满足（1）),,,,()1(-'n y y y x f 在1+n R 上连续；（2）),,,,()1(-'n y y y x f 在1+n R 上关于变量)1(,,,-'n y y y 满足李普希兹条件，即存在正数N ，使得对于任何一对点1)1(222)1(111),,,,(),,,,,(+--∈''n n n R y y y x y y y x ，总有)(),,,,(),,,,()1(2)1(12121)1(222)1(111-----++'-'+-≤'-'n n n n y y y y y y N y y y x f y y y x f 则初值问题⎩⎨⎧='='='=---)1(00)1(0000)1()()(,,)(,)(),,,,(n n n n y x y y x y y x y y y y x f y 在0000h x x h x +≤≤-上存在唯一解)(x y ϕ=. 这里),,,,(max },,,,,min{)1(),,,,(11001)1(-∈'-'==+-n R y y y x n y y y x f M Mb M b M b a h n n . 定理3.1的证明和第二章中的解的存在唯一性定理的证明是相仿的，读者可以模仿定理2.1的证明，完成定理3.1的证明.n 阶线性微分方程（3.2）只是一般的n 解微分方程的一种特殊形式，和一阶线性微分方程类似，有如下的解的存在唯一性定理.定理3.2 如果方程（3.2）的系数),,2,1)((n i x p i =以及右端函数)(x f 在区间],[b a 上有定义而且都连续，则初值问题⎩⎨⎧='='==+'+++----)1(00)1(00001)1(1)()(,,)(,)()()()()(n n n n n n y x yy x y y x y x f y x p y x p y x p y 在],[b a 上存在唯一解)(x y ϕ=.该定理的证明可以利用定理3.1，只要),,2,1)((n i x p i =以及)(x f 在区间],[b a 上连续，则该初值问题就满足了定理3.1的两个条件，从而定理3.2是成立的.如无特别声明，在本章的讨论中，总假定方程（3.2）的系数),,2,1)((n i x p i =以及右端函数)(x f 在区间],[b a 上有定义而且都连续，从而，方程（3.2）满足初值条件（3.3）的解在闭区间],[b a 上存在且唯一.特别地，如果0)(=x f ，则方程（3.2）变为0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n （3.4）方程（3.4）称为n 阶线性齐次微分方程；如果0)(≠x f ，则称方程（3.2）为n 阶线性非齐次微分方程.这时，称方程（3.4）为方程（3.2）所对应的n 阶线性齐次微分方程.3.2 n 阶线性齐次微分方程的一般理论方程（3.4）称为n 阶线性齐次微分方程，对于这类方程应该如何处理呢？我们先研究方程（3.4）的解的性质（1）如果1y 是方程（3.4）的解，则对任意常数C ，1Cy 也是方程（3.4）的解.证明 因为1y 是方程（3.4）的解，所以0)()()(111)1(11)(1=+'+++--y x p y x p y x p y n n n n 从而有))(())(())(()(111)1(11)(1Cy x p Cy x p Cy x p Cy n n n n +'+++--00))()()((111)1(11)(1=⨯=+'+++=--C y x p y x p y x p y C n n n n 因此，1Cy 也是方程（3.4）的解.（2）如果21,y y 是方程（3.4）的解，则21y y +也是方程（3.4）的解.证明 因为21,y y 是方程（3.4）的解，所以0)()()(111)1(11)(1=+'+++--y x p y x p y x p y n n n n 0)()()(221)1(21)(2=+'+++--y x p y x p y x p y n n n n 从而有))(())(())(()(21211)1(211)(21y y x p y y x p y y x p y y n n n n ++'++++++--))(())((])())[(()()(21211)1(2)1(11)(2)(1y y x p y y x p y y x p y y n n n n n n ++'+'+++++=--- ))()()(())()()((221)1(21)(2111)1(11)(1y x p y x p y x p y y x p y x p y x p y n n n n n n n n +'+++++'+++=---- 000=+=因此，21y y +也是方程（3.4）的解.推论 如果n y y y ,,,21 是方程（3.4）的解，则对任意n 个常数n C C C ,,,21 ，线性组合n n y C y C y C +++ 2211也是方程（3.4）的解.该推论可由性质（1）和性质（2）直接推出.并且根据性质（1）和性质（2），我们可以得出，n 阶线性齐次微分方程（3.4）的解构成一个线性空间，称为解的线性空间.例1 易于验证函数x x e x y e x y -==)(,)(21是方程0=-''y y的解，因此，函数x x e C e C x y -+=21)(也是原方程的解.反过来，方程0=-''y y 的通解是不是x x e C e C x y -+=21)(呢？同样地，给出了方程（3.4）的n 个解)(,),(),(21x y x y x y n 后，含有n 个任意常数n C C C ,,,21 的函数)()()()(2211x y C x y C x y C x y n n +++=是否就是线性齐次微分方程的通解呢？为了解决这个问题，我们首先给出函数组的线性相关和线性无关的概念.定义3.1 函数组)(,),(),(21x y x y x y n 在区间I 上有定义，如果其中的某个函数可由其余的1-n 个函数线性表出，则称函数组)(,),(),(21x y x y x y n 在区间I 上是线性相关的.如果任何一个函数都不能由其余的1-n 个函数线性表出，则称函数组)(,),(),(21x y x y x y n 在区间I 上是线性无关的.定理3.3（判定定理）函数组)(,),(),(21x y x y x y n 在区间I 上有定义，如果存在一组不全为零的常数n k k k ,,,21 ，使得0)()()(2211=+++x y k x y k x y k n n则函数组)(,),(),(21x y x y x y n 在区间I 上线性相关.如果只有n k k k ,,,21 全为零时，0)()()(2211=+++x y k x y k x y k n n 才成立，则称函数组)(,),(),(21x y x y x y n 在区间I 上线性无关.证明 如果存在一组不全为零的常数n k k k ,,,21 ，使得0)()()(2211=+++x y k x y k x y k n n不妨假定01≠k ，则)()()(12121x y k k x y k k x y n n ---= 即)(1x y 可由)(,),(2x y x y n 线性表出，因此函数组)(,),(),(21x y x y x y n 在区间I 上线性相关.由线性相关与线性无关的定义，我们可以很容易的得出下面的结论：（1）在函数组)(,),(),(21x y x y x y n 中，如果含有一个零函数，比如0)(=x y i ，则函数组)(,),(),(21x y x y x y n 线性相关.事实上，)(x y i 可由其余1-n 个函数线性表出，即)(0)(0)(0)(0)(00)(1121x y x y x y x y x y x y n i i i ⨯++⨯+⨯++⨯+⨯=+-（2）如果函数组只含有两个函数)(),(21x y x y ，则它们线性相关等价于它们之比)()(21x y x y 为常数. 证明 如果)(),(21x y x y 线性相关，则其中的一个可以由另一个线性表出，不妨设)()(21x y x y →，则存在常数k ，使得)()(21x ky x y =，所以k x y x y =)()(21为常数. 反过来，如果)()(21x y x y 为常数，设k x y x y =)()(21，则)()(21x ky x y =， 即)()(21x y x y →，所以)(),(21x y x y 线性相关.例2 函数组x x e x y e x y -==)(,)(21在任意区间是线性无关的.证明 因为x e x y x y 221)()(=不是常数，所以x x e x y e x y -==)(,)(21线性无关. 例3 函数组x x y x x y x x y 23221sin )(,cos )(,2cos )(===在任意区间上是线性相关的.证明 因为x x x 22sin cos 2cos -=，即)()()(321x y x y x y -=，所以函数组 x x y x x y x x y 23221sin )(,cos )(,2cos )(===线性相关.对于一般的函数组，直接用线性相关和线性无关的定义，或者用定理3.3去判定是非常困难的.比如函数组1)(,cos )(,sin )(321===x y x x y x x y 是线性相关的，还是线性无关的呢？为了解决这个问题，我们下面给出朗斯基（Wronski ）行列式的定义.定义3.2 设函数组)(,),(),(21x y x y x y n 在区间I 上有定义，而且都存在1-n 阶导数，我们称行列式)()()()()()()()()()()1()1(2)1(12121x y x y x y x y x y x y x y x y x y x W n n n n n n ---'''=为函数组)(,),(),(21x y x y x y n 的朗斯基行列式.函数组的朗斯基行列式有如下的性质：定理3.4 函数组)(,),(),(21x y x y x y n 在区间I 上有定义而且线性相关，又对每个),,2,1)((n k x y k =存在1-n 阶导数，则它们的朗斯基行列式恒等于零.证明 因为函数组)(,),(),(21x y x y x y n 线性相关，所以其中的一个函数可由其余1-n 个函数线性表出，不妨设)(x y i 可由其余1-n 个函数线性表出，即存在1-n 个常数n i i k k k k k ,,,,,,1121 +-，使得)()()()()()(11112211x y k x y k x y k x y k x y k x y n n i i i i i ++++++=++--而且有)()()()()()(11112211x y k x y k x y k x y k x y k x y n n i i i i i '++'+'++'+'='++-- …………………………………………………………………………………………)()()()()()()1()1(11)1(11)1(22)1(11)1(x y k x y k x y k x y k x y k x y n nn n i i n i i n n n i --++------++++++= 则在朗斯基行列式)()()()()()()()()()()()()()1()1()1(2)1(12121x y x y x y x y x y x y x y x y x y x y x y x y x W n n n n n i n n i i ----''''=的计算中，第一列乘以1k -，第二列乘以2k -，…，第1-i 列乘以1--i k ，第1+i 列乘以1+-i k ，…，第n 列乘以n k -后，全加至第i 列，则第i 列中的每个元素全为零，所以0)()()(0)()(0)()(0)()()()1()1(2)1(12121='''=---x y x y x y x y x y x y x y x y x y x W n nn n n n推论1 函数组)(,),(),(21x y x y x y n 在区间I 上有定义，而且都存在1-n 阶导数，如果存在某一点I x ∈0，使得0)(0≠x W ，则函数组)(,),(),(21x y x y x y n 线性无关.该推论是定理3.4的逆否命题，所以显然是成立的.例4 函数组1)(,cos )(,sin )(321===x y x x y x x y 在任意区间上是线性无关的.证明 函数组1)(,cos )(,sin )(321===x y x x y x x y 的朗斯基行列式为01cos sin sin cos 0cos sin 0sin cos 1cos sin )(≠-=---=---=x x x x x x x x x x x W 所以函数组1)(,cos )(,sin )(321===x y x x y x x y 是线性无关的.例5 函数组x n x x n e x y e x y e x y λλλ===)(,,)(,)(2121 （其中n λλλ ,,21两两互异）在任意区间上是线性无关的.证明 函数组x n x x n e x y e x y e x y λλλ===)(,,)(,)(2121 的朗斯基行列式为1121121)(1121121111)(21212121---+++---==nnn n nx x nn x n x n x n x x x x xn n n n e e e e e e e e e e x W λλλλλλλλλλλλλλλλλλλλλλλλ∏≤<≤+++-=n j i i j xn e 1)()(21λλλλλ因为n λλλ ,,21两两互异，所以0)()(1)(21≠-=∏≤<≤+++n j i i j xn e x W λλλλλ ，因此，函数组x n x x n e x y e x y e x y λλλ===)(,,)(,)(2121 线性无关.由定理3.4知，函数组)(,),(),(21x y x y x y n 线性相关，则它们的朗斯基行列式恒为零；反过来，如果函数组)(,),(),(21x y x y x y n 的朗斯基行列式恒为零，能不能得出函数组)(,),(),(21x y x y x y n 线性相关呢？即朗斯基行列式恒为零是不是线性相关的充分必要条件呢？下面的例子给出了答案.例6 函数组,0,0)(21<≥⎩⎨⎧=x x x x y 0,0,0)(,22<≥⎩⎨⎧=x x x x y 显然是线性无关的，因为)()(21x y x y 不是常数.但是，它们的朗斯基行列式为 当0≥x 时，0020)(2==x x x W ，当0<x 时，0200)(2==xx x W . 即，对所有的R x ∈，朗斯基行列式恒为零.通过这个例子，我们可以看出即使函数组的朗斯基行列式恒为零，该函数组也有可能是线性无关的，所以函数组的朗斯基行列式恒为零是判定函数组线性相关的必要条件，而不是充分条件.我们只能用朗斯基行列式在某点处不为零，判定该函数组线性无关.但是，如果函数组)(,),(),(21x y x y x y n 是方程0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n的n 个解，这时它们的朗斯基行列式恒为零是判定该函数组线性相关的充分必要条件.这可由下面的定理得到.定理3.5 如果函数组)(,),(),(21x y x y x y n 是方程0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n定义在区间),(b a 上的n 个线性无关的解，则它们的朗斯基行列式0)(≠x W在区间),(b a 上恒成立.证明 假设0)(≠x W 在区间),(b a 上不恒成立，即存在),(0b a x ∈，使得0)(0=x W .构造n C C C ,,,21 的方程组⎪⎪⎩⎪⎪⎨⎧=+++='++'+'=+++---0)()()(0)()()(0)()()(0)1(0)1(220)1(1100220110022011x y C x y C x y C x y C x y C x y C x y C x y C x y C n n n n n n n n n因为方程组的系数行列式0)(0=x W ，所以方程组有非零解，设为)0()0(2)0(1,,,n C C C （)0()0(2)0(1,,,nC C C 不全为零） 以)0()0(2)0(1,,,nC C C 为系数，构造)(,),(),(21x y x y x y n 的线性组合)()()()()0(2)0(21)0(1x y C x y C x y C x y n n +++=根据齐次线性微分方程解的性质知，它是方程（3.4）的解.而且满足初始条件⎪⎪⎩⎪⎪⎨⎧=+++=='++'+'='=+++=----0)()()()(0)()()()(0)()()()(0)1()0(0)1(2)0(20)1(1)0(10)1(0)0(02)0(201)0(100)0(02)0(201)0(10x y C x y C x y C x yx y C x y C x y C x y x y C x y C x y C x y n n n n n n n n n n 而0)(≡x y 也是方程（3.4）的解，也满足0)(,,0)(,0)(0)1(00=='=-x y x y x y n .因为初值问题⎩⎨⎧=='==+'+++---0)(,,0)(,0)(0)()()(0)1(001)1(1)(x yx y x y y x p y x p y x p y n n n n n 的解是唯一的，所以0)()()()()0(2)0(21)0(1≡+++=x y C x y C x y C x y n n又)0()0(2)0(1,,,nC C C 不全为零，由定理3.3知，函数组)(,),(),(21x y x y x y n 是线性相关，这与函数组)(,),(),(21x y x y x y n 线性无关矛盾. 所以，假设不成立，因此，0)(≠x W 在区间),(b a 上恒成立. 推论2如果函数组)(,),(),(21x y x y x y n 是方程0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n定义在区间),(b a 上的n 个解，如果存在),(0b a x ∈，使得它们的朗斯基行列式0)(0=x W则该解组在),(b a 上线性相关.该推论是定理3.5的逆否命题，所以显然是成立的. 推论3方程0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n的n 个解)(,),(),(21x y x y x y n 在其定义区间),(b a 上线性无关的充要条件是，存在),(0b a x ∈，使得它们的朗斯基行列式0)(0≠x W . 这样，我们就可以得出下面的结论：设函数组)(,),(),(21x y x y x y n 是方程0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n在区间),(b a 上的n 个解，则)(,),(),(21x y x y x y n 线性相关⇔)),((0)(0)(00b a x x W x W ∈=⇔≡ )(,),(),(21x y x y x y n 线性无关⇔)),((0)(0)(00b a x x W x W ∈≠⇔≠.这样，我们判定方程的n 个解是线性相关还是线性无关就可以看),,(0b a x ∈∀ )(0x W 是否为零.解决了线性相关与线性无关的判定问题后，我们继续解决线性齐次微分方程的通解问题，首先给出基本解组的概念.定义3.3 方程（3.4）的定义在区间),(b a 上的n 个线性无关的解，称为方程（3.4）的基本解组.例7 在例1中，我们验证了函数x x e x y e x y -==)(,)(21是方程0=-''y y的解，而且x x e x y e x y -==)(,)(21是线性无关的两个解，即)(),(21x y x y 是方程0=-''y y的基本解组，那么x x e C e C x y -+=21)(是否为方程的通解呢？我们需要证明，方程的任一解是否可由基本解组线性表出？定理 3.6 如果函数组)(,),(),(21x y x y x y n 是方程（3.4）的一个基本解组，则对于（3.4）的任一解)(x y ，均可由函数组)(,),(),(21x y x y x y n 线性表出，即，存在一组数)0()0(2)0(1,,,nC C C ，使得 )()()()()0(2)0(21)0(1x y C x y C x y C x y n n +++= .证明 设)(x y 是方程（3.4）的任一解，并且满足初始条件)1(00)1(0000)(,,)(,)(--='='=n n y x y y x y y x y 构造n C C C ,,,21 的方程组⎪⎪⎩⎪⎪⎨⎧=+++'='++'+'=+++----)1(00)1(0)1(220)1(110002201100022011)()()()()()()()()(n n n n n n n n n n y x y C x y C x y C y x y C x y C x y C y x y C x y C x y C因为)(,),(),(21x y x y x y n 是方程（3.4）的基本解组，即)(,),(),(21x y x y x y n 是线性无关的，所以它们的朗斯基行列式在0x x =的值不为零，即方程组的系数行列式0)(0≠x W ，因此方程组存在唯一的解，设其解为)0()0(2)0(1,,,n C C C ，这样，我们可以用)0()0(2)0(1,,,nC C C 构造函数 )()()()()0(2)0(21)0(1~x y C x y C x y C x y n n +++=根据线性齐次微分方程解的性质知，)(~x y 是方程（3.4）的解，而且)(~x y 满足：⎪⎪⎩⎪⎪⎨⎧=+++='='++'+'='=+++=-----)1(00)1()0(0)1(2)0(20)1(1)0(10~)1(00)0(02)0(201)0(10~00)0(02)0(201)0(10~)()()()()()()()()()()()(n n n n n n n n n n n y x y C x y C x y C x yy x y C x y C x y C x y y x y C x y C x y C x y 又)(x y 是方程（3.4）的任一解，并且满足初始条件)1(00)1(0000)(,,)(,)(--='='=n n y x y y x y y x y 因为初值问题⎩⎨⎧='='==+'+++----)1(00)1(00001)1(1)()(,,)(,)(0)()()(n n n n n n y x y y x y y x y y x p y x p y x p y 的解是唯一的，所以)()()()()()0(2)0(21)0(1~x y C x y C x y C x y x y n n +++==即，)(x y 可由函数组)(,),(),(21x y x y x y n 线性表出. 由定理3.6，可以得出下面的基本定理.定理3.7（基本定理）如果函数组)(,),(),(21x y x y x y n 是方程（3.4）的一个基本解组，则)()()()(2211x y C x y C x y C x y n n +++=是方程（3.4）的通解，其中n C C C ,,,21 是n 个任意常数.证明 首先由线性齐次微分方程解的性质知，对任意的n C C C ,,,21)()()()(2211x y C x y C x y C x y n n +++=是方程（3.4）的解.其次，由定理3.6知，方程（3.4）的任一解均可由)(,),(),(21x y x y x y n 线性表出，即，任一解都可以表示成)()()()(2211x y C x y C x y C x y n n +++= 的形式. 因此，)()()()(2211x y C x y C x y C x y n n +++= 是方程（3.4）的通解. 由基本定理可知，方程（3.4）的求解，关键是找到方程（3.4）的一个基本解组，即，方程（3.4）的n 个线性无关的解)(,),(),(21x y x y x y n ，这样就可以很容易的写出方程（3.4）的通解：)()()()(2211x y C x y C x y C x y n n +++= .例8 求方程0=-''y y的通解.解 我们验证了函数x x e x y e x y -==)(,)(21是方程0=-''y y的两个线性无关的两个解，因此方程的通解为：x x e C e C x y -+=21)(.根据基本定理，我们可以得到下面的推论.推论4 n 阶线性齐次微分方程（3.4）的线性无关解的个数不超过n 个. 证明 设)(),(,),(),(121x y x y x y x y n n + 是方程（3.4）的任意1+n 个解.如果前n个解)(,),(),(21x y x y x y n 是线性相关的，则上述1+n 个解是线性相关的. 如果前n 个解)(,),(),(21x y x y x y n 是线性无关的，由定理3.6知，)(,),(),()(211x y x y x y x y n n →+所以)(),(,),(),(121x y x y x y x y n n + 是线性相关的. 因此，方程（3.4）的线性无关解的个数不超过n 个.由推论4知，n 阶线性齐次微分方程（3.4）的线性无关解的个数不超过n 个，那么方程（3.4）的线性无关解的个数是不是有且只有n 个呢？其基本解组存在吗？下面的定理可以回答这个问题.定理3.8 方程（3.4）总存在定义在),(b a 上的基本解组，即，总存在n 个线性无关的解.证明 在),(b a 上任取一点0x x =，由解的存在唯一性定理，方程（3.4）在),(b a 上必存在n 个解)(,),(),(),(321x y x y x y x y n ，它们分别满足下列初始条件：0)(,,0)(,0)(,1)(0)1(1010101==''='=-x y x y x y x y n 0)(,,0)(,1)(,0)(0)1(2020202==''='=-x y x y x y x y n 0)(,,1)(,0)(,0)(0)1(3030303==''='=-x y x y x y x y n………………………………………………………1)(,,0)(,0)(,0)(0)1(000==''='=-x y x y x y x y n n n n n由于这n 个解在0x x =点的朗斯基行列式的值011000010100001)(0≠==x W所以)(,),(),(),(321x y x y x y x y n 是线性无关的解，从而它们是方程（3.4）定义在),(b a 上的基本解组.由定理3.8知，线性齐次微分方程（3.4）的基本解组一定存在，且含有n 个线性无关的解.而且方程（3.4）的解与它的系数之间满足如下的刘维尔（Liouville ）公式.定理3.9 设)(,),(),(21x y x y x y n 是方程0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n的任意n 个解，)(x W 是这n 个解的朗斯基行列式，则对),(b a 上的任一点0x x =，总有⎰=-xx dtt p ex W x W 01)(0)()(.在给出定理3.9的证明之前，先给出行列式函数求导法则 设n 阶行列式函数为)()()()()()()()()()(212222111211x a x a x a x a x a x a x a x a x a x D nn n n n n n =则+'''+'''=)()()()()()()()()()()()()()()()()()())((212222111211212222111211x 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 dxx D d nn n n n n nn n n n n n)()()()()()()()()(212222111211x a x a x a x a x a x a x a x a x a nn n nn n '''+证明 当2=n 时，这时)()()()()()()()()(21122211222112112x a x a x a x a x a x a x a x a x D -===dxx D d ))((2)]()()()([)]()()()([2112221121122211x a x a x a x a x a x a x a x a '-'+'-' )()()()()()()()(2221121122211211x a x a x a x a x a x a x a x a ''+''= 当3=n 时，)()()()()()()()()()()()()()()()(3231222113333123211233322322113x 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 D +-==dx x D d ))((3)()()()()()()()()()()()()()()(323122211333312321123332232211x 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 '+'-' ))()()()()(())()()()()(())()()()()((323122211333312321123332232211'+'-'+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)()()()()()()()()()()()()()()()()()()()()()()()()()()(333231232221131211333231232221131211333231232221131211x 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 a x a x a x a x a x a x a '''+'''+'''= 假设对1-n 阶方阵成立，则对于n 的情形有：)()()()()()()(1112121111x A x a x A x a x A x a x D n n n +++==dxx D d n ))(()]()()()()()([1112121111x A x a x A x a x A x a n n '++'+' )()([1111x A x a '+ )]()()()(111212x A x a x A x a n n '++'+ )()()()()()()()()()()()()()()()()()()()()()()()()()()(212222111211212222111211212222111211x 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 a x a x a x a x a x a x a nnn nn n nn n n n n nn n n n n '''++'''+'''=下面，我们给出定理3.9的证明.证明 因为)()()()()()()()()()()1()1(2)1(12121x y x y x y x y x y x y x y x y x y x W n n n n n n ---'''=对朗斯基行列式求导得：)()()()()()()()()()()()()()()(2)(1)2()2(2)2(12121x y x y x y x y x y x y x y x y x y x y x y x y dxx dW n n n n n n n n n n---'''=分别用)(,),(),(21x p x p x p n n -乘以上述行列式的第一行，第二行，…，第1-n 行后，全部加至第n 行，这时第n 行元素为：),,2,1(,)()()(1)2(2)(n i y x p y x p y x p y i n i n n i n i =+'+++--因为)(,),(),(21x y x y x y n 均是方程（3.4）的解，即),,2,1(,0)()()()(1)2(2)1(1)(n i y x p y x p y x p y x p y i n i n n i n i n i ==+'++++--- 所以),,2,1(,)()()()()1(11)2(2)(n i y x p y x p y x p y x p y n i i n i n n i n i =-=+'+++--- 所以第n 行元素可以换为：),,2,1(,)()1(1n i y x p n i =--因此)()()()()()()()()()()()()()()()()()(1)1(1)1(21)1(11)2()2(2)2(12121x W x p x y x p x y x p x y x p x y x y x y x y x y x y x y x y x y dxx dW n n n n n n n n n n -=---'''=------即)()()(1x W x p dxx dW -= 即dx x p x W x dW )()()(1-= 从0x 到x 积分得：⎰=-xx dtt p ex W x W 01)(0)()(.在前面我们已经得出：如果函数组)(,),(),(21x y x y x y n 是方程0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n在区间),(b a 上的n 个解，则)(,),(),(21x y x y x y n 线性相关⇔)),((0)(0)(00b a x x W x W ∈=⇔≡ )(,),(),(21x y x y x y n 线性无关⇔)),((0)(0)(00b a x x W x W ∈≠⇔≠. 现在刘维尔公式直接给出了n 个解)(,),(),(21x y x y x y n 的朗斯基行列式的值与它在某一点0x x =处的值之间的关系：⎰=-xx dtt p ex W x W 01)(0)()(.下面给出刘维尔公式的一个简单应用：对于二阶线性齐次方程0)()(=+'+''y x q y x p y如果已知它的一个非零解)(1x y ，则由刘维尔公式可以求得与)(1x y 线性无关的另一个解，从而可求得方程的通解.设)(x y 是二阶线性齐次方程0)()(=+'+''y x q y x p y 的与)(1x y 线性无关的解，则由刘维尔公式得：)0()()(11≠⎰=''=-C Ce y y y y x W dxx p即⎰='-'-dxx p Ce y y y y )(11两边同时乘以211y 得： ⎰='-'-dx x p ey C y y y y y )(212111 积分得：⎰⎰⎰=⎰=--dx e y C dx e y C y y dxx p dx x p )(21)(2111 即⎰⎰=-dx e y Cy y dxx p )(2111 因此方程的通解为：⎰⎰+=-dx e y Cy y C x y dxx p )(2111*1)(. 例9 已知方程011)ln 1(2=-'+''-y xy x y x的一个解x y ln 1=，试求其通解.解 这里)1(ln 1)(--=x x x p ，由公式可得通解为：]ln 1[ln ]1[)1(ln 1ln 12*)1(ln 121*1⎰⎰⎰+=⎰+=---dx e xC C x dx e y C C y y x d x dx x x)]ln 1ln 1([ln ]ln 1ln [ln 2*2*⎰⎰⎰-+=-+=dx x dx x C C x dx xx C C x Cx x C xCx C x dx x x xd x x C C x +=+=--+=⎰⎰ln )ln (ln )]ln 1ln 1ln ([ln **2*. 上述例子表明，一个二阶的线性齐次微分方程，如果能得到其一个非零解，利用刘维尔公式可求得该线性齐次微分方程的通解.我们也可以采用换元法，将二阶线性齐次微分方程降阶为一阶微分方程，从而求得该方程的通解，下面介绍一下换元法.对于二阶线性齐次方程0)()(=+'+''y x q y x p y如果已知它的一个非零解)(1x y ，下面做变量代换，令z y y 1=则z y z y z y y z y z y y ''+''+''='''+'='111112,，代入原方程得： 0)())((2111111=+'+'+''+''+''z y x q z y z y x p z y z y z y 即0))(2())()((111111=''+'+'++'+''z y z y x p y z y x q y x p y 因为)(1x y 是原方程的解，所以0)()(111=+'+''y x q y x p y ，从而上式可化为： 0))(2(111=''+'+'z y z y x p y 令u z ='，则u z '=''，代入得：0))(2(111='++'u y u y x p y 降为了一阶微分方程，而且是变量可分离的方程，整理得：dx y y x p y u du111)(2+'-=102积分得：⎰=+'-dxy y x p y Ceu 111)(2.从而得到⎰⎰⎰=⎰=+'-+'-dx Cey y dx Cez dxy y x p y dxy y x p y 111111)(21)(2,.因此方程的通解为：⎰⎰+=+'-dx eCy y C x y dx y y x p y 111)(211*)(.例10 求方程066323=-'+''-'''y y x y x y x的通解，已知它的两个特解221,x y x y ==.解 令xz y =，则z z x y z z x y z z x y ''+'''=''''+''=''+'='3,2, 代入066323=-'+''-'''y y x y x y x 得：06)(6)2(3)3(23=-+'+'+''-''+'''xz z z x x z z x x z z x x即04='''z x ，所以,,,12321x z x z z ===因此33x y =故原方程的通解为：33221x C x C x C y ++=.习 题 3.21.试讨论下列各函数组在它们的定义区间上是线性相关的还是线性无关的？（1）;sin ,cos ,2sin t t t （2）;tan ,cos ,sin x x x （3）;42,2,322+++-x x x x x （4）.,,2t t t e t te e2.设在方程0)()(=+'+''y x q y x p y 中，)(x p 在某区间I 上连续且恒不为零，103试证它的任意两个线性无关的解的朗斯基行列式是区间I 上的严格单调函数.3.试证明：二阶线性齐次方程的任意两个线性无关解的朗斯基行列式之比是一个不为零的常数.4.已知方程022)1(2=+'-''-y y x y x 的一个解x y =1，试求其通解.5.已知方程0)1(=+'-''-y y x y x 的一个解x y =1，试求其通解.6.在方程0)()(=+'+''y x q y x p y 中，当系数满足什么条件时，其基本解组的朗斯基行列式等于常数.7.设)(1x y 是n 阶线性齐次方程0)()()(1)1(1)(=+'+++--y x a y x a y x a y n n n n的一个非零解.试证明：利用线性变换z x y y )(1=可将已知方程化为1-n 阶的齐次方程.3.3 n 阶线性非齐次微分方程的一般理论对于线性非齐次微分方程)()()()(1)1(1)(x f y x p y x p y x p y n n n n =+'+++-- （3.2）而言，我们首先研究其解与对应的线性齐次微分方程解之间的关系.（1）如果)(x y 是线性齐次微分方程（3.4）的解，)(*x y 是线性非齐次微分方程（3.2）的解，则)()(*x y x y +是线性非齐次微分方程（3.2）的解.证明 因为)(x y 是线性齐次微分方程（3.4）的解， 所以0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n因为)(*x y 是线性非齐次微分方程（3.2）的解， 所以)()())(())(()(**1)1(*1)(*x f y x p y x p y x p y n n n n =+'+++--因此))(())(())(()(**1)1(*1)(*y y x p y y x p y y x p y y n n n n ++'++++++--104))(())()(())()(()(**1)1(*)1(1)(*)(y y x p y y x p y y x p y y n n n n n n ++'+'+++++=--- )1(*1)(*1)1(1)())(()(())()()((---+++'+++=n n n n n n y x p y y x p y x p y x p y )()(0))())((**1x f x f y x p y x p n n =+=+'++- . 即，)()(*x y x y +是线性非齐次微分方程（3.2）的解.（2）如果)(),(*2*1x y x y 是线性非齐次微分方程（3.2）的解，则)()(*2*1x y x y -是线性齐次微分方程（3.4）的解.证明 因为)(),(*2*1x y x y 是线性非齐次微分方程（3.2）的解，所以)()())(())(()(*1*11)1(*11)(*1x f y x p y x p y x p y n n n n =+'+++-- )()())(())(()(*2*21)1(*21)(*2x f y x p y x p y x p y n n n n =+'+++--因此))(())(())(()(*2*1*2*11)1(*2*11)(*2*1y y x p y y x p y y x p y y n n n n -+'-++-+---))(())())((())())((()()(*2*1*2*11)1(*2)1(*11)(*2)(*1y y x p y y x p y y x p y y n n n n n n -+'-'++-+-=--- )1(*21)(*2*1*11)1(*11)(*1))(()(())())(())(()((---+-+'+++=n n n n n n y x p y y x p y x p y x p y 0)()())())((*2*21=-=+'++-x f x f y x p y x p n n . 即，)()(*2*1x y x y -是线性齐次微分方程（3.4）的解.根据上述两条性质，我们可以得到下面的定理. 定理3.10 n 阶线性非齐次微分方程)()()()(1)1(1)(x f y x p y x p y x p y n n n n =+'+++--的通解等于它对应的齐次方程0)()()(1)1(1)(=+'+++--y x p y x p y x p y n n n n的通解与它本身的一个特解之和.即)()()()()(*2211x y x y C x y C x y C x y n n ++++=证明 因为)()()(2211x y C x y C x y C n n +++ 是齐次方程（3.4）的解，)(*x y 是105非齐次方程（3.2）的解，由性质（1）知，)()()()()(*2211x y x y C x y C x y C x y n n ++++=是非齐次方程（3.2）的解.下面证明非齐次方程（3.2）的任一解)(x y 都可以表示成)()()()(*2211x y x y C x y C x y C n n ++++事实上，因为)(x y 和)(*x y 都是非齐次方程（3.2）的解，由性质（2）知，)()(*x y x y -是齐次方程（3.4）的解，所以可表示成方程（3.4）的通解形式，即)()()()()(2211*x y C x y C x y C x y x y n n +++=-因此)()()()()(*2211x y x y C x y C x y C x y n n ++++= .故非齐次方程（3.2）的通解可以表示为齐次方程的通解与非齐次方程的特解之和.由定理3.10知，求解一个线性非齐次方程（3.2）的关键是先找到对应的线性齐次方程（3.4）的通解，再找到一个非齐次方程的特解就可以了，这时有非齐通解=齐通解+非齐特解.假定我们已经求得了线性齐次微分方程（3.4）的齐通解)()()()(2211x y C x y C x y C x y n n +++=这时，可以采用常数变易法求非齐次方程的特解)(*x y .下面来看一下n 阶线性非齐次微分方程的常数变易法.已知)()()()(2211x y C x y C x y C x y n n +++= 是齐次方程（3.4）的通解，设)()()()()()()(2211*x y x C x y x C x y x C x y n n +++=是线性非齐次方程（3.2）的一个特解，为了将)(*x y 代入方程（3.2），我们需要首先计算)(*x y 的一阶，二阶，…，n 阶导数.则)()()()([)]()()()()()([))((22112211*x y x C x y x C x y x C x y x C x y x C x y n n '+'+'++'+'='106)]()(x y x C nn '++ 在求二阶导数之前，我们先研究一下特解)(*x y ，要想得到)(*x y ，必须求得)(,),(),(21x C x C x C n ，但是代入方程（3.2）只能得到一个等式，所以我们必须构造)(,),(),(21x C x C x C n 满足的另外1-n 个等式，因此在))((*'x y 中，令0)()()()()()(2211='++'+'x y x C x y x C x y x C n n 则)()()()()()())((2211*x y x C x y x C x y x C x y n n '++'+'=' 这时，再求))((*''x y ，有)()()()([)]()()()()()([))((22112211*x y x C x y x C x y x C x y x C x y x C x y n n ''+''+''++''+''='' )]()(x y x C n n ''++ 再令0)()()()()()(2211=''++''+''x y x C x y x C x y x C n n 则)()()()()()())((2211*x y x C x y x C x y x C x y n n ''++''+''='' 序行此法，可得)()([)]()()()()()([))(()1(11)2()2(22)2(11)1(*x y x C x y x C x y x C x y x C x y n n n n n n n -----+'++'+'= )]()()()()1()1(22x y x C x y x C n n n n --+++令0)()()()()()()2()2(22)2(11='++'+'---x y x C x y x C x y x C n n n n n则)()()()()()())(()1()1(22)1(11)1(*x y x C x y x C x y x C x y n n n n n n ----++++=这时，再求)(*))((n x y ，有)()([)]()()()()()([))(()(11)1()1(22)1(11)(*x y x C x y x C x y x C x y x C x y n n n n n n n +'++'+'=--- )]()()()()()(22x y x C x y x C n n n n +++107最后将求得的)(*x y 的一阶，二阶，…，n 阶导数代入方程（3.2）得：)()()()([)]()()()()()({[)(22)(11)1()1(22)1(11x y x C x y x C x y x C x y x C x y x C n n n n n n n ++'++'+'--- )]()()()()()()[()]}()()1()1(22)1(111)(x y x C x y x C x y x C x p x y x C n n n n n n n n ---++++++)()()[()]()()()()()()[(1122111x y x C x p x y x C x y x C x y x C x p n n n n +'++'+'++- )()]()()()(22x f x y x C x y x C n n =+++ 即)]()()()()()()()[(111)1(11)(11x y x p x y x p x y x p x y x C n n n n +'+++-- ++'++++--)]()()()()()()()[(221)1(21)(22x y x p x y x p x y x p x y x C n n n n )]()()()()()()()[(1)1(1)(x y x p x y x p x y x p x y x C n n n n n n n n n +'++++-- )()]()()()()()([)1()1(22)1(11x f x y x C x y x C x y x C n n n n n ='++'+'+---因为),,2,1(,0)()()()()()()(1)1(1)(n i x y x p x y x p x y x p x y i n i n n i n i ==+'+++-- 所以)()()()()()()()1()1(22)1(11x f x y x C x y x C x y x C n n n n n ='++'+'---这样，就得到了)(,),(),(21x C x C x C n ''' 的n 个方程 ⎪⎪⎪⎩⎪⎪⎪⎨⎧='++'+'='++'+'=''++''+''='++'+'------)()()()()()()(0)()()()()()(0)()()()()()(0)()()()()()()1()1(22)1(11)2()2(22)2(1122112211x f x y x C x y x C x y x C x y x C x y x C x y x C x y x C x y x C x y x C x y x C x y x C x y x C n n n n n n nn n n n n n n 该方程组的系数行列式为)(,),(),(21x y x y x y n 的朗斯基行列式)(x W ，因为这n 个解是线性无关的，所以0)(≠x W ，因此，该方程组存在唯一解，这样就可以求得)(,),(),(21x C x C x C n ''' ，再逐个积分，求得)(,),(),(21x C x C x C n ，从而得到特解)(*x y .例1 求方程10812-=-''x xe e y y的通解.解 齐次方程0=-''y y的通解为：x x e C e C y -+=21.设xxe x C e x C x y -+=)()()(21*为方程12-=-''x xe e y y 的一个特解，则⎪⎩⎪⎨⎧-='-'='+'--12)()(0)()(2121x x xx x x e e e x C e x C e x C e x C ， 解方程组得：⎪⎩⎪⎨⎧-='-='x xx e e x C e x C 1)(11)(221，所以x e de ee de e e dx e x C xx x x x x x x --=--=-=-=⎰⎰⎰1ln )111()1(111)(1 )1(ln )111(11)(22xx x x x x x x x e e de e de e e dx e e x C +--=+--=-=-=⎰⎰⎰ 因此非齐特解为)11ln (1ln )(*+----=-x x x x x e e xe e e x y故方程的通解为)11ln (1ln )(21+----++=--x x x x x x x e e xe e e e C e C x y .例2 求方程)(2t f x x =+∙∙ω的通解. 解 对应的齐次方程020=+∙∙x x ω的通解为t C t C t x 0201sin cos )(ωω+=.。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第三讲一阶线性非齐次微分方程组的一般理论(2课时)
一、目的与要求：理解一阶线性非齐次方程组的一般理论，掌握一阶线性非齐次方程组的通
解结构，理解常数变易法.
二、重点：一阶线性非齐次方程组的通解结构，常数变易法.
三、难点：常数变易法.
四、教学方法：讲练结合法、启发式与提问式相结合教学法.
五、教学手段：传统板书与多媒体课件辅助教学相结合.
六、教学过程：
1.课题引入
本节研究一阶线性非齐次方程组
dY
—=A(X Y+F x ) dx (3.7
)
的通解结构与常数变易法.
2.通解结构
定理3.8如果Y(x)是线性非齐次方程组(3.7)的解，而Y0(x)是其对应齐次方程组(3.8)
的解，则Y0(x) +Y(x)是非齐次方程组(3.7)的解.资料个人收集整理，勿做商业用途
证明这只要直接代入验证即可
定理 3.9线性非齐次方程组(3.7)的任意两个解之差是其对应齐次方程组(3.8)的解.
证明设Y(X)和Y(x)是非齐方程组(3.7)的任意两个解，即有等式
dY-^ =A(x)Y(x)+F(x) ，■dY(x) = A(X)Y(X)+F(X)dx
dx
于是有
知(x)-Y(x)]=皿-如
dx dx dx
=A( x)Y(X)+ F(X)- A(X)Y(X) — F (X)
= A(x)[ Y(X)-Y(X)]
上式说明Y(X)-Y(X)是齐次方程组(3.8)的解.
定理3.10 线性非齐次方程组(3.7)的通解等于其对应的齐次方程组(3.8)的通解与方程
组(3.7)的一个特解之和.即若Y(x)是非齐次方程组(3.7)的一个特解，YdhKx)，川，YJ X)
是对应齐次方程组(3.8)的一个基本解组，则方程组(3.7)的通解为
Y(X)-CMd) +C2Y2(X)+川+C n Y n(x pH Y'(x)资料个人收集整理，勿做商业用途
这里C i,C2,HI,C n是任意常数.
证明首先由定理3.8，不论C i,C2,川，C n是什么常数，(3.16)都是(3.7)的解.其次对于方程组(3.7)的任何一个解Y(x)，由定理3.9知，是Y(X)-Y(X)对应齐次方程组的解.于是
由基资料个人收集整理，勿做商业用途
本定理3.6，存在常数G'G,川,c n使得
Y(x) -Y(X)=CY I(X)+ CY2(X) +川+ CY n(x)
Y(x) = CY I(X)+C Y2(X)+ 川+ CY n(x) + Y(X)
所以(3.16)是(3.7)的通解.定理证毕.
3.拉格朗日常数变易法
在第一章我们介绍了对于一阶线性非齐次方程，可用常数变易法求其通解.现在，对于
线性非齐次方程组，自然要问，是否也有常数变易法求其通解呢？事实上，定理3.10告诉我们，为了求解非齐次方程组(3.7)，只需求出它的一个特解和对应齐次方程组 (3.8)的一个基本解组.而当(3.8)的基本解组已知时，类似于一阶方程式，有下面的常数变易法可以求得(3.7) 的一个特解.资料个人收集整理，勿做商业用途
为了计算简洁，我们定义(3.8)的基本解矩阵如下:
l_yn1 (x) y n2(X)lli y nn (x).
其中每一列均为(3.8)的解Y i (x)(i =1,2,川，n)，且Y ,(X ),Y 2(X )川，Y n (x)是(3.8)的一个基本 解组.因此det ①(X)=W(x) H0.
由定理3.6知，齐次方程组(3.8)的通解可表为
Y(x)=e(x)C ,
其中C 为列向量
「C J
C 2
L Cn J
它的各个分量C i (i =1,2，川，n)为任意常数.现在求(3.7)的形如
Y(x) =6(x)C(x)
的解，其中
「C i (x)「
LC n (x)J
为待定向量函数.将(3.17)代入(3.7)有
①'(x)C(x) +①(x)C'(x) = A(x)①(X)+F(x)
其中
①(x) =
y 2i (x) y 22(x)川 y 2n (x) (3.17)
C(x) = C 2(x)
Ly ni (x) y n2(x)川 y nnW” 因为①(x)是(3.8)的基本解矩阵，所以有 ①'(X)=A(x)①(X).从而，上式变为
①(x)C'(x) =F(x)
由于①(x)是非奇异矩阵，故①」(x)存在，于是
C(x)=①」(x)F(x)
积分得
x
C(x) = f a(t)F(t)dt
x 0
X 0为I 中任一点代入(3.17)得到
Y(x) = J x ①(X)①'(t)F(t)dt
x 0
显然Y(x)是(3.7)的一个特解，于是得到非齐次方程组 (3.7)的通解公式
x
Y(x) =e(x)C + J ①(X)①'(t)F(t)dt
x 0
例1求解方程组x=y —5cost, y=2x + y
解 由3.3节例4知，向量函数组
是对应齐次方程组的基本解组.现在求非齐次方程组形如
的特解，此时(3.18)的纯量形式为
①(x) = y 2i (x) y 22(x)川 y 2n
(x) (3.18) (3.19)
|C/(t)e-^+C2(t)e2^-5cost
[-C1(t)e」+2C；(t)e2t =0
解之得
10 5
C i'(t) = ——e t cost, C/t) =-e't cost
3 3
从而
5 1
G(t) =e (cost +sint), C2(t^-e^t(2cos t-s int)
3 3
最后可得该方程组的通解为
J x(t) rGe丄+C2e2t —cost —2sin t
[y(t) = -C i e~^ +2C2e2^3cos^sin t
本讲要点：
1 .非齐次通解=对应齐次通解+非齐次一个特解
2.常数变易法适用于：
先求出齐次通解Y =①(x)C ,再令W =0(x)C(x)为非齐次特解代入原方程确定C(x).。