• 検索結果がありません。

A 型笹野系の有理解の分類

N/A
N/A
Protected

Academic year: 2021

シェア "A 型笹野系の有理解の分類"

Copied!
35
0
0

読み込み中.... (全文を見る)

全文

(1)

A(1)4

型笹野系の有理解の分類

三重大学大学院

教育学研究科教科教育専攻数学教育専修

野田 真司

(2)

目次

序文 1

1 Painlev´e方程式に関する序論 2

1.1 ベックルント変換 . . . . 3

1.2 Painlev´e方程式の歴史. . . . 4

1.3 Painlev´e II型方程式のハミルトン系表示とベックルント変換 . . . . 5

1.4 Painlev´e II型方程式の有理解 . . . . 8

1.5 τ関数 . . . . 10

2 A(1)4 型笹野系の有理解の分類 13 3 行列IIPainlev´e方程式 18 3.1 q1= 0およびq2= 0での行列IIPainlev´e方程式の有理解 . . . . 18

3.2 行列IIPainlev´e方程式のベックルント変換 . . . . 19

3.3 q16= 0, q26= 0での行列IIPainlev´e方程式の有理解の構成 . . . . 22

(3)

序文

 この論文では,岡本和夫氏や野海正俊氏など数々の数学者が研究を行なっているPainlev´e方程式を拡張した

A(1)4 型笹野系」および「行列IIPainlev´e方程式」と呼ばれる方程式系の,有理解の構成について述べてい .

 第1章では, Painlev´e方程式の序論として,本論文を読み進めるにあたって必要な用語の定義や,本研究に有 用な定理等の説明, Painlev´e方程式の歴史,そしてPainlev´e II型方程式のハミルトン系表示とベックルント変 換から有理解の構成方法について記述している.さらに, Painlev´e方程式にとって重要な役割を果たすτ関数 を導入して,これを用いた有理解の表し方についても述べている.

 第2章では, Painlev´e II型方程式を拡張した「A(1)4 型笹野系」について,ハミルトニアンの満たす微分方程 式および,H と独立変数zとの関係式を導き出し,その2つの方程式の有理解がYablonskii - Vorob’evの多項 式を用いて記述できることを示した.そして,その事実を利用してA(1)4 型笹野系の有理解を決定している.  第3章では,東京大学大学院数理科学研究科の川上拓志氏が発見した「行列IIPainlev´e方程式」につい , この方程式系が有理解をもつときのパラメータの条件を調べ,三重大学大学院の吉田和史氏が導き出した ベックルント変換を利用して有理解を構成している.

 研究対象に 有理解 を選んだ理由としては,有理解を調べることにより,変換群の構造がある程度分かり, らにいくつかある高階Painlev´e方程式の分類に有用であると考えたからである.

 最後に,本論文の作成にあたり,熱心かつ丁寧に指導してくださった石谷寛教授および川向洋之准教授,そし て本研究に際して 行列IIPainlev´e方程式 という大変貴重な議題を提供してくださった川上拓志氏,並び に研究仲間の吉田和史氏に深い感謝をし,序文とする.

(4)

1 Painlev´e 方程式に関する序論

 この章では,本論文を読み進めるために必要な用語や定理,そしてPainlev´e方程式が発見されるまでの歴史 およびその後の様々な研究を紹介し,さらにPainlev´e方程式の有理解の構成方法を, II型方程式PIIを例に とって紹介していく.

 複素領域で定義された2階の有理関数係数の線形常微分方程式

y00+p(x)y0+q(x)y= 0 (1.1)

について考察していく(ただし,0 =d/dxである) .

定義1.1 微分方程式(1.1)について,p(x), q(x)x=aで極をもち, (1.1)のどの解もx=aで正則であると ,x=aを見かけの特異点であるという.

定義1.2 方程式(1.1)におけるp(x), q(x)x=cでの極の位数がそれぞれ高々1,2であるとき,x=c

(1.1)の確定特異点といい,そうでない特異点を不確定特異点という.

x=cp(x), q(x)の確定特異点になっているとき, (1.1) p(x) = p1

xc +p0+p1(xc) +· · · q(x) = q2

(xc)2 + q1

xc +q0+q1(xc) +· · · y= (xc)ρ{c0+c1(xc) +c2(xc)2+· · · } を代入したときの(xc)ρ2の係数をみると

{ρ(ρ1) +p1ρ+q2}c0

となる.

定義1.3 f(ρ) = ρ(ρ1) +p1ρ+q2としたときの変数ρについての代数方程式f(ρ) = 0, (1.1) x=cにおける決定方程式といい,その根ρ1, ρ2を特性指数という.

