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

함수 f(x)의 Fourier transform:

˜f(k)=1(i(xiϵ))αeikxdx=|k|α1Γ(1α)sin(πα)(1+sgn(k))여기서 ϵ0+은 미소양수이고, 0<α<1는 실수이다. 함수 f(z)=1(i(ziϵ))αz=iϵ가 branch point이므로 그림과 같이 cutline을 잡는다: 3π2arg(z)π2.

k<0일 때는 lower half plane을 도는 적분경로를 잡으면 ˜f=0임을 보일 수 있고(k<0에서 0이 아닌 결과를 얻기 위해서는 xiϵx+iϵ으로 치환한다), k>0인 경우는 그림과 같이 upper half plane에서 cutline을 감싸는 반원경로에서 적분을 고려하자.

이 경로 내에서 analytic 하므로 f(z)의 Fourier transform은 cutline을 감싸는 두 경로 C1C2에서 적분으로 표현할 수 있다. Branch point를 감싸는 미소원호에서 적분은 0으로 수렴한다.

경로 C1에서 z=yeiπ/2 (y:ϵ), dz=idy,C1=ϵeiπα/2eαlogyiπα/2eky(idy)=ieiπα0ekyyαdy 경로 C2에서는 z=yei3π/2 (y:ϵ), dz=idyC2=ϵeiπα/2eαlogy+i3πα/2eky(idy)=ieiπα0ekyyαdy 이 두 경로에서의 적분합은 C1+C2=i(eiπαeiπα)0ekyyαdy=2sin(πα)0ekyyαdy=2kα1Γ(1α)sin(πα)(k>0) 따라서 그림의 경로에서 f(z)eikz가 analytic 하므로 

f(z)eikzdz=0f(x)eikxdx=C1+C2f(z)eikzdz k<0일 때 0임을 고려하면, 1(i(xiϵ))αeikxdx=|k|α1Γ(1α)sin(πα)(1+sgn(k)) f(x)=(i(x+iϵ))α인 경우에는 lower half plane(k<0)에서 적분경로를 선택하면 된다.  1(i(x+iϵ))αeikxdx=|k|α1Γ(1α)sin(πα)(1+sgn(k))Γ(α)Γ(1α)=πsin(πα)임을 이용하면 

1(±i(xiϵ))αeikxdx=2π|k|α1Γ(α)θ(±k)

728x90

'Mathematics' 카테고리의 다른 글

Integration along a branch cut-042  (0) 2024.12.02
Shokhotski-Plemelj theorem  (0) 2024.11.26
Integration along a branch cut-040  (0) 2024.11.24
Integration along a branch cut-039  (0) 2024.11.23
Integration along a branch cut-038  (0) 2024.11.22
,

 

 

I=π/2π/2cos(aθ)cosb(θ)dθ(a>b>1)

먼저 cosθ=(eiθ+eiθ)/2임을 고려하면 

I=Reπ/2π/2eiaθ(eiθ+eiθ2)2dθ로 쓸 수 있다. 이는 복소평면에서 π/2θπ/2인 반원(C0:|z|=1,π/2arg(z)π/2)에서 적분으로 생각할 수 있다. z=eiθ놓으면 dθ=dz/iz이므로

I=Reright half circleza(z2+12z)bdziz=12bImright half circlezab1(z2+1)bdz 따라서 f(z)=zα(1+z2)β(α>1, β>1)를 그림과 같은 폐경로를 반시계방향으로 순환하는 선적분을 고려하자.


z=0,±if(z)의 branch point이므로 그림에 표시된 cutline을 선택한다. 그러면 πarg(z)π, π2arg(z+i), arg(zi)3π2로 잡을 수 있다. z=0 (z=ϵeiθ)z=±i (z=±i+ϵeiθ)을 감싸는 원호에서 적분의 기여는 없음을 쉽게 알 수 있다. C1에서 적분은

z=yeiπ/2  (y:10)

z+i=(1+y)eiπ/2, zi=(1y)eiπ/2이므로 

C1=01yαeiπα/2(1y2)β(eiπ/2dy)=eiπ(α+1)/210yα(1y2)βdy

C2에서는 z=yeiπ/2 (y:01)z+i=(1+y)eiπ/2, zi=(1y)eiπ/2이므로

C2=10yαeiπα/2(1y2)β(eiπ/2dy)=eiπ(α+1)/210yα(1y2)βdy

그림의 경로에서 f(z)가 analytic하므로 

f(z)=0

  C0f(z)dz=C1+C2f(z)dz=2isinπ(α+1)210yα(1y2)βdy이고 α=ab1, β=b이므로

I=12bImC0f(z)dz=12bsinπ(ab)210yab1(1y2)bdy

y2=t로 치환하면 

I=sinπ(ab)22b10t(ab)/21(1t)bdt이고 Γ(z)의 정의를 사용하면 I=sinπ(ab)22bΓ(ab2)Γ(b+1)Γ(a+b2+1)=π2bΓ(b+1)Γ(a+b2+1)Γ(ba2+1)=π2b(ba+b2)여기서 중간 등호는 Γ(m)Γ(1m)=πsin(mπ/2)임을 사용했다. Check

a=1,b=0:  π/2π/2cosθdθ=2I=πΓ(1)Γ(1+12)Γ(12)=π112(π)2=2

728x90

'Mathematics' 카테고리의 다른 글

Shokhotski-Plemelj theorem  (0) 2024.11.26
Integration along a branch cut-041  (0) 2024.11.25
Integration along a branch cut-039  (0) 2024.11.23
Integration along a branch cut-038  (0) 2024.11.22
Integration along a branch cut-037  (0) 2024.11.21
,

복소함수 F(z)의 inverse Z-transform은

Z1[F(z)]=12πiF(z)zndz    (n=0,1,2,...)

로 정의되는데, F(z)=(z1z)α    (α=Real)인 경우에 대해 구하자. z=0z=1F(z)의 branch point이므로 그림과 같은 dog bone 모양의 경로(반시계 방향)를 선택한다. 

그러면 πarg(z),arg(z+1)π로 선택할 수 있다.  z=0z=1을 감싸는 미소원에서 적분의 기여가 없고, C1에서는 z=xeiπ, z+1=(1x)ei0,  (x:10)이므로 

C1=01(1xx)αeiπαxn1eiπ(n1)eiπdx=eiπ(nα)10xn1α(1x)αdx 그리고 C2에서는 z=xeiπ, z+1=(1x)eiπ,  (x:01)이므로

C2=10(1xx)αeiπαxn1eiπ(n1)eiπdx=eiπ(nα)10xnα1(1x)αdx 따라서

12πiF(z)zndz=sin[π(nα)]π10xnα1(1x)αdx그리고 Gamma 함수 Γ(z)를 사용하면,

12πiF(z)zndz=sin[π(nα)]πΓ(nα)Γ(α+1)Γ(n+1)로 쓸 수 있고, 또 다음 관계

Γ(m)Γ(1m)=πsin(πm)를 이용하면

12πiF(z)zndz=Γ(α+1)Γ(n+1)Γ(αn+1)=(αn)

728x90

'Mathematics' 카테고리의 다른 글

Integration along a branch cut-041  (0) 2024.11.25
Integration along a branch cut-040  (0) 2024.11.24
Integration along a branch cut-038  (0) 2024.11.22
Integration along a branch cut-037  (0) 2024.11.21
Integration along a branch cut-036  (0) 2024.11.20
,