• 沒有找到結果。

解: y = 0是解

N/A
N/A
Protected

Academic year: 2021

Share "解: y = 0是解"

Copied!
5
0
0

加載中.... (立即查看全文)

全文

(1)

浙浙浙江江江大大大学学学2014——2015学学学年年年春春春夏夏夏学学学期期期

《常常微微分分方方程程》》课课程程期期末末考考试试试试卷卷答答案 一. 求解下列方程(25分)

1. dx dy

=

x+y 2y

2

解: y = 0是解

当y = 0时,

dxdy

=

2y1

x +

y2

x = Cp|y| + 1

3 y

2

, C 为任意常数 所以解是y = 0或x = Cp|y| +

13

y

2

, C 为任意常数

2.

2y

2

+ 5



dy 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.

x

2 d dx

2

y

2 − 3x

dy dx

+ 5y = x

2

sin(ln x)

解: 令x = e

t

d

2

y

dt

2

− 4 dy

dt + 5y = e

2t

sin t 特征方程

λ

2

− 4λ + 5 = 0, λ = 2 ± i.

齐次方程解Y = c

1

e

2t

cos t + c

2

e

2t

sin t

设非齐次方程解为y

= Ate

2t

cos t + Bte

2t

sin t, 得 到 A = − 1

2 , B = 0.

y = c

1

e

2t

cos t + c

2

e

2t

sin t − t

2 e

2t

cos t

= c

1

x

2

cos ln x + c

2

x

2

sin ln x − ln x

2 x

2

cos ln x

二. 求解下列方程(组)(25分)

1.

用幂级数法求解y

00

+ 4xy = 0, y(0) = 1, y

0

(0) = 0

解:设y = P

n=0

a

n

x

n

(2)

2a

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)

k

4

k

(3k)(3k − 1) . . . (3 · 2) a

0

= (−1)

k

4

k

(3k − 2)!!!

(3k)! a

0

a

3k+1

= (−1)

k

4

k

(3k + 1)(3k) . . . (4 · 3) a

1

= (−1)

k

4

k

(3k − 1)!!!

(3k + 1)! a

1

a

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)

k

4

k

(3k − 2)!!!

(3k)! x

3k

2.

(

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

1

e

−4t

+ c

2

e

−7t

设x

00

+ 11x

0

+ 28x = 5e

−t

解为x

1

= Ae

−t

,则A =

185

, x

1

=

185

e

−t

设x

00

+ 11x

0

+ 28x = 2e

−2t

解为x

2

= Be

−2t

,则B =

15

, x

2

=

15

e

−2t

x = c

1

e

−4t

+ c

2

e

−7t

+ 5

18 e

−t

+ 1 5 e

−2t

y = c

1

2 e

−4t

− c

2

e

−7t

+ 1

18 e

−t

+ 3 10 e

−2t

3.

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页

(3)

(2)+(3)得2x

0

+ 2x + 2y = 0 (6) (5)-(6)得z

0

+ z = 0,

z = c

1

e

−t

(4)-(5)得−2y

0

+ 2z

0

− 2y = 0

y

0

+ y = −c

1

e

−t

y = c

2

e

−t

− c

1

te

−t

代入(6)得x

0

+ x + c

2

e

−t

− c

1

te

−t

= 0

x = c

3

e

−t

+ c

1

2 t

2

e

−t

− c

2

te

−t

三.(20分)对于系统

 x

0

= 1 − x + y − x

2

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±25

, 鞍点

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

3

y

0

= −4x − y

3

(2) ( x

0

= −x

4

y

y

0

= x

3

y

2

(4)

构造

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

4

y) + x(x

3

y

2

) = 0 不稳定或者:区域A = {x > 0, y > 0},构造

V (x, y) = xy

2

V

0

(x, y) = y

2

(−x

4

y) + 2xy(x

3

y

2

) = x

4

y

3

不稳定或者:区域A = {x > 0, y < 0},构造

V (x, y) = −x

2

y

V

0

(x, y) = −2xy(−x

4

y) − x

2

(x

3

y

2

) = x

5

y

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)).

考虑如下问题



dx

dt

= 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

t

0

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

t

0

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

t

0

|x(s − u(s))) − y(s − u(s))|ds ≤ Lαkx − yk ≤ 1/2kx − yk

4. 唯一性。

(5)

需要证明区间[−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

0

Z

t

β

kx − yk(s)ds ⇒ kx − yk(t) ≤ L

0

Z

t

β

kx − yk(s)ds ⇒ kx − yk(t) ≡ 0

(Gronwall不等式)证毕。

參考文獻

相關文件

INFORMAÇÃO GLOBAL SOBRE AS ASSOCIAÇÕES DE SOLIDARIEDADE SOCIAL E OS SERVIÇOS SUBSIDIADOS REGULARMENTE PELO INSTITUTO DE ACÇÃO SOCIAL. STATISTICS ON SOCIAL SOLIDARITY ASSOCIATIONS

EQUIPAMENTO SOCIAL A CARGO DO INSTITUTO DE ACÇÃO SOCIAL, Nº DE UTENTES E PESSOAL SOCIAL SERVICE FACILITIES OF SOCIAL WELFARE BUREAU, NUMBER OF USERS AND STAFF. 數目 N o

INFORMAÇÃO GLOBAL SOBRE AS ASSOCIAÇÕES DE SOLIDARIEDADE SOCIAL E OS SERVIÇOS SUBSIDIADOS REGULARMENTE PELO INSTITUTO DE ACÇÃO SOCIAL. STATISTICS ON SOCIAL SOLIDARITY ASSOCIATIONS

產」五大社區發展面向,而產生出來的政策性名詞。居住在同一地理範圍內的居民,持續以

„ FPGA –現場可規劃邏輯陣列 (field- programmable

z 可規劃邏輯區塊 (programmable logic blocks) z 可規劃內部連接

(In Section 7.5 we will be able to use Newton's Law of Cooling to find an equation for T as a function of time.) By measuring the slope of the tangent, estimate the rate of change

[r]