跳到主要内容

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

Author​

Casablanca, 祭音Myyura

Description​

日本語版​

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

dxdt=f(x),x∈Rn\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,y∈Rn(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)=(∂f1∂x1(x)…∂f1∂xn(x)⋮⋱⋮∂fn∂x1(x)…∂fn∂xn(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+=lim⁡t→+∞ϕ(t)a_+ = \lim_{t \to +\infty} \phi(t) と a−=lim⁡t→−∞ϕ(t)a_- = \lim_{t \to -\infty} \phi(t) が存在するとき,x=a+x = a_+ と a−a_- が式 (1) の定数解であることを示せ。また,変分方程式 (2) が lim⁡t→±∞ψ(t)=0\lim_{t \to \pm \infty} \psi(t) = 0 かつ R\mathbb{R} 上で有界な解 y=ψ(t)y = \psi(t) をもつことを示せ。

(ii) 次式を満たす C1C^1 級関数 u:Rn→Rnu : \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) 次式を満たす n−1n - 1 個の C1C^1 級関数 vj:Rn→Rn (j=1,2,…,n−1)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,…,n−1)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,…,n−1)v_j(\phi(0)) \ (j = 1, 2, \ldots, n - 1) が線形独立であるとき,変分方程式 (2) の一般解を求めよ。

English Version​

题目描述​

设整数 n>1n>1,f:Rn→Rnf:\mathbb R^n\to\mathbb R^n 为 C1C^1 函数,考虑自治系统

dxdt=f(x),x∈Rn.(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,y∈Rn,(2)\frac{d\boldsymbol y}{dt}=Df(\phi(t))\boldsymbol y, \qquad \boldsymbol y\in\mathbb R^n, \tag{2}

其中 DfDf 是 ff 的 Jacobian 矩阵。回答:

  1. 若极限 a+=lim⁡t→+∞ϕ(t)a_+=\lim_{t\to+\infty}\phi(t) 与 a−=lim⁡t→−∞ϕ(t)a_-=\lim_{t\to-\infty}\phi(t) 存在,证明 x=a+\boldsymbol x=a_+ 和 x=a−\boldsymbol x=a_- 都是 (1) 的常值解;并证明 (2) 存在一个在 R\mathbb R 上有界且满足 lim⁡t→±∞ψ(t)=0\lim_{t\to\pm\infty}\psi(t)=0 的解 y=ψ(t)\boldsymbol y=\psi(t)。

  2. 假设存在 C1C^1 函数 u:Rn→Rnu:\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. 假设存在 n−1n-1 个 C1C^1 函数 vj:Rn→Rnv_j:\mathbb R^n\to\mathbb R^n 满足

    Dvj(x)f(x)−Df(x)vj(x)=0(j=1,…,n−1),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)),…,vn−1(ϕ(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 lim⁡t→+∞ϕ(t)=a+\lim_{t\to+\infty}\phi(t)=a_+ and continuity, we have

lim⁡t→+∞dϕ(t)dt=f(a+).\lim_{t\to+\infty}\frac{d\phi(t)}{dt}=f(a_+).

If L=f(a+)≠0L=f(a_+)\ne0, put q=L/∥L∥q=L/\lVert L\rVert. For all sufficiently large tt,

q⊤ϕ′(t)>∥L∥2.q^\top\phi'(t)>\frac{\lVert L\rVert}{2}.

Integrating this inequality contradicts the convergence of q⊤ϕ(t)q^\top\phi(t). Hence f(a+)=0f(a_+)=0, so x=a+x=a_+ is a constant solution. The same argument at −∞-\infty gives f(a−)=0f(a_-)=0, so x=a−x=a_- is also a constant solution.

Notice that

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

Set

ψ(t)=ϕ′(t)=f(ϕ(t)).\psi(t)=\phi'(t)=f(\phi(t)).

Then

ψ′(t)=Df(ϕ(t))ψ(t),lim⁡t→±∞ψ(t)=f(a±)=0.\psi'(t)=Df(\phi(t))\psi(t),\qquad \lim_{t\to\pm\infty}\psi(t)=f(a_\pm)=0.

Since ϕ(R)\phi(\mathbb R) is bounded, its closure is compact; continuity of ff therefore implies that ψ(t)=f(ϕ(t))\psi(t)=f(\phi(t)) is bounded on R\mathbb R. Since ϕ\phi is nonconstant, ψ\psi is not identically zero.

(ii)​

The two solutions are

y1(t)=f(ϕ(t)),y2(t)=u(ϕ(t)).y_1(t)=f(\phi(t)),\qquad y_2(t)=u(\phi(t)).

Indeed,

y1′=Df(ϕ)y1,y2′=Du(ϕ)f(ϕ)=Df(ϕ)u(ϕ)=Df(ϕ)y2.y_1'=Df(\phi)y_1,\qquad y_2'=Du(\phi)f(\phi)=Df(\phi)u(\phi)=Df(\phi)y_2.

They are linearly independent because their values at t=0t=0 are linearly independent.

(iii)​

By the same calculation, f(ϕ(t))f(\phi(t)) and vj(ϕ(t))v_j(\phi(t)) (j=1,…,n−1)(j=1,\ldots,n-1) are solutions. Their values at t=0t=0 form a basis of Rn\mathbb R^n, hence the general solution is

y(t)=c0f(ϕ(t))+∑j=1n−1cjvj(ϕ(t)),c0,c1,…,cn−1∈R.y(t)=c_0f(\phi(t))+\sum_{j=1}^{n-1}c_jv_j(\phi(t)), \qquad c_0,c_1,\ldots,c_{n-1}\in\mathbb R.