《特殊函数概论》Chapter 3习题解答

众所周知,王竹溪、郭敦仁所著的《特殊函数概论》是一本对初学特殊函数的同学非常友好的书,特别是对我这种英语不好的人来讲,不用一边翻字典一边看Whittaker&Waston了

但是据我所知,特殊函数概论应该是没有完整的习题解析(b站有大佬上传过第一章的解析)

笔者在去年写过第三章部分习题的解析,这段时间我补全了全部解答(花了好几节体育课)并将其码成文章发出来,希望能给学习特殊函数的您一点帮助~

为方便起见,创作过程中我将在不改变题意的前提下对原题进行一定的改写

现在让我们开始吧!

正文

1.证明在 \sum_{1\leq m\leq k}(a_m-b_m)=0 的条件下 \prod_{n=1}^\infty\left\{\frac{(n-a_1)\cdots(n-a_k)}{(n-b_1)\cdots(n-b_k)}\right\}=\prod_{1\leq m\leq k}\frac{\Gamma(1-b_m)}{\Gamma(1-a_m)}

依 \mathbf{Weierstrass} 公式有:

\frac1{\Gamma(1+x)}=\mathrm e^{\gamma x}\prod_{n=1}^\infty\left(\frac{n+x}{n}\right)\mathrm e^{-\frac xn}\\ {\Gamma(1+x)}=\mathrm e^{-\gamma x}\prod_{n= 1}^\infty\left(\frac{n}{n+x}\right)\mathrm e^{\frac xn}\\

所以

1mkΓ(1bm)Γ(1am)=n=1{(na1)(nak)(nb1)(nbk)}exp{1mk(ambm)(1n+γ)}=n=1{(na1)(nak)(nb1)(nbk)}1mkΓ(1bm)Γ(1am)=n=1{(na1)(nak)(nb1)(nbk)}exp{1mk(ambm)(1n+γ)}=n=1{(na1)(nak)(nb1)(nbk)}(1)(2)\\

 

2.证明: \frac{(a+b)!}{a!b!}=\prod^\infty_{s=1}\frac{(s+a)(s+b)}{s(s+a+b)}

(a+b)!a!b!=Γ(a+b+1)Γ(a+1)Γ(b+1)=eγ(a+b)s=1ss+a+bea+bseγas=1s+aseaseγbs=1s+bsebs=s=1(s+a)(s+b)s(s+a+b)(a+b)!a!b!=Γ(a+b+1)Γ(a+1)Γ(b+1)=eγ(a+b)s=1ss+a+bea+bseγas=1s+aseaseγbs=1s+bsebs=s=1(s+a)(s+b)s(s+a+b)(3)(4)(5)\\

 

3.证明:若 \omega:=\mathrm e^{2\pi\mathrm i/n},n\in\mathbb{Z}_{\geq 1} ,有 x^n\prod^\infty_{k=1}\left(1-\frac{x^n}{k^n}\right)=-\prod^{n-1}_{m=0}\frac1{\Gamma(-\omega^mx)}

依 \mathbf{Weierstrass} 公式有:

\frac1{\Gamma(-\omega^mx)}=-\omega^mx\mathrm e^{-\omega^mx}\prod^\infty_{k=1}\left(1-\frac {\omega^mx}{k}\right)\mathrm e^{{\omega^mx}/{k}}\\

所以

n1m=01Γ(ωmx)=(1)n+1ωn(n1)/2exp{xn1j=0ωj}×n1m=0k=1(1ωmxk)eωmx/k=(1)n+1eπ(n1)ik=1n1m=0(1ωmxk)eωmx/k=k=1(1xnkn)n1m=01Γ(ωmx)=(1)n+1ωn(n1)/2exp{xn1j=0ωj}×n1m=0k=1(1ωmxk)eωmx/k=(1)n+1eπ(n1)ik=1n1m=0(1ωmxk)eωmx/k=k=1(1xnkn)(6)(7)(8)(9)\\

 

4.证明 \Gamma(z)=\int^\infty_0t^{z-1}\left\{\mathrm e^{-t}+\sum^k_{j=0}(-1)^{k+1}\frac{t^k}{k!}\right\}\mathrm dt ,其中 k\in\mathbb{N}^+,-k-1<\Re(z)<-k

