常微分方程数值解实验报告

《常微分方程数值解》

实验报告

姓名：李想

班级：数学091

学号：0909300114

运用Euler 方法及其改进方法求解常微分方程初值问题，具体问题在教材课后习题中学生任意选取。

实验内容：

dy/dx=(5*(y^2)+2*(x^2))/(x*y)

实验原理：

?????

?

??????

?????++=++=+

+==++++==+)

,()2,2()2,2()

,()22(6

)0(342312143211

0hk y h x f k hk y h x f k hk y h

x f k y x f k k k k k h y y y y n n n n n n n

n n n ?????

?

???

???

?????+-++=+-+=+

+==++++==+)

,()3,32()3,3()

,()33(8

)

0(321421312143211

0hk hk hk y h x f k hk hk y h x f k hk y h

x f k y x f k k k k k h y y y y n n n n n n n

n n n

图1

图1

表1

实验体会：

学习并了解欧拉方法和欧拉改进方法，初步练习MATLAB软件的方法，使我深深的感受到了它在线性代数的矩阵，向量等各方面领域的重要性！

常微分方程练习题及答案复习题)

常微分方程练习试卷 一、 填空题。 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 第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程： ①、形如 )()(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

常微分方程期末试题B答案

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 +-= （*） ，两端同时关于求导，

一阶常微分方程的奇解

摘要.................................................... 错误!未定义书签。 1.何谓奇解.............................................. 错误!未定义书签。 2.奇解的产生............................................ 错误!未定义书签。 3.包络跟奇解的关系...................................... 错误!未定义书签。 4.理论上证明C-判别曲线与P-判别曲线方法................. 错误!未定义书签。 克莱罗微分方程 ..................................... 错误!未定义书签。 5.奇解的基本性质........................................ 错误!未定义书签。 定理１ ............................................. 错误!未定义书签。 定理２ ............................................. 错误!未定义书签。 定理３ ............................................. 错误!未定义书签。 6.小结.................................................. 错误!未定义书签。参考文献：.............................................. 错误!未定义书签。

常微分方程初值问题的数值解法

第七章 常微分方程初值问题的数值解法 --------学习小结 一、本章学习体会 通过本章的学习，我了解了常微分方程初值问题的计算方法，对于解决那些很难求解出解析表达式的，甚至有解析表达式但是解不出具体的值的常微分方程非常有用。在这一章里求解常微分方程的基本思想是将初值问题进行离散化，然后进行迭代求解。在这里将初值问题离散化的方法有三种，分别是差商代替导数的方法、Taylor 级数法和数值积分法。常微分方程初值问题的数值解法的分类有显示方法和隐式方法，或者可以分为单步法和多步法。在这里单步法是指计算第n+1个y 的值时,只用到前一步的值，而多步法则是指计算第n+1个y 的值时，用到了前几步的值。通过对本章的学习，已经能熟练掌握如何用Taylor 级数法去求解单步法中各方法的公式和截断误差，但是对线性多步法的求解理解不怎么透切，特别是计算过程较复杂的推理。 在本章的学习过程中还遇到不少问题，比如本章知识点多，公式多，在做题时容易混淆，其次对几种R-K 公式的理解不够透彻，处理一个实际问题时，不知道选取哪一种公式，通过课本里面几种方法的计算比较得知其误差并不一样，，这个还需要自己在往后的实际应用中多多实践留意并总结。 二、本章知识梳理 常微分方程初值问题的数值解法一般概念 步长h ，取节点0,(0,1,...,)n t t nh n M =+=，且M t T ≤，则初值问题000 '(,),()y f t y t t T y t y =≤≤?? =?的数值解法的一般形式是 1(,,,...,,)0,(0,1,...,)n n n n k F t y y y h n M k ++==-

常微分方程数值解实验报告

常微分方程数值解实验报告 学院:数学与信息科学 专业:信息与计算科学 姓名：郑思义 学号:2０121６52４ 课程：常微分方程数值解 实验一：常微分方程得数值解法 1、分别用Eｕleｒ法、改进得Euler法(预报校正格式)与S—K法求解初值问题。(h＝0、1)并与真解作比较。 1、1实验代码: %欧拉法 functioｎ [ｘ，y]=naeuler(dｙｆｕｎ，xspan,y０,ｈ) %ｄyfun就是常微分方程,xｓpan就是x得取值范围，ｙ０就是初值,h就是步长 x＝ｘsｐan(1）:h:xspan(2)； y(1）=y０； ｆｏrｎ=1:length(x)-１ ｙ（n+１）=ｙ(n)+h*feval（dｙｆun,ｘ(ｎ),ｙ(ｎ））; eｎd ％改进得欧拉法 ｆunction [x,ｍ,ｙ]=naｅuler2(ｄyfun,xspaｎ,y０,h) %ｄyfｕn就是常微分方程,xｓpan就是ｘ得取值范围,y0就是初值，h就是步长。 %返回值ｘ为x取值,ｍ为预报解,ｙ为校正解 x＝xspan(1)：h:xspaｎ(2); y（1)=y0; m=zerｏs（lenｇth(x)－1，１); foｒｎ=1:leｎgｔh(x）-1 k1＝fｅｖal（dｙfun,x(n),y(n）)； y(n＋1)=y（n)+h*k1; ｍ（n)=y（n+1); k2＝ｆevａｌ(dyfuｎ,x(n＋1）,y(n+1）); y（ｎ+1）=y(n)+h*（k１+k2)／2;

end %四阶S—K法 functioｎ [ｘ，ｙ]=rk(ｄyfun,xspａn,y0,h） %dｙｆun就是常微分方程,xｓpaｎ就是ｘ得取值范围,y0就是初值，h就是步长。 ｘ＝xsｐan(1):h：xsｐan(2)； ｙ(１)=y0; for n=1:lｅｎgth(x)-１ ｋ1=fｅｖａl(dyfｕn，x(n），y(ｎ）)； k2=fｅｖal(dyfun,x(ｎ)+ｈ/2，y(ｎ）+(h＊ｋ1)/2); k3=feval(ｄyfuｎ,ｘ(n)+ｈ/２，y(n)+(h*ｋ2)/2); k4=ｆevａl(dｙfun，x(ｎ)+h,y(n)+ｈ＊k3）； y(n+１)＝ｙ(n)+（h/6）*(k1+２*k2+2*k3+k４）; enｄ ％主程序 x=［０:0、1:１］;y＝ｅxp(-ｘ)+x； dyfun=inline('-ｙ+x+1'）； [x1，y1]=naeuler(ｄyｆｕn,[0,１]，1,0、1); ［x2,m，y2]＝naeuler2（dｙfuｎ,[０，1］,1,０、1）; [x3,y3]=rｋ(dyｆuｎ,[0,１］,1,０、１); ploｔ(x，y,'r',x1,y１,'+'，x2,y2,＇*',x3，y３,'o'); xlabｅl('ｘ');yｌabel('y＇); legｅnd(＇y为真解',＇y1为欧拉解＇,'ｙ2为改进欧拉解','y3为S—K解','Locatiｏn','ＮorthWeｓｔ＇)； 1、２实验结果： x 真解y 欧拉解y1 预报值ｍ校正值y2 S—K解y3 0、０１、0000 1、0000 1、０000 １、000０ 0、1 1、0048 1、000０1、００00 １、005０１、0０4８ 0、2 １、0187 １、0100 1、０145 1、０190 1、0187 ０、3 1、0408 １、02９0 1、０３7１1、0412 1、0４0８ ０、4 １、0７03 1、0５６1 1、06７1 1、0708 １、0７03 0、5 1、10６5 1、09０５1、1０3７1、107１1、１0６5 0、6 1、１488 1、１３1４１、1464 １、149４1、１488 0、7 1、1966 １、17８3 1、194５1、19７2 1、１966 0、８1、2493 １、2３05 1、2475 1、2500 1、２４９3 0、9 1、306６1、２87４1、３050 1、3072 １、３066 1、0 1、3679 1、34８７1、36６5 1、３６85 1、36７9

(整理)常微分方程试题及参考答案

常微分方程试题 一、填空题(每小题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分)解方程

总结一阶常微分方程奇解的求法

总结一阶微分方程奇解的求法 摘要：利用有关奇解的存在定理，总结出求一阶微分方程奇解的几种方法，并通过一些具体的例题说明这几种方法的应用 Using relevant theorems to develop several methods of finding singular solution of ordinary differential equation. In addition, illustrate the application of these methods through the concrete examples. 关键词：常微分方程 奇解 c-判别式 p-判别式 方法一：利用c-判别式求奇解 设一阶微分方程0, ,=?? ? ?? dx dy y x F ① 可求出方程①的通解为()0,,=c y x φ ② 如果()()???==0 ,,0,,' c y x c y x c φφ ③ 是微分方程①的解，且对③式满足：()()02 '2 '≠+y x φφ ④ 则③是微分方程①的奇解，且是通解②的包络。 例1：方程() 2 2 2 x x y dy dx dy dx + -= 的奇解 解：首先，本具题意求出该微分方程的通解为2 2 2 c cx y x ++= 与4 2 x y = 其中c 为任意常数 当时2 2 2 c cx y x ++= ， ()y c cx x c y x -++= 2 2 2 ,,φ 其相应的c -判别式为 ? ??=+=-++02022x 2 c x y c cx 易得到： ? ??=-=2 2c y c x

代入原微分方程，可知? ??=-=2 2c y c x 不是原微分方程的解； 当4 2 x y = 时，易求出2 ,1''x y x ==φφ，则有()()02 '2 '≠+y x φφ 故4 2 x y = 为原微分方程的奇解 例2：试求微分方程() () y y dy dx 9 42 2 1= -的奇解 解：首先，根据题意求出微分方程的通解为：()()0322=---y y c x 其中c 为任意常数 再由相应的c-判别式： ()()()? ??=--=---020 322c x y y c x 易求出：? ??==0y c x 或 ???==3y c x 当???==0y c x 时，代入原微分方程成立； 所以? ??==0y c x 为原微分方程的解 且有()02'=--=c x x φ；()()93232 '-=---=y y y y φ 满足(Φ‘ x )2 +(Φ‘ y )2≠0 易验证???==3y c x 不是原微分方程的解 故x=c, y=0 是元微分方程的奇解。 方法二：利用p-判别法求奇解 在微分方程①中，设y ′=p,则此方程的p-判别式为： ()()?????==0,,0 ,,' p y x F p y x F p ⑤ 消去p 之后得到的函数y=?(x)是微分方程①身为解，

常微分方程习题及答案.[1]

第十二章 常微分方程 (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 *?=

常微分方程数值解

第四章常微分方程数值解 [课时安排]6学时 [教学课型]理论课 [教学目的和要求] 了解常微分方程初值问题数值解法的一些基本概念，如单步法和多步法，显式和隐式，方法的阶数，整体截断误差和局部截断误差的区别和关系等；掌握一阶常微分方程初值问题的一些常用的数值计算方法，例如欧拉（Euler）方法、改进的欧拉方法、龙贝－库塔（Runge-Kutta）方法、阿达姆斯（Adams）方法等，要注意各方法的特点及有关的理论分析；掌握构造常微分方程数值解的数值积分的构造方法和泰勒展开的构造方法的基本思想，并能具体应用它们导出一些常用的数值计算公式及评估截断误差；熟练掌握龙格－库塔（Ｒ－Ｋ）方法的基本思想，公式的推导，Ｒ－Ｋ公式中系数的确定，特别是能应用“标准四阶Ｒ－Ｋ公式”解题；掌握数值方法的收敛性和稳定性的概念，并能确定给定方法的绝对稳定性区域。[教学重点与难点] 重点：欧拉方法，改进的欧拉方法，龙贝－库塔方法。 难点：R—K方法，预估-校正公式。 [教学内容与过程] 4.1 引言 本章讨论常微分方程初值问题 (4.1.1) 的数值解法，这也是科学与工程计算经常遇到的问题，由于只有很特殊的方程能用解析方法求解，而用计算机求解常微分方程的初值问题都要采用数值方法.通常我们假定(4.1.1)中 f(x,y)对y满足Lipschitz条件，即存在常数L＞0，使对，有 (4.1.2) 则初值问题(4.1.1)的解存在唯一. 假定(4.1.1)的精确解为，求它的数值解就是要在区间上的一组离散点 上求的近似.通常取 ，h称为步长，求(4.1.1)的数值解是按节点的顺序逐步 推进求得.首先，要对方程做离散逼近，求出数值解的公式，再研究公式的局部截

[vip专享]2013春数学实验基础 实验报告(1)常微分方程

1. 分别用Euler 法和ode45解下列常微分方程并与解析解比较： (1) 30,1)0(,<<=+='x y y x y 编写Euler 法的matlab 函数，命名为euler.m function [t,y]=euler(odefun,tspan,y0,h)t=tspan(1):h:tspan(2);y(1)=y0; for i=1:length(t)-1 y(i+1)=y(i)+h*feval(odefun,t(i),y(i));end t=t';y=y'; 下面比较三者的差别：% ode45 odefun=inline('x+y','x','y');[x1,y1]=ode45(odefun,[0,3],1);plot(x1,y1,'ko');pause hold on ;% Euler·¨ [x2,y2]=euler(odefun,[0,3],1,0.05);plot(x2,y2,'r+');pause hold on ; % 解析解 y0=dsolve('Dy=t+y','y(0)=1');ezplot(y0,[0,3]);pause hold off ; legend('ode45','euler 法','解析解');

Euler 法只有一阶精度，所以实际应用效率比较差，而ode45的效果比较好，很接近真实值。 (2) 2 0.01()2sin(),(0)0,(0)1,05 y y y t y y t ''''-+===<<先写M 文件ex1_2fun.m function f=ex1_2fun(t,y)f(1)=y(2); f(2)=0.01*y(2).^2-2*y(1)+sin(t);f=f(:);% ode45 [t1,y1]=ode45(@ex1_2fun,[0,5],[0;1]);plot(t1,y1(:,1),'ko'); % 解析解 s=dsolve('D2y-0.01*(Dy)^2+2*y=sin(t)','y(0)=0','Dy(0)=1','t') s = [ empty sym ] %由此可知该微分方程无解析解2. 求一通过原点的曲线，它在处的切线斜率等于若上限增为1.58，1.60会(,)x y 2 2,0 1.57.x y x +<

2018常微分方程考研复试真题及答案

常微分方程计算题 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

一阶常微分方程的奇解

摘要 (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 称为奇积分曲线，他所对应的解就是奇 解，这就是奇解的产生。 我们现在给出曲线族包络的定义 某些微分方程，存在一些特殊的积分 曲线，会存在一些特殊的积分曲线，他并 不属于这方程的积分曲线族，但是，在这些特殊的积分曲线上的每一点处，都有积分曲线族中的一条曲线和他在此处相切。在几何学里，这些特殊的积分曲线称为上述积分曲线族的包络，在微分方程里，这些特殊的积分曲线所对应的解酒称为方程的奇解。

常微分方程数值解法

第七章 常微分方程数值解法 常微分方程中只有一些典型方程能求出初等解(用初等函数表示的解)，大部分的方程是求不出初等解的。另外，有些初值问题虽然有初等解，但由于形式太复杂不便于应用。因此，有必要探讨常微分方程初值问题的数值解法。本章主要介绍一阶常微分方程初值问题的欧拉法、龙格-库塔法、阿达姆斯方法，在此基础上推出一阶微分方程组与高阶方程初值问题的 数值解法；此外，还将简要介绍求解二阶常微分方程值问题的差分方法、试射法。 第一节 欧拉法 求解常微分方程初值问题 ?????==0 0)() ,(y x y y x f dx dy （1） 的数值解，就是寻求准确解)(x y 在一系列离散节点 <<<<

(完整版)计算物理实验报告常微分方程

常微分方程的边值问题和本征值问题 一、问题描述 利用搜索法和弦割法，得到该常微分方程的本征值，再利用打靶法计算多个本征值。 二、解决方法 （一）搜索法 1.先随便猜测k的一个试验值,程序中令k=1 2.由Numerov算法 根据本题的条件，kn+1=kn=kn-1=k,s=0,得到yn+2,yn+1,yn间的迭代公式 令con=(k*h)^2/12 yn+2=2*(1-5*con)*yn+1/(1+con)-yn 3自己给定φ的初始条件，然后利用公式得到边界值φ(1) 4.然后以小的步长dk增加k值，这里令dk=1,每当φ(1)改变符号时，就将步长减半后倒退回来重复 5.当步长小于所要求的容许误差时终止程序，此时的k值即为所求。 （二）弦割法 1.随便猜测两个k值，这里令k0=1,k1=2 2.自己给定φ的初始条件，对两个k值分别利用上述公式进行迭代，得到边界值y1(1)和y2(1)。 3.比较y1(1)和y2(1)的绝对值大小。若绝对值大，说明对应的k值距离本征值距离较远。 4.将（k0+k1）/2赋给k2,边界值绝对值小的对应的k值保持不变，边界值绝对值大的对应k值重新定位k2的值。 5.重复进行实验，当y1(1)和y(2)的差的绝对值小于容许误差时终止程序。此时k1的值即为所求。 当搜索法和弦割法大致求出了一个本征值后，利用打靶法，调整k值再度进行搜索，得到多个本征值，绘出其中一个本征值对应的函数图像，观察其性质。 三、程序实现 1.搜索法 subroutine add(t,y0,y1) !利用子程序表示函数值的迭代 implicit none real(8)::t,h,con,y0,y1,y2

二阶常微分方程解

二阶常微分方程解

————————————————————————————————作者: ————————————————————————————————日期:

第七节 二阶常系数线性微分方程 的解法 在上节我们已经讨论了二阶线性微分方程解的结构，二阶线性微分方程的求解问题,关键在于如何求二阶齐次方程的通解和非齐次方程的一个特解。本节讨论二阶线性方程的一个特殊类型，即二阶常系数线性微分方程及其求解方法。先讨论二阶常系数线性齐次方程的求解方法。 §7.1 二阶常系数线性齐次方程及其求解方法 设给定一常系数二阶线性齐次方程为 ?? 22 dx y d +p dx dy ＋ｑy=０ (7.1) 其中p 、q 是常数,由上节定理二知,要求方程(7．1)的通解，只要求出其任意两个线性无关的特解y 1,ｙ2就可以了，下面讨论这样两个特解的求法。 我们先分析方程(7.1)可能具有什么形式的特解, 从方程的形式上来看，它的特点是22dx y d ，dx dy ,y 各乘以 常数因子后相加等于零，如果能找到一个函数ｙ,其

22dx y d ，dx dy ，ｙ之间只相差一个常数因子，这样的函数有可能是方程(７.1)的特解,在初等函数中,指数函数e rx ,符合上述要求，于是我们令 y=e r ｘ (其中r 为待定常数）来试解 将y ＝e rx ,dx dy =re r ｘ,22dx y d ＝r ２e r ｘ 代入方程(7.1) 得 r 2e ｒx ＋ｐre rx ＋ｑｅrx ＝０ 或 e r ｘ(r ２+pr+q ）＝0 因为e rx ≠０,故得 ? r 2 +pr ＋q=0 由此可见，若r 是二次方程 ?? r 2＋pr ＋q=０ (7．2) 的根，那么e r ｘ就是方程(7.１)的特解,于是方程（７.1)的求解问题，就转化为求代数方程(7.2）的根问题。称（7.2)式为微分方程(７.１)的特征方程。 特征方程(7.2)是一个以r 为未知函数的一元二次代数方程。特征方程的两个根r 1，r 2，称为特征根,由代数知识，特征根r １，r 2有三种可能的情况,下面我们分别进行讨论。 (1)若特证方程（7．2）有两个不相等的实根r 1， r 2，此时e r １ｘ ,e r2x 是方程（７.1)的两个特解。

常微分方程应用题和答案

应 用 题（每题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 轴总相交，交点记为

常微分方程的数值解

实验4 常微分方程的数值解 【实验目的】 1．掌握用MATLAB软件求微分方程初值问题数值解的方法； 2．通过实例用微分方程模型解决简化的实际问题； 3．了解欧拉方法和龙格-库塔方法的基本思想和计算公式，及稳定性等概念。 【实验内容】 题3 小型火箭初始重量为1400kg，其中包括1080kg燃料。火箭竖直向上发射时燃料燃烧率为18kg/s，由此产生32000N的推力，火箭引擎在燃料用尽时关闭。设火箭上升时空气阻力正比于速度的平方，比例系数为m，求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点的时的高度和加速度，并画出高度、速度、加速度随时间变化的图形。 模型及其求解 火箭在上升的过程可分为两个阶段，在全过程中假设重力加速度始终保持不变，g=s2。 在第一个过程中，火箭通过燃烧燃料产生向上的推力，同时它还受到自身重力（包括自重和该时刻剩余燃料的重量）以及与速度平方成正比的空气阻力的作用，根据牛顿第二定律，三个力的合力产生加速度，方向竖直向上。因此有如下二式： a=dv/dt=/m=/(1400-18t) dh/dt=v 又知初始时刻t=0，v=0，h=0。记x(1)=h，x(2)=v，根据MATLAB 可以求出0到60秒内火箭的速度、高度、加速度随时间的变化情况。程序如下： function [ dx ] = rocket( t,x ) a=[*x(2)^2)/(1400-18*t)]; dx=[x(2);a]; end ts=0:1:60;

x0=[0,0]; [t,x]=ode45(@rocket,ts,x0); h=x(:,1); v=x(:,2); a=[*(v.^2))./(1400-18*t)]; [t,h,v,a]; 数据如下： t h v a 000