跳到主要内容

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

Author

Casablanca

Description

日本語版

a(t),b(t)a(t), b(t)tt のある有理式として次の実微分方程式を考える.

d2xdt+a(t)dxdt+b(t)x=0\frac{d^2 x}{dt} + a(t) \frac{dx}{dt} + b(t)x = 0

以下の問いに答えよ.

(i) k1k \geqq 1 をある整数として,x=tkx = t^k が式 (1) の解であるための a(t),b(t)a(t), b(t) に関する必要十分条件を求めよ.

以下では,ある整数 k1k \geqq 1 に対して (i) で求めた条件が成り立つものとし,ϕ(t)\phi(t)tkt^k と線形独立な解として,

p(t)=tdϕdt(t)kϕ(t)p(t) = t\frac{d \phi}{dt} (t) - k \phi(t)

とおく.

(ii) a(t),b(t)a(t), b(t)p(t)p(t) を用いて表わせ.

(iii) p(t)=tp(t) = t のとき a(t),b(t)a(t), b(t) を定めよ.

(iv) 式 (1) のすべての解が定数でない多項式のとき,a(t),b(t)a(t), b(t) は多項式でないことを示せ.

English Version

题目描述

a(t),b(t)a(t),b(t) 为关于 tt 的有理函数,考虑实微分方程

d2xdt2+a(t)dxdt+b(t)x=0.(1)\frac{d^2x}{dt^2}+a(t)\frac{dx}{dt}+b(t)x=0. \tag{1}

回答:

  1. 对给定整数 k1k\ge1,求 x=tkx=t^k 为 (1) 的解时 a(t),b(t)a(t),b(t) 所满足的充要条件。

以下假设该条件对某个 k1k\ge1 成立。令 ϕ(t)\phi(t) 为与 tkt^k 线性无关的另一解,并定义

p(t)=tdϕdt(t)kϕ(t).p(t)=t\frac{d\phi}{dt}(t)-k\phi(t).

继续回答:

  1. p(t)p(t) 表示 a(t),b(t)a(t),b(t)
  2. p(t)=tp(t)=t 时确定 a(t),b(t)a(t),b(t)
  3. 若 (1) 的每一个解都是非常值多项式,证明 a(t),b(t)a(t),b(t) 不是多项式。

考点

  • 二阶线性方程的多项式解条件:将 tkt^k 代入方程,得到两个有理系数之间的必要充分关系。
  • Wronskian 与第二解:由 tkt^kϕ(t)\phi(t) 的组合 p(t)p(t) 提取 Wronskian 信息,反求方程系数。
  • 多项式解的系数反证:分析两线性无关多项式解所决定的有理系数,证明它们不可能同时为多项式。

Kai

(i)

If k=1k = 1, then a(t)=tb(t)=0a(t) = tb(t) = 0.

If k2k \geq 2, then k(k1)+kta(t)+t2b(t)=0k(k-1) + kta(t) + t^2b(t) = 0.

Easy to see k(k1)+kta(t)+t2b(t)=0k(k-1) + kta(t) + t^2b(t) = 0 is neccessary and sufficient.

(ii)

Let Φ(t)=u(t)tk,u(t)≢Constant\Phi(t) = u(t)t^k, u(t) \not\equiv \text{Constant},

Φ(t)=ktk1u(t)+tku(t)\Phi '(t) = kt^{k-1}u(t) + t^ku'(t)
Φ(t)=tku(t)+2ktk1u(t)+k(k1)tk2u(t)\Phi''(t) = t^ku''(t) + 2kt^{k-1}u'(t) + k(k-1)t^{k-2}u(t)

then

p(t)=tk1u(t),p(t)=(k1)tk2u(t)+tk1u(t)p(t) = t^{k-1}u'(t), p'(t) = (k-1)t^{k-2}u'(t) + t^{k-1}u''(t)

And we can obtain:

tu(t)+(2k+a(t)t)u(t)=0t u''(t) + (2k + a(t)t)u'(t) = 0

Therefore

a(t)=u(t)u(t)2kt=(3k1)p(t)tp(t)a(t) = -\frac{u''(t)}{u'(t)} - \frac{2k}{t} = -\frac{(3k-1)p'(t)}{tp(t)}
b(t)=2k2t2kp(t)tp(t)b(t) = \frac{2k^2}{t^2} - \frac{kp'(t)}{tp(t)}

(iii)

a(t)=1t,b(t)=1t2a(t) = -\frac{1}{t}, b(t) = \frac{1}{t^2}

(iv)

Let x1x_1, x2x_2 be 2 independent particular solutions, then

x1=a(t)x1b(t)x1,x2=a(t)x2b(t)x2x_1'' = -a(t)x_1' - b(t)x_1 , x_2'' = -a(t)x_2' - b(t)x_2

and we have

(x1x2x1x2)=x1x2x1x2=a(t)x1x2+a(t)x1x2=a(t)(x1x2x1x2)\begin{aligned} (x_1'x_2 - x_1 x_2')' &= x_1''x_2 - x_1 x_2''\\ &= -a(t)x_1'x_2 + a(t)x_1x_2'\\ &= -a(t)(x_1'x_2 - x_1 x_2') \end{aligned}

so we get

x1x2x1x2=Cea(t)dtx_1'x_2 - x_1x_2' = C e^{\int a(t)dt}

Let A(t)A(t) denote a(t)dt\int a(t)dt. Since a(t)a(t) is a polynomial, and CeA(t)Ce^{-A(t)} is a polynomial for x1x2x1x2x_1'x_2 - x_1x_2' is a polynomial, we obtain a(t)0a(t) \equiv 0.

Then, consider

d2xdt2+b(t)x=0\frac{d^2x}{dt^2} + b(t)x = 0

If x1,x2x_1, x_2 are independent polynomials, and b(t)b(t) is polynomials, suppose that x(t)x(t) is mm times and b(t)b(t) is nn times, then b(t)0b(t) \equiv 0. Thus a(t),b(t)a(t), b(t) can't be polynomials.