由 \mathbf{Gamma} 函数的定义 \Gamma(z)=\int^\infty_0t^{z-1}\mathrm e^{-t}\mathrm dt=\int^1_0t^{z-1}\mathrm e^{-t}\mathrm dt+\int^\infty_1t^{z-1}\mathrm e^{-t}\mathrm dt

\Re(z)<0 时,第一项积分不收敛,考虑把 \mathrm e^{-t} 泰勒级数的前 k 项放进去,使之变为

\int^1_0t^{z-1}\mathrm e^{-t}\mathrm dt=\int^1_0t^{z-1}\color{blue}{\left[\mathrm e^{-t}+\sum_{j=0}^k(-1)^{k+1}\frac{t^k}{k!}\right]}\mathrm dt+\color{red}{\sum^k_{j=0}\frac{(-1)^{j}}{j!(t+j)}}\\

蓝色部分 =\sum_{j\geq k+1}\frac{(-1)^j}{j!}t^{j+z-1}=\frac{(-1)^{k+1}}{(k+1)!}t^{k+z}+o(x^{k+z})

此时收敛域变为 \Re(z)>-k-1

红色部分 =\sum^k_{j=0}\frac{(-1)^{j+1}}{j!}\int^\infty_1t^{j+z-1}\mathrm dt=\int^\infty_1t^{z-1}\sum^k_{j=0}\frac{(-1)^{j+1}}{j!}t^j\mathrm dt

收敛域变为 \Re(z)<-k

所以 \Gamma(z)=\int^\infty_0t^{z-1}\left\{\mathrm e^{-t}+\sum^k_{j=0}(-1)^{k+1}\frac{t^k}{k!}\right\}\mathrm dt

 

6.证明 \mathbf{Binet} 公式的另一形式,对于 \Re(z)>0\log\Gamma(z)=\left(z-\frac12\right)\log z-z+\frac12\log(2\pi)+\int^\infty_0\left(\frac12-\frac1t+\frac1{\mathrm e^{t}-1}\right)\frac{\mathrm e^{-zt}}{t}\mathrm dt\\

 (z12)logzz+12log(2π)+0(121t+1et1)ezttdt=(z12)logzz+12log(2π)+0(111et+1et1)ezttdt0(121t+1et1)ezttdt=(z12)logzz+12log(2π)0(121t+1et1)ezttdt=logΓ(z) (z12)logzz+12log(2π)+0(121t+1et1)ezttdt=(z12)logzz+12log(2π)+0(111et+1et1)ezttdt0(121t+1et1)ezttdt=(z12)logzz+12log(2π)0(121t+1et1)ezttdt=logΓ(z)(10)(11)(12)(13)(14)

 

7. 证明:\int^\infty_0\left(\frac12-\frac1t+\frac1{\mathrm e^t-1}\right)\frac{\mathrm e^{-zt}}{t}\mathrm dt=1-\log(2\pi)

前一问的简单推论罢了

 

8.证明: \gamma=\int^\infty_0\left\{\frac1{1+t}-\mathrm e^{-t}\right\}\frac{\mathrm dt}{t}

 0{11+tet}dtt=limr0+{rdtt(1+t)rettdt}=limr0+{logr[logrγ+O(r)]}=γ 0{11+tet}dtt=limr0+{rdtt(1+t)rettdt}=limr0+{logr[logrγ+O(r)]}=γ(15)(16)(17)(18)\\

 

9.由 \log\Gamma(z)=\int^\infty_0\left\{z-1-\frac{1-\mathrm e^{-(z-1)t}}{1-\mathrm e^{-t}}\right\}\frac{\mathrm e^{-t}}{t}\mathrm dt\quad,(\Re(z)>0) 证明在 0<\Re(z)<1 的条件下 2\log\Gamma(z)=\log\pi-\log\sin(\pi z)+\int^\infty_0\left\{\frac{\sinh(1/2-z)t}{\sinh(t/2)}-(1-2z)\mathrm e^{-t}\right\}\frac{\mathrm dt}{t}

2logΓ(z)=logΓ(z)+logπlogsin(πz)logΓ(1z)=logπlogsin(πz)+0{z11e(z1)t1et+z+1ezt1et}ettdt=logπlogsin(πz)+0{2z1+e(z1)tezt1et}ettdt=logπlogsin(πz)+0{sinh(1/2z)tsinh(t/2)(12z)et}dtt2logΓ(z)=logΓ(z)+logπlogsin(πz)logΓ(1z)=logπlogsin(πz)+0{z11e(z1)t1et+z+1ezt1et}ettdt=logπlogsin(πz)+0{2z1+e(z1)tezt1et}ettdt=logπlogsin(πz)+0{sinh(1/2z)tsinh(t/2)(12z)et}dtt(19)(20)(21)(22)(23)

 

