常微分方程第5章答案
常微分方程及其应用
第5章常微分方程及其应用习题5.22(6) X dy + ( 2xy-x +1 Hx = 0 , y(1)=0 .5.3可降阶微分方程及二阶常系数线性微分方程案例引入 求微分方程y" = 6x 的通解. 解 两边积分,得y ' = J6xdx = 3x 2+Gy = /(3x 2+0 dx =x 3t C r X +C 2所以,原方程的通解为y=x'+C i x+C 2,其中C i 、C 2为任意常数.5.3.1可降阶微分方程 1.形如y(n)= f(X)的微分方程特点:方程右端为已知函数f (x).(1)2X dx + ydy = 0 ; (2) x y'-yin y =0 ; (3) 22(xy +x)dx +(y-x y)dy =( );(4) y ‘ - 3xy =0 ;(5)c p x 2y —y =e ;(6) y' = ytanx+cosx .2.求下列各微分方程满足所给初始条件的特解: (1) 八e 23, y(0) =0;(2)dx - ydy = 0 , y(0)=1 + y 1+x (3) y ,- y =cosx , y(0) = 0; (4) y ‘ 一 y tan x = secx ,y(0) = 0 ;1 .求下列各微分方程的通解:1;y ・+Y=沁,y ⑴=1 ; x x两边再积分,得解法:对y(n)= f(x)连续积分n 次,即可得含有n 个任意常数的通解.2.形如y " = f(X, y )的微分方程 特点:方程右端不显含未知函数y .解法:令y = p(x),则y = p'(x).于是,原方程可化为 p' = f (x, p).这是关于p,p'的一阶微分方程•设其通解为 p(x)=®(x,C i ),即y'=®(x,C i ) •两边积分,即可得原方程通解y =代(X'CJdx +C 2,其中C 1、C 2为任意常数.3.形如yJf(y,y)的微分方程 特点:方程右端不显含自变量 x .= x+C 2,其中C i 、C2为任意常数. 屮(y,C i)例5-7 求微分方程y w = sin X - cosx 的通解.解 两边积分,得 y"=J(si n X-cosx) dx =-cosx-si nx + 2C i两边再积分,得 y = J(—COSX —sin ^2C 1 dx = -sin x + cosx +2C 1x + C ?解法:令/ = p(y),则y=如理*生=pd p.于是,原方程可化为dy dx dy dyPp = f (y, p).这是关于 p, p的一阶微分方程.设其通解为p(y)=屮(y,C i ),即吐=屮(y,C i ).分离变量,得dxdy 屮(y,C i )=dx•然后两边积分,即可得原方程通解dy2第三次积分,得 y = J(—sin X+C0SX+2C 1X +C 2 dx =cosx + sinx+Gx +C 2X + C 3所以,原方程的通解为 y =cosx+sin X+0x 2+C 2X+C 3,其中c i例5-8求微分方程xy”_y'=0的通解.解 令y ,= p(x),则y = p'(x).原方程可化为xp' — P = 0 , 是关于p, p'的一阶线性齐次微分方程•其通解为:y = f2C 1xd^C 1x^C 2,其中 C 1、C 2为任意常数.例5-9 求微分方程y"-一 y'xe 」的通解.X解 令y = p(x),则y " = p(X)•于是,原方程可化为P,P '的一阶线性非齐次微分方程.其通解为一dx IV一 -dxp(x)=e'x If xe e 'x dx + 2C 1=x ( fe^dx +2Ci )= xjr +2Ci )即y'=x(-ej +2C 1)•两边积分,即得原方程通解y = J x (-e " +2G dx = [(—xe * +2C 1xdx = f xd (e ^ )+C 1x 2= xe^ - fe^dx +C 1x^(^1)e^ +C 1x^C 2、C 2、C 3为常数.P(X)=2C ieRx =2C 1e ln^2C 1x ,即 /=2C 1x •两边积分, 即得原方程通解p = xe* •这是关于-In x , , c —dx +2C i其中C i、C2为任意常数.则y " = PP '( y ).于是,原方程可化为ypp' 一 P 2= 0 ,即例5-10 求微分方程2yy"—(y j =0 的通解.p'—丄p = 0 •这是关于yp, p'的一阶线性齐次微分方程.其通解为 fldylp(y) = C i e y = Ge ny = C i y ,即 y ' = G y .所以原方程通解为y=C 2e 卩"'=C z e&x ,其中C i 、C 2为任意常数.5.3.2二阶常系数齐次线性微分方程 定义5.4形如y" + py’ + qy =0 , p 、q 为常数(5-5)的微分方程,称为二阶常系数齐次线性微分方程1.二阶常系数齐次线性微分方程解的结构定理5.1如果函数y^x )和y 2(x )是方程(5-5 )的两个解,那么y =C i y i (x)+C 2y 2(x), G 、C 2为任意常数(5-6)也是方程(5-5)的解.(证明略)定理5.1表明,二阶常系数齐次线性微分方程的解具有叠加性.那么叠加起来的解y =C i y i (x ) + C 2y 2(x )就是通解吗?不一定.例如,设函数 y i (x )是方程(5-5)的一个解,则函数 y 2(x )=2y i (x )也是方程(5-5)的一个解.由定理 5.i 可知,y =C i y i (x)+2C 2y i (x) =(C i +2C 2)y i (x)是方程(5-5)的解•但C^2C 2=C 仍是一个任意常数, 所以y =(C i +2C 2)y i (x ) =Cy i (x )不是方程(5-5)的通解.那么在什么条件下才能保证y^Gyjx ) +C 2y 2(x )就是通解呢?定义5.5设y i (x )和y 2(x )是定义在某区间I 上的两个函数,如果存在两个不全为零 的常数k i 和k 2,使k i y i (x )+k 2y 2(x ) =0在区间I 上恒成立,则称函数 % (x )与y 2(x )在区间I 上线性相关,否则称线性无关.由定义5.5可知,判断函数y^x )与y 2(x )线性相关或线性无关的方法:当y^X )= 一鱼=常数时,y i (x )与y 2(x )线性相关•当 y i (x ) k 2y 2(x )线性无关.定理5.2 如果函数y i (x )和y 2(x )是方程(5-5)的两个线性无关的特解,那么 (5-6)是方程(5-5)的通解.(证明略)2.二阶常系数齐次线性微分方程的解法(5-5)的两个线性无关的特解.所以只要r 满足方程P 、q 为常数即当r 是方程(5-7)的根时,函数y =e rx就是方程(5-5)的解.定义5.6方程(5-7)称为方程(5-5)的特征方程.特征方程的根称为 特征根.y 2(x) y i (x)H 常数时,y i (x )与由上述关于解的结构分析可知,欲求方程(5-5)的通解,首先需讨论如何求出方程猜想方程(5-5)有形如y =e rx的解,其中r 为待定常数•将y=e rx代入该方程,得(e")” + 卩(「)’初(「)+ pre" r/2rx=(r 中 pr 中q)e =0 , rx由于e H0,(5-7)设r,、r2为特征方程(5-7)的两个特征根.根据特征根的不同情形,确定方程(5-5)的通解有以下三种情况:(1)若方程(5-7)有两个不相等的实根r i H「2,则y i = e rix和y^ e r2x是方程(5-5)的两个线性无关的特解,故方程(5-5)的通解为^C1e r1^C2e r2x,其中C1、C2为任、、八Nr、匸意吊数.(2)若方程(5-7)有两个相等实根r, =「2 = r =-卫则仅得到一个特解y^e rx2利用常数变易法可得到与y1 =e rx线性无关的另一个特解y= xe rx,故方程(5-5)的通解为y =C1e r x +C2xe r x,其中C,、C2为任意常数.(3)若方程(5-7)有一对共轭复根» =a +i P与『2 =a - i p,则y^e^*^x和y2 是方程(5-5)的两个复数特解.为便于在实数范围内讨论问题,在此基础上可找到两个线性无关的实数特解e^cos Px和e^x sin Px .故方程(5-5 )的通解为y =e g C i cos P x +C2 sin P x),其中C,、C2 为任意常数.由定理5.1可知,以上两个函数e Q X cos P x和^x si n P x均为方程(5-5)的解,且它们线性无关.上述依据特征根的不同情形来求二阶常系数齐次线性微分方程通解的方法, 称为特征根法.一般步骤:第一步写出所给微分方程的特征方程;第二步求出特征根;第三步根据特征根的三种不同情形,写出通解. (特征根与通解的关系参见表5-1)表5-1特征根与通解的关系特征方程r2+pr +q =0的两个根r, , q 微分方程y"+ py' + qy=O的通解例5-11 求微分方程y“-2y'-3y=0的通解.解该方程的特征方程r — 2r —3 = 0的特征根为匚=—1, D = 3 ( 口H j).所以,方程的通解为y - De- +C2e3x.例5-12 求微分方程y”+2y' + y=0满足初始条件y(0)=0, y'(0) =1的特解.解该方程的特征方程r2 +2r +1 =0的特征根为R =「2 = -1.所以方程的通解为y =(C1 +C2X)e」上式对x求导,得: y,=C2e 一-(G +C2x)e 一将y(0) =0, y'(0) =1代入上两式,解得G =0 , C2 =1 .因此,所求特解为y = xe^ 例5-13 求微分方程y”—2y,+5y =0 的通解.解该方程的特征方程r2 -2r +5 =0的特征根为匚=1 +2i ,心=1 -2i .所以,方程的通解为y = e x(C1 cos2x +C2Sin2x).5.3.3二阶常系数非齐次线性微分方程定义5.7形如y" + p yyqy = f(x) , p 、q 为常数的微分方程,称为二阶常系数非齐次线性微分方程1.二阶常系数非齐次线性微分方程解的结构如果函数y^x)是方程(5-8)的一个特解,Y(x)是该方程所对应的线性齐 其中Q n (x)是与P n (x)同次待定多项式.(5-8)定理5.3 次方程(5-5)的通解,那么y =Y(X)+y "x)(5-9)是方程(5-8) 的通解.定理 5.4 如果函数y 1(x)是方程y" + Py'+qy = fjx)的特解,函数y 2(x)是方程y +py +qy =f 2(x)的特解,那么/ = y i”(x) +y 2(x)(5-10)就是方程y" + py' + qy = f 1(x) + f 2(x)的特解.2.二阶常系数非齐次线性微分方程的解法二阶常系数齐次线性微分方程的通解问题已经解决,根据定理 5.3,求二阶常系数非齐次线性微分方程的通解的关键在于求其自身的一个特解.以下介绍当自由项 f(x)为几类特殊函数时求特解的方法: (1) f(X)=P n (x)e/x, P n (x)是x 的n 次多项式,A 是常数 微分方程的特解可设为沖k. /Xy =x Q n (x)e",不是特征根时,k =/是单特征根时,k = 1(2) f (X) =P n(x)cos©x (或P n(x)s in ©x), P n(x)是x 的n 次多项式,©是常数微分方程的特解可设为… ks /、 C /、. ■, 血i非特征根时,k=0y=X[Qn(X)沁x+R n(X)s gX] '{勿j是特征根时,心其中Q n(x)和R n(x)是与P n(x)同次待定多项式.(3) f (X) =ecos^x (或e赵sin w x), A与©均为常数微分方程的特解可设为k7 几+切)非特征根时,k=0y =x e那[Acoseox+Bsin^x], —,iA+©i是特征根时,k= 1(4)当f (x)为上述任意两类函数之和时,根据定理5.4处理即可.例5-14 求微分方程y"-2y' = 3x+1的通解.解方程y"—2y'=0的特征方程r -2r=0的特征根为r, = 2 , “=0 •于是方程y"-2y'=0的通解为又因为P n(x) =3x +1,几=0是单特征根,所以原方程的特解可设为y* = xQ n(X)= x( Ax + B)3 5代入原方程,解得 A = ——, B = - 一.故原方程的通解为4 4C 2x 土小3 2 5^C1e 中C2 一一X —-x .4 4例5-15 求微分方程y"中y ' + y=3e2x的一个特解.y =C i cosx +C 2 sin xy * = x( Acosx + Bsin x)解 方程y " + y + y =0的特征方程 r 2+r +1=0的特征根为r 1r2 = 一一1 "*^i . f(x^3e 2x , Z =2非特征根,所以原方程的特解可设为2 2代入原方程,解得T .故所求特解为几尹例5-16 求微分方程y" + 3/ +2y=xe 上X 的一个特解. 解 方程y ”+3y ' + 2y =0的特征方程r2+3r +2=0的特征根为口 = 一2 ,「2 = —1 . f(X)=xe'x, P n (x) =x ,扎=-2是单特征根,所以原方程的特解可设为/ =x(Ax +B)e41 X 2代入原方程,解得A = —— ,B = —1.故所求特解为y*=x(—— —1)e2 2例5-17 求微分方程y " + y = Sin X 的通解.=0的特征方程r 2 +1 =0的特征根为 山=i ,「2 =-i .于是方程¥丄y +y =0的通解为又因为f(X)=sin X , 几+ Q i =i 是特征根,所以原方程的特解可设为代入原方程,解得 A=—1, B=0 •故原方程的通解为 21y = G cosx +C 2Sin x —一 xcosx .2例5-18 求微分方程y" + y=xcos2x 的一个特解.解 方程y" + y =0的特征方程r 2+1=0的特征根为r , = i , q = —i .f (X )= XCOS2X ,乙+ © i = 2i 不是特征根,所以原方程的特解可设为y * =(Ax + B) cos2x + (Cx + D)sin 2x1 4 代入原方程,解得 A = — — , B=0 , C=0, D=—.故所求特解为39* 1 4y =——xcos2x + — sin 2x .3 9例5-19 求微分方程y "+3y '-y =e X cos2x 的一个特解.解 方程y "+3y ' —y=0的特征方程r2+3r —1=0的特征根为r 1设为y * =eX ( Acos2x + Bsin 2x)解 方程y -2y ' + y =0的特征方程r 2 -2r +1=0的特征根为n = q = 1 .1f 1(x ) =—e X , f 2(x ) =sinx , =1是二重特征根,=i 不是特征根,所以两个分解r 23 V 13f (x ) =e x cos2x ,几+ eo i =i+2i 不是特征根,所以原方程的特解可代入原方程,例 5-201 10得 A =———,B=——.故所求特解为101101 y * =e X ( -^^cos2x + n 2x).101 1011求微分方程y "—2y' + y=—e +sinx 的一个特解. cos2x +2方程的特解可分别设为y r = Ax 2e x 与 y2 = B cosx +Csin x1分别代入两个分解方程,解得A= — ,B= — ,C=0 .故所求特解为4X 12 xy = — x e4习题1.求下列各微分方程的通解:5. 求下列各微分方程满足所给初始条件的特解:(1) y" —4y' + 3y =0, y(0) =6, y'(0) =10 ; (2)y" —4y' + 4y=0, y(0) =1 , y'(0) = 4 .1 +- cosx . 25.3(1) y " = X +sin x ;(2)利xy =(3) xy" + y'=0 ;(4)H 1, xy — 一 y = xe ;(5)八1 +(y )2;+ _2 1 2 (y) =0. -y2.求下列各微分方程满足所给初始条件的特解:⑴ yJe 2x, y(1) =y(1) =y"(1) =0 ;2(2) y” —3(y) =0 , y(0) =0 , y(0) =—3.判断下列各函数组是线性相关还是线性无关: (1) x与 X 2; (2) e2x与 6e 2x; X(3) x与xe ; / 、 x. x .(4) e cosx 与e sin4. 求下列各微分方程的通解:(1) y” —y'=0;(2) y “+ 4y=0 ; (3) y” —10y'+25y =0 ;(4) y "+y ' + y= 0 .6. 求下列各微分方程的一个特解:7. 求下列各微分方程的通解:&求下列各微分方程满足所给初始条件的特解:(1) yJ3yq2y=5,y(0) =1,y'(0) =2 ; (2) y"-y =4xe X , y(0)=0, "(0)=1.5.4微分方程应用举例微分方程在实践中有着广泛的应用.在实际应用中,常常需要应用微分方程寻求实际 问题中的未知函数.而要建立微分方程,除了需要数学知识外,往往还需要许多专业方面 的知识•本节通过举例来介绍微分方程在几何学、电工学及力学方面的一些简单应用.例5-21曲线L 上点M(x,y)处的法线与x 轴的交点为 N ,且线段 MN 被y 轴平 分.求曲线L 的方程.解 如图5-2,设曲线的方程为 y =y(x).先建立法线 MN 的方程.设法线上的动点坐标为(X ,Y),由于法线MN 的的斜率为 1k 法=-二,于是法线MN 的方程为y'1Y-y = -rXy(1) y"_2y ,—3y =3x +1 ; yN_4y ,+ 4y =e 2x ; (3)y "-2y ,+ 2y =e 」sinx ;y "+ 4y = X +1 +sin X .(1) y "-2y y = X 2y" + 2y'-3y =e X ;(3) y "+y =e x +cosx ;yH-y ,-2y = x + cos2x .-x)又因为线段 MN 被y 轴平分,从而MN 与y 轴交点坐标为p(0丄),代入上式,得,2^-y =—丄(0-X),即 yy' = —2x例5-22 设有一 RC 电路如图5-3所示,电阻 R = 1O 0,电容C =0.1F ,电源电 压U =10si nt(V),开关K 闭合前,电容电压u c =0,求开关K 闭合后电容电压随时图5-3解 设开关K 闭合后电路中的电流为i(t),电容极板上的电荷为q(t),则有用分离变量法解得X 2其中c 为任意正数.CRc . dq d(Cu c) c du e q = cu c, i =— = ------------------- = c ---dt dt dt根据回路电压定律:电容电压与电阻电压之和等于电源电压,即U c + Ri = U,于是有u^Rcdu c- .将R = 10 , C = 0.1 , U =10sint代入,得uC +u c =10sint .又因为开关K闭合前, 电容电压u c =0,即u c (0) =0 .从而问题转化为初值问题:(u C +U c =10si nti u c(0) =0用通解公式求得通解u c = Ae,+5(s i ti-cos)将初始条件u c (0) =0代入通解,求得A=5•所以,所求特解为U c =5e 丄+5(s in t-cost)此即为所求规律U c(t)的表达式.例5-23 设跳伞员开始跳伞后所受的空气阻力与其下落的速度成正比(比例系数为常数k > 0),起跳时的速度为0 .求跳伞员下落的速度与时间之间的函数关系.解这是一个运动问题,可利用牛顿第二定律 F =ma建立微分方程.设跳伞员下落的速度与时间之间的函数关系为v =v(t),则加速度a =v'(t).由于跳伞员在下落过程中所受外力只有重力和空气阻力,于是有 F =mg -kv,由牛顿第二定律F = ma可得速度V =v(t)应满足的微分方程为mg -kv =mv .又因为起跳时的速度为0 ,即其初始条件为v(0) =0 •所以,这个运动问题可化为初值问题:何g - kv = mv’ 卜(0) = 0上tmg -kv =Ce m .将初始条件为 v(0) =0代入通解,解得kv^^O-e 韦t) , 0<t<T (T 为降落伞着地时间),此 k即为所求函数关系.在时刻t =0时,测得其温度为150C ,10分钟后测得温度为100£ .已知牛顿冷却定律: 物体冷却速率与物体和介质的温差成正比. 求物体的温度与时间的函数关系,并计算20分钟后该物体的温度.解 设物体的温度与时间的函数关系为T =T(t).因为热量总是从温度高的物体向温度低的物体传导,从而物体随时间增加而逐渐冷却,所以冷却速率(温度的变化速度)t =20代入,得 T(20) =24 +12604051 过0疋 64(°C).例5-25弹簧振动问题.设有一弹簧上端固定, 下端挂着一个质量为 m 的物体.簧处于平衡位置时,物体所受的重力与弹簧恢复力大小相等,方向相反•设给物体一个初 始位移X o ,初速度V o ,则物体便在其平衡位置附近上下振动•已知阻力与其速度成正比,用分离变量法求出通解为C =mg .因此,所求特解为例5-24 物体冷却过程. 将某高温物体置于空气中冷却,假定空气温度恒为 249 , T (t) <0 ,而物体和空气的 温差恒为正.所以,根据牛顿冷却定律可得¥十(一24).又因为在时刻而问题转化为初值问题:t=0时,测得其温度为1509,即有T(0)=150 .从其中k > 0为比例常数.T(0) =150用分离变量法或通解公式解得T =24+1260山.将T(10) =100代入,求得k.^^^0051.故物体的温度与时间的函数关系为T =24+126eq 051•将当弹求振动过程中位移 x 的变化规律.////////解 建立坐标系如图5-4所示,平衡位置为原点. 位移x 是时间t 的函数X = x(t).物体在 振动过程中受到弹簧恢复力f 与阻力R 的作用•由虎克定律,有 f「1,2 = -n ±v n 2 -豹2.具体情况讨论如下:弹性系数,负号表示弹簧恢复力与位移方向相反;R = -P v ,其中 卩> 0为比例系数(或称阻尼系数),负号表示阻力与速度方向相反.根据牛顿第二定律F = ma ,可得ma = —kx - 出.又因为 a = x"(t),V =x'(t),记 2n =巴八2m所以上述弹簧振动问题化为初值问题:亍2|d X , _ dx ." --- +2n —— dt dt1x(0) =X 0,x'(0)这是一个二阶常系数齐次线性方程, 其特征方程为r 2 +2nr 2=0,特征根为=—kx ,其中k >0为(1) 大阻尼情形,即n X .这时, 特征根是两个不相等实根,所以方程的通解为x = GeG 』2^2^ +c e -(E n 2看)t(2) 临界阻力情形,即 n =©O图5-4停止,称为弹簧的阻尼自由振动.习题5.41 •设过点(1,1)的曲线L 上任意点M(x,y)处的切线分别与x 轴、y 轴交于点A 、B , 且线段AB 被点M 平分.求曲线L 的方程.2 •在如图5-5所示的RC 电路中,已知开关S 闭合前,电容上没有电荷,电容两端C ,电源电压为E •把开关S 合上,电源对电容充电,电图5-5 3.将温度为100心的沸水注入杯中,放在室温为20°C 的环境中自然冷却,5min 后 测得温度为60乜.求水温与时间的函数关系, 并计算水温自100七降至30°C 所需时间.4. 设有一弹簧上端固定,下端挂着一个质量这时,特征根r ,=『2 = -n , 所以方程的通解为X = (C i + C 2t)e 』t(3)大阻尼情形,即n X这时,特征根是一对共轭复根 「1,2 = -n ±晶2 _n 2i ,所以方程的通解为X =e 』t (C 1 cos JJ -n 21 + C2 sin上述三种情形中的任意常数均可由初始条件确定. 这类振动问题均会因阻尼的作用而容电压U c 逐渐升高.求电容电压U c 随时间t 变化的规律.电压为零,电阻为 R ,电容为为0.025kg 的物体.先将物体用手拉到离平衡位置0.04m 处,然后放手,让物体自由振动.若物体所受的阻力大小与运动速度成 弹簧的弹性系数 k=0.625N/m ,阻尼系数4 =0.2N ‘s/m .求物体 的运动规律. 知识拓展:马尔萨斯(Malthus )模型Malthus )模型是最简单的生态学模型.给定一个种群,我们的目的是确 定种群的数量是如何随着时间发展变化的.为此,我们作出如下假设:模型假设:1. 初始种群规模已知 x (0) =X o ,种群数量非常大,世代互相重叠,因此种群的数量 可以看作是连续变化的;2. 种群在空间分布均匀,没有迁入和迁出(或迁入和迁出平衡)3.种群的出生率和死亡率为常数, 即不区分种群个体的大小、年龄、性别等;4.环境资源是无限的.确定变量和参数: 为把问题转化为数学问题,我们首先确定建模中所需变量和参数:t :时间(自变量),x (t ) : t 时刻的种群密度,b :瞬时出生率,d :瞬时死亡率.模型的建立与求解:考察时间段[t,t +A t ](不失一般性,设 M>0),由物质平衡原理,在此时间段内种 群的数量满足:t + A t 时刻种群数量-t 时刻种群数量 =也t 内新出生个体数- A t 内死亡个体数,正比,方向相反, 马尔萨斯(x(t +A t) -x(t) =bx(t)At -dx(t)At=(b-d)x(t):=Ax(t) dt满足初始条件x(0) =X o 的解为X(t) ^00(5 = X 。
《常微分方程》第五章练习题
x
y
C1
e3t 2e3t
C2
et 2et
3、满足初值条件的解为
~
(t )
et e t
4、方程组的通解为
x y
C1e2t
4 5
C2e7t
1 1
。
4
5、所求基解矩阵为 (2 e
3t
3)e
3t
e 3t (2 3)r
3t .
6、 (t )
e3t [E
t(A
3E)]
A1 (t)
A2 (t)
,t
(a,b) .
部分参考答案 一、填空题
1、 (t) (t)C
2、(t) exp[(t t0 )A]
t t0
exp[(t s)A] f (s)ds
3、必要
t t0
1 (s) f
(s)ds
三、计算题
1、
A
4 3
3
4
2、原方程组的通解为
x ' Ax ce mt 有一解形如(t) pemt ,其中 c , p 是常数向量.
3
4、证明:如果 φ(t) 是方程组 x Ax 满足初始条件 φ(t0 ) η 的解,那么
φ(t) [exp A(t t0 )]η 。
5、证明:如果 Φ(t),Ψ (t) 在区间 a t b 上是 n 阶线性方程组
1、向量
X1
(t)
2et 0
,
X
2
(t)
t 2et et
的伏朗斯基行列式
W (t) =(
).
A 、0 ; B 、 tet ; C 、2 e t ; D 、2 e2t .
2、有关矩阵指数 exp A 的性质,以下说法正确的是( )
常微分方程智慧树知到答案章节测试2023年内蒙古民族大学
第一章测试1.下列方程不是常微分方程的是()A: .B: .C: .D: .答案:D2.下面方程中不是线性微分方程的是()A: .B:C: .D: .答案:D3.下列微分方程不是驻定的是()A:B:C:D:答案:C4.下面是微分方程的特解的是()A: .B: .C: .D:答案:D5.微分方程的阶数是().A:1;B:3;C:4.D:2;答案:B6.下列方程中的线性微分方程是().A: .B: ;C: ;D: ;答案:D第二章测试1.下列微分方程中,可分离变量的是( )。
A:B:C:D:答案:D2.下列函数中,哪个是微分方程的解( )。
A:y=2xB:y=-xC:y=-2xD:y=x2答案:D3.微分方程的一个特解是( )。
A:B:C:D:答案:C4.满足的特解是( )。
A:B:C:D:答案:C5.方程的通解是( )。
A:B:C:D:答案:B第三章测试1.利用唯一性充分条件,在平面上微分方程有唯一解的区域是()A:B:C:D: .答案:A2.微分方程的第二次近似解是()A:B:C:D:答案:D3.按存在唯一解定理,微分方程第一次近似解在区域中的误差估计是()A:0.375B:0.625C:0.125D:0.325答案:A4.方程存在唯一解的区域是()A:除了外均存在唯一解B:除了外均存在唯一解C:除了外均存在唯一解D:除了外均存在唯一解答案:D5.方程存在唯一解的区域是以下选项中的()A:B:C:D:答案:C6.方程的第二次近似解在解的存在区间的误差估计是()A:B: .C: .D: .答案:A7.李普希兹条件是保证一阶微分方程初值问题解唯一的()A:充分条件B:不充分不必要条件C:充要条件D:必要条件答案:D8.方程过点共有()A:二个解B:无数个解C:一个解D:没有解答案:B9.方程组的任何一个解的图象是()A: 维的B: 维的C: 维的D: 维的答案:B10.连续是保证方程初值唯一的()A:不充分不必要条件B:必要条件C:充分条件D:充要条件答案:C11.阶线性非齐次微分方程的所有解()A:不能构成一个线性空间B:构成一个维线性空间C:构成一个维线性空间D:构成一个线性空间答案:A12.利普希茨条件是一阶微分方程存在唯一解的()条件.A:充分条件B:必要条件C:充分必要条件D:既非充分也非必要条件.答案:A13.函数对是否满足李普希兹条件()A:满足B:可能不满足C:不满足D:可能满足答案:A14.如果存在常数使得不等式()对于所有都成立,称为利普希兹常数,函数称为在上关于满足利普希兹条件。
常微分方程(第三版) 王高雄等编 高等教育出版社 课后习题答案
1常微分方程习题答案2.11.xy dx dy2=,并求满足初始条件:x=0,y=1的特解.解:对原式进行变量分离得。
故它的特解为代入得把即两边同时积分得:e e xx y c y x x c y c y xdx dy y22,11,0,ln ,212=====+==,0)1(.22=++dy x dx y 并求满足初始条件:x=0,y=1的特解.解:对原式进行变量分离得:。
故特解是时,代入式子得。
当时显然也是原方程的解当即时,两边同时积分得;当xy c y x y x c y c y x y dy dx x y++=====++=+=+≠=+-1ln 11,11,001ln 1,11ln 0,11123.yxy dx dy x y 321++=解:原式可化为:x x y xxyxyx yyxyc c c c x dx x dy y yx ydxdy2222222232232)1(1)1)(1(),0(ln 1ln 21ln 1ln 2111,0111=++=++≠++-=++=+≠+∙+=+)故原方程的解为(即两边积分得故分离变量得显然.0;0;ln ,ln ,ln ln 0110000)1()1(4===-==-+=-++=-=+≠===-++x y c y x xy c y x xy c y y x x dy y y dx x x xy x y xdy y ydx x 故原方程的解为即两边积分时,变量分离是方程的解,当或解:由:10ln 1ln ln 1ln 1,0ln 0)ln (ln :931:8.cos ln sin ln 07ln sgn arcsin ln sgn arcsin 1sgn 11,)1(,,,6ln )1ln(21111,11,,,0)()(:53322222222222c dxdy dx dy xycy ud uudx x x y u dx xydy x y ydx dy y x x c dy yy yydxdy c x y tgxdx ctgydy ctgxdy tgydx cx x xycx x u dxx x du xdxdudxdux u dx dy ux y u x y y dx dy xc x arctgu dx x du u u u dx du x u dxdu xu dx dy ux y u x y x y x y dx dy dx x y dy x y ee ee ee eexy uu xy x uu xyxyyx xx+===+=+-===-∙-=--+-=-=+-===-=+∙=+∙=∙=--=+===-+=+-=++=++-++=++===+-==-++-+--两边积分解:变量分离:。
常微分方程知到章节答案智慧树2023年齐鲁师范学院
常微分方程知到章节测试答案智慧树2023年最新齐鲁师范学院第一章测试1.二阶微分方程的含有两个任意常数的解一定是通解。
()参考答案:错2.满足初值条件的解称为是微分方程的特解。
()参考答案:对3.一阶微分方程的通解表示平面上的一条曲线。
( )参考答案:错4.不是线性微分方程的方程一定是非线性微分方程。
( )参考答案:对5.函数为任意常数是方程的通解。
( )参考答案:对第二章测试1.一阶非齐次线性微分方程的任意两个解之差必为相应的齐次线性微分方程的解。
()参考答案:对2.微分方程()参考答案:二阶线性微分方程3.微分方程的满足的特解为()参考答案:4.微分方程的通解为()参考答案:5.若一阶微分方程有积分因子,则积分因子一定是唯一的。
()参考答案:错第三章测试1.所有的微分方程都可以通过初等积分法求得其通解。
()参考答案:错2.要求得一阶微分方程的特解,应该给定一个初值条件。
()参考答案:对3.李普希兹条件是一阶微分方程初值问题解存在唯一的充要条件。
()参考答案:错4.存在唯一性定理中解的存在区间是唯一的。
()参考答案:错5.微分方程初值问题的解只要存在就一定唯一。
()参考答案:错第四章测试1.若函数在区间上线性相关,则在上它们的伏朗斯基行列式。
()参考答案:错2.如果方程的解在区间上线性无关,则在这个区间的任何点上都不等于零,即()参考答案:对3.由n阶齐线性方程的n个解构成的伏朗斯基行列式或者恒等于零。
( )参考答案:对4.n阶齐线性方程可以有n+1个线性无关的解。
()参考答案:错5.是方程的通解。
()参考答案:对第五章测试1.如果矩阵,维列向量是可微的,则()参考答案:对2.向量是初值问题在区间上的解。
()参考答案:对3.设是矩阵,则。
()参考答案:对4.如果向量函数在区间线性相关,则它们的伏朗斯基行列式,。
( )参考答案:对5.如果,在区间上是的两个基解矩阵,那么,存在一个非奇异常数矩阵,使得在区间上。
第五章 常微分方程初值问题数值解法
则有
yn 1 yn hf ( xn , yn )
( 5.2 ) Euler格式
例5.1 用Euler格式解初值问题
2x y y y y (0) 1
取步长h=0.1.
(0 x 1)
Euler格式的具体形式为
y n 1 y n hf ( x n , y n ) 2 xn yn 0.1( yn ) yn 0.2 xn 1.1 yn yn
计算公式的精度 常以Taylor展开为工具来分析计算公式的精度. 为简化分析,假定yn是准确的,即在 yn y ( xn ) 的前提下估计误差 y ( xn 1 ) yn 1 Euler格式的局部截断误差 由 从而 局部截断误差
f ( xn , yn ) f ( xn , y ( xn )) y '( xn ) y ( xn 1 ) yn 1 y ( xn 1 ) ( yn hf ( xn , yn )) y ( xn 1 ) y ( xn ) hy '( xn )
y ( xn ), y ( xn 1 ), 的近似值 y1 , y2 , , yn , yn 1 ,
相邻两个节点的间距 h xi 1 xi 称为步长,步 长可以相等,也可以不等.本章总是假定h为定数, 称为定步长,这时节点可表示为
xn x0 nh , n 0,1, 2,
由f ( xn 1 , yn 1 ) f ( xn 1 , y ( xn 1 )) f y ( xn 1 , )( yn 1 y ( xn 1 )) f ( xn 1 , y ( xn 1 )) y '( xn 1 )(在xn点Taylor展开) h2 y '( xn ) hy ''( xn ) y '''( xn ) ... 2 3 2 h h 因此yn 1 y ( xn ) hy '( xn ) y ''( xn ) y '''( xn ) 2 4 h f y ( xn 1 , )( yn 1 y ( xn 1 )) 2 h2 h3 y ( xn 1 ) y ( xn ) hy '( xn ) y ''( xn ) y '''( xn ) 2 3!
常微分方程第5章答案
习题1.给定方程组x = x x= (*)a)试验证u(t)= ,v(t)= 分别是方程组(*)的满足初始条件u(0)= , v(0)= 的解.b)试验证w(t)=c u(t)+c v(t)是方程组(*)的满足初始条件w(0)= 的解,其中是任意常数.解:a) u(0)= =u (t)= = u(t)又 v(0)= =v (t)= = = v(t)因此 u(t),v(t)分别是给定初值问题的解.&b) w(0)= u(0)+ u(0)= + =w (t)= u (t)+ v (t)= +=== w(t)因此 w(t)是给定方程初值问题的解.2. 将下面的初值问题化为与之等价的一阶方程组的初值问题:a) x +2x +7tx=e ,x(1)=7, x (1)=-2b) x +x=te ,x(0)=1, x (0)=-1,x (0)=2,x (0)=0;c)x(0)=1, x (0)=0,y(0)=0,y (0)=1解:a)令 x =x, x = x , 得即又 x =x(1)=7 x (1)= x (1)=-2于是把原初值问题化成了与之等价的一阶方程的初值问题:x = x(1)=其中 x= .b) 令=x ===则得:/且 (0)=x(0)=1, = (0)=-1, (0)= (0)=2,(0)= (0)=0于是把原初值问题化成了与之等价的一阶方程的初值问题:= x(0)= , 其中 x= .c) 令w =x, w =,w =y,w =y ,则原初值问题可化为:且即 ww(0)= 其中 w=3. 试用逐步逼近法求方程组】= x x=满足初始条件x(0)=的第三次近似解.解:\0241201 杨素玲习题02412—02 02412—031.试验证 =是方程组x = x,x= ,在任何不包含原点的区间a 上的基解矩阵。
解:令的第一列为 (t)= ,这时 (t)= = (t)故 (t)是一个解。
第5章常微分方程数值解法
2hyxn
2h3 3!
y
yxn1
yxn1 2hyxn
h3 3
y ②
将①、②两式相减:
y
xn1
h3 yn1 3
y
——两步法局部截断误差
18
第19页/共24页
2024年8月7日
(2)梯形公式
yn1
yn
h 2
f
xn , yn
f
xn1 , yn1
yxn
hf
2
xn ,
yn1 yxn1 2hyxn ①
2024年8月7日
17
第18页/共24页
2024年8月7日
将函数用泰勒级数展开:( h 较小, 相差不大)
yxn1
yxn
h 1!
yxn
h2 2!
yxn
h3 3!
y
yxn1
yxn
h
1!
yxn
h2
2!
yxn
h3
3!
y
yxn1
yxn1
yxn
f
xn1 , yxn1
yn1
yxn
h 2
y
xn
yxn1 ①
将函数用泰勒级数展开:
yxn1
yxn
h 1!
yxn
h2 2!
yxn
h3 3!
y ②
yxn1
yxn
h 1!
yxn
h2 2!
y ③
( h 较小, 相差不大)
19
第20页/共24页
2024年8月7日
①、②两式相减,并代入③式:
(图示表示梯形法计算结果)
第5章_常微分方程数值解法
(5.2.6)
由于方程关于 uk +1 是隐式形式,所以式(5.2.6)称为隐式 Euler 公式。前面显式和隐式 Euler 公式在计
u '(tk ) ≈
得到的递推公式:
u (tk +1 ) − u (tk −1 ) 2h
(5.2.7)
uk +1 ≈ uk −1 + 2hf (tk , uk )
在计算 uk +1 时,需要用到前两步结果 uk −1 , uk ,称为两步法公式。 (2)积分近似方法 将(5.2.1)式的微分方程写成 du = f (t , u )dt ,在区间 [tk , tk +1 ] 上积分,有:
5.2.2 Runge-Kutta 方法 Euler 方法比较简单,但它的收敛阶数低。可以利用 Taylor 展开式构造高阶的单步方法。Euler 公式 可以看成是由一阶 Taylor 展开式得到的,所以应用高阶 Taylor 展开就可以得到高阶单步法。例如:将 u (tk +1 ) 在 tk 处作 q 阶 Taylor 展开:
dy = a − by (t ) dt
是一阶常微分方程,而
2 ∂ 2 u ( x, t ) 2 ∂ u ( x, t ) a = ∂t 2 ∂x 2
(5.1.1)
(5.1.2)
是二阶偏微分方程。 所有使微分方程成为等式的函数,都是微分方程的解;在 n 阶微分方程中,将微分方程的含有 n 个任 意常数的解称为该微分方程的通解。为确定微分方程通解中的任意常数而需要的条件称为定解条件;定解 条件可以分为初始条件和边界条件两类。由微分方程和定解条件一起构成的问题称为微分方程定解问题。 根据定解条件的不同,常微分方程分为初值问题和边值问题;若定解条件是描述函数在一点(或初始 点)处状态的,则称为初值问题,一阶常微分方程初值问题的一般形式为:
常微分课后答案第五章
常微分课后答案第五章第五章 线性微分方程组§5.1 存在唯一性定理习题5.11.给定方程组x x ⎥⎦⎤⎢⎣⎡-='0110,⎪⎪⎭⎫ ⎝⎛=21x x x . (*))a 试验证⎪⎪⎭⎫ ⎝⎛-=t t t u sin cos )(,⎪⎪⎭⎫ ⎝⎛=t t t v cos sin )(分别是方程组(*)的满足初始条件⎪⎪⎭⎫ ⎝⎛=01)0(u ,⎪⎪⎭⎫⎝⎛=10)0(v 的解;)b 试验证)()()(21t v c t u c t w +=是方程组(*)的满足初始条件⎪⎪⎭⎫⎝⎛=21)(c c t w 的解,其中21,c c 是任意常数.证明)a ⎪⎪⎭⎫ ⎝⎛=01)0(u ,⎪⎪⎭⎫ ⎝⎛=10)0(v 显然.)(0110sin cos 0110cos sin )(t u t t t t t u ⎪⎪⎭⎫⎝⎛-=⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛--=',)(0110cos sin 0110sin cos )(t v t t t t t v ⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛-=',所以⎪⎪⎭⎫ ⎝⎛-=t t t u sin cos )(,⎪⎪⎭⎫⎝⎛=t t t v cos sin )(分别是方程组(*)的满足初始条件⎪⎪⎭⎫ ⎝⎛=01)0(u ,⎪⎪⎭⎫ ⎝⎛=10)0(v 的解.)b ⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛=+=2121211001)0()0()0(c c c c v c u c w ,又)(0110)(0110)()()(2121t v c t u c t v c t u c t w ⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛-='+'=')(0110))()((011021t w t v c t u c ⎪⎪⎭⎫⎝⎛-=+⎪⎪⎭⎫ ⎝⎛-=,所以)()()(21t v c t u c t w +=是方程组(*)的满足初始条件⎪⎪⎭⎫⎝⎛=21)(c c t w 的解,其中21,c c 是任意常数.2.将下面的初值问题化为与之等价的一阶方程组的初值问题:)a t e tx x x -=+'+''72,7)1(=x ,2)1(-='x ;)b tte x x =+)4(,1)0(=x ,1)0(-='x ,2)0(=''x ,0)0(='''x ;)c ⎩⎨⎧=-'+-''=+-'+''tx y y y e y x y x t cos 15132,675,1)0(=x ,0)0(='x ,0)0(=y ,1)0(='y .(提示:令y w y w x w x w '=='==4321,,,)解 )a 设x x x x '==21,,则21x x x ='=',te tx xx x -+--=''='12272,即与该初值问题等价的一阶方程组的初值问题为⎪⎩⎪⎨⎧-==+--='='-.2)1(,7)1(,27,2121221x x e x tx x x x t)b 设x x x x x x x x'''=''='==4321,,,,则21x x x ='=',32x x x =''=',43x x x ='''=',tte xx +-='14,则得等价的一阶方程组的初值问题为⎪⎪⎩⎪⎪⎨⎧+-='='='='tte x x x x x x x x 14433221,,,,⎪⎪⎪⎪⎪⎭⎫⎝⎛-=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=0211)0()0()0()0()0(4321x x x x x .)c 令y w y w x w x w'=='==4321,,,,有⎪⎪⎩⎪⎪⎨⎧+-+='='+--='='tw w w w w w e w w w w w w t cos 13215,,567,431443431221 ,⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=1001)0()0()0()0()0(4321w w w w w ,为与原初值问题等价的一阶方程组的初值问题. 3.试用逐步逼近法求方程组xx ⎪⎪⎭⎫⎝⎛-='0110,⎪⎪⎭⎫⎝⎛=21x x x满足初始条件⎪⎪⎭⎫ ⎝⎛=10)0(x 的第三次近似解.解 ⎪⎪⎭⎫⎝⎛=10)(0t ϕ,⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫ ⎝⎛=⎰110011010)(01t ds t tϕ, ⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫⎝⎛-+⎪⎪⎭⎫⎝⎛=⎰2210211011010)(t t ds s t tϕ,第三次近似解为 ⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛=⎰2213610221*********)(t t t ds s s t t ϕ.§5.2 线性微分方程组的一般理论习题5.21.试验证⎥⎦⎤⎢⎣⎡=Φ12)(2t t t t是方程组x t tx ⎥⎥⎦⎤⎢⎢⎣⎡-='22102,⎥⎦⎤⎢⎣⎡=21x x x在任何不包含原点的区间b t a ≤≤上的基解矩阵. 证明 设⎪⎪⎭⎫⎝⎛=t t t 2)(21ϕ,⎪⎪⎭⎫ ⎝⎛=1)(2t t ϕ,则由于)(22102221022)(12221t t t t t t t t t ϕϕ⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛=',)(22101221001)(2222t t t t t t t ϕϕ⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛=',所以)(,)(21t t ϕϕ都是方程组的解,因而[])()()(21t t t ϕϕ=Φ是所给方程组的解矩阵.又由于在任何不包含原点的区间],[b a 上,0)(det 2≠-=Φt t (],[b a t ∈),故)(t Φ是所给方程组的基解矩阵. 2.考虑方程组xt A x )(=', (5.15)其中)(t A 是区间b t a ≤≤上的连续n n ⨯矩阵,它的元素为)(t a ij,n j i ,,2,1, =.)a 如果)(,,)(,)(21t x t x t x n是(5.15)的任意n 个解,那么它们的Wronsky 行列式)](,,)(,)([21t x t x t x W n满足下面的一阶线性微分方程Wt a t a t a W nn )]()()([2211+++=' .(提示:利用行列式的微分公式,求出W '的表达式);)b 解上面的一阶线性微分方程,证明下面的公式:⎰=+++tt nn dss a s a s a e t W t W 02211)]()()([0)()( ,],[,0b a t t∈.证明 )a)()()()()()()()()()()()()()()()()()()(212222111211212222111211t 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 nn n nn n nn n n n n '''++'''='+=∑∑∑===)()()()()()()()()()()()(212222111121111t x t x t x t x t x t x t x t at x t at x t ann n n n nk kn knk k knk k k∑∑∑===+nk kn nknk k nknk k nkn n t x t at x t at x t at x t x t x t x t x t x 112112222111211)()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()(21222211121121222211121111t x t x t x 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 x t x t x t x t a nn n n n n nn nn n n n n++=)()]()([11t W t a t a nn ++= ,所以)(t W 是一阶线性微分方程Wt a t a t a W nn )]()()([2211+++=' 的解.)b 由)a 知,Wt a t a t aW nn )]()()([2211+++=' ,分离变量后两边积分求解得⎰=+++tt nn dss a s a s a cet W 02211)]()()([)( ,t t =时就得到)(0t W c =,所以⎰=+++tt nn dss a s a s a et W t W 02211)]()()([0)()( ,],[,0b a t t ∈.3.设)(t A 为区间],[b a 上的连续n n ⨯实矩阵,)(t Φ为方程x t A x )(='的基解矩阵,而)(t x ϕ=为其一解.试证:)a 对于方程yt Ay T)(-='的任一解)(t ψ必有=)()(t t Tϕψ常数;)b )(t ψ为方程yt Ay T)(-='的基解矩阵的充要条件是存在非奇异的常数矩阵C ,使Ct t T=Φψ)()(.证明)a 由于)(t ϕ是方程x t A x )(='的解,故有)()()(t t A t ϕϕ=',)(t ψ为方程yt A y T )(-='的解,故)()()(t t A t T ψψ-='.所以[][])()()()]([)()()()()()(t t t t t t t t t t TTTTTϕψϕψϕψϕψϕψ'+'='+'=')()()()()]()([t t A t t t t A TT T ϕψϕψ+-=)()()()()()(=+-=t t A t t t A t T T ϕψϕψ,所以=)()(t t Tϕψ常数.)b “⇒” )(t Φ是方程x t A x )(='的基解矩阵,因此)()()(t t A t Φ=Φ',)(t ψ是方程yt Ay T)(-='的基解矩阵,故)()()(t t A t T ψ-=ψ',且0)(det ≠Φt 和0)(det ≠t ψ.所以[][])()()()]([)()()()()()(t t t t t t t t t t TTTTTΦ'ψ+Φψ'=Φ'ψ+Φ'ψ='Φψ)()()()()]()([t t A t t t t A TTTΦψ+Φψ-=)()()()()()(=Φψ+Φψ-=t t A t t t A t T T , 故)()(t t TΦψ是常数矩阵,设Ct t T=Φψ)()(,则)(det )(det )(det )(det )]()(det[det ≠Φ⋅ψ=Φ⋅ψ=Φψ=t t t t t t C T T ,因此存在非奇异常数矩阵C ,使Ct t T=Φψ)()(.“⇐”若存在非奇异常数矩阵C ,使Ct t T=Φψ)()(,则有)(det )(det )(det )(det )]()(det[det 0t t t t t t C T T Φ⋅ψ=Φ⋅ψ=Φψ=≠,所以0)(det ≠ψt ,即)(t ψ是非奇异矩阵或说)(t ψ的各列是线性无关的.又[])()()()()]([)()()(])([)()(0t t A t t t t t t t t t T T T t T Φψ+Φψ'=Φ'ψ+Φ'ψ='Φψ=,并注意到)(det ≠Φt ,有)()()]([t A t t T T ψ-=ψ',即)()()(t t A t T ψ-=ψ'.从而)(t ψ是方程yt Ay T)(-='的基解矩阵.4.设)(t Φ为方程Ax x ='(A 为n n ⨯常数矩阵)的标准基解矩阵(即E =Φ)0(),证明)()()(001t t t t -Φ=ΦΦ-,其中0t 为某一值.证明 由于A 为n n ⨯常数矩阵,故A 在),(∞+-∞有定义、连续,从而它的解也在),(∞+-∞连续可导.由)(t Φ为方程Ax x ='的基解矩阵,故),(∞+-∞∈∀t ,有0)(det ≠Φt ,并且有)()(t A t Φ=Φ',从而对某个0t ,有)(det 0≠-Φt t ,且)()()()(])([00000t t A t t t t t t t t -Φ=-Φ'='-⋅-Φ'='-Φ,即)(0t t -Φ亦为方程Ax x ='的基解矩阵.由推论2*,存在一个非奇异常数矩阵G ,使得在区间),(∞+-∞上,G t t t )()(0Φ=-Φ.又因为Gt t tE )()()0(000Φ=-Φ=Φ=,所以)(01t G -Φ=.因此)()()(001t t t t -Φ=ΦΦ-,其中0t 为某一值.5.设)(,)(t f t A 分别为在区间],[b a 上连续的n n ⨯矩阵和n 维列向量.证明方程组)()(t f x t A x +='存在且最多存在1+n 个线性无关解. 证明 设方程组xt A x )(='的基解矩阵为)](,,)(,)([)(21t t t t n ϕϕϕ =Φ,而)(~t ϕ是方程组)()(t f x t A x +='的一个特解,则其通解为)(~)(t c t x ϕ+Φ=,其中c 是任意的常数列向量.若)(t f 不恒为0,则)(~t ϕ必与)(,,)(,)(21t t t n ϕϕϕ 线性无关,从而)(~t ϕ,)(~)(1t t ϕϕ+,)(~)(2t t ϕϕ+,)(~)(,2t t ϕϕ+ 线性无关,即方程组)()(t f x t A x +='存在1+n 个线性无关解.又假若)(t x 是方程组)()(t f x t A x +='的任意一个解,则一定有确定的常数列向量c ,使得)(~)()(t c t t x ϕ+Φ=,将其加入)(~t ϕ,)(~)(1t t ϕϕ+,)(~)(2t t ϕϕ+,)(~)(,2t t ϕϕ+ 这一组向量就线性相关,故方程组)()(t f x t A x +='的任何2+n 个解必线性相关.从而方程组)()(t f x t A x +='存在且最多存在1+n 个线性无关解.6.试证非齐线性微分方程组的叠加原理:设)(,)(21t x t x 分别是方程组)()(1t f x t A x +=',)()(2t fx t A x +='的解,则)()(21t x t x +是方程组)()()(21t f t f x t A x ++='的解. 证明 因为)(,)(21t x t x 分别是方程组)()(1t f x t A x +=',)()(2t fx t A x +='的解,故)()()()(111t f t x t A t x +=',)()()()(222t f t x t A t x +=',所以有)]()()([)]()()([)()(])()([22112121t f t x t A t f t x t A t x t x t x t x +++='+'='+)()()]()()[(2121t f t f t x t x t A +++=,所以)()(21t x t x +是方程组)()()(21t f t f x t A x ++='的解. 7.考虑方程组)(t f Ax x +=',其中⎪⎪⎭⎫⎝⎛=2012A ,⎪⎪⎭⎫ ⎝⎛=21x x x ,⎪⎪⎭⎫⎝⎛=t t t f cos sin )(. )a 试验证⎪⎪⎭⎫ ⎝⎛=Φt t te te e t 2220)(是Ax x ='的基解矩阵;)b 试求)(t f Ax x +='的满足初始条件⎪⎪⎭⎫⎝⎛-=11)0(ϕ的解)(t ϕ.证明)a 00)(det 4222≠==Φtt t te ete e t ,),(∞+-∞∈∀t 成立.而)(0201220)12(2)(222222t A e te e e e t e t t t tt t tΦ=⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛+=Φ',所以)(t Φ是Ax x ='的基解矩阵.)b ⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛-=Φ--10101)(222241s e e se e es s s s s s,这样,由定理8,方程组满足初始条件⎪⎪⎭⎫ ⎝⎛=00)0(ψ的解就是⎰⎰⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛=ΦΦ=--t s t t ttds s s s e e te e ds s f s t t 0222201cos sin 1010)()()()(ψ⎰⎪⎪⎭⎫⎝⎛-⎪⎪⎭⎫ ⎝⎛=-t s t t tds s s s s e e te e 02222cos cos sin 0⎪⎪⎪⎪⎭⎫ ⎝⎛+-++++--⎪⎪⎭⎫ ⎝⎛=--52)cos 2(sin 51252)cos 2sin 14sin 5cos 10(251022222t t e t t t t t t e e te e t tt tt⎪⎪⎪⎪⎭⎫ ⎝⎛-+--+=)cos 2sin 2(51)cos sin 75(252222t t e t t e te t tt ,对应的齐线性方程组满足初始条件⎪⎪⎭⎫⎝⎛-=11)0(h ϕ的解就是⎪⎪⎭⎫⎝⎛--=⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛=ΦΦ=--t t t t th h e t e E e te e t t 2212221)1(110)0()0()()(ϕϕ,所以,所求方程组)(t f Ax x +='的满足初始条件⎪⎪⎭⎫⎝⎛-=11)0(ϕ的解为⎪⎪⎪⎪⎭⎫⎝⎛-+-+--=+=)cos 2(sin 5153)cos sin 7(252)1527(251)()()(22t t e t t t e t t t t t h ψϕϕ.8.试求)(t f Ax x +=',其中⎪⎪⎭⎫ ⎝⎛=2012A ,⎪⎪⎭⎫ ⎝⎛=21x x x ,⎪⎪⎭⎫⎝⎛=t t t f cos sin )( 满足初始条件⎪⎪⎭⎫⎝⎛-=11)0(ϕ的解)(t ϕ.解 由上题知⎪⎪⎭⎫⎝⎛--=t t h e t e t 22)1()(ϕ,且这里⎰⎰⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛=ΦΦ=--t s s t t ttds e s e e te e ds s f s t t 0222220101010)()()()(ψ⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛=⎪⎪⎭⎫⎝⎛-⎪⎪⎭⎫ ⎝⎛=⎰t t t t t t t t tte e t t t e te e ds s e te e 222222202222121010,所以,所求方程组)(t f Ax x +='的满足初始条件⎪⎪⎭⎫⎝⎛-=11)0(ϕ的解为⎪⎪⎪⎭⎫ ⎝⎛-+-=+=t t h e t e t t t t t 222)1()211()()()(ψϕϕ.9.试求下列方程的通解:)a t x x sec =+'',22ππ<<-t ; )b te x x 28=-'''; )c te x x x =+'-''96.解 )a 易知对应的齐线性方程0=+''x x 的基本解组为t t x cos )(1=,t t x sin )(2=,用公式(5.31)来求方程的一个解.这时1cos sin sin cos )](,)([21=-=tt t t t x t x W ,取0=t,有 ⎰⎰-=-=t t t sdss t s t ds s f s x s x W s x t x s x t x t 0212112sec )sin cos cos (sin )()](,)([)()()()()(0ϕtt t t sds t ds t tt cos ln cos sin tan cos sin 0+=-=⎰⎰所以方程的通解为tt t t t c t c x cos ln cos sin sin cos 21+++=. )b 由于特征方程083=-λ的根是21=λ,i313,2±-=λ,故对应的齐线性方程的基本解组为te t x 21)(=,te t x t 3cos )(2-=,tet x t3sin )(3-=.原方程的一个特解由公式(5.29)有(取0=t),∑⎰==313213210)()](,)(,)([)](,)(,)([)()(k tt k k dss f s x s x s x W s x s x s x W t x t ϕ,其中)](,)(,)([)(321t x t x t x W t W =)3sin 3cos 3(2)3sin 33(cos 24)3sin 3cos 3()3sin 33(cos 23sin 3cos 222t t e t t e e t t e t t e e te te e t t tt t tt t t +----+-=------312=,)](,)(,)([)(3211t x t x t x W t W =)3sin 3cos 3(2)3sin 33(cos 21)3sin 3cos 3()3sin 33(cos 03sin 3cos 0t t e t t e t t e t t e te te t t t t t t +----+-=------te 23-=,)](,)(,)([)(3212t x t x t x W t W =)3cos 33sin 3()3sin 3cos 3(214)3sin 3cos 3(023sin 0222t t e t t e e t t e e te e t t tt tt t -=+--=---,)](,)(,)([)(3213t x t x t x W t W =)3sin 33cos 3(1)3sin 33(cos 240)3sin 33(cos 203cos 222t t e t t e e t t e e te e t t tt tt t +-=--+-=---.所以⎰⎰-+⋅=--ts s tts stdse s s e t e ds e eet 020222312)3cos 33sin 3(3cos 3123)(ϕ⎰+-+-ts s tdse s s e t e 02312)3sin 33cos 3(3sin)3cos 33(sin 324124112122t t e e te t t t ++-=-,故通解tt tte t c t c e ec t x 23221121)3sin 3cos ()(+++=-.)c 特征方程0962=+-λλ,得到特征根32,1=λ,故对应的齐线性方程的基本解组为te t x 31)(=,tte t x 32)(=,tttt tee t ete e t W 63333)31(3)(=+=.取0=t,由(5.31),得特解⎰⎰⋅-=-=t sss t st tt dse e se e e te ds sf s W s x t x s x t x t 06333321120)()()()()()()(ϕtt t ts t e te e ds e s t e 33023412141)(++=-=⎰-,所以得到通解tt e et c ct x 41)()(321++=.10.给定方程)(78t f x x x =+'+'',其中)(t f 在+∞<≤t 0上连续,试利用常数变易公式,证明:)a 若)(t f 在+∞<≤t 0上有界,则上面方程的每一个解在+∞<≤t 0上有界;)b 若当∞→t 时,0)(→t f ,则上面方程的每一个解)(t ϕ,满足0)(→t ϕ(当∞→t 时). 证明 对应的特征方程0782=++λλ有特征根7,1--,故对应的齐线性方程的基本解组te t x -=)(1,tet x 72)(-=,ttt t tee e e e t W 87767)(------=--=.由公式(5.31)得原方程的一个特解(0=t)为⎰⎰-------=-=t s st st tt dss f e e e e e ds s f s W s x t x s x t x t 08772112)(6)()()()()()()(~0ϕ⎰⎰---=t s t t s t dss f e e ds s f e e 0770)(61)(61,所以方程的任一解可写为⎰⎰-----++=t st t s t ttdss f e e ds s f e e ec e c t 0770721)(61)(61)(ϕ.)a 由于)(t f 在+∞<≤t 0上有界,故0>∃M ,),0[∞+∈∀t ,有M t f ≤)(.又由于10≤<-te ,107≤<-te,从而当),0[∞+∈t 时,⎰⎰⋅+⋅++≤--ts t ts t ds e M e ds e M e c c t 0770216161)(ϕ=)1(42)1(67721-+-++--tt t t e e M e e M c c)1(42)1(6721t t e M e M c c ---+-++=M c c 21421++<,即方程的每一个解在+∞<≤t 0上有界.)b 当∞→t 时,0)(→t f ,故由⎰⎰-----++=ts t ts t t t ds s f e e ds s f e e e c e c t 0770721)(61)(61)(ϕ知,若⎰t sdss f e)(有界,则)(0)(610∞→→⎰-t ds s f e e t st ,若⎰t sdss f e)(无界,由于)(s f 在),0[∞+连续,故⎰t s dss f e 0)(为无穷大量,因此0)(lim 616)(lim 6)(lim )(61lim 00====∞→∞→∞→-∞→⎰⎰t f et f e e ds s f e ds s f e e t t t t t tst t s t t ,即总有)(0)(610∞→→⎰-t ds s f e e t st .同理)(0)(61077∞→→⎰-t ds s f e e t st .从而对方程的每一个解)(t ϕ,有)(0)(∞→→t t ϕ.11.给定方程组x t A x )(=',这里)(t A 是区间],[b a 上的连续n n ⨯矩阵.设)(t Φ是它的一个基解矩阵,n 维向量函数),(x t F 在∞<≤≤x b t a ,上连续,],[0b a t∈.试证明初值问题:⎩⎨⎧=+='ηϕ)(,),()(0t x t F x t A x(*)的唯一解)(t ϕ是积分方程组⎰--ΦΦ+ΦΦ=tt dss x s F s t t t t x 0))(,()()()()()(101η (**)的连续解.反之,(**)的连续解也是初值问题(*)的解. 证明)(t ϕ是初值问题(*)的解,故))(,()()()(t t F t t A t ϕϕϕ+=',这说明),(x t F 是t 的向量函数,于是由公式(5.27)得⎰--ΦΦ+ΦΦ=t t ds s s F s t t t t 0))(,()()()()()(101ϕηϕ,即)(t ϕ是积分方程组(**)的连续解.反之,设)(t ϕ是积分方程组(**)的连续解,则有⎰--ΦΦ+ΦΦ=t t ds s s F s t t t t 0))(,()()()()()(101ϕηϕ,两端对t 求导,就有))(,()()())(,()()()()()(11010t t F t t ds s s F s t t t t t t ϕϕηϕ---ΦΦ+ΦΦ'+ΦΦ'='⎰))(,(]))(,()()()[(0101t t F ds s s F s t t tt ϕϕη+Φ+ΦΦ'=⎰-- ))(,(]))(,()()()[()(0101t t F ds s s F s t t t A t t ϕϕη+Φ+ΦΦ=⎰-- ))(,(]))(,()()()()()[(0101t t F ds s s F s t t t t A t t ϕϕη+ΦΦ+ΦΦ=⎰--))(,()()(t t F t t A ϕϕ+=,即)(t ϕ也是初值问题(*)的解.§5.3 常系数线性微分方程组习题5.31.假设A 是n n ⨯矩阵,试证:)a 对任意的常数21,c c 都有A c A c A c A c 2121exp exp )exp(⋅=+;)b 对任意整数k ,都有kAA kexp )(exp =.(当k是负整数时,规定kk A A --=])[(exp )(exp 1.证明 )a 因为))(())((1222121A c A c A c c A c A c ==,所以矩阵Ac 1与A c 2可交换,故Ac A c A c A c 2121exp exp )exp(⋅=+.)b ①先证明N k ∈∀,有kAA kexp )(exp =,这只须对k 施以数学归纳法. 当1=k 时,)1exp(exp )(exp 1A A A ⋅==成立,设当k 时,kAA k exp )(exp =,则当1+k 时,有Ak A kA A A A k k )1exp(exp exp exp )(exp )(exp 1+===+,故对一切自然数k ,kAA kexp )(exp =.②)0exp(0exp )(exp 0A E A ===.③若k 是负整数,则N k ∈-,注意到)exp()(exp 1A A -=-,并由以上证明应用于矩阵A -,就有kAA k A A A k k k exp )](exp[)][exp(])[(exp )(exp 1=--=-==---,由①②③,对一切整数k ,均有kAA kexp )(exp =.2.试证:如果)(t ϕ是Ax x ='满足初始条件ηϕ=)(0t 的解,那么ηϕ)]([exp )(0t t A t -=.证明 由于 ηηϕ⋅⋅-='-='A t t A t t A t )]([exp ])([exp )(0,)(})]({[exp 0t A t t A A ϕη=-=,又ηηηϕ==⋅=E A t )]0[exp()(0,故ηϕ)]([exp )(0t t A t -=是方程组Axx ='满足初始条件ηϕ=)(0t 的解.由解的唯一性,命题得证.3.试计算下列矩阵的特征值及对应的特征向量.)a ⎪⎪⎭⎫ ⎝⎛3421; )b ⎪⎪⎪⎭⎫⎝⎛---244354332;)c ⎪⎪⎪⎭⎫⎝⎛-102111121;)d ⎪⎪⎪⎭⎫ ⎝⎛---6116100010.解 )a 特征方程0543421)det(2=--=----=-λλλλλA E ,特征值11-=λ,52=λ,对应于特征值11-=λ的特征向量⎪⎪⎭⎫ ⎝⎛=21u u u 必须满足方程组0)(1=+-u E A λ,得到0≠∀α,⎪⎪⎭⎫ ⎝⎛-=11αu 是对应于特征值11-=λ的特征向量.类似地可求得对应于特征值52=λ的特征向量为⎪⎪⎭⎫ ⎝⎛=21βv ,其中0≠β的任意常数.)b 特征方程0)2)(1)(2(244354332)det(=++-=---+---=-λλλλλλλA E ,特征值21-=λ,12-=λ,23=λ.对应于特征值21-=λ的特征向量u 必须满足方程组0)(1=+-u E A λ,得到≠∀α,⎪⎪⎪⎭⎫ ⎝⎛=110αu 是对应于特征值21-=λ的特征向量.类似地,可以求出对应于特征值12-=λ以及23=λ的特征向量分别为⎪⎪⎪⎭⎫ ⎝⎛=011βv (0≠β的任意常数)和⎪⎪⎪⎭⎫ ⎝⎛=111γw (0≠γ的任意常数).)c 特征方程0)1)(3(12111121)det(2=+-=---+----=-λλλλλλA E ,特征值12,1-=λ,33=λ.对应于特征值12,1-=λ的特征向量⎪⎪⎪⎭⎫ ⎝⎛=321u u u u 必须满足方程组0)(1=+-u E A λ,得0≠∀α,⎪⎪⎪⎭⎫ ⎝⎛--=212αu 是对应于特征值12,1-=λ的特征向量.类似地,可以求出对应于特征值33=λ的特征向量为⎪⎪⎪⎭⎫ ⎝⎛=212βv (0≠β的任意常数).)d 特征方程0)3)(2)(1(61161001)det(=+++=+--=-λλλλλλλA E ,特征值11-=λ,22-=λ,33-=λ.由0)(1=+-u E A λ,推出0≠∀α,⎪⎪⎪⎭⎫ ⎝⎛-=111αu 是对应于特征值11-=λ的特征向量.同样可求得对应于特征值22-=λ和33-=λ的特征向量分别为⎪⎪⎪⎭⎫⎝⎛-=421βv (0≠β的任意常数)和⎪⎪⎪⎭⎫ ⎝⎛-=931γw (0≠γ的任意常数).4.试求方程组Ax x ='的一个基解矩阵,并计算Atexp ,其中A 为:)a ⎪⎪⎭⎫⎝⎛--2112;)b ⎪⎪⎭⎫⎝⎛3421;)c ⎪⎪⎪⎭⎫⎝⎛---244354332;)d ⎪⎪⎪⎭⎫⎝⎛--115118301.解)a 特征方程032112)det(2=-=--+=-λλλλA E ,得32,1±=λ是特征值.对应的特征向量分别为⎪⎪⎭⎫⎝⎛-=3211αu ,⎪⎪⎭⎫ ⎝⎛+=3212βu ,0,0≠≠βα为任意常数.所以方程组Axx ='的一个基解矩阵为⎪⎪⎭⎫ ⎝⎛+-=Φ--t ttt e e ee t 3333)32()32()(.133331323211)32()32()0()(exp ----⎪⎪⎭⎫ ⎝⎛+-⎪⎪⎭⎫ ⎝⎛+-=ΦΦ=t ttt e e ee t At⎪⎪⎭⎫⎝⎛--+----+=----t ttttt tt e eee e ee e 33333333)32()32()32()32(63.)b 由第3题)a 立即得到方程组Ax x ='的一个基解矩阵为⎪⎪⎭⎫⎝⎛-=Φ--t tt te e e e t 552)(. 155121112)0()(exp ----⎪⎪⎭⎫⎝⎛-⎪⎪⎭⎫ ⎝⎛-=ΦΦ=t tt t e e e e t At⎪⎪⎭⎫ ⎝⎛+--+=----t t t t t t tt e e e e e e e e 55552)(2231.)c 由第3题)b 立即得到方程组Ax x ='的一个基解矩阵为⎪⎪⎪⎭⎫ ⎝⎛=Φ----t t t t tt t e e e e ee e t 222220)(.12222211011111100)0()(exp ------⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛=ΦΦ=t t tt tt t e e e e ee e t At⎪⎪⎪⎭⎫⎝⎛----+---=--------t tt t t tt tt t t tt t t t te e e e e e e e e e e e e e e e e 2222222222222. )d 特征方程)34)(3(11511831)det(2=--+=+------=-λλλλλλλA E ,特征值为31-=λ,723,2±=λ.对应的特征向量分别为⎪⎪⎪⎭⎫ ⎝⎛-=4731αu ,⎪⎪⎪⎭⎫ ⎝⎛++-=7174532βu ,⎪⎪⎪⎭⎫⎝⎛+-+-=7174533γu ,γβα,,均为不等于零的任意常数.故方程组Ax x ='的一个基解矩阵为⎪⎪⎪⎪⎭⎫⎝⎛-++---=Φ-+--+--+-t tt tt tttt e e e ee e e e e t )72()72(3)72()72(3)72()72(3)17()17(4)574()574(7333)(.由)0()(exp 1-ΦΦ=t At 立即可得[])()()(exp 321t t t At ψψψ=,其中列向量函数⎪⎪⎪⎪⎭⎫⎝⎛-+++--+++--+++=-+--+--+-t t t t t t t t t e e e ee e e e e t )72()72(3)72()72(3)72()72(31)7514(2)7514(256)71349()49713(98)737(3)737(342841)(ψ, ⎪⎪⎪⎪⎭⎫⎝⎛++-+-++-+-+-+-+=-+--+--+-t t t tt t t t t e e e e e e e e e t )72()72(3)72()72(3)72()72(31)714(2)714(256)753175()753175(98)757(3)757(3422521)(ψ,⎪⎪⎪⎪⎭⎫⎝⎛-++++--+++-+-=-+--+--+-t t t tt t t t t e e e e e e e e e t )72()72(3)72()72(3)72()72(31)7137()7137(112)98761()98761(196)714(3)714(3841261)(ψ.(该题计算量太大,作为该法的习题不是太好!)5.试求方程组Ax x ='的一个基解矩阵,并求满足初始条件ηϕ=)0(的解)(t ϕ:)a ⎪⎪⎭⎫⎝⎛=3421A ,⎪⎪⎭⎫⎝⎛=33η;)b ⎪⎪⎪⎭⎫⎝⎛--=115118301A ,⎪⎪⎪⎭⎫⎝⎛--=720η;)c ⎪⎪⎪⎭⎫ ⎝⎛-=102111121A ,⎪⎪⎪⎭⎫ ⎝⎛=001η.解 )a 由上题)b 知⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-=--1112231exp 55t tt te e e e At ,所以所求解为⎪⎪⎭⎫⎝⎛+-+==--t t t t e e e e At t 5542)(exp )(ηϕ.)b 由上题)d 知)0()(exp 1-ΦΦ=t At ,其中⎪⎪⎪⎪⎭⎫⎝⎛-++---=Φ-+--+--+-t tt tt tttte e e ee e e e e t )72()72(3)72()72(3)72()72(3)17()17(4)574()574(7333)(.所以所求解为⎪⎪⎪⎭⎫⎝⎛--⎪⎪⎪⎭⎫ ⎝⎛++---+--⋅Φ==720)714(2775)773(3)714(2775)773(33214422521)()(exp )(t At t ηϕ ⎪⎪⎪⎪⎭⎫⎝⎛-++--+--+-++-+=-+--+--+-t t t tt t t t t e e e e e e ee e )72()72(3)72()72(3)72()72(3)7317(3)78977(728)7160289(3)7374511(1274)7435(9)9172(35461261.)c 由第3题)c 知,矩阵A 的特征值为12,1-=λ,33=λ.对应于特征值33=λ的特征向量⎪⎪⎪⎭⎫ ⎝⎛=212αv (0≠α的任意常数).又由648324648)(32121=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛=-u u u u A E λ,得到⎪⎪⎪⎭⎫⎝⎛+-=)24(3331γβγβu (γβ,是任意常数),由⎪⎪⎪⎭⎫⎝⎛+-+⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=)24(3331212001γβγβαη解出41,21,41-===γβα.依公式(5.52),得满足初始条件ηϕ=)0(的解为⎪⎪⎪⎭⎫⎝⎛--⎪⎪⎪⎭⎫⎝⎛+++⎪⎪⎪⎭⎫ ⎝⎛=+++=--212120212124121241)]([)(33t tt t tt t e e u E A t E e Ev e t t t t t ϕ⎪⎪⎪⎭⎫⎝⎛--+=---)(2)(241333t t tt t t e e e e e e6.试求方程组)(t f Ax x +='的解)(t ϕ:)a ⎪⎪⎭⎫⎝⎛-=11)0(ϕ,⎪⎪⎭⎫ ⎝⎛=3421A ,⎪⎪⎭⎫⎝⎛=1)(t e t f ;)b ⎪⎪⎪⎭⎫ ⎝⎛=000)0(ϕ,⎪⎪⎪⎭⎫ ⎝⎛---=6116100010A ,⎪⎪⎪⎭⎫⎝⎛=-t e t f 00)(;)c ⎪⎪⎭⎫⎝⎛=21)0(ηηϕ,⎪⎪⎭⎫⎝⎛--=1234A ,⎪⎪⎭⎫⎝⎛-=t t t f cos 2sin )(.解 )a 由第4题)b 知,⎪⎪⎭⎫⎝⎛-⎪⎪⎭⎫ ⎝⎛-=--1112231exp 55t tt te e e e At ,由公式(5.61)得⎰-+=t ds s f A s t At t 0)(])exp[()(exp )(ηϕ⎰⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-=--------t s s t s t s t s t t t t tds e e e e e e e e e 0)(5)()(5)(5511112231111112231⎪⎪⎪⎪⎭⎫ ⎝⎛+++-++-=--53109235420934355t t t t tt e e e e e e .)b 由第3题)d 知A 的特征值11-=λ,22-=λ,33-=λ,对应的特征向量分别为⎪⎪⎪⎭⎫⎝⎛-=111αu ,⎪⎪⎪⎭⎫ ⎝⎛-=421βv ,⎪⎪⎪⎭⎫ ⎝⎛-=931γw ,其中γβα,,均是不为零的任意常数.Ax x ='的一个基解矩阵为⎪⎪⎪⎭⎫⎝⎛---==Φ---------t tt tt tt t ttt te e e e e ee e e w e v e u et 3232329432][)(321λλλ.⎪⎪⎪⎭⎫⎝⎛---=⎪⎪⎪⎭⎫ ⎝⎛---=Φ--13228615621941321111)0(11,而)0()(exp 1-ΦΦ=t At .由公式(5.61)得⎰-+=t ds s f A s t At t 0)(])exp[()(exp )(ηϕ⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎭⎫ ⎝⎛---⎪⎪⎪⎭⎫ ⎝⎛---=---------000132286156943221323232t tt tt t t t te e ee e e e e e⎰⎪⎪⎪⎭⎫⎝⎛⎪⎪⎪⎭⎫ ⎝⎛---⎪⎪⎪⎭⎫ ⎝⎛---+-------------------t s s t s t s t s t s t s t s t s t s t dse e e e e e e e e e 0)(3)(2)()(3)(2)()(3)(2)(00132286156943221⎪⎪⎪⎭⎫⎝⎛-+-+---+-=⎪⎪⎪⎭⎫ ⎝⎛+--+-+-=------------------⎰t t t tt t t t t t t s t s t t s t s t t s t s t e e e t e e e t e e e t ds e e e e e e e e e 3232320322322322916)72(38)25(4)32(419834221.)c A的特征方程0)2)(1(1234)det(=--=+--=-λλλλλA E ,求解得特征值11=λ,22=λ,对应的特征向量分别是⎪⎪⎭⎫ ⎝⎛=11αu ,⎪⎪⎭⎫ ⎝⎛=23βv ,其中βα,是不为零的任意常数.所以方程组Axx ='的一个基解矩阵为⎪⎪⎭⎫⎝⎛==Φt tt t tte e e e v eu e t 2223][)(21λλ,从而,⎪⎪⎭⎫⎝⎛--Φ=ΦΦ=-1132)()0()(exp 1t t At .由公式(5.61)得⎰-+=t ds s f A s t At t 0)(])exp[()(exp )(ηϕ⎰⎪⎪⎭⎫⎝⎛-⎪⎪⎭⎫ ⎝⎛--⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛=-----t s t s t s t s t t t t tds s s e ee e e e e e 0)(2)()(2)(2122cos 2sin 113223113223ηη⎪⎪⎭⎫⎝⎛+-+-+-+-+⎪⎪⎭⎫ ⎝⎛-+--+-=t t e e t t e e e e e e t t t t t t t t cos 2sin 224cos sin 234)(2)23()(3)23(222211222112ηηηηηηηη⎪⎪⎭⎫⎝⎛+--+--+--+--=t t e e t t e e t t t t cos 2sin 2)(2)423(cos sin 2)(3)423(2211222112ηηηηηηηη.7.假设m 不是矩阵A 的特征值,试证非齐线性方程组mtce Ax x +='有一解形如mte t ρϕ=)(,其中ρ,c 是常数向量.证明 设方程组有形如mte t ρϕ=)(的解,代入方程得m tm t m t ce e A e m +=ρρ,由此得cA m +=ρρ,即cA mE =-ρ)(.因为m 不是矩阵A 的特征值,故0)det(≠-A mE ,即矩阵A mE -可逆,得到c A mE 1)(--=ρ唯一确定.所以方程组有一解m tm t e ce A mE t ρϕ=-=-1)()(8.给定方程组⎩⎨⎧=+'+-'=-'++'-''.02,023221122111x x x x x x x x x)a 试证上面方程组等价于方程组Au u =',其中⎪⎪⎪⎭⎫⎝⎛'=⎪⎪⎪⎭⎫ ⎝⎛=211321x x x u u u u ,⎪⎪⎪⎭⎫ ⎝⎛---=112244010A ;)b 试求)a 中的方程组的基解矩阵;)c 试求原方程组满足初始条件0)0(1=x ,1)0(1='x ,)0(2=x 的解.解 )a 设11x u=,12x u'=,23x u=,则原方程组化为⎪⎩⎪⎨⎧--='=''-+-=''='='=',2,23,32123331212211u u u x u u u u u x u u x u或⎪⎩⎪⎨⎧--='++-='='32133212212,244,uu u u u u u u u u ,即u u ⎪⎪⎪⎭⎫ ⎝⎛---='112244010或Au u ='.反之,设11u x =,21u x =',32u x=,则方程组Au u ='化为⎩⎨⎧-'-='+'+-=''.211221112,244x x x x x x x x即⎩⎨⎧=+'+-'=-'++'-''.02,023221122111x x x x x x x x x)b 由0)2)(1(11224401)det(=--=+----=-λλλλλλλA E ,得矩阵A的特征值01=λ,12=λ,23=λ.对应的特征向量分别为⎪⎪⎪⎭⎫⎝⎛=201αu ,⎪⎪⎪⎭⎫ ⎝⎛=122βv ,⎪⎪⎪⎭⎫ ⎝⎛=021γw ,其中γβα,,均为不等于零的任意常数.由此得Au u ='的一个基解矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==Φ0222021][)(22321t t tt t t t t e e e e e w e v e u e t λλλ.)c 求与之等价的方程组Au u =',满足初始条件η=⎪⎪⎪⎭⎫⎝⎛=010)0(u 的解ηη)0()()(exp )(1-ΦΦ==t At t u⎪⎪⎪⎭⎫ ⎝⎛-+-+-=⎪⎪⎪⎭⎫⎝⎛⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=-t tt t t t t t t t e e e e e e e e e e 226434121010012220121022202122122,所以,原方程组满足初始条件0)0(1=x ,1)0(1='x ,0)0(2=x 的解为⎪⎪⎭⎫⎝⎛-+-=t t t e e e t 2234121)(2ϕ.9.试用Laplace 变换法解第5题和第6题. 解 5.)a 方程组两边取Laplace 变换,有)()(s AX s sX =-η,即η=-)()(s X A sE ,由具体数值代入得方程组⎪⎪⎭⎫⎝⎛=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛----33)()(342121s X s X s s ,根据Gramer 法则得 5211)(1-++=s s s X ,5411)(2-++-=s s s X,所以tte et -+=512)(ϕ,tte et --=524)(ϕ,故初值问题5.)a 的解为⎪⎪⎭⎫⎝⎛-+=⎪⎪⎭⎫ ⎝⎛=--t t t t e e e e t t t 552142)()()(ϕϕϕ.5.)b 对方程组两边施行Laplace 变换,并化简有η=-)()(s X A sE ,用具体数值代入得方程组⎪⎪⎪⎭⎫ ⎝⎛--=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛+------720)()()(115118301321s X s X s X s s s ,根据Gramer 法则得)72(427291)72(4272913313)34)(3(1521)(21--+-+---+=--++-=s s s s s s s s X ,)72(1267376511)72(12673765113991)34)(3(14372)(222--+++--++-=--+-+-=s s s s s s s s s X ,)72(12678977)72(126789773952)34)(3(5127)(223----+-+-+-=--+-+-=s s s s s s s s s X ,所以ttt e e e t )72()72(31427291427291313)(-+-+---=ϕ,ttt ee e t )72()72(3212673765111267376511991)(-+-++-+-=ϕ,ttt ee e t )72()72(331267897712678977952)(-+---+--=ϕ,故初值问题5.)b 的解为⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛--+--++-+-+---=⎪⎪⎪⎭⎫ ⎝⎛=-+--+--+-t t t t t t t t t e e e e e e ee e t t t t )72()72(3)72()72(3)72()72(3321126789771267897795212673765111267376511991427291427291313)()()()(ϕϕϕϕ.5.)c 对方程组两边施行Laplace 变换,并化简有η=-)()(s X A sE ,用具体数值代入得方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫ ⎝⎛---+----001)()()(102111121321s X s X s X s s s ,根据Gramer 法则得31211121)(1-++=s s s X ,31411141)(2--+=s s s X,31211121)(1-++-=s s s X , 所以)(21)(31t te e t -+=ϕ,)(41)(32t te e t ---=ϕ,)(21)(33t te e t --=ϕ,故初值问题5.)a 的解为。
常微分方程第五章测试题及参考答案
常微分方程第五章测试题班级__________姓名__________学号________得分__________一、 填空(30分)1、 在用皮卡逐步逼近法求方程组η=+=')(),()(0t x x f x t A x 的近似解时,若取ηϕ=)(0t ,则=)(t k ϕ( )。
2、 如果)(t A 是n n ⨯矩阵,)(t f 是n 维列向量,则它们在b t a ≤≤上满足( )时,方程组)()(t f x t A x +='满足初始 条件η=)(0t x 的解在b t a ≤≤上存在唯一。
3、 若)(),(),(21t f t a t a 是[b a ,]上的连续函数,)(),(21t x t x 是方程0)()(21=+'+''x t a x t a x 的两个线性无关解,则的通解为( )。
4、 若)(t Φ和)(t ψ都是x t A x )(='的基解矩阵,则)(t Φ与)(t ψ具有关系( )。
5、 若A 是n n ⨯常数矩阵,则矩阵指数exPA=( )。
6、若A 矩阵具有n 个线性无关的特征向量n v v v ,,21,她们对应的特征值分别为n λλλ ,,21,那么矩阵)(t Φ=( )是常系数线性方程组Ax x ='的一个基解矩阵。
7、 若)(t Φ是x t A x )(=' 的基解矩阵,则)()(t f x t A x +='满足的解=)(t ϕ()。
8、 若)(t Φ是x t A x )(=' 的基解矩阵,则向量函数=)(t ϕ( )是)()(t f x t A x +='的满足初始条件 0)(0=t ϕ的解;向量函数=)(t ϕ( )是)()(t f x t A x +='的满足初始条件ηϕ=)(0t 的解。
9、 方程组x t A x )(='的n 个解)(,),(),(21t x t x t x n 线性无关的充要条件是( )。
常微分方程第四、第五章部分习题参考答案
常微分方程习题4.2 2、解下列方程 (1)045)4(=+''-x x x解:特征方程1122045432124-==-===+-λλλλλλ,,,有根故通解为x=t t t te c e c e c e c --+++432221(2)03332=-'+''-'''x a x a x a x解:特征方程0333223=-+-a a a λλλ有三重根a =λ故通解为x=at at at e t c te c e c 2321++ (3)04)5(=''-x x解:特征方程0435=-λλ有三重根0=λ,=4λ2,=5λ-2故通解为54232221c t c t c e c e c x t t ++++=-(4)0=+'+''x x x解:特征方程012=++λλ有复数根=1λ,231i +-=2λ,231i-- 故通解为t e c t ec xt t 23sin 23cos 212211--+=(5) 12+=-''t s a s解:特征方程022=-a λ有根=1λa,=2λ-a当0≠a 时,齐线性方程的通解为s=atat e c e c -+21Bt A s +=~代入原方程解得21aB A -== 故通解为s=atat e c e c -+21-)1(12-t a当a=0时,)(~212γγ+=t t s 代入原方程解得21,6121==γγ故通解为s=t c c 21+-)3(612+t t (6) 32254+=-'+''-'''t x x x x解:特征方程025423=-+-λλλ有根=1λ2,两重根=λ 1 齐线性方程的通解为x=t t t te c e c e c 3221++又因为=λ0不是特征根,故可以取特解形如Bt A x +=~代入原方程解得A=-4,B=-1 故通解为x=t t t te c e c e c 3221++-4-t (7) 322)4(-=+''-t x x x解:特征方程121201224-===+-λλλλ重根,重根有 故齐线性方程的通解为x=t t t t te c e c te c e c --+++4321 取特解形如c Bt At x ++=2~代入原方程解得A=1,B=0,C=1 故通解为x=t t t t te c e c te c e c --+++4321+12+t (8)t x x cos =-'''解:特征方程013=-λ有复数根=1λ,231i +-=2λ,231i--13=λ 故齐线性方程的通解为t t t e c t e c t ec x 321221123sin 23cos ++=--取特解形如t B t A x sin cos ~+=代入原方程解得A=21,21-=B 故通解为t t t e c t e c t ec x 321221123sin 23cos ++=--)sin (cos 21t t +-(9) t x x x 2sin 82=-'+''解:特征方程022=-+λλ有根=1λ-2,=2λ 1 故齐线性方程的通解为x=tte c e c 221-+因为+-2i 不是特征根取特解形如t B t A x 2sin 2cos ~+=代入原方程解得A=56,52-=-B 故通解为x=tte c e c 221-+t t 2sin 562cos 52--(10)t e x x =-'''解:特征方程013=-λ有复数根=1λ,231i +-=2λ,231i--13=λ 故齐线性方程的通解为t t t e c t e c t ec x 321221123sin 23cos ++=-- =λ1是特征方程的根,故t Ate x =~代入原方程解得A=31 故通解为t t t e c t e c t ec x 321221123sin 23cos ++=--+t te 31(11)t e s a s a s =+'+''22解:特征方程0222=++a a λλ有2重根=λ-a 当a=-1时,齐线性方程的通解为s=t t te c e c 21+,=λ1是特征方程的2重根,故t e At x 2~=代入原方程解得A=21通解为s=22121t te c e c t t ++, 当a ≠-1时,齐线性方程的通解为s=at at te c e c --+21,=λ1不是特征方程的根,故t Ae x =~代入原方程解得A=2)1(1+a故通解为s=at at te c e c --+21+te a 2)1(1+ (12)t e x x x 256=+'+''解:特征方程0562=++λλ有根=1λ-1,=2λ-5 故齐线性方程的通解为x=tte c ec 521--+=λ2不是特征方程的根,故t Ae x 2~=代入原方程解得A=211故通解为x=t te c ec 521--++te 2211 (13)t e x x x t cos 32-=+'-''解:特征方程0322=+-λλ有根=1λ-1+2i,=2λ-1-2i 故齐线性方程的通解为t e c t e c x t t 2sin 2cos 21+=i ±-1 不是特征方程的根, 取特解行如t e t B t A x -+=)sin cos (~代入原方程解得A=414,415-=B 故通解为t e c t e c x t t 2sin 2cos 21+=+t e t t --)sin 414cos 415( (14) t t x x 2cos sin -=+''解:特征方程012=+λ有根=1λi,=2λ- i 故齐线性方程的通解为t c t c x sin cos 21+= 对于t x x sin =+'',=1λi,是方程的解, 设)sin cos (~t B t A t x +=代入原方程解得A=21-B=0 故t t x cos 21~-=对于t x x 2cos -=+'' ,设t B t A x 2sin 2cos ~+=代入原方程解得A=31 B=0 故t x 2cos 31~= 故通解为t c t c x sin cos 21+=t t cos 21-t 2cos 31+ 15)1442++=+'-''ttee x x x解:0442=+-λλ,22,1=λ,齐次方程的通解为)()(212t C C e t x t +=。
常微分方程知到章节答案智慧树2023年东北师范大学
常微分方程知到章节测试答案智慧树2023年最新东北师范大学绪论单元测试1.常微分方程的发展按研究内容可分为几个历史阶段?( )参考答案:定性稳定性理论阶段。
;解析理论阶段;经典阶段;适定性理论阶段2.本课程的主要教学内容有哪些?()参考答案:定性和稳定性理论简介等。
;初等积分法;基本定理;一阶线性微分方程组,n阶线性微分方程3.常微分方程的研究方法主要有哪些?()参考答案:各项均正确第一章测试1.下面方程中是线性方程的有()参考答案:2.下面方程中是齐次方程的是()参考答案:3.方程是常数解()参考答案:4.不是所有的方程都可以用初等积分法求解。
( )参考答案:对5.通解不一定包含微分方程的所有解。
( )参考答案:对第二章测试1.存在且连续是保证方程初值解唯一的必要条件。
( )参考答案:错2.线素场中的线素不能等于0。
( )参考答案:错3.奇解也是方程的解。
( )参考答案:对4.方程满足解的存在唯一性定理条件的区域是()参考答案:除去轴的平面5.方程任意解的存在区间是()参考答案:第三章测试1.函数在区间的朗斯基行列式恒为零是它上线性相关的()参考答案:必要条件2.一阶线性非齐次微分方程组的任两个非零解之差()参考答案:是其对应齐次微分方程组的解3.若的解,为其对应的齐次线性微分方程组的解,则()的解参考答案:是4.线性齐次微分方程组的解组为基本解组的充要条件是它们的朗斯基行列式()参考答案:错5.齐次线性微分方程的基本解组不是唯一的。
()参考答案:对第四章测试1.阶线性齐次微分方程的所有解构成一个()维线性空间.参考答案:2.微分方程的通解中应含的独立常数的个数为().参考答案:33.微分方程的特解具有形式().参考答案:4.若和是二阶线性齐次方程的基本解组,则它们没有共同零点。
()参考答案:对5.只要给出阶线性微分方程的个特解,就能写出其通解.()参考答案:错6.下列方程是二阶线性微分方程的是()。
[理学]5常微分方程_OK
dy ay f ( x)的一切解y( x),均有 lim y( x) b .
dx
x
a
例5-16 若曲线过点N(1,1), 曲线上任意一点P(x,y)处
的切线与 Oy 轴交于Q, 经PQ为直径做的圆过A(1,0) ,
求此曲线方程.
20
(三)特殊的高阶微分方程
例5-17
求方程
y
ln x x2
则 μ(x,y) 应该满足一阶偏微分方程
Q P (P Q ).
x y y x
5.准齐次方程
dy f ( ax by c ).
dx
a1x b1 y c1
当c = c1 = 0时,就是齐次方程;
当c , c1中至少有一个不为零时,总可以做变换,使
之转化为齐次方程;
当ab1 – a1b 0时, 令x = + , y = + ;
特征方程的根
单实根 r k重实根 r 一对虚根 r1,2= i
一对 k 重虚根
r1,2= i
方程通解中对应的项
Cerx ,
erx (C1 C2 x Ck xk1),
ex (C1 cos x C2 sin x),
e x[(C1 C2 x (D1 D2 x
Ck xk1)cos x
8.若 y1 (x)是微分方程
y p( x) y q( x) y 0
的一个解,则该微分方程的通解为
y ( x)[C1 C2
e P( x )dx
2 ( x) dx].
16
二、归类解析
(一)可分离变量的微分方程
例5-1 求微分方程
dy 3xy xy2 的通解.
dx
例5-2 求满足方程 ( x2 y2 1)dx 2x2dy 0
第5章_常微分方程
将 y 视为自变量,可以变成关于 x 的线性方程: dx 1 1 − x= y P( y ) = − , Q( y ) = y dy y y
∴x = e
1 − − dy y
∫
[ ∫ ye
−
∫
1 dy y
dy + C ]
= y( y + C )
由 y | x =3 = 1 得: C = 2 故所求特解为: x = y ( y + 2)
解方程(2x-5y+3)dx-(2x+4y-6)dy=0. 例 解方程
a b 2 解: = a1 b1 2 -5
2 x - 5 y + 3 = 0, ≠ 0 令 4 2 x + 4 y - 6 = 0,
解得x 解得 0=1, y0=1
dy 2 X − 5Y 2 − 5 Y x = X + 1, X 则 = = 令 dx 2 X + 4Y 2 + 4 Y y = Y + 1, X Y dY du 令u = , 有 =u+ X X dX dX du 2 − 5u 4u + 2 1 方程变为u + X = ,即 2 du = − dX dX 2 + 4u 4u + 7u − 2 X 4u + 2 2 1 4 1 1 du = ∫ ( ⋅ + ⋅ )du = ln | (u + 2) 2 (4u − 1) | +c ' ∫ 4u 2 + 7u − 2 3 u + 2 3 4u − 1 3
二.齐次方程 齐次方程 如果方程(1)可化成: 令u=
y 解法: 化成可分离变量方程. x dy du y = xu =u+x dx dx du 1 du = dx ∴u + x = ϕ (u ) ϕ (u ) − u x u) dx
江苏省专转本高等数学第五章常微分方程核心知识点例题讲解(含答案)
第五章 常微分方程(简记ODE )本章主要知识点● 可分离变量的ODE● 一阶线性非齐次常微分方程及推广● 二阶常系数线性齐次与非齐次常微分方程● 一些特殊类方程一、可分离变量的ODE1.基本型的解法 基本型:()()dy G x H y dx= 基本解法: ()()dy G x dx H y = ()()dy G x dx H y =⎰⎰例5.1.1)0(,==-y e dx dy y x 解:dx e dy e x y =⎰⎰=dx e dy e x y通解为:c e e x y += 将1,0==y x 得:1-=e c 得 1-+=e e e x y例5.2.(1)ln y y y xdx '+= 解:(1)ln y dy xdx y+= 1(1)ln dy xdx y +=⎰⎰,得:ln ||ln y y x x x C +=-+例5.3.dx y x dy y x )1()1(122+=+-解:dx x x y dy y 2211)1(-=++,2(1)1y dy y +=+⎰ 得:()21arctan ln 12y y C ++= 例5.4.已知()f x 满足0()(1)()1x f t dt x f x +-=⎰,求()f x 。
解:由0()(1)()1xf t dt x f x +-=⎰知(0)1f =-。
方程两边对x 求导得()()(1)()0f x f x x f x '++-=,分离变量求得2()(1)c f x x =-, 将(0)1f =-代入得1c =-,21()(1)f x x =--。
2.可转化的可分离变量的齐次方程 ()x y f y'= 方法:令()y p y p x x y p xp x''=⇒=⇒=+ xdx p p f dp p f dx dp x p =-⇒=+⇒)()(。
例5.5.y x y x dx dy +-= 解:xyx ydx dy +-=11 令p p dx dp x p xp p y px y x y p +-=+⇒+=⇒=⇒=11'', pp p p p p dx dp x +--=-+-=⇒121112 xdx p p dp p =--+⇒221)1( x dx p dp p =+-+⇒⎰2)1(2)1( C x p p +=---⇒ln 21ln 212,将xy p =代入即可。
常微分方程习题答案(第五章定性与稳定性理论简介)
常微分方程习题答案第五章定性与稳定性理论简介教材习题同步解答习题5.21. 对于方程组41114221,,xx x x x x ⎧=-⎨=⎩ 试说明 441212(,)V x x x x =+是正定的,而dVdt是常负的。
证:易知(0,0)0V =,当22120x x +≠时,12(,)0V x x > 正定。
34344444121122211212124()4()440dV V V x x x x x x x x x x x x dt x x ∂∂=+=-+-=-+=∂∂ ,故dV dt是常负。
(0,0)0V =。
2. 讨论方程组312132124,3,xx x x x x ⎧=--⎨=-⎩ 零解的稳定性。
证:取 221212(,)V x x x x =+, 易知(0,0)0V =,当22120x x +≠时, 12(,)0V x x >即正定。
334411221212121212222(4)2(3)22()0dV x x x x x x x x x x x x x x dt=+=--+-=---< ,故方程的零解是渐进稳定的。
3. 讨论自治系统2111222212,,x Ax x x x Ax x x ⎧=-⎨=-⎩ 零解的稳定性。
证:证:取 221212(,)V x x x x =+, 易知(0,0)0V =,当22120x x +≠时,12(,)0V x x >即正定。
222211221112221212222()2()2()dV x x x x x Ax x x x Ax x x A x x dt=+=-+-=+ ,故方程的0A >,则零解是不稳定的;若0A <,则零解是渐进稳定的。
习题5.3通过求解,确定下列各方程的奇点类型,画出相图,并确定奇点的稳定性:(1)2,3;dx x dt dy y dt ⎧=-⎪⎪⎨⎪=-⎪⎩(2)3,3;dx x dt dy x y dt⎧=⎪⎪⎨⎪=+⎪⎩(3),;dx y dt dy x dt ⎧=⎪⎪⎨⎪=-⎪⎩(4)23,3;dxx y dtdy x y dt ⎧=+⎪⎪⎨⎪=+⎪⎩解:(1)方程的奇点为(0,0)O ,方程所对应的系数矩阵为2003A -⎡⎤=⎢⎥-⎣⎦,系数矩阵所对应的特征方程为20003λλ--=-- 或2560λλ++= ,特征根为 1220,30,λλ=-<=-<奇点(0,0)O 为稳定结点。
常微分方程第5章答案
又因为=-A (t) (t),所以=- (t) A(t)
[ (t) (t)] =- (t) (t)A(t)+ (t) A(t) (t)=0,
所以对于方程y =-A (t)y的任一解y= (t)必有(t) (t)=常数
b)“”假设为方程y =-A (t)y的基解矩阵,则
又v(0)= =
v (t)= = = v(t)
因此u(t),v(t)分别是给定初值问题的解.
b) w(0)= u(0)+ u(0)= + =
w (t)= u (t)+ v (t)
= +
=
=
= w(t)
因此w(t)是给定方程初值问题的解.
2.将下面的初值问题化为与之等价的一阶方程组的初值问题:
a) x +2x +7tx=e ,x(1)=7, x (1)=-2
习题5.1
1.给定方程组
x = x x=(*)
a)试验证u(t)= ,v(t)=分别是方程组(*)的满足初始条件u(0)= , v(0)=的解.
b)试验证w(t)=c u(t)+c v(t)是方程组(*)的满足初始条件w(0)=的解,其中是任意常数.
解:a) u(0)= =
u (t)= = u(t)
6、试证非齐线性微分方程组的叠加原理:
的解,则是方程组
的解。
证明:(1)(2)
分别将代入(1)和(2)
则
则
令
即证
7.考虑方程组,其中
a)试验证是的基解矩阵;
b)试求的满足初始条件的解。
证明:a)首先验证它是基解矩阵
国家开放大学电大本科《常微分方程》网络课形考任务1-6试题及答案
国家开放大学电大本科《常微分方程》网络课形考任务1-6试题及答案国家开放大学电大本科《常微分方程》网络课形考任务1-6试题及答案100%通过考试说明:2022年秋期电大把该网络课纳入到“国开平台”进行考核,该课程共有6个形考任务,针对该门课程,本人汇总了该科所有的题,形成一个完整的标准题库,并且以后会不断更新,对考生的复习、作业和考试起着非常重要的作用,会给您节省大量的时间。
做考题时,利用本文档中的查找工具,把考题中的关键字输到查找工具的查找内容框内,就可迅速查找到该题答案。
本文库还有其他网核及教学考一体化答案,敬请查看。
课程总成绩=形成性考核某50%+终结性考试某50%形考任务1题目1本课程的教学内容共有五章,其中第三章的名称是().选择一项: A.一阶线性微分方程组B.定性和稳定性理论简介C.初等积分法D.基本定理题目2本课程安排了6次形成性考核任务,第2次形成性考核作业的名称是().选择一项: A.第一章至第四章的单项选择题B.第二章基本定理的形成性考核书面作业C.初等积分法中的方程可积类型的判断D.第一章初等积分法的形成性考核书面作业题目3网络课程主页的左侧第3个栏目名称是:().选择一项: A.课程公告B.自主学习C.课程信息D.系统学习题目4网络课程的“系统学习”栏目中第一章初等积分法的第4个知识点的名称是().选择一项: A.一阶隐式微分方程B.分离变量法C.全微分方程与积分因子D.常数变易法题目5网络课程的“视频课堂”栏目中老师讲课的电视课共有()讲.选择一项: A.18B.20C.19D.17题目6网络课程主页的左侧“考试复习”版块中第二个栏目名称是:().选择一项: A.考核说明B.复习指导C.模拟测试D.各章练习汇总题目7请您按照课程的学习目标、学习要求和学习方法设计自己的学习计划,并在下列文本框中提交,字数要求在100—1000字.答:常微分方程是研究自然现象,物理工程和工程技术的强有力工具,熟练掌握常微分方程的一些基本解法是学习常微分方程的主要任务,凡包含自变量,未知函数和未知函数的导数的方程叫做微分方程。
常微分方程课后答案
d2 y dx2
−
1 1−x2
y
=
1
+
x,
y(0) = 1;
(3) y = ex +
x 0
y(t)
dt;
(4)
dy dx
=
; x4 +y 3
xy2
(5) 2xydy − (2y2 − x)dx = 0;
(6) (y ln x − 2)ydx = xdy;
(7)
3xy2
dy dx
+
y3
+
x3
=
0;
(8)
y
=
0(其中c是任意常数);
(5)
y = ecx,
(
dy dx
)2
−
y
d2 y dx2
= 0(其中c是任意常数);
2
−
(x−C1 4
)2
,
−∞ < x < C1;
(6)
y =
0, C1 < x < C2,
(x−C1 4
)2
,
C2 < x < +∞,
dy dx
=
|y|.
答:将解代入验证就可得知是否为微分方程的解:
dy dx
=
y x+y
3
.
3
解: (1)方程两边同时乘以因子e−x,由此得到方程的通解为
y
=
C ex
−
sin
x
+ 2
cos x
其中C为任意常数;
R
(2)方程两边同时乘以因子e−
1 1−x2
dx,由此得到方程的通解为
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1.给定方程组x = x x= (*)a)试验证u(t)= ,v(t)= 分别是方程组(*)的满足初始条件u(0)= , v(0)= 的解.b)试验证w(t)=c u(t)+c v(t)是方程组(*)的满足初始条件w(0)= 的解,其中是任意常数.解:a) u(0)= =u (t)= = u(t)又v(0)= =v (t)= = = v(t)因此u(t),v(t)分别是给定初值问题的解.b) w(0)= u(0)+ u(0)= + =w (t)= u (t)+ v (t)= +=== w(t)因此w(t)是给定方程初值问题的解.2. 将下面的初值问题化为与之等价的一阶方程组的初值问题:a) x +2x +7tx=e ,x(1)=7, x (1)=-2b) x +x=te ,x(0)=1, x (0)=-1,x (0)=2,x (0)=0c)x(0)=1, x (0)=0,y(0)=0,y (0)=1解:a)令x =x, x = x , 得即又x =x(1)=7 x (1)= x (1)=-2于是把原初值问题化成了与之等价的一阶方程的初值问题:x =x(1)=其中x=.b) 令=x ===则得:且(0)=x(0)=1, = (0)=-1, (0)= (0)=2,(0)= (0)=0于是把原初值问题化成了与之等价的一阶方程的初值问题:= x(0)= , 其中x= .c) 令w =x,w =,w =y,w =y ,则原初值问题可化为:且即ww(0)= 其中w=3. 试用逐步逼近法求方程组=x x=满足初始条件x(0)=的第三次近似解.0241201 杨素玲习题5.202412—02 02412—031.试验证=是方程组x = x,x= ,在任何不包含原点的区间a 上的基解矩阵。
解:令的第一列为(t)= ,这时(t)= = (t)故(t)是一个解。
同样如果以(t)表示第二列,我们有(t)= = (t)这样(t)也是一个解。
因此是解矩阵。
又因为det =-t 故是基解矩阵。
2.考虑方程组x =A(t)x (5.15)其中A(t)是区间a 上的连续n n矩阵,它的元素为a (t),i ,j=1,2,…,na) 如果x (t),x (t),…,x (t)是(5.15)的任意n个解,那么它们的伏朗斯基行列式W[x (t),x (t),…,x (t)] W(t)满足下面的一阶线性微分方程W =[a (t)+a (t)+…+a (t)]Wb) 解上面的一阶线性微分方程,证明下面公式:W(t)=W(t )e t ,t [a,b]解:w (t)= + +…+= +…+ = +…+ 整理后原式变为(a +…+a )=(a +…+a )w(t)=(a (t)+…+a (t))w(t)b)由于w (t)=[ a (t)+…+a (t)] w(t),即=[ a (t)+…+a (t)]dt两边从t 到t积分ln -ln = 即w(t)=w(t )e ,t [a,b]3.设A(t)为区间a 上的连续n n实矩阵,为方程x =A(t)x的基解矩阵,而x= (t)为其一解,试证:a) 对于方程y =-A (t)y的任一解y= (t)必有(t) (t)=常数;b) (t)为方程y =-A (t)y的基解矩阵的充要条件是存在非奇异的常数矩阵C,使(t) (t)=C.解a)[ (t) (t)] = (t)+ (t)= (t)+ (t)A(t)又因为=-A (t) (t),所以=- (t) A(t)[ (t) (t)] =- (t) (t)A(t)+ (t) A(t) (t)=0,所以对于方程y =-A (t)y的任一解y= (t)必有(t) (t)=常数b) “”假设为方程y =-A (t)y的基解矩阵,则[ (t) (t)] = [ (t)] + (t) (t)=[- A (t) (t)] + (t) A (t) ) + (t)[ A(t) (t)]=- (t) A (t) + (t) A (t) =0,故(t) (t)=C“”若存在非奇异常数矩阵C,detc 0,使(t) (t)=C,则[ (t) (t)] = (t)+ (t)=0,故(t) (t)=- (t) (t)A(t) (t)=- (t) A(t) 所以(t)=- (t) A(t), (t)=- (t) A (t)即(t)为方程y =-A (t)y的基解矩阵4.设为方程x =Ax(A为n n常数矩阵)的标准基解矩阵(即(0)=E),证明:(t )= (t- t )其中t 为某一值.证明:(1), (t- t )是基解矩阵。
(2)由于为方程x =Ax的解矩阵,所以(t )也是x =Ax的解矩阵,而当t= t 时,(t ) (t )=E, (t- t )= (0)=E. 故由解的存在唯一性定理,得(t )= (t- t )5.设A(t),f(t)分别为在区间a 上连续的n n矩阵和n维列向量,证明方程组x =A(t)x+f(t)存在且最多存在n+1个线性无关解。
证明:设x ,x ,…x 是x =A(t)x的n个线性无关解,是x =A(t)x+f(t)的一个解,则x + , x + ,…, x + , 都是非齐线性方程的解,下面来证明它们线性无关,假设存在不全为零的常数C ,(I=1,2,…,n)使得+c =0,从而x + , x + ,…, x + , 在a 上线性相关,此与已知矛盾,因此x + , x + ,…, x + , 线性无关,所以方程组x =A(t)x+f(t)存在且最多存在n+1个线性无关解。
6、试证非齐线性微分方程组的叠加原理:的解,则是方程组的解。
证明:(1)(2)分别将代入(1)和(2)则则令即证7.考虑方程组,其中a)试验证是的基解矩阵;b)试求的满足初始条件的解。
证明:a)首先验证它是基解矩阵以表示的第一列则故是方程的解如果以表示的第二列我们有故也是方程的解从而是方程的解矩阵又故是的基解矩阵;b)由常数变易公式可知,方程满足初始条件的解而8、试求,其中满足初始条件的解。
解:由第7题可知的基解矩阵则若方程满足初始条件则有若则有9、试求下列方程的通解:a)解:易知对应的齐线性方程的基本解组为这时由公式得通解为b)解:易知对应的齐线性方程的基本解组为是方程的特征根故方程有形如的根代入得故方程有通解c)解:易知对应的齐线性方程对应的特征方程为故方程的一个基本解组为因为是对应的齐线性方程的解故也是原方程的一个解故方程的通解为10、给定方程其中f(t)在上连续,试利用常数变易公式,证明:a)如果f(t)在上有界,则上面方程的每一个解在上有界;b)如果当时,,则上面方程的每一个解(当时)。
证明:a) 上有界存在M>0,使得又是齐线性方程组的基本解组非齐线性方程组的解又对于非齐线性方程组的满足初始条件的解x(t),都存在固定的常数使得从而故上面方程的每一个解在上有界b) 时,当t>N时由a)的结论故时,原命题成立11、给定方程组(5.15)这里A(t)是区间上的连续矩阵,设是(5.15)的一个基解矩阵,n维向量函数F(t,x)在,上连续,试证明初值问题:(*)的唯一解是积分方程组(**)的连续解。
反之,(**)的连续解也是初值问题(8)的解。
证明:若是(*)的唯一解则由非齐线性方程组的求解公式即(*)的解满足(**)反之,若是(**)的解,则有两边对t求导:即(**)的解是(*)的解习题5.31、假设A是n n矩阵,试证:a) 对任意常数、都有exp( A+ A)=exp A•exp Ab) 对任意整数k,都有(expA) =expkA(当k是负整数时,规定(expA) =[(expA) ] )证明:a) ∵(A)•(A)=(A)•(A)∴exp( A+ A)= exp A•exp Ab) k>0时,(expA) =expA•expA……expA=exp(A+A+……+A)=expkAk<0时,-k>0(expA) =[(expA) ] =[exp(-A)] = exp(-A)•exp(-A)……exp(-A)=exp[(-A)(-k)]=expkA故k,都有(expA) =expkA2、试证:如果是=Ax满足初始条件=的解,那么=[expA(t-t )]证明:由定理8可知=Ф(t)Ф-1(t0) +Ф(t)又因为Ф(t)= expAt , Ф-1(t0)=( expAt0)-1= exp(-At0), f(s)=0,又因为矩阵(At)•(- At0)=(- At0)•(At)所以=[expA(t-t )]3、试计算下面矩阵的特征值及对应的特征向量a)b)c)d)解:a)det(E-A)= =( -5)( +1)=0∴=5, =-1对应于=5的特征向量u= , ( )对应于=-1的特征向量v= , ( )b) det(E-A)=( +1)( +2)( -2)=0∴=-1,=2,=-2对应于=-1的特征向量u1=,(0 )对应于=2的特征向量u2=,()对应于=-2的特征向量u3=,()c)det(E-A)= =( +1)2( -3)=0∴=-1(二重),=3对应于=-1(二重)的特征向量u=,(0 )对应于=3的特征向量v=, ()d) det(E-A)= =( +3)( +1)( +2)=0∴=-1,=-2,=-3对应于=-1的特征向量u1=,(0 )对应于=-2的特征向量u2=,()对应于=-3的特征向量u3=,()4、试求方程组=Ax的一个基解矩阵,并计算expAt,其中A为:a)b)c)d)解:a)det(E-A)=0得=,=-对应于的特征向量为u=,(0 )对应于的特征向量为v=,()∴u=,v=是对应于,的两个线性无关的特征向量Ф(t)= 是一个基解矩阵ExpAt=b) 由det(E-A)=0得=5,=-1解得u=,v=是对应于,的两个线性无关的特征向量则基解矩阵为Ф(t)=Ф(0)=Ф-1(0)=则expAt=Ф(t) Ф-1(0)=c)由det(E-A)=0得=2,=-2,=-1解得基解矩阵Ф(t)=Ф-1(0)=则expAt=Ф(t) Ф-1(0)=d)由det(E-A)=0得=-3,=2+,=2-解得基解矩阵Ф(t)=则expAt=Ф(t) Ф-1(0)=5、试求方程组=Ax的基解矩阵,并求满足初始条件解:a)由第4题(b)知,基解矩阵为所以b)由第4题(d)知,基解矩阵为Ф(t)=所以c) 由3(c)可知,矩阵A的特征值为=3,=-1(二重)对应的特征向量为u1=,u2=∴=+解得=6、求方程组=Ax+f(t)的解:解:a)令=Ax的基解矩阵为Ф(t)解得Ф(t)=,则Ф-1(t)=Ф-1(0)=求得=b)由det(E-A)=0得=-1,=-2,=-3设对应的特征向量为v1,则(E-A)v1=0,得v1=取v1=,同理可得v2 =,v3=则Ф(t)=从而解得c)令=Ax的基解矩阵为Ф(t)由det(E-A)=0得=1,=2解得对应的基解矩阵为Ф(t)=∴Ф-1(t)=从而Ф-1(0)=∴7、假设m不是矩阵A的特征值。