常微分方程练习试卷 一、 填空题。 1. 方程23 2 10d x x dt +=是 阶 （线性、非线性）微分方程. 2. 方程 ()x dy f xy y dx =经变换_______，可以化为变量分离方程 . 3. 微分方程 3230d y y x dx --=满足条件(0)1,(0)2y y '==的解有 个. 4. 设常系数方程 x y y y e αβγ'''++=的一个特解*2()x x x y x e e xe =++，则此方程的系数α= ，β= ，γ= . 5. 朗斯基行列式 ()0W t ≡是函数组12(),(),,()n x t x t x t 在a x b ≤≤上线性相关的 条件. 6. 方程 22(2320)0xydx x y dy ++-=的只与y 有关的积分因子为 . 7. 已知 ()X A t X '=的基解矩阵为()t Φ的，则()A t = . 8. 方程组 20'05??=???? x x 的基解矩阵为 ． 9.可用变换 将伯努利方程 化为线性方程. 10 .是满足方程 251y y y y ''''''+++= 和初始条件 的唯一解. 11.方程 的待定特解可取 的形式: 12. 三阶常系数齐线性方程 20y y y '''''-+=的特征根是 二、 计算题 1.求平面上过原点的曲线方程, 该曲线上任一点处的切线与切点和点(1,0)的连线相互垂直. 2．求解方程13 dy x y dx x y +-=-+. 3. 求解方程 222()0d x dx x dt dt += 。 4．用比较系数法解方程. . 5．求方程 sin y y x '=+的通解. 6．验证微分方程 22(cos sin )(1)0x x xy dx y x dy -+-=是恰当方程，并求出它的通解.
摘要.................................................... 错误!未定义书签。 1.何谓奇解.............................................. 错误!未定义书签。 2.奇解的产生............................................ 错误!未定义书签。 3.包络跟奇解的关系...................................... 错误!未定义书签。 4.理论上证明C-判别曲线与P-判别曲线方法................. 错误!未定义书签。 克莱罗微分方程 ..................................... 错误!未定义书签。 5.奇解的基本性质........................................ 错误!未定义书签。 定理１ ............................................. 错误!未定义书签。 定理２ ............................................. 错误!未定义书签。 定理３ ............................................. 错误!未定义书签。 6.小结.................................................. 错误!未定义书签。参考文献：.............................................. 错误!未定义书签。
2005——2006学年第二学期 常微分方程课程试卷（B） 一、填空题（每空2 分，共16分）。 1．李普希滋条件是初值问题存在唯一解的充分条件． 2. 一阶微分方程的一个特解的图像是二 维空间上的一条曲线． 3．线性齐次微分方程组Y A Y ) ( d d x x =的一个基本解组的个数不能多于n个，其中R ∈ x，n R Y∈． 4．二阶线性齐次微分方程的两个解) ( 1 x y? =,) ( 2 x y? =成为其基本解组的充要条件是线性无关． 5．方程2 sin() y xy y '' =+的通解是 6．变量可分离方程()()()()0= +dy y q x p dx y N x M的积分因子是()() x P y N 1 7．性齐次微分方程组的解组) ( , ), ( ), ( 2 1 x x x n Y Y Y 为基本解组的充分必要条件是它们的朗斯基行列式0 ) (≠ x W． 8．方程540 y y y ''' ++=的基本解组是x x e e4 ,- - 二、选择题（每小题3 分，共15分）。 9．两个不同的线性齐次微分方程组（ D ）的基本解组． (A) 一定有相同(B) 可能有相同 (C) 一定有相似(D) 没有相同 10．方程组 ? ? ? ?? ? ? + = + = y x t y y x t x 4 3 d d 2 d d 的奇点)0,0(的类型是（D ）． （A）稳定焦点（B）不稳定焦点（C）鞍点（D）不稳定结点11．方程x(y2－1)d x+y(x2－1)d y=0的所有常数解是（ C ）． (A) 1± = x(B)1± = y
(C )1±=y , 1±=x (D )1=y , 1=x 12．n 阶线性非齐次微分方程的所有解（ D ）． （A ）构成一个线性空间 （B ）构成一个1-n 维线性空间 （C ）构成一个1+n 维线性空间 （D ）不能构成一个线性空间 13．方程4d d +-=x y x y （ A ）奇解． (A) 无 (B) 有一个 (C) 有两个 (D) 可能有 三、计算题（每小题8分，共48分） 。 14．求方程 x y x y x y tan d d +=的通解 解：令x y u =，则u x u y '+='， u x u x tan d d = 当0tan ≠u 时，等号两边积分 1d tan d C x x u u +=?? C x u ln ln sin ln += 0≠C Cx x y =sin 15．求方程0d d )1(2=+--y x x y x 的通解 解：积分因子21)(x x =μ， 则 0d 1d 122=+--y x x x y x 为全微分方程．取10=x ，00=y ，于是通积分为 1012 2d d 1C y x x y x y x =+--?? 即 C x x x y =++1 16．求方程2221)(x y x y y + '-'=的通解 解:令 p y ='，得到2 2 2x xp p y +-= （*） ，两端同时关于求导，
P77 31.利用改进欧拉方法计算下列初值问题，并画出近似解的草图：dy + =t = t y y ≤ ≤ ,2 ;5.0 0,3 )0( )1(= ,1 ? dt 代码： %改进欧拉法 function Euler(t0,y0,inv,h) n=round(inv(2)-inv(1))/h; t(1)=t0; y(1)=y0; for i=1:n y1(i+1)=y(i)+h*fun(t(i),y(i)); t(i+1)=t(i)+h; y(i+1)=y(i)+1/2*h*(fun(t(i),y(i))+ fun(t(i+1),y1(i+1))) end plot(t,y,'*r') function y=fun(t,y); y=y+1; 调用：Euler(0,3,[0,2],0.5) 得到解析解：hold on; y=dsolve('Dy=y+1','(y(0)=3)','t'); ezplot(y,[0,2]) 图像：
dy y =t - t y ;2.0 t = ≤ )0( 0,5.0 ,4 )2(2= ≤ ? ,2 dt 代码： function Euler1(t0,y0,inv,h) n=round(inv(2)-inv(1))/h; t(1)=t0; y(1)=y0; for i=1:n y1(i+1)=y(i)+h*fun(t(i),y(i)); t(i+1)=t(i)+h; y(i+1)=y(i)+1/2*h*(fun(t(i),y(i))+ fun(t(i+1),y1(i+1))) end plot(t,y,'*r') function y=fun(t,y); y=y^2-4*t; 调用： Euler1(0,0.5,[0,2],0.2) 图像：
常微分方程试题 一、填空题(每小题3分，共39分) 1.常微分方程中的自变量个数是________. 2.路程函数S(t)的加速度是常数a，则此路程函数S(t)的一般形式是________. 3.微分方程=g( )中g(u)为u的连续函数，作变量变换________，方程可化为变 量分离方程. 4.微分方程F(x,y′)=0中令P=y′,若x、P平面上的曲线F(x,P)=0的参数形式 为x= (t),P=ψ(t),t为参数，则方程参数形式的通解为________. 5.方程=(x+1)3的通解为________. 6.如果函数f(x,y)连续，y= (x)是方程=f(x,y)的定义于区间x0≤x≤x0+h上，满 足初始条件 (x0)=y0的解.则y= (x)是积分方程________定义于x0≤x≤x0+h 上的连续解. 7.方程=x2+xy，满足初始条件y(0)=0的第二次近似解是________. 8.方程+a1(t) +…+a n-1(t) +a n(t)x=0 中a i(t) i=1,2,…，n是〔a,b〕上的连续函数，又x1(t),x2(t),…，x n(t)为方程n 个线性无关的解，则其伏朗斯基行列式W(t) 应具有的性质是：________. 9.常系数线性方程x(4)(t)-2x″(t)+x(t)=0的通解为________. 10.设A(t)是区间a≤t≤b上的连续n×n矩阵，x1(t),x2(t),…，x n(t)是方程组 x′=A(t)x的n个线性无关的解向量.则方程组的任一解向量x(t)均可表示为：x(t)=________的形式. 11.初值问题(t)+2x″(t)-tx′(t)+3x(t)=e-t,x(1)=1,x′(1)=2,x″(1)=3 可化为与之 等价的一阶方程组________. 12.如果A是3×3的常数矩阵，-2为A的三重特征值，则方程组x′=Ax的基 解矩阵exp A t=________. 13.方程组 的奇点类型是________. 二、计算题(共45分) 1.(6分)解方程 = . 2.(6分)解方程 x″(t)+ =0. 3.(6分)解方程 (y-1-xy)dx+xdy=0. 4.(6分)解方程
生物信息技术0801 徐聪U200812594 #include
师大学本科毕业论文（设计） 常微分方程的初等解法与求解技巧 姓名娟 院系数学与计算机科学学院 专业信息与计算科学 班级12510201 学号1251020126 指导教师王晓锋 答辩日期 成绩
常微分方程的初等解法与求解技巧 容摘要 常微分方程在数学中发挥着举足轻重的作用，同时它的应用在日常生活里随处可见，因此掌握常微分方程的初等解法与求解技巧是非常必要的.本论文主要论述了其发展、初等解法与求解技巧，前者主要有变量分离、积分因子、一阶隐式微分方程的参数表示，通过举例从中总结出其求解技巧，目的是掌握其求解技巧. 【关键词】变量分离一阶隐式微分方程积分因子求解技巧
Elementary Solution and Solving Skills of Ordinary Differential Equation Abstract Ordinary differential equations take up significant position in mathematics, and at the same time, the application of it can be seen everywhere in our daily life, therefore, it’s necessary to grasp the elementary solution of ordinary differential equations and solving skills. This paper mainly introduced the definition of ordinary differential equations, elementary solution method and solving skills, the former mainly included the separation of variables, integral factor, a parameter-order differential equations implicit representation, by way of examples to sum up their solving skills, the purpose is to master the skills to solve. 【Key Words】the separation of variables the first order implicit differential equation integrating factor solution techniques
页脚内容1 第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程： ①、形如 )()(y g x f dx dy = 当0)(≠y g 时，得到 dx x f y g dy )()(=，两边积分即可得到结果； 当0)(0=ηg 时，则0)(η=x y 也是方程的解。 例1.1、xy dx dy = 解：当0≠y 时，有xdx y dy =，两边积分得到)(2ln 2为常数C C x y += 所以)(11212 C x e C C e C y ±==为非零常数且 0=y 显然是原方程的解； 综上所述，原方程的解为)(1212 为常数C e C y x = ②、形如0)()()()(=+dy y Q x P dx y N x M 当0)()(≠y N x P 时，可有dy y N y Q dx x P x M ) ()()()(=，两边积分可得结果； 当0)(0=y N 时，0y y =为原方程的解，当0(0=） x P 时，0x x =为原方程的解。 例1.2、0)1()1(22=-+-dy x y dx y x
页脚内容2 解：当0)1)(1(22≠--y x 时，有dx x x dy y y 1 122-=-两边积分得到 )0(ln 1ln 1ln 22≠=-+-C C y x ，所以有)0()1)(1(22≠=--C C y x ； 当0)1)(1(22=--y x 时，也是原方程的解； 综上所述，原方程的解为)()1)(1(22为常数C C y x =--。 ⑵可化为变量可分离方程的方程： ①、形如)(x y g dx dy = 解法：令x y u = ，则udx xdu dy +=，代入得到)(u g u dx du x =+为变量可分离方程，得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x x y f =。 ②、形如)0(),(≠+=ab by ax G dx dy 解法：令by ax u +=，则b du adx dy +=，代入得到)(1u G b a dx du b =+为变量可分离方程，得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x by ax f =+。 ③、形如 )(222111c y b x a c y b x a f dx dy ++++= 解法：01、02211 =b a b a ，转化为)(by ax G dx dy +=，下同①； 02、0221 1 ≠b a b a ，???=++=++00222111c y b x a c y b x a 的解为),(00y x ，令???-=-=00y y v x x u
第十二章 常微分方程 (A) 一、是非题 1．任意微分方程都有通解。( ) 2．微分方程的通解中包含了它所有的解。( ) 3．函数x x y cos 4sin 3-=是微分方程0=+''y y 的解。( ) 4．函数x e x y ?=2是微分方程02=+'-''y y y 的解。( ) 5．微分方程0ln =-'x y x 的通解是()C x y += 2 ln 2 1 (C 为任意常数)。( ) 6．y y sin ='是一阶线性微分方程。( ) 7．xy y x y +='33不是一阶线性微分方程。( ) 8．052=+'-''y y y 的特征方程为0522=+-r r 。( ) 9． 2 2 1xy y x dx dy +++=是可分离变量的微分方程。( ) 二、填空题 1．在横线上填上方程的名称 ①()0ln 3=-?-xdy xdx y 是 。 ②()()022=-++dy y x y dx x xy 是 。 ③x y y dx dy x ln ?=是 。 ④x x y y x sin 2+='是 。 ⑤02=-'+''y y y 是 。 2．x x y x y cos sin =-'+'''的通解中应含 个独立常数。 3．x e y 2-=''的通解是 。 4．x x y cos 2sin -=''的通解是 。 5．124322+=+'+'''x y x y x y x 是 阶微分方程。 6．微分方程()06 ='-''?y y y 是 阶微分方程。
7．x y 1 =所满足的微分方程是 。 8．x y y 2='的通解为 。 9． 0=+ x dy y dx 的通解为 。 10． ()25 11 2+=+- x x y dx dy ，其对应的齐次方程的通解为 。 11．方程()012=+-'y x y x 的通解为 。 12．3阶微分方程3x y ='''的通解为 。 三、选择题 1．微分方程()043='-'+''y y y x y xy 的阶数是( )。 A ．3 B ．4 C ．5 D ． 2 2．微分方程152=-''-'''x y x y 的通解中应含的独立常数的个数为( )。 A ．3 B ．5 C ．4 D ． 2 3．下列函数中，哪个是微分方程02=-xdx dy 的解( )。 A ．x y 2= B ．2x y = C ．x y 2-= D ． x y -= 4．微分方程32 3y y ='的一个特解是( )。 A ．13+=x y B ．()3 2+=x y C ．()2 C x y += D ． ()3 1x C y += 5．函数x y cos =是下列哪个微分方程的解( )。 A ．0=+'y y B ．02=+'y y C ．0=+y y n D ． x y y cos =+'' 6．x x e C e C y -+=21是方程0=-''y y 的( )，其中1C ，2C 为任意常数。 A ．通解 B ．特解 C ．是方程所有的解 D ． 上述都不对 7．y y ='满足2|0==x y 的特解是( )。 A ．1+=x e y B ．x e y 2= C ．22x e y ?= D ． x e y ?=3 8．微分方程x y y sin =+''的一个特解具有形式( )。 A ．x a y sin *= B ．x a y cos *?=
常微分方程解题方法总结 来源：文都教育 复习过半，课本上的知识点相信大部分考生已经学习过一遍 . 接下来，如何将零散的知识点有机地结合起来，而不容易遗忘是大多数考生面临的问题 . 为了加强记忆，使知识自成体系，建议将知识点进行分类系统总结 . 著名数学家华罗庚的读书方法值得借鉴，他强调读 书要 “由薄到厚、由厚到薄 ”，对同学们的复习尤为重要 . 以常微分方程为例， 本部分内容涉及可分离变量、 一阶齐次、 一阶非齐次、 全微分方程、 高阶线性微分方程等内容， 在看完这部分内容会发现要掌握的解题方法太多， 遇到具体的题 目不知该如何下手， 这种情况往往是因为没有很好地总结和归纳解题方法 . 下面以表格的形 式将常微分方程中的解题方法加以总结，一目了然，便于记忆和查询 . 常微分方程 通解公式或解法 ( 名称、形式 ) 当 g( y) 0 时，得到 dy f (x)dx ， g( y) 可分离变量的方程 dy f ( x) g( y) 两边积分即可得到结果； dx 当 g( 0 ) 0 时，则 y( x) 0 也是方程的 解 . 解法：令 u y xdu udx ，代入 ，则 dy 齐次微分方程 dy g( y ) x dx x u g (u) 化为可分离变量方程 得到 x du dx 一 阶 线 性 微 分 方 程 dy P ( x)dx P ( x) dx Q(x) y ( e Q( x)dx C )e P( x) y dx
伯努利方程解法：令 dy P( x) y Q( x) y n（n≠0,1） 代入得到dx —u y1 n，有 du(1 n) y n dy ， du(1 n) P(x)u(1 n)Q(x) dx 求解特征方程： 2pq 0三种情况： 二阶常系数齐次线性微分方程 y p x y q x y0 二阶常系数非齐次线性微分方程y p x y q x y f ( x) (1)两个不等实根： 1 ,2 通解： y c1 e 1x c2 e 2x (2)两个相等实根：12 通解： y c1c2 x e x (3)一对共轭复根：i , 通解： y e x c1 cos x c2 sin x 通解为y p x y q x y 0 的通解与 y p x y q x y f ( x) 的特解之和. 常见的 f (x) 有两种情况： x （ 1）f ( x)e P m ( x) 若不是特征方程的根，令特解y Q m ( x)e x；若是特征方程的单根，令特 解 y xQ m ( x)e x；若是特征方程的重根， 令特解 y*x2Q m (x)e x； （2）f (x) e x[ P m ( x) cos x p n ( x)sin x] 当i不是特征值时，令 欢迎下载2
第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程： ①、形如 )()(y g x f dx dy = 当0)(≠y g 时，得到 dx x f y g dy )() (=，两边积分即可得到结果； 当0)(0=ηg 时，则0)(η=x y 也是方程的解。 例1.1、 xy dx dy = 解：当0≠y 时，有 xdx y dy =，两边积分得到)(2ln 2为常数C C x y += 所以)(112 12 C x e C C e C y ±==为非零常数且 0=y 显然是原方程的解； 综上所述，原方程的解为)(12 12 为常数C e C y x = ②、形如0)()()()(=+dy y Q x P dx y N x M 当0)()(≠y N x P 时，可有 dy y N y Q dx x P x M ) () ()()(=，两边积分可得结果； 当0)(0=y N 时，0y y =为原方程的解，当0(0=） x P 时，0x x =为原方程的解。 例1.2、0)1()1(2 2 =-+-dy x y dx y x 解：当0)1)(1(2 2 ≠--y x 时，有 dx x x dy y y 1 122-=-两边积分得到 )0(ln 1ln 1ln 22≠=-+-C C y x ，所以有)0()1)(1(22≠=--C C y x ； 当0)1)(1(2 2 =--y x 时，也是原方程的解； 综上所述，原方程的解为)()1)(1(2 2 为常数C C y x =--。 ⑵可化为变量可分离方程的方程： ①、形如 )(x y g dx dy = 解法：令x y u =，则udx xdu dy +=，代入得到)(u g u dx du x =+为变量可分离方程，得到
常微分方程计算题 2.指出下列方程中的阶数，是线性方程还是非线性方程，并说明理由； （1） t 2 2 2dt u d +t dt du +( t 2 -1)u=0 （2） dx dy =x 2+y 2 ; （3）dx dy + 2 x y =0 3.求曲线族y=C 1e x +C 2x e x 所满足的微分方程 4．验证函数y= C 1e x 2+ C 2e x 2-是微分方程y `` -4y=0的解，进一步验证它是通解。 5.试用一阶微分方程形式不变性求解方程dx dy =2x 6.什么叫积分一个微分方程 7．什么是求解常微分方程的初等积分法 8．分离变量一阶方程的特征是什么 9．求下列方程的通解 （1） y ` ＝sinx （2） x 2 y 2 y ` +1=y （3） tgx dx dy =1+y （4） dx dy =exp(2x-y) （5） dx dy =21y 2- （6） x 2 ydx=(1- y 2 +x-2 x 2 y 2 )dx （7）( x 2 +1)( y 2 -1)dx+xydy=0 10.叙述齐次函数的定义 11．试给出一阶方程y ` ＝f(x,y)或p(x,y)dx+ q(x,y)dy=0为齐次方程的特征。说明二
个方程的关系。 12．求解齐次方程通常用什么初等变换，新旧函数导数关系如何 13．求解下列方程 dx dy ＝2 22y x xy - 14．求解下列方程 （1）（x+2y ）dx —xdy=0 (2) dx dy =x y +y x 2 15. dx dy =22y x xy + 16(x 2 +y 2 )dx —2xydy=0 17. dx dy =5 242+---y x x y 18―――――19 20―――――――27
第八章 常微分方程 【教学要求】 一、了解微分方程的基本概念：微分方程，微分方程的阶、解、特解、通解、初始条件和初值问题，线性微分方程。 二、熟练掌握一阶可分离变量微分方程的解法。 三、熟练掌握一阶线性非齐次微分方程)()(x q y x p y =+' 的解法——常数变易法和公式法。 四、理解线性微分方程解的性质和解的结构。 五、熟练掌握二阶线性常系数齐次微分方程0=+'+''qy y p y 的解法——特征根法。 会根据特征根的三种情况，熟练地写出方程的通解，并根据定解的条件写出方程特解。 六、熟练掌握二阶线性常系数非齐次微分方程qy y p y +'+'' )(x f =，当自由项f (x )为某些特殊情况时的解法——待定系数法。 所谓f (x )为某些特殊情况是指f (x )为多项式函数，指数函数 或它们的和或乘积形式、三角函数x x x ββαsin cos ,e 。 关键是依据f (x )的形式及特征根的情况，设出特解y *,代入原方程，定出y *的系数。 【教学重点】 一阶可分离变量微分方程、一阶线性微分方程、二阶线性常系数微分方程的解法。 【典型例题】 。的阶数是微分方程例)(e )(12x y y y =-'+'' 2.1.B A 4. 3.D C 解：B 。的特解形式是微分方程例)( e 232x x y y y +=+'-'' x x x b ax B b ax A e )(.e ).(++ x x c b ax D cx b ax C e ).(e ).(++++ 解：C 是一阶线性微分方程。下列方程中例)( ,3 x x y y x B y A y x cos sin 1.e .2=+'='+ y x y D y y x y C ='=+'+''.0 . 解：B ???=='++1)1(0)1(4y y x y y 求解初值问题例 ??-=+x x y y y d )1(d 解：由变量可分离法得 c x y y ln ln 1ln +-=+∴ 代入上式得通解为由21ln ln 1)1(=?=c y x y y 211=+ 的特解。满足求解微分方程例1)0(e 252==-'y x y y x 解：由公式法得 ]d e e 2[e d 12d 1c x x y x x x +???=---?
应 用 题（每题10分） 1、设()f x 在(,)-∞∞上有定义且不恒为零，又()f x '存在并对任意,x y 恒有 ()()()f x y f x f y +=，求()f x 。 2、设()()()F x f x g x =，其中函数(),()f x g x 在(,)-∞∞内满足以下条件 ()(),()(),(0)0,()()2x f x g x g x f x f f x g x e ''===+= （1）求()F x 所满足的一阶微分方程； （2）求出()F x 的表达式。 3、已知连续函数()f x 满足条件320 ()3x x t f x f dt e ??=+ ??? ?，求()f x 。 4、已知函数()f x 在(0,)+∞内可导，()0,lim ()1x f x f x →+∞ >=，且满足 1 1 0()lim ()h x h f x hx e f x →? ?+ ?= ? ?? ? ，求()f x 。 5、设函数()f x 在(0,)+∞内连续，5 (1)2 f =，且对所有,(0,)x t ∈+∞，满足条件 1 1 1 ()()()xt x t f u du t f u du x f u du =+? ??，求()f x 。 6、求连续函数()f x ，使它满足10 ()()sin f tx dt f x x x =+?? 。 7、已知可微函数()f t 满足 31() ()1()x f t dt f x t f t t =-+?，试求()f x 。 8、设有微分方程 '2()y y x ?-=， 其中21 ()01x x x ?? 。试求在(,)-∞∞内的连续函 数()y y x =使之在(,1)-∞和()1,+∞内部满足所给方程，且满足条件(0)0y =。 9、设位于第一象限的曲线()y f x = 过点122?? ? ? ?? ，其上任一点(,)P x y 处的法线与y 轴的交点为Q ，且线段PQ 被x 轴平分。 （1）求曲线()y f x =的方程； （2）已知曲线sin y x =在[0,]π上的弧长为l ，试用l 表示曲线()y f x =的弧长s 。 10、求微分方程(2)0xdy x y dx +-=的一个解()y y x =，使得由曲线()y y x =与直线 1,2x x ==以及x 轴所围成的平面图形绕x 轴旋转一周的旋转体体积最小。 11、设曲线L 位于xOy 平面的第一象限内，L 上任一点M 处的切线与y 轴总相交，交点记为
第8章 常微分方程数值解法 本章主要内容： 1．欧拉法、改进欧拉法. 2．龙格-库塔法。 3．单步法的收敛性与稳定性。 重点、难点 一、微分方程的数值解法 在工程技术或自然科学中，我们会遇到的许多微分方程的问题，而我们只能对其中具有较简单形式的微分方程才能够求出它们的精确解。对于大量的微分方程问题我们需要考虑求它们的满足一定精度要求的近似解的方法，称为微分方程的数值解法。本章我们主要 讨论常微分方程初值问题?????==00 )() ,(y x y y x f dx dy 的数值解法。 数值解法的基本思想是：在常微分方程初值问题解的存在区间[a,b]内，取n+1个节点a=x 0＜x 1＜…＜x N =b （其中差h n = x n –x n-1称为步长，一般取h 为常数，即等步长），在这些节点上把常微分方程的初值问题离散化为差分方程的相应问题，再求出这些点的上的差分方程值作为相应的微分方程的近似值（满足精度要求）。 二、欧拉法与改进欧拉法 欧拉法与改进欧拉法是用数值积分方法对微分方程进行离散化的一种方法。 将常微分方程),(y x f y ='变为() *+=?++1 1))(,()()(n x n x n n dt t y t f x y x y 1．欧拉法（欧拉折线法） 欧拉法是求解常微分方程初值问题的一种最简单的数值解法。 欧拉法的基本思想：用左矩阵公式计算（＊）式右端积分，则得欧拉法的计算公式为：N a b h N n y x hf y y n n n n -= -=+=+)1,...,1,0(),(1 欧拉法局部截断误差 11121 )(2 ++++≤≤''=n n n n n x x y h R ξξ或简记为O （h 2）。
《常微分方程》模拟练习题及参考答案 一、填空题（每个空格4分，共80分） 1、n 阶线性齐次微分方程基本解组中解的个数恰好是 n 个。 2、一阶微分方程 2=dy x dx 的通解为 2=+y x C （C 为任意常数） ，方程与通过点（2，3）的特解为 2 1=-y x ，与直线y=2x+3相切的解是 2 4=+y x ，满足条件3 3ydx =?的解为 22=-y x 。 3、李普希兹条件是保证一阶微分方程初值问题解惟一的 必要 条件。 4、对方程 2()dy x y dx =+作变换 =+u x y ，可将其化为变量可分离方程，其通解为 tan()=+-y x C x 。 5、方程 21d d y x y -=过点)1,2 (π 共有 无数 个解。 6、方程 ''2 1=-y x 的通解为 42 12122=-++x x y C x C ，满足初始条件13|2,|5====x x y y 的特解为 4219 12264 =-++x x y x 。 7、方程 x x y x y +-=d d 无 奇解。 8、微分方程2260--=d y dy y dx dx 可化为一阶线性微分方程组 6?=??? ?=+??dy z dx dz z y dx 。 9、方程 y x y =d d 的奇解是 y=0 。 10、35323+=d y dy x dx dx 是 3 阶常微分方程。 11、方程 22dy x y dx =+满足解得存在唯一性定理条件的区域是 xoy 平面 。 12、微分方程22450d y dy y dx dx --=通解为 512-=+x x y C e C e ，该方程可化为一阶线性微分方程组 45?=??? ?=+??dy z dx dz z y dx 。 13、二阶线性齐次微分方程的两个解12(),()y x y x ??==成为其基本解组的充要条件是 线性无关 。
利用MATLAB求解常微分方程数值解
目录 1. 内容简介 (1) 2. Euler Method(欧拉法)求解 (1) 2.1. 显式Euler法和隐式Euler法 (2) 2.2. 梯形公式和改进Euler法 (3) 2.3. Euler法实用性 (4) 3. Runge-Kutta Method(龙格库塔法)求解 (5) 3.1. Runge-Kutta基本原理 (5) 3.2. MATLAB中使用Runge-Kutta法的函数 (7) 4. 使用MATLAB求解常微分方程 (7) 4.1. 使用ode45函数求解非刚性常微分方程 (8) 4.2. 刚性常微分方程 (9) 5. 总结 (9) 参考文献 (11) 附录 (12) 1. 显式Euler法数值求解 (12) 2. 改进Euler法数值求解 (12) 3. 四阶四级Runge-Kutta法数值求解 (13) 4.使用ode45求解 (14)
1.内容简介 把《高等工程数学》看了一遍，增加对数学内容的了解，对其中数值解法比较感兴趣，这大概是因为在其它各方面的学习和研究中经常会遇到数值解法的问题。理解模型然后列出微分方程，却对着方程无从下手，无法得出精确结果实在是让人难受的一件事情。 实际问题中更多遇到的是利用数值法求解偏微分方程问题，但考虑到先从常微分方程下手更为简单有效率，所以本文只研究常微分方程的数值解法。把一个工程实际问题弄出精确结果远比弄清楚各种细枝末节更有意思，因此文章中不追求非常严格地证明，而是偏向如何利用工具实际求解出常微分方程的数值解，力求将课程上所学的知识真正地运用到实际方程的求解中去，在以后遇到微分方程的时候能够熟练运用MATLAB得到能够在工程上运用的结果。 文中求解过程中用到MATLAB进行数值求解，主要目的是弄清楚各个函数本质上是如何对常微分方程进行求解的，对各种方法进行MATLAB编程求解，并将求得的数值解与精确解对比，其中源程序在附录中。最后考察MATLAB中各个函数的适用范围，当遇到实际工程问题时能够正确地得到问题的数值解。 2.Euler Method(欧拉法)求解 Euler法求解常微分方程主要包括3种形式，即显式Euler法、隐式Euler法、梯形公式法，本节内容分别介绍这3种方法的具体内容，并在最后对3种方法精度进行对比，讨论Euler法的实用性。 本节考虑实际初值问题 使用解析法，对方程两边同乘以得到下式
摘要 (2) 1.何谓奇解 (2) 2.奇解的产生 (3) 3.包络跟奇解的关系 (4) 4.理论上证明C-判别曲线与P-判别曲线方法 (5) 4.1 克莱罗微分方程 (9) 5.奇解的基本性质 (12) 5.1 定理１ (12) 5.2 定理２ (14) 5.3 定理３ (14) 6.小结 (14) 参考文献： (15)
一阶常微分方程的奇解 摘要 在常微分方程中，我们知道方程的解可以有多种，现在我们来讨论求奇解的方法。我们看到某些微分方程，会存在一些特殊的积分曲线，他并不属于这方程的积分曲线族，但是，在这些特殊的积分曲线上的每一点处，都有积分曲线族中的一条曲线和他在此处相切。从而我们引出了积分曲线族的包络，而为了求微分方程的奇解，，我们应先求出他的通解，然后求通解的包络。 关键词：奇解，包络，C-判别式，P-判别式 1.何谓奇解 设一阶隐式方程),,(,y y x F =0有一特解
)(:x y ψ=Γ，j x ∈ 如果对每一点Γ∈P ，在P 点的任何一个领域，方程),,(,y y x F =0都有一个不同于Γ的解在P 点与Γ相切，则称Γ是微分方程的),,(,y y x F =0的奇解 定义：如果一个一阶微分方程的一个特解的积分曲线上的每一点都至少和这个微分方程的不同的积分曲线相切，并且这相切的积分曲线在切点的任何邻域都不重合，则称这个特解为这个微分方程的奇解 2.奇解的产生 先看一个例子，求方程 033=-?? ? ??y dx dy （1） 或与它等价的方程 3y dx dy = 的解。 经分离变量后，可得（1）的通解 3)(27 1c x y += 容易看出，y=0也是原方程的一个解。现在来研究这个解y=0有什么特殊的地方。由图我们看到，在解y=0上的每一 点)0,(0x 处相切，这种特殊的积分曲线y=0 称为奇积分曲线，他所对应的解就是奇 解，这就是奇解的产生。 我们现在给出曲线族包络的定义 某些微分方程，存在一些特殊的积分 曲线，会存在一些特殊的积分曲线，他并 不属于这方程的积分曲线族，但是，在这些特殊的积分曲线上的每一点处，都有积分曲线族中的一条曲线和他在此处相切。在几何学里，这些特殊的积分曲线称为上述积分曲线族的包络，在微分方程里，这些特殊的积分曲线所对应的解酒称为方程的奇解。
常微分方程练习题及答案(复习题)
常微分方程练习试卷 一、 填空题。 1. 方程23 2 10d x x dt +=是 阶 （线性、非线性）微分方程. 2. 方程 ()x dy f xy y dx =经变换_______，可以化为变量分离方程 . 3. 微分方程 3230d y y x dx --=满足条件(0)1,(0)2y y '==的解有 个. 4. 设常系数方程 x y y y e αβγ'''++=的一个特解*2()x x x y x e e xe =++，则此方程的系数α= ，β= ，γ= . 5. 朗斯基行列式 ()0W t ≡是函数组12(),(),,()n x t x t x t L 在a x b ≤≤上线性相关的 条件. 6. 方程 22(2320)0xydx x y dy ++-=的只与y 有关的积分因子为 . 7. 已知 ()X A t X '=的基解矩阵为()t Φ的，则()A t = . 8. 方程组 20'05??=???? x x 的基解矩阵为 ． 9.可用变换 将伯努利方程 化为线性方程. 10 .是满足方程 251y y y y ''''''+++= 和初始条件 的唯一解. 11.方程 的待定特解可取 的形式: 12. 三阶常系数齐线性方程 20y y y '''''-+=的特征根是 二、 计算题 1.求平面上过原点的曲线方程, 该曲线上任一点处的切线与切点和点(1,0)的连线相互垂直. 2．求解方程13 dy x y dx x y +-=-+. 3. 求解方程 222()0d x dx x dt dt += 。 4．用比较系数法解方程. . 5．求方程 sin y y x '=+的通解. 6．验证微分方程 22(cos sin )(1)0x x xy dx y x dy -+-=是恰当方程，并求出它的通解.