定義1.4 リーマン球面P1上の特異点がすべて確定特異点である微分方程式をFuchs型であるという.

定義1.5 方程式(1.1)に含まれるパラメータのうちで,特性指数と無関係なものをアクセサリー・パラメータ

という.

(1.1)の1次独立な解をy1(x), y2(x)とすると,y1(x), y2(x)は解空間の基底になる.,x=cを正の向きに 一周する道Lに沿ってy1(x), y2(x)を解析接続したものをLy1(x), Ly2(x)と表すと,Ly1(x), Ly2(x) y1(x), y2(x)の線形結合で表すことができる.すなわち

L

( y1(x) y2(x)

)

=ML

( y1(x) y2(x)

)

, ML=

( α β γ δ

)

(1.2) と表される(α, β, γ, δC).

(5)

定義1.6 関係式(1.2)における行列MLをモノドロミー行列という.

定義1.7 f(z)0<|za|< Rにおいて正則な関数とし,z=aにおけるローラン展開を f(z) =

n=1

cn

(za)n +

n=0

cn(za)n とする.このとき,係数c1z=aにおける留数といい, Res

z=af(z)dzと表す.

定義1.8 f(z)R <|z|<において正則な関数とする.このときf(1/w)0<|w|<1/Rで正則な関数 となる.f(1/w)w= 0におけるローラン展開を

f ( 1

w )

=

n=1

cn

wn +

n=1

cnwn としたときのc1,z=におけるf(z)の留数といい, Res

z=f(z)dzと表す. 定義1.9 1変数関数f(x), g(x)に対して,Dnf·g,次が成り立つものとして定める.

Dnf ·g=

n j=0

(1)jnCj

dnjf dxnj

djg dxj . この作用素(f, g) Dnf·gを広田微分という.

定理1.1 Cを単純閉曲線,関数f(z)Cの内部にある有限個の特異点a1, a2, · · ·, anを除いて正則である とする.このとき次が成り立つ.

C

f(z)dz= 2πi

n k=1

z=aReskf(z)dz .

定理1.2 f(z)C上の点a1, a2,· · ·, an z =で極をもち,それ以外のC上の点で正則であるとする. このとき,次が成立する.

n k=1

z=aResk

f(z)dz+ Res

z=f(z)dz= 0.

1.1 ベックルント変換

 1階の連立微分方程式

q0 =pq2 t 2 p0= 2qp+α+ 1

2

(1.3)

を考える(ただし,0 =d/dt,αは複素パラメータである).この方程式に対して,次の変換

Q=q , P=p+ 2q2+t (1.4)

を施すと, (1.3)p, q, αP, Q,αにしたものに等しくなる.実際 Q0 =q0 =

(

pq2 t 2

)

=PQ2 t 2 P0 =p0+ 4qq0+ 1 =

(

2qp+α+ 1 2

) + 4q

(

pq2 t 2

)

+ 1 = 2QP α+ 1 2

(6)

となり,確かにパラメータがαに変換されている.

定義1.10 (1.4)のような,パラメータを動かすことは許して方程式を不変に保つ有理的な変換をベックルント

変換という.

1.2 Painlev´e方程式の歴史 (1) Painlev´e方程式が発見された経緯

 物理的な現象を数式で表すと,微分方程式が現れることが多い.しかし,大抵の場合,その微分方程式の一般解 は初等的な関数で記述することはできない.そこで, 19世紀の解析学の問題意識として

微分方程式で定義される新しい特殊関数を見つける

というものが生まれた.そこで, 「極以外の解の特異点の位置が初期条件に依存しない」という性質(このよ うな性質をPainlev´e性という)を持つ代数的常微分方程式を考えることになった.線形微分方程式は常に Painlev´e性を持つが,非線形微分方程式でPainlev´e性を持つものは稀である. Painlev´e性を持つ1階方程式は 必ずリッカチ方程式,もしくは楕円関数の満たす微分方程式に変形できることがポワンカレとフックスによっ て示された.

P.Painlev´e, Painlev´e性をもつ2階微分方程式を研究していて,そのような方程式でtの解析関数を係数 とする有理関数Rを用いて

y00=R(y0, y , t)

の形に表されるものは,適当な変形を加える違いを除いて50個の「標準形」に直すことができることを発見し .さらに,このうちの44個は既知の関数を用いて解けることを示し,解として新たな特殊関数の導入が必要 である方程式は6個であることが分かった.

定義1.11 極以外の解の特異点の位置が初期条件に依存しないというPainlev´e性を備えた次の6種類の2階 非線形常微分方程式をPainlev´e方程式といい,それらの解をPainlev´e超越関数という.

