【典型例题】 第五章 线性微分方程组
- 格式:doc
- 大小:809.50 KB
- 文档页数:22
第五章:线性微分方程组本章教学目的和要求:使学生掌握线性微分方程组解的结构。
要求学生熟练掌握求解常系数线性问粉方程组。
熟练掌握常数变易法。
本章重点:解的性质与结构,常系数方程组的解法,常数变易法。
本章难点:向量函数组的线性相关性,一般理论中的定理证明。
本章课时安排:讲16学时,习题及总结测验2学时第五章:线性微分方程组说明:本章所讨论的线性微分方程组仅限与一阶微分方程,从讲义的开头所说的,方程组不仅能在实际中应用广泛,而且她对高阶方程的求解具有不可忽视的作用。
不仅如此,方程组的有关定理在近代微分方程理论中也占有重要地位。
本章内容:一.一阶微分线性方程组及其解的概念;初值问题解的存在和唯一性定理。
二.线性方程组及其解的一般理论/包括解的线线性相关性,基本解组和解的结构定理。
三.方程组的具体解法。
§5.1 存在唯一性定理5.1.1 记号和定义①引言:在第二章我们研究了含有一个未知函数的微分方程的解法以及它们的性质。
但是,在很多实际问题与理论问题中,还要求我们去求解含有多个未知数函数的微分方程组,或者研究它们的解的性质。
如空间运动质点P 的速度与t 以及坐标(,,)x y z 的关系式为:112232(,,,)(,,,)(,,,)x y z v f t x y z x f v f t x y z y f z f v f t x y z ⎧==⎧⎪⎪=⇒=⎨⎨⎪⎪==⎩⎩ 又如: 22sin d dt l θθθ=-令 sin d dtd dtl θωωθθ⎧=⎪⎪⎨⎪=-⎪⎩化成一阶微分方程组。
用类似的方法,如果在 n 阶微分方程 ()(1)(,,...,)n n y x y y y -'=中,令(1)121.,,...,n n y y y y y y --'''=== 它就可以化成方程组 1212(1)121()(1),........(,,...,)n n n n n n y y y y y y y y y yy x y y y -----⎧'=⎪'''==⎪⎪⎨⎪'==⎪⎪'=⎩共同点:出现的未知函数的导数都是一阶的 它 们都是一阶微分方程组。
常微分方程1 .( 05,4 分)微分方程xy 2yxln x 满足y(1)22x y)= x ln x.2 .( 06,4 分) 微分方程 y= y(1 x)的通解为 ———— x分析:这是可变量分离的一阶方程,分离变量得dy( 11)dx.积分得 ln y ln x x C 1,即 y e C1xe x yxy Cxe x, 其中C 为任意常数 .(二)奇次方程与伯努利方程1 .( 97,2,5 分) 求微分方程 (3x2 2xy y 2)dx (x 22xy)dy 0的通解解:所给方程是奇次方程 . 令 y=xu, 则 dy=xdu+udx. 代入原方程得 3 ( 1+u- u 2) dx+x(1-2 u) du=0. 分离变量得1-2u2 du 3dx, 1uu x积分得 ln 1 u u 2 3ln x C 1,即 1 u u 2=Cx 3. 以 u y代入得通解 x 2xy y 2.xx( y x 2y 2)dx xdy 0(x 0),2 .(99,2,7 分 ) 求初值问题 的解 .y x1 0分析:这是一阶线性微分方程原方程变形为 . dy +2y dx x 2 dx lnx, 两边乘 e x=x 得积分得y(1)x 2y=C+ x 2 ln xdx C 1 ln xdx 3 3 1 11 得 C 0 y xln x x.9 39 C 1 x 3 ln x 3 13 x. 9 1 的解解:所给方程是齐次方程 (因 dx, dy 的系数 (y+ x 2 y 2)与 (-x)都是一次齐次函数)令 dy xdu udx,带入得x(u 1 u 2dx x( xdu udx) 0, 化简得 12u 2dx xdu 0.分离变量得dx- du=0. x 1 u 2积分得 ln x ln(u 1 u 2) C 1,即 u 1 u 2Cx. 以 u y代入原方程通解为y+ x 2 y 2 Cx 2.x 再代入初始条件 y x 1 0,得 C=1.故所求解为 y+x 2y2x 2,或写成y 12 (x 2 1).(三)全微分方程 练习题(94,1,9 分)设 f ( x)具有二阶连续导数, f (0) 0, f (0) 1,且 [xy(x+y)- f(x)y]dx+[ f (x)+x 2y]dy=0为一全微分方程,求 f(x)以及全微分方程的通解先用凑微分法求左端微分式的原函数:122 122( y dx x dy ) 2( ydx xdy ) yd (2sin x cos x) (2sin x cos x)dy 0, 22 122d [ x y 2xy y (cos x 2sin x)] 0. 2其通解为 1x 2y 2 2xy y (cos x 2sin x) C.4.( 98,3分) 已知函数y y(x)在任意点x 处的增量 y= y2 x ,当 x0时 ,1x是 x 的高阶无穷小,y(0)= ,则 y(1)等于 ( )解:由全微分方程的条件,有 即 x22xy f (x) f (x)y因而 f (x)是初值问题y x 2[xy(x y) f(x)y] y 2xy, 亦即 f (x) f (x) x 2.2yx的解,从而解得0, y x 0 12.22[ f (x) xy], x 2sin x cosx)dy 0.(A)2 .(B) .(C)e 4 .(D) e 4 .分析:由可微定义,得微分方程 y y. 分离变量得21x1y dx2,两边同时积分得 ln y arctan x C ,即 y Ce arctanx.y1x代入初始条件y(0) ,得 C= ,于是 y(x) earctanx,由此, y(1) e 4.应选 ( D)二、二阶微分方程的可降阶类型5( . 00,3分) 微分方程 x y 3y 0的通解为分析:这是二阶微分方程的一个可降阶类型,令 y =P( x),则 y =P ,方程可化为一阶线性方程xP 3P 0,标准形式为 P+3P=0,两边乘 x 3得 (Px 3) =0. 通解为 y P C 30 .xx再积分得所求通解为 y C 22C 1.x216 .( 02,3分)微分方程 yy y 2=0满足初始条件y x 01, y x 0 2的特解是分析:这是二阶的可降阶微分方程 .令 y P(y)(以 y 为自变量 ),则 y dy dP P dP.dx dx dy代入方程得 yP dP +P 2=0,即 y dP+P=0(或 P=0, ,但其不满足初始条件y x 0 1)dy dy2分离变量得 dP dy 0,PyC积分得 ln P +ln y =C ,即 P= 1(P=0对应 C 1=0); y11由 x 0时 y 1, P=y , 得 C 1 ,于是221 y P ,2 ydy dx, 积分得 y x C 2 2y .又由 y x 0 1 得 C 2. 1,所求特解为 y 1 x.三、二阶线性微分方程(一)二阶线性微分方程解的性质与通解结构7 .( 01,3分)设 y e x(C 1sin xC 2cosx)(C 1,C 2为任意常数 )为某二阶常系数线性齐次微分方程的通解,则该方程为 ___ .r1,r2 1 i,从而得知特征方程为分析一:由通解的形式可得特征方程的两个根是22(r r1 )(r r2) r (r1 r2 )r r1r2 r 2r 2 0.由此,所求微分方程为y 2y 2y 0.分析二:根本不去管它所求的微分方程是什么类型(只要是二阶),由通解y e x(C1sinx C2 cosx)求得y e x[( C1 C2 )sin x (C1 C2)cos x], y e x( 2C2 sin x 2C1 cos x),从这三个式子消去C1与C2,得y 2y 2y 0.(二)求解二阶线性常系数非齐次方程9.( 07,4分) 二阶常系数非齐次线性微分方程y 4y 3y 2e2x的通解为y=分析:特征方程24 3 ( 1)( 3) 0的根为1, 3.非齐次项 e x, 2不是特征根,非齐次方程有特解y Ae2x.代入方程得(4A 8A 3A)e2x2e2x A 2.因此,通解为y C1e x C2e3x2e2x..10.(10,10分 )求微分方程y 3y 2y 2xe x的通解.分析:这是求二阶线性常系数非齐次方程的通解.1由相应的特征方程2 3 2 0, 得特征根 1 1, 2 2 相应的齐次方程的通解为y C1e x C2e2x.2非齐次项 f ( x) 2xe x , 1是单特征根,故设原方程的特解xy x(ax b)e .代入原方程得ax2 (4a b)x 2a 2b 3[ax2 (2a b)x b] 2(ax2 bx) 2x,即 2ax 2a b 2x, a 1,b 2.3原方程的通解为y C1e x C2e2x x(x 2)e x,其中 C1,C2为两个任意常数.04, 2, 4分)微分方程y y x2 1 sin x的特解形式可设为( )22(A)y ax bx c x(Asin x B cosx).(B)y x(ax bx c Asin x B cos x).22(C)y ax bx c Asin x.(D )y ax bx c Acosx.分析:相应的二阶线性齐次方程的特征方程是2 1 0,特征根为i .y y x2 1L()与 1 y y sin xL( 2)方程 (1) 有特解 y ax2 bx c,方程(2)的非齐次项 f (x) e x sin x sin x( 0, 1,i 是特征根), 它有特解y x(Asin x B cosx).y ax2 bx c x(Asin x Bbcosx).应选 (A).(四)二阶线性变系数方程与欧拉方程12.(04, 4分 )欧拉方程x2 d2y 4x dy 2y 0(x 0)的通解为dx dx分析:建立 y 对 t 的导数与y 对 x 的导数之间的关系 .222dy dy dx dyd y d y 2 dy 2 d y dy( sin x), 2 2 sin t cost (1 x ) 2 x .dt dx dt dx dt dx dx dx dxd 2y于是原方程化为 2 y 0,其通解为 y C 1 cost C 2sint.dt 2 回到 x 为自变量得 y C 1x C 2 1 x 2.x由 y (0) C 2 1 C 2 1.y(0) C 1x 02 C 1 2.1 x 2因此 特解为 y 2x 1 x 2 .四、高于二阶的线性常系数齐次方程13.( 08, 4分)在下列微分方程中,以 y C 1e xC 2cos2x C 3 sin 2x(C 1, C 2, C 3为任意常数)为通 解的是()(A)y y 4y 4y 0.(B)y y 4y 4y 0. (C)y y 4y 4y 0.(D ) y y 4y 4y 0.分析:从通解的结构知,三阶线性常系数齐次方程相应的三个特征根是: 1, 2i(i 1),对 应的特征方程是 ( 1)( 2i)( 2i) ( 1)( 24) 3244 0,因此所求的微分方程是 y y 4y 4y 0,选(D).(00,2,3分 ) 具有特解 y 1 e x , y 2 2xe x ,y 3 3e x的三阶常系数齐次线性微分方程是( )(A)y y y y 0.(B)y y y y 0. (C)y 6y 11y 6y 0.(D)y2y y 2y 0.分析:首先,由已知的三个特解可知特征方程的三个根为 r 1 r 21,r 3 1,从而特征方程为(1)求导数 f (x); (2)证明:当 x 0时 ,成立不等式 e分析:求解欧拉方程的方法是:作自变量22d y dy d y dy 2 (4 1) 2y 0,即 2 3 2y xe t(t l n x),将它化成常系数的情形: 0.1, 2 2, 通解为 yC 1e t C 2e 2t. y C 1 x C 22,其中C 1,C 2为任意常数(05,2,12分 )用变量代换 xcost (0 t)化简微分方程 (1 x 2)y xy y 0,并求其(r 1)2(r 1) 0,即r3r 2r 1 0,由此,微分方程为y y y y 0.应选(D).五、求解含变限积分的方程00, 2,8分) 函数y=f(x)在0, 上可导,f (0) 1,且满足等式1xf (x) f (x) 1 f (t)dt 0,x10f(x) 1.求解与证明()首先对恒等式变形后两边求导以便消去积分: 1x(x 1)f (x) (x 1)f(x) 0f (t)dt 0,(x 1)f (x)(x 2)f (x)0.在原方程中令变限 x 0得 f (0) f (0) 0,由 f (0) 1,得 f (0) 1.现降阶:令 u f (x),则有 u x 2u 0,解此一阶线性方程得x1x e f (x) u C eu 0x1 x e 由 f (0) 1,得 C 1,于是 f (x) e. x1xe (2)方法 1 用单调性 . 由f (x) e0(x 0), f (x)单调减 , f(x) f(0) 1(x );x1x 又设 (x) f (x) e x ,则 (x) f (x) e x x e x0(x 0), (x)单调增,因此 (x)x1 (0) 0(x 0),即 f(x) e x(x 0) . 综上所述,当 x 0时 ,e x f (x) 1.方法 2 用积分比较定理 . 由 牛顿 -莱布尼茨公式,有六、应用问题 (一)按导数的几何应用列方程 练习题 1 .( 96,1,7分)设对任意 x 0,曲线 y f(x)上点 (x, f(x))处的切线在 y 轴上的截距等于1 xf (t)dt,求 f ( x)的一般表达式 . x 0解:曲线 y f (x)上点 (x, f ( x))处的切线方程为 Y f ( x) f ( x)( X x).令 X 0得 y 轴上的截距 Y f(x) xf (x).由题意 1x1f(t)dt f(x) xf (x) x 0x, 得x 2f(t)dt xf (x) x 2f (x)( ) 恒等式两边求导,得 f (x) f (x) xf (x) 2xf (x) x 2f ( x),即 xf (x) f (x) 0 在 ( )式中令 x 0得 0 0,自然成立 . 故不必再加附加条件. 就是说f (x)是微分方程 xy y 0的通解 . 令 y P(x),则 y P ,解 xP P 0,得 y P C 1.xf ( x) f (0) x0 f (t)dt, f(x) t 由于 0 e t1从而有 e x e t (t 0),有 0 f (x) 1. 0t e t d t 1 dt . 1 x t e t dt x e (x再积分得 y f ( x) C1 ln x C2.12( . 98,2,8分) 设 y y(x)是一向上凸的连续曲线 ,其上任意一点 (x, y)处的曲率为 1,1 y 2y P tan( x).(二 )按定积分几何应用列方程3.(97,2,8分 )设曲线 L 的极坐标方程为 r r( ), M (r, )为 L 上任一点 ,M 0(2,0)为 L 上一定点 ,若极径 OM 0,OM 与曲线 L 所围成的曲边扇形面积值等于 L 上 M 0、 M 两点间弧长值的一半, 求曲线L 的方程 .且此曲线上点 (0,1)处的切线方程为 y x 1, 求该曲线的方程,并求函数 y y( x)的极值 .解:由题设和曲率公式有y( x)向上凸 , y 0, y令 y P(x),则 y P ,方程化为 y) ,化简得 y 12. yP1 P 21, dP 分离变量得 2 dx,积分得C 1.y (0) 1即 P(0) 1,代入可得 C 1,故再积分得 y ln cos( x) C 2 又由题设可知y(0)1,代入确定 C 2 11ln 2,1y ln cos( x) 1 ln 2x , 即当 4 2,3时 ,cos( x) 0, 而3 或 时, 44cos( x)y ln cos( 40,ln cos( x)1 x) 12 ln2( 4 x34 )显然,当 x 时 ,ln cos( x) 4410, y 取最大值 1 1ln 2,显然 y 在 (3),没有极小值解:由已知条件得r 2d r 2 r 2d , 2020 两边对 求导 ,,得 r 2 r 2 r (隐式微分方程)2 ,解出 r r r 2 1,从而, L 的直角坐标方程为 x m 3y 2.1 arccos r 分离变量,得 dr r r 2 dr r r 2 1 d 1 1 d( )1 r (r 1)2 arccos 1 , 或 r dr r r 2 1d tarccos 1(r sect ) 两边积分,得 代入初始条件 r(0) 2,得 1arccos 2 1arccos r3L 的极坐标方程为 1 r cos( ) 31 co s 3si。
第五章 线性微分方程组5-1 考虑方程组x A x )(t dtd = (1)其中)(t A 是区间b t a ≤≤上的连续n n ⨯矩阵,它的元素为n j i t a ij ,,2,1,),( =,1)如果)(,),(),(21t t t n x x x 是(1)的任意n 个解,那么它们的朗斯基行列式)()](,),(),([21t W t x t x t x W n ≡ 满足下面的一阶线性微分方程W t a t a t a W nn )]()()([2211+++=' (2); 2)解上面的一阶线性微分方程,证明下面的公式:],[,,)()(0)]()([0011b a t t et W t W tt nn dss a s a ∈=⎰++ 。
证 1)根据行列式的微分公式)()()()()()()()()()()()()()()()()()()(122111112211111221111t x t x t x t x t x t x t x t x t x t x t x t x t x t x t x t x t x t x t W nnnn n nn n n n nn n n n ''++''+''='(3)由于)(,),(),(21t t t n x x x 是(1)的解,所以⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛=⎪⎪⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎪⎪⎭⎫⎝⎛='∑∑∑===nj jk nj nj jk j n j jk j nk k k nn n n n k t x t a t x t a t x t a t x t x t x t a t a t a t a t a t a t 11211211221111)()()()()()()()()()()()()()()()(x , 所以∑==='nj jk ijikn k i t x t at x 1),,2,1,(),()( )(,把这些等式代入(3)的右端,化简计算每个行列式,如(3)式右端第一项等于)()()()()()()()()()()()()()()()()(11122111111122111111t W t a t x t x t x t x t x t x t a t x t x t x t x t x t at x t ann n n n nn n n nj jn jn j j j==∑∑==类似地可以算出(3)式右端其它各项分别为)()(,),()(22t W t a t W t a nn ,代入(3)得W t a t a t a W nn )]()()([2211+++=' (2)2)方程(2)是关于)(t W 的一阶线性微分方程,分离变量可求得通解为 ⎰++=tt nn dss a s a Cet W 011)]()([)( ,C 为任意常数。
若)(,00t W W t t ==,则 )(0t W C =, 于是 ⎰++=tt n n dss a s a et W t W 011)]()([0)()( 。
评注:公式 ⎰++=tt nn dss a s a et W t W 011)]()([0)()( 称为刘维尔公式,反映了线性齐次方程组的解与系数矩阵)(t A 的关系。
)()()()(22111t a t a t a t a nn ni ii +++=∑= 称为矩阵)(t A 的迹,记为)(t tr A ,所以刘维尔公式又可表示为⎰=tt dss tr et W t W 0)(0)()(A 。
从公式中可以看出,线性齐次方程组(1)的n 个解构成的朗斯基行列式)(t W 或者恒为零,或者恒不为零。
5-2 设)(t A 为区间b t a ≤≤上连续的n n ⨯实矩阵,)(t Φ为方程x A x )(t ='的基本解矩阵,而)(t φx =为其一解。
试证:1) 对于方程y A y )(t T-='的任一解)(t ψy =必有=)()(t t T φψ常数;2) )(t Ψ为方程y A y )(t T -='的基本解矩阵的充要条件是存在非奇异的常数矩阵C ,使C ΦΨ=)()(t t T。
证 1) 由于)(t ψy =为方程y A y )(t T-='的解,则)()()(t t t TψA ψ-=',两边转置,得())()()(t t t TT A ψψ-=',即())()()(t t t T T A ψψ-='。
因为[]())()()()()()(t t t t dtt t d TTTφψφψφψ'+'=()())()()()()()(t t t t t t TTφA ψφA ψ+-=0= ,所以必有=)()(t t T φψ常数。
2) 必要性。
由于)(t Ψ为方程y A y )(t T -='的基本解矩阵,则 )()()(t t t T ΨA Ψ-=', 转置后,得 ())()()(t t t T T A ΨΨ-='。
因为[]())()()()()()(t t t t dtt t d TTTΦΨΦΨΦΨ'+'=()())()()()()()(t t t t t t TTΦA ΨΦA Ψ+-=0= (零矩阵)。
所以 C ΦΨ=)()(t t T (常数矩阵),而)(t Ψ和)(t Φ都是基本解矩阵,因而C 还为非奇异矩阵。
充分性。
由于存在非奇异的常数矩阵C ,使C ΦΨ=)()(t t T,两边关于t 求导数,有[]())()()()()()(t t t t dtt t d TTTΦΨΦΨΦΨ'+'=()0=+'=)()()()()(t t t t t TT ΦA ΨΦΨ即 ())()()()()(t t t t t TT ΦA ΨΦΨ-=',而)(t Φ是基本解矩阵,则)(t Φ为非奇异矩阵,故有())()()(t t t TTA ΨΨ-=',即())()()(t t t TTA ΨΨ-=',两边再转置,得)()()(t t t TΨA Ψ-=',即证明了)(t Ψ为方程y A y )(t T-='的基本解矩阵。
评注:由证明过程可以看出,方程y A y )(t T-='和x A x )(t ='的解曲线之间满足=)()(t t Tφψ常数。
5-3 设)(t Φ是n 阶线性方程组Ax x =dtd (A 是n n ⨯的常数矩阵)的标准基本解矩阵,(即E Φ=)0()证明)()()(001t t t t -=-ΦΦΦ其中0t 为某一值。
证 因)(t Φ为基本解矩阵,则有)()(t dtt d A ΦΦ=,0)(det ≠t Φ)()()(000t t t t d t t d -=--A ΦΦ,即)()(00t t dtt t d -=-A ΦΦ,所以)(0t t -Φ也是基本解矩阵。
由于线性齐次方程组任意两个基本解矩阵可以互相线性表示,故C ΦΦ)()(0t t t =-,由条件E Φ=)0(得,E ΦC Φ==)0()(0t ,即得 )(01t -=ΦC ,所以有)()()(001t t t t -=-ΦΦΦ。
评注:这是标准基本解矩阵的一个性质,即)exp()exp(])exp[(00A A A t t t t -=-。
5-4 试求下列方程的通解 1)22,sec πt πt x x <<-=+'',2)te x x 28=-'''。
解 1)i λλ±==+2,12,01,齐次方程的基本解组为t t x t t x sin )(,cos )(21== 所以1cos sin sin cos )](),([21=-=ttt t t x t x W ,取00=t ,利用常数变易公式ds s f s x s x W s x t x s x t x t φtt )()](),([)]()()()()(~0212112⎰-=可得原方程的特解为t t t t ds ss t s t t tcos ln cos sin cos 1)sin cos cos (sin )(~0⋅+=-=⎰ϕ,原方程的通解为t C t C t t t t x sin cos cos ln cos sin 21++⋅+=。
2)083=-λ,i 31,23,21±-==λλ,齐次方程基本解组为t et x t e t x e t x tt t 3sin )(,3cos)(,)(3221--===。
利用常数变易公式,原方程满足初始条件的特解为:ds s f s x s x s x W s x s x s x W t x t φk tk k )()](),(),([)](),(),([)()(~31321321∑⎰==,其中)](),(),([321s x s x s x W k 是在朗斯基行列式)](),(),([321s x s x s x W 中的第k 列代以()T 1,,0,0 后得到的行列式。
经计算可得),3cos33sin3()(),3cos33sin 3()(,3)(,312)(3221t t e t W t t e t W et W t W ttt+-=-===-可得原方程的特解为t et eetet tttt3sin 57633cos 1925241121)(~22---+-=ϕ,原方程的通解为 tttteeC et C t C x 22321121)3sin3cos (+++=-。
评注:此题主要是常数变易公式的应用。
常数变易公式表明线性非齐次方程的特解可以由对应齐次方程的基本解组的朗斯基行列式表示。
当然,此题中的2)用待定系数方法求特解会更简单。
5-5 给定方程)(78t f x x x =+'+''其中)(t f 在∞<≤t 0上连续,试利用常数变易公式,证明:1)如果)(t f 在∞<≤t 0上有界,则上面方程的每一个解在∞<≤t 0上有界; 2) 如果当∞→t 时0)(→t f ,则上面方程的每一解)(t φ,满足)(0)(∞→→t t φ当。
证 1)1,7,078212-=-==++λλλλ,齐次方程有基本解组t t e e 7,--tttt t eee eet W 87767)(------=--=。
利用常数变易公式:ds s f s x s x W s x t x s x t x t φtt )()](),([)]()()()()(~0212112⎰-=可得原方程的一个特解ds s f ee eeet φststst)()(61)(~7780------=⎰⎰⎰---=tsttstds s f eeds s f e e770)(61)(61 ,所以原方程的任一解为⎰⎰-----++=tsttstttds s f eeds s f e eeC eC t 0770721)(61)(61)(ϕ。