浙浙浙江江江大大大学学学2014——2015学学学年年年春春春夏夏夏学学学期期期
《《《常常常微微微分分分方方方程程程》》》课课课程程程期期期末末末考考考试试试试试试卷卷卷答答答案案案 一. 求解下列方程(25分)
1. dx dy = x+y 2y
2
解: y = 0是解
当y = 0时,
dxdy=
2y1x +
y2x = Cp|y| + 1
3 y
2, C 为任意常数 所以解是y = 0或x = Cp|y| +
13y
2, C 为任意常数
2.
2y2
+ 5dy dx
2
= 4
解: 设y = √
2 sin t, y
0=
2√ 5
5
cos t, x = x. 则 2 √
5
5 cos t = y
0= d( √ 2 sin t) dx = √
2 cos t dt dx t =
√ 10
5 x + C或 cos t = 0 y = √
2 sin(
√ 10
5 x + C), C 为任意常数 或特解y = ± √
2
3.
x2 d dx
2y
2 − 3xdy dx
+ 5y = x2
sin(ln x)解: 令x = e
td
2y
dt
2− 4 dy
dt + 5y = e
2tsin t 特征方程
λ
2− 4λ + 5 = 0, λ = 2 ± i.
齐次方程解Y = c
1e
2tcos t + c
2e
2tsin t
设非齐次方程解为y
∗= Ate
2tcos t + Bte
2tsin t, 得 到 A = − 1
2 , B = 0.
y = c
1e
2tcos t + c
2e
2tsin t − t
2 e
2tcos t
= c
1x
2cos ln x + c
2x
2sin ln x − ln x
2 x
2cos ln x
二. 求解下列方程(组)(25分)
1.
用幂级数法求解y00
+ 4xy = 0, y(0) = 1, y0
(0) = 0解:设y = P
∞n=0a
nx
n2a
2+
∞
X
n=1
[a
n+2(n + 2)(n + 1) + 4a
n−1]x
n= 0 所以a
2= 0, a
n+2=
(n+2)(n+1)−4an−1, n = 1, 2, . . . 则
a
3k= (−1)
k4
k(3k)(3k − 1) . . . (3 · 2) a
0= (−1)
k4
k(3k − 2)!!!
(3k)! a
0a
3k+1= (−1)
k4
k(3k + 1)(3k) . . . (4 · 3) a
1= (−1)
k4
k(3k − 1)!!!
(3k + 1)! a
1a
3k+2= 0, k = 1, 2, . . .
因为y(0) = 1, y
0(0) = 0, 所 以a
0= 1, a
1= 0 y = 1 +
∞
X
k=1
(−1)
k4
k(3k − 2)!!!
(3k)! x
3k2.
(dx
dt
= 2y − 5x + e−t
(1)dy
dt
= x − 6y + e−2t
(2)解:(1)得到
y = 1
2 (x
0+ 5x − e
−t) 代入(2),
x
00+ 11x
0+ 28x = 5e
−t+ 2e
−2t特征方程λ
2+ 11λ + 28 = 0
λ
1= −4, λ
2= −7 齐次方程解
X = c
1e
−4t+ c
2e
−7t设x
00+ 11x
0+ 28x = 5e
−t解为x
∗1= Ae
−t,则A =
185, x
∗1=
185e
−t设x
00+ 11x
0+ 28x = 2e
−2t解为x
∗2= Be
−2t,则B =
15, x
∗2=
15e
−2tx = c
1e
−4t+ c
2e
−7t+ 5
18 e
−t+ 1 5 e
−2t则
y = c
12 e
−4t− c
2e
−7t+ 1
18 e
−t+ 3 10 e
−2t3.
dx
dt
−dy dt
−dz dt
+ x − 2z = 0 (1)dx
dt
−dy dt
+dz dt
+ x = 0 (2)dx
dt
+dy dt
−dz dt
+ x + 2y = 0 (3)解:(1)+(2)得2x
0− 2y
0+ 2x − 2z = 0 (4) (1)+(3)得2x
0− 2z
0+ 2x + 2y − 2z = 0 (5)
第2页
(2)+(3)得2x
0+ 2x + 2y = 0 (6) (5)-(6)得z
0+ z = 0,
z = c
1e
−t(4)-(5)得−2y
0+ 2z
0− 2y = 0
y
0+ y = −c
1e
−ty = c
2e
−t− c
1te
−t代入(6)得x
0+ x + c
2e
−t− c
1te
−t= 0
x = c
3e
−t+ c
12 t
2e
−t− c
2te
−t三.(20分)对于系统
x0
= 1 − x + y − x2
y
0
= x(x − y)找出所有平衡点(奇点),写出关于这些 平衡点所相应的线性化系统,判断平衡点的类型,并画出平衡点附近相图的草图。
解: 奇点: M
1(0, −1), M
2(1, 1), M
3(−1, −1) M
1(0, −1)所 相应的线性化系统 x
0= −x + y
y
0= x
A
1= −1 1 1 0
特征λ =
−1±2√5, 鞍点
k = 1
−1 + k , k = 1 ± √ 5 2 (1,
1+√ 5
2
) 点x
0> 0, y
0> 0
M
2(1, 1)所 相应的线性化系统 x
0= −3x + y y
0= x − y
A
2= −3 1 1 −1
特征λ = −2 ± √
2, 稳 定结点
k = 1 − k
−3 + k , k = 1 ± √ 2
dy
dx
|
y=x= 0 < 1
M
3(−1, −1)所 相应的线性化系统 x
0= x + y y
0= −x + y
A
3= 1 1
−1 1
特征λ = 1 ± i, 不稳定焦点
xy
0− yx
0= −x
2− y
2< 0 ,顺时针
四.(15分)讨论下面2个方程组零解的李雅普诺夫稳定性 (1) ( x
0= 4y
3− x
3y
0= −4x − y
3(2) ( x
0= −x
4y
y
0= x
3y
2构造
V (x, y) = y
4+ 2x
2,
V
0(x, y) = 4x(4y
3− x
3) + 4y
3(−4x − y
3) = −4x
4− 4y
6渐近稳定 首次积分V = xy
V
0(x, y) = y(−x
4y) + x(x
3y
2) = 0 不稳定或者:区域A = {x > 0, y > 0},构造
V (x, y) = xy
2V
0(x, y) = y
2(−x
4y) + 2xy(x
3y
2) = x
4y
3不稳定或者:区域A = {x > 0, y < 0},构造
V (x, y) = −x
2y
V
0(x, y) = −2xy(−x
4y) − x
2(x
3y
2) = x
5y
2不稳定
五.(15分)给定区间I = [0, a],非负连续函数u(t) ≤ 1,u(0) = 0,连续可微函数f : (t, x) ∈ I × R → R,以及区间[−2, 0]中的一个连续可微函数φ(t),并满足φ
0(0−) = f (0, φ(0)).
考虑如下问题
dxdt
= f (t, x(t − u(t))) t ∈ [0, a]
x(t) = φ(t) t ∈ [−2, 0]
(1) 试证明存在一个α > 0使得该问题在t ∈ [0, α]至少存在一个解。
(2) 更进一步,这样的解是否有唯一性,给出充足的理由。
1. 方程在[−2, T ]的解与积分方程的等价性 x(t) =
φ(0) + R
t0
f (s, x(s − u(s)))ds t ∈ [0, T ]
φ(t) t ∈ [−2, 0]
2. 基本设定。记
M = max
t∈[−2,0]
|φ(t)|
N = max
t∈[−2,a],x∈[−2M,2M ]
|f (t, x)|, L = max
t∈[−2,a],x∈[−2M,2M ]
| ∂f
∂x (t, x)|
B
α= {x ∈ C
1([−2, α]) : |x(t)| ≤ 2M, x(t) = φ(t), t ≤ 0}
3. 存在性。可以用欧拉折线法,或者Picard迭代,或者压缩映射。这里用压缩映射来说明。
对于0 < α ≤ min(
MN,
2L1, a) ,我们定义映射
T : B
α→ C
1([−2, α]) y(t) = (T x)(t) :=
φ(0) + R
t0
f (s, x(s − u(s)))ds, t > 0,
φ(t), t ≤ 0 .
所以
|y(t)| ≤ |y(t) − φ(0)| + |φ(0)| ≤ N α + M ≤ 2M, T : B
α→ B
α另外,如果t ≥ 0,x, y ∈ B
α,kxk = max
t∈[−2,α]|x(t)| ,
|T x − T y(t)| ≤ L Z
t0
|x(s − u(s))) − y(s − u(s))|ds ≤ Lαkx − yk ≤ 1/2kx − yk
4. 唯一性。
需要证明区间[−2, β](0 < β ≤ a)中的两个解x(t), y(t)相等即可。记 M
0= max
t∈[−2,β]
|x(t)| + |y(t)|, L
0= max
t∈[−2,β],x∈[−M0,M0]
| ∂f
∂x (t, x)|
则如果t ∈ [0, β],记非负非减函数kx − yk(t) := max
s∈[−2,t]|x − y(s)| ,则kx − yk(0) = 0,并且
|x − y(t)| ≤ L
0Z
tβ
kx − yk(s)ds ⇒ kx − yk(t) ≤ L
0Z
tβ