§E13.19Liénard 系と極限閉軌道

最終更新

van der Pol 方程式は、外力を周期的に加えなくても振動が持続する自己励振の基本的な例である。原点近傍でのエネルギーの増加と遠方での散逸が、どのような周期運動につながるかを調べるために、平面上の Liénard 系を考える。本記事では、Liénard の条件のもとで周期軌道がただ一つ存在し、原点以外のすべての点から出発する解が、長時間の後に同じ周期軌道へ近づくことを扱う。

1 Liénard 系と半周の軌道

定義 1.1.Φ,g:R→R\Phi,g:\R\to\RをC1C^1級関数とする。F:R2→R2F:\R^2\to\R^2、F(x,y)=(y−Φ(x), −g(x))F(x,y)=(y-\Phi(x),\,-g(x))を(Φ,g)(\Phi,g)のLiénard 系 (Liénard system) という。(Φ,g)(\Phi,g)が次の条件をすべて満たすとき、(Φ,g)(\Phi,g)はLiénard の条件 (Liénard conditions) を満たすという。

  1. Φ\Phiとggは奇関数である。
  2. x≠0x\ne0ならばxg(x)>0xg(x)>0である。
  3. Φ′(0)<0\Phi'(0)<0である。
  4. Φ\Phiは(0,∞)(0,\infty)にただ一つの零点aaをもち、[a,∞)[a,\infty)上で狭義単調増加であり、x→∞x\to\inftyのときΦ(x)→∞\Phi(x)\to\inftyである。

補題 1.2.(Φ,g)(\Phi,g)を Liénard の条件を満たす組、aaをΦ\Phiの正の零点、FFを(Φ,g)(\Phi,g)の Liénard 系、φ\varphiをFFの流れとする。z∈R2z\in\R^2とt∈I(z)t\in I(z)についてφt(z)=(x(t),y(t))\varphi_t(z)=(x(t),y(t))と書く。このとき、次が成り立つ。

  1. 0<x<a0<x<aならばΦ(x)<0\Phi(x)<0、x>ax>aならばΦ(x)>0\Phi(x)>0である。FFの平衡点は原点だけである。z∈R2z\in\R^2についてI(−z)=I(z)I(-z)=I(z)であり、t∈I(z)t\in I(z)についてφt(−z)=−φt(z)\varphi_t(-z)=-\varphi_t(z)である。
  2. z=(x0,y0)z=(x_0,y_0)がx0>0x_0>0またはx0=0<y0x_0=0<y_0を満たすとする。このとき、[0,σ]⊆I(z)[0,\sigma]\subseteq I(z)を満たし、0<t<σ0<t<\sigmaならばx(t)>0x(t)>0であり、x(σ)=0x(\sigma)=0、y(σ)<0y(\sigma)<0を満たすσ>0\sigma>0がただ一つ存在する。このσ\sigmaをσ(z)\sigma(z)と書く。yyは[0,σ(z)][0,\sigma(z)]上で狭義単調減少である。
  3. α>0\alpha>0についてz=(0,α)z=(0,\alpha)とし、σ(α)=σ(z)\sigma(\alpha)=\sigma(z)、β(α)=−y(σ(α))>0\beta(\alpha)=-y(\sigma(\alpha))>0、A(α)={φt(z):0≤t≤σ(α)}A(\alpha)=\{\varphi_t(z):0\le t\le\sigma(\alpha)\}とおく。このとき、0≤t<tm0\le t<t_mならばx′(t)>0x'(t)>0、tm<t≤σ(α)t_m<t\le\sigma(\alpha)ならばx′(t)<0x'(t)<0を満たすtm∈(0,σ(α))t_m\in(0,\sigma(\alpha))がただ一つ存在する。xmax⁡(α)=x(tm)x_{\max}(\alpha)=x(t_m)とおく。η∈[−β(α),α]\eta\in[-\beta(\alpha),\alpha]についてy(t)=ηy(t)=\etaを満たすt∈[0,σ(α)]t\in[0,\sigma(\alpha)]はただ一つであり、そのttでのx(t)x(t)をXα(η)X_\alpha(\eta)と書く。ξ∈[0,xmax⁡(α)]\xi\in[0,x_{\max}(\alpha)]についてx(t)=ξx(t)=\xiを満たすt∈[0,tm]t\in[0,t_m]とt∈[tm,σ(α)]t\in[t_m,\sigma(\alpha)]はそれぞれただ一つであり、それらのttでのy(t)y(t)をYα+(ξ)Y^+_\alpha(\xi)、Yα−(ξ)Y^-_\alpha(\xi)と書く。XαX_\alphaとYα±Y^\pm_\alphaは連続であり、Yα+Y^+_\alphaは狭義単調減少、Yα−Y^-_\alphaは狭義単調増加で、0≤ξ<xmax⁡(α)0\le\xi<x_{\max}(\alpha)ならばYα−(ξ)<Φ(ξ)<Yα+(ξ)Y^-_\alpha(\xi)<\Phi(\xi)<Y^+_\alpha(\xi)である。
  4. σ\sigmaとβ\betaは(0,∞)(0,\infty)上の連続関数である。

証明.主張 (1)を示す。Φ\Phiは奇関数であるからΦ(0)=0\Phi(0)=0であり、Φ′(0)<0\Phi'(0)<0から十分小さいx>0x>0でΦ(x)<0\Phi(x)<0である。Φ\Phiは(0,a)(0,a)に零点をもたないから、§E2.43 系 2.3により(0,a)(0,a)上でΦ<0\Phi<0である。Φ\Phiは[a,∞)[a,\infty)上で狭義単調増加であるから(a,∞)(a,\infty)上でΦ>0\Phi>0である。F(x,y)=0F(x,y)=0ならばg(x)=0g(x)=0であり、xg(x)>0xg(x)>0(x≠0x\ne0)からx=0x=0、y=Φ(0)=0y=\Phi(0)=0である。Φ,g\Phi,gは奇関数であるからF(−w)=−F(w)F(-w)=-F(w)である。t↦−φt(z)t\mapsto-\varphi_t(z)はI(z)I(z)上で(−φt(z))′=−F(φt(z))=F(−φt(z))(-\varphi_t(z))'=-F(\varphi_t(z))=F(-\varphi_t(z))を満たす古典解であり、t=0t=0で−z-zをとるから、§E13.15 補題 1.2によりI(z)⊆I(−z)I(z)\subseteq I(-z)であり、φt(−z)=−φt(z)\varphi_t(-z)=-\varphi_t(z)である。zzを−z-zに取り替えてI(−z)⊆I(z)I(-z)\subseteq I(z)を得る。

主張 (2)を示す。x0>0x_0>0ならば連続性により、x0=0<y0x_0=0<y_0ならばx′(0)=y0−Φ(0)=y0>0x'(0)=y_0-\Phi(0)=y_0>0により、十分小さいt>0t>0でx(t)>0x(t)>0である。T∗=sup⁡{t∈I(z):t>0, 0<s≤t ならば x(s)>0}T^*=\sup\{t\in I(z):t>0,\ 0<s\le t\text{ ならば }x(s)>0\}とおき、β∗=sup⁡I(z)\beta^*=\sup I(z)とおく。0<t<T∗0<t<T^*でx(t)>0x(t)>0であるからy′(t)=−g(x(t))<0y'(t)=-g(x(t))<0であり、yyは[0,T∗)[0,T^*)上で狭義単調減少である。X>max⁡{x0,a}X>\max\{x_0,a\}をΦ(X)>y0\Phi(X)>y_0を満たすようにとる。0≤t<T∗0\le t<T^*ならばx(t)<Xx(t)<Xである。実際、x(t1)≥Xx(t_1)\ge Xを満たすt1∈(0,T∗)t_1\in(0,T^*)があるとすると、t0=min⁡{t∈[0,t1]:x(t)≥X}t_0=\min\{t\in[0,t_1]:x(t)\ge X\}についてt0>0t_0>0、x(t0)=Xx(t_0)=X、[0,t0)[0,t_0)上でx<Xx<Xであるからx′(t0)≥0x'(t_0)\ge0であり、一方x′(t0)=y(t0)−Φ(X)<y0−Φ(X)<0x'(t_0)=y(t_0)-\Phi(X)<y_0-\Phi(X)<0である。M=max⁡[0,X]gM=\max_{[0,X]}gとおくと、0≤t<T∗0\le t<T^*でy(t)≥y0−Mty(t)\ge y_0-Mtである。

