跳到主要内容

京都大学 情報学研究科 数理工学専攻 2014年8月実施 力学系数学

Author

Casablanca

Description

日本語版

a(t)a(t) を半無限区間 [1,)[1, \infty) で定義された連続関数として、微分方程式

td2xdt2+(t21)a(t)dxdtta(t)x=0,t1\begin{align} t \frac{d^2 x}{d t^2} + (t^2 - 1) a(t) \frac{dx}{dt} - t a(t) x = 0, \qquad t \geqq 1 \tag{1} \end{align}

を考える。x=t3t2+1x = \frac{t^3}{t^2 + 1} を式 (1) のひとつの解とする。以下の問いに答えよ。

(i) 関数 a(t)a(t) を求めよ。

(ii) 式 (1) で x=t3t2+1yx = \frac{t^3}{t^2 + 1} y とおく。yy が満たす微分方程式を求めよ。

(iii) (ii) で得られた微分方程式を解き、式 (1) の一般解を求めよ。

(iv) 式 (1) の解 x(t)x(t) が半無限区間 [1,)[1, \infty) において有界となるための、t=1t=1 における初期値 (x0,v0)=(x(1),dxdt(1))(x_0, v_0) = \left( x(1), \frac{dx}{dt} (1) \right) の必要十分条件を求めよ。

English Version

Kai

(i)

dxdt=t4+3t2t4+2t2+1,d2xdt2=2t(3t2)(t2+1)3\frac{dx}{dt} = \frac{t^4 + 3t^2}{t^4 + 2t^2 + 1},\quad \frac{d^2x}{dt^2} = \frac{2t(3-t^2)}{(t^2 + 1)^3}

and we get

a(t)=2t2+1a(t) = \frac{2}{t^2 + 1}

(ii)

for solution x1,x2x_1, x_2

(x1x2x1x2)=x1x2x1x2=(t1t)a(t)(x1x2x1x2)\begin{aligned} (x_1'x_2 - x_1 x_2')' &= x_1 ''x_2 - x_1 x_2'' \\ &= -(t-\frac 1t )a(t)(x_1'x_2 - x_1 x_2') \end{aligned}

x2=x1yx_2 = x_1 y , x2=x1y+x1y1x_2 ' = x_1 ' y + x_1 y_1 ', then

(x1x1yx1(x1y+x1y))=2(1t2)t(t2+1)x12y(x_1' x_1 y - x_1 (x_1 ' y + x_1 y'))' = -\frac{2(1-t^2)}{t(t^2 + 1)} x_1^2 y'

and we obtain

ty+4y=0ty'' + 4y' = 0

(iii)

Let y=μy' = \mu, then tμ=4μt\mu ' = -4\mu, we get μ=Ct4\mu = \frac {C}{t^4}, y=C1t3+C2y = \frac{C_1}{t^3} + C_2.

Let C1=1C_1 = 1 , C2=0C_2 = 0, we have

t3t2+1y=1t2+1\frac{t^3}{t^2 + 1}y = \frac{1}{t^2 + 1}

since 1t2+1\frac{1}{t^2 + 1} , t3t2+1\frac{t^3}{t^2 + 1} is linear independent,

C1+C2t3t2+1\frac{C_1 + C_2 t^3}{t^2 + 1}

is a general solution

(iv)

x(t)=C1+C2t3t2+1x(t) = \frac{C_1 + C_2 t^3}{t^2 + 1}
x(t)=C2t4+3C2t22C1t(t2+1)2x'(t) = \frac{C_2 t^4 + 3C_2 t^2 - 2C_1 t}{(t^2 + 1)^2}

C2=0C_2 = 0 is a necessary and sufficient condition. And by solving:

[x(1)dxdt(1)]=[C1+C222C1+4C24]=[x0v0]\begin{bmatrix} x(1)\\ \frac{dx}{dt}(1) \end{bmatrix} = \begin{bmatrix} \frac{C_1 + C_2}{2} \\ \frac{-2C_1 + 4C_2}{4} \end{bmatrix} = \begin{bmatrix} x_0 \\ v_0 \end{bmatrix}

we get

C1=23(2x0v0)C2=23(x0+v0)\begin{aligned} &C_1 = \frac 23 (2x_0 - v_0) \\ &C_2 = \frac 23 (x_0 + v_0) \end{aligned}

obviously x0=v0x_0 = -v_0 is necessary and sufficient.