常微分方程初值问题答案

常微分方程初值问题答案
常微分方程初值问题答案

1.(10分)对常微分方程初值问题(0)1(01)

dy

y

dx y x ?=-???=≤≤?

取步长0.1,h = 分别用改进的Euler 法和标准的四阶Runge-Kutta 法作数值计算,写出公式和简要推导过程,并把结果填入表内。

解:(1) 改进的Euler 方法: 代入公式得10.905n n y y +=,即0.905n n y = …2分 (2)标准的四阶Runge-Kutta 方法:

1

12341213

2430.1(22)0.90483756

(0.05)0.95(0.05)0.9525(0.1)0.90475n n n n n n

n n

n n y y k k k k y k y k y k y k y k y

k y k y +?=++++=??

=-??

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

即0.9048375n n y = ……(4分)

2. 对常微分方程初值问题12

(0)1(01)

dy

y

dx y x ?=-???=≤≤? 取步长0.1,h = 分别用改进的Euler 法和标准的四阶Runge-Kutta 法作数值计算,写出公式和推导过程,并把结果填入表内。

解:(1) 改进的Euler 方法: 代入公式得10.95125n n y y +=,即0.95125n n y =

……………….(2分)

(2)标准的四阶Runge-Kutta 方法:

1

12341213

2430.1(22)0.9512196

/2(0.05)/20.4875(0.05)/20.4878125(0.1)/20.47622n n n n n n

n n

n n y y k k k k y k y k y k y k y k y

k y k y +?=++++=??

=-??

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

即0.95145314n n y =……(4分)

《数值分析》复习题

一、填空题

1.绝对误差限=末位的一半+单位,相对误差限=绝对误差限/原值*100%

1. 度量一根杆子长250厘米,则其绝对误差限为 ,相对误差限是 。

2. 测量一支铅笔长是16cm , 那么测量的绝对误差限是 ,测量的相对误差限是 。

3. 称量一件商品的质量为50千克,则其绝对误差限为 ,相对误差限是 。 2.利用平方差的方法

4. 在数值计算中,当a _____________

5. 在数值计算中,计算356-应变成 来计算。

6. 在数值计算中,计算1cos3-

应变为 来计算。 3.f 的位数与f (x )的最高次相同的话,就是最高位的常数,大于的话为0

7. 若543()2792100f x x x x x =-+-+,则12345[1,4,4,4,4,4]f =______________,

123456[1,3,3,3,3,3,3]f = 。

8. 函数()f x 关于三个节点012,,x x x 的拉格朗日二次插值多项式为

3.f(x)=f(x0)[(x-x1)(x-x2)/(x0-x1)(x0-x2)] ,

4(,)n B f x =∑f (k/n )Pk(x)=x

9. 当()f x x =时,(,)n B f x = 。

10. 代数式222236()66

x x

R x x x +=++ ______________,

323222122

()23

x x R x x x ++=++ __________________.

11. 已知方程组1231231

23103127322115

x x x x x x x x x --=-??

-++=??+-=-?,那么收敛的Jacobi 迭代格式为:

,收敛的G S -迭代格式为:

收敛理由是 严格对角占优矩阵 ,

12. 已知线性方程组1233111193234184x x x -????????????-=????????????-??????

,那么收敛的Jacobi 迭代格式:

12.化为线性方程2.调整排序

收敛的G-S 迭代格式: 。

收敛理由是 严格对角占优矩阵 ,

13. 求积公式0

()n

n k

k k I A

f x ==∑至少有n 次代数精度的充要条件是________它是插值型

____________;

当n 是偶数时,牛顿-柯特斯公式()0

()()n

n n k

k k I b a C

f x ==-∑至少有___n+1__P103_____次代

数精度;

高斯求积公式

()()()n

b

k k a

k f x x dx A f x ρ=≈∑?

至少有_____2n+1_P116____次代数精度。

14. 设7227227227n n A R ???????

??=∈??

??

????

,则矩阵A 的特征值的界为 (2.2)与7的和、差为界 ,矩阵1

A -的特征值的界为 界的倒数 。

A ∞=max (1<=i<=n )∑(j(1,n))|aij|等价于每一列中最大值的和 1A = max (1<=j<=n )∑(i(1,n))|aij|等价于每一行中最大值的和 2A =作业第五章12 P116

15. 已知1235A -??=??-??,314x ????=-??????