10&11.证明:
\log\Gamma(x)=\left(\frac12-x\right)(\gamma+\log 2)+(1-x)\log\pi-\frac12\log\sin(\pi x)+\frac1\pi\sum^\infty_{n=1}\frac{\log n}{n}\sin(2\pi nx)

将 \log \Gamma \left( x \right) 在 x\in \left( 0,1 \right) 上作傅里叶级数展开:

\log \Gamma \left( x \right) =\frac{a_0}{2}+\sum_{n=1}^\infty{a_n\cos \left( 2\pi nx \right) +b_n\sin \left( 2\pi nx \right)} \quad , x\in \left( 0,1 \right)

a0=210logΓ(x)dx=10[logΓ(x)+logΓ(1x)]dx=10logπsin(πx)dx=lnπ10lnsin(πx)dx=ln(2π)a0=210logΓ(x)dx=10[logΓ(x)+logΓ(1x)]dx=10logπsin(πx)dx=lnπ10lnsin(πx)dx=ln(2π)(24)(25)

an=210logΓ(x)cos(2πnx)dx=10[logΓ(x)+logΓ(1x)]cos(2πnx)dx=10log2π2sin(πx)cos(2πnx)dx=10log(2sinπx)cos(2πnx)dx=k=11k10cos(2πkx)cos(2πnx)dx=12k=11k[k=n]=12nan=210logΓ(x)cos(2πnx)dx=10[logΓ(x)+logΓ(1x)]cos(2πnx)dx=10log2π2sin(πx)cos(2πnx)dx=10log(2sinπx)cos(2πnx)dx=k=11k10cos(2πkx)cos(2πnx)dx=12k=11k[k=n]=12n(26)(27)(28)

bn=210logΓ(x)sin(2πnx)dx=1πn10logΓ(x)d[1cos(2πnx)]=1πn10[cos(2πnx)1]ψ(x)dx=1πn10[cos(2πnx)1](γ+10tx11t1dt)dx=γπn+1πn101t1{10[cos(2πnx)1](tx11)dx}dt=γπn+1πn10[1t1+1tlogtlog(t)2+(2πn)21tlogt]dt=γ+log(2πn)πnbn=210logΓ(x)sin(2πnx)dx=1πn10logΓ(x)d[1cos(2πnx)]=1πn10[cos(2πnx)1]ψ(x)dx=1πn10[cos(2πnx)1](γ+10tx11t1dt)dx=γπn+1πn101t1{10[cos(2πnx)1](tx11)dx}dt=γπn+1πn10[1t1+1tlogtlog(t)2+(2πn)21tlogt]dt=γ+log(2πn)πn

所以

logΓ(x)=a02+n=1ancos(2πnx)+bnsin(2πnx)=12log(2π)+n=1cos(2πnx)2n+n=1γ+log(2πn)πnsin(2πnx)=12log(2πn)12log(2sinπx)+n=1γ+log(2π)πnsin(2πnx)+1πn=1lognnsin(2πnx)=(12x)(γ+log2)+(1x)logπ12logsin(πx)+1πn=1lognnsin(2πnx)logΓ(x)=a02+n=1ancos(2πnx)+bnsin(2πnx)=12log(2π)+n=1cos(2πnx)2n+n=1γ+log(2πn)πnsin(2πnx)=12log(2πn)12log(2sinπx)+n=1γ+log(2π)πnsin(2πnx)+1πn=1lognnsin(2πnx)=(12x)(γ+log2)+(1x)logπ12logsin(πx)+1πn=1lognnsin(2πnx)(29)(30)(31)(32)(33)

 

12.在 \lambda>0,x>0,-\pi/2<\alpha<\pi/2 的条件下,证明:
\int^\infty_0t^{x-1}\mathrm e^{-\lambda\cos \alpha t}cossincossin(\lambda t\sin \alpha)=\lambda^{-x}\Gamma(x)cossincossin(\alpha x)

 0tx1eλcosαteiλtsinαdt=0tx1e(isinαcosα)λtdtλtt=1λx0tx1exp(eiαt)dteiαtt=eiαxλxeiα0tx1etdt=eiαλxΓ(x) 0tx1eλcosαteiλtsinαdt=0tx1e(isinαcosα)λtdtλtt====1λx0tx1exp(eiαt)dteiαtt======eiαxλxeiα0tx1etdt=eiαλxΓ(x)(34)(35)(36)(37)(38)\\

 

