Calculate Laplace Transform with Im(k)=0, Im(p)=0

Click For Summary

Homework Help Overview

The discussion revolves around the analytical calculation of a complex Laplace transform involving square roots and integrals. The parameters k and p are specified to have positive real parts, and the original poster expresses difficulty in obtaining results using Mathematica.

Discussion Character

  • Exploratory, Mathematical reasoning, Assumption checking

Approaches and Questions Raised

  • Participants explore various substitutions and transformations to simplify the integrals involved in the Laplace transforms. There is mention of potential connections to Bessel functions and the use of residue theory, although some express skepticism about the applicability of residues due to the nature of the functions involved.

Discussion Status

Several participants are actively engaging with the problem, sharing their thoughts on possible substitutions and transformations. There is no explicit consensus on a solution, but productive lines of inquiry are being explored, including numerical evaluations and theoretical approaches.

Contextual Notes

Participants note that the original problem's constraints and the nature of the integrals may limit the applicability of certain mathematical techniques, such as residue theory, due to the absence of poles in the functions being analyzed.

anv
Messages
4
Reaction score
0
I have a proplem to analitic calculate this :

(k/p^2)*LaplaceTransform[Sqrt[z^2+2p*z]/(z+p)],z,k]+
(1/p^2)*LaplaceTransform[-Sqrt[z^2+2p*z]/(z+p)^2],z,k]

Im(k)=0, k>0
Im(p)=0, p>0

The Mathematica 5 doesn't calculate this.

Very glad to help.
 
Physics news on Phys.org
anv said:
I have a proplem to analitic calculate this :

(k/p^2)*LaplaceTransform[Sqrt[z^2+2p*z]/(z+p)],z,k]+
(1/p^2)*LaplaceTransform[-Sqrt[z^2+2p*z]/(z+p)^2],z,k]

Im(k)=0, k>0
Im(p)=0, p>0

The Mathematica 5 doesn't calculate this.

Very glad to help.

Let me make sure I've got it straight.I'll think about it later.So u want to compute
[tex]I=\frac{k}{p^{2}}L(\frac{\sqrt{z^{2}+2pz}}{z+p})+\frac{1}{p^{2}}L[-\frac{\sqrt{z^{2}+2pz}}{(z+p)^{2}}][/tex]
??That's an ugly "animal"... :-p

Did u try other version of Mathematica ?? :-p

Daniel.

PS.I'll work on it...I smell some Bessel functions...
 
You undestand me right!

I didn't use other version, but I found the new release 5.1 at http://www.wolfram.com
 
[tex]I=L(\frac{\sqrt{z^{2}+2pz}}{z+p})=\int_{0}^{+\infty} \frac{\sqrt{z^{2}+2pz}}{z+p}e^{-sz} dz[/tex](1)

[tex]J=L[-\frac{\sqrt{z^{2}+2pz}}{(z+p)^{2}}]=-\int_{0}^{+\infty} \frac{\sqrt{z^{2}+2pz}}{(z+p)^{2}}e^{-sz} dz[/tex](2)

The first integral,i complete the square under the sq.root and make a substitution:
[tex]I=\int_{0}^{+\infty} \frac{\sqrt{(\frac{z+p}{p})^{2}-1}}{\frac{z+p}{p}} e^{-sz} dz[/tex](3)

And now i make the substitution
[tex]\frac{z+p}{p}\rightarrow \cosh t[/tex](4)
[tex]dz=p\sinh t dt[/tex](5)
The limits of integration are the same ([itex]\arg\cosh 1 =0[/itex]).
The exponential becomes:
[tex]e^{-sz}=e^{sp}e^{-sp\cosh t}[/tex](6)

Therefore
[tex]I=p e^{sp}\int_{0}^{+\infty} e^{-sp\cosh t}\frac{\sinh^{2} t}{\cosh t} dt[/tex] (7)

Can u convince 'Mathematica' to evaluate this integral??Maybe numerically...

Using the same kind of substitution,for evaluating the second integral (L.transform),u get
[tex]J=e^{sp} \int_{0}^{+\infty} e^{-sp\cosh t}\frac{\sinh^{2} t}{\cosh^{2} t} dt[/tex] (8)

Again,i don't know what to do to it.

Daniel.
 
Last edited:
Thx for help, but I don't simplify your answer.

I'm going use the answer (if it possible) in math model.
Numerical answer is suit, but if it in the form:
f(k,p)*numerical_answer

I try to apply the teory residues for a solve.
What do you think about?

P.S. I should note that the root of initial problem is


\int {0} {\inf} { exp(-rk)/r*(1/r+k) }dt

r = sqrt ( p^2 + (vt)^2 )
 
without ambiguity:

\int {0} {\inf} { (exp(-rk)/r)*(1/r+k) }dt

r = sqrt ( p^2 + (vt)^2 )


p.s. my calculation is big, that is why I don't insert its.
 
anv said:
I have a proplem to analitic calculate this :

(k/p^2)*LaplaceTransform[Sqrt[z^2+2p*z]/(z+p)],z,k]+
(1/p^2)*LaplaceTransform[-Sqrt[z^2+2p*z]/(z+p)^2],z,k]

Im(k)=0, k>0
Im(p)=0, p>0

The Mathematica 5 doesn't calculate this.

Very glad to help.


I'm sorry,but it doesn't work with residues,because the functions doesn't have poles."Im(p)=0, p>0",and the integration is not on entire R,but only on its positive semiaxis.,where the denominator is never zero.

Daniel.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 5 ·
Replies
5
Views
1K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
6
Views
2K
  • · Replies 11 ·
Replies
11
Views
2K
Replies
3
Views
3K
Replies
10
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 14 ·
Replies
14
Views
2K