0
大学数学基礎解説
文献あり

RiemannのP関数と超幾何関数で遊んでみようなのだ!

121
0

はじめに

超幾何関数のあのどっから来たのかよく分からない変換公式がどうやって導かれたのか?そのモヤモヤの解決に少しでもつなげるために本稿を書きましたのだ!
すっごくなっがーい記事だけど、ずんだもちでも食べながらゆっくり咀嚼して読んでくれると嬉しいのだ!

RiemannのP関数の定義と定理なのだ!

RiemannのP微分方程式

d2ydx2+[1ααxa+1ββxb+1γγxc]dydx+[αα(ab)(ac)xa+ββ(bc)(b1)xb+γγ(ca)(cb)xc]y(xa)(xb)(xc)=0
ここで、α+α+β+β+γ+γ=1
また、この微分方程式の解を次のように書きRiemannのP関数というのだ!
y=P{abcαβγxαβγ}

RiemannのP関数の特性指数

d2ydx2+[1ααxa+1ββxb+1γγxc]dydx+[αα(ab)(ac)xa+ββ(bc)(ba)xb+γγ(ca)(cb)xc]y(xa)(xb)(xc)=0
x=a近傍の級数解をu=n=0an(xa)n+sとおくのだ。また下記のことが成り立つので、RiemannのP微分方程式に(xa)2を両辺にかけて(xa)sの係数を考えるのだ。{dudx=n=0(n+s)(xa)n+sd2udx2=n=0(n+s)(n+s1)(xa)n+s
s(s1)+(1αα)s+αα=(sα)(sα)=0
このことから、x=aでの特性指数はα,αであることが分かる。x=b,x=cでも同様なのだ。このことからRiemannのP関数は確定特異点と特性指数を書いたものだとわかるのだ。

メビウス変換とRiemannのP関数

x=rx+spx+qとすると次のように書き直せるのだ。
x=rp(px+q)+sqrppx+q=rp+sqrppx+q=rp+1p(sqrp)x+qp
xで行っていることを細かく分けるのだ。

  1. x1x+qp()
  2. x21x1()
  3. x31p(sqrp)x2()
  4. x4x3+rp()
  1. 平行移動:x1=x+δとすると、v(x1)=u(x1δ),ddx=ddx1,d2dx2=d2dx12よりRiemannのPBe分方程式は次のように書けるのだ。d2vdx12+[1ααx1(a+δ)+1ββx1(b+δ)+1γγx1(c+δ)]dvdx1+[αα(ab)(ac)x1(a+δ)+ββ(bc)(ba)x1(b+δ)+γγ(ca)(cb)x1(c+δ)]v{x1(a+δ)}{x1(b+δ)}{x1(c+δ)}=0より、平行移動による解は次のように書けるのだ。P{abcαβγxαβγ}P{a+δb+δc+δαβγx+δαβγ}
  2. 拡大・回転:x2=pxとするとv(x2)=u(x2p),ddx=pddx2,d2dx2=p2d2dx22より、RiemannのP微分方程式は次のように書き直せるのだ。p2d2vdx22+p[1ααx2pa+1ββx2pb+1γγx2pc]dvdx2+[αα(ab)(ac)x2pa+ββ(bc)(ba)x2pb+γγ(ca)(cb)x2px]v{x2pa}{x2pb}{x2pc}=0d2vdx22+[1ααx2pa+1ββx2pb+1γγx2pc]dvdx2+[αα(papb)(papc)x2pa+ββ(pbpc)(pbpa)x2pb+1γγx2pc]v(x2pa)(x2pb)(x2pc)=0より回転・拡大によるRiemannのP関数は次のように変換されるのだ。P{abcαβγxαβγ}P{papbpcαβγxαβγ}
  3. 反転の場合:x3=1xとすると次のように書けるのだ。v(x3)=u(1x3),ddx=x32ddx3,d2dx2=x34d2dx32+2x33ddx3この事を用いるとRiemannのP微分方程式は下記のように変換されることが分かるのだ。x34d2vdx32+2x33dvdx3x32[1αα1x3a+1ββ1x3b+1γγ1x3c]dvdx3+[αα(ab)(ac)1x3a+ββ(bc)(ba)1x3b+γγ(ca)(cb)1x3c]v(1x3a)(1x3b)(1x3c)=0x3d2vdx32+[21αα1ax31ββ1bx31γγ1cx3]dvdx3+[αα(ab)(ac)1ax3+ββ(bc)(ba)1bx3+γγ(ca)(cb)1cx3]x3v(1ax3)(1bx3)(1cx3)=0ここでA=1αα,B=1ββ,C=1γγ,D=αα,E=ββ,F=γγとおくと次のように書けるのだ。A+B+C=2x3d2vdx32+[A+B+CA1ax3B1bx3C1cx3]dvdx3+[D(ab)(ac)1ax3+E(bc)(ba)1bx3+F(c1)(cb)1cx3]x3v(1ax3)(1bx3)(1cx3)=x3d2vdx32x3[aA1ax3+bB1bx3+cC1cx3]dvdx3+[D(ab)(ac)1ax3+E(bc)(ba)1bx3+F(c1)(cb)1cx3]x3v(1ax3)(1bx3)(1cx3)=0d2vdx32+[Ax31a+Bx31b+Cx31c]dvdx3+[D(1a1b)(1a1c)x31a+E(1b1c)(1b1a)x31b+F(1c1a)(1c1b)x31c]v(x31a)(x31b)(x31c)=0
    この計算より、反転によってRiemannのP関数は下記のように変形されることが示されたのだ。P{abcαβγxαβγ}P{1a1b1cαβγxαβγ}
  4. まとめると、x=rx+spx+qなる変換についてRiemannのP関数は下記のようになるのだ。P{abcαβγxαβγ}P{a+qpb+qpc+qpαβγxαβγ}P{1a+qp1b+qp1c+qpαβγxαβγ}P{1p(sqrp)a+qp1p(sqrp)b+qp1p(sqrp)c+qpαβγxαβγ}P{rp+1p(sqrp)a+qprp+1p(sqrp)b+qprp+1p(sqrp)c+qpαβγrx+spx+qαβγ}よって、x=T(x)=rx+spx+qとすれば下記のように変換されたことが分かったのだ。P{abcαβγxαβγ}P{T(a)T(b)T(c)αβγT(x)αβγ}

(xaxb)μ(xcxb)νP{abcαβγxαβγ}=P{abcα+μβμνγ+νxα+μβμνγ+ν}
が成り立つ。
実際、x=aでの近傍の解で特性指数αの解をu=n=0an(xa)n+αとおくのだ。またv=(xaxb)μ(xcxb)νuと書くのだ。
v=(xaxb)μ(xcxb)νn=0an(xa)n+α=(xc)ν(xb)μ+νn=0an(xa)n+α+μ
また次の様に書けるのだ。
(xc)ν(xb)μ+ν=(ac)ν(ab)μ+νm=0(νm)(xa)mn=0(1)n(xa)n=(ac)ν(ab)μ+νm=0(xa)mn=0m(1)n(νmn)このことから、vがRiemannのP微分方程式とすれば、x=aの近傍で、特性指数α+μの解であることが分かるのだ。
同様にして、vがRiemannのP微分方程式であるとすれば
v=P{abcα+μβμνγ+νxα+μβμνγ+ν}
またvは次のRiemannのP微分方程式を満たすことを示すのだ。
d2vdx2+[1αα2μxa+1ββ+2(μ+ν)xb+1γγ2νxc]dvdx+[(α+μ)(α+μ)(ab)(ac)xa+(βμν)(βμν)(bc)(b1)xb+(γ+ν)(γ+ν)(ca)(cb)xc]v(xa)(xb)(xc)=0
これを満たすことを確認するために下記の事実を用いるのだ。
dvdx=ddx((xa)μ(xc)ν(xb)μ+νu)=(xa)μ(xc)ν(xb)μ+ν(dudx+(μxaμ+νxb+νxc)u)
d2vdx2=(xa)μ(xc)ν(xb)μ+ν(d2udx2+(μxaμ+νxb+νxc)dudx)+μ(xa)μ1(xc)ν(xb)μ+ν(dudx+(μ1xaμ+νxb+νxc)u)(μ+ν)(xa)μ(xc)ν(xb)μ+ν+1(dudx+(μxaμ+ν+1xb+νxc)u)+ν(xa)μ(xc)ν1(xb)μ+ν(dudx+(μxaμ+νxb+ν1xc)u)=(xa)μ(xc)ν(xb)μ+νd2udx2+2(xa)μ(xc)ν(xb)μ+ν(μxaμ+νxb+νxc)dudx+(xa)μ(xc)ν(xb)μ+ν(μ(μ1)(xa)2μ(μ+ν)(xa)(xb)+μν(xa)(xc)μ(μ+ν)(xa)(xb)+(μ+ν)(μ+ν+1)(xb)2ν(μ+ν)(xb)(xc)+μν(xa)(xc)ν(μ+ν)(xb)(xc)+ν(ν1)(xc)2)u=(xa)μ(xc)ν(xb)μ+νd2udx2+2(xa)μ(xc)ν(xb)μ+ν(μxaμ+νxb+νxc)dudx+(xa)μ(xc)ν(xb)μ+ν(μ(μ1)(xa)2+(μ+ν)(μ+ν+1)(xb)2+ν(ν1)(xc)22(μ(μ+ν)(xa)(xb)μν(xa)(xc)+ν(μ+ν)(xb)(xc)))u
これらの事実を用いてRiemannのP微分方程式に代入すると下記の様になるのだ。
d2udx2+2(μxaμ+νxb+νxc)dudx+(μ(μ1)(xa)2+(μ+ν)(μ+ν+1)(xb)2+ν(ν1)(xc)22(μ(μ+ν)(xa)(xb)μν(xa)(xc)+ν(μ+ν)(xb)(xc)))u+[1αα2μxa+1ββ+2(μ+ν)xb+1γγ2νxc]×(dudx+(μxaμ+νxb+νxc)u)+[(α+μ)(α+μ)(ab)(ac)xa+(βμν)(βμν)(bc)(b1)xb+(γ+ν)(γ+ν)(ca)(cb)xc]u(xa)(xb)(xc)=d2udx2+[1ααxa+1ββxb+1γγxc]dudx+(μ(μ1)(xa)2+(μ+ν)(μ+ν+1)(xb)2+ν(ν1)(xc)22(μ(μ+ν)(xa)(xb)μν(xa)(xc)+ν(μ+ν)(xb)(xc)))u+[1αα2μxa+1ββ+2(μ+ν)xb+1γγ2νxc](μxaμ+νxb+νxc)u+[(α+μ)(α+μ)(ab)(ac)xa+(βμν)(βμν)(bc)(b1)xb+(γ+ν)(γ+ν)(ca)(cb)xc]u(xa)(xb)(xc)
よって次の事を示せばよいのだ。
(μ(μ1)(xa)2+(μ+ν)(μ+ν+1)(xb)2+ν(ν1)(xc)22(μ(μ+ν)(xa)(xb)μν(xa)(xc)+ν(μ+ν)(xb)(xc)))u+[1αα2μxa+1ββ+2(μ+ν)xb+1γγ2νxc](μxaμ+νxb+νxc)u+[μ(α+α+μ)(ab)(ac)xa+(μ+ν)(β+βμν)(bc)(ba)xb+ν(γ+γ+ν)(ca)(cb)xc]u(xa)(xb)(xc)=0
μ(μ1)(xb)2(xc)2+(μ+ν)(μ+ν+1)(xa)2(xc)2+ν(ν1)(xa)2(xb)22(μ(μ+ν)(xc)μν(xb)+ν(μ+ν)(xa))(xa)(xb)(xc)+((1αα2μ)(xb)(xc)+(1ββ+2μ+2ν)(xa)(xc)+(1γγ2ν)(xa)(xb))×((μa(μ+ν)b+νc)x+μbc(μ+ν)ca+νab)+μ(α+α+μ)(ab)(ac)(xb)(xc)(μ+ν)(β+βμν)(bc)(ba)(xa)(xc)+ν(γ+γ+ν)(ca)(cb)(xa)(xb)=0
白状するとこれを手計算で行うのはあきらめたのだ(´;ω;`)。
ただ、どうやってこの式が成り立つかを確認したかは伝えておくのだ。それは下記の様なのだ。

  1. pythonによるコーディング:from sympy import symbols, expand, collect

    # 定数と変数の定義
    x, a, b, c, mu, nu = symbols('x a b c mu nu')
    alpha, alpha, beta, beta, gamma, gamma_ = symbols('alpha alpha_ beta beta_ gamma gamma')

    # 与えられた式
    expr = (mu(mu-1)(x-b)2*(x-c)2
            + (mu+nu)(mu+nu+1)(x-a)2*(x-c)2
            + nu(nu-1)(x-a)2*(x-b)2
            - 2(mu(mu+nu)(x-c) - munu(x-b) + nu(mu+nu)(x-a))(x-a)(x-b)(x-c)
            + ((1-alpha-alpha
    - 2mu)(x-b)(x-c)
               + (1-beta-beta_ + 2
    mu + 2nu)(x-a)(x-c)
               + (1-gamma-gamma_ - 2
    nu)(x-a)(x-b))
            * ((mua - (mu+nu)b + nuc)x + mubc - (mu+nu)ca + nuab)
            + mu(alpha + alpha_ + mu)(a-b)(a-c)(x-b)(x-c)
            - (mu+nu)
    (beta + beta_ - mu - nu)(b-c)(b-a)(x-a)(x-c)
            + nu(gamma + gamma_ + nu)(c-a)(c-b)(x-a)*(x-b))

    # 展開と整理
    expanded_expr = expand(expr)

    # xの多項式として整理
    polynomial = collect(expanded_expr, x)

    print(polynomial)
  2. 出力された結果をchatgptに入力して整形
  3. α+α+β+β+γ+γ=1であることを用いてなんやかんやするのだ!

ただし、4次の項が0である事を示すのは容易なのだ!
実は仕事中に気づいた方法があったので紹介するのだ(2024/07/22追記)
結論から述べると、先の式が3次多項式であることを用いてx=0,a,b,cで0になることを示すことでxによらず恒等的に0になる事を示す。
  1. x=aの場合:μ(μ1)(ab)2(ac)2+(1αα2μ)μ(ab)2(ac)2+μ(α+α+μ)(ab)2(ac)2=0
  2. x=cの場合:a,c0
  3. x=bの場合:(μ+ν)(μ+ν+1)(ba)2(bc)2(μ+ν)(1ββ+2μ+2ν)(ba)2(bc)2(μ+ν)(β+βμν)(ba)2(bc)2=0
  4. x=0の場合:μ(μ1)b2c2+(μ+ν)(μ+ν+1)a2c2+ν(ν1)a2b22μ(μ+ν)abc2+2μνab2c2ν(μ+ν)a2bc+((1αα2μ)bc+(1ββ+2μ+2ν)ac+(1γγ2ν)ab)(μbc(μ+ν)ca+νab)+μbc(α+α+μ)(ab)(ac)(μ+ν)(β+βμν)ac(ba)(bc)+νab(γ+γ+ν)(ca)(cb)=0

4.の計算は面倒なので、さらに分解して考えるのだ。4.の式はa,b,c4次の斎次方程式であることは見たら分かるのだ。またa=0とした場合、これはb2c2の項のみが残り、3.から与えられた0であることが分かるのだ。同様にして、c2a2,a2b2の項も0であることが分かるので、結局a2bc,ab2c,abc2の項にのみ着目すればよいのだ。
  1. abc2の係数について:(abc2)=ν(α+α+β+β+γ+γ1)=0
  2. a,cの対称性より(a2bc)=0
  3. ab2cの係数について:(ab2c)=(μ+ν)(α+α+β+β+γ+γ1)=0
  4. 以上より、証明が完了したのだ!

Gaussの超幾何微分方程式

皆が大好きなGauss超幾何微分方程式は次のようなのだ!

z(1z)d2udz2+{c(a+b+1)z}dudxabu=0

実はこの超幾何微分方程式は実はRiemannの微分方程式の一つなのだ!
実際、RiemannのP微分方程式を下記の様に計算していくことで分かるのだ。
d2udz2+(1ααza+1ββzb+1γγzc)dudz+(αα(ab)(ac)za+ββ(bc)(ba)zb+γγ(ca)(cb)zc)u(za)(zb)(zc)d2udz2+(1ααza+1ββzb)dudz+(αα(ab)za+ββ(ba)zb+γγ)u(za)(zb)d2udz2+(1ααz+1ββz1)dudz+(ααz+ββz1+γγ)uz(z1)
また,超幾何微分方程式の両辺をz(1z)で割る事で下記のように書けるのだ。
d2udz2+(cz+a+bc+1z1)dudz+abz(z1)u=0
これらを比較すると超幾何微分方程式の解が次の様に書けることが分かるのだ。
u=P{0100az1ccabb}=F(a,b;c;z)

T(z)=rz+spz+qとすると
v(z)=u(T(z))=P{sqr+sp+qrp00aT(z)1ccabb}
w(z)=(zsq)μ(zrp)ν(zr+sp+q)μ+νv(z)=P{sqr+sp+qrpμμνa+νT(z)1c+μcabμnub+ν}
(zsq)μ(zrp)ν(zr+sp+q)μ+νF(a,b;c;rz+spz+q)=P{sqr+sp+qrpμμνa+νT(z)1c+μcabμνb+ν}
ここで次の様に定めると
{μ=c1ν=b
次のような結果が得られるのだ。
(zsq)c1(zrp)b(zr+sp+q)cb1F(a,b;c;rz+spz+q)=P{sqr+sp+qrpc1bc+1abrz+spz+q01a0}
ここで次の事実を用いるのだ。
T(z)=rz+spz+qとし、T(a)=0,T(b)=1,T(c)=となるようになるp,q,r,sは次の様に書けるのだ。
{p=babcq=babccr=1s=a
そこで、a=sq,b=rp,c=r+sp+qとおき、次のような変換を考えるのだ。
U(z)=(bc)(za)(ba)(z+c)=(rpr+sp+q)(zsq)(rpsq)(zr+sp+q)=p(qrps)(qzs)(qrps)((p+q)z(r+s))=p(qzs)(p+q)z(r+s)
このUを用いると最終的に次のような結論を得るのだ。
(U(z)sq)c1(U(z)rp)b(U(z)r+sp+q)cb1F(a,b;,c;U(rz+spz+q))=P{01001aU(rz+spz+q)c1abbc+1}=F(1a,bc+1;2c;U(rz+spz+q))
長くなったのでまとめるのだ!

{T(z)=rz+spz+qU(z)=pqzs(p+q)z(r+s)(U(z)T(0))c1(U(z)T())b(U(z)T(1))cb1F(a,b;c;U(T(z)))=F(1a,bc;2c;U(T(z)))

参考文献

[1]
WHITTAKER AND WATSON, MODERN ANALYSIS, Cambridge University Press, 1935
投稿日:2024721
更新日:2024722
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. RiemannのP関数の定義と定理なのだ!
  3. Gaussの超幾何微分方程式
  4. 参考文献