跳到主要内容

京都大学 情報学研究科 数理工学専攻 2024年8月実施 常微分方程式

Author

思齐塾, 祭音Myyura

Description

r(t) を tt のある関数として次の微分方程式を考える.

d2xdt2=r(t)x(1)\frac{d^2x}{dt^2} = r(t)x \qquad (1)

ω(t)\omega(t)tt のある関数として、

ϕ(t)=exp(ω(t)dt)\phi(t) = \exp\left(\int \omega(t) dt\right)

としたとき, x=ϕ(t)x = \phi(t)r(t)r(t) の定義域上で式(1)の解になるものとする。以下の問いに答えよ。ただし、関数 r(t),ω(t),ϕ(t)r(t), \omega(t), \phi(t)C\mathbb{C} から高々有限個の点を除いた集合を定義域とし、定義域上で少なくとも連続であるとする.

(i) ω(t)\omega(t) が満たす1階微分方程式を求めよ.

(ii) ϕ(t)\phi(t)tt の有理関数となるとき, ω(t)\omega(t)r(t)r(t)tt の有理関数であることを示せ.

以下では、 ϕ(t)\phi(t)tt の有理関数であり, r(t)r(t) は恒等的には零ではないものと仮定する.

(iii) limtω(t)=0\lim_{t \to \infty} \omega(t) = 0 であることを示せ.

(iv) r1r_1r2r_2 を、それぞれ、 r(t)r(t) の分子と分母の次数とする。 r2>r1+1r_2 > r_1 + 1 であることを示せ.

题目描述

r(t)r(t)tt 的函数,考虑微分方程

d2xdt2=r(t)x.(1)\frac{d^2x}{dt^2}=r(t)x. \tag{1}

再设 ω(t)\omega(t)tt 的函数,并定义

ϕ(t)=exp ⁣(ω(t)dt).\phi(t) = \exp\!\left(\int\omega(t)\,dt\right).

假设 x=ϕ(t)x=\phi(t)r(t)r(t) 的定义域上是方程 (1) 的解。函数 r(t)r(t)ω(t)\omega(t)ϕ(t)\phi(t) 的定义域均为复平面 C\mathbb{C} 去掉至多有限个点后的集合,并且它们在定义域上至少连续。回答:

  1. ω(t)\omega(t) 所满足的一阶微分方程。
  2. 证明:若 ϕ(t)\phi(t)tt 的有理函数,则 ω(t)\omega(t)r(t)r(t) 也都是 tt 的有理函数。

以下再假设 ϕ(t)\phi(t)tt 的有理函数,并且 r(t)r(t) 不恒等于零。

  1. 证明
limtω(t)=0.\lim_{t\to\infty}\omega(t)=0.
  1. 分别以 r1r_1r2r_2 表示有理函数 r(t)r(t) 的分子次数和分母次数。证明
r2>r1+1.r_2>r_1+1.

Kai

(i) x=ϕ(t)x = \phi(t)d2xdt2=r(t)x\frac{d^2x}{dt^2} = r(t)x の解であるから、 ϕ(t)=r(t)ϕ(t)\phi''(t) = r(t)\phi(t) が成り立つ。

ϕ(t)=exp(ω(t)dt)\phi(t) = \exp\left(\int \omega(t) dt\right) より、

ϕ(t)=ω(t)ϕ(t)\phi'(t) = \omega(t) \phi(t)

ϕ(t)=ω(t)ϕ(t)+ω(t)ϕ(t)=ω(t)ϕ(t)+ω(t)2ϕ(t)=(ω(t)+ω(t)2)ϕ(t)\phi''(t) = \omega'(t) \phi(t) + \omega(t) \phi'(t) = \omega'(t) \phi(t) + \omega(t)^2 \phi(t) = (\omega'(t) + \omega(t)^2)\phi(t)

したがって、 ω(t)+ω(t)2=r(t)\omega'(t) + \omega(t)^2 = r(t) が成り立つ。

(ii) ϕ(t)\phi(t)tt の有理関数であるとする。このとき、 ϕ(t)=exp(ω(t)dt)\phi(t) = \exp\left(\int \omega(t) dt\right) より、 ln(ϕ(t))=ω(t)dt\ln(\phi(t)) = \int \omega(t) dt が成り立つ。 ϕ(t)\phi(t) が有理関数であるから、 ln(ϕ(t))\ln(\phi(t)) の微分も有理関数である。したがって、 ω(t)=ddtln(ϕ(t))=ϕ(t)ϕ(t)\omega(t) = \frac{d}{dt} \ln(\phi(t)) = \frac{\phi'(t)}{\phi(t)}tt の有理関数である。

ω(t)=r(t)ω(t)2\omega'(t) = r(t) - \omega(t)^2 より、 r(t)=ω(t)+ω(t)2r(t) = \omega'(t) + \omega(t)^2 である。 ω(t)\omega(t) が有理関数であるから、 ω(t)\omega'(t) も有理関数であり、 ω(t)2\omega(t)^2 も有理関数である。したがって、 r(t)r(t)tt の有理関数である。

(iii) ϕ(t)=P(t)/Q(t)\phi(t)=P(t)/Q(t) とし、 P,QP,Q を零でない多項式とする。 m=degPdegQm=\deg P-\deg Q とおけば、ある零でない定数 cc に対して

ϕ(t)=ctm(1+O(t1))(t)\phi(t)=ct^m\left(1+O(t^{-1})\right)\qquad(t\to\infty)

と書ける。(ii) より ω=ϕ/ϕ\omega=\phi'/\phi なので、この式を対数微分して

ω(t)=mt+O(t2)\omega(t)=\frac{m}{t}+O(t^{-2})

を得る。特に m=0m=0 のときは ω(t)=O(t2)\omega(t)=O(t^{-2}) である。いずれの場合も

limtω(t)=0\boxed{\lim_{t\to\infty}\omega(t)=0}

が成り立つ。

(iv) 上の漸近式から

ω(t)=O(t2),ω(t)2=O(t2)\omega'(t)=O(t^{-2}),\qquad \omega(t)^2=O(t^{-2})

である。したがって、(i) の関係式より

r(t)=ω(t)+ω(t)2=O(t2)(t)r(t)=\omega'(t)+\omega(t)^2=O(t^{-2})\qquad(t\to\infty)

となる。 rr は恒等的には零でないと仮定されているので、既約な形で r(t)=R(t)/S(t)r(t)=R(t)/S(t) と書き、 r1=degRr_1=\deg Rr2=degSr_2=\deg S とすれば、この評価から

r2r12r_2-r_1\ge 2

でなければならない。次数は整数なので、これは

r2>r1+1\boxed{r_2>r_1+1}

と同値である。