PI y00= 6y2+t PII y00= 2y3+ty+α PIII y00= 1

y(y0)2 1 ty0+ 1

t(αy2+β) +γy3+ δ y PIV y00= 1

2y(y0)2+ 3

2y3+ 4ty2+ 2(t2α)y+ β y PV y00=

( 1 2y + 1

y1 )

(y0)2 1

ty0+(y1)2 t2

( αy+ β

y )

+γy

t +δy(y+ 1) y1 PVI y00= 1

2 (1

y + 1

y1 + 1 yt

) (y0)2

(1 t + 1

t1+ 1 yt

) y0

+y(y1)(yt) t2(t1)2

{

α+β t

y2 +γ t1

(y1)2 +δt(t1) (yt)2

}

ここで,0 =d/dtであり,α, β, γ, δは複素パラメータである.

(7)

(2) Painlev´e方程式発見後の進展(モノドロミー保存変形,ベックルント変換,ハミルトン構造)

Painlev´e方程式の発見後, R.Fuchs,P1 上に4つの確定特異点と1つの見かけの特異点をもつ2階の

Fuchs型方程式を研究していた.そして1905年に,この方程式のモノドロミー表現を不変にするためには,アク

セサリー・パラメータの1つがPVIを満たさなければならないことを発見した. L.ShlesingerはこのR.Fuchs の結果をうけ,ベックルント変換を発見した.また,確定特異点だけでなく不確定特異点に対して, R.Fuchs 同様の計算がR.Garnierによって行なわれ,PIからPVが導き出された.しかし, R.Garnierは形式的に計算 をしただけであったが,上野喜三雄氏によって正確な意味づけがなされた.さらに岡本和夫氏は,これらPI PVIのハミルトン構造を発見した.

(3) Painlev´e方程式の規約性に関する問題

Painlev´e方程式の一般解は本当に新しい関数になっているのかという問題があった.この問題に現代的な

定式化を与えたのは,梅村浩氏である.梅村氏は古典関数を

有理関数から出発して,既知関数の加減乗除と微分,既知関数を係数とする代数方程式を解く,既知関数を係数 とする線形常微分方程式を解く,アーベル関数に既知関数を代入する以上の操作を有限回繰り返して得られ るもの

と定義し,古典関数でないものを 新しい関数 とした.さらに梅村氏は,Painlev´e方程式の解が古典関数な らば,代数関数かまたは線形方程式の解を使って表すことができる」ということを証明した.このことにより

Painlev´e方程式の古典解をすべて決定する」という問題が生まれた.この問題は2008年にPVIまでの代数

解が完全に発見されたことにより解決し, Painlev´e方程式の一般解は本当に新しい関数になっているというこ とが分かった.

(4)高階Painlev´e方程式

Painlev´e方程式の研究が盛んになるにつれ,この方程式を拡張する動きも盛んになってきた. Painlev´e性を もつ方程式について,2階の場合は定義1.11に挙げた6種類のみであることが分かっているので,それよりも 高階の方程式が研究されるようになってきた.その中で,東京大学准教授の坂井秀隆氏は,高階のPainlev´e方程 式にはどれほどの種類があるかを考えるようになった.東京大学教授の大島利雄氏の分類によると,アクセサ リー・パラメータを4つ持つFuchs型方程式は13個あり,そのうち(連続的な)モノドロミー保存変形を許す ものは4個ある.そのうち1つは変形の次元が2次元であり,2変数Garnier系に対応している.坂井氏は, の3つの変形方程式を計算し,多項式ハミルトン系としての記述を得た.それらは藤・鈴木系,笹野系,行列型の Painlev´e方程式と呼ばれている.

1.3 Painlev´e II型方程式のハミルトン系表示とベックルント変換

 この論文では, Painlev´e方程式を拡張した高階化微分方程式の有理解をすべて決定することを研究の目的と している.ここでは,拡張する以前のPainlev´e II型方程式PIIのハミルトン系表示と,有理解を求めるために重 要な役割を果たすベックルント変換について再び触れることにする.

Painlev´e II型方程式

PII: y00= 2y3+ty+α (1.5)

について考察していく(ここで,0 =d/dt,αは複素パラメータである).高階の非線形微分方程式は,それと 等価な1階の連立形の方程式に書き直しておく方が都合の良いことも多い.そこで,

q=y , p=y0+y2+ t 2

(8)

とおいて,q, pの方程式に書き改めると, (1.5)は次の連立形の方程式と等価である.

HII:

q0=pq2 t 2 p0 = 2qp+α+ 1

2

(1.6)

,3つの変数q, p, tについての多項式H=H(q, p;t) H = 1

