二阶线性常微分方程的幂级数解法 从微分方程学中知道，在满足某些条件下，可以用幂级数来表示一个函数。因此，自然想到，能否用幂级数来表示微分方程的解呢？ 例1、求方程 ''0y xy -=的通解 解：设2012n n y a a x a x a x =+++++…… 为方程的解，这里(0,1,2,,,)i a i n =……是待定常系数，将它对x 微分两次，有 ''212312132(1)(1)n n n n y a a x n n a x n na x --+=?+?++-+++ 将y ，'y 的表达式代入方程，并比较的同次幂的系数，得到 x -∞<<∞2210a ?=，30320,a a ?-= 41430,a a ?-= 52540,a a ?-= 或一般的可推得 32356(31)3k a a k k = ?????-? ， 1 3134673(31) k a a k k += ??????+ ， 320k a += 其中1a ，2a 是任意的，因而代入设的解中可得： 36347 01[1][] 2323562356(31)33434673(31) n x x x x x y a a x n n n n =+++++++++?????????-????????+ 这个幂级数的收敛半径是无限大的，因而级数的和（其中包括两个任意常数0a 及1a ）便是所要求的通解。
例6 求方程'''240y xy y --=的满足初值条件(0)0y =及'(0)1y =的解。 解 设级 2012n n y a a x a x a x =+++++……为方程的解。首先，利用初值 条件，可以得到 00a =， 11a =， 因而 2323'2123''223123232(1)n n n n n n y x a x a x a x y a x a x na x y a a x n n a x --=+++++=+++++=+?++-+ 将y ，'y ，''y 的表达式带入原方程，合并x 的各同次幂的项，并令各项系数等于零，得到 21422 0,1,0,,,1 n n a a a a a n -==== - 因而 567891111 ,0,,0,,2!63!4! a a a a a = ===== 最后得 21111 (1)!! k a k k k += ?=- , 20k a =, 对一切正整数k 成立。 将i a (0,1,2,)i = 的值代回2012n n y a a x a x a x =+++++……就得到 521 3 2!! k x x y x x k +=+++++ 2 422 (1),2!! k x x x x x xe k =++++ += 这就是方程的满足所给初值条件的解。 是否所有方程都能按以上方式求出其幂级数解？或者说究竟方程应该满足什么条件才能保证它的解可用幂级数来表示呢？级数的
目录 1 引言 (1) 2 二阶常系数常微分方程的几种解法 (1) 2.1 特征方程法 (1) 2.1.1 特征根是两个实根的情形 (2) 2.1.2 特征根有重根的情形 (2) 2.2 常数变异法 (4) 2.3 拉普拉斯变化法 (5) 3 常微分方程的简单应用 (6) 3.1 特征方程法 (7) 3.2 常数变异法 (9) 3.3 拉普拉斯变化法 (10) 4 总结及意义 (11) 参考文献 (12)
二阶常微分方程的解法及其应用 摘要：本文通过对特征方程法、常数变易法、拉普拉斯变换法这三种二阶常系数常微分方程解法进行介绍，特别是其中的特征方程法分为特征根是两个实根的情形和特征根有重根的情形这两种情况，分别使用特征值法、常数变异法以及拉普拉斯变换法来求动力学方程，现今对于二阶常微分方程解法的研究已经取得了不少成就，尤其在二阶常系数线性微分方程的求解问题方面卓有成效。应用常微分方程理论已经取得了很大的成就，但是，它的现有理论也还远远不能满足需要，还有待于进一步的发展，使这门学科的理论更加完善。 关键词：二阶常微分方程;特征分析法；常数变异法；拉普拉斯变换
METHODS FOR TWO ORDER ORDINARY DIFFERENTIAL EQUATION AND ITS APPLICATION Abstract:This paper introduces the solution of the characteristic equation method, the method of variation of parameters, the Laplasse transform method the three kind of two order ordinary differential equations with constant coefficients, especially the characteristic equation method which is characteristic of the root is the two of two real roots and characteristics of root root, branch and don't use eigenvalue method, method of variation of constants and Laplasse transform method to obtain the dynamic equation, the current studies on solution of ordinary differential equations of order two has made many achievements, especially in the aspect of solving the problem of two order linear differential equation with constant coefficients very fruitful. Application of the theory of ordinary differential equations has made great achievements, however, the existing theory it is still far from meeting the need, needs further development, to make the discipline theory more perfect. Keywords:second ord er ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform 1 引言 数学发展的历史告诉我们，300年来数学分析是数学的首要分支，而微分方程
第4章 一阶线性微分方程组 一 内容提要 1． 基本概念 一阶微分方程组：形如 ??? ????? ???===) ,,,,( ),,,,(),,,,(2121222111 n n n n n y y y x f dx dy y y y x f dx dy y y y x f dx dy ΛΛΛΛΛ （3.1） 的方程组，（其中n y y y ,,,21Λ是关于x 的未知函数）叫做一阶微分方程组。 若存在一组函数)(,),(),(21x y x y x y n Λ使得在[a,b]上有恒等式 ),,2,1))((,),(),(,() (21n i x y x y x y x f dx x dy n i i ΛΛ==成立，则 )(,),(),(21x y x y x y n Λ称为一阶微分方程组(3.1)的一个解 含有n 任意常数n C C C ,,,21Λ的解 ?????? ?===) ,,,,( ),,,,(),,,,(21321222111n n n n C C C x y C C C x y C C C x y ΛΛΛΛΛ??? 称为(3.1)通解。如果通解满方程组 ???????=Φ=Φ=Φ0 ),,,,,,,,( 0),,,,,,,,(0),,,,,,,,(21212121221211n n n n n n n C C C y y y x C C C y y y x C C C y y y x ΛΛΛΛΛΛΛΛ 则称这个方程组为（3.1）的通积分。 满足初始条件,)(,,)(,)(0020021001n n y x y y x y y x y ===Λ的解，叫做初值问题的解。
二阶常系数线性微分方程 一、二阶常系数线形微分方程的概念 形如 )(x f qy y p y =+'+'' (1) 的方程称为二阶常系数线性微分方程.其中p 、q 均为实数,)(x f 为已知的连续函数. 如果0)(≡x f ,则方程式 (1)变成 0=+'+''qy y p y (2) 我们把方程(2)叫做二阶常系数齐次线性方程,把方程式(1)叫做二阶常 系数非齐次线性方程. 本节我们将讨论其解法. 二、二阶常系数齐次线性微分方程 1．解的叠加性 定理1 如果函数1y 与2y 是式(2)的两个解, 则2211y C y C y +=也是 式(2)的解,其中21,C C 是任意常数. 证明 因为1y 与2y 是方程(2)的解,所以有 0111 =+'+''qy y p y 0222 =+'+''qy y p y 将2211y C y C y +=代入方程(2)的左边,得 )()()(22112211221 1y C y C q y C y C p y C y C ++'+'+''+'' =0)()(2222111 1=+'+''++'+''qy y p y C qy y p y C 所以2211y C y C y +=是方程(2)的解. 定理1说明齐次线性方程的解具有叠加性. 叠加起来的解从形式看含有21,C C 两个任意常数,但它不一定是方程式(2)的通解. 2.线性相关、线性无关的概念
设,,,,21n y y y 为定义在区间I 内的n 个函数,若存在不全为零的常数 ,,,,21n k k k 使得当在该区间内有02211≡+++n n y k y k y k , 则称这n 个函数在区间I 内线性相关,否则称线性无关. 例如 x x 22sin ,cos ,1在实数范围内是线性相关的,因为 0sin cos 12 2≡--x x 又如2,,1x x 在任何区间(a,b)内是线性无关的,因为在该区间内要使 02321≡++x k x k k 必须0321===k k k . 对两个函数的情形,若=21y y 常数, 则1y ,2y 线性相关,若≠2 1y y 常数, 则1y ,2y 线性无关. 3.二阶常系数齐次微分方程的解法 定理 2 如果1y 与2y 是方程式(2)的两个线性无关的特解,则 212211,(C C y C y C y +=为任意常数)是方程式(2)的通解. 例如, 0=+''y y 是二阶齐次线性方程,x y x y cos ,sin 21==是它的 两个解,且≠=x y y tan 2 1常数,即1y ,2y 线性无关, 所以 x C x C y C y C y cos sin 212211+=+= ( 21,C C 是任意常数)是方程0=+''y y 的通解. 由于指数函数rx e y =(r 为常数)和它的各阶导数都只差一个常数因子, 根据指数函数的这个特点,我们用rx e y =来试着看能否选取适当的常数r , 使rx e y =满足方程(2).
第七章 一阶线性偏微分方程 研究对象 一阶线性齐次偏微分方程 0),,,(),,,() ,,,(2122121211=??++??+??n n n n n x u x x x X x u x x x X x u x x x X 1基本概念 1） 一阶线性齐次偏微分方程 形如 0),,,(),,,(),,,(2122121211=??++??+??n n n n n x u x x x X x u x x x X x u x x x X （7.1） 的方程，称为一阶线性齐次偏微分方程，其中n x x x ,,,21 是自变量，u 是n x x x ,,,21 的未知函数，n X X X ,,,21 是域n R D ?内的已知函数，并设n X X X ,,,21 在域D 内不同时为零。 2） 一阶拟线性偏微分方程 形如 );,,,();,,,();,,,(21211211z x x x Z x z z x x x Y x z z x x x Y n n n n n =??++?? （7.2） 的方程，称为一阶拟线性偏微分方程，其中Z Y Y Y n ;,,,21 是1+n 个变元z x x x n ;,,,21 的已知函数。n Y Y Y ,,,21 在其定义域1+?'n R D 内不同时为零。 所谓“拟线性”是指方程仅对未知函数的各个一阶偏导数是线性的，以下总设n Y Y Y ,,,21 和Z 在域D '内连续可微。 3） 特征方程组 常微分方程组 n n X dx X dx X dx === 2211 （7.3） 称为一阶线性齐次偏微分方程（7.1）的特征方程组。 常微分方程组
习题2.2 求下列方程的解 1．dx dy =x y sin + 解： y=e ?dx (?x sin e ?-dx c dx +) =e x [- 2 1e x -(x x cos sin +)+c] =c e x -21 (x x cos sin +)是原方程的解。 2．dt dx +3x=e t 2 解：原方程可化为： dt dx =-3x+e t 2 所以：x=e ?-dt 3 (?e t 2 e -? -dt 3c dt +) =e t 3- (5 1e t 5+c) =c e t 3-+5 1e t 2 是原方程的解。 3．dt ds =-s t cos +21t 2sin 解：s=e ?-tdt cos (t 2sin 2 1?e dt dt ?3c + ) =e t sin -(?+c dt te t t sin cos sin ) = e t sin -(c e te t t +-sin sin sin ) =1sin sin -+-t ce t 是原方程的解。 4． dx dy n x x e y n x =- ， n 为常数. 解：原方程可化为：dx dy n x x e y n x += )(c dx e x e e y dx x n n x dx x n +??=?- )(c e x x n += 是原方程的解.
5． dx dy +1212--y x x =0 解：原方程可化为：dx dy =-1212+-y x x ?=-dx x x e y 21 2(c dx e dx x x +?-221) )21(ln 2+=x e )(1 ln 2?+--c dx e x x =)1(1 2 x ce x + 是原方程的解. 6． dx dy 234xy x x += 解：dx dy 234xy x x += =23y x +x y 令 x y u = 则 ux y = d x d y =u dx du x + 因此：dx du x u +=2u x 21u dx du = dx du u =2 c x u +=33 1 c x x u +=-33 （*） 将x y u =带入 （*）中 得：3433cx x y =-是原方程的解.
二阶非线性常微分方程的打靶法 1.问题： 试用打靶法求二阶非线性常微分方程亮点边值的数值解： 要求用Matlab 编程计算，请给出一些例子，验证你的算法与程序的正确性。 2.打靶法分析： 非线性打靶法: 非线性打靶法的基本原理是将两点边值问题(1)转化为下面形式的初值问题 令z = y ′,将上述二阶方程降为一阶方程组
3.Matlab 源代码： 创建M 文件： function ys=dbf(f,a,b,alfa,beta,h,eps) ff=@(x,y)[y(2),f(y(1),y(2),x)]; xvalue=a:h:b;%x取值范围 n=length(xvalue) s0=a-0.01;%选取适当的s的初值 x0=[alfa,s0];%迭代初值 flag=0;%用于判断精度 y0=rk4(ff,a,x0,h,a,b); if abs(y0(1,n)-beta)<=eps flag=1; y1=y0; else s1=s0+1; x0=[alfa,s1]; y1=rk4(ff,a,x0,h,a,b); if abs(y1(1,n)-beta)<=eps flag=1; end end if flag~=1 while abs(y1(1,n)-beta)>eps s2=s1-(y1(1,n)-beta)*(s1-s0)/(y1(1,n)-y0(1,n)); x0=[alfa,s2]; y2=rk4(ff,a,x0,h,a,b); s0=s1; s1=s2;
y0=y1; y1=y2; end end xvalue=a:h:b; yvalue=y1(1,:); ys=[xvalue',yvalue']; function x=rk4(f,t0,x0,h,a,b)%rung-kuta法求每个点的近似值（参考大作业一）t=a:h:b;%迭代区间 m=length(t);%区间长度 t(1)=t0; x(:,1)=x0;%迭代初值 for i=1:m-1 L1=f(t(i),x(:,i)); L2=f(t(i)+h/2,x(:,i)'+(h/2)*L1); L3=f(t(i)+h/2,x(:,i)'+(h/2)*L2); L4=f(t(i)+h,x(:,i)'+h*L3); x(:,i+1)=x(:,i)'+(h/6)*(L1+2*L2+2*L3+L4); end 4.举例 求二阶非线性方程的边值问题： 在matlab 控制台中输入： f=@(x,y,z)(x^2+z*x^2); x0l=0; x0u=2*exp(-1); alfa=0; beta=2; h=0.01 dbf(f,x0l,x0u,y0l,y0u,h,1e-6); >> y=ans(:,2); x=ans(:,1); >> plot(x,y,'-r') >> 结果：
一阶常系数线性齐次微分方程组的求解 【模型准备】一只虫子在平面直角坐标系内爬行. 开始时位于点P 0(1, 0)处. 如果知道虫子在点P (x , y )处沿x 轴正向的速率为4x - 5y , 沿y 轴正向的速率为2x - 3y . 如何确定虫子爬行的轨迹的参数方程? 图31 虫子爬行的轨迹 【模型假设】设t 时刻虫子所处位置的坐标为(x (t ), y (t )). 【模型构成】由已知条件和上述假设可知 d 45,d d 23,d x x y t y x y t ?=-????=-??而且(x (0), y (0)) = (1, 0). 现要由此得出虫子爬行的轨迹的参数方程. 【模型求解】令A =4523-?? ?-?? , 则|λE -A | =4523λλ--+= (λ+1)(λ-2). 可见A 的特征值为λ1 = -1, λ2 = 2. (-E -A )x = 0的一个基础解系为: ξ1 = (1, 1)T ; (2E -A )x = 0的一个基础解系为: ξ2 = (5, 2)T . 令P = (ξ1, ξ2), 则P -1AP =1002-?? ??? . 记X =x y ?? ???, Y =u v ?? ??? , 并且作线性变换X = PY , 则Y = P -1X , d d t Y = P -1d d t X = P -1AX = P -1APY =1002-?? ??? Y , 即 d d d d u t v t ?? ???=1002-?? ???u v ?? ??? , 故u = c 1e -t , v = c 2e 2t , 即Y =122t t c e c e -?? ??? . 因而 12c c ?? ??? = Y |t =0 = P -1X |t =0 =2/35/31/31/3-?? ?-??10?? ???=2/31/3-?? ???. 于是 x y O 1 何去何从?
第六章 二阶线性偏微分方程的分类与小结 一 两个自变量的二阶线性方程 1 方程变换与特征方程 两个自变量的二阶线性偏微分方程总表示成 f cu u b u b u a u a u a y x yy xy xx =+++++212212112 ① 它关于未知函数u 及其一、二阶偏导数都是线性的，其中f u c b b a a a ，，，，，，，21221211都是自变量y x ,的已知函数，假设它们的一阶偏 导数在某平面区域D 内都连续，而且221211a a a ，，不全为0 。 设),(000y x M 是D 内给定的一点，考虑在0M 的领域内对方程进行简化。取自变量变换 ),(y x ξξ=,),(y x ηη= 其中它们具有二连续偏导数，而且在0M 处的雅可比行列式。 = ??),(),(y x ηξy x y x ηηξξ =x y y x ηξηξ- 根据隐函数存在定理，在0M 领域内存在逆变换, ),(ηξx x =,),(ηξy y = 因为 x x x u u u ηξξξ+=，y y y u u u ηξξξ+=
xx xx x x x x xx u u u u u u ηξηηξξηξηηξηξξ++++=222 yy yy y y y y yy u u u u u u ηξηηξξηξηηξηξξ++++=222 xy xy y x x y y x x x xy u u u u u u ηξηηηξηξξξηξηηξηξξ+++++=)( 将代入①使其变为 F Cu u B u B u A u A u A =+++++ηξηηξηξξ212212112 经过变换后，方程的阶数不会升高，由变换的可逆性，方程的阶数也不会降低，所以221211,,A A A 不全为0。并可验证 222112122211212))((x y y x a a a A A A ηξηξ--=- 这表明，在可逆变换下2 22112 12A A A -与22112 12 a a a -保持相同的正负号。 定理 在0M 的领域内，不为常数的函数),(y x ?是偏微分方程022*******=++y y x x a a a ????之解的充分必要条件是： C y x ≡),(?是常微分方程的 0)(2)(22212211=++dx a dxdy a dy a 通解。 2 方程的类型及其标准形式 根据以上结论简化方程的问题归结为寻求其特征曲线。为此将特征方程分解成两个方程： 11 22 11 2 12 12 a a a a a dx dy -+=，11 22 11 2 12 12 a a a a a dz dy --= (1) 若在0M 的邻域内022112 12>-a a a 时，方程可以化为
南京师范大学泰州学院 毕业论文（设计） （一六届） 题目：二阶常微分方程的解法 院（系、部）：数学科学与应用学院 专业：数学与应用数学 姓名：潘陆 学号08120146 指导教师：刘陆军 南京师范大学泰州学院教务处制
摘要：本文主要是介绍了二阶常微分方程众多解法中的三种，分别为特征方程法，拉普拉斯变换法和常数变易法，研究并讨论了二阶常微分方程在特征方程法中特征方程根为实根，复根和重根的情形。我们选用了弹簧振子系统的振子运动，用这三种不同的方法来解决该问题。 关键词：二阶常微分方程；特征根法；常数变易法；拉普拉斯变换 Abstract:The main purpose of this paper is the second-order ordinary many differential equation solution of three, respectively as the characteristic equation method, Laplace transform method and variation of constants method, study and discuss the second-order often differential equation in the characteristic equation of the roots of the characteristic equation for real roots, complex roots and root weight. We choose the spring oscillator the oscillator motion, these three different methods to solve the problem. Keywords: second order ordinary differential equation; Characteristic analysis; constant variation method; Laplasse transform
附录A 线性常微分方程 本课程的研究内容与常微分方程理论有非常密切的联系，因此在本附录里，我们将对线性常微分方程的知识——包括解的存在性、解的结构和求解方法做一些回顾和总结。 把包含未知函数和它的j 阶导数()j y (的方程称为常微分方程。线性常微分方程的标准形式 ()(1)110()()'()()n n n y p x y p x y p x y f x --++++=L （A.1） 其中n 称为方程的阶数，()j p x 和()f x 是给定的函数。可微函数()y y x =在区间 I 上满足方程（A.1），则称其为常微分方程（A.1）在 I 上的一个解。，()f x 称为方程（A.1）的自由项，当自由项()0f x ≡时方程（A.1）称为是齐次方程，否则称为非齐次方程。一般来说常微分方程的解是不唯一的，我们将方程的全部解构成的集合称为解集合，解集合中全部元素的一个通项表达式称为方程的通解，而某个给定的解称为方程的特解。 在本附录里，我们重点介绍一阶和二阶常微分方程的相关知识。 A.1 一阶线性常微分方程 一阶线性常微分方程表示为 '()()y p x y f x x I +=∈，. （A.2） 当()0f x ≡，方程退化为 '()0y p x y +=, （A.3） 假设()y x 不恒等于零，则上式等价于 而()'ln 'y y y =，从而（A.3）的通解为 ()d ()p x x y x Ce -?= （ A.4） 对于非齐次一阶线性常微分方程（A.2），在其两端同乘以函数()d p x x e ?
注意到上面等式的左端 因此有 两端积分 其中C 是任意常数。进一步有 综上有如下结论 定理A.1 假设()()p x f x I 和在上连续，则一阶线性非齐次常微分方程（A.1）的通解具有如下形式 ()d ()d ()d ()()d p x x p x x p x x y x Ce e e f x x --? ??=+?‘ （A.5） 其中C 是任意常数。 观察（A.4）式和（A.5）式，我们发现一阶线性非齐次常微分方程（A.1）的解等于 一阶线性齐次常微分方程（ A.2）的通解()d p x x Ce -?加上函数()d ()d *()()d p x x p x x y x e e f x x -??=?。容易验证，*()y x 是方程（A.1）的一个特解。这符合线性方程解的结构规律。 例1 求解一阶常微分方程 解 此时()2()1p x f x =-=，,由（A.5）式，解为 其中C 是任意常数。 A.2 二阶线性常微分方程 将具有以下形式的方程 "()'()()y p x y q x y f x x I ++=∈，， （A.6） 称为二阶线性常微分方程，其中(),(),()p x q x f x 都是变量x 的已知连续函数。称 "()'()0y p x y q x y x I ++=∈，， （A.7） 为与（A.6）相伴的齐次方程. A ．2．1 二阶线性微分方程解的结构 首先讨论齐次方程（A.7）解的结构。
一阶偏微分方程基本知识 这一章我们来讨论一阶线性偏微分方程和一阶拟线性偏微分方程的解法，因为它们都可以化为常微分方程的首次积分问题，所以我们先来介绍常微分方程的首次积分。 1一阶常微分方程组的首次积分 1.1首次积分的定义 从第三章我们知道，n 阶常微分方程 ()()() 1,,'',',-=n n y y y x f y ， （ 1.1） 在变换 ( ) 1'12,,,,n n y y y y y y -=== （ 1.2） 之下，等价于下面的一阶微分方程组 ()()()1 112221212,,,,,,,,,,,,,,. n n n n n dy f x y y y dx dy f x y y y dx dy f x y y y dx ?=?? ?=???? ?=? ? （ 1.3） 在第三章中，已经介绍过方程组（ 1.3）通解的概念和求法。但是除了常 系数线性方程组外，求一般的（ 1.3）的解是极其困难的。然而在某些情况下，可以使用所谓“可积组合”法求通积分，下面先通过例子说明“可积组合”法，然后介绍一阶常微分方程组“首次积分”的概念和性质，以及用首次积分方法来求解方程组（ 1.3）的问题。先看几个例子。 例1 求解微分方程组 ()()22221,1.dx dy y x x y x y x y dt dt =-+-=--+- （ 1.4） 解：将第一式的两端同乘x ，第二式的两端同乘y ，然后相加，得到 ()() 12222-++-=+y x y x dt dy y dt dx x ， ()()()2222221 12 d x y x y x y dt +=-++-。 这个微分方程关于变量t 和()22x y +是可以分离，因此不难求得其解为 122 2221C e y x y x t =+-+， （ 1.5） 1C 为积分常数。（ 1.5）叫做（ 1.4）的首次积分。
习 题 6 —— 1 1．求出齐次线性微分方程组 y t A dt dy )(=的通解，其中A （t ）分别为：（1）???? ??=1011)(t A ；（2）???? ??-=0110)(t A ；（3）???? ? ??=000010100)(t A 。 解 （1）方程组的分量形式为： 211y y dt dy += ，22y dt dy = 从后一式容易求出2y 的通解为 t ke y =2 ，其中K 为任意常数，可分别取02=y 和 t e y =2，代入前一式得到两个相应的特解，t e y =1和 t te y =2这样就求得方程组的一个解矩阵为 ()0t t t e te t e ??Φ= ??? 又 2det ()0t t e Φ=≠ 。因此，)(t Φ是方程组的一个基解矩阵，根据定理6.1 ,方程的通解为 ??? ? ??+???? ??=???? ??t t t e te c e c y y 21210 （2）方程的分量形式为 ?????-==1221y dt dy y dt dy 由①、②可和 21120d y y dt += 由观察法知，t y cos 1=，t y sin 1=为此方程的两个特解，将其代入②式可得两个相应的特解，将其代入②式可得两个相应的特解：2sin y t =-，2cos y t =。这样就求得方 程组的一个解矩阵为 cos int ()int cos t s t s t ??Φ= ?-?? 又 []01)(det ≠=Φ=t ，因此)(t Φ中方程组的一个基解矩阵。故方程组的通解为 1122cos int int cos y t s c c y s t ??????=+ ? ? ?-???? ?? （3）程组的分量形式为：?????='='='13 2231y y y y y y ① ② ① ② ③
习题 1 求方程dx dy =x+y 2通过点(0,0)的第三次近似解； 解： 取0)(0=x ? 20020012 1)()(x xdx dx y x y x x x ==++=??? 522200210220 121])21([])([)(x x dx x x dx x x y x x x +=+=++=???? dx x x x y x x ])20 121([)(252003+++=?? = 118524400 1160120121x x x x +++ 2 求方程dx dy =x-y 2通过点(1,0)的第三次近似解； 解： 令0)(0=x ? 则 20020012 1)()(x xdx dx y x y x x x ==-+=??? 522200210220 121])21([])([)(x x dx x x dx x x y x x x -=-=-+=???? dx x x x y x x ])20 121([)(252003--+=?? =118524400 1160120121x x x x -+- 3 题 求初值问题： ?????=-=0 )1(2y x dx dy R ：1+x ≤1，y ≤1 的解的存在区间，并求解第二次近似解，给出在解的存在空间的误差估计； 解： 因为 M=max{22y x -}=4 则h=min(a,M b )=4 1 则解的存在区间为0x x -=)1(--x =1+x ≤4 1 令 )(0X ψ=0 ; )(1x ψ=y 0+?-x x x 0)0(2dx=31x 3+31;
)(2x ψ =y 0+])3131([2132?-+-x x x dx=31x 3-9x -184x -637x +4211 又 y y x f ??),(2≤=L 则：误差估计为：)()(2x x ψ-ψ≤32 2 )12(*h L M +=2411 4 题 讨论方程：31 23y dx dy =在怎样的区域中满足解的存在唯一性定理的条件， 并求通过点（0，0）的一切解； 解：因为y y x f ??),(=3221-y 在y 0≠上存在且连续； 而312 3y 在y 0φσ≥上连续 由 3123y dx dy =有：y =（x+c ）23 又 因为y(0)=0 所以：y =x 2 3 另外 y=0也是方程的解； 故 方程的解为：y =?????≥00023πx x x 或 y=0; 6题 证明格朗瓦耳不等式： 设K 为非负整数，f(t)和g(t)为区间βα≤≤t 上的连续非负函数，
第二章 二阶线性偏微分方程的分类 1．把下列方程化为标准形式： （1）02=+++++u cu bu au au au y x yy xy xx 解：因为 02 22112 12=?-=-a a a a a a 所以该方程是抛物型方程，其特征方程为 12 2 =-± =a a a a dx dy 。 它只有一族实的特征线 c x y =- 在这种情况下，我们设x y -=ξ，x =η（或令y =η，总之，此处η是与ξ无关的任一函数，当然宜取最简单的函数形式x =η或y =η）。 方法一：用抛物型方程的标准形式 ][12122 F Cu u B u B A +++- =ηξηηη 先算出： ? ??? ? ? ?? ? ? ?-====?+?+?+?+?=++++=?+-+?+?+?=++++==?+?+=++=b c C b c b a a a b b a a a B c b a a a b b a a a B a a a a a a a A y x yy xy xx y x yy xy xx y y x x 0F ,1010020 2 1)1(0020 2 002 2212212112 2122121112 221221122ηηηηηξξξξξηηηη ∴])[(1 u bu u c b a u +++--=ηξηη 即 01=+ + -+ u a u a b u a b c u ηξηη 方法二：应用特征方程，作自变量变换，求出 ??? ??=+-=+-=+--==+-= ,2 ,ξξηξξξηηξηξξηηηξξηξξξηξu u u u u u u u u u u u u u u u u u yy xy xx y x 代入原方程得，0)(=++-+u bu u b c au ηξξη
常微分方程试卷及答案 Company Document number：WTUT-WT88Y-W8BBGB-BWYTT-19998
2010-2011 学年第 二 学期常微分方程考试 AB 卷答案 理学 院 年级 信息与计算科学 专业 填空题（每题4分，共20分） 1. 形如)()('x Q y x P y += （)(),(x Q x P 连续）的方程是 一阶线性微分 方程，它的通解为?? ? ???+?-? =c dx dx x P e x Q dx x P e y )()()( . 2. 形如0y y '''-=的方程是 3 阶__齐次__(“齐次”还是”非齐次”)___常__系数的微分方程,它的特征方程为310λ-=. 3. 形如1 11 111 0n n n n n n n n d y d y dy x a x a x a y dx dx dx ----++++=的方程为 欧拉 方程, 可通过变换t x e =把它转化成常系数方程. 4. 2 (1)0,y dx x dy ++= 满足初始条件：x =0, y =1的特解1 1ln 1y x = ++ 5．5.微分方程0000(,),(),:,dy f x y y x y R x x a y y b dx ==-≤-≤满足的解存在且唯一的条件是: (,)f x y 在R 上连续且满足利普希茨条件 一、下列微分方程的解（每题5分，共30分） 1． dx dy =2 )(1y x + 解：令x+y=u ,则dx dy =dx du -1 (3) dx du -1=21 u u-arctgu=x+c y-arctg(x+y)=c. (5) 2．()()053243 =+++xdy ydx y xdy ydx x 解：两边同乘以y x 2得： ()() 0532******* =+++ydy x dx y x ydy x dx y x (3)
第八章 8.4讲 第四节 二阶常系数线性微分方程 一、二阶常系数线形微分方程的概念 形如 )(x f qy y p y =+'+'' (1) 的方程称为二阶常系数线性微分方程.其中p 、q 均为实数,)(x f 为已知的连续函数. 如果0)(≡x f ,则方程式 (1)变成 0=+'+''qy y p y (2) 我们把方程(2)叫做二阶常系数齐次线性方程,把方程式(1)叫做二阶常系数非齐次线性方程. 本节我们将讨论其解法. 二、二阶常系数齐次线性微分方程 1．解的叠加性 定理1 如果函数1y 与2y 是式(2)的两个解, 则2211y C y C y +=也是式(2)的解,其中21,C C 是任意常数. 证明 因为1y 与2y 是方程(2)的解,所以有 0111 =+'+''qy y p y 0222 =+'+''qy y p y 将2211y C y C y +=代入方程(2)的左边,得 )()()(22112211221 1y C y C q y C y C p y C y C ++'+'+''+'' =0)()(2222111 1=+'+''++'+''qy y p y C qy y p y C 所以2211y C y C y +=是方程(2)的解. 定理1说明齐次线性方程的解具有叠加性. 叠加起来的解从形式看含有21,C C 两个任意常数,但它不一定是方程式(2)
的通解. 2.线性相关、线性无关的概念 设,,,,21n y y y 为定义在区间I 内的n 个函数,若存在不全为零的常数 ,,,,21n k k k 使得当在该区间内有02211≡+++n n y k y k y k , 则称这n 个函数在区间I 内线性相关,否则称线性无关. 例如 x x 2 2 sin ,cos ,1在实数范围内是线性相关的,因为 0sin cos 12 2≡--x x 又如2,,1x x 在任何区间(a,b)内是线性无关的,因为在该区间内要使 02321≡++x k x k k 必须0321===k k k . 对两个函数的情形,若 =21y y 常数, 则1y ,2y 线性相关,若≠2 1y y 常数, 则1y ,2y 线性无关. 3.二阶常系数齐次微分方程的解法 定理 2 如果1y 与2y 是方程式(2)的两个线性无关的特解,则 212211,(C C y C y C y +=为任意常数)是方程式(2)的通解. 例如, 0=+''y y 是二阶齐次线性方程,x y x y cos ,sin 21==是它的两个解,且 ≠=x y y tan 2 1 常数,即1y ,2y 线性无关, 所以 x C x C y C y C y cos sin 212211+=+= ( 21,C C 是任意常数)是方程0=+''y y 的通解. 由于指数函数rx e y =(r 为常数)和它的各阶导数都只差一个常数因子,
一阶线性微分方程及伯 努利介绍 Company number：【0089WT-8898YT-W8CCB-BUUT-202108】
第三节 一阶线性微分方程 内容要点 一、一阶线性微分方程 形如 )()(x Q y x P dx dy =+ 的方程称为一阶线性微分方程. 其中函数)(x P 、)(x Q 是某一区间I 上的连续函数. 当,0)(≡x Q 方程成为 0)(=+y x P dx dy 这个方程称为一阶齐次线性方程. 相应地，方程称为一阶非齐次线性方程. 方程的通解 .)(?-=dx x P Ce y 其中C 为任意常数. 求解一阶非齐次线性微分方程的常数变易法：即在求出对应齐次方程的通解后，将通解中的常数C 变易为待定函数)(x u ，并设一阶非齐次方程通解为 一阶非齐次线性方程的通解为 [] ?-?+=?dx x P dx x P e C dx e x Q y )()()( 二、伯努利方程：形如 n y x Q y x P dx dy )()(=+ 的方程称为伯努利方程，其中n 为常数，且1,0≠n . 伯努利方程是一类非线性方程，但是通过适当的变换，就可以把它化为线性的. 事实上，在方程两端除以n y ，得 或 ),()()(1111x Q y x P y n n n =+'?--- 于是，令n y z -=1，就得到关于变量z 的一阶线性方程
)()1()()1(x Q n z x P n dx dz -=-+. 利用线性方程的求解方法求出通解后，再回代原变量，便可得到伯努利方程的通解 雅各布.伯努利（Jacob Bermoulli ，1654~1705) 伯努利瑞士数学、力学、天文学家。1654年12月27日生于瑞士巴塞尔；1705年8月16日卒于巴塞尔。 雅各布.伯努利出生于一商人世家。他的祖父是一位药商，1662年移居巴塞尔。他的父亲接过兴隆的药材生意，并成了市议会的一名成员和地方行政官。他的母亲是市议员兼银行家的女儿。雅格布在1684年一位富商的女儿结婚，他的儿子尼古拉，伯努得是艺术家，巴塞尔市议会的议员和艺术行会会长。 雅格布毕业于巴塞尔大学，1671年获艺术硕士学位。这里的艺术是指“自由艺术”，它包括算术、几何、天文学、数理音乐的基础，以及方法、修辞和雄辩术等七大门类。遵照他父亲的愿望，他又于1676年得硕士学位。同时他对数学有着浓厚的兴趣，但是他在数学上的兴趣遭到父亲的反对，他违背父亲的意愿，自学了数学和天文学。1676年，他到日内瓦做家庭教师。从1677年起，他开始在这里写内容丰富的《沉思录》。1678年雅格布进行了他第一次学习旅行，他到过法国、荷兰、英国和德国，与数学家们建立了广泛的通信联系。然后他又在法国度过了两年时光，这期间他开始研究数学问题。起初他还不知道牛顿和莱布尼兹的工作，他首先熟悉了笛卡尔的《几何学》、活利斯的《无穷的算术》以及巴罗的《几何学讲义》。他后来逐渐地熟悉了莱布尼兹的工作。1681-1682年间，他做了第二次学习旅行，接触了许多数学家和科学家。通过访问和阅读文献，丰富了他的知识，拓宽了个人的兴趣。这次旅行，他在科学上