主張 1.2.1.T∗<β∗T^*<\beta^*である。

証明.T∗=β∗T^*=\beta^*と仮定する。β∗<∞\beta^*<\inftyならば、t∈I(z)t\in I(z)、t≥0t\ge0についてφt(z)\varphi_t(z)はコンパクト集合[0,X]×[y0−Mβ∗,y0][0,X]\times[y_0-M\beta^*,y_0]に属するから、§E13.15 補題 3.1 (1)により[0,∞)⊆I(z)[0,\infty)\subseteq I(z)であり、β∗<∞\beta^*<\inftyに反する。したがってβ∗=∞\beta^*=\inftyであり、t>0t>0でx(t)>0x(t)>0、yyは狭義単調減少である。y∞=lim⁡t→∞y(t)∈[−∞,y0)y_\infty=\lim_{t\to\infty}y(t)\in[-\infty,y_0)とおく。y∞=−∞y_\infty=-\inftyならば、P=max⁡[0,X]∣Φ∣P=\max_{[0,X]}\lvert\Phi\rvertとおくと、十分大きいt2t_2についてt≥t2t\ge t_2でy(t)≤−P−1y(t)\le-P-1、x′(t)≤−1x'(t)\le-1であり、x(t)≤x(t2)−(t−t2)x(t)\le x(t_2)-(t-t_2)はt→∞t\to\inftyで負になってx>0x>0に反する。y∞y_\inftyが有限ならば、前方軌道はコンパクト集合K=[0,X]×[y∞,y0]K=[0,X]\times[y_\infty,y_0]に含まれる。§E13.15 定理 5.3 (2)によりω(z)\omega(z)は空でなくKKに含まれ、w∈ω(z)w\in\omega(z)についてI(w)=RI(w)=\Rであり、§E13.15 定理 5.3 (1)によりγ(w)⊆ω(z)\gamma(w)\subseteq\omega(z)である。ω(z)\omega(z)の点の第二座標はy∞y_\inftyであるから、φt(w)\varphi_t(w)の第二座標はttによらずy∞y_\inftyであり、その微分−g-gが零になることからφt(w)\varphi_t(w)の第一座標は常に00である。第一座標の微分y∞−Φ(0)=y∞y_\infty-\Phi(0)=y_\inftyも零であるから、w=0w=0であり、ω(z)={0}\omega(z)=\{0\}、y∞=0y_\infty=0である。§E13.15 定理 5.3 (2)によりφt(z)→0\varphi_t(z)\to0である。yyは狭義単調減少で極限00をもつから、すべてのt≥0t\ge0でy(t)>0y(t)>0である。t3t_3をt≥t3t\ge t_3で0<x(t)<a0<x(t)<aとなるようにとると、t≥t3t\ge t_3でx′(t)=y(t)−Φ(x(t))>0x'(t)=y(t)-\Phi(x(t))>0であり、x(t)≥x(t3)>0x(t)\ge x(t_3)>0となってx(t)→0x(t)\to0に反する。▨

主張 1.2.1によりT∗∈I(z)T^*\in I(z)である。xxは(0,T∗)(0,T^*)上で正であるからx(T∗)≥0x(T^*)\ge0であり、x(T∗)>0x(T^*)>0ならばT∗T^*の少し先までx>0x>0となってT∗T^*の定義に反するから、x(T∗)=0x(T^*)=0である。(0,T∗)(0,T^*)上でx>0x>0であるからx′(T∗)≤0x'(T^*)\le0であり、x′(T∗)=y(T∗)−Φ(0)=y(T∗)x'(T^*)=y(T^*)-\Phi(0)=y(T^*)である。y(T∗)=0y(T^*)=0ならばφT∗(z)=0\varphi_{T^*}(z)=0は平衡点であり、§E13.15 命題 4.3 (1)と§E13.15 命題 4.3 (2)によりz=0z=0となって仮定に反するから、y(T∗)<0y(T^*)<0である。σ=T∗\sigma=T^*は条件を満たし、条件を満たすσ\sigmaはxxの最初の正の零点であるからT∗T^*に等しい。(0,T∗)(0,T^*)上でy′<0y'<0であるから、yyは[0,T∗][0,T^*]上で狭義単調減少である。