2p2 (

q2+ t 2

) p

( α+ 1

2 )

q (1.7)

で定義する.このとき,HII

q0= ∂H

∂p , p0 =∂H

∂q (1.8)

と表すことができる.

定義1.12 1個の関数H =H(q, p;t)を用いて(1.8)のような形に表される連立形の方程式をハミルトン系と いい,H をそのハミルトニアンという.

 以上のことをまとめると,次の命題が得られる. 命題1.1 Painlev´e II型方程式PII

H = 1 2p2

( q2+ t

2 )

p (

α+ 1 2

)

q (1.9)

をハミルトニアンとする次のハミルトン系と等価である.

HII:

q0= ∂H

∂p =pq2 t 2 p0=∂H

∂q = 2qp+α+ 1 2

(1.10)

 ここで,PIIの基本的なベックルント変換を,構成方法にも触れていくつか紹介しておく.まず,H0, H00 q, p, tで表すと

H0=1

2p , H00=1 2

(

2qp+α+ 1 2

)

となる.また,これよりq, pH0, H00で表したものをそれぞれϕ(H0, H00, α), ψ(H0, H00, α)とすると q=ϕ(H0, H00, α) = 2H00+α+ 1/2

4H0 , p=ψ(H0, H00, α) =2H0 (1.11) となる.このq, p(1.9)に代入すると,ハミルトニアンが満たす微分方程式

(H00)2+ 4 (H0)3+ 2H0(tH0H) 1 4

( α+ 1

2 )2

= 0 (1.12)

が得られる.逆に, (1.12)を満たすH に対してq, p(1.11)で定めると,これらは(1.10)を満たすことが言え .ところで,Hの満たす微分方程式(1.12),αα1に置き換えても不変である.従って,

Q=ϕ(H0, H00,α1), P =ψ(H0, H00,α1) (1.13)

(9)

とおくと,これらは(1.10)q, p, αQ, P,α1 に置き換えた方程式系を満たすことが分かる. よっ , (1.11), (1.13)からH0, H00を消去して得られるq, p, Q, Pの関係式

Q=q+ α+ 1/2

p , P=p .

,PIIαα1にする変換を与えていることが分かる. このことから,次の定理が得られる. 定理1.3 次のベックルント変換

Q=q+ α+ 1/2

p , P =p , H¯ =H (1.14)

を行なって得られるハミルトン系表示 dQ

dt = H¯

∂P , dP

dt =H¯

∂Q

, (1.10)q, p, αQ, P,α1にしたものに等しい.また,このときのハミルトニアンH¯ は次のとおり である.

H¯ = 1 2P2

( Q2+ t

2 )

P+ (

α+ 1 2

) Q . HIIには他にも基本的なベックルント変換があるので,それを以下にまとめておく.

定理1.4 Painlev´e II型方程式のハミルトン系HII,次の基本的なベックルント変換をもつ.

r1:

Q=q α+ 1/2 p P =p+ 2

(

q+ α+ 1/2 p

)2

+t H¯ =H+q+ α+ 1/2

p αα+ 1

r2:

Q=q+ α1/2 p2q2t P =p+ 2q2+t H¯ =H+q αα1

r3:

Q=q

P =p+ 2q2+t H¯ =H+q

α→ −α

ただし,p= 0のときはベックルント変換r1,p2q2t= 0のときはベックルント変換r2が定義できない ことに注意しておく.そこで,ベックルント変換r1, r2が定義できないときのパラメータαの条件を確認して おく.

命題1.2 Painlev´e II型方程式のハミルトン系HIIについて,以下が成り立つ. (1)p= 0ならば,α=1/2である.

(2)p2q2t= 0ならば,α= 1/2である. 証明 (1)p= 0のときは,HIIの第2式より明らか. (2)p2q2t= 0のとき,両辺をtで微分すると

p2q2t= 0 −→ p04qq01 = 0

−→

(

2qp+α+ 1 2

)

4q (

pq2 t 2

)

1 = 0

−→ −2q(p2q2t) +α 1 2 = 0

(10)

となる.よって,p2q2t= 0ならばα= 1/2である. (証明終) この命題は,次節で有理解を構成していくときに使われる.

1.4 Painlev´e II型方程式の有理解

 前節で紹介したベックルント変換を用いて,PIIの有理解を構成していく.その前に,PIIが有理解をもつとき のパラメータαの条件を確認することにする.そのために,いくつかの補題を用意しておく.

補題1.1 PIIの有理解yが有界な点t=cで極をもつとき,極の位数は1で, Res

t=cy dt=±1である. 証明 t=cでの有理解yの展開を