13.若 b>0 ,证明:
\int^\infty_0x^{-z}\sin(bx)\mathrm dx=\frac\pi2\frac{b^{z-1}}{\Gamma(z)\sin(\pi z/2)}\quad,(0<z<2)\\ \int^\infty_0x^{-z}\cos(bx)\mathrm dx=\frac\pi2\frac{b^{z-1}}{\Gamma(z)\cos(\pi z/2)}\quad,(0<z<1)

利用上题结论变形即得结论

 

14.证明: \mathrm B(np,nq)=n^{-nq}\prod_{i=0}^{n-1}\mathrm B\left(p+\frac{i}{n},q\right)\left[\prod^{n-1}_{i=1}\mathrm B(iq,q)\right]^{-1}

 nnqn1i=0B(p+in,q)[n1i=1B(iq,q)]1=nnqΓ(q)nn1i=0Γ(p+in)Γ(p+q+in)[n1i=1Γ(iq)Γ(q)Γ((i+1)q)]1=nnqΓ(q)n(2π)12(n1)n12npΓ(np)(2π)12(n1)n12n(p+q)Γ(n(p+q))Γ(nq)Γ(q)n=Γ(nq)Γ(np)Γ(n(p+q))=B(np,nq) nnqn1i=0B(p+in,q)[n1i=1B(iq,q)]1=nnqΓ(q)nn1i=0Γ(p+in)Γ(p+q+in)[n1i=1Γ(iq)Γ(q)Γ((i+1)q)]1=nnqΓ(q)n(2π)12(n1)n12npΓ(np)(2π)12(n1)n12n(p+q)Γ(n(p+q))Γ(nq)Γ(q)n=Γ(nq)Γ(np)Γ(n(p+q))=B(np,nq)(39)(40)(41)(42)\\

 

15.证明:
\int^1_0\int^1_0f(xy)(1-x)^{\mu-1}y^\mu(1-y)^{\nu-1}\mathrm dx\mathrm dy\\=\frac{\Gamma(\mu)\Gamma(\nu)}{\Gamma(\mu+\nu)}\int^1_0f(z)(1-z)^{\mu+\nu-1}\mathrm dz

不妨设 f 任意阶可微,故 f 可展开成如下的幂级数

f(x)=\sum^\infty_{n=0}a_nx^n\\

所以

 1010f(xy)(1x)μ1yμ(1y)ν1dxdy=1010n=0ann!xnyn(1x)μ1yμ(1y)ν1dxdy=n=0ann!10xn(1x)μ1dx10yn+μ(1y)ν1dy=n=0ann!B(n+1,μ)B(n+μ+1,ν)=n=0ann!Γ(n+1)Γ(μ)Γ(ν)Γ(n+μ+1)Γ(n+μ+1)Γ(n+μ+ν+1)=Γ(μ)Γ(ν)Γ(μ+ν)n=0ann!B(n+1,ν+μ)=Γ(μ)Γ(ν)Γ(μ+ν)n=0ann!10zn(1z)μ+ν1dz=Γ(μ)Γ(ν)Γ(μ+ν)10(1z)μ+ν+1n=0ann!zndz=Γ(μ)Γ(ν)Γ(μ+ν)10(1z)μ+ν+1f(z)dz 1010f(xy)(1x)μ1yμ(1y)ν1dxdy=1010n=0ann!xnyn(1x)μ1yμ(1y)ν1dxdy=n=0ann!10xn(1x)μ1dx10yn+μ(1y)ν1dy=n=0ann!B(n+1,μ)B(n+μ+1,ν)=n=0ann!Γ(n+1)Γ(μ)Γ(ν)Γ(n+μ+1)Γ(n+μ+1)Γ(n+μ+ν+1)=Γ(μ)Γ(ν)Γ(μ+ν)n=0ann!B(n+1,ν+μ)=Γ(μ)Γ(ν)Γ(μ+ν)n=0ann!10zn(1z)μ+ν1dz=Γ(μ)Γ(ν)Γ(μ+ν)10(1z)μ+ν+1n=0ann!zndz=Γ(μ)Γ(ν)Γ(μ+ν)10(1z)μ+ν+1f(z)dz(43)(44)(45)(46)(47)(48)(49)(50)(51)\\

 

