跳到主要内容

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

Author

Casablanca

Description

日本語版

n>1n > 1 を整数,f:RnRnf : \mathbb{R}^n \rightarrow \mathbb{R}^nC1C^1 級関数として,R\mathbb{R} 上の微分方程式系

dxdt=f(x),xRn\begin{align} \frac{dx}{dt} = f(x), \quad x \in \mathbb{R}^n \tag{1} \end{align}

を考える。x=ϕ(t)x = \phi(t)R\mathbb{R} 上で有界な式 (1) の非定数解とする。 次式を式 (1) の解 x=ϕ(t)x = \phi(t) のまわりの変分方程式という:

dydt=Df(ϕ(t))y,yRn(2)\frac{dy}{dt} = Df(\phi(t))y, \quad y \in \mathbb{R}^n \tag{2}

ここで,Df(x)Df(x)f(x)f(x) のヤコビ行列で,各 j=1,2,,nj = 1, 2, \ldots, n に対して,fj(x)f_j(x)xjx_j をそれぞれ,f(x)f(x)xx の第 jj 成分として,

Df(x)=(f1x1(x)f1xn(x)fnx1(x)fnxn(x))Df(x) = \begin{pmatrix} \frac{\partial f_1}{\partial x_1}(x) & \dots & \frac{\partial f_1}{\partial x_n}(x) \\ \vdots & \ddots & \vdots \\ \frac{\partial f_n}{\partial x_1}(x) & \dots & \frac{\partial f_n}{\partial x_n}(x) \end{pmatrix}

で与えられる nn 次正方行列である。以下の問いに答えよ。

(i) 極限 a+=limt+ϕ(t)a_+ = \lim_{t \to +\infty} \phi(t)a=limtϕ(t)a_- = \lim_{t \to -\infty} \phi(t) が存在するとき,x=a+x = a_+aa_- が式 (1) の定数解であることを示せ。また,変分方程式 (2) が limt±ψ(t)=0\lim_{t \to \pm \infty} \psi(t) = 0 かつ R\mathbb{R} 上で有界な解 y=ψ(t)y = \psi(t) をもつことを示せ。

(ii) 次式を満たす C1C^1 級関数 u:RnRnu : \mathbb{R}^n \rightarrow \mathbb{R}^n が存在するものとする。

Du(x)f(x)Df(x)u(x)=0Du(x)f(x) - Df(x)u(x) = 0

2 個のベクトル f(ϕ(0))f(\phi(0))u(ϕ(0))u(\phi(0)) が線形独立であるとき,変分方程式 (2) の線形独立な解を2個求めよ。

(iii) 次式を満たす n1n - 1 個の C1C^1 級関数 vj:RnRn (j=1,2,,n1)v_j : \mathbb{R}^n \rightarrow \mathbb{R}^n \ (j = 1, 2, \ldots, n - 1) が存在するものとする。

Dvj(x)f(x)Df(x)vj(x)=0(j=1,2,,n1)Dv_j(x)f(x) - Df(x)v_j(x) = 0 \quad (j = 1, 2, \ldots, n - 1)

nn 個のベクトル f(ϕ(0))f(\phi(0))vj(ϕ(0)) (j=1,2,,n1)v_j(\phi(0)) \ (j = 1, 2, \ldots, n - 1) が線形独立であるとき,変分方程式 (2) の一般解を求めよ。

English Version

Kai

(i)

Since ϕ(t)\phi (t) is a solution, we have

dϕ(t)dt=f(ϕ(t))\frac{d \phi(t)}{dt} = f(\phi(t))

From limt+ϕ(t)=a+\lim_{t\rightarrow +\infty}\phi(t) = a_+ and continuity, we have

limt+dϕ(t)dt=f(a+)\lim_{t \rightarrow +\infty} \frac{d\phi(t)}{dt} = f(a_+)

suppose that

f(a+)>0f(a_+) > 0

then we have

limtdϕ(t)dt>0limtϕ(t)=+\lim_{t \rightarrow \infty} \frac{d\phi(t)}{dt} > 0 \Rightarrow \lim_{t\rightarrow \infty} \phi(t) = +\infty

which is conflict with

limtϕ(t)=a+\lim_{t\rightarrow \infty} \phi(t) = a_+

and therefore f(a+)=0f(a_+) = 0. Thus ϕ(t)=a+\phi(t) = a_+ is a constant solution. Similarly, ϕ(t)=a\phi(t) = a_{-} is a constant solution,too.

Notice that

df(ϕ(t))dt=Df(ϕ(t))dϕ(t)dt,f(ϕ(t))=dϕ(t)dt\frac{df(\phi(t))}{dt} = Df(\phi(t))\frac{d\phi(t)}{dt}, f(\phi(t)) = \frac{d\phi(t)}{dt}

hence

ddϕ(t)dtdt=Df(ϕ(t))dϕ(t)dt,ψ(t)=dϕ(t)dt=f(ϕ(t))\frac{d\frac{d\phi(t)}{dt}}{dt} = Df(\phi(t))\frac{d\phi(t)}{dt}, \psi(t) = \frac{d\phi(t)}{dt} = f(\phi(t))

and we konw

limtψ(t)=f(a+)=0\lim_{t\rightarrow \infty} \psi(t) = f(a_+) = 0

since fC1f \in C^1, ϕ(t)\phi(t) is bounded, then ψ(t)\psi(t) is bounded

(ii)

f(ϕ(t))f(\phi(t)) is a solution, then

v(ϕ(t))dt=Dv(ϕ(t))dϕ(t)dt=Dv(x)f(x)=Df(x)v(x)=Df(ϕ(t))v(ϕ(t))\frac{v(\phi(t))}{dt} = Dv(\phi(t))\frac{d\phi(t)}{dt} = Dv(x)f(x) = Df(x)v(x) = Df(\phi(t))v(\phi(t))

and we see that v(ϕ(t))v(\phi(t)) is a solution to (2) and we have

v(ϕ(t))=v(ϕ(0))+Df(ϕ(t))v(ϕ(t))dtv(\phi(t)) = v(\phi(0)) + \int Df(\phi(t))v(\phi(t)) dt
f(ϕ(t))=f(ϕ(0))+Df(ϕ(t))v(ϕ(t))dtf(\phi(t)) = f(\phi(0)) + \int Df(\phi(t))v(\phi(t)) dt

thus, v(ϕ(0))v(\phi(0)), f(ϕ(0))f(\phi(0)) are independent \Rightarrow v(ϕ(t))v(\phi(t)), f(ϕ(t))f(\phi(t)) are independent.

(iii)

Similar to (ii)\boldsymbol{(\text{ii})}, omitted