y= c0

(tc)m + c1

(tc)m1 +· · · · (c06= 0, m >0 )

と表し,これをPIIに代入すれば,最低次の項からm= 1, c0=±1が分かる. (証明終) 補題1.2 PIIの有理解y,無限遠t=で次のように展開される.

y=α

t + 2α(α21) t4 + λ5

t5 +· · ·. 証明 t= 1/sと変換すると,PII

s3d2y

dt2 + 2s3dy

dt = 2y3+ y

s +α (1.15)

となる.よって,有理解ys= 0における展開を

y=λ0sn+λ1sn+1+λ2sn+2+· · · (λ06= 0, nZ)

と表し,これを(1.15)に代入すれば,λ0=α, λ1=λ2= 0, λ3= 2α(α21),· · · が得られる. (証明終) これらの補題および定理1.2から,次が成り立つ.

定理1.5 Painlev´e II型方程式PIIが有理解をもつとき,αZである.

命題1.2より,α Zであればベックルント変換r1, r2はいつでも定義できる.従って,ベックルント変換 r1, r2を施していけば, Painlev´e II型方程式の有理解はα= 0での有理解を調べればすべて決定されることに なる.

定理1.6 α= 0のとき, Painlev´e II型方程式PIIの有理解はy= 0のみである. 証明 y= 0以外に有理解が存在するとして,その有理解を

y= P(t)

Q(t) (1.16)

と表す(P(t), Q(t)tの多項式).補題1.2より,α= 0ならばt=でのローラン展開は y= λ5

t5 + λ6

t6 +· · · ·

表 2 q 1 = 0 のときの有理解の生成 q 1 p 1 q 2 p 2 θ 1 θ 2 .. . ... ... ... ... ... 0 3t(t 3 + 2) (t 3 − 1) 2 t(t 6 − 5t 3 − 5)(t3−1)2 − 3t 2 (t 6 − 2t 3 + 10)2(t3−1)(t6−5t3− 5) 12 − 3 0 1 t 2 t 3 − 1t2 − 2t 3 + 12t(t3 − 1) 12 − 2 0 0 t − 1 2t 12 − 1 0 0 t 0 1 2 0 0 1
表 3 ベックルント変換 σ 5 もしくは σ 6 が定義できないときの有理解の生成 ( θ 1 ∈ Z ≥0 + 3/2, θ 2 = 0 ) q 1 p 1 q 2 p 2 θ 1 θ 2 1 2t 0 4t 3 − 14t2 0 32 0 1 t 0 t 0 52 0 3t 2 2t 3 + 1 0 t(4t 6 + 7t 3 + 7)(2t3+ 1)2 0 72 0 2t 3 + 1 t(t 3 + 2) 0 t 9 + 6t 6 + 12t 3 − 1t2(t3+ 2)2 0 92 0 5(t 3
表 5 ベックルント変換 σ 5 もしくは σ 6 が定義できないときの有理解の生成 ( θ 1 ∈ Z ≤0 + 1/2, θ 2 = −1 ) q 1 p 1 q 2 p 2 θ 1 θ 2 − 1 2t 1t 2 4t 3 − 14t2 − 2t 3 + 1t(4t3 + 1) − 12 − 1 − 1 t 2t 2 t − 1 2t − 32 − 1 − 3t 2 2t 3 + 1 12t(t 3 − 1)(2t3+ 1)2 t(4t 6 + 7t 3 + 7)(2t3+ 1)2 − 8t 9 −

参照

関連したドキュメント

Tskhovrebadze, On two-point boundary value problems for systems of higher- order ordinary differential equations with singularities, Georgian Mathematical Journal 1 (1994),

Tskhovrebadze, On two-point boundary value problems for systems of higher- order ordinary differential equations with singularities, Georgian Mathematical Journal 1 (1994),

Shen, “A note on the existence and uniqueness of mild solutions to neutral stochastic partial functional differential equations with non-Lipschitz coefficients,” Computers

The basic idea is that, due to (2), if a Fuchsian system has finite linear monodromy group then the solution to the isomonodromy equations, controlling its deformations, will only

Wong, “Oscillation criteria for second order forced ordinary differential equations with mixed nonlinearities,” Journal of Mathematical Analysis and Applications, vol.

Stavroulakis, Oscillations of first-order delay differential equations in a critical state, Applicable Analysis, 61(1996), 359-377..

A new direct operational inversion method is introduced for solving coupled linear systems of ordinary fractional differential equations.. The solutions so-obtained can be

In the case of constant growth rates and homogeneous measure chains, that is, for ordinary differential equations and ordinary difference equations, the above gap condition (4.4)