f(t)=L1[1sexp(s1+s2)](t)f(t)=L1[1sexp(s1+s2)](t)

역 Laplace 변환을 Bromwich 적분으로 표현하면 f(t)=12πiγ+iγiF(s)estds=12πiC1zexp(z1+z2+zt)dzf(t)=12πiγ+iγiF(s)estds=12πiC1zexp(z1+z2+zt)dz 적분경로 C={C1+C2+C3+C4+Cϵ+half circle}C={C1+C2+C3+C4+Cϵ+half circle}는 그림의 폐경로에서 오른쪽  z=γz=γ 을 통과하는 수직선을 제외한 나머지 부분이다. 제곱근 함수 때문에 branch cut을 z=±iz=±i을 연결하는 선분으로 잡는다. 이 경우 위상은 π2arg(zi), arg(z+i)3π2π2arg(zi), arg(z+i)3π2 로 잡을 수 있다. 

그러면 각 경로에서,

C1:  z=yeiπ/2, zi=(1y)ei3π/2, z+i=(1+y)eiπ/2, y:10C1:  z=yeiπ/2, zi=(1y)ei3π/2, z+i=(1+y)eiπ/2, y:10

C2:  z=yeiπ/2, zi=(1+y)eiπ/2, z+i=(1y)eiπ/2, y:01C2:  z=yeiπ/2, zi=(1+y)eiπ/2, z+i=(1y)eiπ/2, y:01

C3:  z=yeiπ/2, zi=(1+y)eiπ/2, z+i=(1y)eiπ/2, y:01C3:  z=yeiπ/2, zi=(1+y)eiπ/2, z+i=(1y)eiπ/2, y:01

C4:  z=yeiπ/2, zi=(1y)ei3π/2, z+i=(1+y)eiπ/2, y:10C4:  z=yeiπ/2, zi=(1y)ei3π/2, z+i=(1+y)eiπ/2, y:10

따라서 

C1=01exp(iy1y2+iyt)dyyC1=01exp(iy1y2+iyt)dyy

C2=10exp(iy1y2+iyt)dyy

C3=01exp(iy1y2iyt)dyy

C4=10exp(iy1y2iyt)dyy

그리고  simple pole z=0을 감싸는 적분은 (z=ϵeiθ)

Cϵ=2πi

이므로 f(t)는 다음과 같이 쓸 수 있다.

f(t)=12π10cos(yt)sin(y1y2)dyy

 

728x90
,

f(t)=L1[s2+1s2](t)

Inverse Laplace 변환을 Bromwich 적분으로 표현하기 위해 그림과 같은 폐경로를 고려하면 (C=Ci) ΓF(z)eszdz=(γ+iγiC)F(z)eszdz=2πi×Res(z=0)인데, z=0이 double pole이므로, f(t)=ddz(etzz2+1)|z=0+12πiCz2+1z2etzdz Branch cut을 z=±i을 연결하는 선분 대신 그림과 같이 각각 z=±i에서 출발하는 두 개의 반직선으로 선택한다: πarg(zi), arg(z+i)π

경로 C1에서 z=x+i, zi=xeiπ, z+i=x2+4eiθ

여기서 tanθ=x/2, x:0이다.

cosθ=x/x2+4임을 이용하면

z2+1=(zi)(z+i)=x4x2+4eiπ/2eiθ/2=ix4x2+4(cosθ2+isinθ2)=ix[x2/4+1x/2+ix2/4+1+x/2]ix(g1+ig2)

그리고 C2에서는 zi=xeiπ이므로 (note: z2=x212ix

C1=C2=0ixg1+ig2x212ixexteitdx

마찬가지로 C5,C6에 대한 적분도 (zxi, θθ)

C5=C6=0ixg1ig2x21+2ixdxexteitdx

이므로

f(t)=t2π0dxxext(x21)(g1costg2sint)2x(g1sint+g2cost)(x2+1)2

mathematica V10
mathematica V14

 

728x90
,

I=11(1x)α(1+x)1+αdx5+x2=π30sin[πα+(12α)tan1(5)]sin(πα)(0<α<1)

f(z)=(1z)α(1+z)1+α5+z2 위상은 0arg(1z)2π,πarg(1+z)π로 선택할 수 있다. 그리고 z±=±i5f(z)의 simple pole이다.

C1에서 z1=(1x)eiπ1z=(1x)ei2π, 1+z=(1+x)i0, (x:11)이므로

C1f(z)dz=11(1x)αei2πα(1+x)1+αdx5+x2=ei2παI

C2에서는 z1=(1x)eiπ1z=(1x)i0, 1+z=(1+x)i0, (x:11)이어서

C2f(z)dz=11(1x)α(1+x)1+αdx5+x2=I

따라서 

C1+C2f(z)dz=(ei2πα+1)I=2isin(πα)eiπαI 그리고 C에서는 적분은 0에 수렴한다. z+=i5에서 residue을 구하기 위해서 tanφ=5로 놓으면

arg(z+1)=πφ  arg(1z+)=2πφ arg(1+z+)=φ  αarg(1z+)(1α)arg(1+z+)=[2πα+(12α)φ]

이므로 Resf(i5)=(6)1ei[2πα+(12α)φ]i25=1i230ei[2πα+(12α)φ]

z=i5에서 residue는 

arg(z1)=π+φ  arg(1z)=φ arg(1+z)=φ  αarg(1z)(1α)arg(1+z)=(12α)φ이어서 Resf(i5)=(6)1ei(12α)φi25=1i230ei(12α)φ

두 residue의 합은

Resf(zi)=i230eiπα(ei[πα+(12α)]φei[πα+(12α)φ])=130eiπαsin[πα+(12α)φ]

따라서 residue 정리에 의해서 

2ieiπαsin(πα)I=2πi×(130eiπαsin[πα+(12α)φ])  I=π30sin[πα+(12α)tan1(5)]sin(πα)

728x90
,