16.若 p+q>1,q<1 ,证明 \int^{\pi/2}_0\cos^{p+q-2}\theta\cos((p-q)\theta)\mathrm d\theta=\frac{\pi}{(p+q-1)2^{p+q-1}\mathrm B(p,q)}

此类积分本人在文章里给出了推广形式,只需代入数据运算即可

 

 

17.证明
(i)当 \Re(s)>0 时, \eta(s)=(1-2^{1-s})\zeta(s)=\frac1{\Gamma(s)}\int^\infty_0\frac{x^{s-1}}{\mathrm e^x+1}\mathrm dx
(ii)当 \Re(s)>1 时, (2^s-1)\zeta(s)=\zeta(s,1/2)=\frac{2^{s-1}}{\Gamma(s)}\int^\infty_0\frac{x^{s-1}}{\sinh x}\mathrm dx

1Γ(s)0xs1ex+1dx=1Γ(s)0xs1n=1(1)n+1enxdx=1Γ(s)n=1(1)n+10xs1enxdxnxx=1Γ(s)n=1(1)n+1ns0xs1exdx=n=1(1)n+1ns=n=11ns2n=11(2n)s=(121s)ζ(s)=η(s)1Γ(s)0xs1ex+1dx=1Γ(s)0xs1n=1(1)n+1enxdx=1Γ(s)n=1(1)n+10xs1enxdxnxx=====1Γ(s)n=1(1)n+1ns0xs1exdx=n=1(1)n+1ns=n=11ns2n=11(2n)s=(121s)ζ(s)=η(s)(52)(53)(54)(55)(56)\\

2s1Γ(s)0xs1sinhxdx=2sΓ(s)0exxs11e2xdx=2sΓ(s)0xs1n=0e(2n+1)xdx=2sΓ(s)n=00xs1e(2n+1)xdx(2n+1)=1Γ(s)n=01(n+1/2)s0xs1exdx=ζ(s,1/2)=(2s1)ζ(s)2s1Γ(s)0xs1sinhxdx=2sΓ(s)0exxs11e2xdx=2sΓ(s)0xs1n=0e(2n+1)xdx=2sΓ(s)n=00xs1e(2n+1)xdx(2n+1)=====1Γ(s)n=01(n+1/2)s0xs1exdx=ζ(s,1/2)=(2s1)ζ(s)(57)(58)(59)(60)(61)\\

 

18.证明: \zeta(s)=\frac{\Gamma(1-s)}{2^{s-1}-1}\frac1{2\pi\mathrm i}\int^{(0+)}_{\infty}\frac{(-z)^{s-1}}{\mathrm e^{z}+1}\mathrm dz

利用Qn21的结论化简即可,不多赘述

 

19.证明:
\log\Gamma(z)=\left(z-\frac12\right)\log z-z+\frac12\log(2\pi)\\+\frac12\sum^\infty_{s=1}\frac{s}{(s+1)(s+2)}\zeta(s+1,z+1)

 (z12)logzz+12log(2π)+12s=1s(s+1)(s+2)ζ(s+1,z+1)=(z12)logzz+12log(2π)+12s=1sΓ(s+3)0tseztet1dt=(z12)logzz+12log(2π)+120eztet1s=1stsΓ(s+3)dt=(z12)logzz+12log(2π)+120ezt[et(t2)+t+2]t2(et1)dt=(z12)logzz+12log(2π)+0(121t+1et1)ezttdt=logΓ(z) (z12)logzz+12log(2π)+12s=1s(s+1)(s+2)ζ(s+1,z+1)=(z12)logzz+12log(2π)+12s=1sΓ(s+3)0tseztet1dt=(z12)logzz+12log(2π)+120eztet1s=1stsΓ(s+3)dt=(z12)logzz+12log(2π)+120ezt[et(t2)+t+2]t2(et1)dt=(z12)logzz+12log(2π)+0(121t+1et1)ezttdt=logΓ(z)(62)(63)(64)(65)(66)(67)\\

 

20.证明: \Phi(z,s,a)=\frac1{\Gamma(s)}\int^\infty_0\frac{t^{s-1}\mathrm e^{-at}}{1-z\mathrm e^{-t}}\mathrm dt

