How to get an expression for a Hypergeometric function with symbiotic variable

43 views (last 30 days)
I want to calculate the inverse Laplace transform of ,
where G(r)= , s is the variable.
My code can't compute a Hypergeometric function function with symbolic variable, and it can't get a concrete inverse Laplace transformation.
here is my code
phi=4;alpha=2.2;beta=4.7e-05;rho=25;r1=500;r2=400;lambda=4.0e-05;
syms s r t
G_r=0.5*r^2*(hypergeom([phi,-2/alpha],(alpha-2)/alpha,-r^(-alpha)*beta*rho*s )-1);
L_IBd00_s=exp(2*pi*lambda*( vpa(subs(G_r,r,r1),5) -vpa(subs(G_r,r,r2),5) ));
L_IBd00_inv=ilaplace(1/s*L_IBd00_s);
the result is
>> vpa(subs(G_r,r,radius),5)
ans =
125000.0*hypergeom([-0.90909, 4.0], 0.090909, -1.342e-9*s) - 125000.0
>> vpa(L_IBd00_inv,5)
ans =
0.000012253*heaviside(t)*ilaplace(exp(10.0*pi*hypergeom([-0.90909, 4.0], 0.090909, -1.342e-9*s) - 6.4*pi*hypergeom([-0.90909, 4.0], 0.090909, -2.1926e-9*s))/s, s, t)
So how can i get an expression of L_IBd00_inv with respect to the variable t?

Accepted Answer

Torsten
Torsten on 9 Apr 2024 at 17:18
Moved: Torsten on 9 Apr 2024 at 17:18
MATLAB is not able to find the Inverse Laplace transform of such a complicated function as an analytic expression.
  3 Comments
Torsten
Torsten on 10 Apr 2024 at 9:48
What kind of simplified representation of 125000.0*hypergeom([-0.90909, 4.0], 0.090909, -1.342e-9*s) - 125000.0 do you expect to get ? "hypergeom" is defined by an infinite series expression. Do you think this will be easier to handle for "ilaplace" ?
Drawing
Drawing on 11 Apr 2024 at 6:32
Hello, Torsten! I thought i can get a simplified representation like " hypergeom([4.2,-1],2.2,1.1*s ) ans =1 - (21*s)/10", so that this will be easier to handle for "ilaplace" .
I read this opinion from a paper, the author mentioned "the inverse Laplace transform of 1/s*" can be computed directly in MATLAB.

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!