ordinarykuma的blog

高斯超幾何函數(Gauss' Hypergeometric)

  今天來分享暑假學到的數學,順便來練習LaTex輸入公式(我之前都是使用LibreOffice的語法來輸入公式,接著再用AI幫我轉換,但這次想來點不一樣的)。今天的這一個ODE的解法是來自Frobenius展開。那麼,就來開始上課吧。(順帶一提,這學期應數A+實在是很開心,本來還想說要完蛋了,畢竟段考的分數實在不可口)

Def. Analytic

    We consider a function, f(x). If exists ρ>0 such that xC(aρ,a+ρ), then we say that f(x) is Analytic at x=a.

Thm. Analytic

    If f(x) is analytic at x=a, f(x) can be expanded to be Taylor's Polynomial.

Frobenius Series Solution

    We consider an ODE: x2y+xP(x)y+Q(x)y=0 .

    Frobenius discovered that the series solutions of an ODE are usually like _n=0a_nxn+r where r is the Indicial Root(指標根). Now, We consider two indicial roots, r_1 and r_2, and we let r_1r_20 .

Gauss' Hypergeometric ODE

    Consider x(1x)y+(c(a+b+1)x)yaby=0

Indicial Root

  1. We rewrite the ODE first.(寫成Frobenius的標準形式) So, we have
    x2y+xP(x)y+Q(x)y=x2y+x(c(a+b+1)x1x)y+abx1xy=0. We can note that both P(x) and Q(x) have a Regular Singular Point at x=0.
  2. Solve Indicial Equation
    r(r1)+P(0)r+Q(0)=0r_1=0 and r_2=1c

General Solution for (y_1)

    Now, we consider that r_1=0 and r_2=1c where cc>0

  1. Set y1=l=0alxl+r1=l=0alxl So, we obtain
    x(1x)l=2all(l1)xl2+(c(a+b+1)x)l=1allxl1abl=0alxl=
l=2all(l1)xl1l=2all(l1)xl+l=1callxl1l=1(a+b+1)allxll=0abalxl=l=1al+1l(l+1)xll=2all(l1)xl+l=0cal+1(l+1)xll=1(a+b+1)allxll=0abalxl=2a2x+l=2al+1l(l+1)xll=2all(l1)xl+ca1+ca2x+l=2cal+1(l+1)xl(a+b+1)a1x=(aba0+ca1)+((a+b+1+ab)a1+2(c+1)a2)x+l=2(al+1(l+1)(l+c)al(l+a)(l+b))xlal+1=m=0l(a+m)(b+m)(c+m)(1+m)a0=1(l+1)!m=0l(a+m)(b+m)(c+m)
  1. We take a_0=1 y_1=F(a,b,c;x)=1+_l=1(_m=0l1(a+m)(b+m)(c+m)(1+m))xl

3.推廣

F(1,1,1;x)=1+l=1(m=0l1(1+m)(1+m)(1+m)(1+m))xl=1+l=1xl=11x

for |x|<1

F(1,b,b;x)=1+l=1(m=0l1(1+m)(b+m)(b+m)(1+m))xl=1+l=1xl=11xF(a,1,a;x)=1+l=1(m=0l1(a+m)(1+m)(a+m)(1+m))xl=1+l=1xl=11x
  1. By Ratio Test, we have liml|al+1al|=liml|x(a+l)(b+l)(c+l)(1+l)|=liml|x|<1

5.Common Expansion F(n,b,b;x)=1+l=1(m=0l1(n+m)(b+m)(b+m)(1+m))(x)l=1+l=1(m=0l1nm1+m)xl =1+l=1n(n1)...(n(l1))(nl)(n(l+1))...2*1(1*2*...*(l1)l)((nl)(n(l+1))...2*1)xl1+l=1n!l!(nl)!xl=(1+x)n =1+_l=1n(n1)...(n(l1))(nl)(n(l+1))...2\*1(1\*2\*...\*(l1)l)((nl)(n(l+1))...2\*1)xl1+_l=1n!l!(nl)!xl=(1+x)n (可以注意到的是,當l>n時,an+1會是0,因此級數收斂,所以F(n,b,b;x)=(1+x)n

1nxF(1n,1,2;x)=1nx(1+l=1(m=0l1(1n+m)(1+m)(2+m)(1+m))xl) =1nx(1+l=1(1)l(m=0l1(n1)m2+m)xl) =1nx(1+l=1(1)l(m=0l1(n1)m2+m)xl) 1nx(1+l=1n(1)l(n1)!(l+1)!(n(l+1))!xl) =1nx+l=1n(1)l+1n!(l+1)!(n(l+1))!xl+1=1nx+l=2n(1)l(nl)xl=(1x)n (可以注意到的是,當l>n時,an+1會是0,因此級數收斂,所以1nxF(1n,1,2;x)=(1x)n

xF(12,1,32;x2)=x(1+l=1(m=0l1(12+m)(1+m)(32+m)(1+m))(x2)l) =x(1+l=112(12+2)......(12+(l1))(12+1)(12+2)......(12+(l1))(12+(l1)+1)(x2)l) =x(1+l=11212+l(x2)l)=x+l=1(1)lx2l+12l+1=tan1x for |x|<1

xF(12,12,32;x2)=x(1+l=1(m=0l1(12+m)(12+m)(32+m)(1+m))(x2)l) =x(1+l=112l+112l12l1*(2*1)*2*(2*2)*3*(2*3)......(2l1)(2*l)1*2*...*l*(1*2*...*l)x2l) =x(1+l=1x2l4l(2l+1)(l!)2)=sin1x for |x|<1

xF(1,1,2;x)=x(1+l=1(m=0l1(1+m)(1+m)(2+m)(1+m))(x)l) =x(1+l=1(1*2*......*l2*3*...*l*(l+1))(x)l)=x+l=1(1)lx2l+1l+1=ln(1+x)

2xF(12,1,32;x2)=2x(1+l=1(m=0l1(12+m)(1+m)(2+m)(1+m))(x2)l)=2x+l=12x2l+12l+1=ln1+x1x

General Solution for (y_2)

    Now, we consider that r_1=0 and r_2=1c where cc<0

    Set y_2=_l=0a_lxl+(1c). So, We have x(1x)l=0al(l+1c)(lc)xl+(1c)2+(c(a+b+1)x)l=0al(l+1c)xl+(1c)1abl=0alxl+(1c)= l=0al(l+1c)xlc+l=0al((l+1c)(labc1)ab)xl+(1c) =l=0al(l+1c)xlc+l=1al1((1c)(labc2)ab)xlc =a0(1c)cc+l=1(al(l+1c)+al1((1c)(labc2)ab))xlc a0=0,al=m=1l1(ac+m1)(bc+m1)(1+m)((2c)+m) y2=x1c(1+l=1(m=1l1(ac+m1)(bc+m1)(1+m)((2c)+m))xl)=x1cF(ac+1,bc+1,2c;x)

#物理