,那么A ∞

1A = ________,

2A ∞12x = ________,

其中相等的范数有_______A ∞=____________1x =__________. 二、判断题

1. 如果插值节点01,,...,n x x x 互不相同,则满足插值条件的n 次插值多项式是存在且唯一。( x )

2. 迭代改善法能够解决一切方程组的病态问题。 ( x )

3. 区间[,]a b 上的三次样条插值函数()S x 在[,]a b 上具有直到三阶的连续函数。( )

4. 已知1 2.53 3.5A -??=?

?--??

,51x -??

=????,那么1A =1x 。 ( 1 )

5. 法来完成。 ( 1 )

6. 插值法是函数逼近、数值微分和微分方程数值解的基础。 ( 1 )

7.。 ( )

8. 在使用松弛法(SOR )解线性代数方程组AX b =时,若松弛因子ω满足11ω-≥,则迭代法一定不收敛。 ( 1 ) 9. 求解单变量非线性方程()0f x =,弦截法具有1.618阶收敛,抛物线法具有1.840阶收敛,牛顿法具有2阶收敛。 ( 1) 10. 常微分方程初值问题数值解法的理论根据是函数的泰勒展开。 ( 1 ) 11. 解单变量非线性方程()0f x =,牛顿法在单根附近具有2阶收敛,若再用Steffensen 迭代法,则为3阶收敛。 ( 1 ) 三、计算解答题和证明题

构造差分表,用三点牛顿插值多项式,求12.0e和72.0e的近似值。

1.列出牛顿的插值表

2.Px=f(x0)+……P32

3、试用最小二乘法求一次多项式拟合以上数据,并求出均方误差,某实验数据如下:P75

(2)

4、二分法求根 作业第七章1

(1) 方程034

=-+x x 在[1,2]附近的根,使绝对误差不超过0.01(绝对误差估计式:

1

1*2

++-≤

-n n a

b x x ); (2) 方程32()33f x x x x =+--在[1,2]附近的根,使绝对误差不超过0.01;

(3) 方程4

210x x +-=,在[-2,-1]附近的根,使绝对误差不超过0.01。

第六章

5、用适当的方法解方程组:(1)1231231

23424421051145219

x x x x x x x x x ++=??

++=??++=-?;

(2)123310413150134x x x ??????

??????=??????

????????????

; (3)123210315210231x x x -????????????-=????????????-??????. 作业第四章14

6、写出复合梯形公式、复合辛普生公式、复合柯特斯公式及龙贝格公式之间的关系,并用龙贝格方法计算积分

dx x

?

3

1

1

,误差限不超过310-。 7、写出复化梯形公式、复化辛普生公式、复化柯特斯公式及龙贝格公式关系式,并计算积分

cos x e xdx π

?

,已知134.778519T =-,217.389259T =-,

413.336023T =-,812.382162T =-

8、设方程组12341234

123412347239210683032910

x x x x x x x x x x x x x x x x -+-=-??++-=??--+-=??-+-=?,写出Jacobi 迭代法和G S -迭代法的迭代格式,

并证明它们是收敛的。

9、对常微分方程初值问题 (1)(0)1(01)dy

y

dx y x ?=-???=≤≤?

(2)2(0)1(01)dy y dx y x ?=-???=≤≤?(3)12

(0)1(01)

dy

y dx y x ?=-???=≤≤? 取步长0.1,h = 分别用Euler 法和标准的四阶Runge-Kutta 法作数值计算,列表写出结果,并与准确值比较。

10