1Γ(s)0ts1eat1zetdt=1Γ(s)0ts1eatn=0znentdt=1Γ(s)n=0zn0ts1e(a+n)tdt(a+n)tt=1Γ(s)n=0zn(a+n)s0ts1etdt=Φ(z,s,a)1Γ(s)0ts1eat1zetdt=1Γ(s)0ts1eatn=0znentdt=1Γ(s)n=0zn0ts1e(a+n)tdt(a+n)tt=======1Γ(s)n=0zn(a+n)s0ts1etdt=Φ(z,s,a)(68)(69)(70)(71)\\

 

21.证明: \Phi(z,s,a)=-\frac{\Gamma(1-s)}{2\pi\mathrm i}\int^{(0{+})}_{\infty}\frac{(-t)^{s-1}\mathrm e^{-at}}{1-z\mathrm e^{-t}}\mathrm dt\quad,|\arg(-t)|<\pi

 Γ(1s)2πi(0+)(t)s1eat1zetdt=Γ(1s)2πi(0+)(t)s1eatn=0znentdt=Γ(1s)2πin=0zn(0+)(t)s1e(a+n)tdttt=Γ(1s)2πin=0zn(0+)ts1e(a+n)tdt(a+n)tt=Γ(1s)2πin=0zn(n+a)s(0+)ts1etdt=n=0zn(n+a)s=Φ(z,s,a) Γ(1s)2πi(0+)(t)s1eat1zetdt=Γ(1s)2πi(0+)(t)s1eatn=0znentdt=Γ(1s)2πin=0zn(0+)(t)s1e(a+n)tdttt====Γ(1s)2πin=0zn(0+)ts1e(a+n)tdt(a+n)tt=======Γ(1s)2πin=0zn(n+a)s(0+)ts1etdt=n=0zn(n+a)s=Φ(z,s,a)(72)(73)(74)(75)(76)(77)\\

 

22.证明:
\Phi(z,s,a)=z^{-a}\Gamma(1-s)\sum_{n\in\mathbb{Z}}(-\log z+2n\pi\mathrm i)^{s-1}\mathrm e^{2n\pi a\mathrm i}\\0<a\leq 1,\Re(s)<0,|\arg(-\log z+2n\pi\mathrm i)|<\pi
\Phi(z,s,a)=\Gamma(1-s)z^{-a}\log^{s-1}(1/z)+z^{-a}\sum^\infty_{r=0}\zeta(s-r,a)\frac{\log(z)^r}{r!}\\|\log z|<2\pi,s≠1,2,3,\cdots;a≠0,-1,-2,\cdots

考虑积分围道 C 以原点为圆心的圆,半径为 (2N+1)\pi,N\in\mathbb{Z}_{>0}

考虑被积函数 f(t):=\frac{(-t)^{s-1}\mathrm e^{-at}}{1-z\mathrm e^{-t}} ,除了一阶极点 \log z+2n\pi\mathrm i(n=-N,\cdots,N)

留数 \mathrm{Res}[f;\log z+2n\pi\mathrm i]=z^{-a}(-\log z-2n\pi\mathrm i)^{s-1}

\frac1{2\pi\mathrm i}\int_C\frac{(-t)^{s-1}\mathrm e^{-at}}{1-z\mathrm e^{-t}}\mathrm dt-\frac1{2\pi\mathrm i}\int^{(0+)}_{(2N+1)\mathrm \pi}\frac{(-t)^{s-1}\mathrm e^{-at}}{1-z\mathrm e^{-t}}\mathrm dt=\\ \sum^N_{n=-N}\mathrm{Res}\left[\frac{(-t)^{s-1}\mathrm e^{-at}}{1-z\mathrm e^{-t}},\log z+2n\pi\mathrm i\right] =z^{-a}\sum^N_{n=-N}(-\log z-2n\pi\mathrm i)^{s-1}\mathrm e^{-2an\pi\mathrm i}

\lim_{N\to\infty}\frac1{2\pi\mathrm i}\int_C\frac{(-t)^{s-1}\mathrm e^{-at}}{1-z\mathrm e^{-t}}\mathrm dt=0\\\lim_{N\to \infty}-\frac1{2\pi\mathrm i}\int^{(0+)}_{(2N+1)\mathrm \pi}\frac{(-t)^{s-1}\mathrm e^{-at}}{1-z\mathrm e^{-t}}\mathrm dt=\frac{\Phi(z,s,a)}{\Gamma(1-s)}\\

