跳到主要内容

京都大学 情報学研究科 数理工学専攻 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

题目描述

设整数 n>1n>1f:RnRnf:\mathbb R^n\to\mathbb R^nC1C^1 函数,考虑自治系统

dxdt=f(x),xRn.(1)\frac{d\boldsymbol x}{dt}=f(\boldsymbol x), \qquad \boldsymbol x\in\mathbb R^n. \tag{1}

x=ϕ(t)\boldsymbol x=\phi(t) 为定义在 R\mathbb R 上、有界且非常值的解。其沿轨道的变分方程为

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

其中 DfDfff 的 Jacobian 矩阵。回答:

  1. 若极限 a+=limt+ϕ(t)a_+=\lim_{t\to+\infty}\phi(t)a=limtϕ(t)a_-=\lim_{t\to-\infty}\phi(t) 存在,证明 x=a+\boldsymbol x=a_+x=a\boldsymbol x=a_- 都是 (1) 的常值解;并证明 (2) 存在一个在 R\mathbb R 上有界且满足 limt±ψ(t)=0\lim_{t\to\pm\infty}\psi(t)=0 的解 y=ψ(t)\boldsymbol y=\psi(t)
  2. 假设存在 C1C^1 函数 u:RnRnu:\mathbb R^n\to\mathbb R^n 满足
    Du(x)f(x)Df(x)u(x)=0.Du(\boldsymbol x)f(\boldsymbol x) -Df(\boldsymbol x)u(\boldsymbol x)=0.
    f(ϕ(0))f(\phi(0))u(ϕ(0))u(\phi(0)) 线性无关,求变分方程 (2) 的两个线性无关解。
  3. 假设存在 n1n-1C1C^1 函数 vj:RnRnv_j:\mathbb R^n\to\mathbb R^n 满足
    Dvj(x)f(x)Df(x)vj(x)=0(j=1,,n1),Dv_j(\boldsymbol x)f(\boldsymbol x) -Df(\boldsymbol x)v_j(\boldsymbol x)=0 \quad(j=1,\ldots,n-1),
    f(ϕ(0)),v1(ϕ(0)),,vn1(ϕ(0))f(\phi(0)),v_1(\phi(0)),\ldots,v_{n-1}(\phi(0)) 线性无关。求 (2) 的一般解。

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