跳到主要内容

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

Author

Casablanca

Description

日本語版

nn を自然数、ijij 成分が

aij(t)={1(for i=j)t(for i=j+1)0otherwisea_{ij}(t) = \begin{cases} &1 &(\text{for } i=j) \\ &t &(\text{for } i=j+1) \\ &0 &\text{otherwise} \end{cases}

nn 次正方行列を A(t)A(t) として、t>0t > 0 において nn 元連立線形微分方程式

dxdt=1tA(t)x,xRn\frac{dx}{dt} = \frac{1}{t} A(t) x, \quad x \in \mathbb{R}^n

を考える。以下の問いに答えよ。

(i) n=1n = 1 のとき一般解を求めよ。

(ii) n=2n = 2 のとき一般解を求めよ。

(iii) 任意の自然数 nn に対して一般解を求めよ。

English Version

Kai

(i)

dxdt=xt\frac{dx}{dt} = \frac{x}{t}
x=kt,k is a constantx = kt, \quad k \text{ is a constant}

(ii)

dxdt=[1t011t]x\frac {dx}{dt} = \begin{bmatrix} \frac 1t & 0 \\ 1 & \frac 1t \end{bmatrix} x

then

dx1dt=1tx1xi=kt\frac{dx_1}{dt} = \frac 1t x_1 \Rightarrow x_i = kt
dx2dt=x2+x2tkt2+mt\frac{dx_2}{dt} = x_2 + \frac{x_2}{t} \Rightarrow kt^2 + mt

thus

x(t)=[kt,kt2+mt]\boldsymbol{x}(t) = [kt, kt^2 + mt]^\top

(iii)

dxdt=[1t00011t00011t00001t]\frac{dx}{dt} = \begin{bmatrix} \frac 1t & 0 &0 & \cdots &0 \\ 1 & \frac 1t &0 &\cdots & 0 \\ 0 & 1 & \frac 1t & \cdots &0 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 0 & 0 &0 & \cdots &\frac 1t \end{bmatrix}
x1=c1t,and dxk+1dt=xk+xk+1tx_1 = c_1 t, \text{and } \frac{dx_{k+1}}{dt} = x_k + \frac {x_{k+1}}{t}

then we get

xk+1=u(t)tx_{k+1} = u(t)t

thus

xk+1=(xktdt+C)tx_{k+1} = (\int \frac{x_k}{t} dt + C)t
xk=i=1kci(i1)!tix_{k} = \sum_{i=1}^{k} \frac{c_i}{(i-1)!} t^i

Therefore

x(t)=[x1,x2,,xn]\boldsymbol{x}(t) = [x_1, x_2, \ldots, x_n]^\top