所以

Φ(z,s,a)=zaΓ(1s)nZ(logz2nπi)s1e2nπai=zaΓ(1s)nZ(logz+2nπi)s1e2nπaiΦ(z,s,a)=zaΓ(1s)nZ(logz2nπi)s1e2nπai=zaΓ(1s)nZ(logz+2nπi)s1e2nπai(78)(79)\\

接下来看第二个

由常用梅林变换公式知: \mathcal M\{x^{-s}\mathrm e^{-\beta x}\}(\nu)=\beta^{s-\nu}\Gamma(\nu-s),\quad\Re(\nu)>\Re(s),\Re(\beta)>0

作梅林逆变换得 \frac{\mathrm e^{-\beta x}}{x^s}=\beta^s\frac1{2\pi\mathrm i}\int^{\sigma+\mathrm i\infty}_{\sigma-\mathrm i\infty}(\beta x)^{-\nu}\Gamma(\nu-s)\mathrm d\nu ,所以

n=1exp{z(n+a)}(n+a)s=zsn=112πiσ+iσiΓ(νs)dνzν(n+a)ν=zs12πiσ+iσiΓ(νs)zνn=11(n+a)νdν=zs12πiσ+iσiΓ(νs)zνζ(ν,a)dν=zsr=0Res[Γ(νs)zνζ(ν,a),sr(r=0,1,)]+zsRes[Γ(νs)zνζ(ν,a),1]=r=0(1)rζ(sr,a)r!zr+Γ(1s)zs1n=1exp{z(n+a)}(n+a)s=zsn=112πiσ+iσiΓ(νs)dνzν(n+a)ν=zs12πiσ+iσiΓ(νs)zνn=11(n+a)νdν=zs12πiσ+iσiΓ(νs)zνζ(ν,a)dν=zsr=0Res[Γ(νs)zνζ(ν,a),sr(r=0,1,)]+zsRes[Γ(νs)zνζ(ν,a),1]=r=0(1)rζ(sr,a)r!zr+Γ(1s)zs1(80)(81)(82)(83)(84)(85)

两个留数,前者是平凡的,只需利用伽马函数在负整数处的留数即可

后者利用

\mathrm{Res}[\zeta(\nu,a),1]=\lim_{\nu\to 1}\zeta(\nu,a)(s-1)=\lim_{\nu\to 1}\left\{1+\sum^\infty_{n=0}\frac{(-1)^n\gamma_n(a)(s-1)^{n+1}}{n!}\right\}=1

即可

变形得到: \Phi(z,s,a)=\Gamma(1-s)z^{-a}\log^{s-1}(1/z)+z^{-a}\sum^\infty_{r=0}\zeta(s-r,a)\frac{\log(z)^r}{r!}

 

23.Lerch的变换公式
\Phi(z,s,a)=\mathrm iz^{-a}(2\pi )^{s-1}\Gamma(1-s)\left\{\mathrm e^{-\mathrm i\pi s/2 }\Phi\left(\mathrm e^{-2\pi a\mathrm i},1-s,\frac{\log z}{2\pi\mathrm i}\right)\\-\mathrm e^{\mathrm i\pi (s/2+2a) }\Phi\left(\mathrm e^{2\pi a\mathrm i},1-s,1-\frac{\log z}{2\pi\mathrm i}\right)\right\}

 eiπs/2Φ(e2πai,1s,logz2πi)=Γ(s)eiπs/2zanZ(2πai+2nπi)szn=Γ(s)eiπsza(2π)snZzn(n+a)s eiπs/2Φ(e2πai,1s,logz2πi)=Γ(s)eiπs/2zanZ(2πai+2nπi)szn=Γ(s)eiπsza(2π)snZzn(n+a)s(86)(87)

 eiπ(s/2+2a)Φ(e2πai,1s,1logz2πi)=Γ(s)eiπs/2zanZ(2πai+2nπi)szn=Γ(s)za(2π)snZzn(na)s eiπ(s/2+2a)Φ(e2πai,1s,1logz2πi)=Γ(s)eiπs/2zanZ(2πai+2nπi)szn=Γ(s)za(2π)snZzn(na)s(88)(89)

