ordinarykuma的blog

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

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

Def. Analytic

    We consider a function, f(x). If exists ρ>0 such that ∀x∈C∞(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=0∞a_nxn+r where r is the Indicial Root(指標根). Now, We consider two indicial roots, r_1 and r_2, and we let r_1−r_2≥0 .

Gauss' Hypergeometric ODE

    Consider x(1−x)y″+(c−(a+b+1)x)y′−aby=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)x1−x)y′+−abx1−xy=0. We can note that both P(x) and Q(x) have a Regular Singular Point at x=0.
  2. Solve Indicial Equation
    r(r−1)+P(0)r+Q(0)=0⟹r_1=0 and r_2=1−c

General Solution for (y_1)

    Now, we consider that r_1=0 and r_2=1−c where c∈ℝ∩c>0

  1. Set y1=∑l=0∞alxl+r1=∑l=0∞alxl So, we obtain
    x(1−x)∑l=2∞all(l−1)xl−2+(c−(a+b+1)x)∑l=1∞allxl−1−ab∑l=0∞alxl=
∑l=2∞all(l−1)xl−1−∑l=2∞all(l−1)xl+∑l=1∞callxl−1−∑l=1∞(a+b+1)allxl−∑l=0∞abalxl=∑l=1∞al+1l(l+1)xl−∑l=2∞all(l−1)xl+∑l=0∞cal+1(l+1)xl−∑l=1∞(a+b+1)allxl−∑l=0∞abalxl=2a2x+∑l=2∞al+1l(l+1)xl−∑l=2∞all(l−1)xl+ca1+ca2x+∑l=2∞cal+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))xl⟹al+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=0l−1(a+m)(b+m)(c+m)(1+m))xl

3.推廣

F(1,1,1;x)=1+∑l=1∞(∏m=0l−1(1+m)(1+m)(1+m)(1+m))xl=1+∑l=1∞xl=11−x

for |x|<1

F(1,b,b;x)=1+∑l=1∞(∏m=0l−1(1+m)(b+m)(b+m)(1+m))xl=1+∑l=1∞xl=11−xF(a,1,a;x)=1+∑l=1∞(∏m=0l−1(a+m)(1+m)(a+m)(1+m))xl=1+∑l=1∞xl=11−x
  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=0l−1(−n+m)(b+m)(b+m)(1+m))(−x)l=1+∑l=1∞(∏m=0l−1n−m1+m)xl =1+∑l=1∞n(n−1)...(n−(l−1))(n−l)(n−(l+1))...2*1(1*2*...*(l−1)l)((n−l)(n−(l+1))...2*1)xl⟹1+∑l=1∞n!l!(n−l)!xl=(1+x)n =1+∑_l=1∞n(n−1)...(n−(l−1))(n−l)(n−(l+1))...2\*1(1\*2\*...\*(l−1)l)((n−l)(n−(l+1))...2\*1)xl⟹1+∑_l=1∞n!l!(n−l)!xl=(1+x)n (可以注意到的是,當l>n時,an+1會是0,因此級數收斂,所以F(−n,b,b;−x)=(1+x)n)

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

xF(12,1,32;−x2)=x(1+∑l=1∞(∏m=0l−1(12+m)(1+m)(32+m)(1+m))(−x2)l) =x(1+∑l=1∞12(12+2)......(12+(l−1))(12+1)(12+2)......(12+(l−1))(12+(l−1)+1)(−x2)l) =x(1+∑l=1∞1212+l(−x2)l)=x+∑l=1∞(−1)lx2l+12l+1=tan−1x for |x|<1

xF(12,12,32;x2)=x(1+∑l=1∞(∏m=0l−1(12+m)(12+m)(32+m)(1+m))(x2)l) =x(1+∑l=1∞12l+112l12l1*(2*1)*2*(2*2)*3*(2*3)......(2l−1)(2*l)1*2*...*l*(1*2*...*l)x2l) =x(1+∑l=1∞x2l4l(2l+1)(l!)2)=sin−1x for |x|<1

xF(1,1,2;−x)=x(1+∑l=1∞(∏m=0l−1(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=0l−1(12+m)(1+m)(2+m)(1+m))(x2)l)=2x+∑l=1∞2x2l+12l+1=ln1+x1−x

General Solution for (y_2)

    Now, we consider that r_1=0 and r_2=1−c where c∈ℝ∩c<0

    Set y_2=∑_l=0∞a_lxl+(1−c). So, We have x(1−x)∑l=0∞al(l+1−c)(l−c)xl+(1−c)−2+(c−(a+b+1)x)∑l=0∞al(l+1−c)xl+(1−c)−1−ab∑l=0∞alxl+(1−c)= ∑l=0∞al(l+1−c)xl−c+∑l=0∞al((l+1−c)(l−a−b−c−1)−ab)xl+(1−c) =∑l=0∞al(l+1−c)xl−c+∑l=1∞al−1((1−c)(l−a−b−c−2)−ab)xl−c =a0(1−c)c−c+∑l=1∞(al(l+1−c)+al−1((1−c)(l−a−b−c−2)−ab))xl−c ⟹a0=0,al=∏m=1l−1(a−c+m−1)(b−c+m−1)(1+m)((2−c)+m) ⟹y2=x1−c(1+∑l=1∞(∏m=1l−1(a−c+m−1)(b−c+m−1)(1+m)((2−c)+m))xl)=x1−cF(a−c+1,b−c+1,2−c;x)

#物理