# Problem computing inverse Laplace transform of the Bessel function

조회 수: 30(최근 30일)
Kevin B 2017년 11월 11일
답변: BeeTiaw 2020년 8월 5일
I am trying to evaluate the inverse Laplace transform of a function that includes the modified Bessel function of the second kind, i.e., my function is f_s = sqrt(a*s)*besselk(1,2*sqrt(a*s)), where a is a constant (say a = 0.5) and 's' is defined as a symbol. To evaluate the inverse Laplace transform of this function, I used the built-in Matlab function 'ilaplace' as transV = ilaplace(f_s,s,t), where I want to evaluate my function at 't' = 2. However, the problem I am getting is that the inverse transformed expression transV still contains 's' variable, which I supposed to transform to 't' after the inverse Laplace transform.
Can anyone help me with what is the problem with this approach? I checked the Matlab documentation and still couldn't figure out the solution. I already tried using 'vpa' function as it was suggested in one of the Matlab forums, but that didn't help to get rid of 's'. I am stuck with this and would really appreciate any help on it.
My code is
a = 0.5; t = 2;
syms s;
f_s = 1/s*sqrt(a*s)*besselk(1,2*sqrt(a*s));
transV = ilaplace(f_s,t); % or transV = vpa(ilaplace(f_s,t));
With this, the final result transV still has an expression which is a function of 's' although it should be transformed after ilaplace.
##### 댓글 수: 2표시숨기기 이전 댓글 수: 1
Kevin B 2017년 11월 11일
Hi David,
Thanks for your suggestion. Yes, the closed form solution seems not to exist. I tried Mathematica to check if it can help but without success. For me, the numerical solution is fine for now and I am not able to get that either. Do you have any further suggestions to get the numerical solution? Thanks a ton!

댓글을 달려면 로그인하십시오.

### 채택된 답변

Walter Roberson 2017년 11월 11일
Maple says that for a > 0, the result is (1/2)*exp(-a/t)
##### 댓글 수: 3표시숨기기 이전 댓글 수: 2
Kevin B 2017년 11월 12일
Thank you so much!

댓글을 달려면 로그인하십시오.

### 추가 답변(1개)

BeeTiaw 2020년 8월 5일

댓글을 달려면 로그인하십시오.

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!