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

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

경로 C1C1에서 z=x+i, zi=xeiπ, z+i=x2+4eiθ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
,