跳到主要内容

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

Author

Casablanca, 祭音Myyura

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\to+\infty}\phi(t)=a_+ and continuity, we have

limt+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/Lq=L/\lVert L\rVert. For all sufficiently large tt,

qϕ(t)>L2.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=ax=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),limt±ψ(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,,n1)(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=1n1cjvj(ϕ(t)),c0,c1,,cn1R.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.