跳到主要内容

東京大学 新領域創成科学研究科 人間環境学専攻 2016年8月実施 微分方程

Author

思齐塾, 祭音Myyura

Description

(1) 直交座標 (x,y)(x, y) を極座標 (r,θ)(r, \theta) へ変換することによって, 微分方程式

[(x2+y2a)xy]dydx=(x2+y2a)y+x(1)\left[(x^2 + y^2 - a)x - y\right] \frac{dy}{dx} = (x^2 + y^2 - a)y + x \quad (1)

を解くことを考える。ただし、 aa は正定数とする。

  1. x=rcosθx = r\cos\thetay=rsinθy = r\sin\theta の関係を利用して, xdx+ydyxdx + ydy 及び xdyydxxdy - ydxrrθ\theta で表せ。

  2. 問1)の結果を用いて, 式(1)を極座標で表した微分方程式 drdθ=f(r,θ)\frac{dr}{d\theta} = f(r, \theta) に変換せよ。

  3. 問2)で得られた微分方程式の一般解 r(θ)r(\theta) を求め, 式(1)の一般解を求めよ。

题目描述

aa 为正常数。考虑通过把直角坐标 (x,y)(x,y) 变换为极坐标 (r,θ)(r,\theta),求解微分方程

[(x2+y2a)xy]dydx=(x2+y2a)y+x.(1)\left[(x^2+y^2-a)x-y\right]\frac{dy}{dx} =(x^2+y^2-a)y+x. \tag{1}

利用 x=rcosθ, y=rsinθx=r\cos\theta,\ y=r\sin\theta,回答下列问题:

  1. r,θr,\theta 表示微分形式 xdx+ydyx\,dx+y\,dyxdyydxx\,dy-y\,dx

  2. 利用第 1 问的结果,把式 (1) 化为极坐标下形如

    drdθ=f(r,θ)\frac{dr}{d\theta}=f(r,\theta)

    的微分方程。

  3. 求第 2 问所得方程的一般解 r(θ)r(\theta),进而给出式 (1) 的一般解。

Kai

(1)

x=rcosθ, y=rsinθx=r\cos\theta,\ y=r\sin\theta を微分すると

dx=cosθdrrsinθdθ,dy=sinθdr+rcosθdθ.dx=\cos\theta\,dr-r\sin\theta\,d\theta,\qquad dy=\sin\theta\,dr+r\cos\theta\,d\theta.

したがって

xdx+ydy=rdr,xdyydx=r2dθ.\boxed{x\,dx+y\,dy=r\,dr},\qquad \boxed{x\,dy-y\,dx=r^2\,d\theta}.

(2)

式 (1) に dxdx を掛けて一方の辺に移すと

(x2+y2a)(xdyydx)(xdx+ydy)=0.(x^2+y^2-a)(x\,dy-y\,dx)-(x\,dx+y\,dy)=0.

問1の結果と x2+y2=r2x^2+y^2=r^2 を用いれば

(r2a)r2dθrdr=0.(r^2-a)r^2\,d\theta-r\,dr=0.

よって r>0r>0 の範囲では

drdθ=r(r2a).\boxed{\frac{dr}{d\theta}=r(r^2-a)}.

(3)

r0r\ne0 かつ r2ar^2\ne a の区間で変数分離すると

drr(r2a)=dθ.\int\frac{dr}{r(r^2-a)}=\int d\theta.

左辺は

1alogr+12alogr2a-\frac1a\log r+\frac{1}{2a}\log|r^2-a|

であるから、定数をまとめて

logr2ar2=2aθ+C1.\log\left|\frac{r^2-a}{r^2}\right|=2a\theta+C_1.

符号も任意定数に含めれば

1ar2=Ce2aθ,r(θ)=a1Ce2aθ1-\frac{a}{r^2}=Ce^{2a\theta},\qquad \boxed{r(\theta)=\sqrt{\frac{a}{1-Ce^{2a\theta}}}}

となる(分母が正である各区間で考える)。 C=0C=0 は定常解 r=ar=\sqrt a も含む。

したがって、角度 θ=atan2(y,x)\theta=\operatorname{atan2}(y,x) の枝を一つ固定した局所表示では、式 (1) の一般解は

1ax2+y2=Cexp ⁣(2aatan2(y,x)).\boxed{1-\frac{a}{x^2+y^2}=C\exp\!\left(2a\,\operatorname{atan2}(y,x)\right)}.

極座標の r=0r=0 は原点一個を表し、通常の関数 y(x)y(x) としての積分曲線ではない。また原点では元の方程式の両係数がともに 0 となり、接線方向は方程式から定まらない。