11、设A 是正交矩阵,证明2()1Cond A =。 12、(1)当()f x x =时,(,)n B f x x =; (2)1()k k k k k k f g f g g f +?=?+?; (3)如果A 是正交阵,则2()1cond A =。 13、证明:适当选取待定参数a , 求积公式

)]()0([)]()0([2

)(''20

h f f ah h f f h

dx x f h

-++≈?

的代数精度可达到3=m 。

14、试证明:适当选取待定参数0A , 1A ,2A ,求积公式

)2()()0()(2130

0h f A h f A f A dx x f h

++≈?

的代数精度可达到2=m 。

15、证明Chebyshev 多项式()n T x 满足微分方程

2'''2

(1)()()()0n n n x T x xT x n T x --+=。

常微分方程习题及答案

第十二章 常微分方程 (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 += 2ln 2 1 (C 为任意常数)。( ) 6.y y sin ='是一阶线性微分方程。( ) 7.xy y x y +='33不是一阶线性微分方程。( ) 8.052=+'-''y y y 的特征方程为0522=+-r r 。( ) 9. 221xy 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.y 1 = 所满足的微分方程是 。

8.x y y 2='的通解为 。 9. 0=+x dy y dx 的通解为 。 10.()2511 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.微分方程3 23y 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 .2 2x e y ?= D . x e y ?=3 8.微分方程x y y sin =+''的一个特解具有形式( )。 A .x a y sin *= B .x a y cos *?= C .()x b x a x y cos sin *+= D . x b x a y sin cos *+= 9.下列微分方程中,( )是二阶常系数齐次线性微分方程。

常微分方程知识点总结

常微分方程知识点总结 常微分方程知识点你学得怎么样呢?下面是的常微分方程知识 点总结,欢迎大家阅读! 微分方程的概念 方程对于学过中学数学的人来说是比较熟悉的;在初等数学中 就有各种各样的方程,比如线性方程、二次方程、高次方程、指数方程、对数方程、三角方程和方程组等等。这些方程都是要把研究的问题中的已知数和数之间的关系找出来,列出包含一个数或几个数的一个或者多个方程式,然后取求方程的解。 但是在实际工作中,常常出现一些特点和以上方程完全不同的 问题。比如:物质在一定条件下的运动变化,要寻求它的运动、变化的规律;某个物体在重力作用下自由下落,要寻求下落距离随时间变化的规律;火箭在发动机推动下在空间飞行,要寻求它飞行的轨道,等等。 物质运动和它的变化规律在数学上是用函数关系来描述的,因此,这类问题就是要去寻求满足某些条件的一个或者几个函数。也就是说,凡是这类问题都不是简单地去求一个或者几个固定不变的数值,而是要求一个或者几个的函数。 解这类问题的基本思想和初等数学解方程的基本思想很相似, 也是要把研究的问题中已知函数和函数之间的关系找出来,从列出的包含函数的一个或几个方程中去求得函数的表达式。但是无论在方程

的形式、求解的具体方法、求出解的性质等方面,都和初等数学中的解方程有许多不同的地方。 在数学上,解这类方程,要用到微分和导数的知识。因此,凡是表示函数的导数以及自变量之间的关系的方程,就叫做微分方程。 微分方程差不多是和微积分同时先后产生的,苏格兰数学家耐普尔创立对数的时候,就讨论过微分方程的近似解。牛顿在建立微积分的同时,对简单的微分方程用级数来求解。后来瑞士数学家雅各布?贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论。 常微分方程的形成与发展是和力学、天文学、物理学,以及其他科学技术的发展密切相关的。数学的其他分支的新发展,如复变函数、李群、组合拓扑学等,都对常微分方程的发展产生了深刻的影响,当前计算机的发展更是为常微分方程的应用及理论研究提供了非常 有力的工具。 牛顿研究天体力学和机械力学的时候,利用了微分方程这个工具,从理论上得到了行星运动规律。后来,法国天文学家勒维烈和英国天文学家亚当斯使用微分方程各自计算出那时尚未发现的海王星 的位置。这些都使数学家更加深信微分方程在认识自然、改造自然方面的巨大力量。 微分方程的理论逐步完善的时候,利用它就可以精确地表述事物变化所遵循的基本规律,只要列出相应的微分方程,有了解方程的方法。微分方程也就成了最有生命力的数学分支。

常微分方程第三版答案

常微分方程第三版答案 Document serial number【KK89K-LLS98YT-SS8CB-SSUT-SST108】

习题 1. dx dy =2xy,并满足初始条件:x=0,y=1的特解。 解: y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2 x +e c =cex 2另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2 x . 2. y 2dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。 解:y 2dx=-(x+1)dy 2 y dy dy=-1 1+x dx 两边积分: - y 1 =-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解:y= | )1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+31 x x + y y 21+dy=31 x x +dx 两边积分:x(1+x 2)(1+y 2)=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1 +dx 两边积分:ln|xy|+x-y=c 另外 x=0,y=0也是原方程的解。

5.(y+x )dy+(x-y)dx=0 解:原方程为: dx dy =-y x y x +- 令 x y =u 则dx dy =u+x dx du 代入有: -112++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x | |-2)(1x y - 则令 x y =u dx dy =u+ x dx du 2 11u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为: tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny= x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32 +=0 解:原方程为:dx dy =y e y 2 e x 3 2 e x 3-3e 2 y -=c.

(完整版)常微分方程的大致知识点

= + ?x = + ?x = + ?x 常微分方程的大致知识点 (一)初等积分法 1、线素场与等倾线 2、可分离变量方程 3、齐次方程(一般含有 x 或 y 的项) y x 4、一阶线性非齐次方程 常数变易法,或 y = e ? a ( x )dx [? b (x )e -? a ( x )dx dx + C ] 5、伯努力方程 令 z = y 1-n ,则 dz = (1 - n ) y -n dy ,可将伯努力方程化成一阶线性非齐次或一阶线性齐次 dx 6、全微分方程 若?M ?y 若 ?M ?y dx = ?N ,则u (x , y ) = C ,(留意书上公式) ?x ≠ ?N ,则找积分因子,(留意书上公式) ?x f (x f ( y , (二)毕卡序列 x y 1 y 0 0 x f (x , y 0 )dx , y 2 y 0 0 x f (x , y 1 )dx , y 3 y 0 0 f (x , y 2 )dx ,其余类推 (三)常系数方程 1、常系数齐次L (D ) y = 0 方法:特征方程 7、可降阶的二阶微分方程 d 2 y = , dy ) ,令 dy = d 2 y p ,则 = dy dx 2 d 2 y = dx dy ) ,令 dx dy = p ,则 dx 2 d 2 y dx = p dp dx 2 dx dx dx 2 dy 8、正交轨线族

? ? dy 单的实根, , y = C e 1x + C e 2 x 1 2 1 2 单的复根1, 2 = ± i , y = e x (C cos x + C 2 sin x ) 重的实根 = = , y = (C + C x )e x 1 2 1 2 重的复根1, 2 = ± i ,3, 4 = ± i , y = e x [(C + C 2 x ) c os x + (C 3 + C 4 x ) sin x ] 2、常系数非齐次L (D ) y = 方法:三部曲。 f (x ) 第一步求L (D ) y = 0 的通解Y 第二步求L (D ) y = f (x ) 的特解 y * 第三步求L (D ) y = f (x ) 的通解 y = Y + y * 如何求 y * ? 当 f (x ) = P m (x )e x 时, y * = x k Q (x )e x 当 f (x ) = P m (x )e ux cos vx + Q (x )e ux sin vx 时, y * = x k e ux (R (x ) cos vx + S m (x ) sin vx ) 当 f (x ) 是一般形式时, y * = ? x W (x ,) f ()d ,其中 W(.)是郎斯基行列式 x 0 W () (四)常系数方程组 方法:三部曲。 第一步求 dX dt = A (t ) X 的通解, Φ(t )C 。利用特征方程 A - I = 0 ,并分情况讨论。 第二步求 dX dt 第三步求 dX dt = A (t ) X + f (t ) 的特解, Φ(t )?Φ-1 (s ) f (s )ds ,(定积分与不定积分等价) = A (t ) X + f (t ) 的通解, Φ(t )C + Φ(t )?Φ-1 (s ) f (s )ds (五)奇点与极限环 ? dx = ax + b y dt ? ? = cx + dy 1、分析方程组? dt 的奇点的性质,用特征方程: A - I = 0 特征方程的根有 3 种情况:相异实根、相异复根、相同实根。第一种情况:相异实根,1 ≠ 2 1 1 m m m

2.5常微分方程课后答案(第三版)王高雄

习题2.5 2.ydy x xdy ydx 2=- 。 解: 2x ,得: ydy x xdy ydx =-2 c y x y d +-=221 即c y x y =+2 2 1 4. xy x y dx dy -= 解:两边同除以x ,得 x y x y dx dy - =1 令u x y = 则dx du x u dx dy += 即 dx du x u dx dy +=u u -=1 得到 ()2ln 2 1 1y c u -=, 即2 ln 21?? ? ??-=y c y x 另外0=y 也是方程的解。 6.()01=-+xdy ydx xy 解:0=+-xydx xdy ydx x d x y x d y y d x -=-2 得到c x y x d +-=??? ? ??2 21

即 c x y x =+2 2 1 另外0=y 也是方程的解。 8. 32 x y x y dx dy += 解:令 u x y = 则: 21u x u dx du x u dx dy +=+= 即2 1u x dx du x = 得到22x dx u du = 故c x u +-=-11 即 21 1x x c y += 另外0=y 也是方程的解。 10. 2 1?? ? ??+=dx dy dx dy x 解:令 p dx dy = 即p p x 2 1+= 而 p dx dy =故两边积分得到 c p p y +-=ln 2 12 因此原方程的解为p p x 21+=,c p p y +-=ln 212 。 12.x y xe dx dy e =?? ? ??+-1 解: y x xe dx dy +=+1

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

常微分方程练习试卷 一、 填空题。 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 -+-=是恰当方程,并求出它的通解.

常微分方程的大致知识点

常微分方程的大致知识点Last revision on 21 December 2020

常微分方程的大致知识点 (一)初等积分法 1、线素场与等倾线 2、可分离变量方程 3、齐次方程(一般含有x y y x 或的项) 4、一阶线性非齐次方程 常数变易法,或])([)()(?+??=-C dx e x b e y dx x a dx x a 5、伯努力方程 令n y z -=1,则dx dy y n dx dz n --=)1(,可将伯努力方程化成一阶线性非齐次或一阶线性齐次 6、全微分方程 若x N y M ??=??,则C y x u =),(,(留意书上公式) 若 x N y M ??≠??,则找积分因子,(留意书上公式) 7、可降阶的二阶微分方程 ),(22dx dy x f dx y d =,令dx dy dx y d p dx dy ==22,则 ),(22dx dy y f dx y d =,令dy dp p dx y d p dx dy ==22,则 8、正交轨线族 (二)毕卡序列 ?+=x x dx y x f y y 0),(001,?+=x x dx y x f y y 0),(102,?+=x x dx y x f y y 0),(203,其余类推 (三)常系数方程 1、常系数齐次0)(=y D L 方法:特征方程 单的实根21,λλ,x x e C e C y 2121λλ+= 单的复根i βαλ±=2,1,)sin cos (21x C x C e y x ββα+= 重的实根λλλ==21,x e x C C y λ)(21+= 重的复根i βαλ±=2,1,i βαλ±=4,3,]sin )(cos )[(4321x x C C x x C C e y x ββα+++=

常微分方程第三版答案2.1

常微分方程习题2.1 1. xy dx dy 2=,并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得 。 故它的特解为代入得 把即两边同时积分得:e e x x y c y x x c y c y xdx dy y 2 2 ,11,0,ln ,21 2 =====+== ,0)1(.22 =++dy x dx y 并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得: 。 故特解是 时,代入式子得。当时显然也是原方程的解当即时,两边同时积分得;当x y c y x y x c y c y x y dy dx x y ++=====++=+=+≠=+- 1ln 11 ,11,001ln 1 ,11ln 0,1112 3 y xy dx dy x y 32 1++ = 解:原式可化为: x x y x x y x y x y y x y c c c c x dx x dy y y x y dx dy 2 22 2 22 2 2 3 22 3 2 )1(1)1)(1(),0(ln 1ln 2 1ln 1ln 2 1 1 1,0111=++ =++ ≠++-=+ +=+≠+ ? + =+)故原方程的解为(即两边积分得故分离变量得显然

.0;0;ln ,ln ,ln ln 0 110000 )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 故原方程的解为即两边积分时,变量分离是方程的解,当或解:由:

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

常微分方程初值问题数值解法 朱欲辉 (浙江海洋学院数理信息学院, 浙江舟山316004) [摘要]:在常微分方程的课程中讨论的都是对一些典型方程求解析解的方法.然而在生产实 际和科学研究中所遇到的问题往往很复杂, 在很多情况下都不可能给出解的解析表达式. 本篇文章详细介绍了常微分方程初值问题的一些数值方法, 导出了若干种数值方法, 如Euler法、改进的Euler法、Runge-Kutta法以及线性多步法中的Adams显隐式公式和预测校正 公式, 并且对其稳定性及收敛性作了理论分析. 最后给出了数值例子, 分别用不同的方法计算出近似解, 从得出的结果对比各种方法的优缺点. [关键词]:常微分方程;初值问题; 数值方法; 收敛性; 稳定性; 误差估计 Numerical Method for Initial-Value Problems Zhu Yuhui (School of Mathematics, Physics, and Information Science, Zhejiang Ocean University, Zhoushan, Zhejiang 316004) [Abstract]:In the course about ordinary differential equations, the methods for analytic solutions of some typical equations are often discussed. However, in scientific research, the problems are very complex and the analytic solutions about these problems can’t be e xpressed explicitly. In this paper, some numerical methods for the initial-value problems are introduced. these methods include Euler method, improved Euler method, Runge-Kutta method and some linear multistep method (e.g. Adams formula and predicted-corrected formula). The stability and convergence about the methods are presented. Some numerical examples are give to demonstrate the effectiveness and accuracy of theoretical analysis. [Keywords]:Ordinary differential equation; Initial-value problem; Numerical method; Convergence; Stability;Error estimate

常微分方程课后答案

习题 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 上的连续非负函数,

2018年电大第三版常微分方程答案知识点复习考点归纳总结参考

习题1.2 1.dx dy =2xy,并满足初始条件:x=0,y=1的特解。 解:y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2x +e c =cex 2另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2x . 2. y 2dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。 解:y 2dx=-(x+1)dy 2y dy dy=-11+x dx 两边积分: -y 1=-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解:y=|)1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+31x x + y y 21+dy=31x x +dx 两边积分:x(1+x 2)(1+y 2)=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1+dx 两边积分:ln|xy|+x-y=c 另外 x=0,y=0也是原方程的解。 5.(y+x )dy+(x-y)dx=0 解:原方程为:

dx dy =- y x y x +- 令x y =u 则dx dy =u+x dx du 代入有: -1 12++u u du=x 1 dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x ||-2)(1x y - 则令x y =u dx dy =u+ x dx du 211u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为:tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny=x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32+=0 解:原方程为:dx dy =y e y 2e x 3 2 e x 3-3e 2y -=c. 9.x(lnx-lny)dy-ydx=0 解:原方程为: dx dy =x y ln x y 令 x y =u ,则dx dy =u+ x dx du

《常微分方程》第三版答案

《常微分方程》第三版答案 习题1.2 1. dx dy =2xy,并满足初始条件:x=0,y=1的特解。解: y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2 x +e c =cex 2 另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时c=1 特解为y= e 2 x . 2. y 2 dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。解:y 2dx=-(x+1)dy 2 y dy dy=-1 1+x dx 两边积分: - y 1=-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解x=0,y=1时c=e 特解:y= | )1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+31 x x + y y 21+dy=31 x x +dx 两边积分:x(1+x 2 )(1+y 2 )=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1 +dx 两边积分:ln|xy|+x-y=c 另外x=0,y=0也是原方程的解。 5.(y+x )dy+(x-y)dx=0 解:原方程为: dx dy =-y x y x +-

令 x y =u 则dx dy =u+x dx du 代入有: -1 12++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即ln(y 2+x 2)=c-2arctg 2 x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x | |-2)(1x y - 则令 x y =u dx dy =u+ x dx du 2 11u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为:tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny= x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32+=0 解:原方程为:dx dy =y e y 2 e x 3 2 e x 3-3e 2 y -=c.

常微分方程(第三版)课后答案

常微分方程 1. xy dx dy 2=,并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得 。 故它的特解为代入得 把即两边同时积分得:e e x x y c y x x c y c y xdx dy y 2 2 ,11,0,ln ,21 2 =====+== ,0)1(.22 =++dy x dx y 并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得: 。 故特解是 时,代入式子得。当时显然也是原方程的解当即时,两边同时积分得;当x y c y x y x c y c y x y dy dx x y ++=====++=+=+≠=+- 1ln 11 ,11,001ln 1 ,11ln 0,1112 3 y xy dx dy x y 32 1++ = 解:原式可化为:

x x y x x y x y x y y x y c c c c x dx x dy y y x y dx dy 2 2 2 2 2 2 2 2 322 32)1(1)1)(1(),0(ln 1ln 21ln 1ln 2 1 1 1,0111=++ =++ ≠++-=+ +=+≠+?+=+) 故原方程的解为(即两边积分得故分离变量得显然 .0;0;ln ,ln ,ln ln 0 110000 )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,0 ln 0 )ln (ln :931:8. cos ln sin ln 0 7ln sgn arcsin ln sgn arcsin 1 sgn 11,)1(,,,6ln )1ln(2 11 11,11,,,0 )()(:5332 2 22 2 22 2 22 2 c dx dy dx dy x y cy u d u u dx x x y u dx x y dy x y ydx dy y x x c dy y y y y dx dy c x y tgxdx ctgydy ctgxdy tgydx c x x x y c x x u dx x x du x dx du dx du x u dx dy ux y u x y y dx dy x c x arctgu dx x du u u u dx du x u dx du x u dx dy ux y u x y x y x y dx dy dx x y dy x y e e e e e e e e x y u u x y x u u x y x y y x x x +===+=+-===-?-=--+-=-=+-===-=+?=+?=?=--=+===-+=+-=++ =++-++=++===+-==-++-+-- 两边积分解:变量分离:。 代回原变量得:则有:令解:方程可变为:解:变量分离,得 两边积分得:解:变量分离,得::也是方程的解。 另外,代回原来变量,得两边积分得:分离变量得:则原方程化为: 解:令:。两边积分得:变量分离,得:则令解:

常微分方程解题方法总结.doc

常微分方程解题方法总结 来源:文都教育 复习过半, 课本上的知识点相信大部分考生已经学习过一遍 . 接下来, 如何将零散的知 识点有机地结合起来, 而不容易遗忘是大多数考生面临的问题 . 为了加强记忆, 使知识自成 体系,建议将知识点进行分类系统总结 . 著名数学家华罗庚的读书方法值得借鉴, 他强调读 书要“由薄到厚、由厚到薄”,对同学们的复习尤为重要 . 以常微分方程为例, 本部分内容涉及可分离变量、 一阶齐次、 一阶非齐次、 全微分方程、 高阶线性微分方程等内容, 在看完这部分内容会发现要掌握的解题方法太多, 遇到具体的题 目不知该如何下手, 这种情况往往是因为没有很好地总结和归纳解题方法 . 下面以表格的形 式将常微分方程中的解题方法加以总结,一目了然,便于记忆和查询 . 常微分方程 通解公式或解法 ( 名称、形式 ) 当 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 一 阶 线 性 微 分 方 程 P ( x)dx P ( x) dx dy Q(x) y ( e Q( x)dx C )e P( x) y dx

伯努利方程 解法:令 u y1 n,有 du (1 n) y n dy , dy P( x) y Q( x) y n(n≠0,1)代入得到du (1 n) P(x)u (1 n)Q(x) dx dx 求解特征方程:2 pq 三种情况: 二阶常系数齐次线性微分方程 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) 两个相等实根:1 2 通解: y c1 c2 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]

常微分方程初值问题

常微分方程初值问题 12.1引言 在数学模型中经常出现的常微分方程在科学的许多分支中同样出现,例如工程和经济学。不幸的是却很少出现这些方程可得到表示在封闭的形式的解的情况,所以通常采用数值方法来寻找近似解。如今,这通常可以非常方便的达到高精度和在解析解和数值逼近之间可靠的误差界。在本节我们将关注一阶微分方程(12.1)形式关于实值函数y的实变 量x的结构和数值分析方法,其中和f是一个给定的实值函数的两个变量。为了从解曲线的无限族选择一个特定的积分构成(12.1)的通解,微分方程将与初始条件一起考虑:给定两个实数和,我们寻求一个(12.1)的解决方案,对于有 (12.2) 微分方程(12.1)与初始条件(12.2)被称为一个初值问题。如果你认为任何(12.1),(12.2)形式的初始值问题具有一个唯一解,看看以下例子。 例12.1考虑微分方程,初始条件,其中α是一个固定的实数,α∈(0,1)。 这是一个关于上述想法的简单验证,对于任何非负实数C, 是初值问题在区间[ 0,∞)上的一个解。因此解的存在性是肯定的,但解不一定唯一;事实上,初始值问题的解有一个无限族,当参数。 我们注意到,在与α∈(0,1)相反的情况下,当α≥1,初值问题,具有唯一解y(x)≡0。 例12.1表明函数f必须遵循相对于它的第二个参数的一定的增长性条件,以保证(12.1),(12.2)有唯一解。精确的保证初始值问题(12.1),(12.2)假设f解的存在惟一基于下面的定理。 定理12.1(Picard theorem)假定实值函数是连续的矩形区域D定义 ;当时;且f 满足Lipschitz条件:存在L>0则 。

常微分方程第三版的课后答案

常微分方程 2.1 1. xy dx dy 2=,并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得 。 故它的特解为代入得 把即两边同时积分得:e e x x y c y x x c y c y xdx dy y 2 2 ,11,0,ln ,21 2 =====+== ,0)1(.22 =++dy x dx y 并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得: 。 故特解是 时,代入式子得。当时显然也是原方程的解当即时,两边同时积分得;当x y c y x y x c y c y x y dy dx x y ++=====++=+=+≠=+- 1ln 11 ,11,001ln 1 ,11ln 0,1112 3 y xy dx dy x y 32 1++ = 解:原式可化为: x x y x x y x y x y y x y c c c c x dx x dy y y x y dx dy 2 2 2 2 22 2 2 3 22 3 2 )1(1)1)(1(),0(ln 1ln 21ln 1ln 2 1 1 1,0111=++ =++ ≠++-=+ +=+≠+ ? + =+) 故原方程的解为(即两边积分得故分离变量得显然 .0;0;ln ,ln ,ln ln 0 110000 )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,0 ln 0 )ln (ln :931:8. cos ln sin ln 0 7ln sgn arcsin ln sgn arcsin 1 sgn 11,)1(,,,6ln )1ln(2 11 11,11,,,0 )()(:5332 2 22 2 22 2 22 2 c dx dy dx dy x y cy u d u u dx x x y u dx x y dy x y ydx dy y x x c dy y y y y dx dy c x y tgxdx ctgydy ctgxdy tgydx c x x x y c x x u dx x x du x dx du dx du x u dx dy ux y u x y y dx dy x c x arctgu dx x du u u u dx du x u dx du x u dx dy ux y u x y x y x y dx dy dx x y dy x y e e e e e e e e x y u u x y x u u x y x y y x x x +===+=+-===-?-=--+-=-=+-===-=+?=+?=?=--=+===-+=+-=++ =++-++=++===+-==-++-+-- 两边积分解:变量分离:。 代回原变量得:则有:令解:方程可变为:解:变量分离,得 两边积分得:解:变量分离,得::也是方程的解。 另外,代回原来变量,得两边积分得:分离变量得:则原方程化为: 解:令:。两边积分得:变量分离,得:则令解:

常微分方程期末试题知识点复习考点归纳总结参考

期末考试 一、填空题(每空2 分,共16分)。 1.方程22d d y x x y +=满足解的存在唯一性定理条件的区域是 . 2. 方程组 n x x x R Y R Y F Y ∈∈=,),,(d d 的任何一个解的图象是 维空间中的一条积分曲线. 3.),(y x f y '连续是保证方程),(d d y x f x y =初值唯一的 条件. 4.方程组???????=-=x t y y t x d d d d 的奇点)0,0(的类型是 5.方程2)(2 1y y x y '+'=的通解是 6.变量可分离方程()()()()0=+dy y q x p dx y N x M 的积分因子是 7.二阶线性齐次微分方程的两个解)(1x y ?=,)(2x y ?=成为其基本解组的充要条件是 8.方程440y y y '''++=的基本解组是 二、选择题(每小题 3 分,共 15分)。 9.一阶线性微分方程 d ()()d y p x y q x x +=的积分因子是( ). (A )?=x x p d )(e μ (B )?=x x q d )(e μ (C )?=-x x p d )(e μ (D )?=-x x q d )(e μ 10.微分方程0d )ln (d ln =-+y y x x y y 是( ) (A )可分离变量方程 (B )线性方程 (C )全微分方程 (D )贝努利方程 11.方程x (y 2-1)d x+y (x 2-1)d y =0的所有常数解是( ). (A) 1±=x (B)1±=y

(C )1±=y , 1±=x (D )1=y , 1=x 12.n 阶线性非齐次微分方程的所有解( ). (A )构成一个线性空间 (B )构成一个1-n 维线性空间 (C )构成一个1+n 维线性空间 (D )不能构成一个线性空间 13.方程222+-='x y y ( )奇解. (A )有一个 (B )有无数个 (C )只有两个 (D )无 三、计算题(每小题8分,共48分)。 14.求方程22 2d d x y xy x y -=的通解 15.求方程0d )ln (d 3=++y x y x x y 的通解 16.求方程2 221)(x y x y y +'-'=的通解

相关文档
最新文档