代回去一顿爆算(

 iza(2π)s1Γ(1s){eiπs/2Φ(e2πai,1s,logz2πi)eiπ(s/2+2a)Φ(e2πai,1s,1logz2πi)}=Φ(z,s,a) iza(2π)s1Γ(1s){eiπs/2Φ(e2πai,1s,logz2πi)eiπ(s/2+2a)Φ(e2πai,1s,1logz2πi)}=Φ(z,s,a)(90)

 

24.证明: \frac1{\zeta(s)}=\sum^\infty_{n=1}\frac{\mu(n)}{n^s}

n=1μ(n)ns=p11μ(p)ps=pn=0μ(p)npns=pn=0μ(pn)pns=p(11ps)=1ζ(s)n=1μ(n)ns=p11μ(p)ps=pn=0μ(p)npns=pn=0μ(pn)pns=p(11ps)=1ζ(s)(91)(92)\\

 

25.证明:若 \Re(s)>1 有 \zeta(s)^2=\sum_{n=1}^\infty\frac{\tau(n)}{n^s}

n=1τ(n)ns=p11τ(p)ps=pn=1τ(p)npns=pn=1τ(pn)pns=p1(1ps)2=ζ(s)2n=1τ(n)ns=p11τ(p)ps=pn=1τ(p)npns=pn=1τ(pn)pns=p1(1ps)2=ζ(s)2(93)(94)\\

 

26.证明:若 \Re(s)>1 有 \log\zeta(s)=\sum_p\sum^\infty_{m=1}\frac1{mp^{ms}},-\frac{\zeta'}{\zeta}(s)=\sum^\infty_{n=1}\frac{\Lambda(n)}{n^s}

首先有 \zeta(s)=\prod_{p}\left(1-\frac{1}{p^s}\right)^{-1} ,取对数容易求得

logζ(s)=plog(1ps)=pm=11mpmslogζ(s)=plog(1ps)=pm=11mpms(95)\\

接下来搞第二个

n=1Λ(n)ns=pn=1logppns=ppslogp1ps=Dsplog(1ps)=Dslogζ(s)=ζζ(s)n=1Λ(n)ns=pn=1logppns=ppslogp1ps=Dsplog(1ps)=Dslogζ(s)=ζζ(s)(96)(97)\\

对解析的一些批注

前三题没啥好说的,利用威尔斯特拉斯公式变形即可

第四题是个解析延拓的题,难度也不算大,属于比较典型的解析延拓题了,有个巧妙的方法可以去看 

 的文章

 

第五题没太明白他想问啥,暂时跳过

第六题基于原有结论变形即可,第七题属于第六题的推论

第八题我用渐近分析的方法做的,当然没必要,其他方法留给读者思考

九、十、十一三题都是基于Kummer's Fourier series出的,只要你敢算你就能做出来

第十二题用到了伽马函数的一个积分式 \Gamma(z)=\int^{\infty\mathrm e^{\mathrm i\alpha}}_0t^{z-1}\mathrm e^{-t}\mathrm dt 证明留给读者思考(提示:构造扇形围道)
第十三题是第十二题的推论,其他解法见[1]

第十四题利用倍元公式即可,会算就行

第十五题我怀疑我的解法不太严谨,这里不作说明

第十六题我挂出了文章里的推广结论,并没有参考王竹溪的提示,十分暴力。 

   两位大佬都给出了巧妙的复方法,可以去看看两人的文章

 

第十七题、第二十题都是简单题,我在以前的文章里也给出了很多例子[2]

第十八、二十一题利用Hankel积分即可

第十九题计算罢了

二十二题略有难度,前者涉及到很多复变函数的理论,后者我给出了一个利用梅林变换的很巧妙的解法

第二十三题爆算[3]

第二十四题往后都是披着解析数论外皮的简单题,只需要利用狄利克雷级数的性质[4] L(s,\chi)=\sum^\infty_{n=1}\chi(n)n^{-s}=\prod_p(1-\chi(p)p^{-s})^{-1} 和积性函数的性质即可

“王治淇,著名文字学家,其父、祖父都是清末秀才,自幼熟读四书五经。代表作《新部首大字典》《特殊函数概论》”

(乐)

参考

  1. ^https://www.zhihu.com/question/528516964/answer/2445478028
  2. ^https://zhuanlan.zhihu.com/p/388378336
  3. ^其实这道题我也不知道算没算出来 乐
  4. ^潘承洞、潘承彪《解析数论基础》p.4