Loading [MathJax]/jax/output/CommonHTML/jax.js

I=10log(1+x)1+x2dz=π8log2

함수 f(z)=log(1+z)1+z2을 그림과 같은 경로을 따라 적분을 하자.

경로 1에서

path-1=I

경로 2에서 z=eiθ (θ:0π/2)

1+e2iθ=2cosθeiθlog(1+eiθ)=12log(2+2cosθ)eiθ=12log(2+2cosθ)+iθ

path-2=i4π/2ϵ0log(2+2cosθ)dθcosθπ/2ϵ0θdθcosθ

경로 3에서 z=i+ϵeiθ이므로

path-3=π/20log(1+i+ϵeiθ)ϵ2e2iθ+2iϵeiθiϵeiθdθπ8log(2)iπ216

경로 4에서는 호의 반지름이 ϵ임을 고려하면 z=iy(y:1ϵ0)이므로

path-4=01ϵlog(1+iy)1y2idy=1ϵ0i2log(1+y2)+tan1(y)1y2dy

폐경로 내부에서 f(z)가 analytic하므로 f(z)dz=0이고, 실수부는

I=π8log(2)+π/2ϵ0θdθ4cosθ1ϵ0tan1(y)dy1y2 두번째와 세번째의 적분은 동일하게 logarithmic하게 발산하므로 그 차이는 유한한 값을 가질 수 있다.

θ=π2ϵ: θ4cosθ=π/24ϵ=π8ϵ

y=1ϵ: tan1(y)1y2=π/42ϵ=π8ϵ 그리고 tanθ/2=y로 치환을 하면 cosθ=1y21+y2이므로 두번째 항과 세번째 항이 상쇄됨을 확인할 수 있다. 허수부는 

14π/20log(2+2cosθ)dθcosθπ2161210log(1+y2)dy1y2=0 마찬가지로 위에서와 같이 치환을 하면 등호가 성립함을 보일 수 있다.

728x90

'Mathematics' 카테고리의 다른 글

Integration along a branch cut-035  (0) 2024.11.19
Integrate [log(1+x^2)/1+x^2, {x, 0, ∞}]  (5) 2024.11.15
Integrate [log(1+x)/x, {x, -1, 1}]  (0) 2024.11.13
Integration along a branch cut-034  (0) 2024.11.10
도넛의 표면적, 부피  (1) 2024.11.10
,

I=11log(1+x)xdx=π24

그림과 같은 경로에서 f(z)=log(1+z)z의 적분을 수행하자. z=1f(z)의 simple pole임과 동시에 log(1+z)의 branch point이다. πarg(1+z)π

폐경로 내부에서 analytic하므로 11+rϵf(x)dx+Γ+γϵ=0

반원 Γ에서 z=eiθ (θ:0πϵ)이므로

Γf(z)dz=πϵ0log(1+eiθ)ieiθdθeiθ=iπϵ0[log(2cosθ2)+iθ2]dθ( 1+eiθ=2cosθ2eiθ/2)=iπϵ0log(2cosθ2)dθ(πϵ)24π24as  ϵ0

여기서 π0log(2cosθ2)dθ=0임을 사용했는데, π/20log(sint)dt=π/20log(cost)dt(=π2log(2))를 이용하면 쉽게 보일 수 있다. z=1 둘레의 미소원은 반지름이 rϵ=2sin(ϵ/2)=ϵ이므로 z+1=rϵeiθ (θ:0πϵ2)로 매개화시키면 

γϵf(z)dz=πϵ20log(rϵeiθ)1+rϵeiθieiθdθ=irϵlogrϵπϵ20eiθdθ1+rϵeiθrϵπrϵ20θeiθdθ1+rϵeiθ    0

따라서 ϵ0일 때 I=π24

728x90
,

Action angle variable 이론에서 action 적분을 구해야 하는 경우가 있다. Kepler 궤도 문제일 때 action integral

I=ba(ra)(br)rdr    (0<a<b)을 contour 적분을 이용해서 구하자. 이를 위해서

f(z)=(az)(bz)z을 그림과 같은 경로에서 적분을 하자.

z=a,b는 branch point이다. Cut line을 그림과 같이 선택하는 경우 위상은 πarg(za)π,0arg(bz)2π로 잡는다. z=0이 simple pole이고, 무한대에서도 residue를 가진다. z=0일 때, z=aeiπ, zb=beiπbz=ei2π이므로

Resf(0)=abei3π/2=iab이고(즉 a,b의 기하평균),

(za)(bz)z=i1az1bz=i(1a+b2z+)이므로Resf()=ia+b2 (a,b의 산술평균). 

C1에서 za=(ra)ei0, zb=(br)eiπbz=(br)ei2π이므로 C1=ba(ra)(br)eiπrdr=I C2에서 za=(ra)ei0, zb=(br)eiπbz=(br)ei0이므로 C2=ab(ra)(br)rdr=I

따라서, 

2I=2πi×(ia+b2iab)I=π(a+b2ab)0

728x90
,