主張 (3)を示す。σ=σ(α)\sigma=\sigma(\alpha)とし、k(t)=x′(t)=y(t)−Φ(x(t))k(t)=x'(t)=y(t)-\Phi(x(t))とおく。k(0)=α>0k(0)=\alpha>0、k(σ)=y(σ)<0k(\sigma)=y(\sigma)<0であり、k(t′)=0k(t')=0、t′∈(0,σ)t'\in(0,\sigma)ならばk′(t′)=y′(t′)−Φ′(x(t′))x′(t′)=−g(x(t′))<0k'(t')=y'(t')-\Phi'(x(t'))x'(t')=-g(x(t'))<0である。tm=min⁡{t∈[0,σ]:k(t)=0}t_m=\min\{t\in[0,\sigma]:k(t)=0\}とおくと、[0,tm)[0,t_m)上でk>0k>0であり、k′(tm)<0k'(t_m)<0からtmt_mの直後でk<0k<0である。(tm,σ](t_m,\sigma]にk(t)≥0k(t)\ge0を満たす点があるとし、その下限をt∗t^*とすると、tm<t∗<σt_m<t^*<\sigma、k(t∗)=0k(t^*)=0で(tm,t∗)(t_m,t^*)上でk<0k<0であるからk′(t∗)≥0k'(t^*)\ge0となり、k′(t∗)<0k'(t^*)<0に反する。したがって(tm,σ](t_m,\sigma]上でk<0k<0である。y∣[0,σ]y|_{[0,\sigma]}は[0,σ][0,\sigma]から[−β(α),α][-\beta(\alpha),\alpha]への狭義単調減少な連続全単射であり、x∣[0,tm]x|_{[0,t_m]}、x∣[tm,σ]x|_{[t_m,\sigma]}はそれぞれ[0,xmax⁡(α)][0,x_{\max}(\alpha)]への狭義単調な連続全単射である。§E2.19 定理 6.8によりこれらの逆写像は連続であり、Xα=x∘(y∣[0,σ])−1X_\alpha=x\circ(y|_{[0,\sigma]})^{-1}、Yα+=y∘(x∣[0,tm])−1Y^+_\alpha=y\circ(x|_{[0,t_m]})^{-1}、Yα−=y∘(x∣[tm,σ])−1Y^-_\alpha=y\circ(x|_{[t_m,\sigma]})^{-1}は連続である。x∣[0,tm]x|_{[0,t_m]}の逆写像は狭義単調増加、x∣[tm,σ]x|_{[t_m,\sigma]}の逆写像は狭義単調減少でyyは狭義単調減少であるから、Yα+Y^+_\alphaは狭義単調減少、Yα−Y^-_\alphaは狭義単調増加である。0≤ξ<xmax⁡(α)0\le\xi<x_{\max}(\alpha)ならば、x(t)=ξx(t)=\xiとなるttは[0,tm)[0,t_m)と(tm,σ](t_m,\sigma]にあり、[0,tm)[0,t_m)上でk>0k>0、(tm,σ](t_m,\sigma]上でk<0k<0であるからYα−(ξ)<Φ(ξ)<Yα+(ξ)Y^-_\alpha(\xi)<\Phi(\xi)<Y^+_\alpha(\xi)である。

主張 (4)を示す。α0>0\alpha_0>0、z0=(0,α0)z_0=(0,\alpha_0)、σ0=σ(α0)\sigma_0=\sigma(\alpha_0)、β0=β(α0)\beta_0=\beta(\alpha_0)とし、α>0\alpha>0についてzα=(0,α)z_\alpha=(0,\alpha)とおく。e=(0,1)e=(0,1)とし、Σ0={(0,−β0)+se:∣s∣≤β0/2}\Sigma_0=\{(0,-\beta_0)+se:\lvert s\rvert\le\beta_0/2\}とおくと、(0,y)∈Σ0(0,y)\in\Sigma_0についてy<0y<0、F(0,y)=(y,0)F(0,y)=(y,0)、det⁡(e,F(0,y))=−y≠0\det(e,F(0,y))=-y\ne0であるから、Σ0\Sigma_0は横断線分である。§E13.18 補題 1.4をz0z_0、t0=σ0t_0=\sigma_0、Σ0\Sigma_0に適用してV,η,τV,\eta,\tauを得る。0<η′<min⁡{η,σ0/2}0<\eta'<\min\{\eta,\sigma_0/2\}についてVVをτ−1((σ0−η′,σ0+η′))\tau^{-1}((\sigma_0-\eta',\sigma_0+\eta'))に縮めると§E13.18 補題 1.4 (1)はη′\eta'について成り立つから、η<σ0/2\eta<\sigma_0/2としてよい。§E13.18 補題 1.4 (1)により、z∈Vz\in V、t∈(σ0−η,σ0+η)t\in(\sigma_0-\eta,\sigma_0+\eta)について、φt(z)\varphi_t(z)の第一座標が00であることは、t=τ(z)t=\tau(z)であるための必要十分条件である。{(t,z)∈D:(φt(z))2−Φ((φt(z))1)>0}\{(t,z)\in D:(\varphi_t(z))_2-\Phi((\varphi_t(z))_1)>0\}は(0,z0)(0,z_0)を含む開集合であるから、η1∈(0,σ0−η)\eta_1\in(0,\sigma_0-\eta)とz0z_0の開近傍V1V_1を、[0,η1]×V1[0,\eta_1]\times V_1がこの集合に含まれるようにとることができる。z∈V1z\in V_1についてx(t)=x(0)+∫0tx′(s) dsx(t)=x(0)+\int_0^tx'(s)\,dsであるから、z=zα∈V1z=z_\alpha\in V_1ならば0<t≤η10<t\le\eta_1でx(t)>0x(t)>0である。m=min⁡{x(t):t∈[η1,σ0−η]}m=\min\{x(t):t\in[\eta_1,\sigma_0-\eta]\}(z=z0z=z_0)とおくとm>0m>0である。R2×[η1,σ0−η]\R^2\times[\eta_1,\sigma_0-\eta]の開集合{(z,t):(t,z)∈D, ∥φt(z)−φt(z0)∥<m}\{(z,t):(t,z)\in D,\ \lVert\varphi_t(z)-\varphi_t(z_0)\rVert<m\}は{z0}×[η1,σ0−η]\{z_0\}\times[\eta_1,\sigma_0-\eta]を含むから、§E2.19 補題 5.1によりz0z_0の開近傍V2V_2でz∈V2z\in V_2、t∈[η1,σ0−η]t\in[\eta_1,\sigma_0-\eta]ならば∥φt(z)−φt(z0)∥<m\lVert\varphi_t(z)-\varphi_t(z_0)\rVert<mを満たすものが存在する。zα∈V∩V1∩V2z_\alpha\in V\cap V_1\cap V_2ならば、φt(zα)\varphi_t(z_\alpha)の第一座標は0<t≤σ0−η0<t\le\sigma_0-\etaで正であり、(σ0−η,σ0+η)(\sigma_0-\eta,\sigma_0+\eta)ではt=τ(zα)t=\tau(z_\alpha)でだけ00になる。したがってσ(α)=τ(zα)\sigma(\alpha)=\tau(z_\alpha)であり、β(α)\beta(\alpha)はφτ(zα)(zα)\varphi_{\tau(z_\alpha)}(z_\alpha)の第二座標の−1-1倍である。α↦zα\alpha\mapsto z_\alphaとτ\tauとφ\varphiは連続であるから、σ\sigmaとβ\betaはα0\alpha_0で連続である。▨

2 半周の軌道の比較とエネルギー差

補題 2.1.(Φ,g)(\Phi,g)を Liénard の条件を満たす組、FFをその Liénard 系とし、σ,β,A,xmax⁡,Xα,Yα±\sigma,\beta,A,x_{\max},X_\alpha,Y^\pm_\alphaを補題 1.2 (3)のとおりとし、0<α1<α20<\alpha_1<\alpha_2とする。このとき、β(α1)<β(α2)\beta(\alpha_1)<\beta(\alpha_2)、xmax⁡(α1)<xmax⁡(α2)x_{\max}(\alpha_1)<x_{\max}(\alpha_2)である。η∈[−β(α1),α1]\eta\in[-\beta(\alpha_1),\alpha_1]ならばXα1(η)<Xα2(η)X_{\alpha_1}(\eta)<X_{\alpha_2}(\eta)である。ξ∈[0,xmax⁡(α1)]\xi\in[0,x_{\max}(\alpha_1)]ならばYα1+(ξ)<Yα2+(ξ)Y^+_{\alpha_1}(\xi)<Y^+_{\alpha_2}(\xi)かつYα2−(ξ)<Yα1−(ξ)Y^-_{\alpha_2}(\xi)<Y^-_{\alpha_1}(\xi)である。

証明.i=1,2i=1,2についてzi=(0,αi)z_i=(0,\alpha_i)、σi=σ(αi)\sigma_i=\sigma(\alpha_i)、βi=β(αi)\beta_i=\beta(\alpha_i)、wi=(0,−βi)w_i=(0,-\beta_i)、Ai=A(αi)A_i=A(\alpha_i)とおく。補題 1.2 (2)によりAi∩{x=0}={zi,wi}A_i\cap\{x=0\}=\{z_i,w_i\}であり、−βi<η<αi-\beta_i<\eta<\alpha_iならばXαi(η)>0X_{\alpha_i}(\eta)>0、Xαi(−βi)=Xαi(αi)=0X_{\alpha_i}(-\beta_i)=X_{\alpha_i}(\alpha_i)=0である。

A1∩A2=∅A_1\cap A_2=\varnothingである。実際、φt1(z1)=φt2(z2)\varphi_{t_1}(z_1)=\varphi_{t_2}(z_2)、ti∈[0,σi]t_i\in[0,\sigma_i]とし、t1≤t2t_1\le t_2とすると、§E13.15 定理 1.3によりz1=φt2−t1(z2)z_1=\varphi_{t_2-t_1}(z_2)、0≤t2−t1≤σ20\le t_2-t_1\le\sigma_2である。z1z_1の第一座標は00であるからt2−t1∈{0,σ2}t_2-t_1\in\{0,\sigma_2\}であり、t2−t1=0t_2-t_1=0ならばz1=z2z_1=z_2、t2−t1=σ2t_2-t_1=\sigma_2ならばz1=w2z_1=w_2となって、いずれもα1≠α2\alpha_1\ne\alpha_2、α1>0>−β2\alpha_1>0>-\beta_2に反する。t1>t2t_1>t_2の場合もz2=φt1−t2(z1)z_2=\varphi_{t_1-t_2}(z_1)から同様である。

I=[max⁡{−β1,−β2},α1]I=[\max\{-\beta_1,-\beta_2\},\alpha_1]とおき、η∈I\eta\in Iについてd(η)=Xα2(η)−Xα1(η)d(\eta)=X_{\alpha_2}(\eta)-X_{\alpha_1}(\eta)とおく。IIはXα1X_{\alpha_1}とXα2X_{\alpha_2}の定義域に含まれ、ddはII上で連続である。d(η)=0d(\eta)=0ならば(Xα1(η),η)∈A1∩A2(X_{\alpha_1}(\eta),\eta)\in A_1\cap A_2であるから、ddはII上で零にならない。α1∈(−β2,α2)\alpha_1\in(-\beta_2,\alpha_2)であるからd(α1)=Xα2(α1)>0d(\alpha_1)=X_{\alpha_2}(\alpha_1)>0であり、§E2.43 系 2.3によりII上でd>0d>0である。

β1>β2\beta_1>\beta_2と仮定すると、−β2∈(−β1,α1)-\beta_2\in(-\beta_1,\alpha_1)であるから−β2∈I-\beta_2\in Iであり、d(−β2)=−Xα1(−β2)<0d(-\beta_2)=-X_{\alpha_1}(-\beta_2)<0となってII上でd>0d>0であることに反する。β1=β2\beta_1=\beta_2ならばw1=w2∈A1∩A2w_1=w_2\in A_1\cap A_2となってA1∩A2=∅A_1\cap A_2=\varnothingに反する。したがってβ1<β2\beta_1<\beta_2であり、I=[−β1,α1]I=[-\beta_1,\alpha_1]である。ゆえにη∈[−β1,α1]\eta\in[-\beta_1,\alpha_1]ならばXα1(η)<Xα2(η)X_{\alpha_1}(\eta)<X_{\alpha_2}(\eta)である。xmax⁡(α1)=Xα1(η1)x_{\max}(\alpha_1)=X_{\alpha_1}(\eta_1)を満たすη1∈[−β1,α1]\eta_1\in[-\beta_1,\alpha_1]をとると、xmax⁡(α1)<Xα2(η1)≤xmax⁡(α2)x_{\max}(\alpha_1)<X_{\alpha_2}(\eta_1)\le x_{\max}(\alpha_2)である。

0≤ξ<xmax⁡(α2)0\le\xi<x_{\max}(\alpha_2)について

{η∈[−β2,α2]:Xα2(η)>ξ}=(Yα2−(ξ), Yα2+(ξ))\{\eta\in[-\beta_2,\alpha_2]:X_{\alpha_2}(\eta)>\xi\}=\bigl(Y^-_{\alpha_2}(\xi),\,Y^+_{\alpha_2}(\xi)\bigr)

である。実際、φt(z2)=(x(t),y(t))\varphi_t(z_2)=(x(t),y(t))と書き、tmt_mを補題 1.2 (3)のとおりとし、x(t+)=ξx(t^+)=\xiを満たすt+∈[0,tm]t^+\in[0,t_m]とx(t−)=ξx(t^-)=\xiを満たすt−∈[tm,σ2]t^-\in[t_m,\sigma_2]をとる。xxは[0,tm][0,t_m]上で狭義単調増加、[tm,σ2][t_m,\sigma_2]上で狭義単調減少であるから、t∈[0,σ2]t\in[0,\sigma_2]についてx(t)>ξx(t)>\xiであることは、t∈(t+,t−)t\in(t^+,t^-)であるための必要十分条件である。yyは[0,σ2][0,\sigma_2]から[−β2,α2][-\beta_2,\alpha_2]への狭義単調減少な全単射で、Xα2(y(t))=x(t)X_{\alpha_2}(y(t))=x(t)であるから、左辺はy((t+,t−))=(y(t−),y(t+))=(Yα2−(ξ),Yα2+(ξ))y((t^+,t^-))=(y(t^-),y(t^+))=(Y^-_{\alpha_2}(\xi),Y^+_{\alpha_2}(\xi))に等しい。

ξ∈[0,xmax⁡(α1)]\xi\in[0,x_{\max}(\alpha_1)]とすると、ξ<xmax⁡(α2)\xi<x_{\max}(\alpha_2)である。η±=Yα1±(ξ)\eta^\pm=Y^\pm_{\alpha_1}(\xi)は[−β1,α1][-\beta_1,\alpha_1]に属し、Xα1(η±)=ξX_{\alpha_1}(\eta^\pm)=\xiを満たすから、Xα2(η±)>Xα1(η±)=ξX_{\alpha_2}(\eta^\pm)>X_{\alpha_1}(\eta^\pm)=\xiである。Xα2(η)>ξX_{\alpha_2}(\eta)>\xiを満たすη\etaの集合の表示により、η±∈(Yα2−(ξ),Yα2+(ξ))\eta^\pm\in(Y^-_{\alpha_2}(\xi),Y^+_{\alpha_2}(\xi))である。したがって、Yα2−(ξ)<Yα1−(ξ)Y^-_{\alpha_2}(\xi)<Y^-_{\alpha_1}(\xi)かつYα1+(ξ)<Yα2+(ξ)Y^+_{\alpha_1}(\xi)<Y^+_{\alpha_2}(\xi)である。▨

補題 2.2.(Φ,g)(\Phi,g)を Liénard の条件を満たす組、aaをΦ\Phiの正の零点、FFをその Liénard 系、φ\varphiをFFの流れとし、σ,β,xmax⁡,Xα,Yα±\sigma,\beta,x_{\max},X_\alpha,Y^\pm_\alphaを補題 1.2 (3)のとおりとする。G(x)=∫0xg(s) dsG(x)=\int_0^xg(s)\,ds、E(x,y)=12y2+G(x)E(x,y)=\tfrac12y^2+G(x)、k(ξ)=−g(ξ)Φ(ξ)k(\xi)=-g(\xi)\Phi(\xi)とおき、α>0\alpha>0についてδ(α)=E(0,−β(α))−E(0,α)=12(β(α)2−α2)\delta(\alpha)=E(0,-\beta(\alpha))-E(0,\alpha)=\tfrac12(\beta(\alpha)^2-\alpha^2)とおく。M={α>0:xmax⁡(α)>a}M=\{\alpha>0:x_{\max}(\alpha)>a\}とし、α∈M\alpha\in Mについて

I1(α)=∫0ak(ξ)Yα+(ξ)−Φ(ξ) dξ,I2(α)=−∫Yα−(a)Yα+(a)Φ(Xα(η)) dη,I3(α)=∫0ak(ξ)Φ(ξ)−Yα−(ξ) dξI_1(\alpha)=\int_0^a\frac{k(\xi)}{Y^+_\alpha(\xi)-\Phi(\xi)}\,d\xi,\qquad I_2(\alpha)=-\int_{Y^-_\alpha(a)}^{Y^+_\alpha(a)}\Phi(X_\alpha(\eta))\,d\eta,\qquad I_3(\alpha)=\int_0^a\frac{k(\xi)}{\Phi(\xi)-Y^-_\alpha(\xi)}\,d\xi

とおく。このとき、次が成り立つ。

  1. z∈R2z\in\R^2、φt(z)=(x(t),y(t))\varphi_t(z)=(x(t),y(t))についてddtE(φt(z))=k(x(t))\tfrac{d}{dt}E(\varphi_t(z))=k(x(t))である。δ\deltaは(0,∞)(0,\infty)上で連続であり、α∈M\alpha\in Mならばδ(α)=I1(α)+I2(α)+I3(α)\delta(\alpha)=I_1(\alpha)+I_2(\alpha)+I_3(\alpha)である。
  2. α>0\alpha>0、α∉M\alpha\notin Mならばδ(α)>0\delta(\alpha)>0である。
  3. α∈M\alpha\in M、α′>α\alpha'>\alphaならばα′∈M\alpha'\in Mである。I1I_1とI3I_3はMM上で正の値をとり、I1,I2,I3I_1,I_2,I_3はいずれもMM上で狭義単調減少である。
  4. MMは空でなく、α→∞\alpha\to\inftyのときδ(α)→−∞\delta(\alpha)\to-\inftyである。

証明.主張 (1)を示す。§D1.19 定理 1.1によりG′=gG'=gであり、§E4.3 定理 1.1によりddtE(φt(z))=y y′+g(x) x′=−g(x)y+g(x)(y−Φ(x))=k(x)\tfrac{d}{dt}E(\varphi_t(z))=y\,y'+g(x)\,x'=-g(x)y+g(x)(y-\Phi(x))=k(x)である。δ\deltaの連続性は補題 1.2 (4)による。α∈M\alpha\in Mとし、z=(0,α)z=(0,\alpha)、σ=σ(α)\sigma=\sigma(\alpha)、tmt_mを補題 1.2 (3)のとおりとする。xmax⁡(α)>ax_{\max}(\alpha)>aであるから、x(ta)=ax(t_a)=aを満たすta∈(0,tm)t_a\in(0,t_m)とx(tb)=ax(t_b)=aを満たすtb∈(tm,σ)t_b\in(t_m,\sigma)がただ一つずつ存在する。§D1.19 定理 2.1により

δ(α)=E(φσ(z))−E(z)=∫0tak(x(t)) dt+∫tatbk(x(t)) dt+∫tbσk(x(t)) dt\delta(\alpha)=E(\varphi_\sigma(z))-E(z)=\int_0^{t_a}k(x(t))\,dt+\int_{t_a}^{t_b}k(x(t))\,dt+\int_{t_b}^{\sigma}k(x(t))\,dt

である。[0,ta][0,t_a]上でy(t)=Yα+(x(t))y(t)=Y^+_\alpha(x(t))であり、x′(t)=Yα+(x(t))−Φ(x(t))>0x'(t)=Y^+_\alpha(x(t))-\Phi(x(t))>0であるから、f1(ξ)=k(ξ)/(Yα+(ξ)−Φ(ξ))f_1(\xi)=k(\xi)/(Y^+_\alpha(\xi)-\Phi(\xi))は[0,a][0,a]上で連続でk(x(t))=f1(x(t))x′(t)k(x(t))=f_1(x(t))x'(t)である。§E4.18 補題 3.1により第一項は∫0af1=I1(α)\int_0^af_1=I_1(\alpha)である。[tb,σ][t_b,\sigma]上でy(t)=Yα−(x(t))y(t)=Y^-_\alpha(x(t))、x′(t)=Yα−(x(t))−Φ(x(t))<0x'(t)=Y^-_\alpha(x(t))-\Phi(x(t))<0であるから、同じ置換により第三項は∫a0k(ξ)/(Yα−(ξ)−Φ(ξ)) dξ=I3(α)\int_a^0k(\xi)/(Y^-_\alpha(\xi)-\Phi(\xi))\,d\xi=I_3(\alpha)である。[ta,tb][t_a,t_b]上でk(x(t))=Φ(x(t))y′(t)k(x(t))=\Phi(x(t))y'(t)、x(t)=Xα(y(t))x(t)=X_\alpha(y(t))であり、y(ta)=Yα+(a)y(t_a)=Y^+_\alpha(a)、y(tb)=Yα−(a)y(t_b)=Y^-_\alpha(a)であるから、§E4.18 補題 3.1により第二項は∫Yα+(a)Yα−(a)Φ(Xα(η)) dη=I2(α)\int_{Y^+_\alpha(a)}^{Y^-_\alpha(a)}\Phi(X_\alpha(\eta))\,d\eta=I_2(\alpha)である。

主張 (2)を示す。補題 1.2 (1)により0<ξ<a0<\xi<aでg(ξ)>0g(\xi)>0、Φ(ξ)<0\Phi(\xi)<0であるからk(ξ)>0k(\xi)>0であり、k(0)=k(a)=0k(0)=k(a)=0である。xmax⁡(α)≤ax_{\max}(\alpha)\le aならば、0≤t≤σ(α)0\le t\le\sigma(\alpha)で0≤x(t)≤a0\le x(t)\le aであるからk(x(t))≥0k(x(t))\ge0であり、0<t<tm0<t<t_mで0<x(t)<xmax⁡(α)≤a0<x(t)<x_{\max}(\alpha)\le aであるからk(x(t))>0k(x(t))>0である。tm/2t_m/2を含む閉区間で、その上で連続関数k(x(⋅))k(x(\cdot))が正の定数以上となるものが存在するから、§D1.17 命題 4.4によりδ(α)=∫0σ(α)k(x(t)) dt>0\delta(\alpha)=\int_0^{\sigma(\alpha)}k(x(t))\,dt>0である。

主張 (3)を示す。α∈M\alpha\in M、α′>α\alpha'>\alphaならば、補題 2.1によりxmax⁡(α′)>xmax⁡(α)>ax_{\max}(\alpha')>x_{\max}(\alpha)>aであり、α′∈M\alpha'\in Mである。0<ξ<a0<\xi<aでk(ξ)>0k(\xi)>0であり、補題 1.2 (3)によりYα+(ξ)−Φ(ξ)>0Y^+_\alpha(\xi)-\Phi(\xi)>0、Φ(ξ)−Yα−(ξ)>0\Phi(\xi)-Y^-_\alpha(\xi)>0であるから、I1(α)>0I_1(\alpha)>0、I3(α)>0I_3(\alpha)>0である。α<α′\alpha<\alpha'をMMの点とする。補題 2.1によりξ∈[0,a]\xi\in[0,a]でYα+(ξ)<Yα′+(ξ)Y^+_\alpha(\xi)<Y^+_{\alpha'}(\xi)、Yα′−(ξ)<Yα−(ξ)Y^-_{\alpha'}(\xi)<Y^-_\alpha(\xi)であるから、0<ξ<a0<\xi<aでI1I_1、I3I_3の被積分関数はα′\alpha'のほうが真に小さく、I1(α′)<I1(α)I_1(\alpha')<I_1(\alpha)、I3(α′)<I3(α)I_3(\alpha')<I_3(\alpha)である。Yα′−(a)<Yα−(a)<Yα+(a)<Yα′+(a)Y^-_{\alpha'}(a)<Y^-_\alpha(a)<Y^+_\alpha(a)<Y^+_{\alpha'}(a)である。η∈[Yα′−(a),Yα′+(a)]\eta\in[Y^-_{\alpha'}(a),Y^+_{\alpha'}(a)]に対応する時刻ではx≥ax\ge aであるからXα′(η)≥aX_{\alpha'}(\eta)\ge a、Φ(Xα′(η))≥0\Phi(X_{\alpha'}(\eta))\ge0である。η∈[Yα−(a),Yα+(a)]\eta\in[Y^-_\alpha(a),Y^+_\alpha(a)]ならばa≤Xα(η)<Xα′(η)a\le X_\alpha(\eta)<X_{\alpha'}(\eta)であるから、Φ\Phiの[a,∞)[a,\infty)上の狭義単調性によりΦ(Xα′(η))>Φ(Xα(η))\Phi(X_{\alpha'}(\eta))>\Phi(X_\alpha(\eta))である。したがって

∫Yα′−(a)Yα′+(a)Φ(Xα′(η)) dη≥∫Yα−(a)Yα+(a)Φ(Xα′(η)) dη>∫Yα−(a)Yα+(a)Φ(Xα(η)) dη\int_{Y^-_{\alpha'}(a)}^{Y^+_{\alpha'}(a)}\Phi(X_{\alpha'}(\eta))\,d\eta\ge\int_{Y^-_\alpha(a)}^{Y^+_\alpha(a)}\Phi(X_{\alpha'}(\eta))\,d\eta>\int_{Y^-_\alpha(a)}^{Y^+_\alpha(a)}\Phi(X_\alpha(\eta))\,d\eta

であり、I2(α′)<I2(α)I_2(\alpha')<I_2(\alpha)である。

主張 (4)を示す。x1>ax_1>aをとり、P=max⁡[0,x1]∣Φ∣P=\max_{[0,x_1]}\lvert\Phi\rvert、α♯=G(x1)+P+2\alpha_\sharp=G(x_1)+P+2とおく。α≥α♯\alpha\ge\alpha_\sharp、z=(0,α)z=(0,\alpha)とし、T={t∈[0,tm]:0≤s≤t ならば x′(s)≥1 かつ x(s)≤x1}\mathcal T=\{t\in[0,t_m]:0\le s\le t\text{ ならば }x'(s)\ge1\text{ かつ }x(s)\le x_1\}とおく。x′(0)=α≥2x'(0)=\alpha\ge2であるから0∈T0\in\mathcal Tであり、t′=sup⁡T∈Tt'=\sup\mathcal T\in\mathcal Tである。0≤t≤t′0\le t\le t'について、g≥0g\ge0([0,∞)[0,\infty)上)とx′≥1x'\ge1から∫0tg(x(s)) ds≤∫0tg(x(s))x′(s) ds=G(x(t))≤G(x1)\int_0^tg(x(s))\,ds\le\int_0^tg(x(s))x'(s)\,ds=G(x(t))\le G(x_1)であり、y(t)≥α−G(x1)y(t)\ge\alpha-G(x_1)、x′(t)≥α−G(x1)−P≥2x'(t)\ge\alpha-G(x_1)-P\ge2である。x′(tm)=0x'(t_m)=0であるからt′<tmt'<t_mであり、x(t′)<x1x(t')<x_1ならば連続性によりt′t'の少し先までT\mathcal Tの条件が成り立ってt′=sup⁡Tt'=\sup\mathcal Tに反する。したがってx(t′)=x1x(t')=x_1であり、x1<xmax⁡(α)x_1<x_{\max}(\alpha)、Yα+(x1)=y(t′)≥α−G(x1)Y^+_\alpha(x_1)=y(t')\ge\alpha-G(x_1)である。とくにα∈M\alpha\in Mであり、M≠∅M\ne\varnothingである。Yα+Y^+_\alphaは狭義単調減少、Yα−Y^-_\alphaは狭義単調増加であるから[Yα−(x1),Yα+(x1)]⊆[Yα−(a),Yα+(a)][Y^-_\alpha(x_1),Y^+_\alpha(x_1)]\subseteq[Y^-_\alpha(a),Y^+_\alpha(a)]であり、η∈[Yα−(x1),Yα+(x1)]\eta\in[Y^-_\alpha(x_1),Y^+_\alpha(x_1)]に対応する時刻ではx≥x1x\ge x_1であるからΦ(Xα(η))≥Φ(x1)>0\Phi(X_\alpha(\eta))\ge\Phi(x_1)>0である。その他のη∈[Yα−(a),Yα+(a)]\eta\in[Y^-_\alpha(a),Y^+_\alpha(a)]ではΦ(Xα(η))≥0\Phi(X_\alpha(\eta))\ge0である。補題 1.2 (3)によりYα−(x1)<Φ(x1)Y^-_\alpha(x_1)<\Phi(x_1)であるから

I2(α)≤−Φ(x1)(Yα+(x1)−Yα−(x1))≤−Φ(x1)(α−G(x1)−Φ(x1))I_2(\alpha)\le-\Phi(x_1)\bigl(Y^+_\alpha(x_1)-Y^-_\alpha(x_1)\bigr)\le-\Phi(x_1)\bigl(\alpha-G(x_1)-\Phi(x_1)\bigr)

である。主張 (3)によりα≥α♯\alpha\ge\alpha_\sharpでI1(α)≤I1(α♯)I_1(\alpha)\le I_1(\alpha_\sharp)、I3(α)≤I3(α♯)I_3(\alpha)\le I_3(\alpha_\sharp)であるから、δ(α)≤I1(α♯)+I3(α♯)−Φ(x1)(α−G(x1)−Φ(x1))→−∞\delta(\alpha)\le I_1(\alpha_\sharp)+I_3(\alpha_\sharp)-\Phi(x_1)(\alpha-G(x_1)-\Phi(x_1))\to-\infty(α→∞\alpha\to\infty)である。▨

3 周期軌道の存在と一意性・吸引

定理 3.1 (Liénard の定理).(Φ,g)(\Phi,g)を Liénard の条件を満たす組、FFをその Liénard 系、φ\varphiをFFの流れとする。このときFFはただ一つの周期軌道Γ∗\Gamma_*をもち、すべてのz∈R2∖{0}z\in\R^2\setminus\{0\}について[0,∞)⊆I(z)[0,\infty)\subseteq I(z)、ω(z)=Γ∗\omega(z)=\Gamma_*である。Γ∗\Gamma_*は極限閉軌道である。

証明.a,σ,β,A,xmax⁡,G,E,δ,Ma,\sigma,\beta,A,x_{\max},G,E,\delta,Mを補題 1.2 (3)と補題 2.2のとおりとする。

c=G(a)/2c=G(a)/2とおく。ggは奇関数で(0,∞)(0,\infty)上で正であるから、GGは偶関数で[0,∞)[0,\infty)上で狭義単調増加であり、E(x,y)≤cE(x,y)\le cならばG(x)≤c<G(a)G(x)\le c<G(a)であり∣x∣<a\lvert x\rvert<aである。補題 1.2 (1)とg,Φg,\Phiが奇関数であることにより、∣x∣<a\lvert x\rvert<aならばk(x)=−g(x)Φ(x)≥0k(x)=-g(x)\Phi(x)\ge0である。z≠0z\ne0とし、mz=min⁡{E(z),c}>0m_z=\min\{E(z),c\}>0とおく。

主張 3.1.1.t∈I(z)t\in I(z)、t≥0t\ge0ならばE(φt(z))≥mzE(\varphi_t(z))\ge m_zである。とくに0∉ω(z)0\notin\omega(z)である。

証明.E(φt2(z))<mzE(\varphi_{t_2}(z))<m_zを満たすt2>0t_2>0があると仮定し、t3=sup⁡{t∈[0,t2]:E(φt(z))≥mz}t_3=\sup\{t\in[0,t_2]:E(\varphi_t(z))\ge m_z\}とおく。連続性によりE(φt3(z))=mzE(\varphi_{t_3}(z))=m_z、t3<t2t_3<t_2であり、(t3,t2](t_3,t_2]上でE(φt(z))<mz≤cE(\varphi_t(z))<m_z\le cであるから、補題 2.2 (1)により(t3,t2](t_3,t_2]上でddtE(φt(z))≥0\tfrac{d}{dt}E(\varphi_t(z))\ge0である。したがってE(φt2(z))≥E(φt3(z))=mzE(\varphi_{t_2}(z))\ge E(\varphi_{t_3}(z))=m_zであり、仮定に反する。ω(z)\omega(z)の点はE≥mzE\ge m_zを満たす点の極限であり、EEは連続でE(0)=0E(0)=0であるから、0∉ω(z)0\notin\omega(z)である。▨

補題 2.2 (4)によりδ(α1)<0\delta(\alpha_1)<0を満たすα1∈M\alpha_1\in Mが存在する。補題 2.2 (3)と補題 2.2 (1)により、α≥α1\alpha\ge\alpha_1ならばα∈M\alpha\in M、δ(α)≤δ(α1)<0\delta(\alpha)\le\delta(\alpha_1)<0であり、β(α)<α\beta(\alpha)<\alphaである。

z≠0z\ne0とする。t1≥0t_1\ge0とh1>0h_1>0で、φt1(z)=(0,±h1)\varphi_{t_1}(z)=(0,\pm h_1)を満たすものが存在する。実際、zzがyy軸上にあればt1=0t_1=0であり、zzの第一座標が正ならば補題 1.2 (2)によりt1=σ(z)t_1=\sigma(z)であり、負ならば補題 1.2 (1)によりφt(z)=−φt(−z)\varphi_t(z)=-\varphi_t(-z)であるからt1=σ(−z)t_1=\sigma(-z)である。k≥1k\ge1についてtk+1=tk+σ(hk)t_{k+1}=t_k+\sigma(h_k)、hk+1=β(hk)h_{k+1}=\beta(h_k)と定める。§E13.15 定理 1.3と補題 1.2 (1)により、φtk(z)=(0,±hk)\varphi_{t_k}(z)=(0,\pm h_k)ならばtk+1∈I(z)t_{k+1}\in I(z)、φtk+1(z)=(0,∓hk+1)\varphi_{t_{k+1}}(z)=(0,\mp h_{k+1})であり、{φt(z):tk≤t≤tk+1}=±A(hk)\{\varphi_t(z):t_k\le t\le t_{k+1}\}=\pm A(h_k)である。H=max⁡{h1,α1}H=\max\{h_1,\alpha_1\}とおくと、すべてのkkでhk≤Hh_k\le Hである。実際、hk≤Hh_k\le Hとすると、hk≥α1h_k\ge\alpha_1ならばhk+1=β(hk)<hk≤Hh_{k+1}=\beta(h_k)<h_k\le Hであり、hk<α1h_k<\alpha_1ならば補題 2.1によりhk+1=β(hk)<β(α1)<α1≤Hh_{k+1}=\beta(h_k)<\beta(\alpha_1)<\alpha_1\le Hである。主張 3.1.1により12hk2=E(φtk(z))≥mz\tfrac12h_k^2=E(\varphi_{t_k}(z))\ge m_zであるから、hk≥h0=2mzh_k\ge h_0=\sqrt{2m_z}である。補題 1.2 (4)によりσ\sigmaは[h0,H][h_0,H]上で正の最小値σ0\sigma_0をとり、tk≥t1+(k−1)σ0→∞t_k\ge t_1+(k-1)\sigma_0\to\inftyである。したがって[0,∞)⊆I(z)[0,\infty)\subseteq I(z)である。補題 2.1によりh≤Hh\le HならばA(h)⊆[0,xmax⁡(H)]×[−β(H),H]A(h)\subseteq[0,x_{\max}(H)]\times[-\beta(H),H]であるから、前方軌道は

Kz={φt(z):0≤t≤t1}∪([−xmax⁡(H),xmax⁡(H)]×[−H′,H′]),H′=max⁡{H,β(H)}K_z=\{\varphi_t(z):0\le t\le t_1\}\cup\bigl([-x_{\max}(H),x_{\max}(H)]\times[-H',H']\bigr),\qquad H'=\max\{H,\beta(H)\}

というコンパクト集合に含まれる。

z≠0z\ne0ならば、FFの平衡点は原点だけであり(補題 1.2 (1))、主張 3.1.1により0∉ω(z)0\notin\omega(z)であるから、§E13.18 定理 3.4をzzとKzK_zに適用してω(z)\omega(z)は周期軌道である。とくにFFは周期軌道をもつ。

Γ\Gammaを周期軌道とする。Γ\Gammaの点zzの軌道は正と負のyy軸へ交互に到達するから、(0,α)∈Γ(0,\alpha)\in\Gammaを満たすα>0\alpha>0が存在する。0<c1<c20<c_1<c_2についてΣ={(0,y):c1≤y≤c2}\Sigma=\{(0,y):c_1\le y\le c_2\}は、F(0,y)=(y,0)F(0,y)=(y,0)とdet⁡((0,1),(y,0))=−y≠0\det((0,1),(y,0))=-y\ne0により横断線分である。(0,α)(0,\alpha)から時間σ(α)+σ(β(α))\sigma(\alpha)+\sigma(\beta(\alpha))だけ進んだ点は(0,β(β(α)))∈Γ(0,\beta(\beta(\alpha)))\in\Gammaである。c1<min⁡{α,β(β(α))}c_1<\min\{\alpha,\beta(\beta(\alpha))\}、c2>max⁡{α,β(β(α))}c_2>\max\{\alpha,\beta(\beta(\alpha))\}ととると、§E13.18 補題 2.3 (3)によりΓ∩Σ\Gamma\cap\Sigmaは高々一点であるからβ(β(α))=α\beta(\beta(\alpha))=\alphaである。β\betaは補題 2.1により狭義単調増加であるから、β(α)>α\beta(\alpha)>\alphaならばβ(β(α))>β(α)>α\beta(\beta(\alpha))>\beta(\alpha)>\alpha、β(α)<α\beta(\alpha)<\alphaならばβ(β(α))<α\beta(\beta(\alpha))<\alphaとなるので、β(α)=α\beta(\alpha)=\alpha、δ(α)=0\delta(\alpha)=0である。(0,α),(0,α′)∈Γ(0,\alpha),(0,\alpha')\in\Gamma(α,α′>0\alpha,\alpha'>0)ならば、c1<min⁡{α,α′}c_1<\min\{\alpha,\alpha'\}、c2>max⁡{α,α′}c_2>\max\{\alpha,\alpha'\}としたΣ\Sigmaに両者が属するから、§E13.18 補題 2.3 (3)によりα=α′\alpha=\alpha'である。したがってΓ\Gammaと正のyy軸の共通部分は一点(0,αΓ)(0,\alpha_\Gamma)であり、δ(αΓ)=0\delta(\alpha_\Gamma)=0である。

補題 2.2 (2)によりδ\deltaの零点はMMに属し、補題 2.2 (1)と補題 2.2 (3)によりδ=I1+I2+I3\delta=I_1+I_2+I_3はMM上で狭義単調減少であるから、δ\deltaの零点は高々一つである。周期軌道Γ,Γ′\Gamma,\Gamma'についてδ(αΓ)=δ(αΓ′)=0\delta(\alpha_\Gamma)=\delta(\alpha_{\Gamma'})=0であるからαΓ=αΓ′\alpha_\Gamma=\alpha_{\Gamma'}であり、Γ\GammaとΓ′\Gamma'は点(0,αΓ)(0,\alpha_\Gamma)を共有するから、§E13.15 命題 4.3 (1)によりΓ=Γ′\Gamma=\Gamma'である。したがってFFの周期軌道はただ一つであり、これをΓ∗\Gamma_*とおく。

z≠0z\ne0ならばω(z)\omega(z)は周期軌道であるからω(z)=Γ∗\omega(z)=\Gamma_*である。R2∖(Γ∗∪{0})\R^2\setminus(\Gamma_*\cup\{0\})の点zzをとると、zzは平衡点でなく、γ(z)≠Γ∗\gamma(z)\ne\Gamma_*であるから周期点でもなく、z∉Γ∗z\notin\Gamma_*、ω(z)=Γ∗\omega(z)=\Gamma_*であるから、Γ∗\Gamma_*は極限閉軌道である。▨

4 負減衰条件の必要性と van der Pol 方程式

例 4.1.Φ(x)=x(x2−1)2\Phi(x)=x(x^2-1)^2、g(x)=xg(x)=xとし、F(x,y)=(y−Φ(x), −x)F(x,y)=(y-\Phi(x),\,-x)を(Φ,g)(\Phi,g)の Liénard 系、φ\varphiをFFの流れとする。このとき、Φ\Phiとggは奇関数であり、x≠0x\ne0ならばxg(x)=x2>0xg(x)=x^2>0である。Φ\Phiの(0,∞)(0,\infty)における零点は11だけであり、Φ′(x)=(x2−1)(5x2−1)\Phi'(x)=(x^2-1)(5x^2-1)は(1,∞)(1,\infty)上で正であるから、Φ\Phiは[1,∞)[1,\infty)上で狭義単調増加である。x≥2x\ge2ならばΦ(x)≥x\Phi(x)\ge xであるから、x→∞x\to\inftyのときΦ(x)→∞\Phi(x)\to\inftyである。したがって(Φ,g)(\Phi,g)は定義 1.1 条件 (a)、定義 1.1 条件 (b)、定義 1.1 条件 (d)を満たし、Φ′(0)=1\Phi'(0)=1であるから定義 1.1 条件 (c)を満たさない。

FFは周期軌道をもたない。実際、z∈R2z\in\R^2を周期点、TTをその最小周期とし、φt(z)=(x(t),y(t))\varphi_t(z)=(x(t),y(t))、E(x,y)=12(x2+y2)E(x,y)=\tfrac12(x^2+y^2)、f(t)=x(t)2(x(t)2−1)2f(t)=x(t)^2(x(t)^2-1)^2とおく。§E13.15 命題 4.3 (3)によりI(z)=RI(z)=\R、φT(z)=z\varphi_T(z)=zである。§E4.3 定理 1.1によりddtE(φt(z))=x(y−Φ(x))−yx=−xΦ(x)=−f(t)\tfrac{d}{dt}E(\varphi_t(z))=x(y-\Phi(x))-yx=-x\Phi(x)=-f(t)であるから、§D1.19 定理 2.1により∫0Tf(t) dt=E(z)−E(φT(z))=0\int_0^Tf(t)\,dt=E(z)-E(\varphi_T(z))=0である。ffは[0,T][0,T]上で連続かつ非負である。f(t0)>0f(t_0)>0を満たすt0∈[0,T]t_0\in[0,T]があるとすると、t0t_0を含む閉区間で、その上でffが正の定数以上となるものが存在するから、§D1.17 命題 4.4により∫0Tf(t) dt>0\int_0^Tf(t)\,dt>0であり、これは∫0Tf(t) dt=0\int_0^Tf(t)\,dt=0に反する。したがって[0,T][0,T]上でf=0f=0、すなわちx(t)∈{−1,0,1}x(t)\in\{-1,0,1\}であり、§E2.43 系 2.3によりxxは[0,T][0,T]上で定数である。[0,T][0,T]上で0=x′=y−Φ(x)0=x'=y-\Phi(x)、Φ(x)=0\Phi(x)=0であるからy=0y=0であり、0=y′=−x0=y'=-xであるからx=0x=0である。したがってz=0z=0であり、F(z)=0F(z)=0はzzが周期点であることに反する。定理 3.1の仮定定義 1.1 条件 (c)を外すと、結論が成り立たない。

例 4.2 (van der Pol 方程式).μ>0\mu>0とし、van der Pol 方程式u′′−μ(1−u2)u′+u=0u''-\mu(1-u^2)u'+u=0をR2\R^2上の一階系

(u,v)′=H(u,v),H(u,v)=(v, μ(1−u2)v−u)(u,v)'=H(u,v),\qquad H(u,v)=\bigl(v,\ \mu(1-u^2)v-u\bigr)

として考える。Φ(x)=μ(x3/3−x)\Phi(x)=\mu(x^3/3-x)、g(x)=xg(x)=xとおく。このとき、Φ\Phiとggは奇関数であり、x≠0x\ne0でxg(x)=x2>0xg(x)=x^2>0、Φ′(0)=−μ<0\Phi'(0)=-\mu<0である。Φ(x)=0\Phi(x)=0かつx>0x>0であることは、x=3x=\sqrt3であるための必要十分条件である。x>1x>1でΦ′(x)=μ(x2−1)>0\Phi'(x)=\mu(x^2-1)>0であるからΦ\Phiは[3,∞)[\sqrt3,\infty)上で狭義単調増加で、x→∞x\to\inftyのときΦ(x)→∞\Phi(x)\to\inftyである。したがって(Φ,g)(\Phi,g)は Liénard の条件を満たす。FFをその Liénard 系とする。

h(x,y)=(x,y−Φ(x))h(x,y)=(x,y-\Phi(x))はR2\R^2からR2\R^2への同相写像で、逆写像は(u,v)↦(u,v+Φ(u))(u,v)\mapsto(u,v+\Phi(u))である。(x(t),y(t))(x(t),y(t))がFFの解ならば、u=xu=x、v=y−Φ(x)v=y-\Phi(x)はu′=y−Φ(x)=vu'=y-\Phi(x)=v、v′=−x−Φ′(x)(y−Φ(x))=μ(1−u2)v−uv'=-x-\Phi'(x)(y-\Phi(x))=\mu(1-u^2)v-uを満たし、HHの解である。(u(t),v(t))(u(t),v(t))がHHの解ならば、x=ux=u、y=v+Φ(u)y=v+\Phi(u)はx′=v=y−Φ(x)x'=v=y-\Phi(x)、y′=μ(1−u2)v−u+μ(u2−1)v=−g(x)y'=\mu(1-u^2)v-u+\mu(u^2-1)v=-g(x)を満たし、FFの解である。§E13.15 補題 1.2によりIH(h(z))=IF(z)I_H(h(z))=I_F(z)であり、h(φtF(z))=φtH(h(z))h(\varphi^F_t(z))=\varphi^H_t(h(z))であるから、hhはFFとHHの流れの位相共役である。定理 3.1によりFFはただ一つの周期軌道Γ∗\Gamma_*をもつ。§E13.15 命題 7.2 (1)によりHHの周期軌道はh(Γ∗)h(\Gamma_*)だけであり、§E13.15 命題 7.2 (2)によりw≠0w\ne0についてωH(w)=h(ωF(h−1(w)))=h(Γ∗)\omega_H(w)=h(\omega_F(h^{-1}(w)))=h(\Gamma_*)である。したがってh(Γ∗)h(\Gamma_*)はHHのただ一つの周期軌道であり、原点以外のすべての点のω\omega極限集合であって、極限閉軌道である。

参考文献

  1. Lawrence Perko, Differential Equations and Dynamical Systems, 3rd ed., Texts in Applied Mathematics 7, Springer, New York, 2001.Liénard 系の半周の軌道とエネルギー差の比較による周期軌道の存在・一意性・吸引の証明を参考にした。

前提記事

10 本の記事・単元を表示