常微分方程王高雄第三版答案3.1

习题3.1

1 求方程dy =x+y2通过点(0,0)的第三次近似解; dx

解: 取ϕ0(x ) =0 12x 002

x x 11152x ϕ2(x ) =y 0+⎰[x +ϕ1(x ) ]dx =⎰[x +(x 2) 2]dx =x 2+002220

x 1152x ) ]dx ϕ3(x ) =y 0+⎰[x +(x 2+0220

115181x +x +x 11 = x 2+2201604400

dy 2 求方程=x-y2通过点(1,0)的第三次近似解; dx ϕ1(x ) =y 0+⎰(x +y 0) dx =⎰xdx =x 2x

解: 令ϕ0(x ) =0

12x 002

x x 12212152ϕ(x ) =y +[x -ϕ(x ) ]dx =[x -(x ) ]dx =x -x 201⎰0⎰02220

x 1152x ) ]dx ϕ3(x ) =y 0+⎰[x -(x 2-0220

115181x +x -x 11 =x 2-2201604400

则 ϕ1(x ) =y 0+⎰(x -y 0) dx =⎰xdx =x 2x

3 题 求初值问题:

dy ⎧⎪=x 2-y 2 R:x +≤1,y ≤1 ⎨dx ⎪⎩y (-1) =0

的解的存在区间,并求解第二次近似解,给出在解的存在空间的误差估计;

b 1解: 因为 M=max{x 2-y 2}=4 则h=min(a,)= M 4

1 则解的存在区间为x -x 0=x -(-1) =x +1≤ 4

令 ψ0(X ) =0 ;

11ψ1(x ) =y0+⎰(x 2-0) dx=x 3+; 33x 0x

ψ2(x ) 13x x 4x 7111312=y0+⎰[x -(x +) ]dx=x ---+ 3942186333-12x

又 ∂f (x , y ) ≤2=L ∂y

M *L 2

311则:误差估计为:ψ2(x ) -ψ(x ) ≤h = 24(2+1) 2

dy 334 题 讨论方程:=y 在怎样的区域中满足解的存在唯一性定理的条件, dx 21

并求通过点(0,0)的一切解;

∂f (x , y ) 13解:因为=y 在y ≠0上存在且连续; ∂y 2-2

3 而y 3在y ≥σ 0上连续 2

dy 33由 =y 有:y =(x+c)2 dx 2131

又 因为y(0)=0 所以:y =x

另外 y=0也是方程的解;

⎧3⎪2故 方程的解为:y =⎨x

⎪⎩0x ≥0 x 032

或 y=0;

6题 证明格朗瓦耳不等式:

设K 为非负整数,f(t)和g(t)为区间α≤t ≤β上的连续非负函数,

且满足不等式:

t

f(t)≤k+⎰f (s ) g (s ) ds , α≤t ≤β

α

t

则有:f(t)≤kexp(⎰g (s ) ds ), α≤t ≤β

α

t

证明:令R (t )=⎰f (s ) g (s ) ds , 则R ' (T)=f(t)g(t)

α

R ' (T)-R(t)g(t)= f(t)g(t)- R(t)g(t) ≤kg(t)R ' (T)- R(t)g(t)≤kg(t);

t

两边同乘以exp(-⎰g (s ) ds ) 则有:

α

t

' t R (T) exp(-⎰g (s ) ds )-R(t)g(t) exp(-⎰g (s ) ds )

αα

t

≤ kg(t) exp(-⎰g (s ) ds )

α

两边从α到t 积分:

t t t

R(t) exp(-⎰g (s ) ds ) ≤-⎰kg (s ) ds exp(-⎰g (r ) dr )ds

ααα

t t

即 R(t) ≤⎰kg (s ) ds exp(-⎰g (r ) dr )ds

αs

t t

又 f(t) ≤1≤k+R(t) ≤k+k⎰g (s ) exp(-⎰g (r ) dr )ds

αs

t s

≤k(1-1+ exp(-⎰g (r ) dr )=k exp(⎰g (r ) dr )

s t

t

即 f(t) ≤k ⎰g (r ) dr ;

α

7题 假设函数f(x,y)于(x 0,y 0)的领域内是y 的 不增函数,试证方程

dy = f(x,y)满足条件y(x0)= y0的解于x ≥ x0一侧最多只有一个解; dx

证明:假设满足条件y(x0)= y0的解于x ≥ x0一侧有两个ψ(x),ϕ(x)

则满足:

x

ϕ(x)= y0+⎰f (x , ϕ(x )) dx

x 0

x

ψ(x)= y0+⎰f (x , ψ(x )) dx

x 0

不妨假设ϕ(x) ψ(x),则ϕ(x)- ψ(x)≥0

x x

而ϕ(x)- ψ(x)= ⎰f (x , ϕ(x )) dx-⎰f (x , ψ(x )) dx

x 0x 0

x

=⎰[f (x , ϕ(x )) -f (x , ψ(x )) dx

x 0

又因为 f(x,y)在(x 0,y 0)的领域内是y 的 增函数,则: f(x, ϕ(x))-f(x, ψ(x))≤0

x

则ϕ(x)- ψ(x)= ⎰[f (x , ϕ(x )) -f (x , ψ(x )) dx ≤0

x 0

则ϕ(x)- ψ(x)≤0

所以 ϕ(x)- ψ(x)=0, 即 ϕ(x)= ψ(x) 则原命题方程满足条件y(x0)= y0的解于x ≥ x0一侧最多 只有一个解;

习题3.1

1 求方程dy =x+y2通过点(0,0)的第三次近似解; dx

解: 取ϕ0(x ) =0 12x 002

x x 11152x ϕ2(x ) =y 0+⎰[x +ϕ1(x ) ]dx =⎰[x +(x 2) 2]dx =x 2+002220

x 1152x ) ]dx ϕ3(x ) =y 0+⎰[x +(x 2+0220

115181x +x +x 11 = x 2+2201604400

dy 2 求方程=x-y2通过点(1,0)的第三次近似解; dx ϕ1(x ) =y 0+⎰(x +y 0) dx =⎰xdx =x 2x

解: 令ϕ0(x ) =0

12x 002

x x 12212152ϕ(x ) =y +[x -ϕ(x ) ]dx =[x -(x ) ]dx =x -x 201⎰0⎰02220

x 1152x ) ]dx ϕ3(x ) =y 0+⎰[x -(x 2-0220

115181x +x -x 11 =x 2-2201604400

则 ϕ1(x ) =y 0+⎰(x -y 0) dx =⎰xdx =x 2x

3 题 求初值问题:

dy ⎧⎪=x 2-y 2 R:x +≤1,y ≤1 ⎨dx ⎪⎩y (-1) =0

的解的存在区间,并求解第二次近似解,给出在解的存在空间的误差估计;

b 1解: 因为 M=max{x 2-y 2}=4 则h=min(a,)= M 4

1 则解的存在区间为x -x 0=x -(-1) =x +1≤ 4

令 ψ0(X ) =0 ;

11ψ1(x ) =y0+⎰(x 2-0) dx=x 3+; 33x 0x

ψ2(x ) 13x x 4x 7111312=y0+⎰[x -(x +) ]dx=x ---+ 3942186333-12x

又 ∂f (x , y ) ≤2=L ∂y

M *L 2

311则:误差估计为:ψ2(x ) -ψ(x ) ≤h = 24(2+1) 2

dy 334 题 讨论方程:=y 在怎样的区域中满足解的存在唯一性定理的条件, dx 21

并求通过点(0,0)的一切解;

∂f (x , y ) 13解:因为=y 在y ≠0上存在且连续; ∂y 2-2

3 而y 3在y ≥σ 0上连续 2

dy 33由 =y 有:y =(x+c)2 dx 2131

又 因为y(0)=0 所以:y =x

另外 y=0也是方程的解;

⎧3⎪2故 方程的解为:y =⎨x

⎪⎩0x ≥0 x 032

或 y=0;

6题 证明格朗瓦耳不等式:

设K 为非负整数,f(t)和g(t)为区间α≤t ≤β上的连续非负函数,

且满足不等式:

t

f(t)≤k+⎰f (s ) g (s ) ds , α≤t ≤β

α

t

则有:f(t)≤kexp(⎰g (s ) ds ), α≤t ≤β

α

t

证明:令R (t )=⎰f (s ) g (s ) ds , 则R ' (T)=f(t)g(t)

α

R ' (T)-R(t)g(t)= f(t)g(t)- R(t)g(t) ≤kg(t)R ' (T)- R(t)g(t)≤kg(t);

t

两边同乘以exp(-⎰g (s ) ds ) 则有:

α

t

' t R (T) exp(-⎰g (s ) ds )-R(t)g(t) exp(-⎰g (s ) ds )

αα

t

≤ kg(t) exp(-⎰g (s ) ds )

α

两边从α到t 积分:

t t t

R(t) exp(-⎰g (s ) ds ) ≤-⎰kg (s ) ds exp(-⎰g (r ) dr )ds

ααα

t t

即 R(t) ≤⎰kg (s ) ds exp(-⎰g (r ) dr )ds

αs

t t

又 f(t) ≤1≤k+R(t) ≤k+k⎰g (s ) exp(-⎰g (r ) dr )ds

αs

t s

≤k(1-1+ exp(-⎰g (r ) dr )=k exp(⎰g (r ) dr )

s t

t

即 f(t) ≤k ⎰g (r ) dr ;

α

7题 假设函数f(x,y)于(x 0,y 0)的领域内是y 的 不增函数,试证方程

dy = f(x,y)满足条件y(x0)= y0的解于x ≥ x0一侧最多只有一个解; dx

证明:假设满足条件y(x0)= y0的解于x ≥ x0一侧有两个ψ(x),ϕ(x)

则满足:

x

ϕ(x)= y0+⎰f (x , ϕ(x )) dx

x 0

x

ψ(x)= y0+⎰f (x , ψ(x )) dx

x 0

不妨假设ϕ(x) ψ(x),则ϕ(x)- ψ(x)≥0

x x

而ϕ(x)- ψ(x)= ⎰f (x , ϕ(x )) dx-⎰f (x , ψ(x )) dx

x 0x 0

x

=⎰[f (x , ϕ(x )) -f (x , ψ(x )) dx

x 0

又因为 f(x,y)在(x 0,y 0)的领域内是y 的 增函数,则: f(x, ϕ(x))-f(x, ψ(x))≤0

x

则ϕ(x)- ψ(x)= ⎰[f (x , ϕ(x )) -f (x , ψ(x )) dx ≤0

x 0

则ϕ(x)- ψ(x)≤0

所以 ϕ(x)- ψ(x)=0, 即 ϕ(x)= ψ(x) 则原命题方程满足条件y(x0)= y0的解于x ≥ x0一侧最多 只有一个解;


    相关文章

    数学毕业论文 浅谈积分因子的求解法

    2014届本科毕业论文(设计) 题目:浅谈求积分因子的方法 学 院:数学科学学院 专业班级:数学09-3班 学生姓名:阿曼姑丽·吐鲁甫 指导教师:伊里夏提 答辩日期:2014年5月8日 新疆师范大学教务处 目 录 1. 引言 . ..... ...

    常微分方程中变量变换方法的探讨

    江西师范大学数学与信息科学学院 学士学位论文 常微分方程中变量变换方法的探讨 The Study on Method of Variable-transformed in Ordinary Differential Equations 姓 ...

    一类微分方程组特解的求法

    科技信息 高校理科研究 一类微分方程组特解硇求法 南通职业大学基础课部 周彤 [摘要]对于一阶微分方程组求特解这类问题,可以利用拉普拉斯变换法将原方程组化为二元一次线性方程组:然 后用克莱姆法则和留数法求出特解. [关键词]微分方程组特解拉 ...

    阶非齐次常系数微分方程解的表达式

    2009年10月 第30卷第5期湘南学院学报JournalofXiangnanUniversityOct.,2009V01.30No.5 二阶非齐次常系数微分方程解的表达式 樊自安 (孝感学院数学与统计学院,湖北孝感432000) 摘要:给 ...

    不良导体热导率的测定

    不良导体热导率的测定 彭志伟(第一作者),贾林江(第二作者) 摘要:在稳态法测量不良导体热导率的实验中,传统的通过作图求斜率测定散热速率的方法存在很大的主观性,对实验结果产生的误差是不可控制的,本文通过对散热规律的研究,利用Matlab 多 ...

    关于积分因子的讨论

    第16卷第5期 2006年10月长春大学学报JOURNAI.OFCHANGCHUNUNIVERSITYV01.16No.5Oct.2()016文章编号:1009-3907(2006)05-0017-04 关于积分因子的讨论 温启军1,张丽静 ...

    2012高考天津卷文综地理_及答案_高清__word版

    邹城一中 高二 地理组 2012,6,11 2012年普通高等学校招生全国统一考试(天津卷) 文科综合地理部分 第1卷 注意事项:在每题给出的四个选项中,只有一项是最符合题目要求的. 读城市水循环示意图(图1),回答第1题. 在城镇建设中, ...

    行列式与齐次线性微分方程通解判定

    伏朗斯基行列式与齐次线性微分方程通解判定 作者: 作者单位: 刊名: 英文刊名: 年,卷(期): 被引用次数:李长江, 李淑华, LI Chang-jiang, LI Shu-hua承德民族师专,数学系,河北,承德,067000承德民族师专 ...

    线性微分方程组初值问题的求解公式

    第! ! 卷第" 期宝鸡文理学院学报#自然科学版$ #$, &-. /0'&120&345&''676&18. 9:0/; (! ! ) &(" %&' (! ** ...